Criando linha (vértice mais próximo da linha) usando o ArcGIS Desktop?

11

Estou usando o ArcInfo 10 SP3.

Estou trabalhando para reorganizar nossos dados de utilidade. Há dois anos, começamos a coletar a linha de água do Serviço Privado. Ainda temos muitos deles para extrair dos antigos Desenhos de Registro.

Eu queria saber se havia uma maneira de criar uma linha que unisse nossas pegadas de construção à linha WaterMain?

Eu gostaria de usar o vértice do edifício mais próximo do cano principal de Água como ponto de partida.

insira a descrição da imagem aqui

— Dom
fonte

10

Se você estiver procurando uma solução que não exija o desenvolvimento de uma ferramenta .NET, use o script python abaixo para realizar exatamente o que deseja. Eu tinha exatamente a mesma necessidade e escrevi o seguinte script como solução. Configure-o como uma ferramenta ArcCatalog com os 4 parâmetros, ou comente os parâmetros e descomente as variáveis codificadas permanentemente e execute-as diretamente.

# 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")

— Jeff Berry
fonte

2

Procure o método "NearestFeature" no IIndexQuery2.

Você pode usar isso para obter o recurso principal de água mais próximo de cada prédio. Eu acho que você precisaria, de alguma forma, percorrer os vértices de cada edifício para descobrir qual é a distância mais próxima desse recurso, construir uma nova polilinha usando os vértices do edifício e o domínio da água como pontos finais. A única vez que fiz isso foi usando classes de dois pontos, desejando poder oferecer mais do que isso em cima da minha cabeça ..: 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);

— ericoneal
fonte