9512.net

# Scale Attributes from Continuous Wavelet Transform

INT 3.1

Scale Attributes from Continuous Wavelet Transform
Satish K. Sinha*, School of Geology and Geophysics, University of Oklahoma, Norman, OK; Partha S. Routh, Dept. of Geosciences, Boise State University, Boise, ID; Phil D. Anno, Seismic Imaging and Prediction, ConocoPhillips., Houston, TX; John P. Castagna, University of Houston, Houston, TX.
Summary Average instantaneous attributes of time-frequency decompositions are useful in revealing the time varying spectral properties of seismic data. In the continuous wavelet transform (CWT), a time signal is decomposed into a time-scale spectrum or a scalogram; unlike a timefrequency spectrum or a spectrogram from the short time Fourier transform (STFT). Although there are various approaches of converting a time-scale spectrum into a timefrequency spectrum we introduce new mathematical formulas to calculate spectral attributes from the scalogram. In this process, we bypass the conversion of a scalogram into a time-frequency spectrum and provide average spectral attributes based on scale. The attributes are: center frequency, dominant frequency, and spectral bandwidth. Since these attributes are based on the CWT, computation of these attributes avoids subjective choice of a window length. Introduction Spectral decomposition of non-stationary signals, like seismic signals, is conventionally achieved by the STFT. Therefore, average instantaneous spectral attributes, such as center frequency, dominant frequency etc., from a spectrogram inherits the properties of the STFT. A relatively new approach of spectral decomposition based on the CWT avoids the subjective choice of a window length and produces a time-scale spectrum also called scalogram. Sinha et al. (2003) converted a scalogram into a timefrequency spectrum called the time-frequency from CWT (TFCWT). One can use the TFCWT to compute instantaneous spectral attributes using the formulas outlined by Barnes (1993). However, we define new formulas to compute the instantaneous spectral attributes directly from a scalogram. In this definition we utilize the fact that the dilating support of the given wavelet, i.e. scale, is inversely proportional to the center frequency of the wavelet (Abry et al., 1993). We use Morlet wavelet for which this mapping is a good approximation. Theory Instantaneous spectral attributes of center frequency, dominant frequency, and spectral bandwidth are defined as various moments of a time-frequency distribution using familiar definitions from probability theory (Barnes, 1993). Mean of an instantaneous power spectrum is called instantaneous center frequency and is defined as the first moment of the time-frequency power spectrum along the frequency axis. Mathematically, it is given by,

f W ( f , t ) df f C (t )
0

2

,

--- (1)

W ( f , t ) df
0

2

where

f C (t ) is instantaneous center frequency. W ( f , t )

is the time-frequency map that can be obtained from STFT, TFCWT or any other spectral decomposition. Square root of the second moment of the time-frequency power spectrum along the frequency axis is called instantaneous dominant frequency.

f 2 W ( f , ) df
2 fD ( ) 0

2

,

--- (2)

W ( f , ) df
0

2

where

f D (t ) is instantaneous dominant frequency.

Standard deviation about the instantaneous center frequency is called instantaneous spectral bandwidth.

(f
2 f BW ( ) 0

f C ( )) 2 W ( f , ) df
, --- (3)

2

W ( f , ) df
0

2

where

f BW (t ) instantaneous spectral bandwidth. All the

above equations with different notations have been adopted from the paper by Barnes (1993) though competing equations for all of them exist in the literature. For a symmetrical wavelet, like Morlet wavelet (Torrence and Compo, 1998), center frequency of the wavelet at each scale can be assumed to be inversely proportional to the scale. Mathematically, it can be written as,

fC

k

,

--- (4)

SEG/Houston 2005 Annual Meeting

779

INT 3.1

Scale Attributes from Continuous Wavelet Transform

where

fC ,

, and

k are center frequency, scale and a

wavelet dependent constant respectively. The center frequency can be considered as a representative frequency for a particular scale. Differentiating equation 4, and replacing f C with f , we obtain

(equations 6, 7, and 8) are computationally efficient at the cost of slight smoother representation of the information. This efficiency translates to saving considerable amount of time depending upon length and number of traces in a large 3D data set. The next logical step is to compare results on a real seismic section and examine if there are any difference in interpretation. Instantaneous dominant frequency from TFCWT using equation 2 is shown in Figure 2. Instantaneous dominant frequency from scalogram using equation 5 is displayed in Figure 3. In both Figures (2 and 3) the frequency range on the colorbar is kept the same. Figure 2 and 3 indicates remarkable similarity and virtually can be considered as identical. Difference in timefrequency resolution is minimal. It is important to note that the computation of Figure 3 is about 5 times faster than Figure 2. Therefore, scale based attributes are more efficient from practical point of view. Figure 4 is produced using the STFT with a “Hanning” window of 200 ms window length.

df

k
2

d .

