Estou usando a (incrível!) Caixa de ferramentas de processamento no QGIS 2.0 para reclassificar alguns valores NoData em uma varredura de -9999 a 0, mas em algum momento do processo há um deslocamento de cerca de 0,2 graus na camada de saída.
Estou usando a ferramenta "Reclassify (Simple)" (acredito que faz parte do SAGA), com alguns parâmetros bastante básicos:
. . . e onde a varredura de entrada estava situada corretamente em comparação com essa grade hexadecimal:
. . . a varredura de saída é colidida a sudeste por duas células hexadecimais, ou cerca de 0,2 graus:
A origem da entrada é -180, 83.623. . . mas a origem da saída é -179.05,83.05. Normalmente, eu diria que este é o resultado de uma conversão WGS84 com defeito, mas a entrada, a saída e o projeto são todos EPSG: 4326.
De onde vem essa compensação de origem e vale a pena registrar um ingresso da SAGA ou é apenas algo que eu esqueci?
Atualização O log do processo indica que uma transformação foi realizada na origem da entrada:
Transformation x' = -179.999989 + x * 0.100000 + y * 0.000000
Transformation y' = 83.623125 + x * 0.000000 + y * -0.100000
e, algumas linhas depois, os novos parâmetros aparecem agora incorretos:
Grid system: 0.1; 3600x 1735y; -179.949989x -89.826875y
saga_gui
comando?Respostas:
Seu deslocamento é de 1/2 pixel (0,1 graus).
Longitude:
-179.999989 - -179.949989 = 0,05
A latitude é menos óbvia porque uma grade é referenciada no canto superior esquerdo e a outra no canto inferior esquerdo:
83,623125 - (1735 * 0,1) - -89,826875 = 0,05
Suponho que isso esteja relacionado ao problema RasterPixelIsPoint ou RasterPixelIsArea:
http://www.remotesensing.org/geotiff/faq.html#PixelIsPoint
Você está usando gdal 1.7 ou 1.8?
fonte