Riproiezioni di dati raster
Da Wikipedia, l'enciclopedia libera.
Torna a HOWTOs e procedure tipo
Riproiezioni di dati raster
L'utile programma r.proj di GRASS permette di riproiettare mappe da una location ad un altra. Sfortunatamente, non sempre l'uso è agevole in quanto le region di destinazione deve essere compatibile con la region di partenza. Quando si proietta da lat/long a una proiezione piana o viceversa, tirare a indovinare non è proprio cosa comoda e r.proj può restituire un fastidioso errore 'Input map is outside current region'. Un utile metodo è quello di creare un vector delle dimensioni della region del raster nella location di partenza, riproiettarlo usando v.proj nella location destinazione - il programma vettoriale non ha problema di contenimento della region - e adattare quindi la region al vector risultante.
Nella location 1:
g.region -p rast=myimage v.mkgrid map=myimage_extent grid=1,1 position=region
Nella location 2:
v.proj location=location1 mapset=PERMANENT in=myimage_extent g.region vect=myimage_extent r.proj location=location1 mapset=PERMANENT in=myimage res=risoluzione
A questo punto il raster è perfettamente riproiettato. Si noti che in generale è necessario specificare esplicitamente la risoluzione nella region destinazione delle location, stando ben attenti a definirla confrontabile con la risoluzione dei dati di partenza, nella nuova proiezione.
--FrancescoLovergine 23:29, Apr 2, 2007 (CEST)
Per approfondimenti sul tema "georeferenziare dati raster" si rimanda alla pagina Georeferenziare raster
