Seismic compression. François G. Meyer Department of Electrical Engineering University of Colorado at Boulder
|
|
- Kathleen Manning
- 5 years ago
- Views:
Transcription
1 Seismic compression François G. Meyer Department of Electrical Engineering University of Colorado at Boulder fmeyer IPAM, MGA 2004 Seismic Compression 1
2 Collaborators Amir Averbuch Raphy Coifman Jan-Olov Strömberg Anthony Vassiliou IPAM, MGA 2004 Seismic Compression 2
3 Introduction seismic data: HUGE ( 100 Tbyte) datasets: large dynamic range, geometric features: important seismic data: highly oscillatory Wavelet compression of seismic data [Ergas et al., 1995, Reiter and Heller, 1994, Vassiliou and Wickerhauser, 1997] But: rapid variations of intensity many fine scale coefficients very diffuse representations in a standard wavelet basis IPAM, MGA 2004 Seismic Compression 3
4 Marine shot gathers original, range=[ , ], size = Compression by 20, SNR=35. Error, range =[ 10 6, 10 6 ]. IPAM, MGA 2004 Seismic Compression 4
5 Common depth gathers original, range=[ 4000, 4000], size= Compression by 20, SNR=46. Error, range =[ 10, 10]. IPAM, MGA 2004 Seismic Compression 5
6 Performance comparison compression wavelet LCT ratio (SNR) (SNR) Common depth gathers compression wavelet LCT ratio (SNR) (SNR) Marine shot gathers SNR = l2 ( f ) l 2 ( f f ) (1) Local cosine outperform wavelets for seismic datasets! IPAM, MGA 2004 Seismic Compression 6
7 Adaptive smooth local cosine transforms b n biorthogonal bases a α a a a a n - n n n + β n a - + β n+1 α n+1 n+1 n+1 n+1 R = n=+ n= [a n, a n+1 [ neighborhood around each point a n : [a n α n, a n + β n ] a n + β n < a n+1 α n+1. (2) b n : bell function lives over the interval [a n α n, a n+1 + β n+1 ]. IPAM, MGA 2004 Seismic Compression 7
8 Dual bells b n (x) = θ n (x)b n 1 (2a n x) if x [a n α n, a n + β n ] 1 b n (x) if x [a n + β n, a n+1 α n+1 ] θ n+1 (x)b n+1 (2a n+1 x) if x [a n+1 α n+1, a n+1 + β n+1 ] 0 otherwise (3) θ n (x) = 1 b n (x) b n 1 (2a n x) + b n (2a n x) b n 1 (x) (4) IPAM, MGA 2004 Seismic Compression 8
9 Biorthogonal local cosine bases DCT-IV: c n,k (x) = [ ] 2 x a n cos (k + 1/2) a n+1 a n a n+1 a n (5) We define the family and the dual family: w n,k = b n (x) c n,k (x) (6) w n,k = b n (x) c n,k (x) (7) IPAM, MGA 2004 Seismic Compression 9
10 Biorthogonal local cosine bases Lemma 1 w n,j and w n,j are Riesz biorthogonal bases: w n,j (x) w k,m (x) dx = δ j,k δ n,m (8) x L 2 (R), x(x) = n,j x(x) = n,j x n,j w n,j (x) with x n,j = x n,j w n,j (x) with x n,j = x(x) w n,j (x)dx (9) x(x)w n,j (x)dx (10) IPAM, MGA 2004 Seismic Compression 10
11 Choice of the bell function biorthogonal bases: more freedom to select the bells b n b n can be optimized seismic data : oscillatory [Matviyenko, 1996] other choices: reproduce constants [Jawerth and Sweldens, 1995], or polynomials [Bittner, 1999] IPAM, MGA 2004 Seismic Compression 11
12 Bells of Matviyenko class of signals: x = cos(ωk + ϕ) minimum number of coefficients x n,j to reconstruct x up to an error ε 0 b(x) 1 0 b(x) ( (11) (2) + 1)/2 IPAM, MGA 2004 Seismic Compression 12
13 Bells of Matviyenko b N (x) = ( ) N g n sin(n + 1/2)πx n=0 ) ( N n=0 ( 1) n g n cos(n + 1/2)πx coefficients g n are calculated numerically if x [ 1/2, 1/2] if x [1/2 : 3/2] 0 otherwise (12) IPAM, MGA 2004 Seismic Compression 13
14 N=1 N=2 N=3 N= Optimized bells. N = 1, 2, 3, 4 IPAM, MGA 2004 Seismic Compression 14
15 N=1 N=2 N=3 N= Optimized dual bells. N = 1, 2, 3, 4 IPAM, MGA 2004 Seismic Compression 15
16 Adaptive segmentation quadtree segmentation preserve the original aspect ratio of the data anisotropic segmentations [Bennett, 2000] depth first approach: minimize memory extension on the borders IPAM, MGA 2004 Seismic Compression 16
17 Choice of a cost function entropy [Coifman and Wickerhauser, 1992] h(x) = k x k 2 x 2 log x k 2 x 2 (13) Rate distortion [Ramchandran and Vetterli, 1993] { { max λ min T T node T min q Q { D node(q) + λr node (q)} λ R Q : set of all quantizers T : set of all bases T with quadtree structure very high computational complexity! first order entropy estimates only }} (14) IPAM, MGA 2004 Seismic Compression 17
18 Choice of a cost function cost function: estimate of the actual rate achieved by each node mimics the actual scalar quantization, and entropy coding much faster to compute composed of two complementary terms: c 1 (x): cost of coding the sign and the magnitude of the non zero output levels of the scalar quantizer, c 2 (x): cost of coding the locations of the non zero output levels (significance map), IPAM, MGA 2004 Seismic Compression 18
19 Choice of a cost function c 1 (x) = k/q(xk ) =0 max (log 2 Q(x k ), 0) fast implementation: representation of floating numbers c 2 (x) = N ( p log 2 (p) + (1 p) log 2 (1 p) ) first order entropy of a Bernoulli process: each coefficient x k is significant with a probability p IPAM, MGA 2004 Seismic Compression 19
20 Fast DCT-IV DCT-IV : FFT of half length DCT-IV coefficients, ˆx(j), j = 0,..., N 1 of the sequence x(n), n = 0,..., N 1 ˆx(j) =Re e ijπ 2N N/2 1 y(n)e 2iπjn N ˆx(N j 1) = Im n=0 e ijπ 2N N/2 1 n=0 y(n)e 2iπjn N (15) with + 1/4)π i(n y(n) = (x(2n) + i x(n 2n 1)) e N (16) IPAM, MGA 2004 Seismic Compression 20
21 Scanning the coefficients scanning coefficients by increasing frequency each LCT block divided into a fixed number of frequency subsets gather from all the LCT blocks all the coefficients that are in the same subset IPAM, MGA 2004 Seismic Compression 21
22 Scanning the coefficients IPAM, MGA 2004 Seismic Compression 22
23 IPAM, MGA 2004 Seismic Compression 23
24 IPAM, MGA 2004 Seismic Compression 24
25 Laplacian based scalar quantization distribution of the cosine coefficients: Laplacian [Birney and Fischer, 1995] near optimal scalar quantizer [Sullivan, 1996] [ + δ, δ], the symmetric dead-zone,, the quantizer step size, δ, the reconstruction offset IPAM, MGA 2004 Seismic Compression 25
26 Laplacian based scalar quantization 2 δ δ -3 +δ - 2 +δ - + δ - δ - δ 2 - δ 3 - δ δ - 2 Scalar quantizer, with a dead zone. IPAM, MGA 2004 Seismic Compression 26
27 Entropy coding significance map: n C order arithmetic coder signs of the output levels: packed magnitude of the output levels: variable length encoded best basis geometry: adaptive arithmetic coder. IPAM, MGA 2004 Seismic Compression 27
28 Experiments Fast Local Cosine Transform (FLCT) coder and decoder actual bit stream: size equal to the targeted budget IPAM, MGA 2004 Seismic Compression 28
29 Comparison of the bells Standard bell: [ π ] b N (x) = sin 4 (1 + x N) x j = sin( π 2 x j 1) x 0 = x (17) Optimized bell: Matviyenko Gaussian b(x) = e α(x a n) 2 IPAM, MGA 2004 Seismic Compression 29
30 1 0.9 alpha=1 alpha= alpha=1 alpha= Gaussian bells. α = 1, 2 IPAM, MGA 2004 Seismic Compression 30
31 N=1 N=2 N=3 N= Standard bells. N = 1, 2, 3, 4 IPAM, MGA 2004 Seismic Compression 31
32 Test images Synthetic image: warped cosines, Seismic data: 2D slice uniform grid: blocks of size IPAM, MGA 2004 Seismic Compression 32
33 IPAM, MGA 2004 Seismic Compression 33
34 IPAM, MGA 2004 Seismic Compression 34
35 Synthetic data : optimal N for optimized bell N=1 N=2 N=3 N=4 N=5 44 PSNR (db) Compression ratio IPAM, MGA 2004 Seismic Compression 35
36 Synthetic data : optimal overlap for optimized bell PSNR (db) Compression ratio alpha=0 alpha=2 alpha=4 alpha=6 alpha=8 alpha=10 alpha=12 alpha=12 alpha=16 IPAM, MGA 2004 Seismic Compression 36
37 Synthetic data : comparison of the bells Matviyenko (alpha=16, N =2) Standard (alpha =16) Gauss (alpha =16) No bell Wavelet packets PSNR (db) Compression ratio IPAM, MGA 2004 Seismic Compression 37
38 Seismic data : comparison of the bells Standard (alpha=16) Matviyenko (alpha=16, N =3) Wavelet packets No bell Gauss (alpha =16) PSNR (db) Compression ratio IPAM, MGA 2004 Seismic Compression 38
39 Compression 50, Matviyenko (N = 3), PSNR =20.56 db IPAM, MGA 2004 Seismic Compression 39
40 Compression 50, no bell, PSNR = db IPAM, MGA 2004 Seismic Compression 40
41 Conclusion compression of highly oscillatory signals: local cosine transforms optimal spatial tiling design of the window full 3-D compression: fast 3-D FFT The future: curvelets? IPAM, MGA 2004 Seismic Compression 41
42 References [Bennett, 2000] Bennett, N. (2000). Fast algorithm for best anisotropic Walsh bases and relatives. Applied and Computational Harmonic Analysis, 8(1): [Birney and Fischer, 1995] Birney, K. and Fischer, T. (1995). On the modeling of DCT and subbdand image data for compression. IEEE Trans. on Image Process., 4(2): [Bittner, 1999] Bittner, K. (1999). Error estimates and reproduction of polynomials for biorthogonal local trigonometric bases. Applied and Computational Harmonic Analysis, 6: [Coifman and Wickerhauser, 1992] Coifman, R. and Wickerhauser, M. (1992). Entropy-based algorithms for best basis selection. IEEE Trans. Information Theory, 38(2): [Ergas et al., 1995] Ergas, R., Donoho, P., and Villassenor, J. (1995). High-performance seismic trace compression. In Soc. Expl. Geophys., 65th Intern. Conv. Soc. Expl. Geophys. [Jawerth and Sweldens, 1995] Jawerth, B. and Sweldens, W. (1995). Biorthogonal local trigonometric bases. J. Fourier Anal. Appl., 2(2): [Matviyenko, 1996] Matviyenko, G. (1996). Optimized local trigonometric bases. Applied and Computational Harmonic Analysis, 3: [Ramchandran and Vetterli, 1993] Ramchandran, K. and Vetterli, M. (1993). Best wavelet packet bases in a ratedistortion sense. IEEE Trans. Image Processing, pages [Reiter and Heller, 1994] Reiter, E. and Heller, P. (1994). Wavelet transformation-based compression of nmocorrected cdp gathers. In Soc. Expl. Geophys., 64th Intern. Conv. Soc. Expl. Geophys. [Sullivan, 1996] Sullivan, G. (1996). Efficient scalar quantization of exponential and Laplacian random variables. IEEE Trans. Information Theory, 42(5): IPAM, MGA 2004 Seismic Compression 42
43 [Vassiliou and Wickerhauser, 1997] Vassiliou, A. and Wickerhauser, M. (1997). Comparison of wavelet image coding schemes for seismic data compression. In Soc. Expl. Geophys., 67th Intern. Conv. Soc. Expl. Geophys. IPAM, MGA 2004 Seismic Compression 43
Wavelet transforms and compression of seismic data
Mathematical Geophysics Summer School, Stanford, August 1998 Wavelet transforms and compression of seismic data G. Beylkin 1 and A. Vassiliou 2 I Introduction Seismic data compression (as it exists today)
More informationA NEW BASIS SELECTION PARADIGM FOR WAVELET PACKET IMAGE CODING
A NEW BASIS SELECTION PARADIGM FOR WAVELET PACKET IMAGE CODING Nasir M. Rajpoot, Roland G. Wilson, François G. Meyer, Ronald R. Coifman Corresponding Author: nasir@dcs.warwick.ac.uk ABSTRACT In this paper,
More informationSCALABLE AUDIO CODING USING WATERMARKING
SCALABLE AUDIO CODING USING WATERMARKING Mahmood Movassagh Peter Kabal Department of Electrical and Computer Engineering McGill University, Montreal, Canada Email: {mahmood.movassagh@mail.mcgill.ca, peter.kabal@mcgill.ca}
More informationImage Analysis and Approximation via Generalized Polyharmonic Local Trigonometric Transforms
Image Analysis and Approximation via Generalized Polyharmonic Local Trigonometric Transforms Naoki Saito with Jean-François Remy, Noel Smith, Ernest Woei, Katsu Yamatani, & Jim Zhao saito@math.ucdavis.edu
More informationDesign of Optimal Quantizers for Distributed Source Coding
Design of Optimal Quantizers for Distributed Source Coding David Rebollo-Monedero, Rui Zhang and Bernd Girod Information Systems Laboratory, Electrical Eng. Dept. Stanford University, Stanford, CA 94305
More informationConverting DCT Coefficients to H.264/AVC
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Converting DCT Coefficients to H.264/AVC Jun Xin, Anthony Vetro, Huifang Sun TR2004-058 June 2004 Abstract Many video coding schemes, including
More informationDiffusion Wavelets and Applications
Diffusion Wavelets and Applications J.C. Bremer, R.R. Coifman, P.W. Jones, S. Lafon, M. Mohlenkamp, MM, R. Schul, A.D. Szlam Demos, web pages and preprints available at: S.Lafon: www.math.yale.edu/~sl349
More informationat Some sort of quantization is necessary to represent continuous signals in digital form
Quantization at Some sort of quantization is necessary to represent continuous signals in digital form x(n 1,n ) x(t 1,tt ) D Sampler Quantizer x q (n 1,nn ) Digitizer (A/D) Quantization is also used for
More informationImage compression using an edge adapted redundant dictionary and wavelets
Image compression using an edge adapted redundant dictionary and wavelets Lorenzo Peotta, Lorenzo Granai, Pierre Vandergheynst Signal Processing Institute Swiss Federal Institute of Technology EPFL-STI-ITS-LTS2
More informationAn Investigation of 3D Dual-Tree Wavelet Transform for Video Coding
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com An Investigation of 3D Dual-Tree Wavelet Transform for Video Coding Beibei Wang, Yao Wang, Ivan Selesnick and Anthony Vetro TR2004-132 December
More informationEstimation-Theoretic Delayed Decoding of Predictively Encoded Video Sequences
Estimation-Theoretic Delayed Decoding of Predictively Encoded Video Sequences Jingning Han, Vinay Melkote, and Kenneth Rose Department of Electrical and Computer Engineering University of California, Santa
More informationDigital Image Processing Lectures 25 & 26
Lectures 25 & 26, Professor Department of Electrical and Computer Engineering Colorado State University Spring 2015 Area 4: Image Encoding and Compression Goal: To exploit the redundancies in the image
More informationStatistical Analysis and Distortion Modeling of MPEG-4 FGS
Statistical Analysis and Distortion Modeling of MPEG-4 FGS Min Dai Electrical Engineering Texas A&M University, TX 77843 Dmitri Loguinov Computer Science Texas A&M University, TX 77843 Hayder Radha Hayder
More informationBasic Principles of Video Coding
Basic Principles of Video Coding Introduction Categories of Video Coding Schemes Information Theory Overview of Video Coding Techniques Predictive coding Transform coding Quantization Entropy coding Motion
More informationEE368B Image and Video Compression
EE368B Image and Video Compression Homework Set #2 due Friday, October 20, 2000, 9 a.m. Introduction The Lloyd-Max quantizer is a scalar quantizer which can be seen as a special case of a vector quantizer
More informationDirectionlets. Anisotropic Multi-directional Representation of Images with Separable Filtering. Vladan Velisavljević Deutsche Telekom, Laboratories
Directionlets Anisotropic Multi-directional Representation of Images with Separable Filtering Vladan Velisavljević Deutsche Telekom, Laboratories Google Inc. Mountain View, CA October 2006 Collaborators
More informationECE533 Digital Image Processing. Embedded Zerotree Wavelet Image Codec
University of Wisconsin Madison Electrical Computer Engineering ECE533 Digital Image Processing Embedded Zerotree Wavelet Image Codec Team members Hongyu Sun Yi Zhang December 12, 2003 Table of Contents
More informationDesign and Implementation of Multistage Vector Quantization Algorithm of Image compression assistant by Multiwavelet Transform
Design and Implementation of Multistage Vector Quantization Algorithm of Image compression assistant by Multiwavelet Transform Assist Instructor/ BASHAR TALIB HAMEED DIYALA UNIVERSITY, COLLEGE OF SCIENCE
More informationThe Polyharmonic Local Sine Transform: A New Tool for Local Image Analysis and Synthesis without Edge Effect
The Polyharmonic Local Sine Transform: A New Tool for Local Image Analysis and Synthesis without Edge Effect Naoki Saito Jean-François Remy Department of Mathematics University of California Davis, CA
More informationA WAVELET BASED CODING SCHEME VIA ATOMIC APPROXIMATION AND ADAPTIVE SAMPLING OF THE LOWEST FREQUENCY BAND
A WAVELET BASED CODING SCHEME VIA ATOMIC APPROXIMATION AND ADAPTIVE SAMPLING OF THE LOWEST FREQUENCY BAND V. Bruni, D. Vitulano Istituto per le Applicazioni del Calcolo M. Picone, C. N. R. Viale del Policlinico
More informationFast Progressive Wavelet Coding
PRESENTED AT THE IEEE DCC 99 CONFERENCE SNOWBIRD, UTAH, MARCH/APRIL 1999 Fast Progressive Wavelet Coding Henrique S. Malvar Microsoft Research One Microsoft Way, Redmond, WA 98052 E-mail: malvar@microsoft.com
More informationECG782: Multidimensional Digital Signal Processing
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Spring 2014 TTh 14:30-15:45 CBC C313 Lecture 05 Image Processing Basics 13/02/04 http://www.ee.unlv.edu/~b1morris/ecg782/
More informationNear Optimal Signal Recovery from Random Projections
1 Near Optimal Signal Recovery from Random Projections Emmanuel Candès, California Institute of Technology Multiscale Geometric Analysis in High Dimensions: Workshop # 2 IPAM, UCLA, October 2004 Collaborators:
More informationThe information loss in quantization
The information loss in quantization The rough meaning of quantization in the frame of coding is representing numerical quantities with a finite set of symbols. The mapping between numbers, which are normally
More informationProgressive Wavelet Coding of Images
Progressive Wavelet Coding of Images Henrique Malvar May 1999 Technical Report MSR-TR-99-26 Microsoft Research Microsoft Corporation One Microsoft Way Redmond, WA 98052 1999 IEEE. Published in the IEEE
More informationWedgelets and Image Compression
Wedgelets and Image Compression Laurent Demaret, Mattia Fedrigo, Hartmut Führ Summer school: New Trends and Directions in Harmonic Analysis, Approximation Theory, and Image Analysis, Inzell, Germany, 20
More informationWavelet Scalable Video Codec Part 1: image compression by JPEG2000
1 Wavelet Scalable Video Codec Part 1: image compression by JPEG2000 Aline Roumy aline.roumy@inria.fr May 2011 2 Motivation for Video Compression Digital video studio standard ITU-R Rec. 601 Y luminance
More informationCompression methods: the 1 st generation
Compression methods: the 1 st generation 1998-2017 Josef Pelikán CGG MFF UK Praha pepca@cgg.mff.cuni.cz http://cgg.mff.cuni.cz/~pepca/ Still1g 2017 Josef Pelikán, http://cgg.mff.cuni.cz/~pepca 1 / 32 Basic
More informationEFFICIENT IMAGE COMPRESSION ALGORITHMS USING EVOLVED WAVELETS
International Journal of ISSN 0974-2107 Systems and Technologies IJST Vol.4, No.2, pp 127-146 KLEF 2011 EFFICIENT IMAGE COMPRESSION ALGORITHMS USING EVOLVED WAVELETS 1.G.CHENCHU KRISHNAIAH* 2.T.JAYACHANDRAPRASAD
More informationLecture 9 Video Coding Transforms 2
Lecture 9 Video Coding Transforms 2 Integer Transform of H.264/AVC In previous standards, the DCT was defined as the ideal transform, with unlimited accuracy. This has the problem, that we have encoders
More informationCHAPTER 3. Transformed Vector Quantization with Orthogonal Polynomials Introduction Vector quantization
3.1. Introduction CHAPTER 3 Transformed Vector Quantization with Orthogonal Polynomials In the previous chapter, a new integer image coding technique based on orthogonal polynomials for monochrome images
More informationWavelets and Image Compression Augusta State University April, 27, Joe Lakey. Department of Mathematical Sciences. New Mexico State University
Wavelets and Image Compression Augusta State University April, 27, 6 Joe Lakey Department of Mathematical Sciences New Mexico State University 1 Signals and Images Goal Reduce image complexity with little
More informationApplications of Polyspline Wavelets to Astronomical Image Analysis
VIRTUAL OBSERVATORY: Plate Content Digitization, Archive Mining & Image Sequence Processing edited by M. Tsvetkov, V. Golev, F. Murtagh, and R. Molina, Heron Press, Sofia, 25 Applications of Polyspline
More information- An Image Coding Algorithm
- An Image Coding Algorithm Shufang Wu http://www.sfu.ca/~vswu vswu@cs.sfu.ca Friday, June 14, 2002 22-1 Agenda Overview Discrete Wavelet Transform Zerotree Coding of Wavelet Coefficients Successive-Approximation
More informationWavelet Footprints: Theory, Algorithms, and Applications
1306 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 51, NO. 5, MAY 2003 Wavelet Footprints: Theory, Algorithms, and Applications Pier Luigi Dragotti, Member, IEEE, and Martin Vetterli, Fellow, IEEE Abstract
More informationMULTIPLEXING AND DEMULTIPLEXING FRAME PAIRS
MULTIPLEXING AND DEMULTIPLEXING FRAME PAIRS AZITA MAYELI AND MOHAMMAD RAZANI Abstract. Based on multiplexing and demultiplexing techniques in telecommunication, we study the cases when a sequence of several
More informationOverview. Analog capturing device (camera, microphone) PCM encoded or raw signal ( wav, bmp, ) A/D CONVERTER. Compressed bit stream (mp3, jpg, )
Overview Analog capturing device (camera, microphone) Sampling Fine Quantization A/D CONVERTER PCM encoded or raw signal ( wav, bmp, ) Transform Quantizer VLC encoding Compressed bit stream (mp3, jpg,
More informationScalar and Vector Quantization. National Chiao Tung University Chun-Jen Tsai 11/06/2014
Scalar and Vector Quantization National Chiao Tung University Chun-Jen Tsai 11/06/014 Basic Concept of Quantization Quantization is the process of representing a large, possibly infinite, set of values
More informationMGA Tutorial, September 08, 2004 Construction of Wavelets. Jan-Olov Strömberg
MGA Tutorial, September 08, 2004 Construction of Wavelets Jan-Olov Strömberg Department of Mathematics Royal Institute of Technology (KTH) Stockholm, Sweden Department of Numerical Analysis and Computer
More informationFast Algorithm for Best Anisotropic Walsh Bases and Relatives
Applied and Computational Harmonic Analysis 8, 86 103 (000) doi:10.1006/acha.1998.076, available online at http://www.idealibrary.com on Fast Algorithm for Best Anisotropic Walsh Bases and Relatives icholas.
More informationEE67I Multimedia Communication Systems
EE67I Multimedia Communication Systems Lecture 5: LOSSY COMPRESSION In these schemes, we tradeoff error for bitrate leading to distortion. Lossy compression represents a close approximation of an original
More informationBASICS OF COMPRESSION THEORY
BASICS OF COMPRESSION THEORY Why Compression? Task: storage and transport of multimedia information. E.g.: non-interlaced HDTV: 0x0x0x = Mb/s!! Solutions: Develop technologies for higher bandwidth Find
More informationSIGNAL COMPRESSION. 8. Lossy image compression: Principle of embedding
SIGNAL COMPRESSION 8. Lossy image compression: Principle of embedding 8.1 Lossy compression 8.2 Embedded Zerotree Coder 161 8.1 Lossy compression - many degrees of freedom and many viewpoints The fundamental
More informationPredictive Coding. Prediction
Predictive Coding Prediction Prediction in Images Principle of Differential Pulse Code Modulation (DPCM) DPCM and entropy-constrained scalar quantization DPCM and transmission errors Adaptive intra-interframe
More informationPredictive Coding. Prediction Prediction in Images
Prediction Prediction in Images Predictive Coding Principle of Differential Pulse Code Modulation (DPCM) DPCM and entropy-constrained scalar quantization DPCM and transmission errors Adaptive intra-interframe
More informationSIMULTANEOUS SEGMENTATION, COMPRESSION, AND DENOISING OF SIGNALS USING POLYHARMONIC LOCAL SINE TRANSFORM AND MINIMUM DESCRIPTION LENGTH CRITERION
SIMULTANEOUS SEGMENTATION, COMPRESSION, AND DENOISING OF SIGNALS USING POLYHARMONIC LOCAL SINE TRANSFORM AND MINIMUM DESCRIPTION LENGTH CRITERION Naoki Saito and Ernest Woei Department of Mathematics University
More informationFuzzy quantization of Bandlet coefficients for image compression
Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research, 2013, 4(2):140-146 Fuzzy quantization of Bandlet coefficients for image compression R. Rajeswari and R. Rajesh ISSN:
More informationMultiscale Geometric Analysis: Thoughts and Applications (a summary)
Multiscale Geometric Analysis: Thoughts and Applications (a summary) Anestis Antoniadis, University Joseph Fourier Assimage 2005,Chamrousse, February 2005 Classical Multiscale Analysis Wavelets: Enormous
More informationCompressing a 1D Discrete Signal
Compressing a D Discrete Signal Divide the signal into 8blocks. Subtract the sample mean from each value. Compute the 8 8covariancematrixforthe blocks. Compute the eigenvectors of the covariance matrix.
More informationObjectives of Image Coding
Objectives of Image Coding Representation of an image with acceptable quality, using as small a number of bits as possible Applications: Reduction of channel bandwidth for image transmission Reduction
More informationHölder regularity estimation by Hart Smith and Curvelet transforms
Hölder regularity estimation by Hart Smith and Curvelet transforms Jouni Sampo Lappeenranta University Of Technology Department of Mathematics and Physics Finland 18th September 2007 This research is done
More informationWavelets in Image Compression
Wavelets in Image Compression M. Victor WICKERHAUSER Washington University in St. Louis, Missouri victor@math.wustl.edu http://www.math.wustl.edu/~victor THEORY AND APPLICATIONS OF WAVELETS A Workshop
More informationIntroduction to Wavelets and Wavelet Transforms
Introduction to Wavelets and Wavelet Transforms A Primer C. Sidney Burrus, Ramesh A. Gopinath, and Haitao Guo with additional material and programs by Jan E. Odegard and Ivan W. Selesnick Electrical and
More informationModelling of produced bit rate through the percentage of null quantized transform coefficients ( zeros )
Rate control strategies in H264 Simone Milani (simone.milani@dei.unipd.it) with the collaboration of Università degli Studi di adova ST Microelectronics Summary General scheme of the H.264 encoder Rate
More informationVector Quantization. Institut Mines-Telecom. Marco Cagnazzo, MN910 Advanced Compression
Institut Mines-Telecom Vector Quantization Marco Cagnazzo, cagnazzo@telecom-paristech.fr MN910 Advanced Compression 2/66 19.01.18 Institut Mines-Telecom Vector Quantization Outline Gain-shape VQ 3/66 19.01.18
More informationIMAGE RECONSTRUCTION FROM UNDERSAMPLED FOURIER DATA USING THE POLYNOMIAL ANNIHILATION TRANSFORM
IMAGE RECONSTRUCTION FROM UNDERSAMPLED FOURIER DATA USING THE POLYNOMIAL ANNIHILATION TRANSFORM Anne Gelb, Rodrigo B. Platte, Rosie Renaut Development and Analysis of Non-Classical Numerical Approximation
More informationAnalysis of Rate-distortion Functions and Congestion Control in Scalable Internet Video Streaming
Analysis of Rate-distortion Functions and Congestion Control in Scalable Internet Video Streaming Min Dai Electrical Engineering, Texas A&M University Dmitri Loguinov Computer Science, Texas A&M University
More informationFast Kirchhoff migration in the wavelet domain
Exploration Geophysics (2002) 33, 23-27 Fast Kirchhoff migration in the wavelet domain Valery A. Zheludev 1 Eugene Ragoza 2 Dan D. Kosloff 3 Key Words: Kirchhoff migration, wavelet, anti-aliasing ABSTRACT
More informationGaussian source Assumptions d = (x-y) 2, given D, find lower bound of I(X;Y)
Gaussian source Assumptions d = (x-y) 2, given D, find lower bound of I(X;Y) E{(X-Y) 2 } D
More informationWAVELETS WITH COMPOSITE DILATIONS
ELECTRONIC RESEARCH ANNOUNCEMENTS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Pages 000 000 (Xxxx XX, XXXX S 1079-6762(XX0000-0 WAVELETS WITH COMPOSITE DILATIONS KANGHUI GUO, DEMETRIO LABATE, WANG-Q
More informationTwo-Dimensional Orthogonal Filter Banks with Directional Vanishing Moments
Two-imensional Orthogonal Filter Banks with irectional Vanishing Moments Jianping Zhou and Minh N. o epartment of Electrical and Computer Engineering University of Illinois at Urbana-Champaign, Urbana,
More informationEBCOT coding passes explained on a detailed example
EBCOT coding passes explained on a detailed example Xavier Delaunay d.xav@free.fr Contents Introduction Example used Coding of the first bit-plane. Cleanup pass............................. Coding of the
More informationInverse Problems in Image Processing
H D Inverse Problems in Image Processing Ramesh Neelamani (Neelsh) Committee: Profs. R. Baraniuk, R. Nowak, M. Orchard, S. Cox June 2003 Inverse Problems Data estimation from inadequate/noisy observations
More informationOn Compression Encrypted Data part 2. Prof. Ja-Ling Wu The Graduate Institute of Networking and Multimedia National Taiwan University
On Compression Encrypted Data part 2 Prof. Ja-Ling Wu The Graduate Institute of Networking and Multimedia National Taiwan University 1 Brief Summary of Information-theoretic Prescription At a functional
More informationDiffusion Geometries, Diffusion Wavelets and Harmonic Analysis of large data sets.
Diffusion Geometries, Diffusion Wavelets and Harmonic Analysis of large data sets. R.R. Coifman, S. Lafon, MM Mathematics Department Program of Applied Mathematics. Yale University Motivations The main
More informationMMSE DECODING FOR ANALOG JOINT SOURCE CHANNEL CODING USING MONTE CARLO IMPORTANCE SAMPLING
MMSE DECODING FOR ANALOG JOINT SOURCE CHANNEL CODING USING MONTE CARLO IMPORTANCE SAMPLING Yichuan Hu (), Javier Garcia-Frias () () Dept. of Elec. and Comp. Engineering University of Delaware Newark, DE
More informationDiffusion Geometries, Global and Multiscale
Diffusion Geometries, Global and Multiscale R.R. Coifman, S. Lafon, MM, J.C. Bremer Jr., A.D. Szlam, P.W. Jones, R.Schul Papers, talks, other materials available at: www.math.yale.edu/~mmm82 Data and functions
More informationon a per-coecient basis in large images is computationally expensive. Further, the algorithm in [CR95] needs to be rerun, every time a new rate of com
Extending RD-OPT with Global Thresholding for JPEG Optimization Viresh Ratnakar University of Wisconsin-Madison Computer Sciences Department Madison, WI 53706 Phone: (608) 262-6627 Email: ratnakar@cs.wisc.edu
More informationOn Scalable Coding in the Presence of Decoder Side Information
On Scalable Coding in the Presence of Decoder Side Information Emrah Akyol, Urbashi Mitra Dep. of Electrical Eng. USC, CA, US Email: {eakyol, ubli}@usc.edu Ertem Tuncel Dep. of Electrical Eng. UC Riverside,
More information2D Wavelets. Hints on advanced Concepts
2D Wavelets Hints on advanced Concepts 1 Advanced concepts Wavelet packets Laplacian pyramid Overcomplete bases Discrete wavelet frames (DWF) Algorithme à trous Discrete dyadic wavelet frames (DDWF) Overview
More informationCompression and Coding
Compression and Coding Theory and Applications Part 1: Fundamentals Gloria Menegaz 1 Transmitter (Encoder) What is the problem? Receiver (Decoder) Transformation information unit Channel Ordering (significance)
More informationCan the sample being transmitted be used to refine its own PDF estimate?
Can the sample being transmitted be used to refine its own PDF estimate? Dinei A. Florêncio and Patrice Simard Microsoft Research One Microsoft Way, Redmond, WA 98052 {dinei, patrice}@microsoft.com Abstract
More informationTransform Coding. Transform Coding Principle
Transform Coding Principle of block-wise transform coding Properties of orthonormal transforms Discrete cosine transform (DCT) Bit allocation for transform coefficients Entropy coding of transform coefficients
More informationL. Yaroslavsky. Fundamentals of Digital Image Processing. Course
L. Yaroslavsky. Fundamentals of Digital Image Processing. Course 0555.330 Lec. 6. Principles of image coding The term image coding or image compression refers to processing image digital data aimed at
More informationencoding without prediction) (Server) Quantization: Initial Data 0, 1, 2, Quantized Data 0, 1, 2, 3, 4, 8, 16, 32, 64, 128, 256
General Models for Compression / Decompression -they apply to symbols data, text, and to image but not video 1. Simplest model (Lossless ( encoding without prediction) (server) Signal Encode Transmit (client)
More informationMultiple Description Transform Coding of Images
Multiple Description Transform Coding of Images Vivek K Goyal Jelena Kovačević Ramon Arean Martin Vetterli U. of California, Berkeley Bell Laboratories École Poly. Féd. de Lausanne École Poly. Féd. de
More informationA Fast Spherical Filter with Uniform Resolution
JOURNAL OF COMPUTATIONAL PHYSICS 136, 580 584 (1997) ARTICLE NO. CP975782 A Fast Spherical Filter with Uniform Resolution Rüdiger Jakob-Chien*, and Bradley K. Alpert *Department of Computer Science & Engineering,
More informationLecture 20: Quantization and Rate-Distortion
Lecture 20: Quantization and Rate-Distortion Quantization Introduction to rate-distortion theorem Dr. Yao Xie, ECE587, Information Theory, Duke University Approimating continuous signals... Dr. Yao Xie,
More informationOn Common Information and the Encoding of Sources that are Not Successively Refinable
On Common Information and the Encoding of Sources that are Not Successively Refinable Kumar Viswanatha, Emrah Akyol, Tejaswi Nanjundaswamy and Kenneth Rose ECE Department, University of California - Santa
More informationRate-Constrained Multihypothesis Prediction for Motion-Compensated Video Compression
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS FOR VIDEO TECHNOLOGY, VOL 12, NO 11, NOVEMBER 2002 957 Rate-Constrained Multihypothesis Prediction for Motion-Compensated Video Compression Markus Flierl, Student
More informationSource Coding for Compression
Source Coding for Compression Types of data compression: 1. Lossless -. Lossy removes redundancies (reversible) removes less important information (irreversible) Lec 16b.6-1 M1 Lossless Entropy Coding,
More informationC.M. Liu Perceptual Signal Processing Lab College of Computer Science National Chiao-Tung University
Quantization C.M. Liu Perceptual Signal Processing Lab College of Computer Science National Chiao-Tung University http://www.csie.nctu.edu.tw/~cmliu/courses/compression/ Office: EC538 (03)5731877 cmliu@cs.nctu.edu.tw
More informationImage Filtering, Edges and Image Representation
Image Filtering, Edges and Image Representation Capturing what s important Req reading: Chapter 7, 9 F&P Adelson, Simoncelli and Freeman (handout online) Opt reading: Horn 7 & 8 FP 8 February 19, 8 A nice
More informationMultimedia Networking ECE 599
Multimedia Networking ECE 599 Prof. Thinh Nguyen School of Electrical Engineering and Computer Science Based on lectures from B. Lee, B. Girod, and A. Mukherjee 1 Outline Digital Signal Representation
More informationVector Quantization Encoder Decoder Original Form image Minimize distortion Table Channel Image Vectors Look-up (X, X i ) X may be a block of l
Vector Quantization Encoder Decoder Original Image Form image Vectors X Minimize distortion k k Table X^ k Channel d(x, X^ Look-up i ) X may be a block of l m image or X=( r, g, b ), or a block of DCT
More informationThe Generalized Haar-Walsh Transform (GHWT) for Data Analysis on Graphs and Networks
The Generalized Haar-Walsh Transform (GHWT) for Data Analysis on Graphs and Networks Jeff Irion & Naoki Saito Department of Mathematics University of California, Davis SIAM Annual Meeting 2014 Chicago,
More informationImage Data Compression
Image Data Compression Image data compression is important for - image archiving e.g. satellite data - image transmission e.g. web data - multimedia applications e.g. desk-top editing Image data compression
More informationEnhanced Stochastic Bit Reshuffling for Fine Granular Scalable Video Coding
Enhanced Stochastic Bit Reshuffling for Fine Granular Scalable Video Coding Wen-Hsiao Peng, Tihao Chiang, Hsueh-Ming Hang, and Chen-Yi Lee National Chiao-Tung University 1001 Ta-Hsueh Rd., HsinChu 30010,
More informationIEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 52, NO. 2, FEBRUARY A Posteriori Quantization of Progressive Matching Pursuit Streams
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 52, NO. 2, FEBRUARY 2004 525 A Posteriori Quantization of Progressive Matching Pursuit Streams Pascal Frossard, Member, IEEE, Pierre Vandergheynst, Member,
More informationBases for time frequency analysis
4 Bases for time frequency analysis In the first part of this chapter we construct various bases having good time frequency localization. The construction of wavelet packets, for example, is presented
More informationRate Bounds on SSIM Index of Quantized Image DCT Coefficients
Rate Bounds on SSIM Index of Quantized Image DCT Coefficients Sumohana S. Channappayya, Alan C. Bovik, Robert W. Heath Jr. and Constantine Caramanis Dept. of Elec. & Comp. Engg.,The University of Texas
More informationAn Overview of Sparsity with Applications to Compression, Restoration, and Inverse Problems
An Overview of Sparsity with Applications to Compression, Restoration, and Inverse Problems Justin Romberg Georgia Tech, School of ECE ENS Winter School January 9, 2012 Lyon, France Applied and Computational
More informationContext Adaptive Space Quantization for Image Coding
Context Adaptive Space Quantization for Image Coding by Jeffrey Erbrecht A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Master of Applied Science
More informationLinear Independence of Finite Gabor Systems
Linear Independence of Finite Gabor Systems 1 Linear Independence of Finite Gabor Systems School of Mathematics Korea Institute for Advanced Study Linear Independence of Finite Gabor Systems 2 Short trip
More informationA Tutorial on Wavelets and their Applications. Martin J. Mohlenkamp
A Tutorial on Wavelets and their Applications Martin J. Mohlenkamp University of Colorado at Boulder Department of Applied Mathematics mjm@colorado.edu This tutorial is designed for people with little
More informationTree Partitioning Compression of Piecewise Polynomial Functions
Tree Partitioning ompression of Piecewise Polynomial Functions Arian Maleki Department of Electrical Engineering Stanford University arianm@stanford.edu Morteza Shahram Department of Statistics Stanford
More informationScalable color image coding with Matching Pursuit
SCHOOL OF ENGINEERING - STI SIGNAL PROCESSING INSTITUTE Rosa M. Figueras i Ventura CH-115 LAUSANNE Telephone: +4121 6935646 Telefax: +4121 69376 e-mail: rosa.figueras@epfl.ch ÉCOLE POLYTECHNIQUE FÉDÉRALE
More informationTHE newest video coding standard is known as H.264/AVC
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS FOR VIDEO TECHNOLOGY, VOL. 17, NO. 6, JUNE 2007 765 Transform-Domain Fast Sum of the Squared Difference Computation for H.264/AVC Rate-Distortion Optimization
More informationOptimizing Motion Vector Accuracy in Block-Based Video Coding
revised and resubmitted to IEEE Trans. Circuits and Systems for Video Technology, 2/00 Optimizing Motion Vector Accuracy in Block-Based Video Coding Jordi Ribas-Corbera Digital Media Division Microsoft
More informationWavelets: Theory and Applications. Somdatt Sharma
Wavelets: Theory and Applications Somdatt Sharma Department of Mathematics, Central University of Jammu, Jammu and Kashmir, India Email:somdattjammu@gmail.com Contents I 1 Representation of Functions 2
More information