Constrained simultaneous algebraic reconstruction technique (C-SART) —a new and simple algorithm applied to ionospheric tomography

A simple and relatively fast method (C-SART) is presented for tomographic reconstruction of the electron density distribution in the ionosphere using smooth fields. Since it does not use matrix algebra, it can be implemented in a low-level programming language, which speeds up applications significantly. Compared with traditional simultaneous algebraic reconstruction, this method facilitates both estimation of instrumental offsets and consideration of physical principles (expressed in the form of finite differences). Testing using a 2D scenario and an artificial data set showed that C-SART can be used for radio tomographic reconstruction of the electron density distribution in the ionosphere using data collected by global navigation satellite system ground receivers and low Earth orbiting satellites. Its convergence speed is significantly higher than that of classical SART, but it needs to be speeded up by a factor of 100 or more to enable it to be used for (near) real-time 3D tomographic reconstruction of the ionosphere.

[1]  Takashi Maruyama,et al.  Three‐dimensional ionospheric tomography using observation data of GPS ground receivers and ionosonde by neural network , 2005 .

[2]  G. Herman,et al.  Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and x-ray photography. , 1970, Journal of theoretical biology.

[3]  Stuart Crozier,et al.  Software developed for obtaining GPS‐derived total electron content values , 2007 .

[4]  Kefei Zhang,et al.  Ionospheric temporal and spatial variations during the 18 August 2003 storm over China , 2007 .

[5]  Giulio Ruffini,et al.  GPS tomography of the ionospheric electron content with a correlation functional , 1998, IEEE Trans. Geosci. Remote. Sens..

[6]  D. Wen.  Imaging the Ionospheric Electron Density Using a Combined Tomographic Algorithm , 2007 .

[7]  Ming Jiang,et al.  Convergence of the simultaneous algebraic reconstruction technique (SART) , 2003, Conference Record of Thirty-Fifth Asilomar Conference on Signals, Systems and Computers (Cat.No.01CH37256).

[8]  M. SIAMJ.,et al.  BLOCK-ITERATIVE ALGORITHMS WITH DIAGONALLY SCALED OBLIQUE PROJECTIONS FOR THE LINEAR FEASIBILITY PROBLEM , 2002 .

[9]  Jim R. Ray,et al.  Geodetic techniques for time and frequency comparisons using GPS phase and code measurements , 2005 .

[10]  HansenPer Christian.  The truncated SVD as a method for regularization , 1987 .

[11]  Pradip Kumar Bhuyan,et al.  Application of generalized singular value decomposition to ionospheric tomography , 2004 .

[12]  S. Ivansson.  Correction to "Seismic borehole tomography—Theory and computational methods" , 1986 .

[13]  Dieter Bilitza,et al.  International reference ionosphere , 1978 .

[14]  J. Ou,et al.  Monitoring the three-dimensional ionospheric electron density distribution using GPS observations over China , 2007 .

[15]  Rita Noumeir,et al.  Deriving the multiplicative algebraic reconstruction algorithm (MART) by the method of convex projection (POCS) , 1993, 1993 IEEE International Conference on Acoustics, Speech, and Signal Processing.

[16]  Karl-Rudolf Koch,et al.  Parameter estimation and hypothesis testing in linear models , 1988 .

[17]  Xingliang Huo,et al.  Three-dimensional ionospheric tomography by an improved algebraic reconstruction technique , 2007 .

[18]  Jaume Sanz,et al.  New approaches in global ionospheric determination using ground GPS data , 1999 .

[19]  Xiaoqing Pi,et al.  Data assimilation of ground GPS total electron content into a physics‐based ionospheric model by use of the Kalman filter , 2004 .

[20]  Bodo W. Reinisch,et al.  International Reference Ionosphere 2000 , 2001 .

[21]  Anthony Skjellum,et al.  Writing libraries in MPI , 1993, Proceedings of Scalable Parallel Libraries Conference.

[22]  E. Andreeva,et al.  Phase-difference radiotomography of the ionosphere , 1992 .

[23]  L. Dagum,et al.  OpenMP: an industry standard API for shared-memory programming , 1998 .

[24]  T. Başar,et al.  A New Approach to Linear Filtering and Prediction Problems , 2001 .

[25]  T. D. Raymund,et al.  Ionospheric tomography algorithms , 1994, Int. J. Imaging Syst. Technol..

[26]  G.E. Moore,et al.  Cramming More Components Onto Integrated Circuits , 1998, Proceedings of the IEEE.

[27]  Joachim Feltens,et al.  The activities of the Ionosphere Working Group of the International GPS Service (IGS) , 2003 .

[28]  A. Kak,et al.  Simultaneous Algebraic Reconstruction Technique (SART): A Superior Implementation of the Art Algorithm , 1984, Ultrasonic imaging.

[29]  P.S.J. Spencer,et al.  Ionospheric data assimilation methods for geodetic applications , 2004, PLANS 2004. Position Location and Navigation Symposium (IEEE Cat. No.04CH37556).

[30]  George L. Siscoe,et al.  The Solar–Terrestrial Environment , 1995 .

[31]  Farzad Kamalabadi,et al.  Localized three‐dimensional ionospheric tomography with GPS ground receiver measurements , 2007 .