Hi Micha,
i fixed my previous problem...i can see the map on the display (it was because
of the wrong resolution setting: in fact i had to add to the r.proj command the
parameter res=<resolution of the raster converted into the map units of the
target location>)
I can see my raster map, but the raster changed shape: all the printed squares
of the utm system became rectangles....mmm
i guess it has got smthing to do with the resolution setting...
----Messaggio originale----
Da: micha@arava.co.il
Data: 11/03/2010 14.47
A: "ciccpaol@libero.it"<ciccpaol@libero.it>
Cc: <grass-user@lists.osgeo.org>
Ogg: Re: R: Re: [GRASS-user] Input raster map is outside current region
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
On 03/11/2010 04:49 PM, ciccpaol@libero.it wrote:
Hi Micha,
i fixed my previous problem...i can see the map on the display (it was because
of the wrong resolution setting: in fact i had to add to the r.proj command the
parameter res=<resolution of the raster converted into the map units of the
target location>)
ׂGood. You probably want to set the res even before r.proj, with g.region res=xxx
I can see my raster map, but the raster changed shape: all the printed squares
of the utm system became rectangles....mmm
That's what reprojection does!
i guess it has got smthing to do with the resolution setting...
----Messaggio originale----
Da: micha@arava.co.il
Data: 11/03/2010 14.47
A: "ciccpaol@libero.it"<ciccpaol@libero.it>
Cc:<grass-user@lists.osgeo.org>
Ogg: Re: R: Re: [GRASS-user] Input raster map is outside current region
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
This mail was received via Mail-SeCure System.
--
Micha Silver
http://www.surfaces.co.il/
Arava Development Co. +972-52-3665918