Solution de contournement pour ce que vous décrivez en l'absence d'analyste spatial:
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")
Exécutez ceci sur un nouveau champ (fixé le 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!)
Ajoutez des champs X et Y à lines_3D et remplissez-les en utilisant:
!XY!.split(" ")[0] and !XY!.split(" ")[1]
Créez une couche XY à partir des 2 champs ci-dessus, copiez-la dans la classe d'entités points et utilisez point à ligne pour créer une ligne tampon à l'intérieur de l'original:
PRODUCTION:
Conclusion:
À moins que votre terrain ne soit extrêmement raide, un simple tampon fera l'affaire, à moins que cela ne semble très probable que vous parliez réellement de la longueur du flux. Pas de chance, car pour cela vous avez besoin d'une boîte à outils d'analyse hydrologique d'analyste spatial