[QGIS-pt] Converter tamanho de célula
João Gaspar
joao.f.r.gaspar gmail.com
Quinta-Feira, 24 de Abril de 2014 - 00:52:53 PDT
Pedro, já consegui.
A situação era mesmo do sistema de coordenadas.
1ºConverti os ficheiros vectoriais para o sistema de coordenadas do raster.
2º Fiz as operações de corte.
3º Finalizei com o warp gdal para o ETRS89 Portugal TM06 com tamanho de
célula 20 e já tenho o tamanho certo 20x20.
Obrigado a todos.
Abraço
No dia 24 de Abril de 2014 às 00:11, Pedro Venâncio <
pedrongvenancio gmail.com> escreveu:
> 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
>
> _______________________________________________
> QGIS-pt mailing list
> QGIS-pt lists.osgeo.org
> http://lists.osgeo.org/cgi-bin/mailman/listinfo/qgis-pt
>
>
-------------- próxima parte ----------
Um anexo em HTML foi limpo...
URL: <http://lists.osgeo.org/pipermail/qgis-pt/attachments/20140424/b001ee3e/attachment.html>
More information about the QGIS-pt
mailing list