--- (5)

Changing from frequency to scale in the equations 1, 2, and 3 and assuming the wavelet dependent constant k 1 , we obtain

1 ~ f C (t )
3 0

FW ( , t ) d
, --- (6)

2

1
2 0

FW ( , t ) d FW ( , t ) d
, --- (7)
2

2

1 ~2 f D (t )
and
4 0

1
2 0

FW ( , t ) d
15

2

(a) Seismic Trace

(1 ~2 f BW ( )
0

~ 2 fC )
4

10

Amplitude

5 0 -5 -10 -15

FW ( , t ) d
, --- (8)
2

2

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

(b) Center Frequency
90

2 0

FW ( , t ) d

Frequency (Hz)

1

80 70 60 50 40 30 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

~ ~ ~ where f C , f D , and f BW are instantaneous center
frequency, instantaneous dominant frequency and instantaneous bandwidth respectively from the scalogram FW ( , t ) . Thus equation 6, 7 and 8 allows us to compute these spectral attributes directly from a scale map. Example Comparisons of instantaneous spectral attributes based on frequency and scale for a seismic trace (Figure 1a) are shown in Figures 1-b, c, and d. We note the remarkable similarity between these attributes although the attributes computed from the scalogram seems to be smoother. This is expected because a scale represents a band of frequency and not a single frequency. Therefore, overlap of frequencies from one scale to another produces smoother attributes, but the overall agreement is excellent. The favorable comparison is a consequence of the scale to center frequency mapping for a Morlet wavelet. Computation of attributes based on scalogram is much faster than those based on TFCWT. Therefore, from practical point of view, attribute definitions based on scale

From TFCWT

From Scalogram

90

(c) Dominant Frequency

Ti

( )

Frequency (Hz)

80 70 60 50 40 30 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

From TFCWT

From Scalogram

(d) Bandwidth
40

Bandwidth (Hz)

35 30 25 20 15 10 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7

From TFCWT

From Scalogram

0.8

0.9

1

Time (s)

Figure 1: a) A seismic trace of length 1sec sampled at 4ms. It was used to compute frequency and scale based attributes. b) Instantaneous center frequency attribute, c) instantaneous dominant frequency attribute, and d) instantaneous spectral bandwidth computed from TFCWT (shown in black) and from scalogram (shown in red).

SEG/Houston 2005 Annual Meeting

780

INT 3.1

Scale Attributes from Continuous Wavelet Transform

Conclusions Scalogram obtained from continuous wavelet transform can provide useful average measures that are directly interpretable in terms of time-frequency attributes. We obtain new expressions for commonly used attributes such as dominant frequency, center frequency and bandwidth directly from a scalogram. These attributes can be efficiently calculated without converting a scalogram into a time-frequency spectrum. Field example presented in this paper shows a factor of 5 computational efficiency using scale based computation. In addition these attributes avoids the subjective choice of a window length commonly used in STFT. This is because the wavelet acts as an effective window via dilatation process.

Figure 4: Instantaneous dominant frequency from STFT using 200 ms window length.

Acknowledgements The authors acknowledge the Conoco research scientists for their valuable inputs. We would also like to thank Conoco Inc. (now ConocoPhillips) for providing the field dataset. We thank Dr. Dale Cox of Conoco for many useful discussions and for his help in the practical implementation in SU. References
Figure 2: Instantaneous dominant frequency from TFCWT.

Abry, P., Goncalves, P., and Flandrin, P., 1993, Waveletbased spectral analysis of 1/f processes: IEEE Internat. Conf. Acoustic, Speech and Signal Processing, 3, 237-240. Barnes, A. E., 1993, Instantaneous spectral bandwidth and dominant frequency with applications to seismic reflection data: Geophysics, 58, 3, 419-428. Sinha, S., Routh, P., Anno, P. and Castagna, J., 2003, Time-frequency attribute of seismic data using continuous wavelet transform, 73rd Ann. Internat. Mtg.: Soc. of Expl. Geophys., Expanded Abstracts, 1481-1484. Torrence, C., and Compo, G., P., 1998, A practical guide to wavelet analysis: Bulletin of the American Meteorological Society, 79, 1, 61-78.

Figure 3: Instantaneous dominant frequency from scalogram.

SEG/Houston 2005 Annual Meeting

781

EDITED REFERENCES Note: This reference list is a copy-edited version of the reference list submitted by the author. Reference lists for the 2005 SEG Technical Program Expanded Abstracts have been copy edited so that references provided with the online metadata for each paper will achieve a high degree of linking to cited sources that appear on the Web. Scale Attributes from Continuous Wavelet Transform REFERENCES Abry, P., P. Goncalves, and P. Flandrin, 1993, Wavelet-based spectral analysis of 1/f processes: IEEE International Conference on Acoustic, Speech and Signal Processing, 3, 237-240. Barnes, A. E., 1993, Instantaneous spectral bandwidth and dominant frequency with applications to seismic reflection data: Geophysics, 58, 419-428. Sinha, S., P. Routh, P. Anno, and J. Castagna, 2003, Time-frequency attribute of seismic data using continuous wavelet transform: 73rd Annual International Meeting, SEG, Expanded Abstracts, 1481-1484. Torrence, C., and G. P. Compo, 1998, A practical guide to wavelet analysis: Bulletin of the American Meteorological Society, 79, 61-78.

