|
1.IntroductionOptical coherence tomography (OCT) has become an important biomedical tissue imaging technique over the last decade.1, 2 Especially, it has proven very promising within the area of ophthalmology, where it is currently being applied both qualitatively and quantitatively for a large variety of diseases. As a result, third-generation commercial systems (third-generation time-domain) already are available for this application area. Ophthalmic OCT systems are beneficial both for diagnosing retinal diseases and also in obtaining novel knowledge about the retinal structure and the development of different pathologies. Because OCT is based on low-coherence interferometry, the measured signals are affected by speckle noise,3 which can significantly lower the obtained signal-to-noise ratio (SNR) of the measured images. As the recorded speckle patterns will vary only in cases when the imaging geometry or setup is modified, either the position of the sample or the optical pathway should be allowed to vary between recorded scans to take advantage of repeated measurements in lowering the noise level. Frequency, spatial, angular compounding, and polarization diversity are possible ways3, 4 to decorrelate the speckle noise from one image to the next. If the speckle is fully uncorrelated between images, the SNR of the resulting average image is increased as the square root of the number of images used.5 For standard OCT commercial systems, such modifications of the measuring setup are not an option, and one is then left with the possibility of denoising the images using image post-processing techniques.6, 7 A simple approach is to apply a median filter or different kinds of low-pass filtering. More advanced schemes use adaptive wavelet schemes combined with nonlinear thresholds.8 Nonlinear anisotropic filtering has recently been suggested for enhancing retinal OCT images, and it was reported to facilitate segmentation of the retinal layers.9 One problem with these more advanced post-processing techniques is that there is a risk of generating artificial structures from the noise when applying the filters. We have found this to be a potential problem when doing image post-processing based on, e.g., coherent enhancing diffusion filters.10 Such problems are avoided if one instead implements an average procedure that calculates the statistical mean of the signal over different realizations of the stochastic process responsible for the noise. When applying OCT for in vivo imaging of the human retina, small movements of the head and saccadic eye movements will lead to small axial movements of the retina relative to the imaging geometry during image recording. This behavior will in reality lead to decorrelation of the speckle noise from one image to the next when recording a series of B-scans. The challenge, under the assumption that mesoscopic landmarks can be correlated, is to register these and shift and distort the consecutive tomograms for successful averaging. This paper illustrates how the disturbing noise including speckle can be significantly suppressed by combining a series of recorded OCT images, typically 8 to 12. The algorithm fuses a technique for aligning a series of one-dimensional (1-D) signals with a regularized dynamic programming algorithm that takes into account the spatial constraints/relationships of the two-dimensional images (i.e., using the fact that the estimated registration shift of a given A-scan should be correlated with the expected registration shifts of the neighboring A-scans). The resulting improvement facilitates quantification and clinical interpretation and has already proven most helpful for a number of pathologic applications. Initial reports on contrast enhancement of ophthalmic OCT images from combination of B-scans (using a less-rigorous and a less-robust approach than the methodology presented here) was reported in Ref. 11. In Sec. 2, the image registering problem is formulated. Next, in Secs. 3, 4, 5, 6, the algorithmic approach based on regularized dynamic programming is outlined and illustrated for a specific case. In Sec. 7, the resulting improvement in image quality is both illustrated visually and quantified using a suitable contrast to noise measure. The clinical value of the image enhancement is illustrated from a number of cases in Sec. 8. The paper concludes with a discussion in Sec. 9. OCT is today a standard imaging technique in ophthalmological clinics focusing on retinal diseases. Accordingly, the results of applying the suggested method were based on a retrospective analysis of digital image information and did not require institutional review. 2.Setup and Problem FormulationThe OCT data was recorded at Glostrup Hospital in Copenhagen using the Stratus/OCT3 instrument (Humphrey-Zeiss). The Stratus system provides as a standard 512 A-scans per B-scan, independent of the length of the physical scan length. The optical transverse resolution is , and the axial resolution is approximately . Accordingly, a standard -long scan has a lateral pixel spacing below the transverse resolution of the optical system. Figure 1a shows a series of successive retinal B-scans corresponding to the ideal scan path illustrated in the autofluorescence fundus photograph [Fig. 1b]. It is clearly observed how small changes in eye position relative to the imaging device result in axial displacements of the retinal structure. Our goal is to enhance the SNR by registering such a series of images in order to calculate an average B-scan image. In order to fuse the individual B-scans into a single average image, we need to register the images mutually. Corresponding A-scans (represented by image columns) must be registered vertically to deal with the axial displacements. As we want to ensure a noise-robust registration procedure, we will use a so-called regularized algorithm, which takes the spatial continuity of the retinal structure into account. In addition to the axial displacements, the inevitable eye movements will cause the sampled scanning points to scatter around the desired scan path. These movements include micro-saccade and tremor, which goes up to of rapid involuntary random motion in transverse and torsional directions. Large transverse shifts are more likely to happen if the patient has difficulties in fixating on the spot of the guiding beam used during recording of a B-scan. The transverse uncertainty naturally relates to the severity of the pathology of the retina being examined. Part of the transverse displacements can be dealt with by combining the vertical registrations (corresponding to axial registration) with a possible rigid horizontal registration of one image to another, i.e., all rows in an image must be shifted the same amount relative to one of the other images to maximize an adequate similarity measure. Even though the horizontal registration can significantly reduce the effect of transverse distortions, the fact that the loci of the A-scans scatter around the ideal scan path will naturally lead to some loss of transverse resolution in the combined average image compared to that of a singe B-scan. In Sec. 6, we estimate the transverse accuracy that can typically be obtained when producing an average image with increased SNR. The process of registering two OCT images will be illustrated with the two images shown in Figs. 2a and 2b . 3.Alignment ProcedureWhen axially registering the two B-scans in Figs. 2a and 2b, the objective is to find the vertical shifts between corresponding columns that will align them optimally according to a similarity measure. Here we apply the cross-correlation between shifted versions of image columns to measure how well they align. Because of the speckle noise (as well as other noise contributants), the shift leading to the maximum correlation does not necessarily specify the optimum position for an ideal alignment. We can, however, take advantage of the fact that axial shifts across the image in most cases is known to change smoothly and continuously, simply because the head/eye movements that have caused these shifts describe a continuous and smooth motion. (If a saccade or blink completely obscures the continuity of a B-scan, it is simply discarded by the user; see Sec. 9). The implication is that the obtained values of neigboring translation shifts (obtained by mutually registering two B-scans) must describe a continous and smooth curve. Such a constraint can be implemented using so-called regularization in combination with a dynamic programming algorithm.12 In our case, this is analogous to formulating a constrained shortest-path strategy for navigating through an energy landscape specifying the problem at hand. The regularization then corresponds to applying a penalty term, adding extra cost to the traversed path through the landscape. The penalty will be dependent on the curvature and smoothness of the path. More precisely, we apply a regularized shortest path (or regularized minimum cost) algorithm to an image map or energy landscape formed by calculating the negative cross-correlation for the possible vertical shifts ( -variable) between corresponding columns ( -variable) in the two considered B-scan images and (each with rows and columns): where the raw cross-correlation values between the A-scan given by column of image and the A-scan given by column of image is defined asand for negative The relative shift belongs to the interval , where is the number of rows in each image. In order to take the continuity and smoothness constraint into account, we will now estimate the column shifts by calculating the shortest directed path across the image by the use of regularized dynamic programming.4.Regularized Dynamic ProgrammingDynamic programming describes a methodology where optimal solutions of subproblems can be used to find the optimal solutions of the overall problem. A well-known application is to determine the shortest path or accumulated cost across a map/image. Its general principles were first introduced by Bellman.13 Dynamic programming is known to be computationally efficient and to find the optimal path, but in its original version, it does not allow for any smoothness constraints on the generated path across the image. An extension was introduced in Ref. 14, which applied a penalty depending on the shape of the contour. To maintain optimality, it was an iterative process. The method was adapted by Buckley in Ref. 12, where a method for finding smooth shortest paths across images was introduced. The authors of this paper refer to the method as “regularized shortest-path extraction.” For every pixel in the image with rows and columns, let be the gray/intensity value of the pixel located at . A path of order from the left side of the image to a point on the right side of the image is a set of pixels: The points have ’th-order connectivity, i.e., the change in row index from one point to the next on the path is less than or equal to . The length or cost of a path is defined as the sum of the pixel values in the pathThe shortest path going from one side of the image to the other is the path that minimizes Eq. 5.In regularized dynamic programming, a penalty that is proportional to a measure of the roughness of the path is added to the length of the path. The roughness of a path is defined as:12 This leads to a new definition of the length or regularized cost of path :where is a regularization constant. The optimal path will correspond to the shortest path length under the constraint of limited roughness. One can now define as the length of the shortest path from the left border to the point via , for . The initialization equation for the dynamic programming algorithm iswhich holds for and . (If the pixels addressed are outside the image map, the cost at these pixels is set to infinity.)A path from the left border to via must pass through for some . This gives a roughness contribution of . The recursion formula for the dynamic programming routine is as follows { specifies the minimum value of obtained for the given range of }: from which it is possible to calculate all values of . The value that leads to a minimum is stored for later backtracking in the following variable ( specifies the value of the argument that minimizes for ):When backtracking, the last point and the direction it came from is found asThe last point on the optimal path is, . The second to last point is also known: , where . The rest of the points on the optimal path can be calculated as , with . The computational cost of the regularized version of dynamic programming is .Due to the discrete spatial nature of the digital image, the possible slope values calculated along the directed graph are rather limited, and this in turn affects the smoothness evaluation. In order to increase the number of quantization values of the slope, it is often necessary to perform a vertical subsampling of the map image.12 As this trick is fairly expensive computationally, one can alternatively reduce the horizontal resolution when forming the map image, as also described in Ref. 12. By applying the shortest-path approach to mutually align the images shown in Figs. 2a and 2b, one obtains the aligned images shown in Figs. 2c and 2d. Taking the average of the two aligned images gives the image shown in Fig. 3 . A satisfactory vertical registering has been obtained, but from inspecting the two sides of the macular hole, one can observe that a horizontal displacement exists between the two combined images, and this needs to be compensated. A likely cause for the horizontal misalignment is the fact that patients with macular holes have their central vision impaired. This in turn causes difficulties for the patients to fixate, which increase the likelihood for larger deviations from the ideal scan path. 5.Horizontal RegistrationWhen we perform the vertical registration between two images, we neglect a possible horizontal mismatch and consider the vertical alignment problem as an isolated one. The reason this is possible is that the retinal images have an inherent horizontal structure, corresponding to a high correlation between most neighboring columns in a B-scan. It is therefore also important to perform a horizontal alignment of all images before we carry out the horizontal registration. The horizontal shift is simply estimated by cross-correlating the two images that have initially been vertically aligned using a candidate series of likely horizontal shifts. The average of the two considered images after vertical and horizontal registration is shown in Fig. 4 . The image depicted in Fig. 2d is shifted to the left corresponding to approximately . This horizontal alignment makes a significant improvement to the resultant image. The edges of the macular hole in the combined image are much sharper, and the shadow from a superficial blood vessel in the left part of the image is now clearly visible. Although the retinal layers have been succesfully registered, we can also observe an artifact/limitaton of the method: the vitreous membrane now has a double line on the left side and the operculum image shows a clear double/ghost image. 6.Aligning Multiple B-ScansThe preceding sections illustrated the procedure for registering two OCT images. In order to reduce the noise significantly, we need to align several images. In this section, the procedure is extended to any number of images by basically combining a method suggested for aligning a series of one-dimensional (1-D) signals15 with the preceding procedure. The entire method is illustrated in the flow chart in Fig. 5 . The first image of our sequence is initially used as a template image to which the rest of the images are being registered vertically. As noted earlier, we need to perform a horizontal alignment of all images before we carry out the horizontal registration. This is performed at the end of the first vertical registration (see Fig. 5). The first time the horizontal registration is carried out, we again use the first image as the template. In the following iteration cycles of the registration procedure (both vertical and horizontal), we use the current average image as the template image and register all individual images to this. Two iterations are in general sufficient for the iteration process to converge (toward a stable average image). In order to estimate and illustrate the transverse resolution of a combined average image, we again consider the scan path from Fig. 1, including retinal vessels. Blood vessels give rise to an intermediate signal level compared to the surrounding retinal tissue, but in many cases, the most outstanding feature is a shadow observed in underlying layers, caused by the light absorbance of the blood. This can be observed in Fig. 6 , showing an average image of 15 registered B-scans. The width of retinal veins can be estimated to be approximately , whereas the width of retinal arteries is around near the optic nerve.16, 17, 18 A rough estimate of the vessel diameter of the arteries visible on the left and right side of the average image (full drawn arrows in Fig. 6) is around as the vessels are located peripheral to the optic nerve. Below the superior vessel in Fig. 6 is a smaller arteriole (dashed arrow), the size being approximately . On the average image, this vessel (or its shadow) is just being resolved, indicating that the lateral resolution for this average image is close to . The external limiting membrane is also visible in the average image, indicating a high-quality picture with an axial resolution as claimed by the manufacturer, i.e., below . 7.ResultsThe registration procedure outlined earlier for combining a series of scans has been tested and used on a large number of patients, as discussed more thoroughly in Sec. 8. An example is illustrated in Fig. 7 , where 14 images corresponding to the example case of Fig. 2 have been combined. To quantify the improvement in SNRs, we have applied19 an adequate contrast-to-noise ratio to two neighboring image areas and with different mean intensities and , respectively, and with corresponding standard deviations and : From the definition, we see that the ratio can be considered as the square of the difference between the SNRs in the two considered regions of interest. If we were combining images with totally independent noise contributions, we would expect to gain a factor of in the ratio. Table 1 lists the improvement in as a function of the number of images being aligned for two series of OCT measurements. The numbers in the rightmost column are obtained for the image sequence used to form the average image in Fig. 8 . One observes that a factor of around 3 to 4 can be gained in the ratio when combining around 10 images. Ideally, we would gain a factor of 10 if the noise contributions had been totally independent. The improvement of using more than 12 images is often marginal, as seen from the listed values.8.Clinical ApplicationsThe registration procedure has been tested and used on both healthy subjects and patients with a wide range of pathologies and has shown very encouraging results. Fig. 8 (below) shows the enhancement obtained for a healthy subject and the succeding examples focus on different pathologies including quantitative analysis of macular hole patients. The abbreviations used throughout the case decriptions are: RNFL—retinal nerve fiber layer; GCL—ganglion cell layer; IPL—inner plexiform layer; INL—inner nuclear layer; OPL—outer plexiform layer; PL—photoreceptor layer (nuclei); EXL—external limiting membrane; IS/OS—inner segment/outer segment junction; OS—outer segments; and RPE—retinal pigment epithelium layer. 8.1.Healthy SubjectIn Fig. 8, a single B-scan from the Stratus system is compared to an average image of 14 B-scans. The single scan is well-centered on the fovea, as can be seen from the red hyperreflective spot in the foveal depression (the umbo of the fovea). In the average image, the well-known retinal layers are clearly defined and also the extreme thin line from the external limiting membrane is visible above the IS/OS junction. In principle, this membrane is hard to resolve with the resolution of the Stratus OCT system, but as it is “floating” as a single layer above the IS/OS layer, the increased SNR makes it visible at a lower resolution. The highly reflecting junction line representing the IS/OS junction of the photoreceptors is less prominent just in the foveal center, above the low-reflective layer of photoreceptor outer segments. 8.2.Diabetic Macular EdemaApproximately 20% of diabetic patients develop macular edema with time, ranging from subtle increased thickness to large, diffuse edema including all retinal layers.20 Incipient edema is hard to detect clinically, and OCT has found widespread use for this purpose both for screening and for diagnosis. The leakage originates from the retinal capillaries in the inner retina and diffuses toward the outer PL layer.21, 22 Retinal capillaries are hard to visualize with OCT due to low reflectance, although they do give a shadow beneath. Figure 9 shows a fundus image and a combined OCT image for a patient with moderate edema in the foveal center. In these images, the vessel trunks from the superficial retinal vessels are seen dipping down in the superficial capillary plexus and the deeper plexus at the level of the IPL/INL. When the horizontal alignment is included in the average image, the vessels are clearly visible, and as mentioned, the fluid accumulation is found away from the capillaries, in the photoreceptor layer. The mechanism behind this redistribution is not known in detail, and until now, no method has been available for evaluation of the relation between the morphological appearance of the edema and the effect of treatment.23, 24 In later stages, with diffuse edema, the edema is also found in the inner retinal layers (see Fig. 10 ). Cystic spaces are seen in the fovea, in the layer of photoreceptor nuclei and the overlying Henles fibers. Cystic spaces are also seen in the surrounding perifoveal region confined to the INL layer, and although the perifoveal cysts are smaller than the central foveal cysts, their presence is probably an indication of a less-favorable visual prognosis compared to patients without cystic spaces in the perifovea. Compared to the average image, the layer of cystic spaces in the INL layer in the single B-scan is less recognizable. On the other hand, the vertical boundaries of the cysts are slightly blurred in the average image due to the decrease in transverse resolution. But most noticeable, the presence of the EXL membrane and the junction of IS/OS segments are not visible in the single B-scan. If the outer retinal layers are absent, the visual prognosis is also unfavorable. In severe cases as this one, a serious detachment is seen in the fovea, just above the RPE. Macular holes can be closed surgically with a large improvement in visual acuity, the prognosis being related to the size of the hole and probably also to other morphological details such as the degree of cystic formation at the edge of the hole. OCT is of major interest here. The enhancement of the OCT images obtained by the registration procedure has found an important clinical applicability in the evaluation of retinal conditions after surgery for idiopathic macular hole where central vision often remains compromised25 (see Fig. 11 ). Since restoration of visual function after macular hole closure clearly depends on a remaining PL layer rather than the clinical observable macular surface configuration, visualization of outer retinal layers is important. Especially, visualization of the IS/OS integrity is an important indicator for photoreceptor survival. Since the denoised average OCT images enhance the visibility of the signals from the PL and RPE layers, quantitative analysis of central photoreceptor defects26 that may be present after macular hole surgery is facilitated. A clinical study is ongoing to evaluate the correlation between visual acuity and photoreceptor defects being quantified from enhanced OCT images. Table 1Contrast-to-noise ratio (CNR) as a function of the number of images used to produce the average image. Here, shown for two series of registrations.
9.DiscussionThis paper has demonstrated a novel image registration procedure applicable for retinal OCT B-scans obtained by standard OCT systems used for ophthalmology. The method combines a regularized shortest-path algorithm with an iterative scheme for aligning multiple signals. From the registered images, an average image can be calculated. The resulting suppression of noise (including speckle) produces images of a quality that otherwise can be obtained only using ultra-high-resolution OCT laboratory equipment. The enhanced contrast in the combined images facilitates image diagnostics and quantization of pathologies. In the vast majority of clinical patients, it is possible and feasible to record 12 to 20 images and then manually discard grossly deviating scans. The scans are typically obtained within . Typically around ten scans are left for the alignment procedure. In patients with a substantial loss of vision, this may not be accomplished. Some of these patients also have high levels of retinal scatter, impairing the subjective sorting part where images have to be selected or deselected, and then the combination technique is less powerful. An assumption of the implemented alignment procedure is that the macroscopic anatomy of the retina does not change significantly with transverse shifts over 30 to 40 microns (due to movement of consecutive tomograms). Naturally a certain loss of resolution cannot be avoided, as also can be observed from the average images. On the other hand, the improvement in SNR is in general so high that more retinal details actually are visible in the average image compared with the single recordings. The present implementation of the registration algorithm has been implemented in and the registration (both horizontal and vertical) of 10 images takes around on a Pentium 4. Optimization of the code with respect to the processing time has not been performed. For macular edema, photocoagulation treatment is the proven treatment, although it is well known that the treatment is not always effective, and new treatments are now in ongoing clinical studies. A possible reason for the slow development in this field is our limited understanding of the stages of severity, leading to poor results in clinical studies and in the evaluation of the individual patient. The methods presented here give improved possibilities for classification based on the retinal morphology by itself. The enhanced OCT images are also useful when visual acuity is below that expected with no obvious cause. Edema is easy visualized with standard software, but intraretinal defects like atrophy of the outer photoreceptor layers and/or increased backscatter from the outer photoreceptor nuclear layers is clearly more visible with higher SNR, as in the enhanced images presented here. Although ultra-high-resolution OCT systems give a better resolution,27 these systems will not be available for most ophthalmological departments in the near future. In clinical setups, OCT systems are now used extensively, and for difficult cases and research, the presented contrast enhancement obtained with images from commercial instruments is highy useful. Recently, an eye-tracking system has been combined with the [Humphrey-Zeiss] OCT II commercial system (second-generation time-domain) in order to significantly reduce the registration error relative to the planned transverse scanning path.28 With such a system, one needs only to perform the axial registration when combining a series of recordings to form the average image. Although the average/registration procedure is likely to have most impact on “slow” speed (time-domain OCT) and low-resolution OCT systems, the technique described here can still be beneficial for enhancing the contrast in images obtained using the new generation of Fourier domain OCT systems. AcknowledgmentsThe authors would like to thank the Danish Eye Research Foundation, the Danish Eye Health Society, the Danish Medical Research Council, and the Danish Research Academy (framework program BIOPHOT) for financial support. ReferencesHandbook of Optical Coherence Tomography, Marcel Dekker, New York
(2002). Google Scholar
M. R. Hee,
J. A. Izatt,
E. A. Swanson,
D. Huang,
J. S. Schuman,
C. P. Lin,
C. A. Puliafito, and
J. G. Fujimoto,
“Optical coherence tomography of the human retina,”
Arch. Ophthalmol. (Chicago), 113 325
–332
(1995). 0003-9950 Google Scholar
J. M. Schmitt,
“Speckle in optical coherence tomography,”
J. Biomed. Opt., 4 95
–105
(1999). https://doi.org/10.1117/1.429925 1083-3668 Google Scholar
A. E. Desjardins,
B. J. Vakoc,
G. J. Tearney, and
B. E. Bouma,
“Speckle reduction in OCT using massively-parallel detection and frequency-domain ranging,”
Opt. Express, 14 4736
–4745
(2006). https://doi.org/10.1364/OE.14.004736 1094-4087 Google Scholar
M. Bashkansky and
J. Reintjes,
“Statistics and reduction of speckle in optical coherence tomography,”
Opt. Lett., 25 545
–547
(2000). https://doi.org/10.1364/OL.25.000545 0146-9592 Google Scholar
D. L. Marks,
T. S. Ralston, and
S. A. Boppart,
“Speckle reduction by I-divergence regularization in optical coherence tomography,”
J. Opt. Soc. Am. A, 22 2366
–2371
(2005). https://doi.org/10.1364/JOSAA.22.002366 0740-3232 Google Scholar
S. Paes,
S. Y. Ryu,
J. Na,
E. Choi,
B. H. Lee, and
I. K. Hong,
“Advantages of adaptive speckle filtering prior to alication of iterative deconvolution methods for optical coherent tomography imaging,”
Opt. Quantum Electron., 37 1225
–1238
(2005). https://doi.org/10.1007/s11082-005-4194-5 0306-8919 Google Scholar
D. C. Adler,
T. H. Ko, and
J. G. Fujimoto,
“Speckle reduction in optical coherence tomography images by use of a spatially adaptive wavelet filter,”
Opt. Lett., 29 2878
–2880
(2004). https://doi.org/10.1364/OL.29.002878 0146-9592 Google Scholar
D. C. Fernandez,
“Delineating fluid-filled region boundaries in optical coherence tomography images of the retina,”
IEEE Trans. Med. Imaging, 24 929
–945
(2005). https://doi.org/10.1109/TMI.2005.848655 0278-0062 Google Scholar
D. C. Fernández,
H. M. Salinas, and
C. A. Puliafito,
“Automated detection of retinal layer structures on optical coherence tomography images,”
Opt. Express, 13 10200
–10216
(2005). https://doi.org/10.1364/OPEX.13.010200 1094-4087 Google Scholar
B. Sander,
M. Larsen,
L. Thrane,
J. L. Hougaard, and
T. M. Jørgensen,
“Enhanced optical coherence tomography imaging by multiple scan averaging,”
Br. J. Ophthamol., 89 207
–212
(2005). https://doi.org/10.1136/bjo.2004.045989 0007-1161 Google Scholar
M. Buckley and
J. Yang,
“Regularized shortest-path extraction,”
Pattern Recogn. Lett., 18 621
–629
(1997). https://doi.org/10.1016/S0167-8655(97)00076-7 0167-8655 Google Scholar
R. E. Bellman, Dynamic Programming, Princeton University Press, NJ
(1957). Google Scholar
A. A. Amini,
T. E. Weymouth, and
R. C. Jain,
“Using dynamic programming for solving variational problems in vision,”
IEEE Trans. Pattern Anal. Mach. Intell., 12 855
–867
(1990). https://doi.org/10.1109/34.57681 0162-8828 Google Scholar
K. J. Coakley and
P. Hale,
“Alignment of noisy signals,”
IEEE Trans. Instrum. Meas., 50 141
–149
(2001). https://doi.org/10.1109/19.903892 0018-9456 Google Scholar
Early Treatment Diabetic Retinopathy Study Research Group,
“Grading diabetic retinopathy from stereoscopic color fundus photographs—an extension of the modified Airlie House classification, ETDRS Report No. 10,”
Ophthalmology, 98 786
–806
(1991). 0161-6420 Google Scholar
G. Weigert,
C. Zawinka,
H. Resch,
L. Schmetterer, and
G. Garhofer,
“Intravenous administration of difhenhydramine reduces histamin-induced vasodilator effects in the retina and choroid,”
Invest. Ophthalmol. Visual Sci., 47 1096
–1100
(2006). https://doi.org/10.1167/iovs.05-1174 0146-0404 Google Scholar
T. I. Metelitsina,
J. E. Grunwald,
J. C. DuPont,
G. S. Ying, and
C. C. Liu,
“Effect of Viagra on retinal vein diameter in AMD patients,”
Exp. Eye Res., 83 128
–132
(2006). 0014-4835 Google Scholar
X. Song,
B. W. Pogue,
S. Jiang,
M. M. Doyley,
H. Dehghani,
T. D. Tosteson, and
K. D. Paulsen,
“Automated region detection based on the contrast-to-noise ratio in near-infrared tomography,”
Appl. Opt., 43 1053
–1062
(2004). https://doi.org/10.1364/AO.43.001053 0003-6935 Google Scholar
Early Treatment Diabetic Retinopathy Study Research Group,
“Photocoagulation for diabetic macular edema, ETDRS Report No. 1,”
Arch. Ophthalmol. (Chicago), 103 1796
–1806
(1985). 0003-9950 Google Scholar
Early Treatment Diabetic Retinopathy Study Research Group,
“Fluorescein angiographic risk factors for progression of diabetic retinopathy, ETDRS Report No. 13,”
Ophthalmology, 98 834
–840
(1991). 0161-6420 Google Scholar
B. Sander,
M. Larsen,
C. Engler,
C. Strom,
B. Moldow,
N. Larsen, and
H. Lund-Andersen,
“Diabetic macular edema: a comparison of vitreous fluorometry, angiography and retinopathy,”
Br. J. Ophthamol., 86 316
–320
(2002). https://doi.org/10.1136/bjo.86.3.316 0007-1161 Google Scholar
B. Y. Kim,
S. D. Smith, and
P. K. Kaiser,
“Optical coherence tomographic patterns of diabetic macular edema,”
Am. J. Ophthalmol., 142 405
–412
(2006). 0002-9394 Google Scholar
W. Soliman,
B. Sander, and
T. M. Jørgensen,
“Enhanced optical coherence patterns of diabetic macular oedema and their correlation with the pathophysiology,”
Acta Ophthalmol. Scand., 85 613
–617
(2007). https://doi.org/10.1111/j.1600-0420.2007.00917.x 1395-3907 Google Scholar
U. Christensen,
K. Kroyer,
J. Thomadsen,
T. M. Jørgensen,
M. la Coer, and
B. Sander,
“Quantification of photoreceptor layer thickness using software image enhancing based on stratus optical coherence tomography images,”
Google Scholar
N. Villate,
J. E. Lee,
A. Venkatraman, and
W. EM. Smiddy,
“Photoreceptor layer features in eyes with closed macular holes: optical coherence tomography findings and correlation with visual outcomes,”
Am. J. Ophthalmol., 139 280
–289
(2005). https://doi.org/10.1016/j.ajo.2004.09.029 0002-9394 Google Scholar
W. Drexler,
“Enhanced visualization of macular pathology with the use of ultrahigh-resolution optical coherence tomography,”
Arch. Ophthalmol. (Chicago), 121 695
–706
(2003). https://doi.org/10.1001/archopht.121.5.695 0003-9950 Google Scholar
H. Ishikawa,
M. L. Gabriele,
G. Wollstein,
R. D. Ferguson,
D. X. Hammer,
L. A. Paunescu,
S. A. Beaton, and
J. S. Schuman,
“Retinal nerve fiber layer assessment using optical coherence tomography with active optic nerve head tracking,”
Invest. Ophthalmol. Visual Sci., 47 964
–967
(2006). https://doi.org/10.1167/iovs.05-0748 0146-0404 Google Scholar
|