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

[Maciek Sieczka]

Michael

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 
luck.

Maciek

----- 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
>
>