Glaucoma, the second leading cause of blindness worldwide, is an optic neuropthy characterized by distinctive changes in the optic nerve head (ONH) and visual field. The detection of glaucomatous progression is one of the most important and most challenging aspects of primary open angle glaucoma (OAG) management. In this context, ocular imaging equipment is increasingly sophisticated, providing quantitative tools to measure structural changes in ONH topography, an essential element in determining whether the disease is getting worse. In particular, the Heidelberg Retina Tomograph (HRT), a confocal scanning laser technology, has been commonly used to detect glaucoma and monitor its progression.
In this paper, we present a new framework for detection of glaucomatous progression using HRT images. In contrast to previous works that do not integrate a priori knowledge available in the images, particularly the spatial pixel dependency in the change detection map, the Markov Random Field is proposed to handle such dependency. To our knowledge, this is the first application of the Variational Expectation Maximization (VEM) algorithm for inferring topographic ONH changes in the glaucoma progression detection framework. Diagnostic performance of the proposed framework is compared to recently proposed methods of progression detection.