Program output is either [input - input * filter] or optionally [input * filter].
The data was then filtered (Fig. 3) using:
Subtraction of the output from the input results in a difference plot (Fig. 4) exposing the data removed during filtering. The effects of the filter are made obvious when the data are viewed in the Fourier domain. This was accomplished using:
fftxy -N TimeSliceFkkstrip -O 800msFkstripFft -rs 200 -re 200
The output records (Figs. 5 & 6) display the 2 dimensional spatial Fourier transforms of the input data both pre and post filter application. The axes are labeled in samples and extend from the positive to negative Nyquist frequency in both the Kx and Ky directions. The data of interest is located between the origin and 0.5 Nyquist in both orientations.
The post fkkstrip spectrum (Fig. 6) shows the level of attenuation achieved.
The filtering may have been a little too severe as what appears to be a coherent event [marked by the arrows in Fig. 5] has been attenuated. This energy corresponds to high frequency events on the input TI slice. Though probably not optimum, for this example no further parameterization of the filter was attempted. In production use it is expected that the processor would quality control the filter design so as not to remove any primary reflections of interest.
Should a time variant implementation be required one could use a combination of fftxy and polymute by picking a suite of polygons to apply to the transformed TI sorted data. In this example however a single fkkstrip time invariant filter was applied to the entire TI volume.
ttds3d -N
timeSliceFkkstrip -ND xyt -OD tyx -OD IsliceFkkstrip
The dramatic reduction in noise seen on this presentation (Fig's 7 & 8) is quite remarkable when considering the minimal smoothing observed on the time slices. The difference plot (Fig. 9) reveals that mainly random noise, but also some coherent energy, had been removed. This coherent energy corresponds to events that dip steeply in time and thus land within the reject frequency range on the TI sorted data.
Comparison of the pre and post fkkstrip LI slice data
(Fig's 1 & 10) reveals that the TI slice processing has been very effective at attenuating both random noise and steeply dipping events. Some attenuation of energy at shallower dips is also evident especially on the difference plot (Fig. 11). This was disconcerting, however, an analysis of the DI slice data (Fig's 7, 8 & 9) revealed that these events are associated with steeply dipping arrivals in the crossline direction which fall outside the pass band of the filter applied to the TI slice data.
The filter may have been more effective if convolved with the original dataset as the original crossline sample interval was 25 meters [the processed volume was sampled using a 50 meters crossline sample interval]. It would also have been instructive so see the fourier spectrum of the data prior to decimation in order to determine if an anti-alias filter was required in advance of resampling.
A comparison of the Fourier transforms of the LI records before and after filtering (Fig's 12 &13) shows that the random noise level has been lowered not just at high wavenumbers [as expected], but also over the pass band of the filter. This both explains the unexpected improvement in signal continuity and illustrates the power of time slice processing.
-l the filter length. This parameter specifies the number of points [N], on one side of the filter. The filter used will contain N2 number of points. This parameter greatly affects the spacial rejection characteristics of the filter. The greater the number of points in the filter the lower the spatial cutoff frequency.
-i the filter intercept. This parameter controls the number of lobes of the Bessel function used in the filter. A value less than 0.35 is recommended. Values greater that this may be used but could result in an aliased filter design.
-ftaper the percent of the filter edge to multiply by a cosine taper. The greater this value the lower the resulting sidelobes on the frequency response of the filter.
The command line entries affecting data conditioning prior to filtering are:
-xtaper,-ytaper the x and y edge tapers. These entries are the number of samples at the edge of the input TI slice over which to use a cosine taper.
-xpad,-ypad the x and y edge pads. These entries specify the number of samples to use as a pad around the input data.
fkkstrip -N
in -O out -l 7
The 2-d amplitude spectrum (Fig's 14, 15 & 16) of the output data illustrates that as the filter length is increased the spatial cutoff frequency is reduced.
To demonstrate the effect of the intercept parameter the following was run:
fkkstrip -N
in -O out -i 0.005 -l 31
In this instance the resulting amplitude spectrums (Fig's 17 & 18) reveal that as the intercept is increased the depth of the first sidelobe in [Kx,Ky] is decreased.
To illustrate the effect of the cosine taper the following was run:
Without a taper there is no guarantee that the filter edges are anywhere close to zero amplitude. As the taper width is increased the filter becomes more stable [the amount of filter ringing in (Kx,Ky) is reduced] (Fig's 19 & 20).
A quick look at the 2-d amplitude spectrum of the output (Fig. 21) will verify that this has occurred.
2. Since fdslice has no way of remembering the original data indexing it is necessary to run hdrswap on the final output accessing the headers of the original input dataset. It is also necessary to restore the line header entry for sample interval SmpInt [using utop] as it will be set to unity by fdslice.
The filter has also been applied to map data in order to provide regional / residual identification. This approach has met with great success.
This field of image processing of geophysical data is quite young although image processing itself has a long history in the medical, EE and space sciences fields. If you can think of some other off the wall use for this type of filter give us a call and we will endeavour to implement your algorithm in a timely fashion.
Return to the USP Home Page.Frame Doc to html conversion by Tony Garossino