[GRASS-user] Interpolation with faultlines?

I remember that some time ago, there was some discussion on
interpolation using faultlines. Any advances in that subject?

cheers

Carlos

--
+-----------------------------------------------------------+
              Carlos Henrique Grohmann - Guano
  Geologist M.Sc - Doctorate Student at IGc-USP - Brazil
Linux User #89721 - carlos dot grohmann at gmail dot com
+-----------------------------------------------------------+
_________________
"Good morning, doctors. I have taken the liberty of removing Windows
95 from my hard drive."
--The winning entry in a "What were HAL's first words" contest judged
by 2001: A SPACE ODYSSEY creator Arthur C. Clarke

Can't stop the signal.

It's the same question I was going to ask. It would be an important
advance in Grass, as I'm obliged to use commercial sw to do it (Esri)
both to interpolate with barriers and to do 3D TINs with constraints.

Giovanni

2008/2/13, Carlos Guâno Grohmann <carlos.grohmann@gmail.com>:

I remember that some time ago, there was some discussion on
interpolation using faultlines. Any advances in that subject?

cheers

Carlos

--
+-----------------------------------------------------------+
              Carlos Henrique Grohmann - Guano
  Geologist M.Sc - Doctorate Student at IGc-USP - Brazil
Linux User #89721 - carlos dot grohmann at gmail dot com
+-----------------------------------------------------------+
_________________
"Good morning, doctors. I have taken the liberty of removing Windows
95 from my hard drive."
--The winning entry in a "What were HAL's first words" contest judged
by 2001: A SPACE ODYSSEY creator Arthur C. Clarke

Can't stop the signal.
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

On 13/02/08 18:40, G. Allegri wrote:

It's the same question I was going to ask. It would be an important
advance in Grass, as I'm obliged to use commercial sw to do it (Esri)
both to interpolate with barriers and to do 3D TINs with constraints.

The answer was (and as far as I know still is) to either run v.surf.rst several times, once for each part within fault lines (using mask) or to use r.surf.nnbathy (AddOns).

Moritz

Moritz Lennert pisze:

On 13/02/08 18:40, G. Allegri wrote:

It's the same question I was going to ask. It would be an important
advance in Grass, as I'm obliged to use commercial sw to do it (Esri)
both to interpolate with barriers and to do 3D TINs with constraints.

The answer was (and as far as I know still is) to either run v.surf.rst several times, once for each part within fault lines (using mask) or to use r.surf.nnbathy (AddOns).

r.surf.nnbathy does not support fault lines either.

FWIW, there is surfit [1] (FOSS). It supports fault lines AIUI - linear features which don't have a value assigned, but only depict a curve at which the interpolation should break. They call break/fault lines "inequalities". I still haven't used surfit, but from examples [2] it seems interesting.

[1]http://surfit.sourceforge.net
[2]http://surfit.sourceforge.net/surfit/examples.html

Maciek

The attached picture is a sketch I made based on what I remember from
a conference, long time ago, about some software that did interpolate
with faultlines and how it worked.

Carlos

On Thu, Feb 21, 2008 at 4:42 PM, Maciej Sieczka <tutey@o2.pl> wrote:

Moritz Lennert pisze:

> On 13/02/08 18:40, G. Allegri wrote:
>> It's the same question I was going to ask. It would be an important
>> advance in Grass, as I'm obliged to use commercial sw to do it (Esri)
>> both to interpolate with barriers and to do 3D TINs with constraints.

> The answer was (and as far as I know still is) to either run v.surf.rst
> several times, once for each part within fault lines (using mask) or to
> use r.surf.nnbathy (AddOns).

r.surf.nnbathy does not support fault lines either.

FWIW, there is surfit [1] (FOSS). It supports fault lines AIUI - linear
features which don't have a value assigned, but only depict a curve at
which the interpolation should break. They call break/fault lines
"inequalities". I still haven't used surfit, but from examples [2] it
seems interesting.

[1]http://surfit.sourceforge.net
[2]http://surfit.sourceforge.net/surfit/examples.html

Maciek

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

--
+-----------------------------------------------------------+
Carlos Henrique Grohmann - Guano
Geologist M.Sc - Doctorate Student at IGc-USP - Brazil
Linux User #89721 - carlos dot grohmann at gmail dot com
+-----------------------------------------------------------+
_________________
"Good morning, doctors. I have taken the liberty of removing Windows
95 from my hard drive."
--The winning entry in a "What were HAL's first words" contest judged
by 2001: A SPACE ODYSSEY creator Arthur C. Clarke

Can't stop the signal.

(attachments)

interp.png

On 21/02/08 20:42, Maciej Sieczka wrote:

Moritz Lennert pisze:

On 13/02/08 18:40, G. Allegri wrote:

It's the same question I was going to ask. It would be an important
advance in Grass, as I'm obliged to use commercial sw to do it (Esri)
both to interpolate with barriers and to do 3D TINs with constraints.

The answer was (and as far as I know still is) to either run v.surf.rst several times, once for each part within fault lines (using mask) or to use r.surf.nnbathy (AddOns).

r.surf.nnbathy does not support fault lines either.

Ah, ok. Sorry, my bad memory :wink:

FWIW, there is surfit [1] (FOSS). It supports fault lines AIUI - linear features which don't have a value assigned, but only depict a curve at which the interpolation should break. They call break/fault lines "inequalities". I still haven't used surfit, but from examples [2] it seems interesting.

It definitely looks interesting and as it is a tcltk library, it should be very easy to integrate this into a grass script module à la r.surf.nnbathy, i.e. export GRASS data to ascii, import into surfit and create surface, export to ascii and import into GRASS.

Moritz