Preguntas etiquetadas con numpy

Soy nuevo en SIG. Tengo un código que convierte las imágenes infrarrojas de Marte en mapas de inercia térmica, que luego se almacenan como matrices numpy 2D. He guardado estos mapas como archivos hdf5, pero realmente me gustaría guardarlos como imágenes ráster para poder procesarlos en QGIS. He...

He estado tratando de verificar mis filtros en el ráster DEM para el reconocimiento de patrones y siempre resulta en la falta de las últimas filas y columnas (como ... 20) . He intentado con la biblioteca PIL, carga de imágenes. Luego con numpy. La salida es la misma. Pensé, algo estaba mal con...

Estoy tratando de volver a proyectar / volver a muestrear con los enlaces de Python GDAL, pero obtengo resultados ligeramente diferentes en comparación con los de la utilidad de línea de comandos gdalwarp. Vea la actualización a continuación para ver un ejemplo más corto Este script ilustra el...

Estoy interesado en aprender a utilizar matrices NumPy para optimizar el geoprocesamiento. Gran parte de mi trabajo implica "big data", donde el geoprocesamiento a menudo lleva días para realizar ciertas tareas. No hace falta decir que estoy muy interesado en optimizar estas rutinas. ArcGIS 10.1...

He comenzado a aprender a manipular los datos LAS en python y quería ver cómo otros manejan los archivos LAS. Me gustaría leer los puntos (estoy usando una matriz numpy), y filtrar las clases 1 y 2 (sin clasificar y suelo) a una matriz separada. Tengo el siguiente código, pero parece que no puedo...

Para obtener un solo valor de píxel en un punto de un ráster con rasterio, hay un ejemplo aquí: https://github.com/mapbox/rasterio/pull/275 Sin embargo, ¿hay una API directa dentro de rasterio (y no el cli) que se pueda usar para extraer valor en un solo punto en un ráster? - EDITAR with...

Estoy usando ArcGIS 10.1 y quiero crear un nuevo ráster basado en dos rásteres preexistentes. El RasterToNumPyArray tiene un buen ejemplo que quiero adaptarse. import arcpy import numpy myArray = arcpy.RasterToNumPyArray('C:/data/inRaster') myArraySum = myArray.sum(1) myArraySum.shape =...

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 datos ráster en una tabla de postgres que quiero ingresar a python como una matriz numpy. Estoy usando psycopg2 para conectarme a la base de datos. Puedo descargar los datos pero vuelven como una cadena (probablemente un binario serializado). ¿Alguien sabe cómo tomar esta cadena y...

Estoy usando Python 2.6.5 (32 bits) con Numpy 1.3 y Gdal 1.9.1 instalado en Windows 7 de 64 bits. Estoy tratando de leer un dataset ráster Imagine (.img) de 800 MB en una matriz Numpy para hacer algo de álgebra ráster, pero tan pronto como ejecuto el siguiente código, Python.exe se bloquea. from...

Estoy intentando extraer valores de corte naturales de un ráster usando Python. El proceso previsto va: Objeto de trama Arcpy a NumPy Array (usando RasterToNumPyArray ) a valores de ruptura natural (usando la función PySAL Natural Breaks ). import arcpy, pysal from pysal.esda.mapclassify import...

Estoy tratando de aprender las cuerdas del procesamiento de imágenes de detección remota utilizando enlaces Python GDAL y numpy. Como primer intento, estoy leyendo un archivo geotiff Landsat8, hago una simple manipulación y escribo el resultado en un nuevo archivo. El siguiente código parece...

¿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...

Me gustaría reclasificar un archivo ráster de un ráster con 10 clases a un ráster con 8 clases usando pyhton, gdal y / o numpy. Las clases se representan como enteros. He intentado seguir los pasos de esta publicación Reclasificar rásteres usando GDAL y Python , el numpy.equal doc y también...

Estoy usando el método GDAL ReadAsArray para trabajar con datos ráster usando numpy (específicamente reclasificación). Como mis rásteres son grandes, proceso las matrices en bloques, iterando en cada bloque y procesando con un método similar al ejemplo de GeoExamples . Ahora estoy buscando la...

¿Alguien ha tenido éxito haciendo que los pandas trabajen con la instalación de Python de QGIS en Windows? He usado el instalador OSGeo4W para instalar QGIS2.2, que viene con Python 2.7 y Numpy 1.5.1. Cualquier versión razonablemente nueva de pandas requiere Numpy 1.6.1 o posterior. Así que ese...

Hoy actualicé a ArcGIS Desktop 10.2.2, así como a la herramienta de procesamiento en segundo plano de 64 bits, y noté que tanto Numpy como Matplotlib se instalaron con Python. De hecho, en la sección " ¿Qué hay de nuevo? " Para 10.2.1, dice: ArcGIS 10.2.1 se ha actualizado para incluir Python...

Tengo un problema con la biblioteca GDAL cuando abro un gran archivo GeoTiff, con un tamaño de aproximadamente 32000x32000. No puedo usar la ReadAsArrayfunción porque el tamaño máximo de la matriz numpy en python. Pero me pregunto por qué QGIS puede abrir ese archivo fácilmente. ¿Cuál es la técnica...

El problema consiste en modelar la propagación de una señal (por ejemplo, luz o sonido, etc.) a través de una serie de obstáculos, como en la figura a continuación. La señal no puede pasar a través de la superficie inferior (terreno), pero puede atravesar obstáculos. Quiero contar la cantidad de...