arxiv: v1 [cs.lg] 11 Mar 2016
|
|
- Amanda French
- 5 years ago
- Views:
Transcription
1 A Recursive Born Approach to Nonlinear Inverse Scattering Ulugbek S. Kamilov, Dehong Liu, Hassan Mansour, and Petros T. Boufounos March 15, 2016 arxiv: v1 [cs.lg] 11 Mar 2016 Abstract The Iterative Born Approximation (IBA) is a well-known method for describing waves scattered by semi-transparent objects. In this paper, we present a novel nonlinear inverse scattering method that combines IBA with an edge-preserving total variation (TV) regularizer. The proposed method is obtained by relating iterations of IBA to layers of a feedforward neural network and developing a corresponding error backpropagation algorithm for efficiently estimating the permittivity of the object. Simulations illustrate that, by accounting for multiple scattering, the method successfully recovers the permittivity distribution where the traditional linear inverse scattering fails. 1 Introduction Knowledge of the spatial distribution of the permittivity within an object is important for many applications since it enables the visualization of the internal structure and physical properties of the object. Measurements of the permittivity are typically obtained by first illuminating the object with a known incident wave and recording the resulting scattered waves with sensors located outside the object. The spatial map of permittivity is reconstructed from the measurements, using computational inverse scattering methods that rely on physical models describing the object-wave interaction. Traditional approaches to inverse scattering formulate the task as a linear inverse problem by establishing a linear relationship between the permittivity and the scattered wave. The linear model can be obtained by assuming a straight-ray propagation of waves [1], or by adopting more refined single-scattering models based on the first Born [2] or Rytov approximations [3]. Once linearized, the problem can be efficiently solved using the preferred regularized reconstruction methods, typically based on sparsity and iterative optimization [4 6]. Recent experimental results indicate that the resolution and quality of the reconstructed permittivity is improved when nonlinear physical models are used instead of traditional linear ones [7 10]. In particular, nonlinear physical models can account for multiple scattering and provide a more accurate interpretation of the measured data at the cost of increased computational complexity of the reconstruction. In this paper, we develop a new computational imaging method to reconstruct the permittivity distribution of an object from transmitted or reflected waves. Our method is based on a nonlinear physical model that can account for multiple scattering in a computationally efficient way. Specifically, we propose to interpret the iterations of the iterative Born approximation (IBA) as layers of a feedforward neural network. This formulation leads to an efficient error backpropagation algorithm used to evaluate the gradient of the scattered field with respect to the permittivity, thus enabling the recovery of the latter from a set of measured scattered fields. The quality of the final estimate is further enhanced by regularizing the solution with an edge-preserving total variation (TV) penalty. Our simulations indicate that the proposed method accurately models scattering without prohibitive computational overhead, and successfully recovers the object where the traditional linear approaches fail. U. S. Kamilov ( kamilov@merl.com), D. Liu ( liudh@merl.com), H. Mansour ( mansour@merl.com), and P. T. Boufounos ( petrosb@merl.com) are with Mitsubishi Electric Research Laboratories (MERL), 201 Broadway, Cambridge, MA 02140, USA. 1
2 y u in f(x) x imaging area sensor region Figure 1: Schematic representation of a scattering scenario. An object with a scattering potential f(x), x Ω, is illuminated with an input wave u in, which interacts with the object and results in the scattered wave u sc at the sensor region Γ. The scattered wave is measured and used to computationally form an image of f. 2 Related Work A comprehensive review of nonlinear inverse scattering, with detailed description of standard algorithms, is available in the book-chapter by van den Berg [11]. Three common approaches to nonlinear inverse scattering are iterative Born [12 14], modified gradient [15 17], and contrast-source inversion methods [18]. All these methods attempt to iteratively minimize the discrepancy between the actual and predicted measurements, while enforcing the consistency of the fields inside the object. The actual optimization is performed in an alternating fashion by first updating the permittivity for a fixed field, and then updating the field for a fixed permittivity. The difference between the methods is in the actual computation of the updates. Recently, the beam propagation method (BPM) was proposed for performing nonlinear inverse scattering in transmission [10, 19 21]. BPM-based methods circumvent the need to solve an explicit optimization problem for the internal field, by instead numerically propagating the field slice-by-slice through the object. It was shown that BPM can be related to a neural network where each layer corresponds to a spatial slice of the object. An efficient backpropagation algorithm can thus be derived to reconstruct the object from the measurements [21]. This paper substantially extends the BPM-based method by relying on IBA to formulate the nonlinear physical model. One main advantage of the proposed formulation over BPM is that it allows reconstruction when measuring reflections. Additionally, unlike alternating minimization schemes that assume a linear problem for a given internal field, our method directly optimizes the nonlinear model in a tractable way using error backpropagation. 3 Main Results The method presented here can be generalized to a majority of tomographic experiments in transmission or reflection. For simplicity of derivations, we ignore absorption by assuming a real permittivity; however, the equations can be generalized to handle complex permittivities. We additionally assume coherent measurements, i.e., that both the amplitude and phase of the scattered wave are recorded at the sensor locations. 2
3 3.1 Problem Formulation Consider the two-dimensional (2D) scattering problem illustrated in Fig. 1, where an object of real permittivity distribution ɛ(x), with x = (x, y) Ω, is immersed into the background medium of permittivity ɛ b. The line sources that generate the electromagnetic excitation and the sensors collecting the data are located in the sensor region Γ R 2. By assuming, and ignoring, a time dependence exp(jωt), the incident electric field created by the lth source, located at x l Γ, is given by u in (x) = A j 4 H(2) 0 (k b x x l l2 ), (1) for all x R 2, where A is the strength of the source, H (2) 0 is the zero-order Hankel function of the second kind, k b = k 0 ɛb is the wavenumber in the background medium, k 0 = ω/c 0 = 2π/λ is the wavenumber in free space, λ is the wavelength, and c m/s. In the subsequent derivations, we consider the scenario of a single illumination and drop the indices l. The generalization to an arbitrary number of illuminations L is straightforward. The Lippmann-Schwinger integral equation describes the relationship between the permittivity and the wave-field [11] u(x) = u in (x) + g(x x )f(x )u(x ) dx, (2) Ω for all x Ω, where we define the scattering potential and the Green s function for the homogeneous medium f(x) k 2 b (ɛ b ɛ(x)) (3) g(x) j 4 H(2) 0 (k b x l2 ). (4) Similarly, the scattered field in the sensor region can be expressed as u sc (x) = u(x) u in (x) = g(x x )f(x )u(x ) dx (5) for any x Γ. Note that the integrals (2) and (5) extend only over Ω because the scattering potential f is zero for all x / Ω. The goal of inverse scattering is to estimate f, which is equivalent to ɛ, given M measurements of {u sc (x m )} m [1...M] in the sensor region Γ. 3.2 Iterative Born Approximation At first glance, it might seem that (5) directly provides a linear relationship between f and u sc, which can be used to solve the problem. However, the nonlinear nature of the inverse scattering becomes evident if one realizes that the internal field u = u in + u sc in (5) depends on u sc. We now consider a K-term iterative Born approximation (IBA) [22] of the total field (2) u k (x) = u in (x) + g(x x )f(x )u k 1 (x ) dx, (6) Ω where x Ω, u 0 = 0, and k = 1, 2,..., K. When K = 1, eq. (6) reduces to the well-known first-born approximation, which assumes a single scattering from f by approximating u(x) with the incident field u in (x). For K = 2, the approximation of the total field is improved by taking into account the second scattering due to an additional interaction between the object and the field. For higher values of K the approximation is further improved by accounting for multiple scattering of order K (see also the discussion on multiple scattering by Born and Wolf [22]). Ω 3
4 incident field u 0 1 f 1 u 1 u 2 f 1 f 1 error correction predicted field measured field z 1 y 1 u 0 2 f 2 f 2 f 2 compare error z 2 y 2 u 0 3 f 3 scatter f 3 scatter f 3 Figure 2: Schematic representation of the multiple scattering as a neural network for K = 2, N = 3, and M = 2. The scattering potential f R N plays the role of nonlinearities, while matrices H C M N and G C N N represent weights of the network. The error backpropagation algorithm described here allows to efficiently estimate f by comparing predicted field z C M against the actual measurements y C M. In the context of general theory of integral equations, the iterations of type (6) are known as Liouville- Neumann series [23]. The sufficient condition for convergence states that the norm of the integral operator acting on the field should be less than unity. While this implies that IBA might diverge for large permittivity contrasts, it is still expected to work on a wide range of contrasts where the linearized models fail, which is corroborated by our simulations in Section Network Interpretation We now discretize and combine equations (5) and (6) into the following matrix-vector recursion z H(u K f), u k u 0 + G(u k 1 f), (7a) (7b) for k = 1,..., K. Here, the vector f R N is the discretization of the scattering potential f, z C M is the predicted scattered field u sc at sensor locations {x m } m [1...M], u 0 C N is the discretization of the input field u in inside Ω, H C M N is the discretization of the Green s function at sensor locations, G C N N is the discretization of the Green s function inside Ω, and denotes a component-wise multiplication between two vectors. For every k [1... K], the vector u k C N denotes the discretized version of the internal field after the kth scattering. Figure 2 illustrates the representation of (7) as a feedforward neural network [24], where the edge weights are represented in H and G and the nonlinear nodes are described by the scattering potential f. Note that the linear edge weights correspond to convolution operators and can thus be efficiently implemented with FFTs. Accordingly, the total computational cost of evaluating one forward pass through the network is O(KN log(n)). 3.4 Inverse Scattering We formulate the inverse scattering as the following minimization problem f = arg min {D(f) + τr(f)}, (8) f F where D and R are the data-fidelity and regularization terms, respectively, and τ > 0 is the regularization parameter. The convex set F R N enforces physical constraints on the scattering potential such as its, for example, non-negativity. The data-fidelity term is given by D(f) 1 2 y z(f) 2 l 2, 4
5 where y C M contains measurements of the scattered field and z is the field predicted by the recursion (7). As a regularization term, we propose to use isotropic TV penalty [25] R(f) N [Df] n l2 = n=1 N n=1 [D x f] n 2 + [D y f] n 2, where D : R N R N 2 is the discrete gradient operator with matrices D x and D y denoting the finite difference operations along x and y directions, respectively. The optimization (8) can be performed iteratively using a proximal-gradient scheme or one of its accelerated variants [26 28]. Specifically, the scattering potential can be updated with the following iteration f t prox γτr ( f t 1 γ D(f t 1 ) ), (9) where γ > 0 is a step-size and { } 1 prox τr (g) arg min f F 2 f g 2 l 2 + τr(f) (10) is the proximal operator, which corresponds to the TV regularized solution of the denoising problem. Note that, although, the proximal operator for isotropic TV does not admit a closed form, it can be efficiently computed [27]. The gradient D can be expressed as follows { [ ] H D(f) = Re z(f) (z(f) y)}, (11) where H is the Hermitian transpose of the Jacobian z(f) [ z 1... z N ]. (12) Then, by differentiating equations in (7) with respect to f, and simplifying the resulting expressions, we obtain for any two vectors r C M and b C N [ ] H z r = ( H H r ) [ ] u ū K K H (( + H H r ) f ) [ ] u k H b = ( G H b ) [ ] u ū k 1 k 1 H (( + G H b ) f ), where k = 1,..., K, and vector ū contains complex conjugated elements of u. These relationships lead to the following error backpropagation algorithm g k g k+1 + [ G H v k+1] ū k v k [ G H v k+1] f, (13a) (13b) where k = K 1, K 2,..., 0, with the initialization v K = [H H (z y)] f and g K = [ H H (z y) ] ū K. The final value of the gradient (11) is obtained by returning D(f) = Re{g 0 }. The remarkable feature of our error backpropagation approach is that it allows to efficiently evaluate the gradient of the scattered field with respect to the scattering potential. Due to the convolutional structure of the matrices, its computational complexity is equivalent to running a forward pass, which is of order O(KN log(n)). Equipped with this algorithm, the scattering potential can be optimized via iteration (9). Note that the algorithm does not explicitly evaluate and store the Jacobian (12) by instead computing its product with the residual r = (z(f) y), as indicated in (11). 5
6 a b 1.65 db SNR (db) FB #layers (K) 32 FB c db d db RB Figure 3: Evaluation on Shepp-Logan with 15% permittivity contrast at λ = 7.49 cm. (a) True contrast; (b) first-born approximation; (c) alternating minimization; (d) proposed recursive-born method with K = 32. Scale bar is equal to λ. 4 Experimental Validation To validate our Recursive Born (RB) method, we report results for the tomographic experiment illustrated in Fig. 1, where the scattered wave measurements were obtained by running a high-fidelity Finite-Difference Time-Domain (FDTD) [29] simulator. The object is the Shepp-Logan phantom of size 82 cm 112 cm and the background medium is air with ɛ b = 1. The measurements are collected over 24 transmissions on a circle of radius R = 100 cm with 15 o angle increments and, for each transmission, 360 measurements around the object are recorded. The dimensions of the computational domain for reconstruction are set to L x = L y = 120 cm, with sampling steps δx = δy = 0.6 cm, which implies a measurement ratio of M/N 22%. We define the permittivity contrast as f max (ɛ max ɛ b )/ɛ b, where ɛ max max x Ω {ɛ(x)}. We compare results of our approach against two alternative methods. As the first reference method (denoted first Born (FB)), we consider the TV-regularized solution of linearized model based on the first- Born approximation, which is known to be valid only for weakly scattering objects. In addition to the linearized approach, we consider an optimization scheme (denoted Alternating Minimization ()) that alternates between updating the scattering potential for a fixed field and updating the field for a fixed scattering potential. This method is conceptually similar to the one proposed in [13], but with TV replacing the smooth regularizer, due to the edge-preserving properties of the former. Hence, all three methods minimize the same TV-regularized least-squares error functional; however, each method relies on a distinct physical forward model. Additionally, all the methods were initialized with zero, iterated until convergence by measuring the change in two successive estimates, and their regularization parameters were set to the same value τ = y 2 l 2. Figure 3 compares the quality of the images obtained by all three methods for the permittivity contrast of 15% at λ = 7.49 cm. Note that, due to the object s large size and contrast, FB fails to characterize its structure. Alternatively, both and RB succeed at recovering the object, with RB obtaining an image of significantly superior quality. Figure 4 illustrates the influence of the number of layers K on the quality of the reconstructed image (left), and the evolution of the relative data-fit y z 2 l 2 / y 2 l 2 for every iteration 6
7 b 1.65 db SNR (db) FB data fit FB #layers (K) 32 1 iterations (t) Figure 4: Illustration of the reconstruction performance on Shepp-Logan with 15% contrast. Left: SNR is plotted against FBthe number of layers in the neural network. Right: normalized error between the true and predicted fields is plotted at each iteration for K = 32. Results obtained with FB and are marked with dashed lines db d db Table 1: Comparison of three methods in terms of SNR for λ = 10 cm. The best result for each contrast is highlighted. Permittivity Contrast f max 5% 10% 15% 20% RB First Born Alternating Minimization Recursive Born with K = 32 (right). As can be appreciated from these plots, the proposed method outperforms FB and, both in terms of signal-to-noise ratio (SNR) and data-fit for networks with sufficient number of layers. Additionally, the method converges relatively fast within first few tens of iterations. Table 1 presents the results of quantitative evaluation of the methods for different values of the permittivity contrast at λ = 10 cm. As expected, the performance of all the methods degrades as the contrast value increases, which might be due to growing degree of nonlinearity, and hence, nonconvexity of the inverse scattering problem [30]. However, the solutions computed by the proposed RB approach are substantially better than the two alternatives, FB and, for all values of the contrast. Finally, from a computational perspective, a single iteration of the method requires a number of FFTs proportional to the number of scattering layers. In our simulations, accurate results are obtained using about 100 iterations with networks of approximately 20 layers. More concretely, our basic MATLAB implementation requires about 1.3 seconds per iteration to process a transmission on a 4GHz Intel Core i7 processor with 32 GBs of memory. 5 Conclusion The method developed in this paper reconstructs the distribution of the permittivity in an object from a set of measured scattered waves. In particular, the method accounts for multiple scattering of waves, in both transmission and reflection, and can thus be used when linearized models fail. The method is also computationally tractable due to its convolutional structure and can be further accelerated by parallelizing computations over multiple CPUs. We believe that the approach presented here opens rich perspectives for high-resolution tomographic imaging in a range of practical setups where multiple scattering is an issue. References [1] A. C. Kak and M. Slaney, Principles of Computerized Tomographic Imaging. IEEE, [2] E. Wolf, Three-dimensional structure determination of semi-transparent objects from holographic data, Opt. Commun., vol. 1, no. 4, pp , September/October
8 [3] A. J. Devaney, Inverse-scattering theory within the Rytov approximation, Opt. Lett., vol. 6, no. 8, pp , August [4] M. M. Bronstein, A. M. Bronstein, M. Zibulevsky, and H. Azhari, Reconstruction in diffraction ultrasound tomography using nonuniform FFT, IEEE Trans. Med. Imag., vol. 21, no. 11, pp , November [5] Y. Sung and R. R. Dasari, Deterministic regularization of three-dimensional optical diffraction tomography, J. Opt. Soc. Am. A, vol. 28, no. 8, pp , August [6] T. Kim, R. Zhou, M. Mir, S. Babacan, P. Carney, L. Goddard, and G. Popescu, White-light diffraction tomography of unlabelled live cells, Nat. Photonics, vol. 8, pp , March [7] F. Simonetti, Multiple scattering: The key to unravel the subwavelength world from the far-field pattern of scattered wave, Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys., vol. 73, no. 3, p , March [8] F. Simonetti, M. Fleming, and E. A. Marengo, Illustration of the role of multiple scattering in subwavelength imaging from far-field measurements, J. Opt. Soc. Am. A, vol. 25, no. 2, pp , February [9] G. Maire, F. Drsek, J. Girard, H. Giovaninni, A. Talneau, D. Konan, K. Belkebir, P. C. Chaumet, and A. Sentenac, Experimental demonstration of quantitative imaging beyond Abbe s limit with optical diffraction tomography, Phys. Rev. Lett., vol. 102, p , May [10] U. S. Kamilov, I. N. Papadopoulos, M. H. Shoreh, A. Goy, C. Vonesch, M. Unser, and D. Psaltis, Learning approach to optical tomography, Optica, vol. 2, no. 6, pp , June [11] P. M. van den Berg, Scattering. Academic Press, 2002, ch. Nonlinear Scalar Inverse Scattering: Algorithms and Applications, pp [12] A. G. Tijhuis, Born-type reconstruction of material parameters of an inhomogeneous, lossy dielectric slab from reflected-field data, Wave Motion, vol. 11, no. 2, pp , May [13] Y. M. Wang and W. Chew, An iterative solution of the two-dimensional electromagnetic inverse scattering problem, Int. J. Imag. Syst Tech., vol. 1, pp , [14] W. C. Chew and Y. M. Wang, Reconstruction of two-dimensional permittivity distribution using the distorted Born iterative method, IEEE Trans. Med. Imag., vol. 9, no. 2, pp , June [15] R. E. Kleinman and P. M. van den Berg, A modified gradient method for two-dimensional problems in tomography, J. Comput. Appl. Math., vol. 42, no. 1, pp , [16], An extended range-modified gradient technique for profile inversion, Radio Sci., vol. 28, no. 5, pp , September-October [17] K. Belkebir and A. Sentenac, High-resolution optical diffraction microscopy, J. Opt. Soc. Am. A, vol. 20, no. 7, pp , July [18] P. M. van den Berg and R. E. Kleinman, A contrast source inversion method, Inv. Probl., vol. 13, no. 6, pp , December [19] L. Tian and L. Waller, 3D intensity and phase imaging from light field measurements in an LED array microscope, Optica, vol. 2, pp , [20] L. Waller and L. Tian, Machine learning for 3D microscopy, Nature, vol. 523, no. 7561, pp , July
9 [21] U. S. Kamilov, I. N. Papadopoulos, M. H. Shoreh, A. Goy, C. Vonesch, M. Unser, and D. Psaltis, Optical tomographic image reconstruction based on beam propagation and sparse regularization, IEEE Trans. Comp. Imag., vol. 2, no. 1, pp , March [22] M. Born and E. Wolf, Principles of Optics, 7th ed. Cambridge Univ. Press, 2003, ch. Scattering from inhomogeneous media, pp [23] G. B. Arfken and H. J. Weber, Mathematical Methods for Physicists, 6th ed. Elsevier, 2005, ch. Integral equations, pp [24] C. M. Bishop, Neural Networks for Pattern Recognition. Oxford, [25] L. I. Rudin, S. Osher, and E. Fatemi, Nonlinear total variation based noise removal algorithms, Physica D, vol. 60, no. 1 4, pp , November [26] J. M. Bioucas-Dias and M. A. T. Figueiredo, A new TwIST: Two-step iterative shrinkage/thresholding algorithms for image restoration, IEEE Trans. Image Process., vol. 16, no. 12, pp , December [27] A. Beck and M. Teboulle, Fast gradient-based algorithm for constrained total variation image denoising and deblurring problems, IEEE Trans. Image Process., vol. 18, no. 11, pp , November [28] U. S. Kamilov, Parallel proximal methods for total variation minimization, in IEEE Int. Conf. Acoustics, Speech and Signal Process. (ICASSP 2016), Shanghai, China, March 19-25, 2015, arxiv: [cs.it]. [29] A. Taflove and S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method, 3rd ed. Artech House Publishers, [30] O. M. Bucci, L. Cardace, L. Crocco, and T. Isernia, Degree of nonlinearity and a new solution precudure in scalar two-dimensional inverse scattering problems, J. Opt. Soc. Am. A, vol. 18, no. 8, pp , August
A Recursive Born Approach to Nonlinear Inverse Scattering
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com A Recursive Born Approach to Nonlinear Inverse Scattering Kamilov, U.; Liu, D.; Mansour, H.; Boufounos, P.T. TR2016-057 June 2016 Abstract
More informationSEAGLE: Robust Computational Imaging under Multiple Scattering
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com SEAGLE: Robust Computational Imaging under Multiple Scattering Liu, H.-Y.; Liu, D.; Mansour, H.; Boufounos, P.T.; Waller, L.; Kamilov, U. TR217-81
More informationMinimizing Isotropic Total Variation without Subiterations
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Minimizing Isotropic Total Variation without Subiterations Kamilov, U. S. TR206-09 August 206 Abstract Total variation (TV) is one of the most
More informationarxiv: v1 [cs.cv] 5 May 2017
SEAGLE: Sparsity-Driven Image Reconstruction under Multiple Scattering Hsiou-Yuan Liu, Dehong Liu, Hassan Mansour, Petros T. Boufounos, Laura Waller, and Ulugbek S. Kamilov arxiv:175.4281v1 [cs.cv] 5 May
More informationLearning MMSE Optimal Thresholds for FISTA
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Learning MMSE Optimal Thresholds for FISTA Kamilov, U.; Mansour, H. TR2016-111 August 2016 Abstract Fast iterative shrinkage/thresholding algorithm
More informationSPARSE SIGNAL RESTORATION. 1. Introduction
SPARSE SIGNAL RESTORATION IVAN W. SELESNICK 1. Introduction These notes describe an approach for the restoration of degraded signals using sparsity. This approach, which has become quite popular, is useful
More informationIN THE reconstruction of material shapes and properties, we
1704 IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES, VOL. 46, NO. 11, NOVEMBER 1998 Nonlinear Inversion in TE Scattering Bert Jan Kooij and Peter M. van den Berg Abstract A method for reconstructing
More informationModified 2 gradient method and modified Born method for solving a two-dimensional inverse scattering problem
INSTITUTE OF PHYSICS PUBLISHING Inverse Problems 17 (2001) 1671 1688 INVERSE PROBLEMS PII: S0266-5611(01)24165-3 Modified 2 gradient method and modified Born method for solving a two-dimensional inverse
More informationEE 367 / CS 448I Computational Imaging and Display Notes: Image Deconvolution (lecture 6)
EE 367 / CS 448I Computational Imaging and Display Notes: Image Deconvolution (lecture 6) Gordon Wetzstein gordon.wetzstein@stanford.edu This document serves as a supplement to the material discussed in
More informationVersatile reconstruction framework for diffraction tomography with intensity measurements and multiple scattering
Vol. 6, No. 3 5 Feb 18 OPTICS EXPRESS 749 Versatile reconstruction framework for diffraction tomography with intensity measurements and multiple scattering T HANH - AN P HAM, 1,* E MMANUEL S OUBIES, 1
More informationNew Scaling Factors of 2-D Isotropic-Dispersion Finite Difference Time Domain (ID-FDTD) Algorithm for Lossy Media
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 56, NO., FEBRUARY 8 63 is shown by squares. Each point is computed by averaging reflection coefficient values calculated for each component of the regular
More informationHigh resolution tomographic diffraction microscopy
High resolution tomographic diffraction microscopy J. Girard,Y. Ruan 1, E. Mudry 1, F. Drsek G. Maire 1, P. Chaumet 1, H. Giovannini 1,K. Belkebir 1, A. Talneau 2, A. Sentenac 1 1 Institut Fresnel (Marseille)
More informationPublication II Wiley Periodicals. Reprinted by permission of John Wiley & Sons.
Publication II Ilkka Laakso and Tero Uusitupa. 2008. Alternative approach for modeling material interfaces in FDTD. Microwave and Optical Technology Letters, volume 50, number 5, pages 1211-1214. 2008
More informationAdaptive Primal Dual Optimization for Image Processing and Learning
Adaptive Primal Dual Optimization for Image Processing and Learning Tom Goldstein Rice University tag7@rice.edu Ernie Esser University of British Columbia eesser@eos.ubc.ca Richard Baraniuk Rice University
More informationAccelerated Dual Gradient-Based Methods for Total Variation Image Denoising/Deblurring Problems (and other Inverse Problems)
Accelerated Dual Gradient-Based Methods for Total Variation Image Denoising/Deblurring Problems (and other Inverse Problems) Donghwan Kim and Jeffrey A. Fessler EECS Department, University of Michigan
More informationTRACKING SOLUTIONS OF TIME VARYING LINEAR INVERSE PROBLEMS
TRACKING SOLUTIONS OF TIME VARYING LINEAR INVERSE PROBLEMS Martin Kleinsteuber and Simon Hawe Department of Electrical Engineering and Information Technology, Technische Universität München, München, Arcistraße
More informationA FISTA-like scheme to accelerate GISTA?
A FISTA-like scheme to accelerate GISTA? C. Cloquet 1, I. Loris 2, C. Verhoeven 2 and M. Defrise 1 1 Dept. of Nuclear Medicine, Vrije Universiteit Brussel 2 Dept. of Mathematics, Université Libre de Bruxelles
More informationNovel method for ultrashort laser pulse-width measurement based on the self-diffraction effect
Novel method for ultrashort laser pulse-width measurement based on the self-diffraction effect Peng Xi, Changhe Zhou, Enwen Dai, and Liren Liu Shanghai Institute of Optics and Fine Mechanics, Chinese Academy
More informationPOISSON noise, also known as photon noise, is a basic
IEEE SIGNAL PROCESSING LETTERS, VOL. N, NO. N, JUNE 2016 1 A fast and effective method for a Poisson denoising model with total variation Wei Wang and Chuanjiang He arxiv:1609.05035v1 [math.oc] 16 Sep
More informationA Localized Linearized ROF Model for Surface Denoising
1 2 3 4 A Localized Linearized ROF Model for Surface Denoising Shingyu Leung August 7, 2008 5 Abstract 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 1 Introduction CT/MRI scan becomes a very
More informationRouting of Deep-Subwavelength Optical Beams and Images without Reflection and Diffraction Using Infinitely Anisotropic Metamaterials
Peter B. Catrysse * and Shanhui Fan Routing of Deep-Subwavelength Optical Beams and Images without Reflection and Diffraction Using Infinitely Anisotropic Metamaterials Media that are described by extreme
More informationThe Theory of Diffraction Tomography
The Theory of Diffraction Tomography Paul Müller 1, Mirjam Schürmann, and Jochen Guck Biotechnology Center, Technische Universität Dresden, Dresden, Germany (Dated: October 10, 2016) Abstract arxiv:1507.00466v3
More informationCOLLOCATED SIBC-FDTD METHOD FOR COATED CONDUCTORS AT OBLIQUE INCIDENCE
Progress In Electromagnetics Research M, Vol. 3, 239 252, 213 COLLOCATED SIBC-FDTD METHOD FOR COATED CONDUCTORS AT OBLIQUE INCIDENCE Lijuan Shi 1, 3, Lixia Yang 2, *, Hui Ma 2, and Jianning Ding 3 1 School
More informationAnalysis of second-harmonic generation microscopy under refractive index mismatch
Vol 16 No 11, November 27 c 27 Chin. Phys. Soc. 19-1963/27/16(11/3285-5 Chinese Physics and IOP Publishing Ltd Analysis of second-harmonic generation microscopy under refractive index mismatch Wang Xiang-Hui(
More informationTHE total-field/scattered-field (TFSF) boundary, first proposed
454 IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS, VOL. 5, 2006 Analytic Field Propagation TFSF Boundary for FDTD Problems Involving Planar Interfaces: Lossy Material and Evanescent Fields Kakhkhor Abdijalilov
More information/$ IEEE
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: EXPRESS BRIEFS, VOL. 55, NO. 9, SEPTEMBER 2008 937 Analytical Stability Condition of the Latency Insertion Method for Nonuniform GLC Circuits Subramanian N.
More informationScalar electromagnetic integral equations
Scalar electromagnetic integral equations Uday K Khankhoje Abstract This brief note derives the two dimensional scalar electromagnetic integral equation starting from Maxwell s equations, and shows how
More informationModeling microlenses by use of vectorial field rays and diffraction integrals
Modeling microlenses by use of vectorial field rays and diffraction integrals Miguel A. Alvarez-Cabanillas, Fang Xu, and Yeshaiahu Fainman A nonparaxial vector-field method is used to describe the behavior
More informationDeep convolutional framelets:
Deep convolutional framelets: application to diffuse optical tomography : learning based approach for inverse scattering problems Jaejun Yoo NAVER Clova ML OUTLINE (bottom up!) I. INTRODUCTION II. EXPERIMENTS
More informationFast Angular Synchronization for Phase Retrieval via Incomplete Information
Fast Angular Synchronization for Phase Retrieval via Incomplete Information Aditya Viswanathan a and Mark Iwen b a Department of Mathematics, Michigan State University; b Department of Mathematics & Department
More informationThe Deep Ritz method: A deep learning-based numerical algorithm for solving variational problems
The Deep Ritz method: A deep learning-based numerical algorithm for solving variational problems Weinan E 1 and Bing Yu 2 arxiv:1710.00211v1 [cs.lg] 30 Sep 2017 1 The Beijing Institute of Big Data Research,
More informationMMSE Denoising of 2-D Signals Using Consistent Cycle Spinning Algorithm
Denoising of 2-D Signals Using Consistent Cycle Spinning Algorithm Bodduluri Asha, B. Leela kumari Abstract: It is well known that in a real world signals do not exist without noise, which may be negligible
More informationA discretized Newton flow for time varying linear inverse problems
A discretized Newton flow for time varying linear inverse problems Martin Kleinsteuber and Simon Hawe Department of Electrical Engineering and Information Technology, Technische Universität München Arcisstrasse
More informationPhotonic nanojet enhancement of backscattering of light by nanoparticles: a potential novel visible-light ultramicroscopy technique
Photonic nanojet enhancement of backscattering of light by nanoparticles: a potential novel visible-light ultramicroscopy technique Zhigang Chen and Allen Taflove Department of Electrical and Computer
More informationA far-field based T-matrix method for three dimensional acoustic scattering
ANZIAM J. 50 (CTAC2008) pp.c121 C136, 2008 C121 A far-field based T-matrix method for three dimensional acoustic scattering M. Ganesh 1 S. C. Hawkins 2 (Received 14 August 2008; revised 4 October 2008)
More informationAPPLICATION OF THE MAGNETIC FIELD INTEGRAL EQUATION TO DIFFRACTION AND REFLECTION BY A CONDUCTING SHEET
In: International Journal of Theoretical Physics, Group Theory... ISSN: 1525-4674 Volume 14, Issue 3 pp. 1 12 2011 Nova Science Publishers, Inc. APPLICATION OF THE MAGNETIC FIELD INTEGRAL EQUATION TO DIFFRACTION
More informationComparative Study of Restoration Algorithms ISTA and IISTA
Comparative Study of Restoration Algorithms ISTA and IISTA Kumaresh A K 1, Kother Mohideen S 2, Bremnavas I 3 567 Abstract our proposed work is to compare iterative shrinkage thresholding algorithm (ISTA)
More informationBinary Compressive Sensing via Analog. Fountain Coding
Binary Compressive Sensing via Analog 1 Fountain Coding Mahyar Shirvanimoghaddam, Member, IEEE, Yonghui Li, Senior Member, IEEE, Branka Vucetic, Fellow, IEEE, and Jinhong Yuan, Senior Member, IEEE, arxiv:1508.03401v1
More informationA Novel Single-Source Surface Integral Method to Compute Scattering from Dielectric Objects
SUBMITTED TO IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS ON NOVEMBER 18, 2016 1 A Novel Single-Source Surface Integral Method to Compute Scattering from Dielectric Objects Utkarsh R. Patel, Student
More informationANALYSIS OF DISPERSION RELATION OF PIECEWISE LINEAR RECURSIVE CONVOLUTION FDTD METHOD FOR SPACE-VARYING PLASMA
Progress In Electromagnetics Research Letters, Vol. 22, 83 93, 2011 ANALYSIS OF DISPERSION RELATION OF PIECEWISE LINEAR RECURSIVE CONVOLUTION FDTD METHOD FOR SPACE-VARYING PLASMA X. Ai Science and Technology
More informationAmbiguity of optical coherence tomography measurements due to rough surface scattering
Ambiguity of optical coherence tomography measurements due to rough surface scattering Y. Ashtamker, 1 V Freilikher, 1,* and J C Dainty 2 1 Department of Physics, Bar-Ilan University, Ramat Gan 52900,
More informationarxiv: v1 [physics.optics] 5 Mar 2012
Designing and using prior knowledge for phase retrieval Eliyahu Osherovich, Michael Zibulevsky, and Irad Yavneh arxiv:1203.0879v1 [physics.optics] 5 Mar 2012 Computer Science Department, Technion Israel
More informationConvex Hodge Decomposition of Image Flows
Convex Hodge Decomposition of Image Flows Jing Yuan 1, Gabriele Steidl 2, Christoph Schnörr 1 1 Image and Pattern Analysis Group, Heidelberg Collaboratory for Image Processing, University of Heidelberg,
More informationNumerical Analysis of Electromagnetic Fields in Multiscale Model
Commun. Theor. Phys. 63 (205) 505 509 Vol. 63, No. 4, April, 205 Numerical Analysis of Electromagnetic Fields in Multiscale Model MA Ji ( ), FANG Guang-You (ྠ), and JI Yi-Cai (Π) Key Laboratory of Electromagnetic
More information28th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
FINITE DIFFERENCE METHODS FOR ACOUSTIC AND ACOUSTO-GRAVITY WAVEFIELDS: APPLICATION TO LOW FREQUENCY INFRASOUND PROPAGATION Catherine de Groot-Hedlin University of California, San Diego Sponsored by Army
More informationBackscattering enhancement of light by nanoparticles positioned in localized optical intensity peaks
Backscattering enhancement of light by nanoparticles positioned in localized optical intensity peaks Zhigang Chen, Xu Li, Allen Taflove, and Vadim Backman We report what we believe to be a novel backscattering
More informationGeneralized optical theorem for reflection, transmission, and extinction of power for electromagnetic fields
PHYSICAL REVIEW E 71, 5661 5 Generalized optical theorem for reflection, transmission, and extinction of power for electromagnetic fields D. R. Lytle II Department of Electrical and Computer Engineering,
More information2 Regularized Image Reconstruction for Compressive Imaging and Beyond
EE 367 / CS 448I Computational Imaging and Display Notes: Compressive Imaging and Regularized Image Reconstruction (lecture ) Gordon Wetzstein gordon.wetzstein@stanford.edu This document serves as a supplement
More informationSparse Regularization via Convex Analysis
Sparse Regularization via Convex Analysis Ivan Selesnick Electrical and Computer Engineering Tandon School of Engineering New York University Brooklyn, New York, USA 29 / 66 Convex or non-convex: Which
More informationPAPER Fast Algorithm for Solving Matrix Equation in MoM Analysis of Large-Scale Array Antennas
2482 PAPER Fast Algorithm for Solving Matrix Equation in MoM Analysis of Large-Scale Array Antennas Qiang CHEN, Regular Member, Qiaowei YUAN, Nonmember, and Kunio SAWAYA, Regular Member SUMMARY A new iterative
More informationA memory gradient algorithm for l 2 -l 0 regularization with applications to image restoration
A memory gradient algorithm for l 2 -l 0 regularization with applications to image restoration E. Chouzenoux, A. Jezierska, J.-C. Pesquet and H. Talbot Université Paris-Est Lab. d Informatique Gaspard
More informationFourier Methods in Array Processing
Cambridge, Massachusetts! Fourier Methods in Array Processing Petros Boufounos MERL 2/8/23 Array Processing Problem Sources Sensors MERL 2/8/23 A number of sensors are sensing a scene. A number of sources
More informationNon-negative Quadratic Programming Total Variation Regularization for Poisson Vector-Valued Image Restoration
University of New Mexico UNM Digital Repository Electrical & Computer Engineering Technical Reports Engineering Publications 5-10-2011 Non-negative Quadratic Programming Total Variation Regularization
More informationGeneralized Analysis of Stability and Numerical Dispersion in the Discrete-Convolution FDTD Method
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 48, NO. 6, JUNE 2000 887 Generalized Analysis of Stability and Numerical Dispersion in the Discrete-Convolution FDTD Method William A. Beck, Member,
More informationMagnetoacoustic tomography with current injection
Article Progress of Projects Supported by NSFC Materials Science October 23 Vol.58 No.3: 36366 doi:.7/s434-3-5964-2 Magnetoacoustic tomography with current injection LIU GuoQiang *, HUANG Xin,2, XIA Hui
More informationAnthony J. Devaney Department of Electrical and Computer Engineering, Northeastern University, Boston, Massachusetts 02115
Time-reversal imaging with multiple signal classification considering multiple scattering between the targets Fred K. Gruber a) Department of Industrial Engineering and anagement Systems, University of
More informationRecovery of Sparse Signals from Noisy Measurements Using an l p -Regularized Least-Squares Algorithm
Recovery of Sparse Signals from Noisy Measurements Using an l p -Regularized Least-Squares Algorithm J. K. Pant, W.-S. Lu, and A. Antoniou University of Victoria August 25, 2011 Compressive Sensing 1 University
More informationTHE SCATTERING FROM AN ELLIPTIC CYLINDER IRRADIATED BY AN ELECTROMAGNETIC WAVE WITH ARBITRARY DIRECTION AND POLARIZATION
Progress In Electromagnetics Research Letters, Vol. 5, 137 149, 2008 THE SCATTERING FROM AN ELLIPTIC CYLINDER IRRADIATED BY AN ELECTROMAGNETIC WAVE WITH ARBITRARY DIRECTION AND POLARIZATION Y.-L. Li, M.-J.
More informationDenoising of NIRS Measured Biomedical Signals
Denoising of NIRS Measured Biomedical Signals Y V Rami reddy 1, Dr.D.VishnuVardhan 2 1 M. Tech, Dept of ECE, JNTUA College of Engineering, Pulivendula, A.P, India 2 Assistant Professor, Dept of ECE, JNTUA
More informationNIH Public Access Author Manuscript Proc Soc Photo Opt Instrum Eng. Author manuscript; available in PMC 2013 December 17.
NIH Public Access Author Manuscript Published in final edited form as: Proc Soc Photo Opt Instrum Eng. 2008 March 18; 6913:. doi:10.1117/12.769604. Tomographic Reconstruction of Band-limited Hermite Expansions
More informationA Time Domain Approach to Power Integrity for Printed Circuit Boards
A Time Domain Approach to Power Integrity for Printed Circuit Boards N. L. Mattey 1*, G. Edwards 2 and R. J. Hood 2 1 Electrical & Optical Systems Research Division, Faculty of Engineering, University
More informationTutorial on Tangent Propagation
Tutorial on Tangent Propagation Yichuan Tang Centre for Theoretical Neuroscience February 5, 2009 1 Introduction Tangent Propagation is the name of a learning technique of an artificial neural network
More informationFast Compressive Phase Retrieval
Fast Compressive Phase Retrieval Aditya Viswanathan Department of Mathematics Michigan State University East Lansing, Michigan 4884 Mar Iwen Department of Mathematics and Department of ECE Michigan State
More informationEnforcing Passivity for Admittance Matrices Approximated by Rational Functions
IEEE TRANSACTIONS ON POWER SYSTEMS, VOL. 16, NO. 1, FEBRUARY 2001 97 Enforcing Passivity for Admittance Matrices Approximated by Rational Functions Bjørn Gustavsen, Member, IEEE and Adam Semlyen, Life
More informationSparse Solutions of Systems of Equations and Sparse Modelling of Signals and Images
Sparse Solutions of Systems of Equations and Sparse Modelling of Signals and Images Alfredo Nava-Tudela ant@umd.edu John J. Benedetto Department of Mathematics jjb@umd.edu Abstract In this project we are
More informationRobust Sparse Recovery via Non-Convex Optimization
Robust Sparse Recovery via Non-Convex Optimization Laming Chen and Yuantao Gu Department of Electronic Engineering, Tsinghua University Homepage: http://gu.ee.tsinghua.edu.cn/ Email: gyt@tsinghua.edu.cn
More informationAn explicit time-domain finite-element method for room acoustics simulation
An explicit time-domain finite-element method for room acoustics simulation Takeshi OKUZONO 1 ; Toru OTSURU 2 ; Kimihiro SAKAGAMI 3 1 Kobe University, JAPAN 2 Oita University, JAPAN 3 Kobe University,
More informationAbsorption suppression in photonic crystals
PHYSICAL REVIEW B 77, 442 28 Absorption suppression in photonic crystals A. Figotin and I. Vitebskiy Department of Mathematics, University of California at Irvine, Irvine, California 92697, USA Received
More informationEFFICIENT SIMULATIONS OF PERIODIC STRUC- TURES WITH OBLIQUE INCIDENCE USING DIRECT SPECTRAL FDTD METHOD
Progress In Electromagnetics Research M, Vol. 17, 101 111, 2011 EFFICIENT SIMULATIONS OF PERIODIC STRUC- TURES WITH OBLIQUE INCIDENCE USING DIRECT SPECTRAL FDTD METHOD Y. J. Zhou, X. Y. Zhou, and T. J.
More informationCONTRAST SOURCE INVERSION METHOD: STATE OF ART
Progress In Electromagnetics Research, PIER 34, 189 218, 2001 CONTRAST SOURCE INVERSION METHOD: STATE OF ART P. M. van den Berg and A. Abubakar Centre for Technical Geosciences Delft University of Technology
More informationLecture 1: Numerical Issues from Inverse Problems (Parameter Estimation, Regularization Theory, and Parallel Algorithms)
Lecture 1: Numerical Issues from Inverse Problems (Parameter Estimation, Regularization Theory, and Parallel Algorithms) Youzuo Lin 1 Joint work with: Rosemary A. Renaut 2 Brendt Wohlberg 1 Hongbin Guo
More informationScaled gradient projection methods in image deblurring and denoising
Scaled gradient projection methods in image deblurring and denoising Mario Bertero 1 Patrizia Boccacci 1 Silvia Bonettini 2 Riccardo Zanella 3 Luca Zanni 3 1 Dipartmento di Matematica, Università di Genova
More informationRegularizing inverse problems using sparsity-based signal models
Regularizing inverse problems using sparsity-based signal models Jeffrey A. Fessler William L. Root Professor of EECS EECS Dept., BME Dept., Dept. of Radiology University of Michigan http://web.eecs.umich.edu/
More informationAn iterative hard thresholding estimator for low rank matrix recovery
An iterative hard thresholding estimator for low rank matrix recovery Alexandra Carpentier - based on a joint work with Arlene K.Y. Kim Statistical Laboratory, Department of Pure Mathematics and Mathematical
More informationFull Wave Analysis of RF Signal Attenuation in a Lossy Rough Surface Cave Using a High Order Time Domain Vector Finite Element Method
Progress In Electromagnetics Research Symposium 2006, Cambridge, USA, March 26-29 425 Full Wave Analysis of RF Signal Attenuation in a Lossy Rough Surface Cave Using a High Order Time Domain Vector Finite
More informationECE G: Special Topics in Signal Processing: Sparsity, Structure, and Inference
ECE 18-898G: Special Topics in Signal Processing: Sparsity, Structure, and Inference Sparse Recovery using L1 minimization - algorithms Yuejie Chi Department of Electrical and Computer Engineering Spring
More informationNear Optimal Adaptive Robust Beamforming
Near Optimal Adaptive Robust Beamforming The performance degradation in traditional adaptive beamformers can be attributed to the imprecise knowledge of the array steering vector and inaccurate estimation
More informationA Multilevel Proximal Algorithm for Large Scale Composite Convex Optimization
A Multilevel Proximal Algorithm for Large Scale Composite Convex Optimization Panos Parpas Department of Computing Imperial College London www.doc.ic.ac.uk/ pp500 p.parpas@imperial.ac.uk jointly with D.V.
More informationOptimisation using measured Green s function for improving spatial coherence in acoustic measurements
Ultrasonics 42 (2004) 205 212 www.elsevier.com/locate/ultras Optimisation using measured Green s function for improving spatial coherence in acoustic measurements Matthew Clark *, Steve D. Sharples, Mike
More informationNonparametric Rotational Motion Compensation Technique for High-Resolution ISAR Imaging via Golden Section Search
Progress In Electromagnetics Research M, Vol. 36, 67 76, 14 Nonparametric Rotational Motion Compensation Technique for High-Resolution ISAR Imaging via Golden Section Search Yang Liu *, Jiangwei Zou, Shiyou
More informationPublication I Institute of Physics Publishing (IOPP) Reprinted by permission of Institute of Physics Publishing.
Publication I Ilkka Laakso, Sami Ilvonen, and Tero Uusitupa. 7. Performance of convolutional PML absorbing boundary conditions in finite-difference time-domain SAR calculations. Physics in Medicine and
More informationarxiv: v1 [cs.it] 21 Feb 2013
q-ary Compressive Sensing arxiv:30.568v [cs.it] Feb 03 Youssef Mroueh,, Lorenzo Rosasco, CBCL, CSAIL, Massachusetts Institute of Technology LCSL, Istituto Italiano di Tecnologia and IIT@MIT lab, Istituto
More informationOptimal mass transport as a distance measure between images
Optimal mass transport as a distance measure between images Axel Ringh 1 1 Department of Mathematics, KTH Royal Institute of Technology, Stockholm, Sweden. 21 st of June 2018 INRIA, Sophia-Antipolis, France
More informationVariational methods for restoration of phase or orientation data
Variational methods for restoration of phase or orientation data Martin Storath joint works with Laurent Demaret, Michael Unser, Andreas Weinmann Image Analysis and Learning Group Universität Heidelberg
More informationarxiv: v1 [cs.lg] 20 Feb 2014
GROUP-SPARSE MATRI RECOVERY iangrong Zeng and Mário A. T. Figueiredo Instituto de Telecomunicações, Instituto Superior Técnico, Lisboa, Portugal ariv:1402.5077v1 [cs.lg] 20 Feb 2014 ABSTRACT We apply the
More informationNonlinear Compressive Sensing for Distorted Measurements and Application to Improving Efficiency of Power Amplifiers
Nonlinear Compressive Sensing for Distorted Measurements and Application to Improving Efficiency of Power Amplifiers Hsieh-Chung Chen, H.T. Kung, Marcus Comiter Harvard University Abstract Compressive
More informationPHASE RETRIEVAL OF SPARSE SIGNALS FROM MAGNITUDE INFORMATION. A Thesis MELTEM APAYDIN
PHASE RETRIEVAL OF SPARSE SIGNALS FROM MAGNITUDE INFORMATION A Thesis by MELTEM APAYDIN Submitted to the Office of Graduate and Professional Studies of Texas A&M University in partial fulfillment of the
More informationLarge-Scale L1-Related Minimization in Compressive Sensing and Beyond
Large-Scale L1-Related Minimization in Compressive Sensing and Beyond Yin Zhang Department of Computational and Applied Mathematics Rice University, Houston, Texas, U.S.A. Arizona State University March
More informationGPR prospecting in a layered medium via microwave tomography
ANNALS OF GEOPHYSICS, VOL. 46, N. 3, June 003 GPR prospecting in a layered medium via microwave tomography Lorenzo Crocco and Francesco Soldovieri Istituto per il Rilevamento Elettromagnetico dell Ambiente
More informationInterferometric model for phase analysis in fiber couplers
Interferometric model for phase analysis in fiber couplers Xiaojun Fang, Richard O. Claus, and Guy Indebetouw An interferometric model is proposed to estimate the phase differences in lossless, strongly
More informationarxiv: v1 [physics.optics] 30 Mar 2010
Analytical vectorial structure of non-paraxial four-petal Gaussian beams in the far field Xuewen Long a,b, Keqing Lu a, Yuhong Zhang a,b, Jianbang Guo a,b, and Kehao Li a,b a State Key Laboratory of Transient
More informationContinuous State MRF s
EE64 Digital Image Processing II: Purdue University VISE - December 4, Continuous State MRF s Topics to be covered: Quadratic functions Non-Convex functions Continuous MAP estimation Convex functions EE64
More informationLinear Inverse Problems
Linear Inverse Problems Ajinkya Kadu Utrecht University, The Netherlands February 26, 2018 Outline Introduction Least-squares Reconstruction Methods Examples Summary Introduction 2 What are inverse problems?
More informationFocusing Characteristic Analysis of Circular Fresnel Zone Plate Lens
Memoirs of the Faculty of Engineering, Okayama University, Vol.35, Nos.l&2, pp.53-6, March 200 Focusing Characteristic Analysis of Circular Fresnel Zone Plate Lens Tae Yong KIM! and Yukio KAGAWA 2 (Received
More informationQuadrature based Broyden-like method for systems of nonlinear equations
STATISTICS, OPTIMIZATION AND INFORMATION COMPUTING Stat., Optim. Inf. Comput., Vol. 6, March 2018, pp 130 138. Published online in International Academic Press (www.iapress.org) Quadrature based Broyden-like
More informationAbsorption-Amplification Response with or Without Spontaneously Generated Coherence in a Coherent Four-Level Atomic Medium
Commun. Theor. Phys. (Beijing, China) 42 (2004) pp. 425 430 c International Academic Publishers Vol. 42, No. 3, September 15, 2004 Absorption-Amplification Response with or Without Spontaneously Generated
More informationOWL to the rescue of LASSO
OWL to the rescue of LASSO IISc IBM day 2018 Joint Work R. Sankaran and Francis Bach AISTATS 17 Chiranjib Bhattacharyya Professor, Department of Computer Science and Automation Indian Institute of Science,
More informationLeveraging Machine Learning for High-Resolution Restoration of Satellite Imagery
Leveraging Machine Learning for High-Resolution Restoration of Satellite Imagery Daniel L. Pimentel-Alarcón, Ashish Tiwari Georgia State University, Atlanta, GA Douglas A. Hope Hope Scientific Renaissance
More informationElectromagnetic fields and waves
Electromagnetic fields and waves Maxwell s rainbow Outline Maxwell s equations Plane waves Pulses and group velocity Polarization of light Transmission and reflection at an interface Macroscopic Maxwell
More information2D X-Ray Tomographic Reconstruction From Few Projections
2D X-Ray Tomographic Reconstruction From Few Projections Application of Compressed Sensing Theory CEA-LID, Thalès, UJF 6 octobre 2009 Introduction Plan 1 Introduction 2 Overview of Compressed Sensing Theory
More information