Обходной путь для того, что вы описываете в отсутствие пространственного аналитика:
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 полей выше, скопируйте его в класс объектов точек и используйте пункт в линию, чтобы создать строку буфера внутри оригинала:
ВЫХОД:
Вывод:
Если ваш рельеф местности не очень крутой, простой буфер будет выполнять эту работу, если, скорее всего, вы не говорите о длине потока. Не повезло, потому что для этого вам нужен инструментарий гидрологического анализа от пространственного аналитика