Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis - - PowerPoint PPT Presentation

curvelets contourlets shearlets lets etc multiscale
SMART_READER_LITE
LIVE PREVIEW

Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis - - PowerPoint PPT Presentation

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more) L. Jacques, L. Duval, C. Chaux,


slide-1
SLIDE 1

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example

Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

  • L. Jacques, L. Duval, C. Chaux, G. Peyré

UCL, IFPEN, AMU, Dauphine

2015

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-2
SLIDE 2

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example

In just one slide

Figure : A standard, “dyadic”, separable wavelet decomposition

Where do we go from here? 15 years, 300+ refs in 90 minutes

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-3
SLIDE 3

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example

Wavelets for the eye

Artlets: painting wavelets (Hokusai/A. Unser)

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-4
SLIDE 4

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example

Wavelets for 1D signals

1D scaling functions and wavelets

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-5
SLIDE 5

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example

Wavelets for 2D images

2D scaling functions and wavelets

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-6
SLIDE 6

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example

1D signals

1D and 2D data appear quite different, even under simple:

◮ time shift ◮ scale change ◮ amplitude drift

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-7
SLIDE 7

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example

1D signals

Figure : 1D and 2D → 1D related signals

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-8
SLIDE 8

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example

2D images

Figure : 1D → 2D and 2D related images

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-9
SLIDE 9

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example

1D signals & 2D images

Only time shift/scale change/amplitude drift between:

◮ John F. Kennedy Moon Speech (Rice Stadium, 12/09/1962) ◮ A Man on the Moon: Buzz Aldrin (Apollo 11, 21/07/1969)

Partial conclusion:

◮ There IS a curse of dimensionality. . . between 1 and 2

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-10
SLIDE 10

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 10/34

1.5D signals: motivations for 2D directional "wavelets"

Figure : Geophysics: seismic data recording (surface and body waves)

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-11
SLIDE 11

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 10/34

1.5D signals: motivations for 2D directional "wavelets"

Offset (traces) Time (smpl)

50 100 150 200 250 300 100 200 300 400 500 600 700

Figure : Geophysics: surface wave removal (before)

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-12
SLIDE 12

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 10/34

1.5D signals: motivations for 2D directional "wavelets"

Offset (traces) Time (smpl)

50 100 150 200 250 300 100 200 300 400 500 600 700

(b)

Figure : Geophysics: surface wave removal (after)

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-13
SLIDE 13

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 10/34

1.5D signals: motivations for 2D directional "wavelets"

Issues in geophysics:

◮ different types of waves on seismic "images"

◮ appear hyperbolic [layers], linear [noise] (and parabolic)

◮ not the standard “mid-amplitude random noise problem” ◮ no contours enclosing textures, more the converse ◮ kind of halfway between signals and images (1.5D) ◮ yet local, directional, frequency-limited, scale-dependent

structures to separate

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-14
SLIDE 14

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 11/34

Agenda

◮ To survey 15 years of improvements in 2D wavelets

◮ spatial, directional, frequency selectivity increased ◮ sparser representations of contours and textures ◮ from fixed to adaptive, from low to high redundancy ◮ generally fast, practical, compact (or sparse?), informative ◮ 1D/2D, discrete/continuous hybridization

◮ Outline

◮ introduction + early days ( 1998) ◮ fixed: oriented & geometrical (selected): ◮ ± separable (Hilbert/dual-tree wavelet) ◮ isotropic non-separable (Morlet-Gabor) ◮ anisotropic scaling (ridgelet, curvelet, contourlet, shearlet) ◮ (hidden bonuses): ◮ adaptive, lifting, meshes, spheres, manifolds, graphs ◮ conclusions

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-15
SLIDE 15

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 12/34

Images are pixels (but...):

Figure : Image block as a (canonical) linear combination of pixels

◮ suffices for basic counting, enhancement, filtering ◮ very limited in higher level understanding tasks

