Convertir polilínea con ID izquierda e ID derecha en polígonos con ID

Question

Tengo un archivo de forma polilínea que contiene información de ID izquierda y ID derecha y luchando en la conversión de esta capa polilínea a polígono con ID a través arcpy.

Estoy utilizando la licencia de ArcGIS(editor) por lo que herramienta de conversión de característica en polígono no está disponible para mí. También sé que hay varios métodos disponibles como este y éste en ArcGIS si tienes licencia de editor pero quiero usar arcpy para automatizar mi flujo de trabajo.

¿Alguna idea de cómo puedo hacerlo usando arcpy?

Answer 1

No veo una manera strightforward de hacer esto en ArcPy licencia normal sin reinventar las ruedas..

Acaba de leer sobre Shapely (descargar enlace) un paquete python, que incluye el módulo "polygonize" función de suites JTS/GEOS familiares..

Mi idea es que podemos importar el paquete Shapely en ArcPy y obtener el WKT de las geometrías de polilíneas ESRI y pasarlas al poligonizador. Esto dará los polígonos resultantes. Estos podrían ser convertidos de nuevo a ESRI geometrías polilínea. Aunque debo decir, mi conocimiento de Python es terrible esto es sólo una idea genérica

aquí hay una muestra de polygonizer..

>>> from shapely.geometry import LineString, Point
    >>> from shapely.ops import polygonize
    >>> lines = [
    ... LineString(((0, 0), (1, 1))),
    ... LineString(((0, 0), (0, 1))),
    ... LineString(((0, 1), (1, 1))),
    ... LineString(((1, 1), (1, 0))),
    ... LineString(((1, 0), (0, 0))),
    ... LineString(((5, 5), (6, 6))),
    ... Point(0, 0),
    ... ]
    >>> result = polygonize(lines)
    >>> list(result)
    [<shapely.geometry.polygon.Polygon object at ...>, <shapely.geometry.polygon.Polygon object at ...>]

    >>> lines2 = [
    ... ((0, 0), (1, 1)),
    ... ((0, 0), (0, 1)),
    ... ((0, 1), (1, 1)),
    ... ((1, 1), (1, 0)),
    ... ((1, 0), (0, 0)),
    ... ((5, 5), (6, 6)),
    ... ]
    >>> result2 = polygonize(lines2)
    >>> list(result2)
    [<shapely.geometry.polygon.Polygon object at ...>, <shapely.geometry.polygon.Polygon object at ...>]

Tampoco estoy seguro de lo que podemos hacer con LEFT_ID y RIGHT_ID. Creo que son irrelevantes. Póngame al día si me equivoco.

Answer 2

Si quieres usar Python y el GeoProcessor para construir polígonos a partir de tus líneas tienes que echar un vistazo al archivo de Ayuda y buscar una página llamada trabajar con geometría en Python .

Tendrás que crear un cursor sobre tus polilíneas y encadenarlas, puede que primero tengas que "explotarlas" en puntos De todas formas no estoy seguro de cómo puedes resolver el tema del orden en que visitas las polilíneas. No sé (ya que nunca lo he probado en Python) si la geometría del polígono resultante se puede simplificar, es decir, el orden de los vértices se hace correcto cuando se crea el polígono.

Answer 3

Quizás podrías usar arcpy.FeatureToRaster_conversion() para crear un raster a partir de tus polilíneas y luego arcpy.RasterToPolygon_conversion() para crear el polígono