
Curvelet-based Model for the Generation of Anisotropic Fractional Brownian Fields A PREPRINT
Ponson et al [
Ponson et al.(2006)Ponson, Bonamy, Auradou, Mourot, Morel, Bouchaud, Guillot, and Hulin
]
have studied the self-affine properties of anisotropic fracture surfaces. Jennane et al
[
Jennane et al.(2001)Jennane, Ohley, Majumdar, and Lemineur
] have observed anisotropic correlation in images of
bone X-ray tomographic microscopy projections.
Since the early work of Mandelbrot [
Mandelbrot and Van Ness(1968)
], the Fractional Brownian Motion model (FBM)
is largely being used to treat a wide variety of natural phenomena having self-similarity and long-range correla-
tions. In particular, the fractal analysis study applied to bidimensional fields has proved to be very useful to de-
scribe some physical properties such as roughness and porosity. However, the standard FBM model is only suitable
for describing materials and media that show isotropic symmetry. In that context, generalizations of the FBM to
anisotropic models have been proposed in last years. Kamont [
Kamont(1995)
] studied the Fractional Anisotropic
Wiener Field. Bonami and Estrade [
Bonami and Estrade(2003)
] analysed the properties of various models of anisotropic
Gaussian fields and proposed a new procedure to characterize the anisotropy of such fields. Tavares and Lucena
[
Tavares and Lucena(2003)
] and Heneghan [
Heneghan et al.(1996)Heneghan, Lowen, and Teich
] have implemented
wavelet-based models for anisotropic hypersurfaces. Biermé and Richard [
Biermé and Richard(2008)
] have applyied
Radon Transform to estimate the anisotropy of Fractional Brownian Fields. Xiao [
Xiao(2009)
] has studied sample path
properties of such anisotropic fields. It will be desirable to generalize these concepts and ideas in simple and automatic
ways that can allow sparse representations of highly complex anisotropic fields.
In this paper, we introduce a new method for generating 2-D Anisotropic Fractional Brownian Fields (AFBF)
based on the Curvelet Transform. This is a new multiscale transform with strong directional character that
provides an optimal representation of objects that have discontinuities along edges [
Candes and Donoho(2000)
,
Starck et al.(2002)Starck, Candès, and Donoho
,
Candès and Donoho(2004)
]. The curvelets are localized not only in
the spatial domain (location) and frequency domain (scale), but also in angular orientation, which is a step ahead
compared to Wavelet Transform [
Stephane(1999)
]. This directional feature is the one we use to achieve the anisotropy.
2 The Fractional Brownian Motion
We start with the model of Fractional Brownian Motion proposed by Mandelbrot [
Mandelbrot and Van Ness(1968)
].
An isotropic two-dimensional Fractional Brownian field (FBF)
BH(u)
,
u∈R2
, with Hurst index
H
taking values in
(0,1), is defined by the correlation function
E[BH(u)BH(v)] = CH|u|2H+|v|2H+|u−v|2H,∀u, v ∈R2,(1)
where
E[·]
is the expectation operator and
CH
is a constant depending on
H
.
BH(u)
is not a stationary process, but its
increments form a stationary, zero-mean gaussian process, with variance depending only on the distance ∆u:
Eh|BH(u+ ∆u)−BH(u)|2i∝∆uH.
It follows from (1) that BHis a self-similar field:
BH(λu)d
=λHBH(u),
where
d
=
means equal in distribution, and
λ > 0
is a constant. It can be proved that
BH
has an average spectral density
of the form:
S(ξ)∝|ξ|−2H−2(2)
in which
ξ= (ξ1, ξ2)
are the frequency coordinates in the Fourier domain. As proposed in [
Bonami and Estrade(2003)
],
anisotropic fields with stationary increments can be defined by taking orientation-dependent spectral densities of the
form:
S(ξ)∝|ξ|−2Hθ−2
where the Hurst index
Hθ
now depends on the direction of
ξ
, labeled by the index
θ
. We consider
H
a parameter related
to the roughness and regularity of the surface. So, in the case of anisotropy, we would expect that
Hθ
would be related
to the self-similarity properties of the field on a given direction. A 1-D process obtained by selecting any straight line of
an isotropic FBF
BH
is a FBM with Hurst index
H
. However, contrary to what one would expect, the estimation of the
orientation-dependent Hurst index
H
in an anisotropic Gaussian field cannot be performed simply by analysing sample
lines of the field, as in the isotropic case, since the regularity along a line do not appear to be dependent on the direction.
In order to measure the anisotropy dependence of Hurst exponent
Hθ
, Bonami and Strade [
Bonami and Estrade(2003)
]
2