◮ looking for other (meaningful) linear combinations ◮ what about

41 + 37 + 40 + 36, 41 + 37 − 40 − 36 41 − 37 + 40 − 36, 41 − 37 − 40 + 36?

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-16
SLIDE 16

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 12/34

Images are pixels (but...):

Figure : Image block as a (canonical) linear combination of pixels

◮ suffices for basic counting, enhancement, filtering ◮ very limited in higher level understanding tasks

◮ looking for other (meaningful) linear combinations ◮ what about

41 + 37 + 40 + 36, 41 + 37 − 40 − 36 41 − 37 + 40 − 36, 41 − 37 − 40 + 36?

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-17
SLIDE 17

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 12/34

Images are pixels (but...):

Figure : Image block as a (canonical) linear combination of pixels

◮ suffices for basic counting, enhancement, filtering ◮ very limited in higher level understanding tasks

◮ looking for other (meaningful) linear combinations ◮ what about

41 + 37 + 40 + 36, 41 + 37 − 40 − 36 41 − 37 + 40 − 36, 41 − 37 − 40 + 36?

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-18
SLIDE 18

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 12/34

Images are pixels (but...):

Figure : Image block as a (canonical) linear combination of pixels

◮ suffices for basic counting, enhancement, filtering ◮ very limited in higher level understanding tasks

◮ looking for other (meaningful) linear combinations ◮ what about

41 + 37 + 40 + 36, 41 + 37 − 40 − 36 41 − 37 + 40 − 36, 41 − 37 − 40 + 36?

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-19
SLIDE 19

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 12/34

Images are pixels (but...):

Figure : Image block as a (canonical) linear combination of pixels

◮ suffices for basic counting, enhancement, filtering ◮ very limited in higher level understanding tasks

◮ looking for other (meaningful) linear combinations ◮ what about

41 + 37 + 40 + 36, 41 + 37 − 40 − 36 41 − 37 + 40 − 36, 41 − 37 − 40 + 36?

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-20
SLIDE 20

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 12/34

Images are pixels (but...):

Figure : Image block as a (canonical) linear combination of pixels

◮ suffices for basic counting, enhancement, filtering ◮ very limited in higher level understanding tasks

◮ looking for other (meaningful) linear combinations ◮ what about

41 + 37 + 40 + 36, 41 + 37 − 40 − 36 41 − 37 + 40 − 36, 41 − 37 − 40 + 36?

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-21
SLIDE 21

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 13/34

Images are pixels (but...):

A review in an active research field:

◮ (partly) inspired by:

◮ early vision observations [Marr et al.] ◮ sparse coding: wavelet-like oriented filters and receptive fields

  • f simple cells (visual cortex) [Olshausen et al.]

◮ a widespread belief in sparsity

◮ motivated by first successes (JPEG 2000 compression) ◮ aimed either at pragmatic or heuristic purposes:

◮ known formation model or unknown information content

◮ developed through a legion of *-lets (and relatives)

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-22
SLIDE 22

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 13/34

Images are pixels, wavelets are legion

Room(let) for improvement:

Activelet, AMlet, Armlet, Bandlet, Barlet, Bathlet, Beamlet, Binlet, Bumplet, Brushlet, Caplet, Camplet, Chirplet, Chordlet, Circlet, Coiflet, Contourlet, Cooklet, Craplet, Cubelet, CURElet, Curvelet, Daublet, Directionlet, Dreamlet, Edgelet, FAMlet, FLaglet, Flatlet, Fourierlet, Framelet, Fresnelet, Gaborlet, GAMlet, Gausslet, Graphlet, Grouplet, Haarlet, Haardlet, Heatlet, Hutlet, Hyperbolet, Icalet (Icalette), Interpolet, Loglet, Marrlet, MIMOlet, Monowavelet, Morelet, Morphlet, Multiselectivelet, Multiwavelet, Needlet, Noiselet, Ondelette, Ondulette, Prewavelet, Phaselet, Planelet, Platelet, Purelet, QVlet, Radonlet, RAMlet, Randlet, Ranklet, Ridgelet, Riezlet, Ripplet (original, type-I and II), Scalet, S2let, Seamlet, Seislet, Shadelet, Shapelet, Shearlet, Sinclet, Singlet, Slantlet, Smoothlet, Snakelet, SOHOlet, Sparselet, Spikelet, Splinelet, Starlet, Steerlet, Stockeslet, SURE-let (SURElet), Surfacelet, Surflet, Symmlet, S2let, Tetrolet, Treelet, Vaguelette, Wavelet-Vaguelette, Wavelet, Warblet, Warplet, Wedgelet, Xlet, not mentioning all those not in -let!

