[OpendTect_Users] fault density along 2d line

[mick]

how about search for the pattern 11 and replace with 01 then add up the "1"s
I'm sure someone can figure out how to program that.
How this works in 2D may be tricky if the fault runs parallel to a line.

Mick

---- Laurent.Langhi at csiro.au wrote: 
Hi OD users,
here is a little problem I want to share with you.
I have to capture the fault density along a specific horizon, I am working with 2d seismic. So I have created a similarity attribute capturing the fault along the 2d line, I turned the similarity attribute into an attribute where for each trace "1" is the fault and "0 " is the non-fault.
So, along a specific horizon, this is an example of what I get: 3 faults, a narrow one, a average one and a wide one (from left to right) (there are 44 traces)
00000100000000001110000000000111111111100000

Initially to get a density of the fault along the horizon I thought to just count, along my horizon, the number of "1" that are with a certain range from a sample point (using Volume Statistics and output Sum). But you can see the issue here: if I use an increment of 5 traces on each side of the sample, for the wide fault (far right) I ll count ten "1" and for the narrow fault (far left) I ll count one "1". However the density of the fault is not higher on the right than on the left.
Therefore my idea is to turn each group of "1" between two "0" (ie a fault) into one "1" between two "0" like this:
The initial data <00000100000000001110000000000111111111100000> becomes <00000100000000000010000000000000000000100000>
Then every "1" is a fault and therefore if I count the "1" along my horizon, within a defined trace increment range, I should have a fault density.

Any idea how to turn that <00000100000000001110000000000111111111100000> into that <00000100000000000010000000000000000000100000> ?
Alternatively if you have another way to get the fault density along a 2d line I ll be interested.
Cheers
Laurent



Dr Laurent Langhi| Structural Geology - Exploration Geophysics | Petroleum and Geothermal Research  |
CSIRO Earth Science and Resource Engineering | P +61 8 6436 8741 | F +61 8 6436 8555 | 26 Dick Perry Av, Kensington, WA 6151, AUSTRALIA | www.csiro.au/org/CESRE.html<http://www.csiro.au/org/CESRE.html>