<div dir="ltr"><div>João, fiz aqui uns testes e funcionou das duas formas. Faz assim:<br><br></div><div>1 - Descarrega a CAOP e faz um dissolve all, para ficares com os limites do continente em EPSG 3763;<br></div><div>2 - Transforma essa shape para EPSG 3035;<br>
</div><div><br>Opção a) <br>3.1 - Corta o raster da erodibilidade com a shape em EPSG 3035</div><div>3.2 - Transforma o raster resultante para EPSG 3763<br></div><div>3.3 - Corre o r.resample do Processing com a region cellsize 20.<br>
<br></div><div>Opção b)<br></div><div>3 - Corre este comando:<br>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<br>
<br></div><div>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.<br>
</div><div><br></div><div>Não te envio os resultados porque ficaram pesadotes, mas experimenta aí que vai funcionar!<br><br><br></div><div>Abraço!<br></div><div>Pedro<br></div></div>