Estoy tratando de dibujar en un mapa mapbox, basado en los resultados de suncalc y usando geojson. Primero intenté crear 2 funciones, una para cada línea que intentaba dibujar. Pero cuando hice eso, solo mostraría la última función llamada. Así que concluí que no sé cómo administrar las capas, ya...
yo tengo GeoPandas 0.2.1 py27_0 fiona 1.7.0 np110py27_1 Python 2.7.10 instalado en el entorno Anaconda 2-4.1.1-Windows-x86. Puedo construir GeoDataFrameleyendo conjuntos de datos de entrada y operar los datos, pero guardar un conjunto de datos de salida no conserva el sistema de...
En la documentación de gdal_calc se indica Calculadora ráster de línea de comando con sintaxis numpy . Más adelante hay algunos ejemplos donde en uno de ellos: gdal_calc.py -A input.tif --outfile = result.tif --calc = "A * (A> 0)" --NoDataValue = 0 - significa establecer valores de cero e...
Tengo un shapefile con áreas de la ciudad, pero desafortunadamente el shapefile solo contiene líneas simples. El archivo de forma se ve así: Quiero convertir esas líneas en un archivo de forma de polígono de área en QGIS. Lo probé con la herramienta "Líneas a polígonos" (Vector> Herramientas...
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...
en un intento de importar datos editables a mapbox, no logro importar los datos de una forma que el sistema acepte. Rechaza los datos con un mensaje de error que indica: cities.geojson: error de entrada. no se recomienda el miembro crs de estilo antiguo en la línea 1. Al leer eso, supuse que...
Me pregunto si hay alguna explicación para los parámetros que se pueden establecer al ejecutar el Algoritmo "Join Attributes by Location"en QGIS 2.18. He echado un vistazo aquí: Documentación Pruebas QGIS | Unir atributos por ubicación Documentación para QGIS 2.8 | Complemento fTools Pero no...
Tengo una tabla con 8 columnas y ~ 16.7 millones de registros. Necesito ejecutar un conjunto de ecuaciones if-else en las columnas. He escrito un script usando el módulo UpdateCursor, pero después de algunos millones de registros se queda sin memoria. Me preguntaba si hay una mejor manera de...
Cerrada . Esta pregunta necesita estar más centrada . Actualmente no está aceptando respuestas. ¿Quieres mejorar esta pregunta? Actualice la pregunta para que se centre en un problema solo editando esta publicación . Cerrado hace 7 meses . Tengo una...
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 tratando de hacer coincidir segmentos pequeños con un segmento más grande con el que probablemente estén más relacionados: relativamente cerca, rumbo similar y uno frente al otro. Aquí hay un ejemplo típico de los datos que tengo: Aquí tendría que hacer coincidir el segmento 652 con...
Tengo un gran conjunto de redes lineales y me gustaría encontrar los dos extremos de cada red que están más distantes entre sí a lo largo de la red (en la imagen a continuación, sería D a K). La solución de la fuerza bruta para este problema es calcular la ruta más corta a lo largo de la red para...
Instalé GDAL y verifiqué en la línea de comandos, funciona y cuando ejecuto gdal2tiles.py en la línea de comandos funciona, pero cuando ejecuto el mismo archivo gdal2tiles.py con Python IDLE, obtengo "Ningún módulo llamado 'osgeo '"error. ¿Puedes
Tengo varias capas en el Panel de capas dentro de un grupo, pero no quiero que se muestren (hay muchas, por lo que se ve desordenado). Estas capas también están unidas a varias otras. ¿Es posible hacer que estas capas sean invisibles para el usuario ( es decir, no se muestran en el Panel de capas,...
¿Es posible tener un control deslizante de transparencia que se muestra de forma predeterminada en la ventana Capa, es decir, directamente después de agregar una nueva capa, en lugar de tener que habilitarlo para cada capa manualmente a través de las propiedades de
Tengo un problema al calcular los centroides de ciertos polígonos en QGIS (a través de las herramientas de geometría): para ciertos polígonos, con una forma extraña y bastante larga, los centroides calculados están fuera del polígono. Mire, por ejemplo, el polígono marcado en rojo, en la imagen a...
He creado un complemento para alimentar una base de datos Spatialite con una máscara de entrada. Hay dos campos que representan la dimensión mínima / máxima de un objeto. En el momento en que llamo a la herramienta Medir ( iface.actionMeasure().trigger()) con los dos botones "medir", haga clic en...
Tengo una capa vectorial con millones de polígonos que hacen una cobertura continua. Necesito clasificarlos según su forma. Ya estoy usando varios índices de formas de la ecología del paisaje como la compacidad ( 4piA / P ^ 2 ), el ancho medio ( 2A / P ), el número de forma ( P / sqrt (A) ),...
He estado usando un proceso llamado "Crear retícula de polígono" para crear una cuadrícula dentro de un polígono con atributos de fila y columna asignados. Es una herramienta SAGA que he estado ejecutando en QGIS. Genera una cuadrícula basada en la extensión del polígono de entrada con la celda que...
Me gustaría crear un rectángulo circundante con una relación máxima de ancho / alto ("cuadro de límite oblicuo") a partir de un conjunto de características de entrada como se muestra en la imagen a continuación: Por ejemplo, las funciones PostGIS ST_Envelopey los ST_Box2DTransformadores FME...