Especificar valor de célula de la trama

Question

Estoy buscando automatizar un procedimiento que va a cambiar un único valor de la celda en el ráster. Tengo un conocido x, y y z, y un ráster. Necesito cambiar el valor de la celda (el más cercano a mi x e y) el valor especificado para cada valor de z.

He estado haciendo esto el uso de Saga (cambiar los valores de las celdas [interactivo] módulo), pero necesito una más método automatizado para su uso en un archivo por lotes o secuencia de comandos de python. Pensé que puede ser un módulo para la actualización de una celda ráster con base en el punto (vector) los valores? Nota, no quiero a reclasificar los valores de cuadrícula, en lugar quiero especificar uno.

Tengo acceso a QGiS, Saga y ArcGIS básico (sin espacial o paquetes 3D). Si alguien puede me apunte en una dirección inicial que sería muy apreciada. Gracias,

Answer 1

Puesto que usted está buscando un método automatizado utilizando Python, recomiendo la GDAL bibliotecas. Si estás en windows recomiendo el gisinternals repositorio de los archivos binarios.

Supongo que usted no está familiarizado con el GDAL las bibliotecas y cómo hacer ráster de procesamiento con ellos (de ahí su pregunta), así que aquí está un tutorial. El ráster de procesamiento, de inicio en el tutorial 4, y por el tiempo que usted ha leído a través de tutorial 6 tendrá el conocimiento.

Esencialmente, las librerías GDAL, en conjunto con NumPy leer la trama en una matriz que a continuación, puede manipular como se requiere. Una vez que la trama es una matriz, puede cambiar los valores en las celdas individuales, el trabajo en bloques de celdas, aplicar filtros, etc. Los tutoriales le ayudará a empezar pero para algunas grandes intensos procesos también recomiendo el uso de la ndimage módulos de SciPy - pero para su propósito la NumPy los métodos descritos en el tutorial va a ser suficiente (cue cuando el uso de SciPy es cuando la trama procesamiento se tediosamente lento el uso de NumPy).