Preguntas etiquetadas con shapely

Estoy buscando algunas sugerencias sobre cómo hacer que mi código de Python sea más eficiente. Normalmente, la eficiencia no me importa, pero ahora estoy trabajando con un archivo de texto de ubicaciones de EE. UU. Con más de 1,5 millones de puntos. Con la configuración dada, se tarda unos 5...

¿Hay una manera fácil de transformar objetos Shapely (a saber, polígonos y polígonos múltiples) de una proyección a otra sin tener que cavar y extraer coordenadas a mano? De hecho, ni siquiera me importa si son objetos Shapely en este punto, solo quiero pasar características y una proyección, y...

Estoy tratando de ejecutar Shapely en una nueva máquina, e intenté instalarlo desde la fuente, usar el instalador en su página web y usar pip install shapely. Cada referencia que intento hacer desde la biblioteca, excepto from shapely import *arroja el error WindowsError: [Error 126] The specified...

¿Alguien puede demostrar una manera simple de escribir estructuras de datos de geometría de bien formados a archivos de formas? Estoy particularmente interesado en polígonos con agujeros y cadenas lineales. También sería beneficioso mantenerse alejado de arcpy (por lo que osgeo, pyshp, etc., todo...

Estoy buscando calcular estadísticas focales para cada celda de un ráster, dentro de un vecindario de un criterio específico. Antecedentes: tengo tres rásteres binarios, cada uno de los cuales representa un único tipo de vegetación de interés. Me gustaría calcular el porcentaje de cobertura de...

Antecedentes Desde un punto conocido, requiero establecer el "perímetro visible" circundante más cercano contra una tabla de MultiLineStrings, como se muestra en el diagrama. He buscado en este sitio con una serie de términos (por ejemplo, borde mínimo, perímetro mínimo, vecino más cercano, clip,...

Dado que geopandas GeoDataFramecontiene una serie de polígonos, me gustaría obtener el área en km² de cada entidad en mi lista. Este es un problema bastante común, y la solución habitual sugerida en el pasado ha sido usar shapelyy pyprojdirectamente (por ejemplo, aquí y aquí ). ¿Hay alguna manera...

Estoy tratando de capturar todos los polígonos no superpuestos que se indican a continuación usando Shapely (dados los polígonos A, B y C). Además, espero hacerlo sin iteraciones, pruebas de intersección, etc. La respuesta aceptada a esta pregunta expresa el método PostGIS, pero parece que 'unión'...

He estado tratando de crear una función que básicamente haga lo mismo que la función de "disolución" de QGIS. Pensé que sería súper fácil pero al parecer no. Entonces, por lo que he reunido, el uso de fiona con bien proporcionado debería ser la mejor opción aquí. Empecé a jugar con archivos...

Estoy tratando de ajustar dos líneas entre sí usando Shapely / Geopandas, pero el resultado del ajuste es muy extraño. Lo intenté : import geopandas as gpd from shapely.geometry import * from shapely.ops import snap lines1 = gpd.GeoDataFrame.from_file('lines1.shp') lines1 = lines1.to_crs({'init':...

Tengo polígonos de coordenadas en (Python bien formados) que se ve así POLYGON ((24.8085317 46.8512821, 24.7986952 46.8574619, 24.8088238 46.8664741, 24.8155239 46.8576335, 24.8085317 46.8512821)) Me gustaría calcular el área de este polígono en km ^ 2. ¿Cuál sería la mejor manera de hacer esto...

parece una pregunta tonta, sin embargo, no puedo encontrar una respuesta clara al respecto: ¿qué unidades geopandas / bien proporcionadas usan para calcular la distancia / área entre

He intentado varios ejemplos de código utilizando bibliotecas como shapefile, fiona y ogr para intentar verificar si un punto (x, y) se encuentra dentro de los límites de un multipolígono creado con ArcMap (y, por lo tanto, en formato shapefile). Sin embargo, ninguno de los ejemplos funciona bien...

Estoy usando Shapely en python y me dan un MultiLineStringcon un montón de Linestringobjetos. Puedo garantizar que todos los LineStringobjetos son líneas simples con solo 2 vértices y que todos forman parte de una sola línea (sin ramas). Quiero "conectar los puntos" y crear uno solo LineString....

Estoy buscando una solución de python de código abierto para convertir el ráster a polígono (sin ArcPy). Conocía la función GDAL para convertir el ráster a polígono, y aquí está el manual: http://pcjericks.github.io/py-gdalogr-cookbook/raster_layers.html#polygonize-a-raster-band Sin embargo,...

Tengo problemas para comprender el uso de índices espaciales con RTree. Ejemplo: tengo 300 puntos almacenados en búfer, y necesito saber el área de intersección de cada búfer con un archivo de forma poligonal. El archivo de forma poligonal tiene> 20,000 polígonos. Se sugirió que usara índices...

Tengo un geopandas que GeoDataFrame contiene cientos de formas Polygony MultiPolygongeometrías. Los polígonos se superponen en muchos lugares. Me gustaría hacer una nueva geometría que contenga un recuento de cuántas de ellas se superponen. Algo como esto: ¿Alguien tiene alguna idea sobre cómo...

Cuando intento instalar Shapely en mi computadora con Windows de 64 bits, no puedo hacer que la biblioteca GEOS funcione. Hasta ahora, he ejecutado el instalador OSGeo4W desde el que instalé GDAL (creo que la biblioteca geos está incluida en ese paquete). Después de eso, verifiqué y tengo...

Estoy leyendo un archivo de forma usando Fiona: with fiona.open('data/boroughs/boroughs_n.shp') as source: mpolys = [p for p in source] candidate = polys[0]['geometry'] esto me da un dict que contiene una lista de listas de coordenadas y escribe 'MultiPolygon'. Sin embargo, no estoy seguro de...

[EDITAR: la solución a esto fue simplemente usar OGR para leer archivos de forma. Ver el ejemplo de geographika.] En un archivo de forma ESRI, no hay distinción entre Polígonos y MultiPolygons. Además, no existe una distinción explícita entre los agujeros interiores y los anillos exteriores...