next up previous [pdf]

Next: STEEP-DIP DECON Up: PATCHING TECHNOLOGY Previous: Designing a separate filter

Triangular patches

I have been running patching code for several years and my first general comment is that realistic applications often call for patches of different sizes and different shapes. (Tutorial, non-interactive C code is poorly suited to this need.) Raw seismic data in particular seems more suited to triangular shapes. It is worth noting that the basic concepts in this chapter have ready extension to other shapes. For example, a rectangular shape could be duplicated into two identical patches; then data in one could be zeroed above the diagonal and in the other below; you would have to allow, of course, for overlap the size of the filter. Module pef [*] automatically ignores the zeroed portion of the triangle, and it is irrelevant what mis2() [*] does with a zeroed portion of data, if a triangular footprint of weights is designed to ignore its output.

EXERCISES:

  1. Code the linear operator $ \bold W_{\rm wall} \bold P\T \bold W_{\rm wind} \bold P $ including its adjoint.
  2. Smoothing program. Some familiar operations can be seen in a new light when done in patches. Patch the data. In each patch, find the mean value. Replace each value by the mean value. Reconstruct the wall.
  3. Smoothing while filling missing data. This is like smoothing, but you set window weights to zero where there is no data. Because there will be a different set of weights in each window, you will need to make a simple generalization to mkwallwt [*].
  4. Gain control. Divide the data into patches. Compute the square root of the sum of the squares of all the data in each patch. Divide all values in that patch by this amount. Reassemble patches.


next up previous [pdf]

Next: STEEP-DIP DECON Up: PATCHING TECHNOLOGY Previous: Designing a separate filter

2013-07-26