[GRASSLIST:7034] Re: region related issues with r.proj

Maciek Sieczka werchowyna at epf.pl
Mon Jun 6 16:21:18 EDT 2005


Not sure how to cope with r.proj here, but you could:
1 cut the piece od DEM you want with r.mapcalc
2 r.out.gdal to GeoTiff
3 gdalwarp
4 r.in.gdal

A note: as far I was told on this list before by wise folks and as I have 
confirmed in my experience, reprojecting raster DEM is *not* a good idea for 
most typical purpose (e.g. you won't be happy with the resulting aspect or 
slope). Better t.to.vect, v.proj, re-interpolate (say v.surf.rst). Good 


----- Original Message ----- 
From: "Michael Kowen" <michael.kowen at case.edu>
To: <grasslist at baylor.edu>
Sent: Monday, June 06, 2005 9:19 PM
Subject: [GRASSLIST:7031] region related issues with r.proj

> Hello again,
> I have just stumbled upon a bit of an issue while using r.proj in
> grass6.0.  The situation is as follows:  I have imported a 10m dem of a
> comparatively large scale (zoomed out area) of a 10m dem, and what I
> would like to do is import this dem into a location of a smaller region
> using r.proj (I have done this successfully in a different scenario),
> but it complains that "the input map is outside current region."       I
> figure that the problem is just as grass reports it, but I am not sure
> how to fix it/make it work.  The projection converion is from nad83 to
> nad27. (On a separate note, is this inverse chronical conversion an
> error--is it better to run map projects in the most recent projection
> systems?)
> thanks,
> Michael Kowen

More information about the grass-user mailing list