15.01.2013 Views

U. Glaeser

U. Glaeser

U. Glaeser

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

FIGURE 28.17 The real (top) and imaginary (bottom) parts of a representative 2-D Gabor function.<br />

where ∆ x, ∆ y, ∆ξ x, and ∆ξ y are the effective widths of the functions in the spatial and spatial-frequency<br />

domains. By this measure, then, these functions are optimally local. Their real and imaginary parts also<br />

agree reasonably well with measured receptive field profiles. The basis is not orthogonal, however.<br />

Specifically, the Gabor transform is not equivalent to the finite-support Fourier transform with a Gaussian<br />

window. For a cross-section of the state of the art in Gabor transform-based analysis, see [12].<br />

The Derivative of Gaussian Transform<br />

In 1987, Young [13] proposed a receptive field model based on the Gaussian and its derivatives. These<br />

functions, like the Gabor functions, are spatially and spectrally local and consist of alternating regions<br />

of excitation and inhibition in a decaying envelope. Young showed that Gaussian derivative functions<br />

more accurately model the measured receptive field data than do the Gabor functions [14].<br />

In [15], a spatial/spatial-frequency representation based on shifted versions of the Gaussian and its<br />

derivatives was introduced (the derivative of Gaussian transform (DGT)). As with the Gabor transform,<br />

although this transform is nonorthogonal, with a suitably chosen basis it is invertible. The DGT has<br />

significant practical advantage over the Gabor transform in that both the basis functions and coefficients<br />

of expansion are real-valued.<br />

The family of 2-D separable Gaussian derivatives centered at the origin can be defined as<br />

© 2002 by CRC Press LLC<br />

g0,0( x, y)<br />

= g0( x)g0(<br />

y)<br />

=<br />

e (x2 +y 2 – )<br />

/2σ 2<br />

gm,n ( x, y)<br />

= gm( x)gn(<br />

y)<br />

d m ( )<br />

dx m ( )<br />

( n)<br />

d<br />

dy n ( )<br />

= -----------g 0( x)<br />

----------g 0( y)<br />

(28.35)<br />

(28.36)<br />

This set can then be shifted to any desired location. The variance σ defines the extent of the functions<br />

in the spatial domain. There is an inverse relationship between the spatial and spectral extents, and the<br />

value of this variable may be constant or may vary with context.<br />

The 1-D Gaussian derivative function spectra are bimodal (except for that of the original Gaussian,<br />

which is itself a Gaussian) with modes centered at ±Ω m rad/pixel:<br />

m<br />

Ωm =<br />

-------σ<br />

(28.37)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!