(1) |
iterate {
(2) | ||
(3) | ||
(4) |
} ,where the subscript k denotes the iteration number.
To greatly improve efficiency, we solve equation (3) in the Fourier domain. We apply the divergence to the estimated dips and divide by the Z-transform of the discretized Laplacian in the Fourier domain as
(5) |
Unfortunately, mirroring is required to avoid boundary affects in the Fourier domain. Mirroring is concatenating reversed copies of the data along each dimension so that the input to the tranform has periodic boundaries. In 1D the input is increased by a factor of 2, in 2D the input is increased by a factor of four, and so on. In this case, we mirror the divergence of the residual (). Mirroring is described in more detail by Ghiglia and Pritt (1998) for application to 2D phase unwrapping.
Ghiglia and Pritt (1998) also describe an alternative method using discrete cosine transforms that does not need mirroring. This method does not require mirroring because the cosine transform assumes that image is even (periodic) implicitly. Since it makes this assumption, then only the unmirrored input itself must be transformed. The alternative formulation is
(6) |