Me gustaría saber cómo obtener todos los valores ráster dentro de un polígono usando gdal o pygeoprocessing, sin leer toda la cuadrícula como una matriz. pygeoprocessing y gdal pueden hacer estadísticas zonales, pero solo la función min, max, mean, stdev o count están disponibles desde dicha...
Lo sé, mi pregunta es similar a algunas anteriores en este sitio. Tengo muchos archivos CSV (coordenadas geográficas) para importar a qgis (y luego convertirlos), y la forma habitual no es la mejor manera de hacerlo (demasiado tiempo). Tengo casi 500 archivos CSV (coordenadas wgs84) y esto es lo...
Vi una gran interfaz para FME con Python ¿Qué están haciendo con eso? Quiero
Me gustaría instalar QGIS. Estoy en un sistema que por defecto es Python 2.5x y descargué python 3.7x con homebrew. Creé un entorno virtual con este código. brew install pyenv brew install pyenv-virtualenv exec $SHELL eval "$(pyenv init -)" eval "$(pyenv virtualenv-init -)" pyenv activate...
Estoy tratando de usar ArcObject en Python. Tengo problemas para asignar una interfaz (es decir, IZ.InterpolateZsBetween) a mi objeto con datos de geometría, el código es: from Comtypes.client import CreateObject, GetModule import arcgisscripting gp =
Quiero convertir una tabla PostGIS en un archivo de forma (sin usar pgsql2shp). Para crear una geometría en el archivo de forma, tengo que dar Xmin, Ymin y Xmax, Ymax, y la geometría que tengo en mi tabla PostGIS es de forma irregular (puedo obtener el exterior usando el cuadro delimitador, pero...
Quiero hacer una imagen multiespectral desde cero para hacer algunas pruebas al respecto. Algo realmente simple como 5 bandas completamente uniformes con ruido de sal y pimienta o un cuadrado de valores diferentes en el centro. Claramente, esto sería solo una pila de matrices, una matriz...
Me preguntaba si alguien tiene alguna experiencia en obtener datos de elevación de un ráster sin usar ArcGIS , sino obtener la información como un pitón listo dict? Recibo mis datos XY como una lista de tuplas: xy34 =[perp_obj[j].CalcPnts(float(i.dist), orientation) for j in range...
¿Es posible detectar mensajes de error de salida de PostgreSQL en OGR después de que una consulta haya fallado? Por ejemplo conn.ExecuteSQL('SELECT * FROM non_existing_table;') En Postgres esto devuelve un mensaje de error: "Error: relation 'non_existing_table' does not exist." En ogr...
Al iterar sobre una capa vectorial usando el siguiente código (resumido del ejemplo en el libro de cocina pyqgis ), ¿hay alguna manera de verificar si una función está seleccionada? provider = vlayer.dataProvider() feat = QgsFeature() allAttrs =
Tengo un problema / pregunta con respecto a python y QGIS. ¿Es posible obtener los datos de una característica seleccionada en QGIS usando python, y si es así, cómo? En mi ejemplo, tengo una capa que muestra algunas parcelas para la planificación urbana. Si ahora selecciono una de las parcelas,...
Estoy tratando de hacer una figura que muestre datos azimutales con un rango diferente de incertidumbres en cada punto. Esta figura de la vieja escuela de un artículo de 1991 captura la idea de la "trama de pajarita" a la que apunto: ¿Alguna sugerencia sobre cómo podría hacer una figura similar?...
¿Es posible configurar el CRS del proyecto de un proyecto en QGIS mediante un comando de Python? En mi caso, el usuario debe abrir QGIS, iniciar mi complemento escrito y este complemento cargará datos espaciales que tienen el código epsg 31467 (Gauss Krueger Zone 3). Lamentablemente, el CRS...
Tengo una útil herramienta de script que recorre un espacio de trabajo y renombra y copia archivos de forma en un conjunto de datos de características. Sin embargo, si hay un archivo de forma dañado en algún lugar del espacio de trabajo, el script falla y detiene el procesamiento. ¿Cómo manejas...
¿Hay algún medio para verificar si 2 capas ráster dadas tienen contenido idéntico ? Tenemos un problema con nuestro volumen de almacenamiento compartido corporativo: ahora es tan grande que lleva más de 3 días realizar una copia de seguridad completa. La investigación preliminar revela que uno de...
Estoy tratando de entender cómo puedo usar el módulo csv en python para abrir un archivo csv en la misma carpeta que el script python, y luego crear un shapefile usando el módulo shapefile pyshp. El archivo csv se ve así, pero puede tener un par de miles de filas de
Me gustaría ejecutar un script de Python a través de qgis --code myscript.pyy luego salir inmediatamente. Estoy usando iface.actionExit().trigger(), que mata QGIS cuando lo ejecuto desde la consola de Python, pero no cuando pongo el script al que paso --code. ¿Cuál es la forma correcta de salir...
¿Es posible cambiar el color al agregar una capa vectorial? Tengo una capa vectorial con áreas de agua y quiero cargarla con un estilo de color azul. La capa se agrega con: QgsMapLayerRegistry.instance().addMapLayer(self.vlayer) ¿Puedo cambiar el color de alguna manera después o durante la...
He encontrado datos de elevación de alta precisión (resolución de 1/3 a 1/9 de segundo de arco) del conjunto de datos de elevación nacional proporcionado por el USGS. Viene en formatos IMG, Arcgrid y GridFloat. Sé que ArcGrid corresponde a un paquete de software pago, pero estoy tratando de...
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...