Abstract

Since our algorithm for reconstructing generating currents from magnetic-field measurements uses a vector-Fourier space description of the imaging process, we cast the simulation in these terms. Using the Maxwell equations for static fields, we synthesize the magnetic field which would be generated by a current density distribution. Because we use the discrete Fourier transform, care is taken to configure the current so that the continuity equation is satisfied. We also take into account the spatial-frequency warping and the special form of bounding functions that are required by the discretization. Once the field has been synthesized, a forbidden measurement region is blocked out (including the region occupied by the current source). Multiplicative noise in controlled quantity is mixed into the accessible measurement region. Our reconstruction algorithm is applied. The resulting image or image series is then displayed on a color video display unit.

© 1985 Optical Society of America