POST
|
I'm also having this exact problem in 10.1 SP1. Does anyone know of a fix, or a patch?
... View more
01-07-2014
10:58 PM
|
0
|
0
|
385
|
POST
|
Thanks Chris, hopefully it should make a bit more sense now.
... View more
12-23-2010
01:08 AM
|
0
|
0
|
563
|
POST
|
Hi, I'm new to python and have just made my first script but am very disapointed with the performance acheived, not sure if this is a problem with the script or inherent? First my script...... I'm working on a problem of spatially joining road edge vertices to the nearest road centreline vertex. The potential problem encountered here is where wide road edges are closer to an adjacent roads centreline than their own and hence will be joined to the wrong road. To solve this I draw a line between two candidate vertices and proceed only if this line is completely within the road area (the joining line does not jump between roads). Here is an exerpt of the script which is looping through the road edge vertices: #start the counter. Used to select vertices sequentially
count = 1
#Setup search cursor for the road outline vertices
rows1 = arcpy.SearchCursor("RoadOutline_Points_Layer")
#Loop through each road outline vertex joining each to its nearest intrinsic centreline vertex
for row1 in rows1:
#create a variable used to select outline (OL) vertices by their unique ID
id_sel = "OLID = %f" % count
#Select a point based upon its ID
arcpy.SelectLayerByAttribute_management("RoadOutline_Points_Layer","NEW_SELECTION",id_sel)
#Calculate the distance between each road centreline vertex and the selected outline point
arcpy.Near_analysis(p_centrelinepoints,"RoadOutlin e_Points_Layer")
#Sort the centreline vertices by this calculated distance and export to a new featureclass
arcpy.Sort_management(p_centrelinepoints,"RoadCentreline_Points_Near","NEAR_DIST")
#Delete the near analysis fields from the centre line points
arcpy.DeleteField_management(p_centrelinepoints,["NEAR_DIST","NEAR_FID"])
#Create a layer from the new FC
arcpy.MakeFeatureLayer_management("RoadCentreline_ Points_Near","RoadCentreline_Points_Near_Layer")
#Setup a loop for the centreline vertices so that starting with the nearest point they are interated through until the conditions are met. the conditions being that the line between the two points is intrinsic to the road area
#Variables for loop control and counter respectively
nearloop = 0
nearloop_count = 0
#loop until the conditions are not met
while nearloop == 0:
#Increment the counter
nearloop_count = nearloop_count + 1
#create a varaible to select the centreline vertex by unique ID (CLID)
objectid_sel = "OBJECTID = %f" % nearloop_count
#Select the sorted record
arcpy.SelectLayerByAttribute_management("RoadCentr eline_Points_Near_Layer","NEW_SELECTION",objectid_ sel)
#Merge with the selected outline point to give a new FC. This gives a feature class with the road outline vertex and a candidate centreline vertex
arcpy.Merge_management(["RoadCentreline_Points_Near_Layer","RoadOutline_Po ints_Layer"],"temp_pair_points")
#Draw a line between these two points
arcpy.PointsToLine_management("temp_pair_points"," temp_pair_line")
#Erase the line by the road area. If all of the line is erased then the path between the two points is intrinic to the road and so is ok to spatially join.
arcpy.Erase_analysis("temp_pair_line",p_area,"temp _pair_line_Erase",p_precision)
#Test to see if all of the road has been erased by finding the number of records in the erased layer
erasecount = arcpy.GetCount_management("temp_pair_line_Erase")
erasecount = int(erasecount.getOutput(0))
#Test the erased layer. If there are no erase features the point is good and can be joined with the outline point, otherwise the loop continues.
if erasecount == 0:
#The conditions are met stop the loop
nearloop = 1
#Join the two points
#Create a layer from the feature class containing the two points
arcpy.MakeFeatureLayer_management("temp_pair_point s","temp_pair_points_layer")
#Select the centreline point
arcpy.SelectLayerByAttribute_management("temp_pair _points_layer","NEW_SELECTION","CLID IS NOT NULL")
#Setup a search cursor for this featureclass
rows2 = arcpy.SearchCursor("temp_pair_points_layer")
for row2 in rows2:
#get the unique ID of the selected centreline vertex
joinCLID = row2.getValue("CLID")
#Write this(Join) to the outline vertex
arcpy.CalculateField_management("RoadOutline_Point s_Layer","CLID",joinCLID)
else:
#The conditions are not met continue the loop
nearloop = 0
#This outline vertex is complete, increment the counter by 1 to start the next point
count = count + 1 I know this is pooly constructed code but it does work. My problem is that the performance is terrible. For each outline vertex the above process takes approx. 50 seconds on a decent machine. 1 sqkm of data will have about 9000 road edge vertices which would take around 5 days....I don't understand how it is taking so long it is not to far from the individual process being perfromed manually in real time. Looking at the processes the 'heaviest' is the near analysis, but i timed this (from toolbox) and it took 5 seconds. I have timed all of the processes inddividually and for each point they add up to about 11 seconds...any ideas as to why it is taking 5x longer through python? Incidentally the nested loop which searches through to find a correctly joined centreline almost always only iterates through once as the closest is almost always ok, so this will not be eating up time on a regular occurance. Any help greatly appreciated. Andrew
... View more
12-22-2010
06:12 AM
|
0
|
5
|
3290
|
Online Status |
Offline
|
Date Last Visited |
11-11-2020
02:23 AM
|