[GRASS-user] Reprojection from LLC to UTM (various zones)

[Maciej Sieczka]

Miriam,

Ok, if I got you right you have a volume (3D raster) projected in lcc,
imported into a GRASS location which has it's projection set up
according to this volume. Now you want to visualise your volume over a
DEM. But your DEM is in latlong "projection". If this is right, you
should propably be doing like this:

1. Import DEM into a location that has the same projection as DEM
(propably, it's latlong on WGS84 ellipsoid; the EPSG code 4326); let's
call the location "source" and let's assume you are working in mapset
"miriam".

2. In "source" location:

g.region rast=DEM

3. In the "target" lcc location, the mapset your volume is in, do:

g.region rast3d=the_volume_name
r.proj in=DEM out=DEM_lcc mapset=miriam location=source

This will re-project a subset of the whole DEM which covers your
volume; if you want a bigger piece of DEM, enlarge region before
r.proj. Read r.proj and g.region manual.

Hope it works.

Maciek