fui ao site da união europeia obter raster do factor k (Erodibilidade do Solo), o tiff está com um tamanho de célula 500mx500m e foi produzido com os produtos da ESRI. O sistema de coordenadas é ETRS89 Lambert Azimuthal Equal Area.
Apliquei uma reprojecção do raster para ETRS89 / Portugal TM06. com comando Warp do menu Raster > Projections.
Em seguida, seleccionei a opção do raster gravar como, e definir o tamanho de célula para 20mx20m.
O que acontece, é que ele grava e considera o tamanho de célula 20mx20m mas na realidade, se for medir o pixel ele mantem os 500 metros de tamanho de célula.
Experimentei no SAGA (resampling) e no GRASS GIS (todos os módulos do resample).
Experimentei também com um número de linhas e colunas fixo e deu no mesmo.
Será que seja algum parâmetro interno do tif?
Será que o resample não se aplica para valores que vão do sentido de um tamanho de célula maior para um menor? (faz-me sentido esta última dúvida)
Experimentei o ArcCoiso para despitar e tb acontece o mesmo.
No GRASS (plugin GRASS), usa o r.resample, mas tens que definir o tamanho do pixel na Region. O GRASS do Processing, se não me engano, usa o tamanho do pixel do raster original, porque não se conseguem definir as propriedades da Region.
O que dizem os metadados da imagem? E como estás a “medir”? É que se estás a medir “visualmente”, é possível que te pareça que o pixel ainda tem os 500m, porque apesar de estares a transformar para um tamanho 20x20, todos os 625 pixeis do novo raster 20x20 vão ter o mesmo valor que o pixel correspondente do raster original com 500x500.
No processing podes sim definir o extent (até de forma interactiva no
canvas) da region e a resolução da analise.
Por defeito ele usa o extent e resolução do input.
O r.resample da processing toolbox é (umas das) ferramentas que pode
ser usada para esta tarefa.
O problema poderá ser a projecção ETRS89 Lambert Azimuthal Equal Area
que costuma dar dores de cabeça não só no qgis. Neste casos faço uma
copia do input para um sistema de coordenadas diferente.
-- Giovanni --
2014-04-23 16:33 GMT+01:00 Pedro Venâncio <pedrongvenancio@gmail.com>:
No GRASS (plugin GRASS), usa o r.resample, mas tens que definir o tamanho do
pixel na Region. O GRASS do Processing, se não me engano, usa o tamanho do
pixel do raster original, porque não se conseguem definir as propriedades da
Region.
No processing podes sim definir o extent (até de forma interactiva no
canvas) da region e a resolução da analise.
Por defeito ele usa o extent e resolução do input.
O r.resample da processing toolbox é (umas das) ferramentas que pode
ser usada para esta tarefa.
O problema poderá ser a projecção ETRS89 Lambert Azimuthal Equal Area
que costuma dar dores de cabeça não só no qgis. Neste casos faço uma
copia do input para um sistema de coordenadas diferente.
No GRASS (plugin GRASS), usa o r.resample, mas tens que definir o tamanho do
pixel na Region. O GRASS do Processing, se não me engano, usa o tamanho do
pixel do raster original, porque não se conseguem definir as propriedades da
Region.
Eu já tinha converti para etrs89/ Portugal TM06 (alias nem me deixe ter os
dois no mesmo projecto, o qgis para e diz exceed project limits)
Para passares para PT-TM06, tens previamente de cortar o raster pelos
limites de Portugal Continental, que é para onde o sistema de referência
está definido.
Vê a resposta do Pedro. Esse erro exceed project limits está relacionado com o facto de existirem coordenadas que quando vão ser reprojectadas dão valores fora do extent válido do ETRS-89. Se adoptares a sugestão do Pedro, ou seja, cortares o raster só para o extent de Portugal, tenho a certeza que esse erro desaparece.
As opções vão dar resultados ligeiramente diferentes porque no gdalwarp usei o método de reamostragem bilinear, enquanto o r.resample usa o nearest neighbor. Para usares o nearest neighbor no gdal, usa -r near. O -20000 foi um valor que eu atribui para nodata, podes usar outro.
Não te envio os resultados porque ficaram pesadotes, mas experimenta aí que vai funcionar!
As opções vão dar resultados ligeiramente diferentes porque no gdalwarp usei o método de reamostragem bilinear, enquanto o r.resample usa o nearest neighbor. Para usares o nearest neighbor no gdal, usa -r near. O -20000 foi um valor que eu atribui para nodata, podes usar outro.
Não te envio os resultados porque ficaram pesadotes, mas experimenta aí que vai funcionar!