OSA's Digital Library

Applied Optics

Applied Optics


  • Vol. 51, Iss. 11 — Apr. 10, 2012
  • pp: 1715–1723

Simulation research on improved regularized solution of the inverse problem in spectral extinction measurements

Janusz Mroczka and Damian Szczuczyński  »View Author Affiliations

Applied Optics, Vol. 51, Issue 11, pp. 1715-1723 (2012)

View Full Text Article

Enhanced HTML    Acrobat PDF (1072 KB)

Browse Journals / Lookup Meetings

Browse by Journal and Year


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools



We present further results of the simulation research on the constrained regularized least squares (CRLS) solution of the ill-conditioned inverse problem in spectral extinction (turbidimetric) measurements, which we originally presented in this journal [Appl. Opt. 49, 4591 (2010)]. The inverse problem consists of determining the particle size distribution (PSD) function of a particulate system on the basis of a measured extinction coefficient as a function of wavelength. In our previous paper, it was shown that under assumed conditions the problem can be formulated in terms of the discretized Fredholm integral equation of the first kind. The CRLS method incorporates two constraints, which the PSD sought will satisfy: nonnegativity of the PSD values and normalization of the PSD to unity when integrated over the whole range of particle size, into the regularized least squares (RLS) method. This leads to the quadratic programming problem, which is solved by means of the active set algorithm within the research. The simulation research that is the subject of the present paper is a continuation and extension of the research described in our previous paper. In the present research, the performance of the CRLS method variants is compared not only to the corresponding RLS method variants but also to other regularization techniques: the truncated generalized singular value decomposition and the filtered generalized singular value decomposition, as well as nonlinear iterative algorithms: The Twomey algorithm and the Twomey—Markowski algorithm. Moreover, two methods of selecting the optimum value of the regularization parameter are considered: The L-curve method and the generalized cross validation method. The results of our simulation research provide even stronger proof that the CRLS method performs considerably better with reconstruction of PSD than other inversing methods, in terms of better fidelity and smaller uncertainty.

© 2012 Optical Society of America

OCIS Codes
(290.3200) Scattering : Inverse scattering
(290.4020) Scattering : Mie theory
(290.5820) Scattering : Scattering measurements
(290.5850) Scattering : Scattering, particles
(290.7050) Scattering : Turbid media
(290.2558) Scattering : Forward scattering

ToC Category:

Original Manuscript: November 8, 2011
Manuscript Accepted: December 30, 2011
Published: April 5, 2012

Virtual Issues
Vol. 7, Iss. 6 Virtual Journal for Biomedical Optics

Janusz Mroczka and Damian Szczuczyński, "Simulation research on improved regularized solution of the inverse problem in spectral extinction measurements," Appl. Opt. 51, 1715-1723 (2012)

Sort:  Author  |  Year  |  Journal  |  Reset  


  1. J. Mroczka and D. Szczuczyński, “Improved regularized solution of the inverse problem in turbidimetric measurements,” Appl. Opt. 49, 4591–4603 (2010). [CrossRef]
  2. C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley-Interscience, 1983).
  3. A. R. Jones, “Light scattering for particle characterization,” Progr. Energy Combust. Sci. 25, 1–53 (1992).
  4. P. C. Hansen, “Regularization tools: A Matlab package for analysis and solution of discrete ill-posed problems,” Numerical Algorithms 6, 1–35 (1994).
  5. J. G. Crump and J. H. Seinfeld, “A new algorithm for inversion of aerosol size distribution data,” Aerosol Sci. Technol. 1, 15–34 (1981). [CrossRef]
  6. M. Kandlikar and G. Ramachandran, “Inverse methods for analysing aerosol spectrometer measurements: A Critical Review,” J. Aerosol Sci. 30, 413–437 (1999). [CrossRef]
  7. F. Stout and J. H. Kalivas, “Tikhonov regularization in standardized and general form for multivariate calibration with application towards removing unwanted spectral artifacts,” J. Chemom. 20, 22–33 (2006). [CrossRef]
  8. S. Twomey, “On the numerical solution of Fredholm integral equations of the first kind by the inversion of the linear system produced by quadrature,” J. Assoc. Comput. Mach. 10, 97–101 (1963). [CrossRef]
  9. D. L. Phillips, “A Technique for the numerical solution of certain integral equations of the first kind,” J. Assoc. Comput. Mach. 9, 84–97 (1962). [CrossRef]
  10. C. L. Lawson and R. J. Hanson, Solving Least Squares Problems (Prentice Hall, 1974).
  11. A. R. Roig and J. L. Alessandrini, “Particle size distributions from static light scattering with regularized Non-Negative Least Squares Constraints,” Part. Part. Syst. Charact. 23, 431–437 (2006). [CrossRef]
  12. G. H. Golub and U. von Matt, “Quadratically constrained least squares and quadratic problems,” Numer. Math. 59, 561–580 (1991). [CrossRef]
  13. P. E. Gill, W. Murray, M. A. Saunders, and M. H. Wright, “Procedures for Optimization Problems with a Mixture of Bounds and General Linear Constraints,” ACM Trans. Math. Softw. 10, 282–298 (1984). [CrossRef]
  14. P. C. Hansen, “Numerical tools for analysis and solution of Fredholm integral equations of the first kind,” Inverse Probl. 8, 849–872 (1992). [CrossRef]
  15. P. C. Hansen, “Regularization, GSVD and truncated GSVD,” BIT (Nord. Tidskr. Inf.-behandl.) 29, 491–504 (1989). [CrossRef]
  16. P. C. Hansen, “Relations between SVD and GSVD of discrete regularization problems in standard and general form,” Linear Algebra Appl. 141, 165–176 (1990). [CrossRef]
  17. P. C. Hansen, T. Sekii, and H. Shibahashi, “The modified truncated SVD method for regularization in general form,” SIAM J. Sci. Stat. Comput. 13, 1142–1150 (1992). [CrossRef]
  18. P. C. Hansen and D. P. O’Leary, “The use of the L-curve in the regularization of discrete ill-posed problems,” SIAM J. Sci. Comput. 14, 1487–1503 (1993). [CrossRef]
  19. P. J. McCarthy, “Direct analytic model of the L-curve for Tikhonov regularization parameter selection,” Inverse Probl. 19, 643–663 (2003). [CrossRef]
  20. G. H. Golub, M. Heath, and H. Wahba, “Generalized cross validation as a method for choosing a good ridge parameter,” Technometrics 21, 215–224 (1979). [CrossRef]
  21. S. Twomey, “Comparison of constrained linear inversion and an iterative non-linear algorithm applied to the indirect estimation of particle size distributions,” J. Comput. Phys. 18, 188–200 (1975). [CrossRef]
  22. W. Winklmayr, H. Wang, and W. John, “Adaptation of the Twomey Algorithm to the Inversion of Cascade Impactor Data,” Aerosol Sci. Technol. 13, 322–331 (1990). [CrossRef]
  23. G. R. Markowski, “Improving Twomey’s Algorithm for Inversion of Aerosol Measurement Data,” Aerosol Sci. Technol 7, 127–141 (1987). [CrossRef]
  24. Index of Refraction, Technical Note-007 (Duke Scientific Corporation, December 1, 1996).
  25. “Release on the refractive index of ordinary water substance as a function of wavelength, temperature and pressure” (International Association for the Properties of Water and Steam, 1997).

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited