Me gustaría obtener la intersección de múltiples polígonos. Usando el shapelypaquete de Python , puedo encontrar la intersección de dos polígonos usando la intersectionfunción. ¿Existe una función eficiente similar para obtener la intersección de múltiples polígonos? Aquí hay un fragmento de...
Estoy tratando de encontrar dónde varias líneas se cruzan con un polígono para dos geodataframes diferentes: from shapely.geometry import Polygon, LineString import geopandas as gpd polygon = Polygon([(0, 0), (1, 0), (1, 1), (0, 1), (0, 0)]) line1 = LineString([(0.5, 0.5), (0.7, 0.7)]) line2 =...
Tengo un marco de datos geopandas de un grupo de cadenas de líneas que tienen algunos datos asociados con cada vértice / punto: Point_x = (Lat, Lon, Time, ID, Data1, Data2, Data3) Los puntos se convierten en cadenas de líneas según la ID y se ordenan por tiempo. Quiero romper las cadenas...
Estoy haciendo un flujo de trabajo a pequeña escala en el cuaderno de ipython con geopandas y bien proporcionado, obteniendo un montón de datos geoespaciales, a veces de archivos de forma, a veces de postgis (donde se realiza un procesamiento más costoso). Ahora, luego tiro las tablas postgis a...
Si tiene dos puntos, desde los cuales quiero crear un LineStringobjeto recto : from shapely.geometry import Point, LineString A = Point(0,0) B = Point(1,1) El manual bien proporcionado paraLineString estados: Una secuencia de Pointinstancias no es un parámetro constructor válido. A...
Estoy haciendo un cálculo muy simple de la longitud de una polilínea usando bien: from shapely.geometry import LineString ... xy_list = [map(float,e) for e in xy_intm] line = LineString(xy_list) s = '%s,%s,%s' % (fr,to,line.length) Mis coordenadas están en WGS84. Parece que no puedo encontrar...
Estoy tratando de crear un rectángulo girado en coordenadas lat / lon. Básicamente, tengo las coordenadas del centro y el ángulo de rotación alrededor del centro. ¿Es posible hacer uso, bien proporcionado o gdal, o algo
Para obtener las coordenadas exteriores, necesito convertir un MultiPolygon bien formado en un Polígono. Lo hago así: if poly.geometry.type == 'Polygon': x, y = poly.geometry.exterior.xy elif poly.geometry.type == 'MultiPolygon': allparts = [p.buffer(0) for p in poly.geometry] poly.geometry =...
Estoy tratando de generar polígonos para franjas orbitales satelitales. Hasta ahora tengo un método para generar dos líneas que representan el borde de cada franja en [lat, long]. Algunas de las hileras cruzan la fecha internacional y así se envuelven: Pude resolver esto con ogr2ogr...
Deseo recortar un conjunto de polilíneas (líneas negras en la imagen a continuación) en el límite exterior de un polígono. Cualquier vacío dentro del polígono debe ser ignorado. Mi salida ideal son las líneas amarillas discontinuas. Las líneas iniciales pueden o no ser rectas. La imagen es un...
Sé que este problema se ha abordado anteriormente con diferentes softwares, pero ¿cómo lo hago con Shapely? Tengo un polígono con un espacio muy pequeño en el interior de la forma: Aquí hay un primer plano: El polígono es válido, así que no puedo usar el truco del búfer, y simplificar el...
Tengo un código que estoy usando para determinar qué Polígono Shapely / MultiPolygons se cruzan con una serie de Shapely LineStrings. A través de las respuestas a esta pregunta, el código ha pasado de esto: import fiona from shapely.geometry import LineString, Polygon, MultiPolygon, shape #...
Tengo dos polígonos POLÍGONO ((- 5 -5, -5 5,5 5,5 -5, -5 -5)) POLÍGONO ((1 1, 1 7, 7 7, 7 1, 1 1)) y deseo dividirlos de manera uniforme. El resultado debe ser una división a lo largo de la línea roja marcada. Estoy usando qgis / python bien y deseo hacerlo usando
Estoy ejecutando un análisis específico donde usamos bien para crear buffers alrededor de puntos (ubicaciones de tiendas) y luego verifico si otros puntos (ubicaciones de usuarios) están presentes dentro de ese valor de buffer. Cuando compruebo la distancia usando bien, resulta ser diferente de la...
Estoy trabajando con Python, Shapely y Fiona ... ¿Hay alguna forma en que pueda dar una entrada de un archivo de forma de línea y un archivo de forma de polígono, y obtener un resultado de los puntos de intersección y sus coordenadas? La ilustración proporciona una descripción precisa...
Tengo un marco de datos de coordenadas X, Y que representan puntos a lo largo de los caminos tomados por varias entidades diferentes. Pseudo-datos aquí, pero es más o menos de la forma: entity_id lat lon time 1001 34.5 14.2 4:55 pm 1001 34.7 14.5 4:58 pm 1001 35.0 14.6 5.03 pm 1002 27.1 19.2...
Estoy tratando de obtener un archivo KML ( wijken.kml ) en Shapely . El archivo KML se valida con el esquema XML adecuado, por lo que supongo que la entrada es correcta. Rutas que he probado: 1) Convertir a formato WKT o WKB, y leer con funciones incorporadas Conversión: ogr2ogr -f CSV...
¿Cómo dibujo una elipse con