Now, some reasons behind this quantity

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-23
SLIDE 23

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 13/34

Images are pixels, but altogether different

Figure : Different kinds of images

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-24
SLIDE 24

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 13/34

Images are pixels, but altogether different

Figure : Different kinds of images

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-25
SLIDE 25

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 13/34

Images are pixels, but might be described by models

“Template” image decomposition models:

◮ edge cartoon + texture [Meyer-2001]:

inf

u E(u) =

|∇u| + λv∗, f = u + v

◮ edge cartoon + texture + noise [Aujol-Chambolle-2005]:

inf

u,v,w F(u, v, w) = J(u)+J∗

v µ

  • +B∗ w

λ

  • + 1

2αf −u−v −wL2

◮ heuristically: piecewise-smooth + contours + geometrical

textures + noise (or unmodeled)

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-26
SLIDE 26

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 13/34

Images are pixels, but resolution/scale helps with models

Coarse-to-fine and fine-to-coarse relationships

Figure : Notion of sufficient resolution [Chabat et al., 2004]

◮ discrete 80’s wavelets were “not bad” for: piecewise-smooth

(moments) + contours (gradient-behavior) + geometrical textures (oscillations) + noise (orthogonality)

◮ yet, not enough with noise, complicated images (poor sparsity

decay)

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-27
SLIDE 27

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 13/34

Images are pixels, but decay with regularity

Figure : Piecewise regular image

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-28
SLIDE 28

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 13/34

Images are pixels, but decay with regularity

Compressibility vs regularity: MSE with M-term approximation

◮ 1D

◮ piecewise C α → O(M−2α)

◮ 2D

◮ C α → O(M−α) (standard wavelets) ◮ piecewise C α/C α → O(M−1) (standard wavelets) ◮ piecewise C 2/C 2 → O(M−2) (triangulations)

◮ Notes:

◮ very imprecise statements, many deeper results ◮ piecewise C 2/C 2 → O(M−2f (M)) w/ directional wavelets? ◮ do much better with other regularities (α = 2, BV)?

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-29
SLIDE 29

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 14/34

Images are pixels, but sometimes deceiving

Figure : Real world image and illusions

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-30
SLIDE 30

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 14/34

Images are pixels, but sometimes deceiving

Figure : Real world image and illusions

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-31
SLIDE 31

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 14/34

Images are pixels, but sometimes deceiving

Figure : Real world image and illusions

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-32
SLIDE 32

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 15/34

Images are pixels, but resolution/scale helps

To catch important "objects" in their context

◮ use scales, pyramidal or multiresolution schemes, ◮ combine w/ different description/detection/modeling:

◮ smooth curve or polynomial fit, oriented regularized derivatives

(Sobel, structure tensor), discrete (lines) geometry, parametric curve detectors (e.g. Hough transform), mathematical morphology, empirical mode decomposition, local frequency estimators, Hilbert and Riesz (analytic and monogenic), quaternions, Clifford algebras, optical flow, smoothed random models, generalized Gaussian mixtures, warping operators, etc.

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-33
SLIDE 33

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 16/34

Images are pixels, and need efficient descriptions

