We present a method for imaging small scatterers in a homogeneous medium from polarization measurements of the electric field at an array. The electric field comes from illuminating the scatterers with a point source with known location and polarization. We view this problem as a generalized phase retrieval problem with data being the coherency matrix or Stokes parameters of the electric field at the array. We introduce a simple preprocessing of the coherency matrix data that partially recovers the ideal data where all the components of the electric field are known for different source dipole moments. We prove that the images obtained using an electromagnetic version of Kirchhoff migration applied to the partial data are, for high frequencies, asymptotically identical to the images obtained from ideal data. We analyze the image resolution and show that polarizability tensor components in an appropriate basis can be recovered from the Kirchhoff images, which are tensor fields. A time domain interpretation of this imaging problem is provided and numerical experiments are used to illustrate the theory.
A key aspect of our imaging method is that we do not need to retrieve all phases from the data, since the imaging method we use does not require it. This is a feature that is also used for intensity only imaging in @cite_27 @cite_0 @cite_20 @cite_1 @cite_28 . For full aperture data, phase retrieval is done in @cite_19 @cite_12 @cite_6 using a similar preprocessing of the intensity data as we present here. One novelty in our method is that we image polarizability of scatterers, directly from polarization data. There are other ways of imaging polarizability in different physical setups such as Optical Coherence Tomography from interferometric data @cite_9 and Polarimetric Synthetic Aperture Radar (see e.g. @cite_24 ). Finally we point out that our imaging method can be seen as correlation based imaging, see e.g. @cite_16 @cite_22 @cite_38 @cite_2 .
@cite_38 @cite_22 @cite_28 @cite_9 @cite_1 @cite_6 @cite_0 @cite_19 @cite_27 @cite_24 @cite_2 @cite_16 @cite_12 @cite_20
2963789789 2557428315 2082344075 2036379242
In this paper we consider narrow band, active array imaging of weak localized scatterers when only the intensities are recorded at an array with @math transducers. We assume that the medium is homogeneous and, hence, wave propagation is fully coherent. This work is an extension of our previous paper [A. Novikov, M. Moscoso, and G. Papanicolaou, SIAM J. Imaging Sci., 8 (2015), pp. 1547--1573], where we showed that using linear combinations of intensity-only measurements, obtained from @math illuminations, imaging of localized scatterers can be carried out efficiently using imaging methods based on the singular value decomposition of the time-reversal matrix. Here we show that the same strategy can be accomplished with only @math illuminations, therefore enormously reducing the data acquisition process. Furthermore, we show that in the paraxial regime one can form the images by using six illuminations only. In particular, this paraxial regime includes Fresnel and Fraunhofer diffraction. The key point of this... We propose an illumination strategy for interferometric imaging that allows for robust depth recovery from intensity-only measurements. For an array with colocated sources and receivers, we show that all the possible interferometric data for multiple sources, receivers, and frequencies can be recovered from intensity-only measurements provided that we have sufficient source location and frequency illumination diversity. There is no need for phase reconstruction in this approach. Using interferometric imaging methods we show that in homogeneous media there is no loss of resolution when imaging with intensities only. If in these imaging methods we reduce incoherence by restricting the multifrequency interferometric data to nearby array elements and nearby frequencies we obtain robust images in weakly inhomogeneous background media with a somewhat reduced resolution. Tomographic image reconstruction, such as the reconstruction of computed tomography projection values, of tomosynthesis data, positron emission tomography or SPECT events, and of magnetic resonance imaging data is computationally very demanding. One of the most time-consuming steps is the backprojection. Recently, a novel general purpose architecture optimized for distributed computing became available: the cell broadband engine (CBE). To maximize image reconstruction speed we modified our parallel-beam backprojection algorithm [two dimensional (2D)] and our perspective backprojection algorithm [three dimensional (3D), cone beam for flat–panel detectors] and optimized the code for the CBE. The algorithms are pixel or voxel driven, run with floating point accuracy and use linear (LI) or nearest neighbor (NN) interpolation between detector elements. For the parallel-beam case, 512 projections per half rotation, 1024 detector channels, and an image of size 512 2 was used. The cone-beam backprojection performance was assessed by backprojecting a full circle scan of 512 projections of size 1024 2 into a volume of size 512 3 voxels. The field of view was chosen to completely lie within the field of measurement and the pixel or voxel size was set to correspond to the detector element size projected to the center of rotation divided by 2 . Both the PC and the CBE were clocked at 3 GHz . For the parallel backprojection of 512 projections into a 512 2 image, a throughput of 11 fps (LI) and 15 fps (NN) was measured on the PC, whereas the CBE achieved 126 fps (LI) and 165 fps (NN), respectively. The cone-beam backprojection of 512 projections into the 512 3 volume took 3.2 min on the PC and is as fast as 13.6 s on the cell. Thereby, the cell greatly outperforms today’s top-notch backprojections based on graphical processing units. Using both CBEs of our dual cell-based blade (Mercury Computer Systems) allows to 2D backproject 330 images s and one can complete the 3D cone-beam backprojection in 6.8 s . We address the problem of recovering a scene recorded through a semireflecting medium (i.e. planar lens), with a virtual reflected image being superimposed on the image of the scene transmitted through the semirefelecting lens. Recent studies propose imaging through a linear polarizer at several orientations to estimate the reflected and the transmitted components in the scene. In this study we extend the sparse ICA (SPICA) technique and apply it to the problem of separating the image of the scene without having any a priori knowledge about its structure or statistics. Recent novel advances in the SPICA approach are discussed. Simulation and experimental results demonstrate the efficacy of the proposed methods.© 2005 Wiley Periodicals, Inc. © 2005 Wiley Periodicals, Inc. Int J Imaging Syst Technol, 15, 84–91, 2005; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002 ima.20042