Quantitative 3D-OCT motion correction with tilt and illumination correction, robust similarity measure and regularization

Abstract

Variability in illumination, signal quality, tilt and the amount of motion pose challenges for post-processing based 3D-OCT motion correction algorithms. We present an advanced 3D-OCT motion correction algorithm using image registration and orthogonal raster scan patterns aimed at addressing these challenges. An intensity similarity measure using the pseudo Huber norm and a regularization scheme based on a pseudo L0.5 norm are introduced. A two-stage registration approach was developed. In the first stage, only axial motion and axial tilt are coarsely corrected. This result is then used as the starting point for a second stage full optimization. In preprocessing, a bias field estimation based approach to correct illumination differences in the input volumes is employed. Quantitative evaluation was performed using a large set of data acquired from 73 healthy and glaucomatous eyes using SD-OCT systems. OCT volumes of both the optic nerve head and the macula region acquired with three independent orthogonal volume pairs for each location were used to assess reproducibility. The advanced motion correction algorithm using the techniques presented in this paper was compared to a basic algorithm corresponding to an earlier version and to performing no motion correction. Errors in segmentation-based measures such as layer positions, retinal and nerve fiber thickness, as well as the blood vessel pattern were evaluated. The quantitative results consistently show that reproducibility is improved considerably by using the advanced algorithm, which also significantly outperforms the basic algorithm. The mean of the mean absolute retinal thickness difference over all data was 9.9 um without motion correction, 7.1 um using the basic algorithm and 5.0 um using the advanced algorithm. Similarly, the blood vessel likelihood map error is reduced to 69% of the uncorrected error for the basic and to 47% of the uncorrected error for the advanced algorithm. These results demonstrate that our advanced motion correction algorithm has the potential to improve the reliability of quantitative measurements derived from 3D-OCT data substantially.

OCT angiography and sequential quantitative analysis of type 2 neovascularization after ranibizumab therapy

“…4 Motion correction was performed using registration of two orthogonally captured imaging volumes. 6,7 To delineate the plane to visualize the neovascular membrane, the automated segmentation lines were adjusted to the inner and outer margin of the lesion. En face images of the vasculature were generated by average intensity projection for the identified layer.…”

Projection-resolved optical coherence tomographic angiography

“…All 608 B-scans in each data cube were acquired in 2.9 seconds. Two volumetric raster scans, one x-fast scan and one y-fast scan, were acquired, registered [25,26], and merged into one 3D angiogram. Blood flow is detected using the AngioVue software, a commercial version of the split-spectrum amplitude-decorrelation angiography (SSADA) algorithm [1,27].…”

Automated choroidal neovascularization detection algorithm for optical coherence tomography angiography

“…The SSADA algorithm was applied to detect flow between the 2 consecutive B-scans at the same location [13,14]. The two scans were then registered and merged through an orthogonal registration algorithm [21].…”

“…The end result is that the vascular pattern from the superficial inner retina is replicated on the deeper outer retina. Secondly, while the effect of eye motion during the scan can be minimized by subtracting bulk motion noise [13,20] and using orthogonal registration [21], motion artifacts in the form of horizontal or vertical lines may remain. Finally, the intrinsic complexity of CNV also makes automated detection difficult.…”