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 _______________________________________________ grassuser mailing list [email protected] http://grass.itc.it/mailman/listinfo/grassuser