Depend on application, with sparsity priors:

◮ compression, denoising, enhancement, inpainting, restoration,

contour detection, texture analysis, fusion, super-resolution, registration, segmentation, reconstruction, source separation, image decomposition, MDC, learning, etc.

100 200 300 400 500 600 700 800 900 1000 10 −1 10 10 1 10 2 10 3 10 4 Magnitude Index

Figure : Image (contours/textures) and decaying singular values

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-34
SLIDE 34

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 17/34

Images are pixels: a guiding thread (GT)

Figure : Memorial plaque in honor of A. Haar and F. Riesz: A szegedi matematikai iskola világhírű megalapítói, courtesy Prof. K. Szatmáry

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-35
SLIDE 35

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 18/34

Guiding thread (GT): early days

Fourier approach: critical, orthogonal

Figure : GT luminance component amplitude spectrum (log-scale)

Fast, compact, practical but not quite informative (not local)

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-36
SLIDE 36

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 18/34

Guiding thread (GT): early days

Scale-space approach: (highly)-redundant, more local

Figure : GT with Gaussian scale-space decomposition

Gaussian filters and heat diffusion interpretation Varying persistence of features across scales ⇒ redundancy

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-37
SLIDE 37

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 18/34

Guiding thread (GT): early days

Pyramid-like approach: (less)-redundant, more local

Figure : GT with Gaussian pyramid decomposition

Varying persistence of features across scales + subsampling

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-38
SLIDE 38

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 18/34

Guiding thread (GT): early days

Differences in scale-space with subsampling

Figure : GT with Laplacian pyramid decomposition

Laplacian pyramid: complete, reduced redundancy, enhances image singularities, low-activity regions/small coefficients, algorithmic

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-39
SLIDE 39

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 18/34

Guiding thread (GT): early days

Isotropic wavelets (more axiomatic) Consider Wavelet ψ ∈ L2(R2) such that ψ(x) = ψrad(x), with x = (x1, x2), for some radial function ψrad : R+ → R (with adm. conditions). Decomposition and reconstruction For ψ(b,a)(x) = 1

aψ( x−b a ), Wf (b, a) = ψ(b,a), f with reconstruc-

tion: f (x) =

2π cψ

+∞

  • R2 Wf (b, a) ψ(b,a)(x) d2b da

a3

if cψ = (2π)2

R2 | ˆ

ψ(k)|2/k2 d2k < ∞.

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-40
SLIDE 40

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 18/34

Guiding thread (GT): early days

Wavelets as multiscale edge detectors: many more potential wavelet shapes (difference of Gaussians, Cauchy, etc.)

Figure : Example: Marr wavelet as a singularity detector

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-41
SLIDE 41

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 18/34

Guiding thread (GT): early days

Definition The family B is a frame if there exist two constants 0 < µ♭ µ♯ < ∞ such that for all f µ♭f 2

  • m

|ψm, f |2 µ♯f 2 Possibility of discrete orthogonal bases with O(N) speed. In 2D: Definition Separable orthogonal wavelets: dyadic scalings and translations ψm(x) = 2−jψk(2−jx − n) of three tensor-product 2-D wavelets ψV (x) = ψ(x1)ϕ(x2), ψH(x) = ϕ(x1)ψ(x2), ψD(x) = ψ(x1)ψ(x2)

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-42
SLIDE 42

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 18/34

Guiding thread (GT): early days

1D scaling functions ψ(x1) and wavelets ϕ(x2)

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-43
SLIDE 43

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 18/34

Guiding thread (GT): early days

So, back to orthogonality with the discrete wavelet transform: fast, compact and informative, but... is it sufficient (singularities, noise, shifts, rotations)?

Figure : Discrete wavelet transform of GT

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-44
SLIDE 44

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 19/34

Oriented, ± separable

To tackle orthogonal DWT limitations

◮ 1D, orthogonality, realness, symmetry, finite support (Haar)