Scale Attributes from Continuous Wavelet Transform_图文.pdf
Scale Attributes from Continuous Wavelet Transform - Average instantaneous attributes of time-fre...

...Based on Support Vector Machines and Multi-Scale Wavelet ....pdf
A wavelet transform of the noisy speech is needed to get the wavelet transform coefficient on different scales. Using a threshold function to deal with ...
Medical Image Fusion Based on Wavelet Multi-Scale....pdf
Medical Image Fusion Based on Wavelet Multi-ScaleDecomposition_基础医学_医药...According to the choice of multi-wavelet based on multi-wavelet transform ...
A new method of choosing scales in wavelet transform for ....pdf
A new method of choosing scales in wavelet transform for damping identification_交通运输_工程科技_专业资料。维普资讯 http://www.cqvip.com Junl ab stto...
...emitter signal recognition based on multi-scale wavelet ....pdf
Radar emitter signal recognition based on multi-scale wavelet entropy and feature weighting_能源/化工_工程科技_专业资料。J.Cent.South Univ.(2014)21:4254...
The FBI Wavelet Scalar Quantizat ion Standard for gray-scale ....pdf
3 THE WAVELET TRANSFORM SUBBAND DECOMPOSITION A ...continuous-tone images as far as the compression ... Scale Attributes from ... 4页 免费 Edge ...
Multiscale edge detection of noisy images using wavelets_论文....pdf
Multiscale edge detection of noisy images using wavelets_信息与通信_工程科技_专业资料。维普资讯 http://www.cqvip.com JunlfHri sito cnly(e ee)o.4N.,...
...continuous spherical wavelet transform algorithm....pdf
The Fast Lifting Wavelet... 暂无评价 27页 免费 Scale Attributes from Co... Fast directional continuous spherical wavelet transform algorithms We describe ...
Parallel Wavelet Transform for Large Scale Image.pdf
Parallel Wavelet Transform for Large Scale Image_专业资料。In this paper we discuss several issues relevant to the parallel implementation of a 2-D ...
...detection using the dual-tree complex wavelet transform_....pdf
Multiscale keypoint detection using the dual-tree complex wavelet transform_专业资料。We show that the Dual-Tree Complex Wavelet Transform (DTCWT) [1] ...
WAVELET FOLDING AND DECORRELATION ACROSS THE SCALE.pdf
WAVELET FOLDING AND DECORRELATION ACROSS THE SCALE_专业资料。The discrete wavelet transform (DWT) gives a compact multiscale representation of signals and ...
Multiscale fundamental forms a multimodal image wavelet ....pdf
towards multiscale fundamental forms using the dyadic wavelet transform of ...proposed to orient the eigenvectors, keeping the angle continuous in local ...
...texture segmentation of dip-cube slices using wavelet-....pdf
In this paper, we propose a multiscale texture segmentation algorithm based on the wavelet transform of the image. For many 1 classes of signals and ...
...scale Classification of Document Images Based on Wavelet ....pdf
Context Based Multiscale Classification of Document ...(t) can be estimated by a continuous ...Wavelet Transform," Proceedings of Third ...
...scale parallel implementation of fast wavelet transforms_....pdf
[6] Judith Ford, Ke Chen, and Laurence Scales. A new wavelet transform preconditioner for iterative solution of elastohydrodynamic lubrication problems. Int....
MULTISCALE OBJECT FEATURES FROM CLUSTERED COMPLEX WAVELET ....pdf
scales, derived from dual-tree complex wavelets, ...cients. The ILP transform is a pyramid of ... however, is a set of continuous complex ILP ...
Wavelet Analysis for a New Multiresolution Model for Large-....pdf
Model for Large-Scale Textured Terrains_专业资料。...The Wavelet transform produces a hierarchical ...Real-Time, Continuous Level of Detail Rendering ...
...feature extraction and object recognition with wavelets ....pdf
Multiscale geometrical feature extraction and object recognition with wavelets and morphology_专业资料。MAllat的好文 MULTISCALE GEOMETRICAL RECOGNITION WITH Seema,...
Wavelet transform_图文.ppt
? ? ? Wave and Wavelet Scaling functions Mother wavelet Wavelet functions ...Continuous wavelet transform(CWT) Discrete wavelet transform(DWT) 2012-4-6 ...