วิธีแก้ปัญหาสำหรับสิ่งที่คุณอธิบายเมื่อไม่มีนักวิเคราะห์เชิงพื้นที่:
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 ลงใน lines_3D และเติมโดยใช้:
!XY!.split(" ")[0] and !XY!.split(" ")[1]
สร้างเลเยอร์ XY จากฟิลด์ด้านบน 2 คัดลอกไปยังจุดคุณลักษณะคลาสและใช้ชี้ไปที่บรรทัดเพื่อสร้างบรรทัดบัฟเฟอร์ภายในต้นฉบับ:
เอาท์พุท:
สรุป:
เว้นเสียแต่ว่าภูมิประเทศของคุณจะเป็นบัฟเฟอร์แบบง่าย ๆ ที่ชันมากจะทำงานได้เว้นแต่ว่าคุณจะพูดถึงความยาวของการไหล โชคไม่ดีเพราะการที่คุณต้องการกล่องเครื่องมือการวิเคราะห์ทางอุทกวิทยาจากนักวิเคราะห์เชิงพื้นที่