DOCSLIB.ORG

Wavelet Transform and JPEG2000

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Wavelet Transform and JPEG2000 Wavelet Transform and JPEG2000 Yao Wang Polytechnic School of Engineering, New York University © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 1 Lecture Outline • Introduction • Multi-resolution representation of images: Gaussian and Laplacian pyramids • Wavelet transform through Iterated Filterbank Implementation • Basic Ideas in JPEG2000 Codec • JPEG vs. JPEG2000 • Scalability: why and how © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 2 Multi-Resolution Representation (aka Pyramid Representation) © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 3 Use of Pyramid Representations • Recall HBMA: uses pyramid representation to speed up motion estimation • Many other applications – Feature extraction across scales (SIFT) – Extracting faces of different sizes – … © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 4 © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 5 Pyramid is a redundant representation • A pyramid (either Guassian or Laplacian) includes an image of the original size plus additional smaller images • How many samples in all levels? • Original image (level J-0) NxN • Next level (level J-1): N/2xN/2 • Level J-l (l=0 to J): N/(2^l) x N/(2^l) • Total # samples =N^2 \sum_{l=0}^{J} (1/4^l)=N^2(1-(1/4)^(J-1))/ (1-1/4) \approx 4/3 N^2 • However, with Laplacian pyramid, many samples are close to zero except at the top level. • Wavelet transform is similar to Laplacian pyramid but does not incur additional pixels (non-redundant representation). © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 6 Wavelet vs. Pyramid vs. Subband Decomposition • Pyramid is a redundant transform (more samples than original) • Wavelet is a non-redundant multi-resolution representation • There are many ways to interpret wavelet transform. Here we describe the generation of discrete wavelet transform using the tree-structured subband decomposition (aka iterated filterbank) approach – 1D 2-band decomposition – 1D tree-structured subband decomposition (discrete wavelet transform) – Harr wavelet as an example – Extension to 2D by separable processing © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 7 Two Band Filterbank s(n) u(n) q(n) r(n) t(n) v(n) © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 8 What does the filter bank do? • When H0 and G0 are ideal Low-Pass and H1 and G1 are ideal High-Pass filters h0: Lowpass filter (0-¼), y0: a low-passed and then down-sampled version of x (Sampling theorem tells us we can down-sample after bandlimiting) h1: Highpass filter (1/4-1/2), y1: a high-passed and then down-sampled version of x (Sampling theorem also works in this case) g0: interpolation filter for low-pass subsignal g1: interpolation filter for high-pass subsignal • Can still reach perfect reconstruction even if these filters are not ideal low-pass/high-pass filters! – When the filters h0,h1, g0, g1 are designed appropriately, – X^=X (perfect reconstruction filterbank) © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 9 DTFT of signals after downsampling and upsampling Down&sampling!by!factor!of!2: 1⎛ ⎛ u⎞ ⎛ u 1⎞ ⎞ x (m) x(2m) X (u) X X d = ⇔ d = ⎜ ⎜ ⎟ + ⎜ − ⎟ ⎟ 2⎝ ⎝ 2⎠ ⎝ 2 2⎠ ⎠ Up&sampling!by!factor!of!2: ⎧ ⎛ m⎞ ⎪ x⎜ ⎟ , m = even xu(m)= ⎨ ⎝ 2 ⎠ ⇔ Xu(u)= X(2u) ⎪ 0, otherwise ! ⎩⎪ Conceptual proof by doing sampling on continuous signal under 2 different sampling rates. © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 10 Perfect Reconstruction Conditions x(n)*h0(n)⇔ X(u)H0(u) xl(n)= down(x(n)*h0(n))⇔ Xl(u)= (X(u/2)H0(u/2)+ X(u/2−1/2)H0(u/2−1/2))/2 up(xl(n))⇔ Xl(2u)= (X(u)H0(u)+ X(u−1)H0(u−1))/2 x(n)= up(x (n))* g (n)+up(x (n))* g (n) l 0 h 1 ⇔ X(u)= (X(u)H0(u)G0(u)+ X(u−1)H0(u−1)G0(u))/2 +(X(u)H1(u)G1(u)+ X(u−1)H1(u−1)G1(u))/2 = X(u)(H0(u)G0(u)+ H1(u)G1(u))/2 + X(u−1)(H (u−1)G (u)+ H (u−1)G (u))/2 0 0 1 1 To!guarantee!X(u)= X(u),!we!need H0(u)G0(u)+ H1(u)G1(u)= 2 H (u−1)G (u)+ H (u−1)G (u)= 0!!(To!remove!aliasing!component!) ! 0 0 1 1 © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 11 Perfect!reconstruction!condition: H (u)G (u) H (u)G (u) 2 0 0 + 1 1 = Perfect reconstruction H (u−1)G (u)+ H (u−1)G (u)= 0 0 0 1 1 conditions The!second!equation!(aliasing!cancelation)!can!be!guaranteed!by!requiring n G (u) H (u 1) g (n) 1 h (n) 0 = 1 − ⇔ 0 = (− ) 1 n+1 G (u) H (u 1) g (n) 1 h (n) 1 = − 0 − ⇔ 1 = (− ) 0 (Quadrature!Mirror!Condition) To!guarantee!perfect!reconstruction,!the!filters!must!satisfy!the!biorthogonality!condition: < hi(2n− k), gj (k)>= δ(i − j)δ(n) One!has!freedom!to!design!both!g0(n), g1(n),!which!can!have!different!length. A!more!strict!condition!requires!orthonality!between!g0(n), g1(n)!: < gi(n), gj (n+2m)>= δ(i − j)δ(m) which!yields n g (n) 1 g (L 1 n), 1 = (− ) 0 − − h0(n)= g0(L−1− n) n n h (n) g (L 1 n) 1 g (n) 1 h (L 1 n) 1 = 1 − − = (− ) 0 = (− ) 0 − − One!only!has!freedom!to!design!g (n),!filter!length!L!must!be!even!and!all!filters!have!same!length. ! 0 © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 12 Haar Filter (Simplest Orthogonal Wavelet Filter) h0 : averaging, [1,1]/ 2; h1: difference, [1,-1]/ 2; g0 = [1,1]/ 2; g1 = [-1,1]/ 2 Input sequence : [x1,x2, x3,x4,....] Analysis(Assuming samples outside the boundaries are 0. remember to flip the filter when doing convolution) s = x *h0 = [s0,s1, s2, s3, s4,...], s0 = (x1+ 0)/ 2,s1 = (x2+ x1)/ 2,s2 = (x3+ x2)/ 2,s3 = (x4+ x3)/ 2... y0 = s ↓ 2 = [s1, s3,...,] t = x *h1 = [t0, t1,t2,t3,t4,...], t0 = [x1- 0]/ 2, t1 = [x2- x1]/ 2, t2 = [x3- x2]/ 2, t3 = [x4- x3]/ 2,... y1= t ↓ 2 = [t1, t3,...] Synthesis: u = y0 ↑ 2 = [0,s1,0,s3,...] q = u *g0 = [q1,q2,q3,q4,...],q1 = (s1+ 0)/ 2 = (x1+ x2)/2, q2 = (0+ s1)/ 2 = (x1+ x2)/2, q3 = (s3+ 0)/ 2 = (x3+ x4)/2 v = y1↑ 2 = [0,t1,0, t3,...] r = v*g1 = [r1,r2,r3,r4,...],r1= (-t1+ 0)/ 2 = (x1− x2)/2, r2 = (-0+ t1)/ 2 = (-x1+ x2)/2, r3 = (-t3+ 0)/ 2 = (x3− x4)/2, xˆ = q + r = [q1+ r1,q2 + r2,....] = [x1,x2, x3,...] Note with Haar wavelet, the lowpass subband essentially takes the average of every two samples, L=(x1+x2)/sqrt(2), and the highpass subband takes the difference of every two samples, H=(x1-x2)/ sqrt(2). For synthesis, you take the sum of the lowpass and high pass signal to recover first sample A=(L +H)/sqrt(2), and you take the difference to recover the second sample B=(L-H)/sqrt(2). © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 13 Common Wavelet Filters • Haar: simplest, orthogonal, not very good in energy compaction • Daubechies 8/8: orthogonal • Daubechies 9/7: bi-orthogonal, most commonly used if numerical reconstruction errors are acceptable • LeGall 5/3: bi-orthogonal, integer operation, can be implemented with integer operations only, used for lossless image coding • Differ in energy compaction capability © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 14 © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 15 MATLAB example >> [ca,cd]=dwt(y,'db4'); >> z=idwt(ca,cd,'db4'); >> wy=[ca,cd]; >> subplot(3,1,1),plot(y), title('Original sequence'); >> subplot(3,1,2),plot(wy), title('Wavelet transform: left=low band, righ=high band'); >> subplot(3,1,3),plot(y), title('Reconstructed sequence'); © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 16 Iterated Filter Bank From [Vetterli01] © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 17 MATLAB example >> [ca,cd]=dwt(y,'db4'); [caa,cad]=dwt(ca,'db4'); [caaa,caad]=dwt(caa,'db4'); wy1=[ca,cd]; wy2=[caa,cad,cd]; >> wy3=[caaa,caad,cad,cd]; >> subplot(4,1,1),plot(y),title('Or iginal Signal'); >> subplot(4,1,2),plot(wy1),title(' 1-level Wavelet Transform'); >> subplot(4,1,3),plot(wy2),title(' 2-level Wavelet Transform'); >> subplot(4,1,4),plot(wy3),title(' 3-level Wavelet Transform'); © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 18 Discrete Wavelet Transform = Iterated Filter Bank © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 19 Temporal-Frequency Domain Partition © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 20 Wavelet Transform vs. Fourier Transform • Fourier transform: – Basis functions cover the entire signal range, varying in frequency only • Wavelet transform – Basis functions vary in frequency (called “scale”) as well as spatial extend • High frequency basis covers a smaller area • Low frequency basis covers a larger area • Non-uniform partition of frequency range and spatial range • More appropriate for non-stationary signals © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 21 Haar Wavelet: Analysis Note that the assumed high pass filter in this example has a factor “-1” difference from our previous example. Similarly the synthesis filter is off by the same factor. Both are OK. © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 22 Haar Wavelet: Synthesis © Yao Wang, 2016 EL-GY 6123: Image and Video Processing 23 How to Apply Filterbank to Images? 2D decomposition is accomplished by applying the 1D decomposition along rows of an image first, and then columns.
Load more

