[GRASSLIST:555] raster 3d from DEMs

GRASS GRASS user
1

hello

I created a couple DEMs from drillhole data, DEMs represent a geological boundary (e.g. Triassic/Jurassic), now i want to fill space between this boundary and create a raster 3d. Is it possible to do, i suppose - yes, but how to do it. Can anyone help me.

regards
maciej

Maciej Tomaszczyk
Polish Geological Institute
Geological Mapping Department

2

I created a couple DEMs from drillhole data, DEMs represent a
geological boundary (e.g. Triassic/Jurassic), now i want to fill
space between this boundary and create a raster 3d. Is it possible to
do, i suppose - yes, but how to do it.

Is v.vol.rst what you are looking for?
  http://grass.ibiblio.org/grass61/manuals/html61_user/v.vol.rst.html

see also
  http://grass.ibiblio.org/grass61/manuals/html61_user/raster3dintro.html
  http://grass.ibiblio.org/grass61/manuals/html61_user/raster3D.html

Hamish

3

hi
i need each geological unit to be one value.
v.vol.rst generate in the space beetwen geological bounduary many values, (the slice has many colors)
now i try to use v.to.rast3, maybe this is what i looking for.
thanks for help
maciej

Maciej Tomaszczyk
Polish Geological Institute
Geological Mapping Department

-----Wiadomość oryginalna-----
Od: Hamish [mailto:hamish_nospam@yahoo.com]
Wysłano: Wt 2006-04-04 07:45
Do: Maciej Tomaszczyk
DW: GRASSLIST@baylor.edu
Temat: Re: [GRASSLIST:555] raster 3d from DEMs

I created a couple DEMs from drillhole data, DEMs represent a
geological boundary (e.g. Triassic/Jurassic), now i want to fill
space between this boundary and create a raster 3d. Is it possible to
do, i suppose - yes, but how to do it.

Is v.vol.rst what you are looking for?
  http://grass.ibiblio.org/grass61/manuals/html61_user/v.vol.rst.html

see also
  http://grass.ibiblio.org/grass61/manuals/html61_user/raster3dintro.html
  http://grass.ibiblio.org/grass61/manuals/html61_user/raster3D.html

Hamish

4

Hi Maciej,
to define geological units (or in general a 1D/2D/3D categorical variable in the "space") you have to use algorithms able to compute them.
We generally use stochastic algorithms based on Gaussian Simulations as implemented in gslib and gstat or Uncert packages after a study concerning the variable spatial distribution.

I hope this helps

m.

On 4 Apr 2006, at 10:57 AM, Maciej Tomaszczyk wrote:

hi
i need each geological unit to be one value.
v.vol.rst generate in the space beetwen geological bounduary many values, (the slice has many colors)
now i try to use v.to.rast3, maybe this is what i looking for.
thanks for help
maciej

Maciej Tomaszczyk
Polish Geological Institute
Geological Mapping Department

-----Wiadomość oryginalna-----
Od: Hamish [mailto:hamish_nospam@yahoo.com]
Wysłano: Wt 2006-04-04 07:45
Do: Maciej Tomaszczyk
DW: GRASSLIST@baylor.edu
Temat: Re: [GRASSLIST:555] raster 3d from DEMs

> I created a couple DEMs from drillhole data, DEMs represent a
> geological boundary (e.g. Triassic/Jurassic), now i want to fill
> space between this boundary and create a raster 3d. Is it possible to
> do, i suppose - yes, but how to do it.

Is v.vol.rst what you are looking for?
http://grass.ibiblio.org/grass61/manuals/html61_user/v.vol.rst.html

see also
http://grass.ibiblio.org/grass61/manuals/html61_user/raster3dintro.html
http://grass.ibiblio.org/grass61/manuals/html61_user/raster3D.html

Hamish

*******************************************************
Marino Vetuschi Zuccolini
zucco@dipteris.unige.it
PhD / Geochemist
Laboratory of Geochemistry
http://www.dipteris.unige.it/geochimica

Visit our BAXEICO computing cluster homepage
http://qed.dipteris.unige.it/ganglia -- out of service

DIPartimento per lo studio della TErra e delle sue RISorse - Università di Genova
Tel. ++39 010 3538136 Fax. ++39 010 352169
Corso Europa 26, 16132 - Genova - Italy

5

Hi,

--- Ursprüngliche Nachricht ---
Von: "Maciej Tomaszczyk" <mtom@pgi.gov.pl>
An: <GRASSLIST@baylor.edu>
Betreff: [GRASSLIST:555] raster 3d from DEMs
Datum: Mon, 3 Apr 2006 12:25:05 +0200

hello

I created a couple DEMs from drillhole data, DEMs represent a geological
boundary (e.g. Triassic/Jurassic), now i want to fill space between this
boundary and create a raster 3d. Is it possible to do, i suppose - yes,

but

how to do it. Can anyone help me.

GRASS can only handle raster 3d data with one top and bottom elevation
value.
But if you want to visualize the geological boundaries filled with data, you
can
use a combination of r3.mapcalc, r3.out.vtk and Paraview.

example:
Create for each geological boundary a g3d map with at least 2 depths,

1. g.region t=100 b=0 tbres=50
2. r3.mapcalc "triassic = 1"
3 r3.mapcalc "jurassic = 2"

the use r3.out.vtk to export the data with top and bottom surface

r3.out.vtk -s input=triassic output=/tmp/triassic.vtk bottom=dem1 top=dem2

r3.out.vtk -s input=jurassic output=/tmp/jurassic.vtk bottom=dem2 top=dem3

Now visualize this data with Paraview
available at www.paraview.org

If you have some point data within this geological boundaries, you can
interpolate them within a high resolution g3d map (depth > 40) with
v.vol.rast
and make slices of a specific elevation with r3.cross.rast.
The slices created with r3.cross.rast can be converted to g3d maps with
r.to.rast3.

example:
High res 3d region
g.region b="the lowest dem" t="the highest dem" tbres="a useful resolution"

Now use v.vol.rast to generate the g3d volume data, output is
"generated_volume"

extract the data with r3.cross.rast
r3.cross.rast input=generated_volume elevation=dem1 output=datadem1
r3.cross.rast input=generated_volume elevation=dem2 output=datadem2

Set a 3d region with only two depths
g.region t=100 b=0 tbres=50
r.to.rast3 input=dem1,dem2 output=triassic_data

Now export and visualize this data with Paraview
r3.out.vtk -s input=triassic_data output=/tmp/triassic_data.vtk bottom=dem1
top=dem2
Paraview --data=/tmp/triassic_data.vtk

Maybe this helps a bit
Best
Soeren

regards
maciej

Maciej Tomaszczyk
Polish Geological Institute
Geological Mapping Department

--
Analog-/ISDN-Nutzer sparen mit GMX SmartSurfer bis zu 70%!
Kostenlos downloaden: http://www.gmx.net/de/go/smartsurfer