¿Cómo puedo obtener las coordenadas proyectadas, así como los valores de píxeles reales en esas coordenadas de un archivo GeoTiff y luego guardarlos en una matriz numpy? Tengo el archivo arsenci020l.tif, y sus coordenadas están en metros. A continuación se muestra la salida resumida de gdalinfo que...
Cerrado. Esta pregunta está fuera de tema . Actualmente no está aceptando respuestas. ¿Quieres mejorar esta pregunta? Actualice la pregunta para que sea sobre el tema para el intercambio de pila de sistemas de información geográfica. Cerrado el mes pasado ....
Me gustaría crear un script Python en QGIS que tomará algunos datos (shp + tif, en mi ejemplo), usará Tamplate Map Composer (desde un archivo) y exportará la composición creada a una imagen png. Casi sin experiencia en programación (menos que el conocimiento básico de Python) busco en Google...
Novato aquí luchando con SIG. Estoy tratando de mapear los barrios de la ciudad de Milwuakee usando archivos de formas que se encuentran en el sitio web de su condado . Estoy siguiendo el hilo aquí con cierto éxito. Mi código es da: from pyproj import Proj, transform # wisconsing EPSG:32054 #...
Ocasionalmente, pero con demasiada frecuencia desafortunadamente, mis servicios de AGS simplemente se detienen sin ningún motivo. Por lo general, no me doy cuenta de esto hasta que los usuarios me llaman y me preguntan por qué sus mapas no tienen datos. ¿Hay alguna forma (estaba pensando en...
Utilizo la siguiente línea para informar al usuario sobre el estado: iface.mainWindow().statusBar().showMessage("Status:" + str(i)) El complemento tarda unos 2 minutos en ejecutarse en mi conjunto de datos, pero Windows lo detecta como "no responde" y deja de mostrar las actualizaciones de...
Quiero crear una función que seleccione una función y la amplíe (similar en QGIS). Por lo tanto, existe la siguiente función: QgsMapLayerRegistry.instance().addMapLayer(self.vlayer) def zoomTo(self): layer = self.vlayer atable = self.ui.table selectList=[] for i in
Tengo un conjunto de líneas que representan líneas de autobús. Algunas de las líneas se superponen y toman los mismos caminos. Soy capaz de extraer los nodos. Sin embargo, estoy interesado en extraer solo cruces como este: ¿Cómo puedo hacer esto? Estoy buscando formas con QGIS o...
¿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...
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...
Tengo un geodataframe con puntos y algunos datos asociados. Quiero trazarlo en un mapa usando geopandas y hacer que el tamaño de los puntos corresponda a una de las columnas del geodataframe. Hasta ahora tengo el siguiente código: base = world.plot(color='white',
Acabo de comenzar a usar pycharm con qgis pero no pude conectarlos a ambos. Pycharm siempre permanece en un estado de "espera de conexión". La mayoría de los tutoriales disponibles apuntan a las ventanas, pero estoy usando ubuntu, por lo que no pude encontrar la forma de depurar el código qgis en...
Estoy tratando de construir un complemento para cargar un compositor de impresión desde un archivo, generar un atlas y exportarlo a la imagen. Hasta ahora he tenido éxito al cargar la plantilla y exportarla a la imagen. No he podido agregar ninguna de las capas de la leyenda (que también están en...
Me encargaron crear un modelo usando ArcGIS ModelBuilder, que inicialmente ejecuta un par de tareas de geoprocesamiento, antes de unir una tabla de Excel con el resultado de esos geoprocesos. Sin embargo, el modelo necesita abrir realmente el archivo Excel primero, para que el usuario lo edite...
Quiero crear un conjunto de datos vectoriales de puntos de 10000 puntos (o más grandes) dentro de un ráster binario, donde los puntos deben estar restringidos a áreas donde el valor del ráster es 1. Intenté los siguientes pasos. Poligonalizar ráster QGIS: Vector -> Herramientas de...
Estoy escribiendo una aplicación web de uso intensivo de datos que se entrega a través de apache. Mi pregunta es sobre cómo organizar mejor el procesamiento dado que existen múltiples opciones. Tengo a mi disposición OpenLayers / JQuery / Javascript, PostGIS / Postgresql (con pgsql), python /...
Estoy buscando registros duplicados en archivos dbf basados en el atributo llamado 'ID'. Tengo varios archivos dbf de 500,000 registros a 1.5 millones y sé que hay una gran cantidad de duplicados. Me gustaría agregar un campo 'Duplicar' que dice Sí o No (o 1 o 0 está bien) cuando el atributo ID...
He estado observando un rendimiento inusual con un script de geoprocesamiento de Python. El script (adjunto) realiza las siguientes acciones: Use un cursor de búsqueda para buscar la zona UTM correspondiente a las características del polígono Crear objeto de referencia espacial basado en los...
¿Existe una documentación oficial para el módulo Python GDAL? Hay varios libros de cocina, pero poca información sobre cada función o argumentos para las funciones más