स्थानिक विश्लेषक की अनुपस्थिति में आपके द्वारा वर्णित के लिए वर्कअराउंड:
arcpy.Buffer_analysis("target","../buffer.shp", "100 Meters")
arcpy.FeatureVerticesToPoints_management("buffer","../points.shp","ALL")
arcpy.AddField_management("points", "PID", "LONG")
arcpy.CalculateField_management("points", "PID", "[FID]")
arcpy.Near_analysis("points", "target","LOCATION")
arcpy.MakeXYEventLayer_management("points","NEAR_X","NEAR_Y","points_Layer")
arcpy.CopyFeatures_management("points_Layer","../from_points.shp")
arcpy.Merge_management("from_points;points","../pairs.shp")
arcpy.PointsToLine_management("pairs","../lines_2D.shp", "PID")
arcpy.InterpolateShape_3d("DEM","lines_2D","../lines_3D.shp")
arcpy.AddField_management("lines_3D", "XY", "TEXT")
इसे नए फ़ील्ड पर चलाएँ (05/02/2018 को निर्धारित):
def getPoint(shp):
part=shp.getPart(0)
n=len(part);L=0
for i in xrange(n):
p=part.getObject(i)
x=p.X;y=p.Y;z=p.Z
if i >0:
dX=X-x;dY=Y-y;dZ=Z-z;dL=pow((dX*dX+dY*dY+dZ*dZ),0.5);L+=dL
if L>=100: break
X,Y,Z=x,y,z
xNew=x+dX/dL*(L-100);yNew=y+dY/dL*(L-100)
return '%s %s' %(xNew,yNew)
#-----------------
getPoint( !Shape!)
X और Y फ़ील्ड को लाइनों_3D में जोड़ें और उनका उपयोग करके पॉपुलेट करें:
!XY!.split(" ")[0] and !XY!.split(" ")[1]
2 क्षेत्रों के ऊपर से XY लेयर बनाएं, इसे अंक वर्ग में कॉपी करें और मूल के अंदर बफर लाइन बनाने के लिए पॉइंट टू लाइन का उपयोग करें:
उत्पादन:
निष्कर्ष:
जब तक आपका इलाक़ा बेहद खड़ी साधारण बफर काम नहीं करेगा, जब तक कि बहुत संभावना नहीं लगती कि आप वास्तव में प्रवाह की लंबाई के बारे में बात कर रहे हैं। बुरी किस्मत, क्योंकि इसके लिए आपको स्थानिक विश्लेषक से हाइड्रोलॉजिकल विश्लेषण टूलबॉक्स की आवश्यकता होती है