Onde esse deslocamento de varredura está sendo introduzido na caixa de ferramentas de processamento QGIS?

8

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:

reclassificar parâmetros

. . . e onde a varredura de entrada estava situada corretamente em comparação com essa grade hexadecimal:

raster de entrada

. . . a varredura de saída é colidida a sudeste por duas células hexadecimais, ou cerca de 0,2 graus:

insira a descrição da imagem aqui

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
Bill Morris
fonte
1
você deve testar a ferramenta SAGA nativa usando a GUI ou a CLI e veja se obtém o mesmo resultado.
Giovanni Manghi
1
Muitos becos sem saída com erros desses documentos de instalação na minha máquina - não estou interessado em limpar minha configuração atual do GDAL / PROJ apenas para testar o SAGA nativamente.
Bill Morris
1
Se você estiver usando SAGA no QGIS, o SAGA já está instalado em sua máquina.
Giovanni Manghi
Depois de mais de duas horas examinando os documentos, não tenho a menor idéia de COMO LANÇAR A GUI NO OSX. A documentação é uma reminiscência de R de todas as maneiras ruins. Nenhuma variação na "saga" (-gis, -gui, etc.) produz qualquer coisa na linha de comando e a compilação a partir da fonte retorna erros completamente inúteis. Agradeça aos céus misericordiosos (ou seja, os principais desenvolvedores) que essa bagunça agora está acessível aos humanos em uma barra lateral do QGIS. Espero sinceramente que o SAGA não seja o que está mudando minha origem de saída raster, porque nunca poderei depurá-lo.
Bill Morris
Você é capaz de iniciar a saga com saga_guicomando?
precisa saber é o seguinte

Respostas:

1

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?

Marc Pfister
fonte
Em todo o sistema, tenho o GDAL 1.10.1; por qualquer motivo eu não posso começar o console QGIS python para relatar a versão utilizada lá
Bill Morris
parece muito mais que 1/2 pixel na imagem.
precisa saber é o seguinte
É, mas o tamanho do pixel também está desativado. Aqui está o gdalinfo para ambas as tiffs de entrada e saída
Bill Morris
Ah, vendo o gdalinfo, parece que a mudança é de 0,95 e -0,5731265. O deslocamento de 1/2 pixel acima provavelmente é uma área correta para conversão de pontos. Sua varredura original tem alguns limites estranhos, talvez eles estejam estragando a conversão.
Marc Pfister
Bom ponto. É estranho que esteja pendurado na extremidade do mundo em um grau. A varredura é o resultado de uma conversão do vetor; deixe-me ver se consigo refazer essa conversão com limites corretos.
Bill Morris