Seismic data compression by an adaptive local cosine/sine transform and its effects on migration
|
|
- Ethelbert Morton
- 5 years ago
- Views:
Transcription
1 Geophysical Prospecting, 2000, 48, 1009±1031 Seismic data compression by an adaptive local cosine/sine transform and its effects on migration Yongzhong Wang 1 and Ru-Shan Wu 1 Abstract The local cosine/sine basis is a localized version of the cosine/sine basis with a window function which can have arbitrary smoothness. It has orthogonality and good time and frequency localization properties. The adaptive local cosine/sine basis is a best-basis obtained from an overabundant library of cosine/sine packets based on a costfunctional. We propose a 2D semi-adaptive (time-adaptive or space-adaptive) local cosine transform (referred to as a 2D semi-alct) and apply it to the SEG±EAEG salt model synthetic data set for compression. From the numerical results, we see that most of the important features of the data set can be well preserved even in the high compression ratio (CRˆ40:1) case. Using reconstructed data from the highly compressed ALCT coefficients (CRˆ40:1) for migration, we can still obtain a highquality image including subsalt structures. Furthermore, we find that the window partition, generated by the 2D semi-alct, is well adapted to the characteristics of the seismic data set, and the compression capability of the 2D semi-alct is greater than that of the 2D uniform local cosine transform (2D ULCT). We find also that a (32, 32) or (32, 64) minimum (time, space) window size can generate the best compression results for the SEG±EAEG salt data set. Introduction In seismic data compression, the problem is the selection of suitable bases to represent seismic signals efficiently and therefore to achieve maximum possible compression, while still preserving useful information. This problem has attracted considerable attention (e.g. Donoho, Ergas and Villasenor 1995; Bradley, Fei and Hildebrand 1996; Villasenor, Ergas and Donoho 1996; Wang and Pann 1996; Vassiliou and Wickerhauser 1997; Wu and Wang 1999) and much progress has been made in the past few years. For example, Vassiliou and Wickerhauser (1997) used bi-orthogonal wavelets as the expansion bases for the compression of data sets from Egypt and Trinidad. Their numerical results show that long filters, moderate decomposition depths and frequency-weighted, variance-adjusted quantization yield the best results. Received April 1999, revision accepted March Institute of Tectonics, University of California, Santa Cruz, CA 95064, USA. q 2000 European Association of Geoscientists & Engineers 1009
2 1010 Y. Wang and R.-S. Wu Figure 1. A schematic representation of 2D time-alct, showing uniform LCT along the space axis, and within each fixed space strip, adaptive LCT along the time axis. From their conclusions, we see that compression effectiveness varies considerably for different expansion bases. Local cosine/sine bases are orthogonal bases with good localization in both the space/time and the wavenumber/frequency domains. The adaptive local cosine/sine basis is a best-basis obtained from an overabundant library of cosine/sine packets based on a cost-functional. The adaptability results in a flexible space and/or time segmentation to match the signal characteristics. By contrast, the uniform local cosine transform (ULCT) cannot generate an adaptive matching to a given signal because the width of each window is fixed and equal. We propose a 2D semi-adaptive local cosine transform (either adaptive in the time direction and uniform in the space direction, referred to as 2D time-alct, or adaptive in the space direction and uniform in the time direction, referred to as 2D space-alct) and apply it to seismic data compression of the SEG±EAEG salt model synthetic data set Figure 2. A schematic representation of 2D space-alct, showing uniform LCT along the time axis, and within each fixed time strip, adaptive LCT along the space axis.
3 Adaptive local cosine/sine transform 1011 Figure 3. A sketch drawing of folding on the overlapping zone. b I (x) is a bell function whose nominal support is [a,b]. 1 is the left overlapping radius and 1 0 is the right overlapping radius. Figure 4. A 1D adaptive local cosine multilevel decomposition tree for the 200th trace of the SEG±EAEG salt data set. The preset maximum decomposition level is 4 and Shannon entropy is chosen as the cost-functional.
4 1012 Y. Wang and R.-S. Wu Figure 5. Adaptive window segmentation of the 200th trace of the SEG±EAEG salt data set. The minimum window width is 64 samples, the smoothness parameter n in (9) is 2, the preset maximum decomposition level is 4 and Shannon entropy is chosen as the cost-functional. The dotted lines indicate the edges of the nominal supports of the windows. The solid line shows the input trace. as a test example. The 1D ALCT can be implemented with a fast algorithm of computational complexity O(Nd log 2 N) (Coifman and Wickerhauser 1992), where N is the number of samples in the signal and d is the maximum decomposition level. Adaptive local cosine/sine transform Brief description of the local cosine/sine basis Local cosine/sine bases constructed by Coifman and Meyer (1991) (see also Pascal, Guido and Wickerhauser 1992) consist of cosines/sines multiplied by smooth, compactly supported bell functions. These localized cosine/sine functions remain orthogonal and have small Heisenberg products. The local cosine/sine transform has much in common with the windowed or short-time Fourier transform (WFT or
5 Adaptive local cosine/sine transform 1013 STFT) (Gabor 1946). However, for the latter, the Balian±Low obstruction (Balian 1981) prevents the windowed exponential basis from simultaneously being a frame and having a finite Heisenberg product. Local cosine or sine wavelets, on the other hand, overcome this limitation. The basis element (here, for example, the basis of the fourth type) can be characterized by position a, interval I and frequency/wavenumber index k as follows: or r 2 c aik x ˆ b I x cos p k 1 x 2 a jij 2 jij r 2 c aik x ˆ b I x sin p k 1 x 2 a ; 2 jij 2 jij where b ˆ b I x is a bell function which is a smooth function compactly supported in the interval a 2 1; b 1 0 Š for a 1 # b : This interval contains I ˆ a; bš which is called its nominal support. We can control the window width by I, the left endpoint of the nominal support by the position a and the frequency/wavenumber by the index k. The properties of the bell function, i.e. 8 b I x 2 b I 2a 2 x 2 ˆ 1; x [ a 2 1; a 1Š; >< b I x 2 b I 2b 2 x 2 ˆ 1; x [ b ; b 1 0 Š; 3 >: b I x ˆ1; x [ a 1; b ; ensure the orthonormality of the bases and afford a fast algorithm. Suppose a sequence {a j } is selected to satisfy a j, a j 1 ; lim j!^1 {a j } ˆ ^1 and there also exists an accompanying sequence {1 j } such that a j 1 j # a j j 1 for all j [ Z: Then the functions s 2 c jk x ˆ b aj ;a a j 1 2 a j 1 Š x cos p k 1 x 2 aj ; 4 j 2 a j 1 2 a j j [ Z; k ˆ 0; 1; 2; ¼ ; with positive polarity at a j and negative polarity at a j+1,or s 2 c jk x ˆ b aj ;a a j 1 2 a j 1 Š x sin p k 1 x 2 aj ; 5 j 2 a j 1 2 a j 1 j [ Z; k ˆ 0; 1; 2; ¼ ; with negative polarity at a j and positive polarity at a j+1,form an orthonormal basis for the space L 2 (R). A two-dimensional local cosine/sine basis can be generated by the tensor products
6 1014 Y. Wang and R.-S. Wu Figure 6. Window partition, generated by the 2D time-alct, for the SEG±EAEG data set. The minimum time-window width is 32 samples and the uniform space-window width is 64 traces. Shannon entropy is selected as the cost-functional. c ax I x k x x c ay I y k y y and their nominal supports are the Cartesian product rectangles of the nominal supports of the x and y factors. Adaptive local cosine/sine basis 1D ALCT For binary decomposition, if we fix the window width I jij ˆja j 1 2 a j j ; all wavelets in (4) or (5) can form an orthonormal basis for L 2 (R) and this basis is called level 1; if we then select I /2 as the window width, the wavelets in (4) or (5) can also form an orthonormal basis and this is called level 2; and so on. Thus, a binary-based signal decomposition tree consists of the bases at different levels. However, not all the bases are efficient in matching a given signal. Therefore, we must pick the `best-basis' from all the possible local cosine/sine bases, using a cost-functional. To search for the local cosine/sine best-basis, i.e. an adaptive local cosine/sine basis, in order to achieve the best matching to the signal, a cost-functional is defined based
7 Adaptive local cosine/sine transform 1015 Figure 7. Window partition, generated by the 2D space-alct, for the SEG±EAEG data set. The uniform time-window width is 32 samples and the minimum space-window width is 32 traces. Shannon entropy is selected as the cost-functional. on the entropy of the decomposed signal. There are numerous cost-functionals cited in the literature (see for example Misiti et al. 1996). Here, we use Shannon entropy as the cost-functional, i.e. E X ˆ2 X p n log 2 p n ; 6 n where p n ˆjx n j 2 =kxk 2 ; X is the signal of length N samples, x n is the nth component of X and E(X) denotes the entropy of the signal X. In implementation, we use Coifman and Wickerhauser's (1992) fast algorithm to search for the best-basis based on the Shannon entropy. The main concept of this algorithm is that the full local cosine/sine tree is pruned recursively at each node by comparing its entropy to the summation of the entropy of its corresponding child nodes, i.e. IF Entropy parent node # Entropy child1 Entropy child2 Š THEN cut off the child branches: Initially, a full binary-based decomposition tree with a preset maximum
8 1016 Y. Wang and R.-S. Wu decomposition level is produced. The pruning procedure then starts from the leaf nodes and proceeds towards the root. At the end of this procedure, an optimal pruned tree is obtained for the given signal, i.e. an adaptive local cosine/sine basis is obtained. For more details, see Coifman and Wickerhauser (1992). 2D semi-alct For the 2D case, we propose two semi-adaptive schemes, i.e. 2D time-alct and 2D space-alct. The 2D time-alct is uniform in the space direction and adaptive in the time direction. Its adaptability is accomplished by the above-mentioned 1D ALCT algorithm along the time direction for a strip with fixed width (see Fig. 1). For the 2D space-alct, the segmentation is uniform in the time direction but adaptive in the space direction (see Fig. 2). Implementation Figure 8. SEG±EAEG salt velocity model A±A' (velocity in feet/second). Bell function Generally, the bell function over I ˆ a; bš is defined by b I x ˆS 1 x 2 a C 1 0 x 2 b ; x [ R; 7
9 Adaptive local cosine/sine transform 1017 Figure 9. The zero-offset seismic section from the model generated by exploding reflector simulation. where S 1 x 2 a ; sin u n and x2b C1 0 x 2 b ; cos u n : We select 8 >< u 0 x ˆ and >: 8 >< u n x ˆ >: x2a 1 0; if x,21; p 1 x ; 4 if 2 1 # x # 1; u n21 p ; if x. 1; 2 0; if x,21; sin p 2 x ; if 2 1 # x # 1; p ; if x. 1: Using induction we can show that S 1 x 2 a has 2 n 2 1 vanishing derivatives at x ˆ a 2 1 and x ˆ a 1; and C 1 0 x 2 b has 2 n 2 1 vanishing derivatives at x ˆ
10 1018 Y. Wang and R.-S. Wu Figure 10. Reconstructed data with compression ratio (CR) of 26:1. b and x ˆ b 1 0 : So S 1 x 2 a or C 1 0 x 2 b can be made into an arbitrarily smooth cut-off function. We take n ˆ 2: Folding Rather than calculating inner products with the sequences c a Ik, we can preprocess data so that the standard fast discrete cosine transform of the fourth type (DCT-IV) or discrete sine transform of the fourth type (DST-IV) algorithm may be used. This can be realized by folding the overlapping parts of the bell functions back into the intervals. Suppose we wish to fold a signal f(x) back into the interval I ˆ a; bš across a and b (see Fig. 3). Using the bell function b I (x) defined by (7), we have (using the local cosine transform of the fourth type (LCT-IV) as an example) 8 f x ˆb I x f x b I 2a 2 x f 2a 2 x ; if a # x # a 1; >< f new x ˆ f 2 x ˆb I x f x 2 b I 2b 2 x f 2b 2 x ; if b # x # b; >: f x ; if a 1, x, b : 10 The resultant folded data f new (x) is now defined in the interval [a,b].
11 Adaptive local cosine/sine transform 1019 Figure 11. Residual data between Fig. 9 and Fig. 10 for CR ˆ 26:1 (amplified by a factor of 10). To reconstruct f(x) from f new (x), we can use the following unfolding formulae (again using the LCT-IV as an example): 8 b I x f new x 2 b I 2a 2 x f new 2a 2 x ; if a # x # a 1; >< f x ˆ b I x f new x b I 2b 2 x f new 2b 2 x ; if b # x # b; >: f new x ; if a 1, x, b : 11 Edge extension When the bell shifts to the leftmost endpoint or the rightmost endpoint of the signal, we cannot directly obtain f + (x)orf 2 (x) from the above formulae because of the lack of data in the leftmost and rightmost overlapping zones. Usually, we have four extension schemes: (1) zero-extension; (2) symmetry-extension; (3) smoothness-extension; (4) periodization-extension. In this work, based on the features of seismic signals, we chose the zero-extension.
12 1020 Y. Wang and R.-S. Wu Figure 12. Reconstructed data with CR ˆ 40:1. DCT-IV/DST-IV After the procedures of folding and edge extension, we can apply the fast DCT-IV/ DST-IV to the folded data f new (x) to obtain the local cosine/sine transform coefficients. Searching for the best-basis After we have completed local cosine/sine transforms for all the preset decomposition levels and also calculated the entropy of each transformed subsignal, we can use the Coifman and Wickerhauser (1992) fast algorithm to search for the local cosine/sine best-basis. As an example, the 1D ALCT is applied to the 200th trace of the SEG±EAEG salt data set (the original trace length was 626, but is extended here to 1024 by zero-padding). The minimum window width for this trace is preset at 64 samples, the smoothness parameter n in (9) is chosen as 2, and Shannon entropy is adopted as the cost-functional. Figure 4 indicates its adaptive local cosine multilevel decomposition tree (the preset maximum level is 4). Figure 5 shows its adaptive window segmentation of the trace resulting from Fig. 4. In Fig. 5, the dotted lines
13 Adaptive local cosine/sine transform 1021 Figure 13. Residual data between Fig. 9 and Fig. 12 for CR ˆ 40:1 (amplified by a factor of 10). indicate the edges of the nominal supports of the windows while the solid line shows the input trace. Figures 6 and 7 show the window partitions generated by the 2D semi-alct for the SEG±EAEG data set (in Fig. 6 by the 2D time-alct and in Fig. 7 by the 2D space-alct). In Fig. 6, the minimum time-window width is 32 samples and the uniform space-window width is 64 traces, while in Fig. 7, the uniform time-window width is 32 samples and the minimum space-window width is 32 traces. In Figs 6 and 7, the adaptive grids show the nominal supports of the windows of a 2D best local cosine basis. As can be seen, the adaptive grids match the seismic data well, i.e. along the time or space direction the window widths are shorter where the density of seismic signals is greater. Data compression and imaging using compressed data We now test the performances of our compression schemes (2D time-alct and 2D space-alct) and compare the results with that of 2D ULCT. We also test
14 1022 Y. Wang and R.-S. Wu Figure 14. The coefficient profile with CR ˆ 40:1 after the 2D time-alct of input data in Fig. 9. In this figure, the horizontal axis represents position-wavenumber while the vertical axis represents time-frequency. The minimum (time, space) window size for the 2D time-alct is (16, 16) and Shannon entropy is adopted. the effects of seismic data compression on the image quality obtained by depth migration. Figure 8 shows a 2D profile A±A' of the SEG±EAEG salt velocity model and Fig. 9 shows the synthetic zero-offset data from the salt model using a finite-difference exploding reflector modelling algorithm, generated at AMOCO. This is a very complicated model with a large velocity contrast (salt velocity is as high as three times that of the surrounding media). Figure 10 shows the data reconstructed from the compressed 2D time-alct coefficients with compression ratio CR ˆ 26:1. Compression is achieved by setting a threshold of 1% of the maximum absolute ALCT coefficient for each space-window. Therefore the threshold varies for different spacewindows. Figure 11 shows the corresponding residual data amplified 10 times to show detail. Reconstructed data for the CR ˆ 40:1 case (also using the 2D time-alct) is shown in Fig. 12 and its residual data (amplified by a factor of 10) is plotted in Fig. 13. Here, a threshold of 2% of the maximum absolute ALCT coefficient for each spacewindow is adopted. Figures 10±13 show that the majority of the important features of
15 Adaptive local cosine/sine transform 1023 Figure 15. Hybrid pseudo-screen migration on the original synthetic data. the original synthetic data have been preserved very well even in the case of a high compression ratio (CR ˆ 40:1). In implementing the 2D time-alct with the SEG±EAEG salt data set, we select the Shannon entropy as the cost-functional, the maximum decomposition level in the time direction is 6, the uniform window width in the space direction is 16 traces, i.e. the size of the minimum (time, space) window is (16, 16) (the original seismic data is extended to 1024 in time and 2048 in space by zero-padding), the overlapping radius in the time or space direction is one half of the minimum or uniform window width, i.e. (8, 8), and the fixed folding style is used. Also an automatic gain control (AGC) is applied to the data to preserve low-amplitude seismic signals before applying the 2D time-alct. The AGC is subsequently removed after reconstruction. Figure 14 shows the coefficient profile with CR ˆ 40:1 after the 2D time-alct of input data in Fig. 9 using the (16, 16) minimum (time, space) window. Figures 15±19 show the results of migrating the original data and reconstructed data after compression in order to investigate the influence of data compression on the migrated image. Here, the imaging method used is the hybrid pseudo-screen
16 1024 Y. Wang and R.-S. Wu Figure 16. Hybrid pseudo-screen migration on the reconstructed synthetic data from the CR ˆ 26:1 ALCT coefficients. migration (Jin, Wu and Peng 1998). Figure 15 shows migration on the original synthetic data, and Fig. 16 shows migration on the reconstructed data from the compressed coefficient data with CR ˆ 26:1. Figure 18 is the image obtained from migration on the compressed data with CR ˆ 40:1. Figures 17 and 19 show the corresponding residual images after amplification by a factor of 10. Figures 15±19 show that high-quality images including subsalt structures can still be obtained even in the high compression ratio (CR ˆ 40:1) case. To compare the compression capability of the 2D semi-alct with that of the 2D ULCT, we define the signal-to-noise ratio (SNR) as follows: X SNR ˆ 10 log 10 jc k j 2 = X j1 k j!; 2 12 k k where c k is the coefficient above the threshold in the absolute value sense and is thus retained, and 1 k is the coefficient which is discarded. Figure 20 shows the comparison of compression performance between the 2D time-alct and the 2D ULCT. From Fig. 20, we see that the compression capability of the 2D time-alct is more powerful than that of the 2D ULCT.
17 Adaptive local cosine/sine transform 1025 Figure 17. Residual image between Fig. 15 and Fig. 16 for CR ˆ 26:1 (amplified by a factor of 10). To address the relationship between the SNR (in db) of the transformed data and the compression ratio with different minimum (time, space) windows (here, adaptive along the time axis), we tested six different minimum (time, space) windows: (8, 8), (16, 8), (16, 16), (32, 16), (32, 32) and (32, 64). From the curves of SNR versus compression ratio shown in Fig. 21, we see that either (32, 32) or (32, 64) is the best choice of (time, space) windows for the SEG±EAEG data set. Finally, the compression performances of the 2D time-alct and the 2D space- ALCT were compared. In implementing the 2D space-alct, we selected only the (32, 32) and (32, 64) minimum (time, space) windows as tests because these two windows are already the best choices in the 2D time-alct for the SEG±EAEG salt data set. Figure 22 shows that the compression capabilities of the 2D time-alct and the 2D space-alct are almost the same. Conclusions We proposed a 2D semi-adaptive local cosine/sine transform to achieve data compression and we applied it to the SEG±EAEG salt data set. From the numerical
18 1026 Y. Wang and R.-S. Wu Figure 18. Hybrid pseudo-screen migration on the reconstructed synthetic data from the CR ˆ 40:1 ALCT coefficients. results, it was found that the majority of the important features of the original data set can be preserved very well even in the high compression ratio (CR ˆ 40:1) case. Furthermore, even using the highly compressed data (CR ˆ 40:1) for migration, we can still obtain high-quality seismic images including the subsalt structures. We found that the window partitions, generated by the 2D semi-alct, aligned well with natural breaks in the seismic data, therefore the compression capability of the 2D semi-alct is greater than that of the 2D ULCT. Also, the 2D time-alct and 2D space-alct are equally powerful for data compression. Moreover, we found that a (32, 32) or (32, 64) minimum (time, space) window size generated the best compression result for the SEG±EAEG salt data set. Acknowledgements This work was supported by the WTOPI (Wavelet Transform On Propagation and Imaging for seismic exploration) project at the University of California, Santa Cruz. We are grateful to the sponsors. The help from Dr Shengwen Jin on migration is
19 Adaptive local cosine/sine transform 1027 Figure 19. Residual image between Fig. 15 and Fig. 18 for CR ˆ 40:1 (amplified by a factor of 10). greatly appreciated. The facility support from the W. M. Keck Foundation is also acknowledged. Contribution number 402 of the Institute of Tectonics, University of California, Santa Cruz. References Balian R Un principle d'incertitude en theâorie du signal ou en meâcanique quantique. Comptes Rendus de l'academie des Sciences, Paris, SeÂrie II 292, 1357±1362. Bradley J., Fei T. and Hildebrand S Wavelet compression for 3D depth migration. 66th SEG meeting, Denver, USA, Expanded Abstracts, 1627±1629. Coifman R.R. and Meyer Y Remarques sur l'analyse de Fourier aá feneãtre. Comptes Rendus de l'academie des Sciences, Paris, SeÂrie I 312, 259±261. Coifman R.R. and Wickerhauser M.V Entropy-based algorithms for best basis selection. IEEE Transactions on Information Theory 38, 713±718. Donoho P.L., Ergas R.A. and Villasenor J.D High-performance seismic trace compression. 65th SEG meeting, Expanded Abstracts, Houston, USA, 160±163.
20 1028 Y. Wang and R.-S. Wu Figure 20. Compression performance comparison between the 2D time-alct and the 2D ULCT: minimum (time, space) window (16, 16) for 2D time-alct, uniform (time, space) window (16, 16) for 2D ULCT. Shannon entropy is used as the cost-functional for both schemes.
21 Adaptive local cosine/sine transform 1029 Figure 21. Compression performance comparison of different minimum (time, space) windows by 2D time-alct using Shannon entropy as the cost-functional.
22 1030 Y. Wang and R.-S. Wu Figure 22. Compression performance comparison between the 2D time-alct and the 2D space-alct.
23 Adaptive local cosine/sine transform 1031 Gabor D Theory of communication. Journal of the Institute of Electrical Engineers 93(III), 429±457. Jin S., Wu R.S. and Peng C Prestack depth migration using a hybrid pseudo-screen propagator. 68th SEG meetng, New Orleans, USA, Expanded Abstracts, 1819±1822. Misiti M., Misiti Y., Oppenheim G. and Poggi J.M Wavelet Toolbox for Use with MATLAB, pp. 5±7. The Mathworks, Inc. Pascal A., Guido W. and Wickerhauser M.V Local sine and cosine bases of Coifman and Meyer and the construction of smooth wavelets. In: Wavelets: A Tutorial in Theory and Applications (ed. Charles K. Chui), pp. 237±256. Academic Press, Inc. Vassiliou A. and Wickerhauser M.V Comparison of wavelet image coding schemes for seismic data compression. In: Wavelet Applications in Signal and Image Processing V. (eds Akram Aldroubi, Andrew F. Laine and Michael A. Unser). Proceedings of SPIE 3169, 118± 126. Villasenor J.P., Ergas R.A. and Donoho P.L Seismic data compression using highdimensional wavelet transforms. Proceedings of the Data Compression Conference, Snowbird, UT, USA, pp. 396±405. IEEE Computer Society Press. Wang B. and Pann K Kirchhoff migration of seismic data compressed by matching pursuit decomposition. 66th SEG meeting, Denver, USA, Expanded Abstracts, 1642±1645. Wu R.S. and Wang Y Seismic data compression using adapted local cosine transform and its effects on imaging. 61st EAGE conference, Helsinki, Finland, Extended Abstracts, P102.
SEG Houston 2009 International Exposition and Annual Meeting
Yu Geng* 1, Ru-Shan Wu and Jinghuai Gao 2 Modeling and Imaging Laboratory, IGPP, University of California, Santa Cruz, CA 95064 Summary Local cosine/sine basis is a localized version of cosine/sine basis
More informationSEG/San Antonio 2007 Annual Meeting
Imaging steep salt flanks by super-wide angle one-way method Xiaofeng Jia* and Ru-Shan Wu, Modeling and Imaging Laboratory, IGPP, University of California, Santa Cruz Summary The super-wide angle one-way
More informationSEG/New Orleans 2006 Annual Meeting
Accuracy improvement for super-wide angle one-way waves by wavefront reconstruction Ru-Shan Wu* and Xiaofeng Jia, Modeling and Imaging Laboratory, IGPP, University of California, Santa Cruz Summary To
More informationDownloaded 10/23/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
ACQUISITION APERTURE CORRECTION IN ANGLE-DOMAIN TOWARDS THE TRUE- REFLECTION RTM Rui Yan 1*, Huimin Guan 2, Xiao-Bi Xie 1, Ru-Shan Wu 1, 1 IGPP, Earth and Planetary Sciences Department, University of California,
More informationDistance (*40 ft) Depth (*40 ft) Profile A-A from SEG-EAEG salt model
Proposal for a WTOPI Research Consortium Wavelet Transform On Propagation and Imaging for seismic exploration Ru-Shan Wu Modeling and Imaging Project, University of California, Santa Cruz August 27, 1996
More informationMain Menu. Summary. sampled) f has a sparse representation transform domain S with. in certain. f S x, the relation becomes
Preliminary study on Dreamlet based compressive sensing data recovery Ru-Shan Wu*, Yu Geng 1 and Lingling Ye, Modeling and Imaging Lab, Earth & Planetary Sciences/IGPP, University of California, Santa
More informationA full-wave equation based seismic illumination analysis method
A full-wave equation based seismic illumination analysis method Xiao-Bi Xie and Hui Yang University of California, Santa Cruz, CA 95064, USA Summary We propose a full-wave equation based method for seismic
More informationMulti-azimuth velocity estimation
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 87 Multi-azimuth velocity estimation Robert G. Clapp and Biondo Biondi 1 ABSTRACT It is well known that the inverse problem of estimating interval
More informationSEG Houston 2009 International Exposition and Annual Meeting
Yueming Ye *1, Ru-Shan Wu and Zhenchun Li 2 Modeling and Imaging Laboratory, IGPP, University of California, Santa Cruz, CA 95064 Summary Migration with data acquired on surface with irregular topography
More informationDownloaded 05/09/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
Elastic converted-wave path migration for subsalt imaging Ru-Shan Wu*, Rui Yan, Xiao-Bi Xie, Modeling and Imaging Laboratory, Earth and Planetary Sciences/IGPP, University of California, Santa Cruz, David
More informationSEG/New Orleans 2006 Annual Meeting
and its implications for the curvelet design Hervé Chauris, Ecole des Mines de Paris Summary This paper is a first attempt towards the migration of seismic data in the curvelet domain for heterogeneous
More informationSummary. Introduction
Yaofeng He and Ru-Shan Wu Institute of Geophysics and Planetary Physics, University of California, Santa Cruz, C 9564 Summary We propose an approach to perform migration and imaging using secondary scattered
More informationImage Transformation Techniques Dr. Rajeev Srivastava Dept. of Computer Engineering, ITBHU, Varanasi
Image Transformation Techniques Dr. Rajeev Srivastava Dept. of Computer Engineering, ITBHU, Varanasi 1. Introduction The choice of a particular transform in a given application depends on the amount of
More informationFRACTAL IMAGE COMPRESSION OF GRAYSCALE AND RGB IMAGES USING DCT WITH QUADTREE DECOMPOSITION AND HUFFMAN CODING. Moheb R. Girgis and Mohammed M.
322 FRACTAL IMAGE COMPRESSION OF GRAYSCALE AND RGB IMAGES USING DCT WITH QUADTREE DECOMPOSITION AND HUFFMAN CODING Moheb R. Girgis and Mohammed M. Talaat Abstract: Fractal image compression (FIC) is a
More informationHonors Precalculus: Solving equations and inequalities graphically and algebraically. Page 1
Solving equations and inequalities graphically and algebraically 1. Plot points on the Cartesian coordinate plane. P.1 2. Represent data graphically using scatter plots, bar graphs, & line graphs. P.1
More informationDownloaded 02/17/16 to Redistribution subject to SEG license or copyright; see Terms of Use at
Downloaded 0/7/6 to 8.4.69.89. Redistribution subject to SEG license or copright; see Terms of Use at http://librar.seg.org/ 3D Beamlet migration in VTI medium Jian Mao *, Ru-Shan Wu, Universit of California
More informationAn Improvement in Temporal Resolution of Seismic Data Using Logarithmic Time-frequency Transform Method
Iranian Journal of Oil & Gas Science and Technology, Vol. 4 (2015), No. 2, pp. 27-39 http://ijogst.put.ac.ir An Improvement in Temporal Resolution of Seismic Data Using Logarithmic Time-frequency Transform
More information3D PRESTACK DEPTH MIGRATION WITH FACTORIZATION FOUR-WAY SPLITTING SCHEME
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume, Supp, Pages 183 196 c 5 Institute for Scientific Computing and Information 3D PRESTACK DEPTH MIGRATION WITH FACTORIZATION FOUR-WAY SPLITTING
More informationAnisotropy-preserving 5D interpolation by hybrid Fourier transform
Anisotropy-preserving 5D interpolation by hybrid Fourier transform Juefu Wang and Shaowu Wang, CGG Summary We present an anisotropy-preserving interpolation method based on a hybrid 5D Fourier transform,
More informationMain Menu. Summary. Introduction
Local-angle domain illumination for full-wave propagators Jun Cao* and Ru-Shan Wu Department of Earth and Planetary Sciences/IGPP, University of California, Santa Cruz Summary We propose an efficient split-step
More informationMastery. PRECALCULUS Student Learning Targets
PRECALCULUS Student Learning Targets Big Idea: Sequences and Series 1. I can describe a sequence as a function where the domain is the set of natural numbers. Connections (Pictures, Vocabulary, Definitions,
More information1D internal multiple prediction in a multidimensional world: errors and recommendations
1D internal multiple prediction 1D internal multiple prediction in a multidimensional world: errors and recommendations Pan Pan and Kris Innanen ABSTRACT Internal multiples are more difficult to estimate
More informationP052 3D Modeling of Acoustic Green's Function in Layered Media with Diffracting Edges
P052 3D Modeling of Acoustic Green's Function in Layered Media with Diffracting Edges M. Ayzenberg* (Norwegian University of Science and Technology, A. Aizenberg (Institute of Petroleum Geology and Geophysics,
More informationCompression of RADARSAT Data with Block Adaptive Wavelets Abstract: 1. Introduction
Compression of RADARSAT Data with Block Adaptive Wavelets Ian Cumming and Jing Wang Department of Electrical and Computer Engineering The University of British Columbia 2356 Main Mall, Vancouver, BC, Canada
More informationZ-99 3D Sub-salt Tomography Based on Wave Equation Migration Perturbation Scans
1 Z-99 3D Sub-salt Tomography Based on Wave Equation Migration Perturbation Scans BIN WANG 1, VOLKER DIRKS 1, PATRICE GUILLAUME 2, FRANÇOIS AUDEBERT 1, ANNING HOU 1 AND DURYODHAN EPILI 1 1 CGG Americas;
More informationImage Compression Algorithm for Different Wavelet Codes
Image Compression Algorithm for Different Wavelet Codes Tanveer Sultana Department of Information Technology Deccan college of Engineering and Technology, Hyderabad, Telangana, India. Abstract: - This
More information3D angle decomposition for elastic reverse time migration Yuting Duan & Paul Sava, Center for Wave Phenomena, Colorado School of Mines
3D angle decomposition for elastic reverse time migration Yuting Duan & Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY We propose 3D angle decomposition methods from elastic reverse
More informationCoherent partial stacking by offset continuation of 2-D prestack data
Stanford Exploration Project, Report 82, May 11, 2001, pages 1 124 Coherent partial stacking by offset continuation of 2-D prestack data Nizar Chemingui and Biondo Biondi 1 ABSTRACT Previously, we introduced
More informationWavelet based Keyframe Extraction Method from Motion Capture Data
Wavelet based Keyframe Extraction Method from Motion Capture Data Xin Wei * Kunio Kondo ** Kei Tateno* Toshihiro Konma*** Tetsuya Shimamura * *Saitama University, Toyo University of Technology, ***Shobi
More informationVertical and Horizontal Translations
SECTION 4.3 Vertical and Horizontal Translations Copyright Cengage Learning. All rights reserved. Learning Objectives 1 2 3 4 Find the vertical translation of a sine or cosine function. Find the horizontal
More informationOptimized Progressive Coding of Stereo Images Using Discrete Wavelet Transform
Optimized Progressive Coding of Stereo Images Using Discrete Wavelet Transform Torsten Palfner, Alexander Mali and Erika Müller Institute of Telecommunications and Information Technology, University of
More informationComparative Analysis of Image Compression Using Wavelet and Ridgelet Transform
Comparative Analysis of Image Compression Using Wavelet and Ridgelet Transform Thaarini.P 1, Thiyagarajan.J 2 PG Student, Department of EEE, K.S.R College of Engineering, Thiruchengode, Tamil Nadu, India
More informationCHAPTER 3 WAVELET DECOMPOSITION USING HAAR WAVELET
69 CHAPTER 3 WAVELET DECOMPOSITION USING HAAR WAVELET 3.1 WAVELET Wavelet as a subject is highly interdisciplinary and it draws in crucial ways on ideas from the outside world. The working of wavelet in
More informationMachine-learning Based Automated Fault Detection in Seismic Traces
Machine-learning Based Automated Fault Detection in Seismic Traces Chiyuan Zhang and Charlie Frogner (MIT), Mauricio Araya-Polo and Detlef Hohl (Shell International E & P Inc.) June 9, 24 Introduction
More informationSEG/San Antonio 2007 Annual Meeting
Yaofeng He* and Ru-Shan Wu Modeling and Imaging Laboratory, Institute of Geophysics and Planetary Physics, University of California, Santa Cruz, CA, 95064, USA Summary A one-way and one-return boundary
More informationCHAPTER 3 DIFFERENT DOMAINS OF WATERMARKING. domain. In spatial domain the watermark bits directly added to the pixels of the cover
38 CHAPTER 3 DIFFERENT DOMAINS OF WATERMARKING Digital image watermarking can be done in both spatial domain and transform domain. In spatial domain the watermark bits directly added to the pixels of the
More informationInversion after depth imaging
Robin P. Fletcher *, Stewart Archer, Dave Nichols, and Weijian Mao, WesternGeco Summary In many areas, depth imaging of seismic data is required to construct an accurate view of the reservoir structure.
More informationRobustness of the scalar elastic imaging condition for converted waves
CWP-830 Robustness of the scalar elastic imaging condition for converted waves Yuting Duan & Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT For elastic reverse-time migration, one
More informationSUMMARY ELASTIC SCALAR IMAGING CONDITION
Robust 3D scalar imaging condition for elastic RTM Yuting Duan, presently at Shell International Exploration and Production Inc., formerly at Center for Wave Phenomena, Colorado School of Mines Paul Sava,
More informationComputational Geometry
Windowing queries Windowing Windowing queries Zoom in; re-center and zoom in; select by outlining Windowing Windowing queries Windowing Windowing queries Given a set of n axis-parallel line segments, preprocess
More informationTexture. Outline. Image representations: spatial and frequency Fourier transform Frequency filtering Oriented pyramids Texture representation
Texture Outline Image representations: spatial and frequency Fourier transform Frequency filtering Oriented pyramids Texture representation 1 Image Representation The standard basis for images is the set
More informationKeywords - DWT, Lifting Scheme, DWT Processor.
Lifting Based 2D DWT Processor for Image Compression A. F. Mulla, Dr.R. S. Patil aieshamulla@yahoo.com Abstract - Digital images play an important role both in daily life applications as well as in areas
More informationWriting Kirchhoff migration/modelling in a matrix form
Writing Kirchhoff migration/modelling in a matrix form Abdolnaser Yousefzadeh and John C. Bancroft Kirchhoff migration matrix ABSTRACT Kirchhoff prestack migration and modelling are linear operators. Therefore,
More informationA Nuclear Norm Minimization Algorithm with Application to Five Dimensional (5D) Seismic Data Recovery
A Nuclear Norm Minimization Algorithm with Application to Five Dimensional (5D) Seismic Data Recovery Summary N. Kreimer, A. Stanton and M. D. Sacchi, University of Alberta, Edmonton, Canada kreimer@ualberta.ca
More informationClosing the Loop via Scenario Modeling in a Time-Lapse Study of an EOR Target in Oman
Closing the Loop via Scenario Modeling in a Time-Lapse Study of an EOR Target in Oman Tania Mukherjee *(University of Houston), Kurang Mehta, Jorge Lopez (Shell International Exploration and Production
More informationEdge Equalized Treemaps
Edge Equalized Treemaps Aimi Kobayashi Department of Computer Science University of Tsukuba Ibaraki, Japan kobayashi@iplab.cs.tsukuba.ac.jp Kazuo Misue Faculty of Engineering, Information and Systems University
More informationStatics: the abusive power of trimming
John C. Bancroft, Alan Richards, and Charles P. Ursenbach Statics: the abusive power of trimming ABSTRACT The application of trim statics directly to each seismic trace can be very dangerous as any seismic
More informationImaging diffraction points using the local image matrices generated in prestack migration
GEOPHYSICS, VOL. 75, NO. 1 JANUARY-FEBRUARY 1 ; P. S1 S9, 16 FIGS. 1.119/1.3775 Downloaded 8/3/16 to 18.114.69.87. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/
More informationEquivalence of source-receiver migration and shot-profile migration
Stanford Exploration Project, Report 112, November 11, 2002, pages 109 117 Short Note Equivalence of source-receiver migration and shot-profile migration Biondo Biondi 1 INTRODUCTION At first glance, shot
More informationIntroduction to Wavelets
Lab 11 Introduction to Wavelets Lab Objective: In the context of Fourier analysis, one seeks to represent a function as a sum of sinusoids. A drawback to this approach is that the Fourier transform only
More informationWavelet Transform (WT) & JPEG-2000
Chapter 8 Wavelet Transform (WT) & JPEG-2000 8.1 A Review of WT 8.1.1 Wave vs. Wavelet [castleman] 1 0-1 -2-3 -4-5 -6-7 -8 0 100 200 300 400 500 600 Figure 8.1 Sinusoidal waves (top two) and wavelets (bottom
More informationRecovering the Reflectivity Matrix and Angledependent Plane-wave Reflection Coefficients from Imaging of Multiples
Recovering the Reflectivity Matrix and Angledependent Plane-wave Reflection Coefficients from Imaging of Multiples A. Ordoñez* (PGS), W.F. Sollner (PGS), T. Klüver (PGS) & L.G. Gelius (UiO) SUMMARY A joint
More informationReview of seismic imaging: Prestack
Review of seismic imaging: Prestack John C. Bancroft Published in the CSEG Recorder - November 001 Review of seismic migration: Prestack ABSTRACT The kinematics of prestack data considers an arbitrary
More informationA stochastic approach of Residual Move Out Analysis in seismic data processing
A stochastic approach of Residual ove Out Analysis in seismic data processing JOHNG-AY T.,, BORDES L., DOSSOU-GBÉTÉ S. and LANDA E. Laboratoire de athématique et leurs Applications PAU Applied Geophysical
More informationSummary. Introduction. Source separation method
Separability of simultaneous source data based on distribution of firing times Jinkun Cheng, Department of Physics, University of Alberta, jinkun@ualberta.ca Summary Simultaneous source acquisition has
More informationa) It obeys the admissibility condition which is given as C ψ = ψ (ω)
Chapter 2 Introduction to Wavelets Wavelets were shown in 1987 to be the foundation of a powerful new approach to signal processing and analysis called multi-resolution theory by S. Mallat. As its name
More informationAVA analysis based on RTM angle-domain common image gather Rui Yan* and Xiao-Bi Xie, IGPP, University of California, Santa Cruz
AVA analysis based on RTM angle-domain common image gather Rui Yan* and Xiao-Bi Xie, IGPP, University of California, Santa Cruz Summary We propose an alternative approach of AVA analysis based on RTM angle-domain
More informationSeismic Modeling, Migration and Velocity Inversion
Seismic Modeling, Migration and Velocity Inversion Aliasing Bee Bednar Panorama Technologies, Inc. 14811 St Marys Lane, Suite 150 Houston TX 77079 May 16, 2014 Bee Bednar (Panorama Technologies) Seismic
More informationPractical implementation of SRME for land multiple attenuation
Practical implementation of SRME for land multiple attenuation Juefu Wang* and Shaowu Wang, CGGVeritas, Calgary, Canada juefu.wang@cggveritas.com Summary We present a practical implementation of Surface
More informationVelocity analysis using surface-seismic primaries-only data obtained without removing multiples
Delft University of Technology Velocity analysis using surface-seismic primaries-only data obtained without removing multiples Dokter, E.; Meles, Giovanni; Curtis, A; Wapenaar, Kees DOI.3997/4-469.7 Publication
More informationAngle-gather time migration a
Angle-gather time migration a a Published in SEP report, 1, 141-15 (1999) Sergey Fomel and Marie Prucha 1 ABSTRACT Angle-gather migration creates seismic images for different reflection angles at the reflector.
More informationImprovements in time domain FWI and its applications Kwangjin Yoon*, Sang Suh, James Cai and Bin Wang, TGS
Downloaded 0/7/13 to 05.196.179.38. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/ Improvements in time domain FWI and its applications Kwangjin Yoon*,
More informationA low rank based seismic data interpolation via frequencypatches transform and low rank space projection
A low rank based seismic data interpolation via frequencypatches transform and low rank space projection Zhengsheng Yao, Mike Galbraith and Randy Kolesar Schlumberger Summary We propose a new algorithm
More informationk are within the range 0 k < k < π dx,
Acoustic Wave Equation by Implicit Spatial Derivative Operators Dan Kosloff*, Department of Geophysics, Tel-Aviv University and Paradigm, Reynam Pestana, Department of Geophysics, Federal University of
More informationSummary. Figure 1: Simplified CRS-based imaging workflow. This paper deals with the boxes highlighted in green.
Smoothing and automated picking of kinematic wavefield attributes Tilman Klüver and Jürgen Mann, Geophysical Institute, University of Karlsruhe, Germany Copyright 2005, SBGf Sociedade Brasiliera de Geofísica
More informationReflectivity modeling for stratified anelastic media
Reflectivity modeling for stratified anelastic media Peng Cheng and Gary F. Margrave Reflectivity modeling ABSTRACT The reflectivity method is widely used for the computation of synthetic seismograms for
More informationThe azimuth-dependent offset-midpoint traveltime pyramid in 3D HTI media
The azimuth-dependent offset-midpoint traeltime pyramid in 3D HTI media Item Type Conference Paper Authors Hao Qi; Stoas Alexey; Alkhalifah Tariq Ali Eprint ersion Publisher's Version/PDF DOI.9/segam3-58.
More informationCrosswell Imaging by 2-D Prestack Wavepath Migration
Crosswell Imaging by 2-D Prestack Wavepath Migration Hongchuan Sun ABSTRACT Prestack wavepath migration (WM) is applied to 2-D synthetic crosswell data, and the migrated images are compared to those from
More informationAn imaging technique for subsurface faults using Teleseismic-Wave Records II Improvement in the detectability of subsurface faults
Earth Planets Space, 52, 3 11, 2000 An imaging technique for subsurface faults using Teleseismic-Wave Records II Improvement in the detectability of subsurface faults Takumi Murakoshi 1, Hiroshi Takenaka
More informationTrue Amplitude Extraction Technology and Its Application in Pre-Stack Seismic Data
Journal of Earth Science, Vol. 26, No. 4, p. 522 529, August 25 ISSN 674-487X Printed in China DOI:.7/s2583-5-567-2 True Amplitude Extraction Technology and Its Application in Pre-Stack Seismic Data Xiaohui
More informationSection Graphs of the Sine and Cosine Functions
Section 5. - Graphs of the Sine and Cosine Functions In this section, we will graph the basic sine function and the basic cosine function and then graph other sine and cosine functions using transformations.
More informationIt is widely considered that, in regions with significant
Multifocusing-based multiple attenuation Alex Berkovitch and Kostya Deev, Geomage Evgeny Landa, OPERA It is widely considered that, in regions with significant geologic complexity, methods which work directly
More informationFrom Fourier Transform to Wavelets
From Fourier Transform to Wavelets Otto Seppälä April . TRANSFORMS.. BASIS FUNCTIONS... SOME POSSIBLE BASIS FUNCTION CONDITIONS... Orthogonality... Redundancy...3. Compact support.. FOURIER TRANSFORMS
More informationSeismic modelling with the reflectivity method
Yongwang Ma, Luiz Loures, and Gary F. Margrave ABSTRACT Seismic modelling with the reflectivity method Numerical seismic modelling is a powerful tool in seismic imaging, interpretation and inversion. Wave
More informationWavelet Transform On Propagation and Imaging For Seismic Exploration
Proposal for WTOPI Research Consortium Phase IV Wavelet Transform On Propagation and Imaging For Seismic Exploration Phase IV (4/1/2006-3/31/2009): Imaging and Inversion with Beamlets: Processing in Local
More informationA Short Narrative on the Scope of Work Involved in Data Conditioning and Seismic Reservoir Characterization
A Short Narrative on the Scope of Work Involved in Data Conditioning and Seismic Reservoir Characterization March 18, 1999 M. Turhan (Tury) Taner, Ph.D. Chief Geophysicist Rock Solid Images 2600 South
More informationA Wavelet Tour of Signal Processing The Sparse Way
A Wavelet Tour of Signal Processing The Sparse Way Stephane Mallat with contributions from Gabriel Peyre AMSTERDAM BOSTON HEIDELBERG LONDON NEWYORK OXFORD PARIS SAN DIEGO SAN FRANCISCO SINGAPORE SYDNEY»TOKYO
More informationSeparate CT-Reconstruction for Orientation and Position Adaptive Wavelet Denoising
Separate CT-Reconstruction for Orientation and Position Adaptive Wavelet Denoising Anja Borsdorf 1,, Rainer Raupach, Joachim Hornegger 1 1 Chair for Pattern Recognition, Friedrich-Alexander-University
More informationEffects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines
Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY Velocity models used for wavefield-based seismic
More informationApplication of Daubechies Wavelets for Image Compression
Application of Daubechies Wavelets for Image Compression Heydari. Aghile 1,*, Naseri.Roghaye 2 1 Department of Math., Payame Noor University, Mashad, IRAN, Email Address a_heidari@pnu.ac.ir, Funded by
More informationDownloaded 09/09/15 to Redistribution subject to SEG license or copyright; see Terms of Use at
Recovering the Reflectivity Matrix and Angle-dependent Plane-wave Reflection Coefficients from Imaging of Multiples Alba Ordoñez PGS/UiO*, Walter Söllner PGS, Tilman Klüver PGS and Leiv J. Gelius UiO Summary
More informationDigital Image Processing. Chapter 7: Wavelets and Multiresolution Processing ( )
Digital Image Processing Chapter 7: Wavelets and Multiresolution Processing (7.4 7.6) 7.4 Fast Wavelet Transform Fast wavelet transform (FWT) = Mallat s herringbone algorithm Mallat, S. [1989a]. "A Theory
More informationSeismic data Interpolation in the Continuous Wavenumber Domain, Flexibility and Accuracy
Seismic data Interpolation in the Continuous Wavenumber Domain, Flexibility and Accuracy Ye Zheng Geo-X Exploration Services Inc. Summary Our recently developed algorithm, ASFT (Arbitrarily Sampled Fourier
More informationShort Note. DMO velocity analysis with Jacubowicz s dip-decomposition method. David Kessler and Wai-Kin Chan*
GEOPHYSICS, VOL. 58, NO. 10 (OCTOBER 1993); P. 1517-1524,9 FIGS. Short Note DMO velocity analysis with Jacubowicz s dip-decomposition method David Kessler and Wai-Kin Chan* INTRODUCTION Dip-moveout (DMO)
More informationA hybrid elastic one-way propagator for strong-contrast media and its application to subsalt migration
A hybrid elastic one-way propagator for strong-contrast media and its application to subsalt migration Rui Yan 1, Ru-Shan Wu 2, Xiao-Bi Xie 3, and David Walraven 4 1 University of California, Santa Cruz,
More informationPoint Lattices in Computer Graphics and Visualization how signal processing may help computer graphics
Point Lattices in Computer Graphics and Visualization how signal processing may help computer graphics Dimitri Van De Ville Ecole Polytechnique Fédérale de Lausanne Biomedical Imaging Group dimitri.vandeville@epfl.ch
More informationWave-equation migration from topography: Imaging Husky
Stanford Exploration Project, Report 123, October 31, 2005, pages 49 56 Short Note Wave-equation migration from topography: Imaging Husky Jeff Shragge 1 INTRODUCTION Imaging land seismic data is wrought
More informationFUNCTIONS AND MODELS
1 FUNCTIONS AND MODELS FUNCTIONS AND MODELS 1.3 New Functions from Old Functions In this section, we will learn: How to obtain new functions from old functions and how to combine pairs of functions. NEW
More informationInterpolation using asymptote and apex shifted hyperbolic Radon transform
Interpolation using asymptote and apex shifted hyperbolic Radon transform Amr Ibrahim, Paolo Terenghi and Mauricio D. Sacchi Department of Physics, University of Alberta PGS Abstract The asymptote and
More informationThe Conventional and Non-conventional Seismic Differencing
Summary The Conventional and Non-conventional Seismic Differencing Vanja Vracar* University of Calgary, CREWES Project, Alberta, Canada vmilicev@ucalgary.ca and Robert J. Ferguson University of Calgary,
More informationOptimised corrections for finite-difference modelling in two dimensions
Optimized corrections for 2D FD modelling Optimised corrections for finite-difference modelling in two dimensions Peter M. Manning and Gary F. Margrave ABSTRACT Finite-difference two-dimensional correction
More informationIMAGE DE-NOISING IN WAVELET DOMAIN
IMAGE DE-NOISING IN WAVELET DOMAIN Aaditya Verma a, Shrey Agarwal a a Department of Civil Engineering, Indian Institute of Technology, Kanpur, India - (aaditya, ashrey)@iitk.ac.in KEY WORDS: Wavelets,
More informationAN ALGORITHM FOR BLIND RESTORATION OF BLURRED AND NOISY IMAGES
AN ALGORITHM FOR BLIND RESTORATION OF BLURRED AND NOISY IMAGES Nader Moayeri and Konstantinos Konstantinides Hewlett-Packard Laboratories 1501 Page Mill Road Palo Alto, CA 94304-1120 moayeri,konstant@hpl.hp.com
More informationUnit 1 Quadratic Functions
Unit 1 Quadratic Functions This unit extends the study of quadratic functions to include in-depth analysis of general quadratic functions in both the standard form f ( x) = ax + bx + c and in the vertex
More informationContourlet Based Image Recovery and De-noising Through Wireless Fading Channels
2005 Conference on Information Science and Systems, The Johns Hopkins University, March 16-18, 2005 Contourlet Based Image Recovery and De-noising Through Wireless Fading Channels Yanun Yan and Lisa Osadciw
More informationsignal-to-noise ratio (PSNR), 2
u m " The Integration in Optics, Mechanics, and Electronics of Digital Versatile Disc Systems (1/3) ---(IV) Digital Video and Audio Signal Processing ƒf NSC87-2218-E-009-036 86 8 1 --- 87 7 31 p m o This
More informationCommon-angle processing using reflection angle computed by kinematic pre-stack time demigration
Common-angle processing using reflection angle computed by kinematic pre-stack time demigration Didier Lecerf*, Philippe Herrmann, Gilles Lambaré, Jean-Paul Tourré and Sylvian Legleut, CGGVeritas Summary
More informationFast Implementation of VC-1 with Modified Motion Estimation and Adaptive Block Transform
Circuits and Systems, 2010, 1, 12-17 doi:10.4236/cs.2010.11003 Published Online July 2010 (http://www.scirp.org/journal/cs) Fast Implementation of VC-1 with Modified Motion Estimation and Adaptive Block
More informationSeismic Modeling, Migration and Velocity Inversion
Seismic Modeling, Migration and Velocity Inversion Full Waveform Inversion Bee Bednar Panorama Technologies, Inc. 14811 St Marys Lane, Suite 150 Houston TX 77079 May 18, 2014 Bee Bednar (Panorama Technologies)
More informationSeparation of diffracted waves in TI media
CWP-829 Separation of diffracted waves in TI media Yogesh Arora & Ilya Tsvankin Center for Wave Phenomena, Colorado School of Mines Key words: Diffractions, Kirchhoff, Specularity, Anisotropy ABSTRACT
More information