¿Por qué existe la diferencia entre dos tramas superpuestas?

Question

Tengo dos rásteres de el Monte Vesubio: abarca de 20 kilómetros, y el otro 50 kilómetros. Tanto los rásteres están en WGS84 sistema de coordenadas geográficas, y se superponen entre sí:

Cuando puedo usar QGIS Raster calculator para calcular la diferencia entre ellos, un nuevo archivo raster se crea que es totalmente negro - lo que demuestra que no hay absolutamente ninguna diferencia en los valores de los píxeles en la parte donde se superponen:

Sin embargo, cuando me proyecto, tanto los rásteres a Azimutal equidistante de proyección, y vuelva a comprobar la diferencia con la Trama de la calculadora, la diferencia entre estos dos se vuelve a proyectar los rásteres existe! Aquí está cómo se ve:

¿Por qué es esto así?

La diferencia no existe entre los dos rásteres cuando están en WGS84 sistema de coordenadas geográficas, pero una vez que se proyecta Azimutal equidistante de proyección, una diferencia significativa aparece?

Por qué?

He utilizado el cubic método de remuestreo, pero tengo la similar de la parte superior de diferencia con el bilinear .

Sólo por el bien de la comprobación, he intentado utilizar Transversal de Mercator proyección (en lugar de Azimutal equidistante), y de nuevo existe una diferencia entre los dos primeros proyectado rásteres. ¿Por qué es eso así?

---

Información adicional: El proyecto inicial WGS84 archivos raster, he utilizado la siguiente sintaxis:

gdalwarp -s_srs EPSG:4326 -t_srs "+proj=aeqd +lat_ts=40.81266 +lon_0=14.414252" -r cubic -of GTiff "C:/vesuvius_radius_20KM.tif" "C:/vesuvius_radius_20KM_cubic_aeqd.tif"

Para proyectarlas Transversal de Mercator:

gdalwarp -s_srs EPSG:4326 -t_srs EPSG:32633 -r cubic -of GTiff "C:/vesuvius_radius_50KM.tif" "C:/vesuvius_radius_50KM_cubic_tm.tif"

Aquí están las iniciales de archivos raster (WGS84):

https://www.dropbox.com/s/98hcftjsnmrqs8p/vesuvius_radius_20KM.tif?dl=0 https://www.dropbox.com/s/a68921e6tpszt0f/vesuvius_radius_50KM.tif?dl=0

Aquí están sus Azimutal equidistante proyectado rásteres:

https://www.dropbox.com/s/4mexe0rmieam9ri/vesuvius_radius_20KM_cubic_aeqd.tif?dl=0 https://www.dropbox.com/s/0hcle0b151kbod1/vesuvius_radius_50KM_cubic_aeqd.tif?dl=0

Y la diferencia entre dos Azimutal equidistante proyectado rásteres:

https://www.dropbox.com/s/2w39hgb9kcrzb8h/vesuvius_radius_20-50KM_cubic_aeqd.tif?dl=0

Answer 1

Si usted reproyectar una trama a otra proyección, el tamaño de celda no coincide con el original. Por esta razón, el nuevo valor de la celda será calculado con respecto a las células vecinas, como un valor medio ponderado.

En los archivos originales, ambos tamaños de celda corresponde exactamente (0.000833333 grados), y los valores de la celda eran idénticos. Pero después de reproyección, tamaño de la celda (82.7608 m vs 82.7644 m) y los valores no son idénticos ya.

---

ACTUALIZACIÓN

Para aclarar, yo he creado el vector de las rejillas de su trama archivos con la misma extensión y tamaño de la celda:

La cuadrícula verde es del original de 50 kilómetros - y a 20 km raster, perfectamente superpuestas con un tamaño de celda de 0.000833333° (3 segundos). Las líneas rojas son de los 50km-aeqd de la cuadrícula, y el azul de los 20km-aeqd-grid.

Para cada uno de los rojo y azul de la célula, los valores de los tocados células verdes serán tomados para calcular el nuevo valor de la celda. Desde el rojo y el azul de la célula no se alinean, se obtienen diferentes valores de celda para ambas redes.

Después de reproyección, los rásteres no son exactamente rectangular más. GDAL intenta crear un mapa de bits de tamaño casi pegada al original. Para obtener un número de filas y columnas de la nueva medida, se obtiene de las diferentes espaciamiento de los 50 km y 20 km de ceñir. Puede forzar el tamaño de la celda con un valor determinado con -tr opción en gdalwarp, pero las rejillas todavía puede ser desplazado el uno contra el otro. Usted tendría que definir el destino de las extensiones con -te así que asegúrese de que tanto los rásteres de alinear.

Así que algo como esto debería funcionar:

gdalwarp -overwrite -te -28000 4500000 15000 4540000 -tr 100 100 -s_srs EPSG:4326 -t_srs "+proj=aeqd +lat_0=0 +lon_0=14.414252 +x_0=0 +y_0=0 +ellps=WGS84 +units=m +no_defs" -of GTiff vesuvius_radius_20KM.tif 20km-neu.tif
gdalwarp -overwrite -te -60000 4470000 50000 4570000 -tr 100 100 -s_srs EPSG:4326 -t_srs "+proj=aeqd +lat_0=0 +lon_0=14.414252 +x_0=0 +y_0=0 +ellps=WGS84 +units=m +no_defs" -of GTiff vesuvius_radius_50KM.tif 50km-neu.tif