Seismic trace interpolation with approximate message passing Navid Ghadermarzy and Felix Herrmann and Özgür Yılmaz, University of British Columbia
|
|
- Mary Harper
- 5 years ago
- Views:
Transcription
1 Seismic trace interpolation with approximate message passing Navid Ghadermarzy and Felix Herrmann and Özgür Yılmaz, University of British Colmbia SUMMARY Approximate message passing (AMP) is a comptationally effective algorithm for recovering high dimensional signals from a few compressed measrements. In this paper we se AMP to solve the seismic trace interpolation problem. We also show that we can exploit the fast AMP algorithm to improve the recovery reslts of seismic trace interpolation in crvelet domain, both in terms of convergence speed and recovery performance by sing AMP in Forier domain as a preprocessor for the ` recovery in Crvelet domain. INTRODUCTION Reglar sampling along the time axis is common approach for seismic imaging. However, seismic data are often spatially ndersampled de to economical reasons as well as grond srface limitations. Too large spatial sampling interval may reslt in spatial aliasing. Operational conditions might also reslt is noisy traces or even gaps in coverage and irreglar sampling which often leads to image artifacts. Seismic trace interpolation aims to interpolate missing traces in an otherwise reglarly sampled seismic data to a complete reglar grid. Compressed sensing (CS) is a data acqisition techniqe to efficiently acqiring and recovering high-dimensional signals from seemingly incomplete linear measrements, provided that the signal is sparse or compressible in some transform domain (Donoho (006); Candès et al. (006)). Several works have shown that the strctre of seismic wavefronts can be captred by a few significant transform coefficients (e.g., crvelet or Forier coefficients) and the problem of seismic trace interpolation can be cast as a sparse recovery problem once the data is transformed into a sitable domain, e.g., discrete Forier transform (DFT) (Sacchi and Ulrych (996); Li and Sacchi (004)), Radon transform (Thorson and Claerbot (985); Hampson et al. (986); Herrmann et al. (000)), or crvelet transform (Hennenfent and Herrmann (005, 006, 008); Naghizadeh et al. (00)). Hence, the interpolation problem becomes that of finding transform domain coefficients with the smallest ` norm that satisfy the sbsampled measrements. In this paper we first consider sing crvelet transform as the sparsifying transform (see, e.g., Hennenfent and Herrmann (006)). Crvelet transform is a redndant transform which expands the size of the model and has been shown to be very effective in captring the strctre of the data with very few significant coefficients. Then, we recover the reslting compressible signal by solving a convex minimization problem. Next we consider sing Forier transform as the transform operator. The reslting transform is not as effective as crvelet transform (in terms of sparsity of the reslting transform). However, the transformed data can still be approximated by a solving a sparse recovery problem which is mch faster than the crvelet case. This is de to two reasons. First one is the availability of very fast DFT transforms and second is the applicability of approximate message passing (AMP) for this problem. AMP (Donoho et al. (009)) is a fast (needs few matrix-vector mltiplications) iterative ` recovery algorithm which can be sed to recover sparse Forier coefficients (The crrent variants of AMP can not recover crvelet coefficients). Also notice that AMP can be sed to speed p interpolation techniqes that work on small high dimensional windows. In this paper we se AMP with Forier transform coefficients to generate relatively good approximations of the seismic data which is sed to enhance the speed and recovery performance of seismic trace interpolation in the crvelet domain. First we formlate the problem of seismic trace interpolation as a sparse recovery problem. Next we explain the approach we take to solve this problem both from the crvelet transform and Forier transform coefficients and then we present the -stage algorithm which combines these two approaches. PROBLEM FORMULATION Consider a seismic line with N s sorces located on earth srface which send sond waves into the earth and N r receivers record the reflection in N t time samples. Rearranging the seismic line, we have a signal f R N, where N = N sn rn t. Assme x = Sf where x is the sparse representation of f in some transform domain. We want to recover a very high dimensional seismic data volme x by interpolating between a smaller nmber of measrements b = RMf, where RM is a sampling operator composed of the prodct of a restriction matrix R, which specifies the sample locations that have data and a measrement basis matrix M. To overcome the problem of high dimensionality, recovery is performed by first partitioning the seismic data volmes into freqency slices and then solving a seqence of individal sbproblems (Mansor et al. (03)). For each partition j, let b (j) be the sbsampled measrements of data f (j) in partition j and the corresponding measrement matrix is given by A (j) = R (j) S H (Notice that the sperscript H denoted the Hermitian transpose and for simplicity, we have as-
2 smed the measrement matrix M to be the identity). Then for each partition we need to estimate the soltion of the following minimization problem independently: ex (j) := arg min kk 0 s.t. A (j) = b (j) ; () and then f (j) is approximated by S H ex (j). Recovery in crvelet domain with convex optimization In this case S C C P N is the crvelet transform which is sed as the sparsifying transform S. For each partition we se the -D crvelet transform where N P 8 which is highly redndant and allows for sparse representations of the data in the transformed domain. Notice that in this case A (j) = R (j) SC H. The minimization problem () is a combinatorial problem and therefore we approximate the soltion by the following convex relaxation: ex (j) ` := arg min kk s.t. A (j) = b (j) ; () and then f (j) is approximated by SC H ex (j). Here kk ` = NP j i j and () can be efficiently solved by algorithms i= like SPGL (Van Den Berg and Friedlander (007)). Mansor et al. (03) noted that althogh partitioning the data sing and -D crvelets redces the size of the problem, it doesn t tilize the continity across partitions. Therefore they sggest sing spport set of each partition as a spport estimate for the next one. For partition j > assme T j := spp(s C SC H ex jk) which is the locations of the largest k entries of S C SC H ex in magnitde and define the weight vector w R N sch that wj = for j = T j and wj =! < for j T j, where k is chosen to be the nmber of largest transform-domain coefficients that contribte 95% of the signal energy. Mansor et al. (0) prove that solving the following weighted ` minimization reslts in better recovery compared to () if the spport estimate T j is at least 50% accrate: ex (j) := argmin kk ;wj s.t. A (j) = b (j) ; (3) w` where kk ;wj = NP wij i j is the weighted ` norm of. i= Recovery in Forier domain with AMP In this case we se DFT matrix S F C N N as the sparsifying operator. Hence, A (j) = R (j) S F is the sbsampled -D DFT matrix. To solve () the AMP algorithm (Donoho et al. (009)) starts from an initial x 0 amp and an initial threshold ^ 0 = and iteratively goes by x t+ amp = (x t amp + A (j)h z t ; t ); z t = y A (j) x t amp + z t h 0 (x t amp + A (j)h z t ; t )i; (4) and ex (j) tf inal amp = xamp where t final is the nmber of times we iterate the algorithm. Here SF H xt amp is the crrent estimate of f (j), z t is the crrent residal, h i is the arithmetic mean, and is the soft thresholding fnction [(a;s)]j = sign(a j )(a j s) +. This fnction acts on the coefficients of the signal a and zeroes ot any vale which is less than the scalar s in magnitde. Similar to the previos case we can tilize the continity of seismic data across partitions to estimate the spport set of each partition. As before T j := spp(s F SF H x jk) and the weight vector w is defined sch that wj = for j = T j and wj =! < for j T j. Ghadermarzy and Yilmaz (03) incorporate this prior spport information into the AMP algorithm by the following iterations: x t+ wamp = (x t wamp + A (j)h z t ; t w j ); z t = y A (j) x t wamp + z t h 0 (x t wamp + A (j)h z t ; t w j )i: (5) The only difference between AMP and weighted AMP (WAMP) is in the vale of the thresholds which is set to be a smaller vale for the coefficients in the spport estimate which makes them more likely to remain nonzero. In Figres.a-f we compare the performance of AMP and WAMP (in Forier domain) with weighted ` (in crvelet domain) on recovering a seismic line from the Glf of Sez with 50% randomly sbsampled receivers. The seismic line at fll resoltion has N s = 78 sorces, N r = 78 receivers, and Nt = time samples acqired with a sampling interval of 4 ms. Conseqently, the seismic line contains samples collected in a s temporal window with a maximm freqency of 5 Hz. To access the freqency slices we take the -D DFT of the data along time axis and interpolate the missing receivers by solving the seismic trace interpolation with weighted `, AMP and weighted AMP algorithms explained above (Reslts of ` minimization is omitted as the weighted ` minimization is always giving better reslts (Mansor et al. (03)). Notice that reconstrction reslts of weighted AMP is better than those of AMP. -STAGE ALGORITHM As mentioned before AMP algorithm is a significantly fast algorithm and Figre shows that althogh the reslts of weighted ` is better than those of weighted AMP, the reslts of weighted AMP is still reasonably good (considering that weighted AMP is abot 40 times faster than the weighted ` algorithm). Weighted ` minimization in the crvelet domain gives mch better reslts which is de to the better choice of sparsifying transform domain. On the other hand AMP is a very fast algorithm that approximates the tre data reasonably good. Notice that both these algorithms are solving the same problem in different transform domains.
3 Weighted L error image in SR 0. Weighted L minimization in SR (a) Recovery by weighted ` in crvelet domain AMP error image in SR (c) Recovery by AMP in Forier domain (d) AMP error Weighted AMP error image in SR Weighted AMP in SR (b) Weighted ` error AMP in SR (e) Recovery by weighted AMP in Forier domain (f) WAMP error WAMP+W L minimization in SR WAMP+W L error image in SR (g) -stage reconstrction (h) -stage error Figre : Comparison of the different recovery methods. Figres show recovered shotgather nmber 84 from sbsampled seismic line from the Glf of Sez. Figres on the left show the reconstrctions and figres on the right show the recovery errors for weighted ` (a, b), AMP (c, d), WAMP (e,f), and the -stage algorithm (g, h). f is approximated by: Therefore, we can combine these algorithms to se both the fast convergence of AMP algorithm and sperior reslts of weighted `. Notice that at each partition j, e e H weighted ` : f SC xw` weighted AM P 3 : f S H F xwamp (6)
4 Hence the otcome of the weighted ` algorithm can be approximated by the final otcome of the weighted AMP algorithm (which is in Forier domain) throgh: S H C ex (j) SF H w` ex (j) wamp: (7) Hence, we propose a -stage algorithm, that for each partition j first we apply the weighted AMP algorithm (5) in the Forier domain to get the estimate SC H ex (j) w`. Next we se the largest coefficients of S C SF H ex (j) wamp as a spport estimate for ex (j) and solve a weighted ` w` algorithm in the crvelet domain. Figre shows an example of the accracy achieved to predict the spport of the crvelet coefficients in a seismic line when we jst se the adjacent freqency slices (the weighted ` case) and when we combine the information of the previos freqency slice with the approximation gained by WAMP. Here the spport is chosen to be the set of largest transform-domain coefficients that contribte 95% of the signal energy. The spport estimate gained by the -stage algorithm is always better than the estimate gained by the previos slice. We also feed the approximation S C SF H ex (j) wamp into the weighted ` algorithm as a good initial point for the SPGL algorithm. As the crvelet transform is a redndant transform, we can not mathematically jstify why this choice of initial point improves the convergence rate. However, empirical reslts show significant improvement when we feed the reslt of WAMP into the ` algorithm. Figres.gh show the recovery reslts by the -stage algorithm. In Table we have compared the nmber of matrix-vector mltiplications and average time it takes for the above algorithms to recover different freqency slices. The - stage algorithm presented ses fewer crvelet transforms compared to the weighted ` which makes it faster. Figre 3 shows the SNRs achieved by `, weighted `, AMP, Weighted AMP, and the -stage algorithms for different shot gathers. The reslts of the -satge algorithm is always better than the weighted `. recovery method Comparison of different methods # of DFT # of crvelet transforms transforms average rntime weighted ` s AMP s WAMP s -stage WAMP+w` s Table : Comparison of approximate average comptational complexity for different recovery methods in Figre averaged over freqency slices. Fraction of accrate detected spport 0. -stage Weighted l Freqency (Hz) Figre : Example of the accracy achieved in estimating the spport of a freqency slice by sing the -stage algorithm and the weighted ` algorithm. SNR Shot gather nmber W L on freq. in SR L on freq. in SR AMP on freq. in SR Figre 3: Comparison of the SNRs achieved by different methods explained above in recovering shot gathers from the Sez seismic line in the sorce-receiver domain. CONCLUSIONS In this abstract we showed how we can se the approximate message passing algorithm for the problem of seismic trace interpolation. De to fast natre of AMP algorithm and the DFT transform we can get a good approximation of the seismic shot gathers mch faster than doing the interpolation in the crvelet domain. Hence, AMP can be sed as a fast preprocessor for interpolation in crvelet domain. This reslt shows that doing this not only improves the recovery reslts significantly and bt also the -stage algorithms need less SPGL iterations (with crvelets) to get to a good approximation once the fast WAMP algorithms is done in the Forier domain. Weighted AMP on freq. in SR Two stage WAMP+Wspgl 4
5 REFERENCES Candès, E. J., J. Romberg, and T. Tao, 006, Stable signal recovery from incomplete and inaccrate measrements: Commnications on Pre and Applied Mathematics, 59, Donoho, D., 006, Compressed sensing.: IEEE Transactions on Information Theory, 5, Donoho, D. L., A. Maleki, and A. Montanari, 009, Message-passing algorithms for compressed sensing: Proceedings of the National Academy of Sciences, 06, Ghadermarzy, N., and O. Yilmaz, 03, Weighted and reweighted approximate message passing: SPIE Optical Engineering+ Applications, International Society for Optics and Photonics, 88580C 88580C. Hampson, D., et al., 986, Inverse velocity stacking for mltiple elimination: Presented at the 986 SEG Annal Meeting, Society of Exploration Geophysicists. Hennenfent, G., and F. J. Herrmann, 005, Sparsenessconstrained data contination with frames: Applications to missing traces and aliased signals in /3-d: Presented at the SEG International Exposition and 75th Annal Meeting., 006, Seismic denoising with nonniformly sampled crvelets: Compting in Science & Engineering, 8, 6 5., 008, Non-parametric seismic data recovery with crvelet frames: Geophysical Jornal International, 73, Herrmann, P., T. Mojesky, M. Magesan, P. Hgonnet, et al., 000, De-aliased, high-resoltion radon transforms: 70th Annal International Meeting, SEG, Li, B., and M. D. Sacchi, 004, Minimm weighted norm interpolation of seismic records: Geophysics, 69, Mansor, H., M. P. Friedlander, R. Saab, and O. Yılmaz, 0, Recovering compressively sampled signals sing partial spport information: IEEE Transactions on Information Theory, 58, 34. Mansor, H., F. J. Herrmann, and O. Yılmaz, 03, Improved wavefield reconstrction from randomized sampling via weighted one-norm minimization: Geophysics, 78, no. 5, V93 V06. Naghizadeh, M., M. Sacchi, et al., 00, Hierarchical scale crvelet interpolation of aliased seismic data: Presented at the SEG Inter-national Exposition and 80th Annal Meeting.[S..]: SEG. Sacchi, M. D., and T. J. Ulrych, 996, Estimation of the discrete forier transform, a linear inversion approach: Geophysics, 6, Thorson, J. R., and J. F. Claerbot, 985, Velocitystack and slant-stack stochastic inversion: Geophysics, 50, Van Den Berg, E., and M. Friedlander, 007, Spgl: A solver for large-scale sparse reconstrction: Online: cs. bc. ca/labs/scl/spgl. 5
Tu P7 15 First-arrival Traveltime Tomography with Modified Total Variation Regularization
T P7 15 First-arrival Traveltime Tomography with Modified Total Variation Reglarization W. Jiang* (University of Science and Technology of China) & J. Zhang (University of Science and Technology of China)
More informationOptimal Sampling in Compressed Sensing
Optimal Sampling in Compressed Sensing Joyita Dtta Introdction Compressed sensing allows s to recover objects reasonably well from highly ndersampled data, in spite of violating the Nyqist criterion. In
More informationRandomized sampling strategies
Randomized sampling strategies Felix J. Herrmann SLIM Seismic Laboratory for Imaging and Modeling the University of British Columbia SLIM Drivers & Acquisition costs impediments Full-waveform inversion
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 informationSUMMARY. In combination with compressive sensing, a successful reconstruction
Higher dimensional blue-noise sampling schemes for curvelet-based seismic data recovery Gang Tang, Tsinghua University & UBC-Seismic Laboratory for Imaging and Modeling (UBC-SLIM), Reza Shahidi, UBC-SLIM,
More informationA sufficient condition for spiral cone beam long object imaging via backprojection
A sfficient condition for spiral cone beam long object imaging via backprojection K. C. Tam Siemens Corporate Research, Inc., Princeton, NJ, USA Abstract The response of a point object in cone beam spiral
More informationImage Compression Compression Fundamentals
Compression Fndamentals Data compression refers to the process of redcing the amont of data reqired to represent given qantity of information. Note that data and information are not the same. Data refers
More informationDeconvolution with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC
Deconvolution with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC SUMMARY We use the recently introduced multiscale and multidirectional curvelet transform to exploit the
More informationBias of Higher Order Predictive Interpolation for Sub-pixel Registration
Bias of Higher Order Predictive Interpolation for Sb-pixel Registration Donald G Bailey Institte of Information Sciences and Technology Massey University Palmerston North, New Zealand D.G.Bailey@massey.ac.nz
More informationDIRECT AND PROGRESSIVE RECONSTRUCTION OF DUAL PHOTOGRAPHY IMAGES
DIRECT AND PROGRESSIVE RECONSTRUCTION OF DUAL PHOTOGRAPHY IMAGES Binh-Son Ha 1 Imari Sato 2 Kok-Lim Low 1 1 National University of Singapore 2 National Institte of Informatics, Tokyo, Japan ABSTRACT Dal
More informationSUMMARY. In this paper, we present a methodology to improve trace interpolation
Reconstruction of seismic wavefields via low-rank matrix factorization in the hierarchical-separable matrix representation Rajiv Kumar, Hassan Mansour, Aleksandr Y. Aravkin, and Felix J. Herrmann University
More informationp-norm MINIMIZATION OVER INTERSECTIONS OF CONVEX SETS İlker Bayram
p-norm MINIMIZATION OVER INTERSECTIONS OF CONVEX SETS İlker Bayram Istanbl Technical University, Department of Electronics and Telecommnications Engineering, Istanbl, Trkey ABSTRACT We consider the imization
More informationSeismic Data Interpolation With Symmetry
Seismic Data Interpolation With Symmetry James Johnson joint work with Gilles Hennenfent 8 Consortium Meeting Outline Introduction Reciprocity Real Data Example Data Interpolation Problem Statement CRSI
More informationTime-jittered ocean bottom seismic acquisition Haneet Wason and Felix J. Herrmann
Time-jittered ocean bottom seismic acquisition Haneet Wason and Felix J. Herrmann SLIM University of British Columbia Challenges Need for full sampling - wave-equation based inversion (RTM & FWI) - SRME/EPSI
More informationMultiple-Choice Test Chapter Golden Section Search Method Optimization COMPLETE SOLUTION SET
Mltiple-Choice Test Chapter 09.0 Golden Section Search Method Optimization COMPLETE SOLUTION SET. Which o the ollowing statements is incorrect regarding the Eqal Interval Search and Golden Section Search
More informationIncoherent noise suppression with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC
Incoherent noise suppression with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC SUMMARY The separation of signal and noise is a key issue in seismic data processing. By
More informationSparsity-promoting migration accelerated by message passing Xiang LI and Felix J.Herrmann EOS-UBC
Sparsity-promoting migration accelerated by message passing Xiang LI and Feli J.Herrmann EOS-UBC SUMMARY Seismic imaging via linearized inversion requires multiple iterations to the least-squares misfit
More informationSZ-1.4: Significantly Improving Lossy Compression for Scientific Data Sets Based on Multidimensional Prediction and Error- Controlled Quantization
SZ-1.4: Significantly Improving Lossy Compression for Scientific Data Sets Based on Mltidimensional Prediction and Error- Controlled Qantization Dingwen Tao (University of California, Riverside) Sheng
More informationOn the Computational Complexity and Effectiveness of N-hub Shortest-Path Routing
1 On the Comptational Complexity and Effectiveness of N-hb Shortest-Path Roting Reven Cohen Gabi Nakibli Dept. of Compter Sciences Technion Israel Abstract In this paper we stdy the comptational complexity
More informationSUMMARY. Pursuit De-Noise (BPDN) problem, Chen et al., 2001; van den Berg and Friedlander, 2008):
Controlling linearization errors in l 1 regularized inversion by rerandomization Ning Tu, Xiang Li and Felix J. Herrmann, Earth and Ocean Sciences, University of British Columbia SUMMARY Linearized inversion
More informationMaster for Co-Simulation Using FMI
Master for Co-Simlation Using FMI Jens Bastian Christoph Claß Ssann Wolf Peter Schneider Franhofer Institte for Integrated Circits IIS / Design Atomation Division EAS Zenerstraße 38, 69 Dresden, Germany
More informationCOMPOSITION OF STABLE SET POLYHEDRA
COMPOSITION OF STABLE SET POLYHEDRA Benjamin McClosky and Illya V. Hicks Department of Comptational and Applied Mathematics Rice University November 30, 2007 Abstract Barahona and Mahjob fond a defining
More informationAUTOMATIC REGISTRATION FOR REPEAT-TRACK INSAR DATA PROCESSING
AUTOMATIC REGISTRATION FOR REPEAT-TRACK INSAR DATA PROCESSING Mingsheng LIAO, Li ZHANG, Zxn ZHANG, Jiangqing ZHANG Whan Technical University of Srveying and Mapping, Natinal Lab. for Information Eng. in
More informationEvaluating Influence Diagrams
Evalating Inflence Diagrams Where we ve been and where we re going Mark Crowley Department of Compter Science University of British Colmbia crowley@cs.bc.ca Agst 31, 2004 Abstract In this paper we will
More informationFull-waveform inversion with randomized L1 recovery for the model updates
Full-waveform inversion with randomized L1 recovery for the model updates Xiang Li, Aleksandr Aravkin, Tristan van Leeuwen, Felix Herrmann March 5, 11 Abstract Full-waveform inversion (FWI) is a data fitting
More informationOne-norm regularized inversion: learning from the Pareto curve
One-norm regularized inversion: learning from the Pareto curve Gilles Hennenfent and Felix J. Herrmann, Earth & Ocean Sciences Dept., the University of British Columbia ABSTRACT Geophysical inverse problems
More informationImage Enhancement in the Frequency Domain Periodicity and the need for Padding:
Prepared B: Dr. Hasan Demirel PhD Image Enhancement in the Freqenc Domain Periodicit and the need for Padding: Periodicit propert of the DFT: The discrete Forier Transform and the inverse Forier transforms
More informationCS 153 Design of Operating Systems Spring 18
CS 53 Design of Operating Systems Spring 8 Lectre 2: Virtal Memory Instrctor: Chengy Song Slide contribtions from Nael Ab-Ghazaleh, Harsha Madhyvasta and Zhiyn Qian Recap: cache Well-written programs exhibit
More informationAn SVD-free Pareto curve approach to rank minimization
An SVD-free Pareto curve approach to rank minimization Aleksandr Y. Aravkin IBM T.J. Watson Research Center Yorktown Heights, NY 598, USA Rajiv Mittal Department of Earth and Ocean Sciences, University
More informationProf. Kozyrakis. 1. (10 points) Consider the following fragment of Java code:
EE8 Winter 25 Homework #2 Soltions De Thrsday, Feb 2, 5 P. ( points) Consider the following fragment of Java code: for (i=; i
More informationSub-Nyquist sampling and sparsity: getting more information from fewer samples
Sub-Nyquist sampling and sparsity: getting more information from fewer samples Felix J. Herrmann SLIM Seismic Laboratory for Imaging and Modeling the University of British Columbia SLIM Drivers Our incessant
More informationCS 4204 Computer Graphics
CS 424 Compter Graphics Crves and Srfaces Yong Cao Virginia Tech Reference: Ed Angle, Interactive Compter Graphics, University of New Mexico, class notes Crve and Srface Modeling Objectives Introdce types
More informationSeismic data reconstruction with Generative Adversarial Networks
Seismic data reconstruction with Generative Adversarial Networks Ali Siahkoohi 1, Rajiv Kumar 1,2 and Felix J. Herrmann 2 1 Seismic Laboratory for Imaging and Modeling (SLIM), The University of British
More informationImage Denoising Algorithms
Image Denoising Algorithms Xiang Hao School of Compting, University of Utah, USA, hao@cs.tah.ed Abstract. This is a report of an assignment of the class Mathematics of Imaging. In this assignment, we first
More informationG009 Scale and Direction-guided Interpolation of Aliased Seismic Data in the Curvelet Domain
G009 Scale and Direction-guided Interpolation of Aliased Seismic Data in the Curvelet Domain M. Naghizadeh* (University of Alberta) & M. Sacchi (University of Alberta) SUMMARY We propose a robust interpolation
More informationA Frequency-optimized Discontinuous Formulation for Wave Propagation Problems
th Flid Dynamics Conference and Ehibit 8 Jne - Jly, Chicago, Illinois AIAA -55 A Freqency-optimized Discontinos Formlation for Wave Propagation Problems Yi Li and Z. J. Wang Department of Aerospace Engineering
More informationTopic Continuity for Web Document Categorization and Ranking
Topic Continity for Web ocment Categorization and Ranking B. L. Narayan, C. A. Mrthy and Sankar. Pal Machine Intelligence Unit, Indian Statistical Institte, 03, B. T. Road, olkata - 70008, India. E-mail:
More informationSUMMARY. These two projections are illustrated in the equation
Mitigating artifacts in Projection Onto Convex Sets interpolation Aaron Stanton*, University of Alberta, Mauricio D. Sacchi, University of Alberta, Ray Abma, BP, and Jaime A. Stein, Geotrace Technologies
More informationSeismic data interpolation using nonlinear shaping regularization a
Seismic data interpolation using nonlinear shaping regularization a a Published in Journal of Seismic Exploration, 24, no. 5, 327-342, (2015) Yangkang Chen, Lele Zhang and Le-wei Mo ABSTRACT Seismic data
More informationDigital Image Processing Chapter 5: Image Restoration
Digital Image Processing Chapter 5: Image Restoration Concept of Image Restoration Image restoration is to restore a degraded image back to the original image while image enhancement is to maniplate the
More informationFault Tolerance in Hypercubes
Falt Tolerance in Hypercbes Shobana Balakrishnan, Füsn Özgüner, and Baback A. Izadi Department of Electrical Engineering, The Ohio State University, Colmbs, OH 40, USA Abstract: This paper describes different
More informationLecture 10. Diffraction. incident
1 Introdction Lectre 1 Diffraction It is qite often the case that no line-of-sight path exists between a cell phone and a basestation. In other words there are no basestations that the cstomer can see
More informationEECS 487: Interactive Computer Graphics f
Interpolating Key Vales EECS 487: Interactive Compter Graphics f Keys Lectre 33: Keyframe interpolation and splines Cbic splines The key vales of each variable may occr at different frames The interpolation
More informationPARALLEL MRI RECONSTRUCTION USING SVD-AND- LAPLACIAN TRANSFORM BASED SPARSITY REGULARIZATION
PARALLEL MRI RECONSTRUCTION USING SVD-AND- LAPLACIAN TRANSFORM BASED SPARSITY REGULARIZATION JING BIAN, YAMING WANG, HUAXIONG ZHANG, LITAO ZHU, MINGFENG JIANG School of Information Science and Technology,
More informationNormalized averaging using adaptive applicability functions with application in image reconstruction from sparsely and randomly sampled data
Normalized averaging sing adaptive applicability fnctions with application in image reconstrction from sparsely and randomly sampled data Tan Q. Pham, Lcas J. van Vliet Pattern Recognition Grop, Faclty
More informationPOWER-OF-2 BOUNDARIES
Warren.3.fm Page 5 Monday, Jne 17, 5:6 PM CHAPTER 3 POWER-OF- BOUNDARIES 3 1 Ronding Up/Down to a Mltiple of a Known Power of Ronding an nsigned integer down to, for eample, the net smaller mltiple of
More informationABSOLUTE DEFORMATION PROFILE MEASUREMENT IN TUNNELS USING RELATIVE CONVERGENCE MEASUREMENTS
Proceedings th FIG Symposim on Deformation Measrements Santorini Greece 00. ABSOUTE DEFORMATION PROFIE MEASUREMENT IN TUNNES USING REATIVE CONVERGENCE MEASUREMENTS Mahdi Moosai and Saeid Khazaei Mining
More informationPicking and Curves Week 6
CS 48/68 INTERACTIVE COMPUTER GRAPHICS Picking and Crves Week 6 David Breen Department of Compter Science Drexel University Based on material from Ed Angel, University of New Mexico Objectives Picking
More informationXiang Li. Temple University. September 6, 2016
Introdction to Matlab Xiang Li Temple University September 6, 2016 Otline Matlab Introdction Get yor own copy of Matlab Data Types Matrices Operators MAT-LAB Introdction Matlab is a programming langage
More informationMostafa Naghizadeh and Mauricio D. Sacchi
Ground-roll elimination by scale and direction guided curvelet transform Mostafa Naghizadeh and Mauricio D. Sacchi Summary We propose a curvelet domain strategy to subtract ground-roll from seismic records.
More informationSimultaneous source separation via robust time variant Radon operators
Simultaneous source separation via robust time variant Radon operators Summary Amr Ibrahim, University of Alberta, Edmonton, Canada amr.ibrahim@ualberta.ca and Mauricio D. Sacchi, University of Alberta,
More informationStereo Matching and 3D Visualization for Gamma-Ray Cargo Inspection
Stereo Matching and 3D Visalization for Gamma-Ray Cargo Inspection Zhigang Zh *ab, Y-Chi H bc a Department of Compter Science, The City College of New York, New York, NY 3 b Department of Compter Science,
More informationOCIS codes: ( ) Microscopy; ( ) Partial coherence in imaging; ( ) Image reconstruction techniques
Sparsity based sb-wavelength imaging with partially incoherent light via qadratic compressed sensing Yoav Shechtman 1, Yonina C. Eldar, Alexander Szameit 1 and Mordechai Segev 1 1 Physics Department and
More informationConstructing and Comparing User Mobility Profiles for Location-based Services
Constrcting and Comparing User Mobility Profiles for Location-based Services Xihi Chen Interdisciplinary Centre for Secrity, Reliability and Trst, University of Lxemborg Jn Pang Compter Science and Commnications,
More informationPlenoPatch: Patch-based Plenoptic Image Manipulation
1 PlenoPatch: Patch-based Plenoptic Image Maniplation Fang-Le Zhang, Member, IEEE, Je Wang, Senior Member, IEEE, Eli Shechtman, Member, IEEE, Zi-Ye Zho, Jia-Xin Shi, and Shi-Min H, Member, IEEE Abstract
More informationThe LS-STAG Method : A new Immersed Boundary / Level-Set Method for the Computation of Incompressible Viscous Flows in Complex Geometries
The LS-STAG Method : A new Immersed Bondary / Level-Set Method for the Comptation of Incompressible Viscos Flows in Complex Geometries Yoann Cheny & Olivier Botella Nancy Universités LEMTA - UMR 7563 (CNRS-INPL-UHP)
More informationHardware-Accelerated Free-Form Deformation
Hardware-Accelerated Free-Form Deformation Clint Cha and Ulrich Nemann Compter Science Department Integrated Media Systems Center University of Sothern California Abstract Hardware-acceleration for geometric
More informationA RECOGNITION METHOD FOR AIRPLANE TARGETS USING 3D POINT CLOUD DATA
A RECOGNITION METHOD FOR AIRPLANE TARGETS USING 3D POINT CLOUD DATA Mei Zho*, Ling-li Tang, Chan-rong Li, Zhi Peng, Jing-mei Li Academy of Opto-Electronics, Chinese Academy of Sciences, No.9, Dengzhang
More informationAlliances and Bisection Width for Planar Graphs
Alliances and Bisection Width for Planar Graphs Martin Olsen 1 and Morten Revsbæk 1 AU Herning Aarhs University, Denmark. martino@hih.a.dk MADAGO, Department of Compter Science Aarhs University, Denmark.
More informationPlenoPatch: Patch-based Plenoptic Image Manipulation
1 PlenoPatch: Patch-based Plenoptic Image Maniplation Fang-Le Zhang, Member, IEEE, Je Wang, Senior Member, IEEE, Eli Shechtman, Member, IEEE, Zi-Ye Zho, Jia-Xin Shi, and Shi-Min H, Member, IEEE Abstract
More informationCurves and Surfaces. CS 537 Interactive Computer Graphics Prof. David E. Breen Department of Computer Science
Crves and Srfaces CS 57 Interactive Compter Graphics Prof. David E. Breen Department of Compter Science E. Angel and D. Shreiner: Interactive Compter Graphics 6E Addison-Wesley 22 Objectives Introdce types
More informationAn Introduction to GPU Computing. Aaron Coutino MFCF
An Introdction to GPU Compting Aaron Cotino acotino@waterloo.ca MFCF What is a GPU? A GPU (Graphical Processing Unit) is a special type of processor that was designed to render and maniplate textres. They
More informationP257 Transform-domain Sparsity Regularization in Reconstruction of Channelized Facies
P257 Transform-domain Sparsity Regularization in Reconstruction of Channelized Facies. azemi* (University of Alberta) & H.R. Siahkoohi (University of Tehran) SUMMARY Petrophysical reservoir properties,
More informationFast Obstacle Detection using Flow/Depth Constraint
Fast Obstacle etection sing Flow/epth Constraint S. Heinrich aimlerchrylser AG P.O.Box 2360, -89013 Ulm, Germany Stefan.Heinrich@aimlerChrysler.com Abstract The early recognition of potentially harmfl
More informationVideo Denoising using Transform Domain Method
Video Denoising sing Transform Domain Method Lina Desale 1, Prof. S. B. Borse 1 E&TC Dept, LGNSCOE, Nashik E&TC Dept, LGNSCOE, Nashik Abstract Utmost prevailing practical digital video denoising techniqes
More informationFINITE ELEMENT APPROXIMATION OF CONVECTION DIFFUSION PROBLEMS USING GRADED MESHES
FINITE ELEMENT APPROXIMATION OF CONVECTION DIFFUSION PROBLEMS USING GRADED MESHES RICARDO G. DURÁN AND ARIEL L. LOMBARDI Abstract. We consider the nmerical approximation of a model convection-diffsion
More informationEFFECTS OF OPTICALLY SHALLOW BOTTOMS ON WATER-LEAVING RADIANCES. Curtis D. Mobley and Lydia Sundman Sequoia Scientific, Inc. th
EFFECTS OF OPTICALLY SHALLOW BOTTOMS ON WATER-LEAVING RADIANCES Crtis D. Mobley and Lydia Sndman Seqoia Scientific, Inc. th 15317 NE 90 Street Redmond, WA 98052 USA mobley@seqoiasci.com phone: 425-867-2464
More informationRandomized sampling and sparsity: getting more information from fewer samples
Randomized sampling and sparsity: getting more information from fewer samples Felix J. Herrmann 1 ABSTRACT Many seismic exploration techniques rely on the collection of massive data volumes that are subsequently
More informationThe single-cycle design from last time
lticycle path Last time we saw a single-cycle path and control nit for or simple IPS-based instrction set. A mlticycle processor fies some shortcomings in the single-cycle CPU. Faster instrctions are not
More informationSUMMARY. forward/inverse discrete curvelet transform matrices (defined by the fast discrete curvelet transform, FDCT, with wrapping
Seismic data processing with curvelets: a multiscale and nonlinear approach Felix J. Herrmann, EOS-UBC, Deli Wang, Jilin University and Gilles Hennenfent and Peyman Moghaddam SUMMARY In this abstract,
More informationNetworks An introduction to microcomputer networking concepts
Behavior Research Methods& Instrmentation 1978, Vol 10 (4),522-526 Networks An introdction to microcompter networking concepts RALPH WALLACE and RICHARD N. JOHNSON GA TX, Chicago, Illinois60648 and JAMES
More informationB023 Seismic Data Interpolation by Orthogonal Matching Pursuit
B023 Seismic Data Interpolation by Orthogonal Matching Pursuit Y. Hollander* (Paradigm Geophysical), D. Kosloff (Paradigm Geophysical), Z. Koren (Paradigm Geophysical) & A. Bartana (Paradigm Geophysical)
More informationPara-catadioptric camera auto-calibration from epipolar geometry
CENTER FOR MACHINE PERCEPTION CZECH TECHNICAL UNIVERSITY Para-catadioptric camera ato-calibration from epipolar geometry Branislav Mičšík and Tomáš Pajdla micsb1@cmp.felk.cvt.cz, pajdla@cmp.felk.cvt.cz
More informationReview. How to represent real numbers
PCWrite PC IorD Review ALUSrcA emread Address Write data emory emwrite em Data IRWrite [3-26] [25-2] [2-6] [5-] [5-] RegDst Read register Read register 2 Write register Write data RegWrite Read data Read
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 informationCompressed sensing based land simultaneous acquisition using encoded sweeps
Compressed sensing based land simultaneous acquisition using encoded sweeps Rajiv Kumar 1, Shashin Sharan 1, Nick Moldoveanu 2, and Felix J. Herrmann 3 1 University of British Columbia, Vancouver 2 WesternGeco,
More information3-D SURFACE ROUGHNESS PROFILE OF 316-STAINLESS STEEL USING VERTICAL SCANNING INTERFEROMETRY WITH A SUPERLUMINESCENT DIODE
IMEKO 1 TC3, TC5 and TC Conferences Metrology in Modern Context November 5, 1, Pattaya, Chonbri, Thailand 3-D SURFACE ROUGHNESS PROFILE OF 316-STAINLESS STEEL USING VERTICAL SCANNING INTERFEROMETRY WITH
More informationTAKING THE PULSE OF ICT IN HEALTHCARE
ICT TODAY THE OFFICIAL TRADE JOURNAL OF BICSI Janary/Febrary 2016 Volme 37, Nmber 1 TAKING THE PULSE OF ICT IN HEALTHCARE + PLUS + High-Power PoE + Using HDBaseT in AV Design for Schools + Focs on Wireless
More informationSUMMARY METHOD. d(t 2 = τ 2 + x2
Yujin Liu, Zhi Peng, William W. Symes and Wotao Yin, China University of Petroleum (Huadong), Rice University SUMMARY The Radon transform suffers from the typical problems of loss of resolution and aliasing
More informationVessel Tracking for Retina Images Based on Fuzzy Ant Colony Algorithm
Vessel Tracking for Retina Images Based on Fzzy Ant Colony Algorithm Sina Hooshyar a, Rasol Khayati b and Reza Rezai c a,b,c Department of Biomedical Engineering, Engineering Faclty, Shahed University,
More informationMETAMODEL FOR SOFTWARE SOLUTIONS IN COMPUTED TOMOGRAPHY
VOL. 10, NO 22, DECEBER, 2015 ISSN 1819-6608 ETAODEL FOR SOFTWARE SOLUTIONS IN COPUTED TOOGRAPHY Vitaliy ezhyev Faclty of Compter Systems and Software Engineering, Universiti alaysia Pahang, Gambang, alaysia
More informationCautionary Aspects of Cross Layer Design: Context, Architecture and Interactions
Cationary Aspects of Cross Layer Design: Context, Architectre and Interactions Vikas Kawadia and P. R. Kmar Dept. of Electrical and Compter Engineering, and Coordinated Science Lab University of Illinois,
More informationPavlin and Daniel D. Corkill. Department of Computer and Information Science University of Massachusetts Amherst, Massachusetts 01003
From: AAAI-84 Proceedings. Copyright 1984, AAAI (www.aaai.org). All rights reserved. SELECTIVE ABSTRACTION OF AI SYSTEM ACTIVITY Jasmina Pavlin and Daniel D. Corkill Department of Compter and Information
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 informationTowards Understanding Bilevel Multi-objective Optimization with Deterministic Lower Level Decisions
Towards Understanding Bilevel Mlti-objective Optimization with Deterministic Lower Level Decisions Ankr Sinha Ankr.Sinha@aalto.fi Department of Information and Service Economy, Aalto University School
More informationLayered Light Field Reconstruction for Defocus Blur
Layered Light Field Reconstrction for Defocs Blr KARTHIK VAIDYANATHAN, JACOB MUNKBERG, PETRIK CLARBERG and MARCO SALVI Intel Corporation We present a novel algorithm for reconstrcting high qality defocs
More informationSeismic wavefield inversion with curvelet-domain sparsity promotion Felix J. Herrmann, EOS-UBC and Deli Wang, Jilin University
Seismic wavefield inversion with curvelet-domain sparsity promotion Felix J. Herrmann, EOS-UBC and Deli Wang, Jilin University SUMMARY Inverting seismic wavefields lies at the heart of seismic data processing
More informationStatistical Methods in functional MRI. Standard Analysis. Data Processing Pipeline. Multiple Comparisons Problem. Multiple Comparisons Problem
Statistical Methods in fnctional MRI Lectre 7: Mltiple Comparisons 04/3/13 Martin Lindqist Department of Biostatistics Johns Hopkins University Data Processing Pipeline Standard Analysis Data Acqisition
More informationCompressive sensing in seismic exploration: an outlook on a new paradigm
Compressive sensing in seismic exploration: an outlook on a new paradigm Felix J. Herrmann 1, Haneet Wason 1, and Tim T.Y. Lin 1 ABSTRACT Many seismic exploration techniques rely on the collection of massive
More informationTemporal Light Field Reconstruction for Rendering Distribution Effects
Temporal Light Field Reconstrction for Rendering Distribtion Effects Jaakko Lehtinen NVIDIA Research PBRT, 16 spp, 403 s Timo Aila NVIDIA Research Jiawen Chen MIT CSAIL PBRT, 256 spp, 6426 s Samli Laine
More informationSeismic data interpolation beyond aliasing using regularized nonstationary autoregression a
Seismic data interpolation beyond aliasing using regularized nonstationary autoregression a a Published in Geophysics, 76, V69-V77, (2011) Yang Liu, Sergey Fomel ABSTRACT Seismic data are often inadequately
More informationThe extra single-cycle adders
lticycle Datapath As an added bons, we can eliminate some of the etra hardware from the single-cycle path. We will restrict orselves to sing each fnctional nit once per cycle, jst like before. Bt since
More informationh-vectors of PS ear-decomposable graphs
h-vectors of PS ear-decomposable graphs Nima Imani 2, Lee Johnson 1, Mckenzie Keeling-Garcia 1, Steven Klee 1 and Casey Pinckney 1 1 Seattle University Department of Mathematics, 901 12th Avene, Seattle,
More informationReview Multicycle: What is Happening. Controlling The Multicycle Design
Review lticycle: What is Happening Reslt Zero Op SrcA SrcB Registers Reg Address emory em Data Sign etend Shift left Sorce A B Ot [-6] [5-] [-6] [5-] [5-] Instrction emory IR RegDst emtoreg IorD em em
More informationDiscrete Cost Multicommodity Network Optimization Problems and Exact Solution Methods
Annals of Operations Research 106, 19 46, 2001 2002 Klwer Academic Pblishers. Manfactred in The Netherlands. Discrete Cost Mlticommodity Network Optimization Problems and Exact Soltion Methods MICHEL MINOUX
More informationDynamic Maintenance of Majority Information in Constant Time per Update? Gudmund S. Frandsen and Sven Skyum BRICS 1 Department of Computer Science, Un
Dynamic Maintenance of Majority Information in Constant Time per Update? Gdmnd S. Frandsen and Sven Skym BRICS 1 Department of Compter Science, University of arhs, Ny Mnkegade, DK-8000 arhs C, Denmark
More informationThe Impact of Avatar Mobility on Distributed Server Assignment for Delivering Mobile Immersive Communication Environment
This fll text paper was peer reviewed at the direction of IEEE Commnications Society sbject matter experts for pblication in the ICC 27 proceedings. The Impact of Avatar Mobility on Distribted Server Assignment
More informationA personalized search using a semantic distance measure in a graph-based ranking model
Article A personalized search sing a semantic distance measre in a graph-based ranking model Jornal of Information Science XX (X) pp. 1-23 The Athor(s) 2011 Reprints and Permissions: sagepb.co.k/jornalspermissions.nav
More informationMulti-lingual Multi-media Information Retrieval System
Mlti-lingal Mlti-media Information Retrieval System Shoji Mizobchi, Sankon Lee, Fmihiko Kawano, Tsyoshi Kobayashi, Takahiro Komats Gradate School of Engineering, University of Tokshima 2-1 Minamijosanjima,
More informationSketch-Based Aesthetic Product Form Exploration from Existing Images Using Piecewise Clothoid Curves
Sketch-Based Aesthetic Prodct Form Exploration from Existing Images Using Piecewise Clothoid Crves Günay Orbay, Mehmet Ersin Yümer, Levent Brak Kara* Mechanical Engineering Department Carnegie Mellon University
More information