Evaluation of the areal material distribution of paper from its optical transmission image

The goal of this study was to evaluate the areal mass distribution (defined as the X-ray transmission image) of paper from its optical transmission image. A Bayesian inversion framework was used in the related deconvolution process so as to combine indirect optical information with a priori knowledge about the type of paper imaged. The a priori knowledge was expressed in the form of an empirical Besov space prior distribution constructed in a computationally effective way using the wavelet transform. The estimation process took the form of a large-scale optimization problem, which was in turn solved using the gradient descent method of Barzilai and Borwein. It was demonstrated that optical transmission images can indeed be transformed so as to fairly closely resemble the ones that reflect the true areal distribution of mass. Furthermore, the Besov space prior was found to give better results than the classical Gaussian smoothness prior (here equivalent to Tikhonov regularization).

[1]  Matti Lassas. Eero Saksman,et al.  Discretization-invariant Bayesian inversion and Besov space priors , 2009, 0901.4220.

[2]  H. Triebel Function Spaces and Wavelets on Domains , 2008 .

[3]  Sos S. Agaian,et al.  Quantifying image similarity using measure of enhancement by entropy , 2007, SPIE Defense + Commercial Sensing.

[4]  R. Myllylä,et al.  On the estimate of effective optical parameters of close-packed fibrillar media , 2007 .

[5]  Kannappan Palaniappan,et al.  Adaptive Robust Structure Tensors for Orientation Estimation and Image Segmentation , 2005, ISVC.

[6]  E. Somersalo,et al.  Statistical and computational inverse problems , 2004 .

[7]  J. Hecksher-Sørensen,et al.  Optical Projection Tomography as a Tool for 3D Microscopy and Gene Expression Studies , 2002, Science.

[8]  Gene H. Golub,et al.  Tikhonov Regularization and Total Least Squares , 1999, SIAM J. Matrix Anal. Appl..

[9]  Stéphane Mallat,et al.  A Theory for Multiresolution Signal Decomposition: The Wavelet Representation , 1989, IEEE Trans. Pattern Anal. Mach. Intell..

[10]  J A Seibert,et al.  X-ray scatter removal by deconvolution. , 1988, Medical physics.

[11]  J F Williamson,et al.  Monte Carlo evaluation of kerma at a point for photon transport problems. , 1987, Medical physics.

[12]  J A Seibert,et al.  Characterization of the point spread function and modulation transfer function of scattered radiation using a digital imaging system. , 1986, Medical physics.

[13]  Stéphane Mallat,et al.  A Wavelet Tour of Signal Processing - The Sparse Way, 3rd Edition , 2008 .

[14]  S. D. Keller,et al.  Bipolar Spectral Analysis of Paper Formation Using Wavelet Transforms , 2007 .

[15]  Joost van de Weijer,et al.  Adaptive Structure Tensors and their Applications , 2006, Visualization and Processing of Tensor Fields.

[16]  J. Sharpe Optical projection tomography. , 2004, Annual review of biomedical engineering.

[17]  J.-P. Bernie,et al.  Local grammage distribution and formation of paper by light transmission image analysis , 1996 .

[18]  J. Borwein,et al.  Two-Point Step Size Gradient Methods , 1988 .