A quadratic programming approach for joint image reconstruction: mathematical and geophysical examples

Although a comparative analysis of multiple images of a physical target can be useful, a joint image reconstruction approach should provide better interpretative elements for multi-spectral images. We present a generalized image reconstruction algorithm for the simultaneous reconstruction of band-limited images based on the novel cross-gradients concept developed for geophysical imaging. The general problem is formulated as the search for those images that stay within their band limits, are geometrically similar and satisfy their respective data in a least-squares sense. A robust iterative quadratic programming scheme is used to minimize the resulting objective function. We apply the algorithm to synthetic data generated using linear mathematical functions and to comparative geophysical test data. The resulting images recovered the test targets and show improved structural semblance between the reconstructed images in comparison to the results from two other conventional approaches.

[1]  C. Vincent Tao,et al.  An Initial Study on Automatic Reconstruction of Ground DEMs from Airborne IfSAR DSMs , 2004 .

[2]  Benjamin M. W. Tsui,et al.  Principles of Medical Imaging , 1992 .

[3]  A. N. Tikhonov,et al.  Solutions of ill-posed problems , 1977 .

[4]  A. E. Hoerl,et al.  Ridge Regression: Applications to Nonorthogonal Problems , 1970 .

[5]  F. Morgan,et al.  Joint Seismic and Electrical Tomography , 1997 .

[6]  M. Meju,et al.  Characterization of heterogeneous near‐surface materials by joint 2D inversion of dc resistivity and seismic data , 2003 .

[7]  Jun Li,et al.  Radiography of soft tissue of the foot and ankle with diffraction enhanced imaging , 2004, Journal of the American Podiatric Medical Association.

[8]  M. A. Pérez-Flores,et al.  Imaging low-frequency and dc electromagnetic fields using a simple linear approximation , 2001 .

[9]  M. Meju,et al.  Joint two-dimensional DC resistivity and seismic travel time inversion with cross-gradients constraints , 2004 .

[10]  Michael A. Saunders,et al.  User''s guide for NPSOL (Ver-sion 4.0): A FORTRAN package for nonlinear programming , 1984 .

[11]  S. Webb.  The Physics of Medical Imaging , 1990 .

[12]  Alan G. Green,et al.  Discrete tomography and joint inversion for loosely connected or unconnected physical properties: application to crosshole seismic and georadar data sets , 2003 .

[13]  J. Vidale.  Finite‐difference calculation of traveltimes in three dimensions , 1990 .

[14]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[15]  P. Gill,et al.  Fortran package for constrained linear least-squares and convex quadratic programming. User's Guide for LSSOL (Version 1. 0) , 1986 .

[16]  Yun Zhang,et al.  Understanding image fusion , 2004 .

[17]  Enrique Gómez-Treviño,et al.  A versatile algorithm for joint 3D inversion of gravity and magnetic data , 2003 .

[18]  S. Twomey.  Introduction to the Mathematics of Inversion in Remote Sensing and Indirect Measurements , 1997 .

[19]  Enrique Gómez-Treviño.  Nonlinear integral equations for electromagnetic inverse problems , 1987 .

[20]  E. Haber,et al.  Joint inversion: a structural approach , 1997 .

[21]  Martin Rumpf,et al.  A Variational Approach to Nonrigid Morphological Image Registration , 2004, SIAM J. Appl. Math..

[22]  Donald W. Marquaridt.  Generalized Inverses, Ridge Regression, Biased Linear Estimation, and Nonlinear Estimation , 1970 .

[23]  Reinaldo J. Michelena,et al.  Tomographic traveltime inversion using natural pixels , 1991 .

[24]  R. Parker,et al.  Occam's inversion; a practical algorithm for generating smooth models from electromagnetic sounding data , 1987 .