<div dir="ltr"><div><div>Pedro, já consegui.<br><br>A situação era mesmo do sistema de coordenadas.<br><br>1ºConverti os ficheiros vectoriais para o sistema de coordenadas do raster.<br><br>2º Fiz as operações de corte.<br>
<br>3º Finalizei com o warp gdal para o ETRS89 Portugal TM06 com tamanho de célula 20 e já tenho o tamanho certo 20x20.<br><br></div>Obrigado a todos.<br></div>Abraço<br><div class="gmail_extra"><br><br><div class="gmail_quote">
No dia 24 de Abril de 2014 às 00:11, Pedro Venâncio <span dir="ltr"><<a href="mailto:pedrongvenancio@gmail.com" target="_blank">pedrongvenancio@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<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!<span class="HOEnZb"><font color="#888888"><br></font></span></div><span class="HOEnZb"><font color="#888888"><div>
Pedro<br></div></font></span></div>
<br>_______________________________________________<br>
QGIS-pt mailing list<br>
<a href="mailto:QGIS-pt@lists.osgeo.org">QGIS-pt@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/cgi-bin/mailman/listinfo/qgis-pt" target="_blank">http://lists.osgeo.org/cgi-bin/mailman/listinfo/qgis-pt</a><br>
<br></blockquote></div><br></div></div>