the research protocol was approved by the institutional review board at the oregon health & science university ( ohsu ) and performed in accordance with the tenets of the declaration of helsinki . written informed consent was obtained from each participant after explanation of the nature and possible consequences of the study .
participants were recruited at the casey eye institute / ohsu according to the aig study protocol .
the ohsu ancillary site followed the same eligibility and endpoint protocol as the aig study , but used an advanced swept - source oct prototype system instead of commercially available oct instruments .
briefly , normal control participants met the following criteria in both eyes : iop of less than 21 mm hg in both eyes , and a normal humphrey visual field ( hvf ) on achromatic standard automated perimetry by swedish interactive threshold algorithm 24 - 2 testing ( hfa ii ; carl zeiss meditec , inc . , dublin , ca , usa ) with mean deviation ( md ) , glaucoma hemifield test ( ght ) , and pattern standard deviation ( psd ) within normal limits .
in addition , normal subjects had a normal - appearing optic nerve head ( onh ) and nfl on ophthalmoscopic examination , and an open angle by gonioscopy .
glaucoma participants were required to have a glaucomatous onh rim or nfl thinning on ophthalmoscopic examination .
glaucomatous eyes were classified in the pg subgroup if they have visual field ( vf ) psd or ght outside normal limits ( p < 0.05 and p < 1% , respectively ) in a consistent pattern on two qualifying vf exams .
exclusion criteria for all groups included visual acuity less than 20/40 , age < 40 or > 79 years at enrollment , any ocular surgery other than uncomplicated cataract extraction , other diseases that might cause vf or onh abnormality , and factors that might preclude the participant from performing study procedures or completing the study .
a prototype high speed swept - source fourier - domain oct system was used in this study .
the device operated at an axial scan speed of 100 khz using a swept - source cavity laser operating at 1050 nm with a tuning range of 100 nm .
a resolution of 5.3 m axially and 18 m laterally at an imaging depth of 2.9 mm in tissue was achieved .
the ocular light power exposure was 1.9 mw , within the american national standards institute ( ansi ) safety limit .
participants were scanned using a high density 8  8 mm raster scan pattern , centered at the onh . in the fast transverse scan direction ,
the time to acquire each 3d volumetric scan was 4.3 seconds . in the scan protocol ,
one eye of each participant was scanned with 4 scans consisting of 2 horizontal and 2 vertical scans .
an orthogonal registration algorithm was applied to register and merge all 4 scans to reduce eye motion .
concentric , cylindrical , cross - sectional oct images of varying radii were resampled from the volume scan ( fig .
we selected rings with radii of 1.5 , 1.7 , 2 , 2.5 , 3 , and 3.5 mm for analysis .
the disc center was decided manually by one of us ( ll ) by matching the disc boundary with an ellipse on the en face view ( fig .
( a ) en face view ( average projection of log reflectance intensity ) of the 8  8 mm volumetric  cube  scan of the onh region in a glaucomatous eye . the 3-mm radius circle centered on the onh
( b ) the 3-mm radius cylindrical cross section is resampled from the cube scan .
the ilm contour ( white dashed line ) is measured from this circular cross section .
( c ) en face oct showing resampling circles at radii of 1.5 , 1.7 , 2.0 , 2.5 , 3.0 , and 3.5 mm .
the ilm was detected as the most inner boundary of retina in the cylindrical cross - section images .
an automated algorithm was developed to follow the vessel bumps and avoid artifacts , such as points of vitreous macular adhesion . the ilm elevation profile ( contour )
one of us ( ll ) inspected the cross - sectional images with ilm contour overlay and performed manual correction of the contour if necessary .
the ilm elevation profile was transformed to spatial frequency domain using the fast fourier transform ( fig .
the spatial frequency components at lowest frequencies were removed because they simply relate to the average axial position and the tilt of the oct scan beam relative to the onh plane .
the highest frequency components also were nondiagnostic because of high noise and little anatomic information .
the middle frequency band that optimized the detection of retinal vessel relief and nfl bundle defects was determined empirically based on the difference between the average normal and glaucoma ilm spatial frequency spectrum ( fig .
the rscv is defined as the average log value within a middle spatial frequency band of the fourier transform to the elevation profile of the inner retinal surface .
( a ) spatial frequency spectrum of ilm elevation showing the averages of all glaucoma eyes and all normal eyes derived from the 3.0-mm oct cylindrical sections .
( b ) shows the difference spectrum ( glaucoma  normal ) and the optimize frequency band ( between dotted lines ) used to compute the rscv .
optimized passbands for rscv calculation at different radii . on the cylindrical cross - sectional oct image with radius of 1.7 mm
the wilcoxon rank sum test was performed to determine statistical significance between the study groups .
the area under receiver operating characteristic curve ( aroc ) was calculated for the diagnostic accuracy .
pearson correlation was used to determine correlations between rscv , axial length , nflt , and visual field test md . a coefficient of variance ( cv )
the first one compares rscv value at r = 2.8 mm and r = 3.0 mm .
the second test compares results graded independently by two of the authors ( ot and ll ) .
all image processes and statistical analyses were done using matlab ( mathworks , natick , ma , usa ) .
the research protocol was approved by the institutional review board at the oregon health & science university ( ohsu ) and performed in accordance with the tenets of the declaration of helsinki . written informed consent was obtained from each participant after explanation of the nature and possible consequences of the study .
participants were recruited at the casey eye institute / ohsu according to the aig study protocol .
the ohsu ancillary site followed the same eligibility and endpoint protocol as the aig study , but used an advanced swept - source oct prototype system instead of commercially available oct instruments .
briefly , normal control participants met the following criteria in both eyes : iop of less than 21 mm hg in both eyes , and a normal humphrey visual field ( hvf ) on achromatic standard automated perimetry by swedish interactive threshold algorithm 24 - 2 testing ( hfa ii ; carl zeiss meditec , inc . , dublin , ca , usa ) with mean deviation ( md ) , glaucoma hemifield test ( ght ) , and pattern standard deviation ( psd ) within normal limits .
in addition , normal subjects had a normal - appearing optic nerve head ( onh ) and nfl on ophthalmoscopic examination , and an open angle by gonioscopy .
glaucoma participants were required to have a glaucomatous onh rim or nfl thinning on ophthalmoscopic examination .
glaucomatous eyes were classified in the pg subgroup if they have visual field ( vf ) psd or ght outside normal limits ( p < 0.05 and p < 1% , respectively ) in a consistent pattern on two qualifying vf exams .
exclusion criteria for all groups included visual acuity less than 20/40 , age < 40 or > 79 years at enrollment , any ocular surgery other than uncomplicated cataract extraction , other diseases that might cause vf or onh abnormality , and factors that might preclude the participant from performing study procedures or completing the study .
a prototype high speed swept - source fourier - domain oct system was used in this study . the device operated at an axial scan speed of 100 khz using a swept - source cavity laser operating at 1050 nm with a tuning range of 100 nm .
a resolution of 5.3 m axially and 18 m laterally at an imaging depth of 2.9 mm in tissue was achieved .
the ocular light power exposure was 1.9 mw , within the american national standards institute ( ansi ) safety limit .
participants were scanned using a high density 8  8 mm raster scan pattern , centered at the onh . in the fast transverse scan direction ,
the time to acquire each 3d volumetric scan was 4.3 seconds . in the scan protocol
, one eye of each participant was scanned with 4 scans consisting of 2 horizontal and 2 vertical scans .
an orthogonal registration algorithm was applied to register and merge all 4 scans to reduce eye motion .
concentric , cylindrical , cross - sectional oct images of varying radii were resampled from the volume scan ( fig .
we selected rings with radii of 1.5 , 1.7 , 2 , 2.5 , 3 , and 3.5 mm for analysis .
the disc center was decided manually by one of us ( ll ) by matching the disc boundary with an ellipse on the en face view ( fig .
( a ) en face view ( average projection of log reflectance intensity ) of the 8  8 mm volumetric  cube  scan of the onh region in a glaucomatous eye . the 3-mm radius circle centered on the onh
( b ) the 3-mm radius cylindrical cross section is resampled from the cube scan .
the ilm contour ( white dashed line ) is measured from this circular cross section .
( c ) en face oct showing resampling circles at radii of 1.5 , 1.7 , 2.0 , 2.5 , 3.0 , and 3.5 mm .
the ilm was detected as the most inner boundary of retina in the cylindrical cross - section images .
an automated algorithm was developed to follow the vessel bumps and avoid artifacts , such as points of vitreous macular adhesion . the ilm elevation profile ( contour )
one of us ( ll ) inspected the cross - sectional images with ilm contour overlay and performed manual correction of the contour if necessary .
the ilm elevation profile was transformed to spatial frequency domain using the fast fourier transform ( fig .
3a ) . the spatial frequency components at lowest frequencies were removed because they simply relate to the average axial position and the tilt of the oct scan beam relative to the onh plane .
the highest frequency components also were nondiagnostic because of high noise and little anatomic information .
the middle frequency band that optimized the detection of retinal vessel relief and nfl bundle defects was determined empirically based on the difference between the average normal and glaucoma ilm spatial frequency spectrum ( fig .
the rscv is defined as the average log value within a middle spatial frequency band of the fourier transform to the elevation profile of the inner retinal surface .
( a ) spatial frequency spectrum of ilm elevation showing the averages of all glaucoma eyes and all normal eyes derived from the 3.0-mm oct cylindrical sections .
( b ) shows the difference spectrum ( glaucoma  normal ) and the optimize frequency band ( between dotted lines ) used to compute the rscv .
on the cylindrical cross - sectional oct image with radius of 1.7 mm , the lower nfl boundary also was detected using a method described previously .
the wilcoxon rank sum test was performed to determine statistical significance between the study groups .
the area under receiver operating characteristic curve ( aroc ) was calculated for the diagnostic accuracy .
pearson correlation was used to determine correlations between rscv , axial length , nflt , and visual field test md . a coefficient of variance ( cv )
the first one compares rscv value at r = 2.8 mm and r = 3.0 mm .
the second test compares results graded independently by two of the authors ( ot and ll ) .
all image processes and statistical analyses were done using matlab ( mathworks , natick , ma , usa ) .
of 17 participants in the glaucoma group ( 8 in the pg and 9 in the ppg subgroups ) , 15 had early vf damage ( md > 6 db ) and 2 had more severe vf damage .
we enrolled 25 participants in the normal control group , but due to their younger average age , only the 17 oldest control participants were used in the present analysis .
the age matching was not exact , as the control participants still averaged 9 years younger ( table 1 ) .
as expected , the glaucoma group had significantly worse vf md and thinner nfl in the overall area , superior and inferior quadrants .
the image quality of the normal group is significantly higher than that of the glaucoma group based on signal strength index .
participant characteristics the spatial frequency transforms of the ilm profile were averaged for the normal and glaucoma groups ( fig .
optimized bands were detected at frequency components where the glaucoma group was significantly larger than the normal group ( fig .
we averaged rscvs in the optimized bands for each circle and the rscv distribution for normal and glaucoma groups are illustrated in figure 5 .
the rscv was significantly higher at all sampling radii in the glaucoma group compared to the normal control group ( table 2 ) .
distribution of rscv values for normal ( n ) and glaucoma ( g ) groups at sampling radii between 1.5 and 3.5 mm .
retinal surface contour variability the diagnostic accuracy of nfl and rscv were assessed using aroc values at each sampling radius ( table 3 ) .
the diagnostic accuracy of rscv improved with larger radii and the highest aroc value ( 0.90 ) was found at a radius of 3.5 mm .
this was higher than that for nflt parameters , but only significantly in the temporal and nasal quadrants ( < 0.01 ) .
we examined averages of rscv over various ranges of radii and found that the average over 2.5 to 3.5 mm had the highest aroc ( 0.91 ) .
the combination had marginally higher aroc ( 0.91 ) than nfl only ( p = 0.08 ) . diagnostic accuracy comparison venn diagram analysis ( fig .
6 ) showed that rscv had higher sensitivity for glaucoma diagnosis ( 53% ) than nflt ( 29% ) .
it also was notable that all abnormal eyes detected by nflt also were detected by rscv .
however , the difference in sensitivity was not statistically significant ( p = 0.13 ) .
the threshold for rscv was set 2.33 sd above the mean of the normal group ( 99 percentile cutoff assuming normal distribution ) and the threshold of nfl was 2.33 sd below the mean of normal group . in the normal group ,
all eyes were within the cutoff for both parameters , indicating 100% specificity . in the glaucoma groups , rscv detected abnormality in nine participants , including the five detected by nfl .
the correlations between averaged rscv and nflt , or between averaged rscv and vf md , were estimated for glaucomatous participants .
significant correlations were found between rscv and nfl ( r = 0.54 , p = 0.03 ) but not md ( r = 0.32 , p = 0.21 ) .
a significant correlation also was found between nfl and md ( r = 0.68 , p = 0.002 ) .
we also evaluated the correlation between rscv and nfl for normal participants , which was not significant ( r = 0.10 , p = 0.70 ) .
we tested whether the rscv was sensitive to transverse magnification changes ( i.e. , axial eye length variation ) by comparing the values evaluated at the sampling radii of r = 2.8 and 3 mm .
no significant difference was found between the two radii in either the glaucoma or normal group ( p > 0.26 ) .
1.3% for glaucoma eyes . we also tested the between - grader ( ot and ll ) reproducibility of rscv because the onh boundary and center were defined manually .
the reproducibility ( cv ) was 1.6% for normal eyes and 1.1% for glaucoma eyes .
we did not find a significant correlation between axial length and rscv parameter ( r = 3.5 mm ) for the normal ( r = 0.06 , p = 0.83 ) and glaucoma ( r = 0.42 , p = 0.09 ) groups .
we also did not find a significant correlation between axial length and rscv parameters at any radius for the normal group ( p > 0.32 ) .
in this pilot study , we reported , to our knowledge , the first use of spatial frequency analysis to detect small scale retinal surface contour change in glaucoma .
the rscv , measured in the peripapillary area , was significantly larger in glaucoma eyes .
the diagnostic accuracy of rscv was at least as good as nflt , although a larger study would be needed to see if there is a significant advantage .
in early glaucoma , nflt has limited diagnostic sensitivity , with literature values of 33.3% to 67.4% at a fixed specificity of 99% .
this is not surprising as there is a relatively wide range of population variation in the average nflt , which has been reported as between 9.0% and 10.0% sd . because the rscv measures small scale , focal changes that are unlikely to be due to the inborn variation of nflt , it may be able to better detect glaucoma in the earlier stages , when nfl damage is small and variable in the location . in this study ,
nflt was abnormally ( 99% specificity cutoff ) thin in only 5 of 17 participants in the glaucoma group , while rscv was abnormally elevated in 9 of 17 glaucoma participants .
while this advantage is not statistically significant , it shows a promising trend that deserves further study . because the nflt and the rscv measure different aspects of nfl damage 
while we previously have developed another parameter to measure focal peripapillary nfl and macular ganglion cell complex ( gcc ) thinning called focal loss volume ( flv ) , that algorithm was optimized to detect relatively larger areas ( 500 m superpixels ) of focal damage . by contrast
the rscv spatial frequency analysis detects changes in the 2 to 100 radians / mm , which corresponds to spatial features as small as 60 m ( the diameter of a small arteriole ) .
thus , we believe all three parameters will be synergistic based on the fact they measure glaucoma damage on different spatial scales .
we found that the rscv appeared to have higher diagnostic accuracy in circles with large radii ( 2.53.5 mm ) than small radii ( 1.52 mm ) .
firstly , we observed that , in the normal eye , larger retinal vessels submerge into the retina further from the disc margin .
thus , the inner retinal surface is smoother and rscv is lower further from the disc .
this observation also might be due to the fact that glaucoma damages axon bundles nonuniformly within the onh ( laminar pores , and so forth ) and these bundle losses may be more detectable within the nfl where all bundles are more spread out .
another assumption is that the thinner nfl away from the disc is going to be less likely to mask blood vessels and loss of axon bundles may reveal blood vessels more readily in these regions .
the second observation is that the population variances were tighter with larger radii ( table 2 ) .
we suspected this was simply due to the larger number of a - scans in sampling circles of larger radii .
thus , maximizing the oct scan density in the 2.5- to 3.5-mm radius annulus around the optic disc may improve rscv evaluation .
an additional advantage of rscv over nflt is that the rscv value is relatively insensitive to magnification variation or decentration relative to the optic disc center . between the scan radii of 2.5 to 3.5 mm
, the mean rscv values varied less than 2% in the normal and glaucoma groups .
intergrader variation in disc boundary identification also affected rscv by less than 2% . by comparison
furthermore , nfl quadrant thickness averages are known to be exquisitely sensitive to decentration of the sampling circle relative to the optic disc .
therefore , rscv may be a more robust diagnostic parameter that is relatively resistant to the introduction of bias due to patients , oct operators , or graders .
in previous studies , contour variability was accessed by heidelberg retinal tomography ( hrt ) but was considerably inferior in its ability to diagnose glaucoma or progression compared to nfl measured by oct .
we believe there are two reasons why contour variability by hrt was considerably inferior in diagnostic power .
thus , it is difficult for hrt to catch small changes in the contour variability caused by protrusions of retinal vessels .
, we found that the contour variability is less sensitive in location near the disc edge .
one of the limitations of this study is the relatively small number of participants , especially for advanced glaucoma ( vf < 6db ) .
a power calculation shows that our sample size of 17 normal and 17 glaucoma subjects should be able to detect a difference in aroc of 0.67 compared to the aroc of rscv of 0.90 with a probability of 0.80 at a significance threshold of p < 0.05 .
however , rscv was not significantly affected by sex in our control group ( p = 0.62 ) .
second , the glaucoma group also is 9 years older than the normal group . as a consequence , the oct image quality in normal
however , the ssi is good for the normal ( 60  9 ) and glaucoma ( 50  12 ) groups .
the rscv parameters are not sensitive to ssi when the image quality is good enough for accurate ilm boundary detection .
we also tested the correlation between ssi and rscv but did not find any significance for the normal and glaucoma groups ( p > 0.38 ) .
furthermore , rscv was not correlated with age in the normal group ( p = 0.43 ) .
retinal surface contour variability also may be sensitive to focal nfl loss caused by other optic neuropathies , such as anterior ischemic optic neuropathy ( aion ) and optic neuritis .
since rscv is a nonspecific detector of nfl loss on a small spatial scale , specific diagnosis of particular conditions would require the physician to consider the clinical presentation as well as patterns of nfl , gcc , optic disc topography , and vf changes . in summary ,
since the rscv detects small - scale focal damage and the average nflt measures global damage , they provide different diagnostic information that may be synergistic .
 purposewe investigated the feasibility of glaucoma detection by measuring retinal surface contour variability ( rscv ) using optical coherence tomography ( oct).methodsthe peripapillary region in one eye of each participant was scanned over an 8  8 mm area with a swept source oct prototype . 
 the retinal surface contour was sampled at approximately 1.5- to 3.5-mm radius circles centered on the optic nerve head . 
 the rscv is defined as the average log value within a middle spatial frequency band of the fourier transform to the elevation profile of the inner retinal surface . 
 the spatial frequency band was optimized to distinguish glaucoma from normal . 
 nerve fiber layer thickness ( nflt ) was sampled around a 1.7-mm radius circle . 
 glaucoma severity was assessed by automated static perimetry.resultswe enrolled 17 glaucomatous eyes and 17 healthy eyes . 
 a great majority of the glaucoma group were in the early stage ( visual field mean deviation average 2.48  3.73 db ) . 
 significant differences were found for rscv between glaucoma and control eyes ( p < 0.003 ) at all radii . 
 the area under the receiver operating characteristic curve ( aroc = 0.90 ) of rscv was best at the 3.5-mm radius . 
 this was not significantly better than nflt ( aroc = 0.84 ) . with the 99% specificity , 
 the glaucoma detection sensitivity was 53% for rscv and 29% for nflt 
 ( p = 0.13).conclusionsretinal surface contour variability was significantly increased in glaucoma patients . 
 the diagnostic accuracy of rscv was equal to nflt in early glaucoma . 
 since the rscv detects small - scale focal damage and the average nflt measures global damage , they provide different diagnostic information that may be synergistic . 
Method
Participants
OCT and Scanning
ILM Detection
Retinal Surface Contour Variability (RSCV)
Nerve Fiber Layer Thickness (NFLT)
Statistical Analysis
Results
Discussion