Abstract
This paper deals with the modelling of high bit-depth images acquired by astronomical cameras using the discrete wavelet transform and the undecimated discrete wavelet transform for image representation. The probability density function (PDF) model parameters are estimated using the expectation-maximization (EM) algorithm and the method of moments. As proposed in this paper, the task of estimating the overall PDF model parameters can be simplified by so-called mirroring of the initial model which is estimated only for those wavelet coefficients that are greater than or equal to zero. In the case of the EM algorithm, this technique significantly reduces the computational cost of the model fitting algorithm. In our experiments, we achieved a reduction of more than 70%. In the case of the method of moments, this technique simplifies a system of moment equations. Three main PDF models are presented here: firstly, the mirrored mixture of a half-normal distribution and an exponential distribution, secondly, the mirrored mixture of two exponential distributions, and finally, the mirrored mixture of two half-normal distributions. Performance of these models is evaluated on three sets of astronomical images and also on artificial data using the Jeffrey divergence metric. Overall, the mirrored mixture of a half-normal and an exponential distribution overcomes the commonly used GLM (generalized Laplacian model) and also the other studied models.
Similar content being viewed by others
References
Baraniuk, R., Choi, H., Fernandes, F., et al.: Rice Wavelet Tools version 2.4. http://www.ece.rice.edu/dsp-/software/RWT (2002). Online; Accessed 2007
Bishop, C.: Pattern Recognition and Machine Learning. Springer, New York (2006)
Borran, M.J., Nowak, R.D.: Wavelet-based denoising using hidden markov models. In: IEEE International Conference on Acoustics, Speech, and Signal Processing. Proceedings (2001)
Buil, C.: CCD Astronomy. Willman-Bell, Inc. (1991)
Coleman, T., Li, Y.: On the convergence of interior-reflective newton methods for nonlinear minimization subject to bounds. Math. Program. 67(2), 189–224 (1994)
Coleman, T., Li, Y.: An interior, trust region approach for nonlinear minimization subject to bounds. SIAM J. Optim. 6(2), 418–445 (1996)
Crouse, M.S., et al.: Wavelet-based statistical signal processing using hidden Markov models. IEEE Trans. Signal Process. 46(4), 886–901 (1998)
DuBose, T.B., et al.: Statistical models of signal and noise and fundamental limits of segmentation accuracy in retinal optical coherence tomography. IEEE Trans. Med. Imaging 37(9), 1978–1988 (2018)
Forcina, A., Carbone, P.: Modelling dark current and hot pixels in imaging sensors. Statistical Modelling (2018 (Online First))
Holschneider, M., et al.: Wavelets, time-frequency methods and phase space, Real Time Algorithm for Signal Analysis with the Help of the Wavelet Transform. Springer, Berlin (1990)
Izenman, A.J.: Recent developments in nonparametric density estimation. J. Am. Stat. Assoc. 86(413), 205–224 (1991)
Jeffrey, H.: An invariant form for the prior probability in estimation problems. Proc. R. Soc. Lond. 186(1007), 453–461 (1946)
Jelinek, M., Castro-Tirado, A.J., Cunniffe, R.: A decade of GRB follow-up by BOOTES in Spain (2003–2013). Adv. Astron. 2016, 1–12 (2016)
Lyu, S.W., Simoncelli, E.P.: Modeling multiscale subbands of photographic images with fields of gaussian scale mixtures. IEEE Trans. Pattern Anal. Mach. Intell. 31(4), 693–706 (2009)
Mallat, S.G.: A theory for multiresolution signal decomposition: the wavelet representation. Trans. Pattern Anal. Mach. Intell. 11(7), 674–693 (1989)
McLachlan, G., Peel, D.: Finite Mixture Models. Wiley, London (2000)
Papoulis, A., Pillai, S.U.: Probability, Random Variables, and Stochastic Processes, 4th edn. McGraw Hill, Boston (2002)
Pinho, A.: On the impact of histogram sparseness on some lossless image compression techniques. In: Proceedings of International Conference on Image Processing (ICIP 2001), p. 442
Pizurica, A.: Image denoising using wavelets and spatial context modeling. Ph.D. thesis, Univ. Gent, Gent, Belgium (2002)
Sadreazami, H., Ahmad, M.O., Swamy, M.N.S.: Contourlet domain image modeling by using the alpha-stable family of distributions. In: IEEE International Symposium on Circuits and Systems (ISCAS), pp. 1288–1291 (2014)
Shandoosti, H.R., Hazavei, S.M.: Image denoising in dual contourlet domain using hidden Markov tree models. Digital Signal Process. 67, 17–29 (2017)
Simoncelli, E.P., Adelson, E.H.: Noise removal via bayesian wavelet coring. In: Proceedings of 3rd IEEE International Conference on Image Processing, pp. 379–382 (1996)
Sinno, Z., Caramanis, C., Bovik, A.C.: Towards a closed form second-order natural scene statistics model. IEEE Trans. Image Process. 27(7), 3194–3209 (2018)
Svihlik, J., Kukal, J., Fliegel, K., et al.: Estimation of non-Gaussian noise parameters in the wavelet domain using the moment-generating function. J. Electron. Imaging 21(2), 023025–1–023025–15 (2012)
Svihlik, J., et al.: Estimation of poisson noise in spatial domain. In: Proceedings of SPIE Applications of Digital Image Processing XL, vol. 10396, pp. 10396–10396-7 (2017)
Acknowledgements
This work has been supported by the grant GA17-05840S ”Multicriteria Optimization of Shift-Variant Imaging System Models” of the Czech Science Foundation.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Švihlík, J. Mirrored mixture PDF models for scientific image modelling. SIViP 16, 385–393 (2022). https://doi.org/10.1007/s11760-021-01944-z
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11760-021-01944-z