POST
|
Thats a job for network analyst if you are planning to use linear features, but if you can convert the roads to raster you can do a cost path using spatial analyst. Those are the simplest but may require you to make a purchase if your organization doesnt already have one of those.
... View more
08-12-2010
05:07 AM
|
0
|
0
|
193
|
POST
|
I use komodo edit as my ide. I really dont like pythonwin. Couldnt tell you why, just doesnt feel right. Try a different ide and see if it installs.
... View more
08-11-2010
05:26 AM
|
0
|
0
|
142
|
POST
|
Well you could try gp.delete_management(path, "shapefile") which tells the tool explicitly you are deleting a shapefile, or you can try deleting it with os.remove(). With os.remove() you will have to delete all the files that make up the shapefile individually though. os.remove() is pure python and doesnt muck about with the geoprocessor. for ext in [".shp",".shx",".sbn",".dbf",".prj",".sbx"]: os.remove("C:\\Temp\\Contour"+ext)
... View more
08-11-2010
05:19 AM
|
0
|
0
|
796
|
POST
|
The search cursor will only apply to the single feature class given when the object is created. So the script deletes the file but still has an error? Could you post the error? The while shouldnt fall into an infinite loop with such a simple function inside it. 'while row:' is the same as 'while row != none:'. The loop will continue until cursor.Next() returns a result of none. The while row: loop of geoprocessing cursors work like a for row in cursor: on a normal python cursor.
... View more
08-06-2010
04:53 AM
|
0
|
0
|
796
|
POST
|
I cant create an update cursor on a joined table. This is the error given to me. I can create the cursor on an un-joined table though, but that isnt really the point. import arcgisscripting
gp=arcgisscripting.create(9.3)
gp.makefeaturelayer_management(r"C:\GIS Projects\CRW\TrakIt_Source_RS.gdb\Target_Data\parce_test", "parce_test")
gp.addjoin("parce_test", "A1RENUM", r"C:\GIS Projects\CRW\parce_test_spatial_join.shp", "A1RENUM")
cursor=gp.updatecursor("parce_test")
row=cursor.Next()
while row:
row.Flood_zn = row.getvalue("parce_test_spatial_join.FLD_ZONE")
cursor.updaterow(row)
row=cursor.Next()
del row
del cursor
Traceback (most recent call last): File "C:/GIS Projects/CRW/upgradetest", line 6, in <module> cursor=gp.updatecursor("parce_test") RuntimeError: ERROR 999999: Error executing function.
... View more
08-05-2010
08:09 AM
|
0
|
0
|
848
|
POST
|
I think you are going too deep in your workspace. Is '\\Cl ip_Feature' a directory? I thought that was the name of a shapefile. Also you cant do a for loop on at 9.3 cursor. You have to use a while loop. It seems strange but that is changed at 10 so you can use for.
... View more
08-04-2010
05:16 AM
|
0
|
0
|
796
|
POST
|
Thats awesome. If you have two two features, each in a separate feature class, aFeature and bFeature, could this method of checking polygon intersection be used to see if aFeature intersects aFeature and if so write an attribute of bFeature to a field in aFeature?
... View more
08-03-2010
12:28 PM
|
0
|
0
|
2286
|
POST
|
You need a search cursor to look over the shapefile and see if that value occurs. in 9.3: scursor=gp.SearchCursor(Clip_Feature)
row=scursor.Next()
while row:
if row.GetValue("Region_No") == 3:
gp.Delete_management("C:\\Daves_Stuff\\ArcPad_Data \\Standard_Layers\\Lamberts_Projection\\Pts_Lines\\Contour.shp")
del scursor
in 10: for row in arcpy.SearchCursor(Clip_Feature):
if row.Region_No == 3:
arcpy.Delete_management("C:\\Daves_Stuff\\ArcPad_Data \\Standard_Layers\\Lamberts_Projection\\Pts_Lines\\Contour.shp")
... View more
08-03-2010
04:55 AM
|
0
|
0
|
796
|
POST
|
Really? How do you reference the field names ("name",[name],etc)? I have never had an update cursor work on a joined table before.
... View more
08-02-2010
10:44 AM
|
0
|
0
|
848
|
POST
|
Can we use an update cursor on a joined table in arcpy? I dont have 10 yet (very sad) because of a 3rd party plug-in conflict. I should have it soon so I am updating scripts for arcpy and this abilty would help me greatly.
... View more
08-02-2010
08:06 AM
|
0
|
17
|
5054
|
POST
|
Think of arcgisscripting as a com object that calls out to the geoprocessor. When you import arcpy you dont import the geoprocessor just classes that send data into it and get results back from it. The actual geoprocessor is part of the ArcGIS core. If one computer has 10 already why not just run the scripts from there?
... View more
07-23-2010
05:03 AM
|
0
|
0
|
619
|
POST
|
No because arcpy requires ArcGIS 10 geoprocessor to be installed. The 9.3.* geoprocessor interfaces with python 2.5 and the 10 geoprocessor interfaces with python 2.6 so there is some version disconnect there as well.
... View more
07-22-2010
07:25 AM
|
0
|
0
|
619
|
POST
|
I cant imagine why this would be causing a problem, but could you try modifiying the model so that you dont use a feature layer? If the shapefile is already loaded into ArcMap you shouldnt need to make it a feature layer in the model.
... View more
07-22-2010
07:15 AM
|
0
|
0
|
1228
|
POST
|
If it works outside of arcmap I assume you are using ArcGIS 10. Try selecting a feature by hand and then pasting each line from the script into the python window in ArcMap one at a time. There may be a typo in the code that will rear its head when you are using the code interactivly. Are there multiple data frames in the map?
... View more
07-21-2010
11:23 AM
|
0
|
0
|
1228
|
POST
|
Anyone have an idea of a way in python that would let me take the centroid of a feature and then get the values of intersecting features from other feature classes to put into attributes for the centroid? I dont want to use Spatial Join because I dont want a bunch of intermediate files or to deal with making feature layers to join those files back to the original feature class to calc over the data. The way I am doing it now is: parcel centroids are spatially joined to a polygon feature class like zoning join field permanently joins the pertinent fields from centroids_sj to the parcels feature class the added field is calculated over to another field in the parcels schema and the added field is deleted This is wildly inefficient because Join Field takes a really long time to run all 27 times in this script. If I could go over the parcels with a cursor and for each row get the centroid with describe and just intersect it with other features I think it could speed up this code.
... View more
07-20-2010
09:18 AM
|
0
|
0
|
700
|
Title | Kudos | Posted |
---|---|---|
1 | 11-03-2015 11:33 AM | |
1 | 09-23-2015 10:50 AM | |
1 | 12-14-2015 10:42 AM | |
1 | 12-14-2015 01:10 PM | |
1 | 12-04-2015 08:19 AM |
Online Status |
Offline
|
Date Last Visited |
11-11-2020
02:23 AM
|