Approaches used for simple designs (& more involved as well)

◮ relaxing properties: IIR, biorthogonal, complex ◮ M-adic MRAs with M integer > 2 or M = p/q ◮ hyperbolic, alternative tilings, less isotropic decompositions ◮ with pyramidal-scheme: steerable Marr-like pyramids ◮ relaxing critical sampling with oversampled filter banks ◮ complexity: (fractional/directional) Hilbert, Riesz, phaselets,

monogenic, hypercomplex, quaternions, Clifford algebras

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-45
SLIDE 45

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 20/34

Oriented, ± separable

Illustration of a combination of Hilbert pairs and M-band MRA

  • H{f }(ω) = −ı sign(ω)

f (ω)

−4 −3 −2 −1 1 2 3 −0.8 −0.6 −0.4 −0.2 0.2 0.4 0.6 0.8 1

Figure : Hilbert pair 1

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-46
SLIDE 46

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 20/34

Oriented, ± separable

Illustration of a combination of Hilbert pairs and M-band MRA

  • H{f }(ω) = −ı sign(ω)

f (ω)

−4 −3 −2 −1 1 2 3 −0.5 0.5 1

Figure : Hilbert pair 2

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-47
SLIDE 47

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 20/34

Oriented, ± separable

Illustration of a combination of Hilbert pairs and M-band MRA

  • H{f }(ω) = −ı sign(ω)

f (ω)

−4 −3 −2 −1 1 2 3 4 −2 −1.5 −1 −0.5 0.5 1 1.5 2

Figure : Hilbert pair 3

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-48
SLIDE 48

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 20/34

Oriented, ± separable

Illustration of a combination of Hilbert pairs and M-band MRA

  • H{f }(ω) = −ı sign(ω)

f (ω)

−4 −3 −2 −1 1 2 3 −2 −1 1 2 3

Figure : Hilbert pair 4

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-49
SLIDE 49

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 20/34

Oriented, ± separable

Illustration of a combination of Hilbert pairs and M-band MRA

  • H{f }(ω) = −ı sign(ω)

f (ω) Compute two wavelet trees in parallel, wavelets forming Hilbert pairs, and combine, either with standard 2-band or 4-band

Figure : Dual-tree wavelet atoms and frequency partinioning

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-50
SLIDE 50

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 21/34

Oriented, ± separable

Figure : GT for horizontal subband(s): dyadic, 2-band and 4-band DTT

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-51
SLIDE 51

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 21/34

Oriented, ± separable

Figure : GT (reminder)

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-52
SLIDE 52

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 21/34

Oriented, ± separable

Figure : GT for horizontal subband(s) (reminder)

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-53
SLIDE 53

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 21/34

Oriented, ± separable

Figure : GT for horizontal subband(s): 2-band, real-valued wavelet

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-54
SLIDE 54

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 21/34

Oriented, ± separable

Figure : GT for horizontal subband(s): 2-band dual-tree wavelet

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-55
SLIDE 55

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 21/34

Oriented, ± separable

Figure : GT for horizontal subband(s): 4-band dual-tree wavelet

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-56
SLIDE 56

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 22/34

Directional, non-separable

Non-separable decomposition schemes, directly n-D

◮ non-diagonal subsampling operators & windows ◮ non-rectangular lattices (quincunx, skewed) ◮ non-MRA directional filter banks ◮ steerable pyramids ◮ M-band non-redundant directional discrete wavelets ◮ served as building blocks for:

◮ contourlets, surfacelets ◮ first generation curvelets with (pseudo-)polar FFT, loglets,

directionlets, digital ridgelets, tetrolets

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-57
SLIDE 57

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 22/34

Directional, non-separable

Directional wavelets and frames with actions of rotation or similitude groups ψ(b,a,θ)(x) =

1 a ψ( 1 a R−1 θ

  • x − b)
  • ,

