Używam ArcInfo 10 SP3.
Pracuję nad reorganizacją naszych danych o użyteczności. Dwa lata temu zaczęliśmy zbierać linię wodną Private Service. Wciąż mamy wiele z nich do wyciągnięcia ze starych Rysunków.
Zastanawiałem się, czy istnieje sposób na utworzenie linii, która łączyłaby ślady naszego budynku z linią WaterMain?
Jako punkt wyjścia chciałbym wykorzystać wierzchołek budynku znajdujący się najbliżej magistrali wodnej.
Jeśli szukasz rozwiązania, które nie wymaga opracowania narzędzia .NET, możesz użyć poniższego skryptu python, aby osiągnąć dokładnie to, czego szukasz. Miałem dokładnie taką samą potrzebę i napisałem następujący skrypt jako rozwiązanie. Skonfiguruj go jako narzędzie ArcCatalog z 4 parametrami lub skomentuj parametry i odkomentuj na stałe zmienne i uruchom je bezpośrednio.
# CreateLineFromNearestVertexToFeature.py
# Author: Jeff Berry
# Description: Creates a line between the nearest vertext on source features
# to the nearest feature in target feature class.
# ---------------------------------------------------------------------------
# Import arcpy module
import arcpy
from arcpy import env
# Local variables:
# 1. SourceFC - Feature Class
# 2. TargetFC - Feature Class
# 3. Output_gdb - Geodatabase
# 4. Output_fc - String
SourceFC = arcpy.GetParameterAsText(0)
TargetFC = arcpy.GetParameterAsText(1)
Output_gdb = arcpy.GetParameterAsText(2)
Output_fc = arcpy.GetParameterAsText(3)
## Alternatively setup hardcoded variables
##SourceFC = "Buildings"
##TargetFC = "WaterMains"
##Output_gdb = "D:\\New File Geodatabase.gdb"
##Output_fc = "lines_output"
SourceFeaturePoints = "SrcFtrPoints"
arcpy.env.workspace = Output_gdb
# Process: Feature Vertices To Points
arcpy.FeatureVerticesToPoints_management(SourceFC, SourceFeaturePoints, "ALL")
# Process: Near
arcpy.Near_analysis(SourceFeaturePoints, TargetFC, "1000 Feet", "LOCATION", "NO_ANGLE")
# Process: Create Feature Class...
#arcpy.CreateFeatureclass_management(Output_gdb, Output_fc, "POLYLINE", "", "DISABLED", "DISABLED", "", "", "0", "0", "0")
rows = arcpy.SearchCursor(SourceFeaturePoints)
lstIDs = []
for row in rows:
lstIDs.append(row.ORIG_FID)
uniqueOBJIDS = set(lstIDs)
newLineList = []
shapeName = arcpy.Describe(SourceFeaturePoints).shapeFieldName
for objID in uniqueOBJIDS:
rows = arcpy.SearchCursor(SourceFeaturePoints, "\"NEAR_DIST\" = (SELECT MIN( \"NEAR_DIST\") FROM SrcFtrPoints WHERE \"ORIG_FID\" = " + str(objID) + ")")
for row in rows:
arrayLine = arcpy.Array()
ftr = row.getValue(shapeName)
pointStart = ftr.firstPoint
pointEnd = arcpy.Point(row.NEAR_X, row.NEAR_Y)
arrayLine.add(pointStart)
arrayLine.add(pointEnd)
plyLine = arcpy.Polyline(arrayLine)
newLineList.append(plyLine)
arcpy.CopyFeatures_management(newLineList, Output_fc)
arcpy.Delete_management(SourceFeaturePoints, "FeatureClass")
del rows
del row
del SourceFeaturePoints
del Output_fc
del Output_gdb
arcpy.ClearEnvironment("workspace")Zajrzyj do metody „NearestFeature” na IIndexQuery2.
Możesz użyć tego, aby uzyskać najbliższą główną funkcję wody dla każdego budynku. Sądzę, że musiałbyś w jakiś sposób przejść przez wierzchołki każdego budynku, aby znaleźć najbliższą odległość od tej funkcji, a następnie zbudować nową polilinię, używając wierzchołków z budynku i głównego budynku jako punktów końcowych. Jedyny raz, kiedy to zrobiłem, użyłem dwóch punktowych klas funkcji, chciałbym móc zaoferować coś więcej niż z góry mojej głowy ..: D
IFeatureCursor pDepthCursor = pDepthSoundings.Search(null, false);
IFeatureIndex2 pFtrInd = new FeatureIndexClass();
pFtrInd.FeatureClass = pDepthSoundings.FeatureClass;
pFtrInd.FeatureCursor = pDepthCursor;
pFtrInd.Index(null, pCombinedEnvelope);
IIndexQuery2 pIndQry = pFtrInd as IIndexQuery2;
int FtdID = 0;
double dDist2Ftr = 0;
pIndQry.NearestFeature(ppoint, out FtdID, out dDist2Ftr);
IFeature pCloseFeature = pDepthSoundings.FeatureClass.GetFeature(FtdID);