[GRASS-user] Help: Converting a raster map between locations (from wgs84 lat / long to UTM)

[Nikos Alexandris]

On Wed, 2008-10-22 at 15:24 -0200, Daniel Victoria wrote:
> It's because your region in the target location is wrong.
> 
> Here is a little trick:
> 
> In your lat/lon location, after you set your region, run the command
> v.in.region (check correct syntax). This will create a bounding box of
> the region.
> 
> Then, go to your target location (UTM) and use the command v.proj to
> import the vector you just created.
> 
> Now, set your region to the imported vector with g.region
> vect=name_of_vector_here
> 
> Finally, run r.proj!!
> 
> The only problem at which I don't know any exact method is setting the
> target resolution correctly... I know it's done with g.region res=XX
> but, the conversion from lat/long to meters is not a straight one...
> 
> Cheers
> Daniel
> 
> On Wed, Oct 22, 2008 at 2:44 PM, Corrado <ct529 at york.ac.uk> wrote:
> > Dear friends,
> >
> > I am trying to convert a GTOPO30 DEM (raster) map (for UK) from lat / long wgs
> > 84 to UTM for UK.
> >
> > I am running into rserious problems. I am using the procedure as described in
> > the book open source gis a grass gis approach"
> >
> > 1) in the WGS84 lat long location (ukll) I resize the map using d.zoom
> > 2) I save the new region as well as the new sub-raster (for scotland)
> > 3) I open grass in the UTM location
> > 4) I try to import using
> >
> > r.proj in=gtopo30_cs out=gtopo30_cs_utm location=ukll mapset=PERMANENT
> > method=cubic --o
> >
> > All I can get is: ERROR: Input map is outside current region
> >
> > What am I doing wrong? I have tried so many different alternative solution.
> >
> > Best,
> > --
> > Corrado Topi

OK, maybe your way is better Daniel :-)

Greets, Nikos