where Rθ stands for the 2 × 2 rotation matrix Wf (b, a, θ) = ψ(b,a,θ), f inverted through f (x) = c−1

ψ

da a3

2π dθ

  • R2

d2b Wf (b, a, θ) ψ(b,a,θ)(x)

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-58
SLIDE 58

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 22/34

Directional, non-separable

Directional wavelets and frames:

◮ examples: Conical-Cauchy wavelet, Morlet-Gabor frames

Figure : Morlet Wavelet (real part) and Fourier representation

◮ possibility to decompose and reconstruct an image from a

discretized set of parameters; often (too) isotropic

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-59
SLIDE 59

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 23/34

Directional, anisotropic scaling

Ridgelets: 1-D wavelet and Radon transform Rf (θ, t) Rf (b, a, θ) =

  • ψ(b,a,θ)(x) f (x) d2x =
  • Rf (θ, t) a−1/2ψ((t−b)/a) dt

Figure : Ridgelet atom and GT decomposition

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-60
SLIDE 60

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 23/34

Directional, anisotropic scaling

Curvelet transform: continuous and frame

◮ curvelet atom: scale s, orient. θ ∈ [0, π), pos. y ∈ [0, 1]2:

ψs,y,θ(x) = ψs(R−1

θ (x − y))

ψs(x) ≈ s−3/4 ψ(s−1/2x1, s−1x2) parabolic stretch; (w ≃ √ l) Near-optimal decay: C 2 in C 2: O(n−2 log3 n)

◮ tight frame: ψm(x) = ψ2j,θℓ,xn(x)

where m = (j, n, ℓ) with sampling locations: θℓ = ℓπ2⌊j/2⌋−1 ∈ [0, π) and xn = Rθℓ(2j/2n1, 2jn2) ∈ [0, 1]2

◮ related transforms: shearlets, type-I ripplets

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-61
SLIDE 61

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 23/34

Directional, anisotropic scaling

Curvelet transform: continuous and frame

Figure : A curvelet atom and the wegde-like frequency support

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-62
SLIDE 62

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 23/34

Directional, anisotropic scaling

Curvelet transform: continuous and frame

Figure : GT curvelet decomposition

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-63
SLIDE 63

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 23/34

Directional, anisotropic scaling

Contourlets: Laplacian pyramid + directional filter banks

Figure : Contourlet atom and frequency tiling

from close to critical to highly oversampled

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-64
SLIDE 64

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 23/34

Directional, anisotropic scaling

Contourlets: Laplacian pyramid + directional filter banks

Figure : Contourlet GT (flexible) decomposition

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-65
SLIDE 65

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 23/34

Directional, anisotropic scaling

Shearlets

Figure : Shearlet atom in space and frequency, and frequency tiling

Do they have it all?

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-66
SLIDE 66

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 23/34

Directional, anisotropic scaling

Additional transforms

◮ previously mentioned transforms are better suited for edge

representation

◮ oscillating textures may require more appropriate transforms ◮ examples:

◮ wavelet and local cosine packets ◮ best packets in Gabor frames ◮ brushlets [Meyer, 1997; Borup, 2005] ◮ wave atoms [Demanet, 2007]

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-67
SLIDE 67

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 24/34

Lifting representations

Lifting scheme is an unifying framework

◮ to design adaptive biorthogonal wavelets ◮ use of spatially varying local interpolations ◮ at each scale j, aj−1 are split into ao j and do j ◮ wavelet coefficients dj and coarse scale coefficients aj: apply

(linear) operators Pλj

j

and Uλj

j

parameterized by λj dj = do

j − Pλj j ao j

and aj = ao

j + Uλj j dj

It also

◮ guarantees perfect reconstruction for arbitrary filters ◮ adapts to non-linear filters, morphological operations ◮ can be used on non-translation invariant grids to build

wavelets on surfaces

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-68
SLIDE 68

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 24/34

Lifting representations

dj = do

