¿Dónde se introduce este desplazamiento ráster en la caja de herramientas de procesamiento de QGIS?

8

Estoy usando la caja de herramientas de procesamiento (¡increíble!) En QGIS 2.0 para reclasificar algunos valores NoData en un ráster de -9999 a 0, pero en algún punto del proceso se introduce un desplazamiento de aproximadamente 0.2 grados en la capa de salida.

Estoy usando la herramienta "Reclasificar (simple)" (creo que es parte de SAGA), con algunos parámetros bastante básicos:

reclasificar parámetros

. . . y donde el ráster de entrada se situó correctamente en comparación con esta cuadrícula hexadecimal:

ráster de entrada

. . . el ráster de salida es golpeado hacia el sureste por dos celdas hexadecimales, o aproximadamente 0.2 grados

ingrese la descripción de la imagen aquí

El origen de entrada es -180, 83.623. . . pero el origen de salida es -179.05,83.05. Normalmente supondría que esto es el resultado de una conversión WGS84 defectuosa, pero la entrada, la salida y el proyecto son todos EPSG: 4326.

¿De dónde viene este desplazamiento de origen, y vale la pena presentar un boleto SAGA o es algo que pasé por alto?

Actualización El registro del proceso indica que se realiza una transformación en el origen de entrada:

Transformation x' = -179.999989 + x * 0.100000 + y * 0.000000
Transformation y' = 83.623125 + x * 0.000000 + y * -0.100000

y unas pocas líneas más tarde aparecen los nuevos parámetros, ahora incorrectos:

Grid system: 0.1; 3600x 1735y; -179.949989x -89.826875y
qgis raster geoprocessing qgis-2 saga Bill Morris
fuente
1
debe probar la herramienta SAGA nativa utilizando la GUI o la CLI y ver si obtiene el mismo resultado.
Giovanni Manghi
1
Demasiados puntos muertos de error de estos documentos de instalación en mi máquina: no estoy interesado en borrar mi configuración actual de GDAL / PROJ solo para probar SAGA de forma nativa.
Bill Morris
1
Si está utilizando SAGA en QGIS, SAGA ya está instalado en su máquina.
Giovanni Manghi el
Después de más de 2 horas de leer detenidamente los documentos, no tengo la menor idea de CÓMO LANZAR LA GUI EN OSX. La documentación recuerda a R en todos los aspectos negativos. Ninguna variación en "saga" (-gis, -gui, etc.) produce nada en la línea de comando, y la compilación de la fuente devuelve errores completamente inútiles. Gracias a los cielos misericordiosos (desarrolladores centrales, es decir) que este desastre ahora es accesible para los humanos en una barra lateral de QGIS. Espero sinceramente que SAGA no sea lo que está cambiando mi origen de salida de trama, porque nunca podré depurarlo.
Bill Morris
¿Eres capaz de lanzar saga con saga_guicomando en su lugar?
SaultDon

Respuestas:

1

Su desplazamiento es de 1/2 píxel (0.1 grados).

Longitud:

-179.999989 - -179.949989 = 0.05

La latitud es menos obvia porque se hace referencia a una cuadrícula desde la parte superior izquierda, la otra desde la parte inferior izquierda:

83.623125 - (1735 * 0.1) - -89.826875 = 0.05

Supongo que esto se relaciona con el problema RasterPixelIsPoint o RasterPixelIsArea:

http://www.remotesensing.org/geotiff/faq.html#PixelIsPoint

¿Estás usando gdal 1.7 o 1.8?

Marc Pfister
fuente
En todo el sistema tengo GDAL 1.10.1; por alguna razón no puedo hacer que la consola QGIS python informe la versión utilizada allí
Bill Morris
se ve mucho más de 1/2 píxel en la imagen.
radouxju
Lo es, pero el tamaño de píxel también está desactivado. Aquí está el gdalinfo para
Bill Morris
Ah, al ver el gdalinfo parece que el cambio es de .95 y -0.5731265. El desplazamiento de 1/2 píxel anterior es probablemente una conversión correcta de área a punto. Su ráster original tiene algunos límites extraños, tal vez están arruinando la conversión.
Marc Pfister
Buen punto. Es extraño que cuelgue un poco del borde del mundo. El ráster es el resultado de una conversión de vector; déjame ver si puedo rehacer esa conversión con los límites correctos.
Bill Morris