Tomas Lanczos wrote: > I have contour maps converted to rasters (rock temperatures) in > different depths (-500, -1000, -1500, -2000, -2500, -3000, -4000, -6000 > meters below sea level). What is the best way to interpolate them to one > 3d raster (volume)?
for evenly spaced raster slices first set the 3D region, e.g. g.region t=-500 b=-6000 tbres=500 then use the r.to.rast3 module to make the 3D raster volume map. but you would need to interpolate intermediate levels for missing slices. For -3500 you could use r.series method=average, for 4500,5000,5500 you would have to use r.mapcalc manually. v.vol.rst could be used to interpolate from 3D points, but RST+contour lines usually is prone to artifact problems do to the uneven spatial distribution of data points. see also the 3D help page on the wiki site: http://grass.gdf-hannover.de/wiki/Help_with_3D Hamish ____________________________________________________________________________________ Get easy, one-click access to your favorites. Make Yahoo! your homepage. http://www.yahoo.com/r/hs _______________________________________________ grassuser mailing list grassuser@grass.itc.it http://grass.itc.it/mailman/listinfo/grassuser