j − Pλj j ao j

and aj = ao

j + Uλj j dj

Lazy Predict Update n = m − 2j−1 m m + 2j−1 aj−1[n] aj−1[m] ao

j[n]

do

j[m]

dj[m] aj[n]

−1 2 −1 2 1 4 1 4

Gj−1 Gj ∪ Cj = ∪

Figure : Predict and update lifting steps; MaxMin lifting of GT

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-69
SLIDE 69

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 24/34

Lifting representations

Extensions and related works

◮ adaptive predictions:

◮ possibility to design the set of parameter λ = {λj}j to adapt

the transform to the geometry of the image

◮ λj is called an association field, since it links a coefficient of ao

j

to a few neighboring coefficients in do

j

◮ each association is optimized to reduce the magnitude of

wavelet coefficients dj, and should thus follow the geometric structures in the image

◮ may shorten wavelet filters near the edges

◮ grouplets: association fields combined to maintain

  • rthogonality
  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-70
SLIDE 70

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 25/34

Images are colors, not monochrome!

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-71
SLIDE 71

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 25/34

Images are colors, not monochrome!

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-72
SLIDE 72

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 25/34

Images are colors, not monochrome!

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-73
SLIDE 73

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 25/34

Images are colors, not monochrome!

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-74
SLIDE 74

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 25/34

Images are colors, not monochrome!

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-75
SLIDE 75

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 25/34

Images are colors, not monochrome!

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-76
SLIDE 76

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 26/34

One result among many others

Context: multivariate Stein-based denoising, multi-spectral Different spectral bands [Chaux et al., IEEE TSP, 2008]

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-77
SLIDE 77

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 27/34

One result among many others

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-78
SLIDE 78

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 27/34

One result among many others

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-79
SLIDE 79

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 27/34

One result among many others

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-80
SLIDE 80

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 27/34

One result among many others

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-81
SLIDE 81

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 28/34

What else? Images are not (all) flat

Many multiscale designs have been transported, adapted to:

◮ meshes ◮ spheres ◮ two-sheeted hyperboloid and

paraboloid

◮ 2-manifolds (case dependent) ◮ big deal: data on graphs

see 300+ reference list!

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-82
SLIDE 82

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 29/34

Images are not (all) flat

On the sphere, motivated by astronomy and geosciences

◮ filtering (two-angle spherical param. ξ = (θ, ϕ) ∈ S2)

(f ⋆ h)(ρ) =

  • S2 h(ρ ξ)f (ξ) dµ(ξ),

where ρ ∈ SO(3) is a rotation (driven by three angles) applied to the point ξ ∈ S2 and dµ(ξ) = sin θdθdϕ.

◮ Fourier transform

ˆ fℓm = Yℓm, f =

  • S2 Y ∗

ℓm(ξ) f (ξ) dµ(ξ),

f (ξ) =

  • ℓ,m

ˆ fℓm Yℓm(ξ) with respect to orthonormal basis of spherical harmonics Y = {Yℓm(ξ) : ℓ 0, |m| ℓ}, i.e. the eigenvectors of the spherical Laplacian

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-83
SLIDE 83

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 30/34

Images are not (all) flat

◮ spherical Scale-Space: solution at at time τ > 0 initialized to f

f (ξ, τ) =

  • ℓ,m

ˆ fℓm(τ)Yℓm(ξ) with ˆ fℓm(τ) = ˆ fℓm e−ℓ(ℓ+1) τ and f (ξ, 0) = f (ξ)

◮ spectral Wavelets

◮ spherical wavelets ψa(ξ), Wf (a, ξ) = (f ∗ ψa)(ξ)

f (ξ′) = f +

  • R+
  • S2 Wf (a, ξ) ψa(ξ′ · ξ) da

a dξ

with f =

1 4π

  • S2 f (ξ) dµ(ξ)

◮ MRA on the sphere with QMF ◮ needlets

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-84
SLIDE 84

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 30/34

