¿Función de ArcGIS para restar la celda de la siguiente en un raster?

Question

Usando ArcGIS 10.5.1 - Tengo un raster con ciertos valores (elevación en este caso). Me gustaría generar un ráster que se base en la sustracción de cada celda de su celda vecina (la que está directamente a la izquierda). ¿Existe alguna función que pueda hacer esto? Ejemplo:

Raster de entrada:

• 6 16 14
• 6 10 14
• 8 10 14

Trama de salida deseada:

• Sin datos 10 -2
• Sin datos 4 4
• Sin datos 2 4

Raster usado --> Float. Aquí está mi archivo de cabecera

Answer 1

Sólo necesita el bloc de notas para editar el archivo de cabecera (*.hdr) de la trama flotante.

De:

ncols         348
nrows         355
xllcorner     1758048.7257138
yllcorner     5905009.5863752
cellsize      1
NODATA_value  -9999
byteorder     LSBFIRST

Para:

ncols         348
nrows         355
xllcorner     1758049.7257138
yllcorner     5905009.5863752
cellsize      1
NODATA_value  -9999
byteorder     LSBFIRST

es decir, la tercera fila desde arriba.

Flujo de trabajo para una resolución de trama de 1m*1m utilizando herramientas:

arcpy.Warp_management(in_raster="dem", source_control_points="'0 0'", target_control_points="'1 0'", out_raster="../SCRATCH/shifted", transformation_type="POLYORDER0", resampling_type="NEAREST")
arcpy.gp.RasterCalculator_sa(""""dem" - "shifted"""", "../SCRATCH/diff")

Entrada y salida:

Answer 2

Si no quieres usar python, esto se puede hacer con el model builder, desplazando la rejilla a la izquierda en 1 pixel, hacer los cálculos y luego volver a desplazarla. Al hacer esto se introduce una columna de nódulos a la derecha por lo que el paso final es utilizar el clip original para recortarlo de nuevo al tamaño original.

El modelo sería:

Lo importante es establecer la configuración del entorno de la herramienta menos para que sea la unión de las extensiones de la cuadrícula original y la desplazada, esto permite crear la columna de no datos.