Hi ,
thanks Micha for your help! I tried the steps you suggested me, but they didnt
work properly; after trying that, i found an article about that
http://wiki.gfoss.it/index.php/Riproiezioni_di_dati_raster
that is more or less the same you suggested me.
Working those steps, I finally got the right region extent. In fact when I'm
inside the target location and I'm displaying the reprojected raster map, the
coordinates values shown on the bottom of the page are right, but I can't see
the map! The page is all blank.
What you do here is:
- go back to the source location and run
v.in.region out=src_region
This creates a vector polygon of the bounds of the current region.
- now switch to the new, target region and do
v.proj location=<src_location> mapset=<src_mapset> in=src_region
out=tar_region
to reproject that rectangular polygon to the target region . Next run
g.region vect=tar_region
And the region in your target location should match where the raster
wants to be.
- Now try your r.proj comand
--
Micha Silver
http://www.surfaces.co.il/
Arava Development Co. +972-52-3665918
On 03/11/2010 03:34 PM, ciccpaol@libero.it wrote:
Hi ,
thanks Micha for your help! I tried the steps you suggested me, but they didnt
work properly; after trying that, i found an article about that
http://wiki.gfoss.it/index.php/Riproiezioni_di_dati_raster
that is more or less the same you suggested me.
Working those steps, I finally got the right region extent. In fact when I'm
inside the target location and I'm displaying the reprojected raster map, the
coordinates values shown on the bottom of the page are right, but I can't see
the map! The page is all blank.
I guess it's time to ask for more information: Can you post the outputs for *both locations* of:
(Source location)
g.region -p
v.info -g <src vector region as below>
r.info -g <original raster>
(Target location)
g.region -p
v.info -g <projected vector region>
r.info -g <final projected raster which doesn't appear>
What you do here is:
- go back to the source location and run
v.in.region out=src_region
This creates a vector polygon of the bounds of the current region.
- now switch to the new, target region and do
v.proj location=<src_location> mapset=<src_mapset> in=src_region
out=tar_region
to reproject that rectangular polygon to the target region . Next run
g.region vect=tar_region
And the region in your target location should match where the raster
wants to be.
- Now try your r.proj comand
--
Micha Silver
http://www.surfaces.co.il/
Arava Development Co. +972-52-3665918
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
This mail was received via Mail-SeCure System.
--
Micha Silver
http://www.surfaces.co.il/
Arava Development Co. +972-52-3665918