Images are not (all) flat

Some of the other (many) constructions

◮ stereographic wavelets ◮ Haar transforms (embedded spherical triangulations)

j = 5 j = 4 j = 3 j = 2

Figure : Projection on the spherical Haar multiresolution

◮ steerable wavelets the sphere ◮ HEALPix: trans. wavelets, ridgelets, curvelets on the sphere ◮ SOHO wavelets, FLaglets, S2lets

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-85
SLIDE 85

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 30/34

Images are not (all) flat

Some of the other (many) constructions

◮ wavelets on general 2-manifolds ◮ lifting scheme on meshed surfaces

Figure : Multiresolution semi-regular mesh

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-86
SLIDE 86

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 30/34

Images are not (all) flat

Some of the other (many) constructions

◮ wavelets on general 2-manifolds ◮ lifting scheme on meshed surfaces

Figure : Surface approximation by thresholding (100%, 10%, 5%, 2%)

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-87
SLIDE 87

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 31/34

Images are not images anymore

Figure : Image as a graph

Pixels in images may be viewed as loci:

◮ with connectivity information (even non-local) ◮ scalar, vector-valued ◮ labeled

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-88
SLIDE 88

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 32/34

Images are not images anymore

Figure : Wavelets on a swiss roll

Wavelets on graphs [Hammond et al., 2011]

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-89
SLIDE 89

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 33/34

Conclusion: on a (frustrating) panorama

Take-away messages anyway? If you only have a hammer, every problem looks like a nail

◮ Is there a "best" geometric and multiscale transform?

◮ no: intricate data/transform/processing relationships ◮ more needed on asymptotics, optimization, models ◮ maybe: many candidates, progresses awaited: ◮ “so ℓ2”! Low-rank (ℓ0/ℓ1), math. morph. (+, × vs max, +) ◮ yes: those you handle best, or on wishlist ◮ mild redundancy, invariance, manageable correlation, fast

decay, tunable frequency decomposition, complex or more

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-90
SLIDE 90

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 33/34

Conclusion: on a (frustrating) panorama

Postponed references & toolboxes

◮ A Panorama on Multiscale Geometric Representations, Intertwining

Spatial, Directional and Frequency Selectivity Signal Processing, Dec. 2011 Toolboxes, images, and wavelet names

http://www.sciencedirect.com/science/article/pii/S0165168411001356 http://www.laurent-duval.eu/siva-panorama-multiscale-geometric-representations.html http://www.laurent-duval.eu/siva-wits-where-is-the-starlet.html

Cymatiophilic/leptostatonymomaniac acknowledgments to:

◮ the many *-lets. Last picks: Speclets/Gabor

shearlets/α-curvelets

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-91
SLIDE 91

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 34/34

Local applications

Weak interest in Lenna-like images

Figure : Geophysics: seismic Lenna

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-92
SLIDE 92

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 34/34

Local applications

Offset (traces) Time (smpl)

50 100 150 200 250 300 100 200 300 400 500 600 700

Figure : Geophysics: surface wave removal (before)

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-93
SLIDE 93

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 34/34

Local applications

Offset (traces) Time (smpl)

50 100 150 200 250 300 100 200 300 400 500 600 700

(b)

Figure : Geophysics: surface wave removal (after)

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-94
SLIDE 94

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 34/34

Local applications

Figure : Geophysics: paper to seismics

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-95
SLIDE 95

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 34/34

Local applications

Figure : Catalyst measurements

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)

slide-96
SLIDE 96

Motivations Intro. Early days Oriented & geometrical Far away from the plane End Example 34/34

Local applications

Figure : Process: catalyst measurements

  • L. Jacques, L. Duval, C. Chaux, G. Peyré:

UCL, IFPEN, AMU, Dauphine Curvelets, contourlets, shearlets, *lets, etc.: multiscale analysis and directional wavelets for images (and more)