Wavelet Analysis Basic Concepts and Applications 1st Edition Sabrine Arfaoui Download
Wavelet Analysis Basic Concepts and Applications 1st Edition Sabrine Arfaoui Download
https://ebookbell.com/product/wavelet-analysis-basic-concepts-
and-applications-1st-edition-sabrine-arfaoui-51684912
https://ebookbell.com/product/wavelet-analysis-on-the-sphere-
spheroidal-wavelets-sabrine-arfaoui-imen-rezgui-anouar-ben-mabrouk-
knowledge-unlatched-51110958
https://ebookbell.com/product/wavelet-analysis-and-transient-signal-
processing-applications-for-power-systems-he-22035766
https://ebookbell.com/product/wavelet-analysis-on-local-fields-of-
positive-characteristic-biswaranjan-behera-37754402
Wavelet Analysis In Civil Engineering Pranesh Chatterjee
https://ebookbell.com/product/wavelet-analysis-in-civil-engineering-
pranesh-chatterjee-4986866
https://ebookbell.com/product/wavelet-analysis-in-civil-engineering-
pranesh-chatterjee-5059264
https://ebookbell.com/product/wavelet-analysis-and-its-applications-
second-international-conference-waa-2001-hong-kong-china-
december-1820-2001-proceedings-1st-edition-john-daugman-auth-1147700
Wavelet Analysis And Applications 1st Edition Tao Qian Tao Qianmang I
Vaixu Yuesheng
https://ebookbell.com/product/wavelet-analysis-and-applications-1st-
edition-tao-qian-tao-qianmang-i-vaixu-yuesheng-1189502
Sabrine Arfaoui
University of Monastir, Tunisia
University of Tabuk, Saudi Arabia
Carlo Cattani
University of Tuscia, Italy
First edition published 2021
by CRC Press
6000 Broken Sound Parkway NW, Suite 300, Boca Raton, FL 33487-2742
The right of Sabrine Arfaoui, Anouar Ben Mabrouk, and Carlo Cattani to be identified as authors of this work has been
asserted by them in accordance with sections 77 and 78 of the Copyright, Designs and Patents Act 1988.
Reasonable efforts have been made to publish reliable data and information, but the author and publisher cannot
assume responsibility for the validity of all materials or the consequences of their use. The authors and publishers
have attempted to trace the copyright holders of all material reproduced in this publication and apologize to copyright
holders if permission to publish in this form has not been obtained. If any copyright material has not been acknowl-
edged please write and let us know so we may rectify in any future reprint.
Except as permitted under U.S. Copyright Law, no part of this book may be reprinted, reproduced, transmitted, or
utilized in any form by any electronic, mechanical, or other means, now known or hereafter invented, including pho-
tocopying, microfilming, and recording, or in any information storage or retrieval system, without written permission
from the publishers.
For permission to photocopy or use material electronically from this work, access www.copyright.com or contact the
Copyright Clearance Center, Inc. (CCC), 222 Rosewood Drive, Danvers, MA 01923, 978-750-8400. For works that are
not available on CCC please contact [email protected]
Trademark notice: Product or corporate names may be trademarks or registered trademarks and are used only for
identification and explanation without intent to infringe.
List of Figures ix
Preface xi
Chapter 1 Introduction 1
2.1 INTRODUCTION 5
2.2 WAVELETS ON R 6
2.2.1 Continuous wavelet transform 7
2.2.2 Discrete wavelet transform 10
2.3 MULTI-RESOLUTION ANALYSIS 11
2.4 WAVELET ALGORITHMS 13
2.5 WAVELET BASIS 16
2.6 MULTIDIMENSIONAL REAL WAVELETS 21
2.7 EXAMPLES OF WAVELET FUNCTIONS AND MRA 22
2.7.1 Haar wavelet 22
2.7.2 Faber–Schauder wavelet 24
2.7.3 Daubechies wavelets 25
2.7.4 Symlet wavelets 27
2.7.5 Spline wavelets 27
2.7.6 Anisotropic wavelets 29
2.7.7 Cauchy wavelets 30
2.8 EXERCISES 31
v
vi Contents
4.1 INTRODUCTION 51
4.2 DIFFERENT CONSTRUCTIONS OF CLIFFORD ALGEBRAS 52
4.2.1 Clifford original construction 53
4.2.2 Quadratic form-based construction 53
4.2.3 A standard construction 54
4.3 GRADUATION IN CLIFFORD ALGEBRAS 56
4.4 SOME USEFUL OPERATIONS ON CLIFFORD ALGEBRAS 57
4.4.1 Products in Clifford algebras 57
4.4.2 Involutions on a Clifford algebra 58
4.5 CLIFFORD FUNCTIONAL ANALYSIS 60
4.6 EXISTENCE OF MONOGENIC EXTENSIONS 67
4.7 CLIFFORD-FOURIER TRANSFORM 70
4.8 CLIFFORD WAVELET ANALYSIS 76
4.8.1 Spin-group based Clifford wavelets 76
4.8.2 Monogenic polynomial-based Clifford wavelets 82
4.9 SOME EXPERIMENTATIONS 92
4.10 EXERCISES 96
5.1 INTRODUCTION 99
5.2 BESSEL FUNCTIONS 99
5.3 BESSEL WAVELETS 105
5.4 FRACTIONAL BESSEL WAVELETS 107
5.5 QUANTUM THEORY TOOLKIT 119
Contents vii
Bibliography 209
Index 237
List of Figures
6.1 The normal reduced and centered density N (0, 1). 158
6.2 The gaussian density and its wavelet estimator. 159
6.3 Claw density. 159
6.4 Claw density wavelet estimator at the level J = 4. 160
ix
Preface
Nowadays, wavelets are applied almost everywhere in science. Both pure fields, such
as mathematics and theoretical physics, and applied ones, such as signal/image pro-
cessing, finance and engineering, apply wavelets. Although the references and/or the
documentation about wavelets and their applications are wide, it seems that with the
advancement of technology and the appearance of many phenomena in nature and in
life there still exist some places for more efforts and developments to understand the
new problems, as the existing wavelet methods do not provide good understanding
of them. The new COVID-19 pandemic may be one of the challenges that should be
understood.
On the other hand, especially for young researchers, existing references such as
books in wavelet theory are somehow very restricted. The majority are written for
specific communities. This is, in fact, not surprising and may be due to the necessity
of developing such references to overcome the concerned problems in that time.
Next, with the inclusion of wavelet theory in academic studies such as in mas-
ter’s and PhD programs, the scientific and academic communities have had a great
need to develop references in other forms. Students and generally researchers need
sometimes self-containing references responding to their need, to avoid losing time
in redeveloping existing results, which is a necessary step for both the generalization
and the experiments.
The present volume is composed of eight chapters. In the first introductory chap-
ter, a literal introduction is developed discussing generally the topic. Chapter 2 is
concerned with the presentation of the original developments of wavelet theory on
the real Euclidean space. This is also a preliminary chapter that will be of great help
for young researchers. Chapter 3 is more specialized and constitutes a continuation of
the previous one, in which some extending cases of wavelet theory and applications
have been provided. Chapter 4 is a very specialized part that is developed for the
first time to our knowledge. It is concerned with the presentation of wavelet theory in
a general functional framework based on Clifford algebras. This is very important as
these algebras contain all the Euclidean structures and gather them in one structure
to facilitate calculus. Readers will notice clearly that Clifford wavelet theory induces
naturally the Euclidean ones such as real and complex numbers, circles and spheres.
Chapter 5 is a continuation of the development of the theory in specialized fields
such as quantum theory. Next, in Chapter 6, statistical application of wavelets has
been reviewed. Topics such as density estimation, thresholding concepts, variance and
covariance have been detailed. Chapter 7 is devoted to wavelets applied in solving
partial differential equations. Recall that this field needs many assumptions on the
functional bases applied, especially the explicit form of the basis elements and their
xi
xii Preface
regularities. The last chapter is devoted to the link and/or the use of wavelet theory
in characterizing fractal and multifractal functions and their application. Each chap-
ter contains a series of exercises and experimentations to help understand the theory
and also to show the utility of wavelets.
The present book stems, in fact, from lectures and papers on the topics developed,
which have been gathered, re-developed, improved and sometimes completed with
necessary missing developments. However, naturally it is not exhaustive and should
be always criticized, sometimes corrected and improved by readers. So, we accept and
wait for any comments and suggestions.
We also want to stress the fact that we have provided in some chapters, especially
those on preliminary concepts that may be useful to young researchers, some exercises
and applications that are simple to handle with the aim to help the readers understand
the theory. We apologize if there are simpler applications and details that may be
more helpful to the readers but that have been left out from inclusion in this book.
This, in fact, needs more time and may induce delays in the publication of the book.
We hope that with the present form the readers become acquainted with the topics
presented.
The aim of this book is to provide a basic and self-contained introduction to
the ideas underpinning wavelet theory and its diversified applications. Readers of
our proposed book would include master’s degree students, PhD students and se-
nior researchers. It may also serve scientists and research workers from industrial
settings, where modeling real-world phenomena and data needs wavelets such as fi-
nance, medicine, engineering, transport, images and signals. Henceforth, the book
will interest practitioners and theorists alike. For theorists, rigorous mathematical
developments will be presented with necessary prerequisites that make the book self-
containing. For the practitioner, often interested in model building and analysis, we
provide the cornerstone ideas.
As with any scientific production and reference, the present volume could not have
been realized without the help of many persons. We thus owe thanks to many persons
who have helped us in any direction such as encouragements, scientific discussions
and documentation. We thank the Taylor & Francis Publishing Group for giving us
the opportunity to write and publish the present work. We also would like to express
our gratitude to our professors, teachers, colleagues, and universities. Without their
help and efforts, no such work might be realized. We would also like to thank all the
members of the publishing house, especially the editorial staff for the present volume,
Callum Fraser and Mansi Kabra, for their hospitality, cooperation, collaboration and
for the time they have spent on our project.
CHAPTER 1
Introduction
1
2 Wavelet Analysis: Basic Concepts and Applications
the so-called isotropic and anisotropic wavelets (see [11], [12], [20], [205], [218], [219],
[220], [229], [230], [248], [279], [295], [305], [357], [358], [363]).
Wavelet theory provides for functional spaces and time series good bases, allow-
ing their decomposition into spaces associated with different horizons known as the
levels of decomposition. A wavelet basis is a family of functions obtained from one
function known as the mother wavelet, by translations and dilations. Due to the
power of their theory, wavelets have many applications in different domains such as
mathematics, physics, electrical engineering and seismic geology. This tool permits
the representation of L2 -functions in a basis well localized in time and in frequency.
Wavelets are also associated with many special functions such as orthogonal poly-
nomials and hypergeometric series. The most well known may be the Bessel functions
that have been developed in both classic theory of Bessel functional analysis and the
modified versions in fractional and quantum calculus. As its name indicates, Bessel
wavelets are related to Bessel special function. Historically, special functions differ
from elementary ones such as powers, roots, trigonometric and their inverses, mainly
with the limitations that these latter classes are known for. Many fundamental prob-
lems such as orbital motion, simultaneous oscillatory chains and spherical body grav-
itational potential were not best described using elementary functions. This makes
it necessary to extend elementary functions’ classes to more general ones that may
describe well unresolved problems.
Wavelets are also developed and applied in financial time series such as market
indices and exchange rates. In [42], for example, a study of the largest transaction
financial market was carried out. The exchange market gave some high-frequency
data. Compared to other markets, such data can be available at long periods and
with high frequency. The data were detected for very small periods, which means
that the market is also liquid. Until 1990, economists were interested in intra-daily
data because of which the detection of some behaviors did not appear in the daily
analysis of data such as homogeneity.
A well-known hypothesis in finance is the homogeneity of markets where all in-
vestigators have almost the same behavior. The idea of nonhomogeneous markets is
more recent, and it suggests that investigators have different perceptions and differ-
ent laws. For the exchange market, for example, investigators can differ in profiles,
geographic localizations and also in institutional constraints. Another natural sug-
gestion can be done about traders. Naturally, traders investigating at short time in-
tervals allow some high-frequency behaviors in the change market. Long-time traders
are interested in the general tendency and the volatility of the market along a mi-
croscopic greed. Short-time traders, however, are interested in fractional perceptions
and so in macroscopic greed. This leads to the wavelet analysis of financial time
series.
Recently, other models have been introduced in modeling financial time series
by means of fractals, which are in turn strongly related to wavelets. For example,
in Olsen & Associates, operating the largest financial database, has noticed that the
tick frequency has strongly increased in one decade, causing problems in studying the
time series extracted from such a database. Such problems can be due to transmission
delays, input errors and machine damages. So, some filtering procedure has to be done
Introduction 3
before using the data. The first point that one must take into account in filtering time
series is their scaling behavior. Scaling laws were empirically observed by Olsen et al
until 1990 (see [306]). A time series X(t) has a scaling law if its so-called partition
function has the form
N/t
t X
Sq (t) = |X(jt)|q ∼ tτq +1
N j=1
where N stands for the size of the series for some appropriate function τ . This estima-
tion is well understood when merging wavelet tools into fractal models. (See Chapter
8.) In such models, Fourier analysis could not produce good results in estimating such
behavior. Indeed, Fourier transform of time series is generally limited because a single
analysis window cannot detect features in the signals that are either much longer or
shorter than the window size. Moving-window Fourier transform (MWFT) slides a
fixed-size analysis window along the time axis and is able to detect non-stationarities.
The fixed-size window algorithm of MWFT limits the detection of cycles at wave-
lengths that are longer than the analysis windows, and non-stationarities in short
wavelengths (i.e., high frequencies) are smoothed. Use of the wavelet transform solves
this problem, because it uses narrow windows at high frequencies and wide windows
at low frequencies.
This book is devoted to developing the basic concepts of wavelet analysis neces-
sary for young researchers doing their Master’s level in science and researchers doing
doctoral studies in pure mathematical/physical sciences, as well as applied and inter-
acted ones by providing the basic tools required, with simple and rigorous methods.
It also aims to serve researchers at advanced levels by providing them the necessary
tools that will allow them to understand and adapt wavelet theory to their needs
such as supervision and development of research projects.
The book provides some highly flexible methods and ideas that can be manip-
ulated easily by undergraduate students, and thus may be of interest for Bachelors
in science by providing them a clear idea on what wavelets are, and thus permitting
them to decide in their scientific future.
Wavelets on Euclidean
Spaces
2.1 INTRODUCTION
Wavelets were discovered in the eighteenth century in petroleum exploration, and
since their discovery, they have proven to be powerful tools in many fields from
pure mathematics to physics to applied ones such as images, signals, medicine, fi-
nance and statistics. The study of their constructions and their properties, especially
in functional/signals decompositions on functional wavelet bases, has indeed grown
considerably.
In pure mathematics, wavelets constitute a refinement of Fourier analysis as they
compensate and/or resolve some anomalies in Fourier series. The first wavelet basis
has been, in fact, used before the pure mathematical discovery of wavelet bases since
the introduction of Haar system, which by the next has been proved indeed to be
a possible wavelet basis reminiscent of the regularity. Such a system dates back to
the beginning of the 20th century and was precisely discovered in the year 1909. It
was introduced in order to construct a functional basis permitting the representation
of all continuous functions by means of a uniformly convergent series. Recall that in
Fourier series, there are, as usual, many kinds of convergence that may be investigated
such as the point-wise convergence subject of the well-known Dirichlet theorem, the
uniform convergence which needs more assumptions on the series and the function,
the convergence in norm and Carleson’s or almost everywhere convergence. Each kind
of convergence requires special assumptions on the function. Although a Fourier series
converges in sense of Dirichlet, it does not imply that the graph of the partial Fourier
series converges to that of the function. This phenomenon is known as the Gibbs
phenomenon and is related to the presence of oscillations in the Fourier series near
the discontinuity points of the function. This means that the uniform convergence
is not sufficient. One of the challenging concepts in wavelet analysis is its ability to
describe well the behavior of the analyzed function near its singularities and join or
more precisely extend the notion of Littlewood-Palay decomposition.
To overcome some drawbacks of Fourier analysis, mathematicians have intro-
duced a bit of modification called the windowed Fourier transform by computing the
5
6 Wavelet Analysis: Basic Concepts and Applications
original Fourier transform of the analyzed signals on a special localized extra function
called the window. However, some situations remain non resolved especially with the
emergence of irregular signals or high-frequency variations. The major problem in the
windowed Fourier extension is due to the use of fixed window, which may not be well
adapted to other problems such as high fluctuations of non stationary signals. This
led researchers to think about a stronger tool taking into account nonlinear algo-
rithms, nonstationary signals, nonperiodical, volatile and/or fluctuated ones. It holds
that wavelets, since their discovery, have permitted to overcome these obstacles.
These powers are related simultaneously to many properties of wavelets. Indeed,
wavelet decomposition of functions joins Littlewood-Paley decomposition in many
cases. Wavelets provide simultaneous local analyses related to time–frequency. They
can be adapted to study-specific operators, especially differential and stochastic ones.
From the numerical and/or applied point of view, wavelets provide fast and accurate
algorithms, multi-resolution analyses as well as recursive schemes. These are very
important especially in big data analysis, image processing and also in numerical
resolution of partial differential equations.
In this chapter, we propose to review the basic concepts of wavelets as well as
their basic properties.
2.2 WAVELETS ON R
The first wavelet bases and thus analyses have been constructed on the real line R and
have been next extended to the general cases of the real/complex Euclidean spaces
Rm -Cm using different methods such as the natural tensor product.
Mathematically speaking, a wavelet or an analyzing wavelet on the Euclidean
space Rm may be defined in a large way as a function with specific properties that
may or may not be required necessarily as simultaneous assumptions. More precisely,
wavelet analysis is based primarily on the following points:
• An effective representation for standard functions,
• Robustness to the specification models,
• A reduction in the computation time,
• Simplicity of the analysis,
• An easy generalization and efficient, depending on the dimension,
• A location in time and frequency.
A wavelet is a function ψ ∈ L2 (R) that satisfies the following conditions:
• Admissibility,
dω
Z
|ψ̂(ω)|2 = Cψ < ∞. (2.1)
R+ |ω|
• Zero mean, Z +∞
ψ(0)
b = ψ(u)du = 0. (2.2)
−∞
Wavelets on Euclidean Spaces 7
By varying the parameters s and u, we may cover completely all the time–frequency
plane. This gives a full and redundant representation of the whole signal to be an-
alyzed (see [295]). This transform is called continuous because of the nature of the
parameters s and u that may operate at all levels and positions. The original signal
S can be reproduced knowing its CWT by the following relationship:
1 x−u dsdu
Z Z
S(t) = du,s (S)ψ . (2.7)
Cψ R s s2
It remains to notice that CWT is suitable for continuous time signals and those
representing varying singularities.
The function ψ ∈ L2 (R) satisfies some conditions such as the admissibility con-
dition and somehow describes Fourier-Plancherel identity and says that
dω
Z
|ψ̂(ω)|2 = Cψ < ∞. (2.8)
R+ |ω|
8 Wavelet Analysis: Basic Concepts and Applications
The function ψ has to also satisfy a number of vanishing moments, which is related
in wavelet theory to its regularity order. It states that
Z
p = 0, ..., m − 1, ψ(t)tp dt = 0. (2.9)
R
To analyze a signal by wavelets, one passes via the so-called wavelet transforms. A
wavelet transform is a representation of the signal by means of an integral form similar
to Fourier one in which the Fourier sine and/or cosine is replaced by the analyzing
wavelet ψ. In Fourier transform, the complex exponential source function yields the
copies eis. indexed by s ∈ R, which somehow represent frequencies. This transform
is continuous in the sense that it is indexed on the whole line of indices s ∈ R.
In wavelet theory, the situation is more unified. A CWT is also well known.
Firstly, a frequency, scale or a dilation or compression parameter s > 0 and a second
one related to time or position u ∈ R have to be fixed. The source function ψ
known as the analyzing wavelet is next transformed to yield some copies (replacing
the eis. ):
1 x−u
ψs,u (x) = √ ψ . (2.11)
s s
The CWT of a real-valued function f defined on the real line at the position u and
the scale s is defined by
Z ∞
ds,u (f ) = f (t)ψs,u (t)dt, ∀ u, s. (2.12)
−∞
By varying the parameters s and u, we cover completely all the time–frequency plane.
This gives a full and redundant representation of the whole signal to be analyzed (see
[295]). This transform is called continuous because of the nature of the parameters s
and u that may operate at all levels and positions.
So, wavelets operate according two parameters: the parameter u, which permits
translation of the graph of the source mother wavelet ψ, and the parameter s, which
permits compression or dilation of the graph of ψ. Computing or evaluating the
coefficients du,s means analyzing the function f with wavelets.
Theorem 2.1 The wavelet transform ds,u (f ) possesses some properties such as
and where
(τt f )(x) = f (x − t).
The proof of these properties is easy and the readers may refer to [163] for a review.
In wavelet theory, as in Fourier analysis theory, the original function f can be
reproduced via its CWT by an L2 -identity.
1 x − u dsdu
Z Z
f (x) = ds,u (f )ψ( ) 2 .
Cψ s s
dsdu
Z Z Z
ds,u (f )ds,u (g) = Cψ f (x)g(x)dx, ∀ f, g ∈ L2 (R).
s
Proof. We have
1 x−u 1
Z
{ fˆ(y)ψ̂(sy)e−iuy .
ds,u (f ) = f ∗ ψs (u) = √ f (x)ψ( )dx =
s s 2π
Consequently
1
Z Z
ds,u (f )ds,u (g)du = fb(y)ĝ(y)|ψ̂(sy)|2 dy.
u 2π y
Next, using Fubini’s rule, we observe that the last line is equal to
1 dadb
Z Z Z
sup F (x)G(x)dx − Ca,b (F )Ca,b (G)
||G||=1 Cψ 1/A≤a≤A |b|≤B a
1 dadb
Z
= sup Ca,b (F )Ca,b (G)
||G||=1 Cψ (a,b)∈[1/A,A]×[−B,B]
/ a
Now, Lemma 2.3 shows that the last quantity goes to 0 as A, B tends to +∞.
These are often called wavelet coefficients or detail coefficients of the signal S.
It holds that the set (ψj,k )j,k∈Z constitutes an orthonormal basis of L2 (R) and it
is called wavelet basis. A signal S with finite energy (in L2 (R)) is then decomposed
according to this basis into a series
∞ X
X
S(t) = dj,k ψj,k (t) (2.14)
j=0 k
called the wavelet series of S which replaces the reconstruction formula for the CWT.
Wavelets on Euclidean Spaces 11
For each j ∈ Z, Vj is called the approximation space at the scale or the level j. More
precisely, we have the following definition ([295]).
a) ∀j ∈ Z; V0 ⊂ V1 ⊂ .... ⊂ Vj ⊂ Vj+1 .
b) ∀j ∈ Z; f ∈ Vj ⇔ f (2.) ∈ Vj+1
[
e) Vj = L2 (R).
j∈R
f) ∀j ∈ Z; f ∈ Vj ⇐⇒ f (x − k) ∈ Vj
The property (a) reflects that the approximation of a signal at the resolution’s level
j + 1 contains the necessary information to yield the approximation at the level j.
Assertion (b) called dilatation’s property permits passing from a level of resolution
to another. The space Vj+1 contains signals that are coarser than Vj . Assertion (c)
means that a scaling function exists and permits the decomposition of the signal at the
starting level 0. The property (d) means that at very low resolution level (2−j −→ 0
as j −→ +∞) we lose all the details of the signal. At a minimal resolution, we
lose all the information about the signal. Assertion (e) implies that the signal may
be approximated with elements in Vj . At a maximal resolution, we reconstruct all
the whole signal. Finally, the property (f) of translation means that the space Vj is
invariant under integer translation.
Definition 2.5 The source function ϕ is called the scaling function of the MRA or
also the father wavelet.
It holds that this function generates all the subspaces Vj ’s of the MRA by acting
as dilation/translation parameters. Indeed, the property
(d) combined with
(f) in
j/2 j
Definition 2.4 implies that, for all j, ∈ Z, the set ϕj,k (x) = 2 ϕ(2 x − k) is an
k
orthogonal basis of Vj . From assertion (a), it holds that we may complete Vj in Vj+1
12 Wavelet Analysis: Basic Concepts and Applications
in the sense of the direct sum. Let Wj be the orthogonal supplementary of Vj in Vj+1 ,
that is
Vj+1 = Vj ⊕⊥ Wj . (2.16)
We will see that Wj plays a primordial role in representing the details of the analyzed
signal. This is why it is called the detail space at the level j. Iterating the relation
(2.16), we obtain an orthogonal decomposition for all J ∈ Z,
J,⊥
VJ = V0 ⊕⊥
M
Wj . (2.17)
j=0
The function ψ is called the wavelet function or the mother wavelet associated with
the scaling function ϕ of the MRA.
We will now explain the relationship between ϕ and ψ. It holds from (2.18) that any
S ∈ L2 (R) is decomposed into components according to the subspaces Wj ’s in the
sense that
+∞
X
S= S Wj ,
j=−∞
Observing now that the wavelet coefficients are already defined by (2.12) or (2.13) as
Z +∞
dj,k (S) = S(t)ψj,k (t)dt =< S, ψj,k >
−∞
Wavelets on Euclidean Spaces 13
we obtain the wavelet series of S defined already by (2.14). Returning next to (2.16)
or also (2.17), we may write that
−1 X
X +∞
XX
S= < S, ψj,k > ψj,k + < S, ψj,k > ψj,k . (2.19)
j=−∞ k j=0 k
It results from the nesting property of the Vj ’s and the fact that Wj ⊂ Vj+1
that the first part in (2.19) is an element of V0 which is generated by the basis
(ϕk (t) = ϕ(t − k))k . Consequently, we may also write that
−1 X
X X
< S, ψj,k > ψj,k = Ck ϕk .
j=−∞ k k
In fact, it may truncate at any level J and apply the ϕJ,k instead of ϕk to obtain
X ∞ X
X
S(x) = CJ,k ϕJ,k (x) + dj,k ψj,k (x) (2.21)
k j=J k
which is known as the wavelet series decomposition of S at the level J. In fact, the
coefficients that appears in this decomposition are evaluated via the relation
Z +∞
CJ,k = S(t)ϕJ,k (t)dt
−∞
and are called the scaling or approximation coefficients of S. The first component
in VJ reflects the global behavior or the tendency of the whole signal S and the
second component relative to the dj,k ’s represents the details of S and thus reflects
the dynamic behavior of the signal.
gk = (−1)k h1−k .
For more details, we refer to [177], [236], [295]. We will now prove that these relations
allow obtaining all the decomposition of a signal from each other through specific
algorithms. Indeed, consider a signal S and its approximation coefficients Cj,k and
details dj,k . For j ∈ Z, we have
Hence,
X
Cj,k = hl Cj+1,2k+l .
l
This means that the approximation at level j is obtained from the level j + 1 by the
intermediate of a filter. We have, in fact, the following definition.
Definition 2.9 The sequence H = (hk )k is called the discrete low-pass filter.
Hence, X
dj,k = hS, ψj,k i = gl Cj+1,l+2k .
l
This means that the approximation at level j is obtained from the level j + 1 by the
intermediate of a filter. We have here also the following definition.
Definition 2.10 The sequence G = (gk )k is called the discrete high-pass filter.
Figure 2.1 presents the decomposition algorithm due to Mallat [295]. It consists of a
cascade algorithm permitting to obtain all the levels of resolution.
We will now explain the inverse algorithm which permits to obtain the level j + 1
from the level j. Recall that the orthogonal projection of S on the approximation
space Vj is given by (Hardle et al (1997) [236])
X X
SVj = Cj,k ϕj,k + dj,k ψj,k
k k
Wavelets on Euclidean Spaces 15
Cj+1
H G
↓2 ↓2
Cj dj
H G
↓2 ↓2
Cj−1 dj−1
and that
or else X X
Cj+1,n = hn−2k Cj,m + gn−2k dj,m . (2.23)
k k
3. gn = (−1)n hl−n ; ∀ n.
X
4. hn gn+2j = 0, ∀ j ∈ N. (called mutual orthogonality).
n
16 Wavelet Analysis: Basic Concepts and Applications
Cj−1 dj−1
↑2 ↑2
H G
Cj dj
↑2 ↑2
H G
Cj+1
As a result,
1 c 2 1
Z 2π X
kFϕ k2L2 (R) = kFϕ kL2 (R) = | λk eikξ |2 Γϕ (ξ)dξ.
2π 2π 0 k
1 X
Now, observing that the map Φ : λ = (λk ) → Φ(λ) = √ λk eikξ is an iso-
2π
morphism from `2 (R) topL2 ([0, 2π]), to get the proposition, the map K defined on
L2 ([0, 2π]) by K(f ) = f Γϕ should also be an isomorphism. This is the case if and
1
only if Γϕ and are both bounded (in L∞ ([0, 2π])).
Γϕ
Sometimes we need more regularity on the wavelet basis especially in the case of par-
tial differential equations. Besides, we sometimes need compactly supported wavelets
and/or wavelets with fast decay which are more adaptable and suitable for boundary
and/or limit conditions due to partial differential equations.
It is sometimes necessary to assume that the scaling function ϕ and its derivatives
to an order r (usually the order of the MRA) are of fast decay. This permits to
avoid the assumption on the overlap function Γϕ of ϕ to be upper bounded, as this
characterization becomes already holding. Indeed, whenever ϕ is of fast decay, its
Fourier transform ϕb is therefore infinitely continuously differentiable and is also of
fast decay. It satisfies, in particular,
K
|ϕ(ξ)|
b ≤
1 + |ξ|
18 Wavelet Analysis: Basic Concepts and Applications
Proof 2.2 Assume that the trigonometric series m is analytic on the annulus
|Imω| < R. Observe next that m is also periodic with period 2π. As a consequence,
the function ω 7−→ m(ω + 2π) − m(ω) is vanishing on all the annulus. Now applying
the complex integral, we get
Z 2π Z −iR+2π
λk = m(ω)e−ikω dω = m(z)e−ikz dz
0 −iR
Consequently,
|λk | ≤ Ce−kR .
We will now show that in Definition 2.4 some properties and/or assumptions may
be deduced from others, which thus reduces the number of assumptions and gives
more flexibility in the introduction of multiresolution analyses. We have precisely the
following result.
Properties 2.14 Assume that assertions (a), (b) and (c) in Definition 2.4 hold and
denote Φ, the function defined by its Fourier transform
1
b + 2kπ)|2 )− 2 .
X
Φ(ω)
b = ϕ(ω)
b |ϕ(ω
k
It holds that Z +∞ X
ϕ(x)dx = 1 and ϕ(x − k) = 1.
−∞ k
where the coefficients hj,k are the analogues of the filter coefficients hk in the two-scale
relation. More precisely,
Z +∞ 1
hj,k = ϕ(x)2 2 ϕ(2j x − k)dx.
−∞
Z +∞
Denote next Iϕ = ϕ(x)dx. We may write
+∞
k
Z +∞ k j
− 2j
hj,k − Iϕ 2 ϕ( j ) = hj,k − ϕ( )2 2 ϕ(2j x − k)dx
2 +∞ 2j
Z +∞ k j
≤ C |x − |2 2 |ϕ(2j x − k)|dx
+∞ Z 2j
+∞
−3 2j
≤ C2 |x||ϕ(x)|dx
+∞
−3 2j
≤ C2 .
As a result,
j k j
hj,k = Iϕ 2− 2 ϕ( j
) + O(2−3 2 ).
2
Consequently,
k j
)ϕ(2j x − k) + O(2−3 2 ).
X
ϕ(x) = Iϕ ϕ( j
k
2
We now claim that Iϕ = 1. Indeed, let x0 be such that ϕ(x0 ) 6= 0. Let δ, η > 0 be
such that |ϕ(x) − ϕ(x0 )| ≤ whenever |x − x0 | ≤ η. Therefore,
X X
ϕ(x0 ) − Iϕ ϕ(x0 )ϕ(2j x − k) ≤ Iϕ |ϕ(2j x − k)|
k | kj −x0 |≤η
2
k
X
+ |ϕ( )| + |ϕ(x0 )| |ϕ(2j x − k)|
k
2j
| −x0 |>η
2j
j
+C + (2−3 2 ).
η k
Observe next whenever |x − x0 | ≤ then we necessarily get | j − x0 | > η, which
2 2
k η
implies that | j − x| ≥ . Hence,
2 2
C C
|ϕ(2j x − k)| ≤ ≤ .
(1 + |2j x − k|)4 (1 + |2j x − k|)2 (2j η2 )2
As a consequence,
X X 1 −3 2j
ϕ(x0 ) − Iϕ ϕ(x0 )ϕ(2j x − k) ≤ C η 2 + C + O(2 )
k k
(1 + |2j x 2 j
− k|) (2 2 )
j
≤ C2 −2j −2
η + C + O(2−3 2 ).
20 Wavelet Analysis: Basic Concepts and Applications
X 1
ϕ(y − k) − ≤ C.
k
Iϕ
X
Now observe that the function y 7−→ ϕ(y − k) is 1-periodic. Hence, for j large
k
enough, we get that such a function remains constant on the whole line R, which
means that
X 1
ϕ(y − k) = .
k
M
This yields
Z +∞ XZ k+1 Z 1 X 1
M= ϕ(y)dy = ϕ(y)dy = ( ϕ(y − k))dy = .
∞ k∈Z k 0 k∈Z
M
Hence, M 2 = 1 and \
thus may be chosen to be equal to 1.
We now prove that Vj = {0}. Let F be an element in the intersection. Hence, for
j
all j, we may write X
F = aj,k ϕj,k ,
k
where the aj,k are uniformly bounded. Moreover,
X
kF k2 = |aj,k |2 .
k
Consequently,
j
kχ[a,b] − Fj k22 ≤ C2− 2 .
Letting j → ∞ we observe Fj goes to χ[a,b] in L2 .
The adoption of these one-dimensional MRA and wavelets will be the starting point
to construct the multidimensional case. Indeed, let (Vj1 ) be a multiresolution analysis
of L2 (R) with a scaling function ϕ, a wavelet ψ and an orthogonal projections Pj1 .
Consider then orthogonal projection Pjd in L2 (Rd ) defined as the tensor product of d
copies of Pj1
Pjd = Pj1 ⊗ Pj1 ⊗ ... ⊗ Pj1 .
22 Wavelet Analysis: Basic Concepts and Applications
Denote Vjd = Pjd (L2 (Rd )). We have Vjd = Vj1 ⊗ Vj1 ⊗ ... ⊗ Vj1 . The closure in L2 (Rd )
of Vjd has an orthonormal basis
(
ϕdj,k = 2jd/2 ϕd (2j x − k) ; j ∈ R , k ∈ Rd
ϕd = ϕ ⊗ ϕ ⊗ ... ⊗ ϕ
Having in hand the explicit expression of Haar scaling function ϕ = χ[0,1[ , the latter
becomes √ X
χ[0,1[ (x) = 2 hk χ[k/2,(k+1)/2[ (x).
k
Next, because of the supports of the different functions appearing in the relation, we
obtain
hk = 0 , for all k 6= 0, 1
and thus, √ √
χ[0,1[ (x) = 2h0 χ[0,1/2[ (x) + 2h1 χ[1/2,1[ (x)
which yields immediately that
1
h0 = h1 = √ .
2
It results that the filter G is given by
1
g0 = −g1 = √ and gk = 0 , for all k 6= 0, 1.
2
The associated Haar mother wavelet is then expressed by
ψ = χ[0,1/2[ − χ[1/2,1[ .
We hereafter give the graphic illustrations of the Haar scaling function ϕ and the
Haar mother wavelet ψ (see Figure 2.3).
24 Wavelet Analysis: Basic Concepts and Applications
The subspace Vj is composed of square integrable functions that are of the form
j/2 j j j
aj,k 2 (1 − k + 2 x) , x ∈ [(k − 1)/2 , k/2 [,
j/2 j j j
f (x) = aj,k 2 (1 + k − 2 x) , x ∈ [k/2 , (k + 1)/2 [,
0 , others.
We now evaluate the filter coefficients. The two-scale relation (2.22) is written in the
present case as
√ X2
(1 − |x|)χ[−1,1[ (x) = 2 hk (1 − |2x − k|)χ[(k−1)/2,(k+1)/2[ (x).
k=−2
The coefficients hk are vanishing for |k| ≥ 3. Next by choosing for x the values −1,
−1 1
, 0, and 1, we obtain
2 2
1 1
h−2 = h2 = 0 , h−1 = h1 = √ and h0 = √ .
2 2 2
The G filter will be
1 −1
g0 = g2 = √ and g1 = √ .
2 2 2
The Schauder mother wavelet will be expressed by
1 1
ψ(x) = (1 − |2x|)χ[−1/2,1/2[ (x) − (1 − |2x − 1|)χ[0,1[ (x) + (1 − |2x − 2|)χ[1/2,3/2[ (x).
2 2
The following graph (Figure 2.4) represents Schauder ϕ and ψ.
∀ x, Φ(x) = Σ hk Φ(2x − k)
26 Wavelet Analysis: Basic Concepts and Applications
and √ X
∀ x, ψ(x) = 2 g(n)φ(2x − k)
where X
h(n − 2k)h(n − 2l) = δkl ,
n
X √
h(n) = 2,
n
and
g(n) = (−1)n h(1 − n)
If the filter H = (h(n))n∈Z is finite, the support of the scaling function is also finite.
Let
supp(φ) ⊂ [Nmin , Nmax ].
Hence,
For example for N = 2, we obtain the so-called Db4 wavelet, whose low-pass filter
coefficients hk are √ √
1+ 3 3+ 3
h0 = √ , h1 = √
4 2 4 2
√ √
3− 3 1− 3
h2 = √ , h3 = √
4 2 4 2
The next graph illustrates some wavelet function psi due to Daubechies. The repre-
sentation illustrates clearly the strong relation between the regularity, the number of
vanishing moments and the support length N .
Wavelets on Euclidean Spaces 27
We will start by the most simple example of wavelets and multiresolution analysis
due to the Haar case and show that, even though this is a nonregular case, it may serve
to construct higher order multiresolutions. These are the so-called spline wavelets.
So, let ϕ = χ[0,1] be the Haar scaling function, which is equivalently the char-
acteristic function of the interval [0, 1]. Let r > 0 be in N. The spline wavelets are
defined inductively by
ϕ0 = ϕ and ϕr = ϕr−1 ∗ ϕ, ∀r ≥ 1.
ψ(X) = (2 − XAX)e−XAX/2
with m > 0. In one dimension, the positive half-line is a convex cone. Thus, a natural
generalization to two dimensions will be a wavelet whose support in spatial frequency
space is contained in a convex cone with apex at the origin. Let C ≡ C(α, β) be the
convex cone determined by the unit vectors eα , eβ , where α < β, β − α < π and for
all θ, eθ ≡ (cos θ, sin θ). The axis of the cone is ξαβ = e α+β . In other words,
2
C(α, β) = k ∈ R2 , α ≤ arg(k) ≤ β
= k ∈ R2 , k.ξαβ ≥ eα .ξαβ = eβ .> 0 .
C(α,
e β) = C(α
b , β).
b
where
π b π
b=β−
α , β =α+
2 2
and
eα .eαb = eβ .eβb = 0.
where η ∈ Ce and l, m ∈ N∗ . Note that such a wavelet is also supported by the cone
C. It satisfies the admissibility condition
d2 k
Z
C,η
AψC,η ≡ (2π)2 |ψblm (k)|2 < ∞.
lm |k|2
The following result of Antoine et al is proved in [13] and yields an explicit form for
the two-dimensional Cauchy wavelet.
Wavelets on Euclidean Spaces 31
2.8 EXERCISES
Exercise 1.
√
Denote φ = χ[0,1[ and φjk (t) = 2j φ(2j t − k), t ∈ R. Denote further E = L2 ([0, 1[)
the vector space of square-integrable functions on [0, 1[. Next for j fix V j , the vector
k k+1
space of constant functions, on { j ; j } for k ∈ {0; 2j − 1}.
2 2
a. Show that (φjk )k∈{0,2j −1} is an orthonormal basis of V j . < (φjk ); (φjh ) >= 0∀(k, h) ∈
{0, 2j − 1} × {0, 2j − 1}. and that the elements φjk are unitary.
Exercise 2.
Consider the same assumptions as in Exercise 1 above.
a. Find two functions (ψk1 )k∈{0,1} constituting an orthonormal basis of a vector space
W 1 such that V 2 = V 1 ⊕ W 1 .
b. Express the functions (ψk1 )k∈{0,1 } and (φ1k )k∈{0,1 } by means of (φ2k )k∈{0,3 }.
c. Express conversely the functions (φ2k )k∈{0,3 } in terms of (φ1k )k∈{0,1 } and
(ψk1 )k∈{0,1 }.
Exercise 3.
Consider the same assumptions as in Exercises 1 and 2 above.
a. Provide the 2j Haar wavelets (ψk1 )k∈{0,1 } which constitute an orthonormal basis
of W j = V j+1 V j .
d. Write the approximation coefficients ajk and the detail coefficients (djk ) by means
of the (aj+1
k ).
32 Wavelet Analysis: Basic Concepts and Applications
Exercise 4.
Consider the discrete signal S = [2 4 8 12 14 0 2 1] and the Haar multiresolution
analysis on R.
a. Decomposer S in V 0 ⊕2j=0 W j
b. Sketch the graphs of S relatively to V 2 , V 1 and V 0 .
b. Sketch the graphs of S relatively to the details W 2 , W l and W 0
Exercise 5.
Let ψ : R → C be wavelet in the Schwarz class such that ψ(0)
b = 0). Denote
Z
C= |xψ(x)|dx.
R
Exercise 6.
Let f, ψ : R → C be in the Schwarz class such that ψ(0) = 0) and denote
Z +∞ Z
dt
F (t) = ψs (t − u)Wf (u, s)du .
0 R s2
Show that F = Cψ f , where Cψ is a constant depending on ψ and not on f .
Exercise 7.
Let ψ be a wavelet such that ψ(0)
b = 0, where
Z
ψ(w)
b = ψ(t)e−izt dt, ∀w ∈ R.
R
Wavelets on Euclidean Spaces 33
Exercise 8.
Let f ∈ L2 (R) be such that fb(ξ) = 0 for ξ < 0. For z ∈ C such that I, (z) ≥ 0, let
1
Z +∞
fp (z) = (iξ)p fb(ξ)eizξ dξ.
π 0
Exercise 9.
Let for N ∈ N, f : [0, N ] → R be continuous and (Vm )m≥0 be a multiresolution
analysis on [0, N ] with a scaling function ϕ.
1) Verify that for all ε > 0 there exists m ≥ 0 such that for all elements t1 , t2 ∈ [0, N ]
with |t1 − t2 | < 2−m we have |f (t1 ) − f (t2 )| < ε.
2) Denote next tm,n+1/2 = (n + 21 )2−m and
N −1
2X
fε (t) = f (tm,n+1/2 )ϕm,n (t).
n=0
Prove that
|f (t) − fε (t)| < ε, ∀t ∈ [n2−m , (n + 1)2−m ].
3) Show that √
kf − fε k2 ≤ ε N .
Exercise 10.
Part A. Consider the Haar system on R with its scaling function ϕ = χ[0,1[ and
wavelet function ψ(x) = χ[0,1/2[ − χ[1/2,1[ . Let next f (t) = (t2 + 1)χ[−1,1[ and its
projections on Vj and Wj be, respectively,
X X
fj = aj,k ϕj,k and dj = dj,k ψj,k .
k k
a) Evaluate f0 (x, y), f1 (x, y), f2 (x, y), d0 (x, y) and d1 (x, y).
b) Draw the graphs of f (x, y),f0 (x, y), f1 (x, y), f2 (x, y), d0 (x, y) and d1 (x, y).
CHAPTER 3
Wavelets extended
35
36 Wavelet Analysis: Basic Concepts and Applications
In our case, we will focus on the so-called linear or affine group on Rn . These
are groups that can be realized as groups of linear transformations of a vector space
or, equivalently, as matrix groups. We will see that group representation theory may
help explain the essential properties of wavelet transforms.
Consider the connected affine group G+ , also called the ax + b group, consisting
of transformations of R of the type x 7−→ gx = ax + b, x ∈ R, where a > 0, b ∈ R
and where we denoted g = (b, a) ∈ G+ with the multiplication law
gl g2 = (b1 + a1 b2 , a1 a2 ).
dadb
It holds that an invariant measure on G is dµ (a, b) = .
a2
A nontrivial unitary irreducible representations of G on Hilbert space L2 (R, dx)
will be
1 x−b
(U (a, b)f )(x) = √ f ( ), f ∈ L2 (R, dx), (a, b) ∈ G.
a a
For the sequel, we define the following operators on the set of functions defined on
the Euclidean space R:
These operators are easily seen to be unitary. Using these operators, we immediately
deduce that
(U (a, b)f ) = Da−1 Tb f.
As a consequence, the wavelet transform of f at the scale a > 0 and the position
b ∈ R may be expressed by means of the representation U as
The results presented in this section may be extended in a usual way to multidimen-
sional case.
Wavelets extended 37
• ϕj,k K = N, · · · , 2j − N,
where the functions ϕ̃0k and ϕ̃1k are defined as linear combinations of ϕj,k with supports
containing 0 or 1 and with suitable polynomial coefficients such as
N
X −1
• ϕ̃0k (x) = Pk0 (l)ϕ(x − l),
−N +1
N
X −1
• ϕ̃1k (x) = Pk1 (l)ϕ(x − l).
−N +1
where ϕ is the scaling function of an AMR in L2 (R) with compact support. Let Vj∗
be the subspace of Vjloc of functions with coefficients having polynomial behavior on
some neighborhood of 0 or 1. The degrees of polynomials depend possibly on j. More
precisely, they have proceeded as follows. For j ≥ 0, let Nj = 2j − 2N + 2M + 3, Mj =
min(M, Nj − 1) and
( )
fkj ϕj,k ; ∃Pl , Pr PMj , fkj 1, fkj
X
Vj∗ = ∈ = Pl (k), k ≤ N − = Pr (k), k ≥ 2j − N + 1 .
k∈Z
Random documents with unrelated
content Scribd suggests to you:
Thalásio, a Seawater-colour or blew like the Sea-waues.
Thalassagliéne, an hearbe very salt in taste.
Thalasséglie, as Potamánti.
Thalassiárca, an Admirall.
Thalassíno, as Thalásio.
Thalassinténo, a kind of salt wine vsed among the Grecians.
Thalassoméli, a kind of surrop made of salt sea-water.
Thália, one of the nine Muses, as much to say capacity.
Thaliétro, the hearbe Laske-wort or Flax-wort.
Thállo, vsed of the Florentines for the colour Greene. Vsed also for a
bough that still continueth greene.
Thalmudísta, a professor of the Chaldeans or Egiptians old lawes.
Thalofóri, bringers of branches or boughes.
Thápsia, an hearbe like Dill, which being cut yeeldes a iuice as white as
milke.
Thargélie fẻste, certaine feastes and games kept in honour of Apóllo.
Thásia, a kind of wine in Egypt. Also as Thápsia.
Tháspi, a kind of Mustard-seed.
Tháu, an Hebrew letter vsed often misteriously for the crosse of Christ, as
also for the number of 400.
Theamẻde, a stone which contrary to the Adamant putteth away iron from
it.
Theangélide, a kind of hearb.
Theatineríe, hipocriticall fryers trikes.
Theatíni, a kind of order of friars, that much dissemble and faine purity.
Theatrále, of or belonging to a Theater.
Theátrico, a plaier in a Theater, a Stage plaier.
Theátro, a Theater, an open place where all men may sit, see and behold,
or wherein one sheweth, plaieth or declareth. Also a stage or skaffold.
Also any spectacle or common-play.
Thebáide, a kind of poore drie Date.
Theicóne, as Trogéte.
Thelétteri, Female or shee Ferne.
Thélico, a stone like an Oliue-stone, which being beaten to pouder
breaketh the stone in the kidneis.
Thelifóno, the hearbe Libard-bane, the roote whereof doth kill Serpents.
Theligónio, the hearb Culerage, Arse-mart or as some call it the grace of
God.
Thelipteríno, a kind of Fern.
Théma, a theame, an argument, a subiect, a ground or matter supposed
to be written or spoken of. Also a figure set by Astronomers containing
the state of heauen and the planets at any certaine time. Also an
horoscope or consideration of the time of ones birth or natiuity vnder
this or that planet.
Themísia, Oke of Ierusalem or Mille-graine.
Themología, as Ethimología.
Thenáre, the fleshy part of the hand between the thumbe and the middle
finger.
Thendálida, as Sanguinária.
Theodía, Deity, God-head, the word of God, Diuinity. Also Dauids Psalmes.
Theogónia, the Originall, the beginning or generation of the Gods.
Theolẻpto, inspired of God.
Theologáre, to reason of holy things or preach diuinity.
Theología, reasoning of holy things, science of diuinity, profession of holy
things.
Thómbo, a thump.
Thopalóne, a water-rose.
Thoráce, as Toráce.
Thoriphétro, the hearbe Lions-paw.
Thóro, as Tóro, a bull.
Thrácia, a kinde of blacke stone or Agath, which being steeped in water,
burneth, and sprinkled with Oyle, it quencheth.
Thrásia, a wind between the North and the summer-setting of the Sunne.
Thréni di Gieremía, the lamentations of Ieremy.
Threustóne, a kind of sweet gum like Frankincense.
Thriálli, a kind of Primrose or Couslip.
Thridáce, Sper-hawke-hearb or Hawk-weede.
Thridacẻa, the female Mandrágora.
Thríssa, a kind of Purple fish.
Thróno, a throne or seate of maiesty, an Emperours or Kings seate. Vsed
also for a Kingdome.
Thù, Ground-pine or field-cipres. Vsed also for incense or Frankincense.
Thúmbo, a thump.
Thúnnide, as Thínnide.
Thúnno, as Thínno.
Thuríbolo, as Torríbolo.
Tì, to thee, thee. Also from thee. Looke for it in my rules.
Tiáde, as Ptiáde.
Tiáme, any Ouen or Furnace.
Tiáno, any Ouen or Furnace.
Tiára, a turbant or round wreath of linnen as Turks weare on their heads,
but it was properly a round ornament or head-attire that Princes were
wont to weare in Persia.
Tíbia, any kind of flute, pipe or recorder. Also a trumpet. Also the shanke
or shin-bone of a mans legge which, comprehendeth all the space
betweene the knee and the foote, and containeth two bones, of which
the greater is called Tíbia, it hath also beene vsed for a thrashing of
Corne.
Tibiále, a shin-piece, a greaue or leg-harnish. Also a nether-stocke or
knee-hose.
Tibiáre, to sound vpon any flute, a pipe or recorder. Also to thrash corne.
Tibiatóre, a piper vpon pipes, flutes or recorders, a trumpeter. Also a
thrasher.
Tibidrágo, a kinde of strong Spanish wine or sacke, we call it Rubiedauy.
Tibilústri, daies in which men were wont to go with pipes and trumpets, as
it were in procession about their lands.
Tíbis, a birde like a Storke feeding on Serpents egges.
Tibribástio, iargon, gibrish, pedlers french, fusti[an] tongue, rogues
language.
Tíbulo, a kind of Pine-tree.
Tiburíno, a monstruous Sea-fish which hath two rankes of sharpe teeth.
ebookbell.com