The helical coordinate

Next: Improving low-frequency behavior Up: The helical coordinate Previous: FACTORED LAPLACIAN == HELIX

# HELIX LOW-CUT FILTER

Because the autocorrelation of is is a second derivative, the operator must be something like a first derivative. As a geophysicist, I found it natural to compare the operator with by applying the helix derivative to a local topographic map. The result shown in Figure 11 is that enhances drainage patterns whereas enhances mountain ridges.

helocut
Figure 11.
Topography, helical derivative, slope south.

The operator has curious similarities and differences with the familiar gradient and divergence operators. In 2-dimensional physical space, the gradient maps one field to two fields (north slope and east slope). The factorization of with the helix gives us the operator that maps one field to one field. Being a one-to-one transformation (unlike gradient and divergence), the operator is potentially invertible by deconvolution (recursive filtering).

I have chosen the name helix derivative'' or helical derivative'' for the operator . A flag pole has a narrow shadow behind it. The helix integral (middle frame of Figure ) and the helix derivative (left frame) show shadows with an angular bandwidth approaching .

Our construction makes have the energy spectrum , so the magnitude of the Fourier transform is . It is a cone centered at the origin with there the value zero. By contrast, the components of the ordinary gradient have amplitude responses and that are lines of zero across the -plane.

The rotationally invariant cone in the Fourier domain contrasts sharply with the nonrotationally invariant helix derivative in -space. The difference must arise from the phase spectrum. The factorization (13) is nonunique because causality associated with the helix mapping can be defined along either - or -axes; thus the operator (13) can be rotated or reflected.

In practice, we often require an isotropic filter. Such a filter is a function of . It could be represented as a sum of helix derivatives to integer powers.

If you want to see some tracks on the side of a hill, you want to subtract the hill and see only the tracks. Usually, however, you do not have a very good model for the hill. As an expedient, you could apply a low-cut filter to remove all slowly variable functions of altitude. In Chapter we found the Sea of Galilee in Figure to be too smooth for viewing pleasure, so we made the roughened versions in Figure , a 1-dimensional filter that we could apply over the -axis or the -axis. In Fourier space, such a filter has a response function of or a function of . The isotropy of physical space tells us it would be more logical to design a filter that is a function of . In Figure 11 we saw that the helix derivative does a nice job. The Fourier magnitude of its impulse response is . There is a little anisotropy connected with phase (which way should we wind the helix, on or ?), but it is not nearly so severe as that of either component of the gradient, the two components having wholly different spectra, amplitude or .

Subsections
 The helical coordinate

Next: Improving low-frequency behavior Up: The helical coordinate Previous: FACTORED LAPLACIAN == HELIX

2015-03-25