Segmentation of macular fluorescein angiographies. A statistical approach

This paper is concerned with the use of Bayesian methods in the segmentation of macular #uorescein angiographies. Fluorescein angiography is used in ophthalmic practice to evaluate vascular retinopathies and choroidopathies: Sodium #uorescein is injected in the arm’s cubital vein of the patient and its distribution is observed along retinal vessels at certain times. In this task a previous and essential step is the segmentation of the image into its relevant components. In order to obtain this segmentation Bayesian methods can be used because a previous knowledge about the spatial structure of the scene to be segmented is available in this kind of images. The stochastic model assumed for the observed intensities is a simple model with a Gaussian noise process which is statiscally independent between pixels. The process of labels x is modelled as a Markov random "eld with a space-dependent external "eld expressing the anatomy of the ocular fundus and higher order interactions encouraging blood vessels to be thin and large. This procedure is applied to di!erent cases of diabetic retinophaty and vein occlusions. Two algorithms have been used to estimate x, simulated annealing and iterated conditional modes. In order to evaluate the accuracy of the estimation several error measures have been calculated. ( 2001 Pattern Recognition Society. Published by Elsevier Science Ltd. All rights reserved.

[1]  D. M. Titterington,et al.  Bayesian Image Restoration: An Application to Edge-Preserving Surface Recovery , 1993, IEEE Trans. Pattern Anal. Mach. Intell..

[2]  C. Geyer,et al.  Constrained Monte Carlo Maximum Likelihood for Dependent Data , 1992 .

[3]  Walter R. Gilks,et al.  MCMC in image analysis , 1995 .

[4]  B. Ripley,et al.  Pattern Recognition , 1968, Nature.

[5]  J. Besag Statistical Analysis of Non-Lattice Data , 1975 .

[6]  Robert Haining,et al.  Statistics for spatial data: by Noel Cressie, 1991, John Wiley & Sons, New York, 900 p., ISBN 0-471-84336-9, US $89.95 , 1993 .

[7]  Ying Sun,et al.  Recursive tracking of vascular networks in angiograms based on the detection-deletion scheme , 1993, IEEE Trans. Medical Imaging.

[8]  G. S. Watson,et al.  Goodness-of-fit tests on a circle. II , 1961 .

[9]  Stuart Geman,et al.  Statistical methods for tomographic image reconstruction , 1987 .

[10]  Leszek Wojnar,et al.  Image Analysis , 1998 .

[11]  J. Besag On the Statistical Analysis of Dirty Pictures , 1986 .

[12]  V. Johnson A Model for Segmentation and Analysis of Noisy Images , 1994 .

[13]  L. Baxter Random Fields on a Network: Modeling, Statistics, and Applications , 1996 .

[14]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[15]  Mike Rees,et al.  5. Statistics for Spatial Data , 1993 .

[16]  Sylvia Richardson,et al.  Bayesian mapping of disease , 1995 .

[17]  Arnoldo Frigessi,et al.  Stochastic models, statistical methods, and algorithms in image analysis : proceedings of the special year on image analysis held in Rome, Italy, 1990 , 1992 .

[18]  N. Karssemeijer,et al.  Application of Bayesian Methods to Segmentation in Medical Images , 1992 .

[19]  G. B. Smith,et al.  Preface to S. Geman and D. Geman, “Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images” , 1987 .

[20]  Xavier Guyon,et al.  Random fields on a network , 1995 .

[21]  P. Green,et al.  Global and local priors, and the location of lesions using gamma-camera imagery , 1991, Philosophical Transactions of the Royal Society of London. Series A: Physical and Engineering Sciences.

[22]  A. Sim,et al.  Bayesian detection of the fovea in eye fundus angiographies , 1998 .

[23]  A. Molli'e Bayesian mapping of disease , 1996 .

[24]  William H. Press,et al.  Numerical recipes in C , 2002 .

[25]  Amelia Simó,et al.  Bayesian detection of the fovea in eye fundus angiographies , 1999, Pattern Recognit. Lett..

[26]  K. Mardia Statistics of Directional Data , 1972 .

[27]  Michael H. Goldbaum,et al.  Image Understanding for Automated Retinal Diagnosis. , 1989 .

[28]  Rafael Molina,et al.  On the Hierarchical Bayesian Approach to Image Restoration: Applications to Astronomical Images , 1994, IEEE Trans. Pattern Anal. Mach. Intell..

[29]  C. Geyer On the Convergence of Monte Carlo Maximum Likelihood Calculations , 1994 .

[30]  T. Stephenson Image analysis , 1992, Nature.

[31]  I. Weir Fully Bayesian Reconstructions from Single-Photon Emission Computed Tomography Data , 1997 .

[32]  L. Younes Estimation and annealing for Gibbsian fields , 1988 .

[33]  Guillermo Ayala,et al.  A gray-level 2D feature detector using circular statistics , 1997, Pattern Recognit. Lett..

[34]  P. Green Bayesian reconstructions from emission tomography data using a modified EM algorithm. , 1990, IEEE transactions on medical imaging.

[35]  Sebastian Stier,et al.  Extraction of line properties based on direction fields , 1996, IEEE Trans. Medical Imaging.

[36]  Mário A. T. Figueiredo,et al.  Bayesian estimation of ventricular contours in angiographic images , 1992, IEEE Trans. Medical Imaging.

[37]  F. Zana,et al.  A region merging algorithm using mathematical morphology: application to macula detection , 1998 .

[38]  Peter Green,et al.  Markov chain Monte Carlo in Practice , 1996 .