[QGIS-pt] Converter tamanho de célula

[Pedro Venâncio]

João, fiz aqui uns testes e funcionou das duas formas. Faz assim:

1 - Descarrega a CAOP e faz um dissolve all, para ficares com os limites do
continente em EPSG 3763;
2 - Transforma essa shape para EPSG 3035;

Opção a)
3.1 - Corta o raster da erodibilidade com a shape em EPSG 3035
3.2 - Transforma o raster resultante para EPSG 3763
3.3 - Corre o r.resample do Processing com a region cellsize 20.

Opção b)
3 - Corre este comando:
gdalwarp -dstnodata -20000 -q -cutline Cont_AAD_CAOP2013_pt_3035.shp
-crop_to_cutline -s_srs EPSG:3035 -t_srs "+init=pt:pttm06 +wktext" -tr 20
20 -r bilinear K_new_crop.tif K_new_crop_pt_3_3763_20x20_gdal_2.tif

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!


Abraço!
Pedro
-------------- próxima parte ----------
Um anexo em HTML foi limpo...
URL: <http://lists.osgeo.org/pipermail/qgis-pt/attachments/20140424/a13dedcb/attachment.html>