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)?
Many thanks in advance
Tomas
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