Recommended publications
  • An Overview of Wavelet Transform Concepts and Applications
    An overview of wavelet transform concepts and applications Christopher Liner, University of Houston February 26, 2010 Abstract The continuous wavelet transform utilizing a complex Morlet analyzing wavelet has a close connection to the Fourier transform and is a powerful analysis tool for decomposing broadband wavefield data. A wide range of seismic wavelet applications have been reported over the last three decades, and the free Seismic Unix processing system now contains a code (succwt) based on the work reported here. Introduction The continuous wavelet transform (CWT) is one method of investigating the time-frequency details of data whose spectral content varies with time (non-stationary time series). Moti- vation for the CWT can be found in Goupillaud et al. [12], along with a discussion of its relationship to the Fourier and Gabor transforms. As a brief overview, we note that French geophysicist J. Morlet worked with non- stationary time series in the late 1970's to find an alternative to the short-time Fourier transform (STFT). The STFT was known to have poor localization in both time and fre- quency, although it was a first step beyond the standard Fourier transform in the analysis of such data. Morlet's original wavelet transform idea was developed in collaboration with the- oretical physicist A. Grossmann, whose contributions included an exact inversion formula. A series of fundamental papers flowed from this collaboration [16, 12, 13], and connections were soon recognized between Morlet's wavelet transform and earlier methods, including harmonic analysis, scale-space representations, and conjugated quadrature filters. For fur- ther details, the interested reader is referred to Daubechies' [7] account of the early history of the wavelet transform.
    [Show full text]
  • A Low Bit Rate Audio Codec Using Wavelet Transform
    Navpreet Singh, Mandeep Kaur,Rajveer Kaur / International Journal of Engineering Research and Applications (IJERA) ISSN: 2248-9622 www.ijera.com Vol. 3, Issue 4, Jul-Aug 2013, pp.2222-2228 An Enhanced Low Bit Rate Audio Codec Using Discrete Wavelet Transform Navpreet Singh1, Mandeep Kaur2, Rajveer Kaur3 1,2(M. tech Students, Department of ECE, Guru kashi University, Talwandi Sabo(BTI.), Punjab, INDIA 3(Asst. Prof. Department of ECE, Guru Kashi University, Talwandi Sabo (BTI.), Punjab,INDIA Abstract Audio coding is the technology to represent information and perceptually irrelevant signal audio in digital form with as few bits as possible components can be separated and later removed. This while maintaining the intelligibility and quality class includes techniques such as subband coding; required for particular application. Interest in audio transform coding, critical band analysis, and masking coding is motivated by the evolution to digital effects. The second class takes advantage of the communications and the requirement to minimize statistical redundancy in audio signal and applies bit rate, and hence conserve bandwidth. There is some form of digital encoding. Examples of this class always a tradeoff between lowering the bit rate and include entropy coding in lossless compression and maintaining the delivered audio quality and scalar/vector quantization in lossy compression [1]. intelligibility. The wavelet transform has proven to Digital audio compression allows the be a valuable tool in many application areas for efficient storage and transmission of audio data. The analysis of nonstationary signals such as image and various audio compression techniques offer different audio signals. In this paper a low bit rate audio levels of complexity, compressed audio quality, and codec algorithm using wavelet transform has been amount of data compression.
    [Show full text]
  • Image Compression Techniques by Using Wavelet Transform
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by International Institute for Science, Technology and Education (IISTE): E-Journals Journal of Information Engineering and Applications www.iiste.org ISSN 2224-5782 (print) ISSN 2225-0506 (online) Vol 2, No.5, 2012 Image Compression Techniques by using Wavelet Transform V. V. Sunil Kumar 1* M. Indra Sena Reddy 2 1. Dept. of CSE, PBR Visvodaya Institute of Tech & Science, Kavali, Nellore (Dt), AP, INDIA 2. School of Computer science & Enginering, RGM College of Engineering & Tech. Nandyal, A.P, India * E-mail of the corresponding author: [email protected] Abstract This paper is concerned with a certain type of compression techniques by using wavelet transforms. Wavelets are used to characterize a complex pattern as a series of simple patterns and coefficients that, when multiplied and summed, reproduce the original pattern. The data compression schemes can be divided into lossless and lossy compression. Lossy compression generally provides much higher compression than lossless compression. Wavelets are a class of functions used to localize a given signal in both space and scaling domains. A MinImage was originally created to test one type of wavelet and the additional functionality was added to Image to support other wavelet types, and the EZW coding algorithm was implemented to achieve better compression. Keywords: Wavelet Transforms, Image Compression, Lossless Compression, Lossy Compression 1. Introduction Digital images are widely used in computer applications. Uncompressed digital images require considerable storage capacity and transmission bandwidth. Efficient image compression solutions are becoming more critical with the recent growth of data intensive, multimedia based web applications.
    [Show full text]
  • JPEG2000: Wavelets in Image Compression
    EE678 WAVELETS APPLICATION ASSIGNMENT 1 JPEG2000: Wavelets In Image Compression Group Members: Qutubuddin Saifee [email protected] 01d07009 Ankur Gupta [email protected] 01d070013 Nishant Singh [email protected] 01d07019 Abstract During the past decade, with the birth of wavelet theory and multiresolution analysis, image processing techniques based on wavelet transform have been extensively studied and tremendously improved. JPEG 2000 uses wavelet transform and provides an integrated toolbox to better address increasing needs for compression. In this report, we study the basic concepts of JPEG2000, the LeGall 5/3 and Daubechies 9/7 wavelets used in it and finally Embedded zerotree wavelet coding and Set Partitioning in Hierarchical Trees. Index Terms Wavelets, JPEG2000, Image Compression, LeGall, Daubechies, EZW, SPIHT. I. INTRODUCTION I NCE the mid 1980s, members from both the International Telecommunications Union (ITU) and the International SOrganization for Standardization (ISO) have been working together to establish a joint international standard for the compression of grayscale and color still images. This effort has been known as JPEG, the Joint Photographic Experts Group. The process was such that, after evaluating a number of coding schemes, the JPEG members selected a discrete cosine transform( DCT)-based method in 1988. From 1988 to 1990, the JPEG group continued its work by simulating, testing and documenting the algorithm. JPEG became Draft International Standard (DIS) in 1991 and International Standard (IS) in 1992. With the continual expansion of multimedia and Internet applications, the needs and requirements of the technologies used grew and evolved. In March 1997, a new call for contributions was launched for the development of a new standard for the compression of still images, the JPEG2000 standard.
    [Show full text]
  • Comparison of Image Compressions: Analog Transformations P
    Proceedings Comparison of Image Compressions: Analog † Transformations P Jose Balsa P CITIC Research Center, Universidade da Coruña (University of A Coruña), 15071 A Coruña, Spain; [email protected] † Presented at the 3rd XoveTIC Conference, A Coruña, Spain, 8–9 October 2020. Published: 21 August 2020 Abstract: A comparison between the four most used transforms, the discrete Fourier transform (DFT), discrete cosine transform (DCT), the Walsh–Hadamard transform (WHT) and the Haar- wavelet transform (DWT), for the transmission of analog images, varying their compression and comparing their quality, is presented. Additionally, performance tests are done for different levels of white Gaussian additive noise. Keywords: analog image transformation; analog image compression; analog image quality 1. Introduction Digitized image coding systems employ reversible mathematical transformations. These transformations change values and function domains in order to rearrange information in a way that condenses information important to human vision [1]. In the new domain, it is possible to filter out relevant information and discard information that is irrelevant or of lesser importance for image quality [2]. Both digital and analog systems use the same transformations in source coding. Some examples of digital systems that employ these transformations are JPEG, M-JPEG, JPEG2000, MPEG- 1, 2, 3 and 4, DV and HDV, among others. Although digital systems after transformation and filtering make use of digital lossless compression techniques, such as Huffman. In this work, we aim to make a comparison of the most commonly used transformations in state- of-the-art image compression systems. Typically, the transformations used to compress analog images work either on the entire image or on regions of the image.
    [Show full text]
  • Wavelet Transform for JPG, BMP & TIFF
    Wavelet Transform for JPG, BMP & TIFF Samir H. Abdul-Jauwad Mansour I. AlJaroudi Electrical Engineering Department, Electrical Engineering Department, King Fahd University of Petroleum & King Fahd University of Petroleum & Minerals Minerals Dhahran 31261, Saudi Arabia Dhahran 31261, Saudi Arabia ABSTRACT This paper presents briefly both wavelet transform and inverse wavelet transform for three different images format; JPG, BMP and TIFF. After brief historical view of wavelet, an introduction will define wavelets transform. By narrowing the scope, it emphasizes the discrete wavelet transform (DWT). A practical example of DWT is shown, by choosing three images and applying MATLAB code for 1st level and 2nd level wavelet decomposition then getting the inverse wavelet transform for the three images. KEYWORDS: Wavelet Transform, JPG, BMP, TIFF, MATLAB. __________________________________________________________________ Correspondence: Dr. Samir H. Abdul-Jauwad, Electrical Engineering Department, King Fahd University of Petroleum & Mineral, Dhahran 1261, Saudi Arabia. Email: [email protected] Phone: +96638602500 Fax: +9663860245 INTRODUCTION In 1807, theories of frequency analysis by Joseph Fourier were the main lead to wavelet. However, wavelet was first appeared in an appendix to the thesis of A. Haar in 1909 [1]. Compact support was one property of the Haar wavelet which means that it vanishes outside of a finite interval. Unfortunately, Haar wavelets have some limits, because they are not continuously differentiable. In 1930, work done separately by scientists for representing the functions by using scale-varying basis functions was the key to understanding wavelets. In 1980, Grossman and Morlet, a physicist and an engineer, provided a way of thinking for wavelets based on physical intuition through defining wavelets in the context of quantum physics [3].
    [Show full text]
  • The Wavelet Tutorial Second Edition Part I by Robi Polikar
    THE WAVELET TUTORIAL SECOND EDITION PART I BY ROBI POLIKAR FUNDAMENTAL CONCEPTS & AN OVERVIEW OF THE WAVELET THEORY Welcome to this introductory tutorial on wavelet transforms. The wavelet transform is a relatively new concept (about 10 years old), but yet there are quite a few articles and books written on them. However, most of these books and articles are written by math people, for the other math people; still most of the math people don't know what the other math people are 1 talking about (a math professor of mine made this confession). In other words, majority of the literature available on wavelet transforms are of little help, if any, to those who are new to this subject (this is my personal opinion). When I first started working on wavelet transforms I have struggled for many hours and days to figure out what was going on in this mysterious world of wavelet transforms, due to the lack of introductory level text(s) in this subject. Therefore, I have decided to write this tutorial for the ones who are new to the topic. I consider myself quite new to the subject too, and I have to confess that I have not figured out all the theoretical details yet. However, as far as the engineering applications are concerned, I think all the theoretical details are not necessarily necessary (!). In this tutorial I will try to give basic principles underlying the wavelet theory. The proofs of the theorems and related equations will not be given in this tutorial due to the simple assumption that the intended readers of this tutorial do not need them at this time.
    [Show full text]
  • Discrete Generalizations of the Nyquist-Shannon Sampling Theorem
    DISCRETE SAMPLING: DISCRETE GENERALIZATIONS OF THE NYQUIST-SHANNON SAMPLING THEOREM A DISSERTATION SUBMITTED TO THE DEPARTMENT OF ELECTRICAL ENGINEERING AND THE COMMITTEE ON GRADUATE STUDIES OF STANFORD UNIVERSITY IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY William Wu June 2010 © 2010 by William David Wu. All Rights Reserved. Re-distributed by Stanford University under license with the author. This work is licensed under a Creative Commons Attribution- Noncommercial 3.0 United States License. http://creativecommons.org/licenses/by-nc/3.0/us/ This dissertation is online at: http://purl.stanford.edu/rj968hd9244 ii I certify that I have read this dissertation and that, in my opinion, it is fully adequate in scope and quality as a dissertation for the degree of Doctor of Philosophy. Brad Osgood, Primary Adviser I certify that I have read this dissertation and that, in my opinion, it is fully adequate in scope and quality as a dissertation for the degree of Doctor of Philosophy. Thomas Cover I certify that I have read this dissertation and that, in my opinion, it is fully adequate in scope and quality as a dissertation for the degree of Doctor of Philosophy. John Gill, III Approved for the Stanford University Committee on Graduate Studies. Patricia J. Gumport, Vice Provost Graduate Education This signature page was generated electronically upon submission of this dissertation in electronic format. An original signed hard copy of the signature page is on file in University Archives. iii Abstract 80, 0< 11 17 æ 80, 1< 5 81, 0< 81, 1< 12 83, 1< 84, 0< 82, 1< 83, 4< 82, 3< 15 7 6 83, 2< 83, 3< 84, 8< 83, 6< 83, 7< 13 3 æ æ æ æ æ æ æ æ æ æ 9 -5 -4 -3 -2 -1 1 2 3 4 5 84, 4< 84, 13< 84, 14< 1 HIS DISSERTATION lays a foundation for the theory of reconstructing signals from T finite-dimensional vector spaces using a minimal number of samples.
    [Show full text]
  • A Novel Robust Audio Watermarking Algorithm by Modifying the Average Amplitude in Transform Domain
    applied sciences Article A Novel Robust Audio Watermarking Algorithm by Modifying the Average Amplitude in Transform Domain Qiuling Wu 1,2,* and Meng Wu 1 1 College of Telecommunications and Information Engineering, Nanjing University of Posts and Telecommunications, Nanjing 210003, China; [email protected] 2 College of Zijin, Nanjing University of Science and Technology, Nanjing 210046, China * Correspondence: [email protected]; Tel.: +86-139-5103-2810 Received: 18 April 2018; Accepted: 1 May 2018; Published: 4 May 2018 Featured Application: This algorithm embeds a binary image into an audio signal as a marker to prove the ownership of this audio signal. With large payload capacity and strong robustness against common signal processing attacks, it can be used for copyright protection, broadcast monitoring, fingerprinting, data authentication, and medical safety. Abstract: In order to improve the robustness and imperceptibility in practical application, a novel audio watermarking algorithm with strong robustness is proposed by exploring the multi-resolution characteristic of discrete wavelet transform (DWT) and the energy compaction capability of discrete cosine transform (DCT). The human auditory system is insensitive to the minor changes in the frequency components of the audio signal, so the watermarks can be embedded by slightly modifying the frequency components of the audio signal. The audio fragments segmented from the cover audio signal are decomposed by DWT to obtain several groups of wavelet coefficients with different frequency bands, and then the fourth level detail coefficient is selected to be divided into the former packet and the latter packet, which are executed for DCT to get two sets of transform domain coefficients (TDC) respectively.
    [Show full text]
  • A Very Short Introduction to Sound Analysis for Those Who Like Elephant Trumpet Calls Or Other Wildlife Sound
    A very short introduction to sound analysis for those who like elephant trumpet calls or other wildlife sound J´er^ome Sueur Mus´eum national d'Histoire naturelle CNRS UMR 7205 ISYEB, Paris, France July 14, 2021 This document is a very brief introduction to sound analysis principles. It is mainly written for students starting with bioacoustics. The content should be updated regularly. Demonstrations are based on the package seewave. For more details see the book https://www.springer.com/us/book/9783319776453. Contents 1 Digitization 3 1.1 Sampling........................................3 1.2 Quantisation......................................3 1.3 File format.......................................4 2 Amplitude envelope4 3 Discrete Fourier Transform (DFT)6 3.1 Definitions and principle................................6 3.2 Complete sound.....................................7 3.3 Sound section......................................7 3.3.1 Window shape................................. 10 4 Short-time Fourier Transform (STFT) 10 4.1 Principle......................................... 10 4.2 Spectrogram...................................... 11 4.2.1 3D in a 2D plot................................. 11 4.2.2 Overlap..................................... 12 4.2.3 Values...................................... 13 4.3 Mean spectrum..................................... 13 5 Instantaneous frequency 15 5.1 Zero-crossing...................................... 15 5.2 Hilbert transform.................................... 15 5.3 Cepstral transform..................................
    [Show full text]
  • Image Compression Using the Haar Wavelet Transform
    Image compression using the Haar wavelet transform Colm Mulcahy, Ph.D. ABSTRACT The wavelet transform is a relatively new arrival on the mathematical scene. Wegive a brief intro duction to the sub ject byshowing how the Haar wavelet transform allows information to b e enco ded according to \levels of detail." In one sense, this parallels the way in whichwe often pro cess information in our everydaylives. Our investigations enable us to givetwo interesting applications of wavelet metho ds to digital images: compression and progressive transmission. The mathematical prerequesites will b e kept to a minimum; indeed, the main concepts can b e understo o d in terms of addition, subtraction and division bytwo. We also present a linear algebra implementation of the Haar wavelet transform, and mention imp ortantrecent generalizations. This material is suitable for undergraduate researchby b oth mathematics and computer science ma jors, and has b een used successfully in this way at Sp elman College since 1995. 1 So much information In the real world, we are constantly faced with having to decide howmuch information to absorb or impart. Example 1: An election is imp ending, and you and your b est friend haveyet to decide for whom you will cast your votes. Yourfriendiscontent to base her decision solely on the p olitical party aliations of the available candidates, whereas you take the trouble to nd out how the candidates stand on a variety of issues, and to investigate their records for integrity, e ectiveness, and vision. Example 2: Your younger brother is across the ro om, working on a pro ject for Black History Month.
    [Show full text]
  • A Tutorial of the Wavelet Transform
    A Tutorial of the Wavelet Transform Chun-Lin, Liu February 23, 2010 Chapter 1 Overview 1.1 Introduction The Fourier transform is an useful tool to analyze the frequency components of the signal. However, if we take the Fourier transform over the whole time axis, we cannot tell at what instant a particular frequency rises. Short-time Fourier transform (STFT) uses a sliding window to find spectrogram, which gives the information of both time and frequency. But still another problem exists: The length of window limits the resolution in frequency. Wavelet transform seems to be a solution to the problem above. Wavelet transforms are based on small wavelets with limited duration. The translated-version wavelets locate where we concern. Whereas the scaled-version wavelets allow us to analyze the signal in different scale. 1.2 History The first literature that relates to the wavelet transform is Haar wavelet. It was proposed by the mathematician Alfrd Haar in 1909. However, the con- cept of the wavelet did not exist at that time. Until 1981, the concept was proposed by the geophysicist Jean Morlet. Afterward, Morlet and the physi- cist Alex Grossman invented the term wavelet in 1984. Before 1985, Haar wavelet was the only orthogonal wavelet people know. A lot of researchers even thought that there was no orthogonal wavelet except Haar wavelet. For- tunately, the mathematician Yves Meyer constructed the second orthogonal wavelet called Meyer wavelet in 1985. As more and more scholars joined in this field, the 1st international conference was held in France in 1987. 1 In 1988, Stephane Mallat and Meyer proposed the concept of multireso- lution.
    [Show full text]