Prostate cancer consisted 26% of the new cancer diagnoses in 2015, and it was the second most common cause of cancer-related death in the USA.1 Prostate specific antigen (PSA) screening led to an increase in the diagnosis of low- and intermediate risk prostate cancer.2 Focal therapy is emerging as an attractive treatment option for these patients since it only treats the cancerous part of the prostate while leaving the remaining part unharmed. Focal therapy has significantly fewer side effects when compared to radical treatment.188.8.131.52.8.–9 However, accurate imaging on a microscopic level is crucial for case selection and the application of focal therapy.
Imaging technologies that are currently investigated for prostate cancer include multiparametric MRI and transrectal ultrasound (TRUS).10 These technologies provide information on a macroscopic level with voxel sizes of on MRI and pixel sizes of on TRUS. Standard histopathology provides diagnostic information on prostate cancer on a microscopic level. This process, however, is time consuming, moderately reproducible between observers, and produces artifacts on subsequent radiological images.11,12 An imaging modality that provides reproducible instant information on a microscopic level and that can be used intraoperatively and in the outpatient department would make the process of diagnosis, treatment, and monitoring faster (automated image processing) and more accurate (less open for interpretation error).
Optical coherence tomography (OCT) is a high-resolution imaging method in which contrast is based on differences between light scattering of tissue structures.13 Although the imaging depth is limited to , OCT has an important advantage over standard histopathology, because OCT can, either non- or minimally invasively, be applied in vivo, using scan heads or probes in conjunction with balloons, catheters, or needles. Recent developments in needle-based OCT enable minimally invasive puncturing of tissue and therefore in-tissue imaging.14,15 Furthermore, the OCT signal can be quantitatively analyzed, thereby providing tissue-specific parameters such as the optical attenuation coefficient.14,16 The optical attenuation coefficient can function as a measure for tissue density. In previous studies, one to one matching of the OCT images with histology remained challenging,14,16 whereas precise histopathology correlation is essential to draw well-founded conclusions from the measured OCT data and to understand which tissue types in the prostate can be distinguished by OCT. Therefore, a customized prostate measurement and slicing device was developed which facilitates the necessary one to one correlation of OCT images of fresh prostate tissue with histology. It has been demonstrated that OCT is able to characterize tissue structures as seen on histopathology.1718.–19 In addition, the tissue density is expected to translate into a different attenuation coefficient when compared with benign tissue. The objectives of the present study are to identify unique structural characteristics in needle-based OCT images, which by pathology are proven benign or malignant prostatic tissues, based on visual parameters and/or quantitative analyses by means of the optical attenuation coefficients. We hypothesize that the obtained structural characteristics and quantitative parameters can be translated into a diagnostic accuracy that approaches biopsy levels. The study is performed according to the IDEAL guidelines for the validation of medical devices20 and according to the STARD criteria for diagnostic studies.21
A prospective observational ex vivo study was performed in 20 prostates immediately after radical prostatectomy. The hospital’s ethical board waived the need for evaluation.
Patients at least 18 years of age, diagnosed with prostate cancer, and scheduled for radical prostatectomy were eligible for inclusion in this study. Potentially eligible patients were identified at the outpatient clinic of the urology department at the VU University Medical Center in Amsterdam between August and November 2014. Inclusion was on a consecutive basis.
OCT measurements were recorded using a commercially available Imaging System interfaced to a C7 imaging probe (St. Jude Medical, St. Paul, Minnesota) (Fig. 1). The system uses a wavelength of 1300 nm with a bandwidth of 55 nm, with a scanrate of (500 A-lines/frame).22,23 The imaging probe was inserted into a transparent intravascular (IV)-needle catheter (Terumo Surflo© ) allowing tissue puncturing. A detailed description of the device, its adaption for prostate imaging, and the measurement protocol were provided in a previous paper.14 Briefly, the rotating fiber-optic probe with an outer diameter of 2.7 Fr (0.9 mm) produces cross-sectional images with an axial resolution of 10 to and lateral resolution of 20 to .24 A 360-deg probe rotation provides 504 A-scans of . The 504 A-scans were converted into one cross-sectional image or B-scan. In the cylindrical configuration, the B-scan is thus presented with a diameter of 10 mm. The automatic pullback system scans over a trajectory of 54 mm in , producing a 540-frame/B-scan dataset, resulting in a cylindrical scan of the prostate tissue of , with an imaging depth of 2 mm limited by scattering.14
Ex vivo OCT measurements were performed in fresh tissue at the pathology department directly following radical resection. Four (in a small prostate) or six (in a large prostate) intravascular (IV)-catheters were placed in the prostate using the customized tool as described in previous work.25 After removal of the insertion needle, the C7 OCT probe was inserted into the in-situ transparent IV catheter and OCT measurements of the prostate tissue were acquired.
The prostate was sliced after 48 h formaldehyde fixation. For slicing, the knife was guided through the OCT imaging trajectory (Fig. 2).25 Whole mount histopathological slides were produced and stained with hematoxylin and eosin according to the standard VUMC pathology protocol. The histological diagnosis was made by an expert uropathologist who was informed about biopsy and clinically relevant findings but blinded for the OCT results. Afterward, all whole mount slides were digitized. The pathologist performed the marking and annotation of histological structures in the digitized whole mount slides using annotation software designed by our department. Nine histopathological tissue categories were annotated: cystic atrophy, regular atrophy, benign glands, stroma, inflammation, fat and malignancy Gleason pattern 3, Gleason pattern 4, and Gleason pattern 5.
Measurement of the histopathological OCT trajectory length was performed using the pathology annotation software after digital scale calibration to correct for prostate shrinkage after fixation. A starting point was defined at the border of the prostate. Total trajectory length and all distances from the beginning of the trajectory to marked histological structures were measured. This procedure was performed for the histological structures on both sides of the trajectory [ventral (painted red and yellow) and dorsal (painted blue and green)]. The OCT scan creates a 360-deg view of the tissue. Histopathology was two-dimensional (2-D) and thus contained only information in one plane. Rotational matching was ensured by overlapping identical structures seen in histopathology and OCT, as is indicated by the arrows in Fig. 3. The urethra and outer shape of the prostate provided a rough estimation of orientation. Cysts were used for final precise orientation.
The B-scan locations were matched to the distances measured in the histological slides. Based on this information, a virtual overlay of the OCT scan over the histology was made. In this way, the histological classification of every single B-scan was deduced.
Qualitative data analysis
For each of the nine histological categories, locations were selected where opposing sites (ventral and dorsal to the imaging trajectory) were classified in the same category; we refer to these regions as the regions of interest (ROIs). On these locations, we assume that all tissue surrounding the catheter channel has the same histological characteristics. The corresponding OCT B-scans were saved to a database. Subsequently, 110 B-scans across the histological categories were selected by RK, randomized, and offered to two reviewers (BM and AS) blinded for the biopsy histology assessment and further clinical information. Prior to analysis, the two reviewers received a short training based on example images and a short qualitative description of the OCT images. Both reviewers independently assigned the B-scans to one of the nine histological categories for comparison to the histological classification.
Quantitative data analysis
Quantification of the optical attenuation coefficient was performed using custom-written code (Matlab 7.11.0 R2010b, The Mathworks Inc., Natick, Massachusetts). The optical attenuation coefficients from the ROIs were obtained by processing the data according to steps shown in Fig. 4. The region of interest (ROI) is divided into fit regions (FR) with the following dimensions: ( A-lines). Each separate A-line in an FR was laterally averaged, and attenuation coefficient () was fitted according to the method in Appendix yielding 1 per fit region, which correspond to a specific location in the histopathology. These values were used for analysis.
In order to improve the data quality, a margin of error was included for the OCT–histopathology correlation process. Five B-scans, covering 0.5 mm in scan length, at the beginning and end of each ROI were deleted. The remaining ROIs had to contain at least six B-scans for final inclusion. The attenuation coefficients of these regions were stored in a database based on the histologic classification. The determined values for stroma, inflammation, and all malignancy categories were grouped per patient. In case of missing data of whole mount slides or OCT scans, the specific OCT scan/histology combination was excluded from the analysis.
The results of the visually reviewed OCT B-scan were compared to the whole mount histopathology slides. In addition, answers were grouped based on similar histological classification. Finally, all scores were grouped for benign and malignant histology. Sensitivity, specificity, false negative, and false positive calculations were performed manually using a table. The interobserver agreement expressed in kappa was calculated for every category. Kappa () was defined as “poor agreement,” 0 to 0.20 “slight agreement,” 0.21 to 0.4 “fair agreement,” 0.41 to 0.60 “moderate agreement,” 0.61 to 0.8 “substantial agreement,” and above 0.81 “almost perfect agreement.”
Optical attenuation coefficients of malignancy Gleason patterns 3, 4, and 5 were grouped. With a Kolmogorov–Smirnov test, data were tested for a normal distribution. A Mann–Whitney U test was performed for comparison of the attenuation coefficient of stroma and malignancy, using MedCalc v 15.8. A Kruskal Wallis test was performed to compare attenuation coefficients of stroma, inflammation, Gleason 3, and Gleason 4 individually. A -value of stated significance. In addition, receiver operating characteristic (ROC) curve analyses were performed. The area under the curve (AUC) was calculated to determine test accuracy.
Patient characteristics are described in Table 1: 20 consecutive patients were included. On these 20 prostates, 106 three-dimensional (3-D)-OCT measurements were performed. The mean age was 65 years with a mean PSA of at the time of surgery. All patients except one underwent prostate biopsies with Gleason scores ranging from to . One patient was included based on benign prostate hyperplasia transurethral resection results. Histologic specimen examination provided Gleason scores ranging from to and one benign specimen. Seminal vesicle invasion, extracapsular invasion, pelvic lymph node involvement, and positive resection margins were present in two, six, one, and five patients, respectively.
|Age at operation (years)||mean: 64.5 (range: 55 to 76)|
|PSA at operation (ng/ml)||mean: 11.5 (range: 2.9 to 36)a|
|Biopsies (Gleason)||3 + 3|
|3 + 4|
|4 + 3|
|4 + 4|
|Total number of cores||mean: 9.6 (range: 6 to 15)|
|Number of cores positive||mean: 4 (range: 1 to 7)b|
|Specimen (Gleason)||0 + 0|
|3 + 4|
|4 + 3|
|4 + 5|
|Seminal vesicle invasion|
|Pelvic lymph node dissection|
|Lymph node involvement|
|Resection margins positive|
2.9 under Combodart.
One patient underwent a radical prostatectomy based on a histopathological Gleason score 4 + 4 after transurethral resection of the prostate.
OCT and Histology Correlation
A total number of 52 whole mount histopathological slides were available. The OCT trajectories were visible in 50 slides. Malignancy was visible in 19 histology slides in 13 prostates. Twenty-one slides were damaged during production, resulting in areas with data loss. Matching of OCT scans with histology was achieved with high precision in 36 slides. The urethra and outer boundaries of the prostate contributed to the first rotational orientation. Mainly cysts but also other landmarks such as the urethra provided data for precise rotational orientation as can be seen in Fig. 3. Rotational inaccuracy was estimated at 10 deg.
Qualitative Description of B-Scans
The earlier selected representative OCT B-scans that were grouped per histopathological tissue type were reviewed for identification of unique marks on OCT. Figure 6 shows a representative B-scan per category. Review of the B-scans showed that part of the histopathology was well identifiable by OCT based on unique characteristics. Cystic atrophy [Fig. 6(1)] was visually identified by cavities (), divided by septae. The content of the cavities appears opaque due to backscattered light. Regular atrophy [Fig. 6(2)] has smaller (0.1 to 0.3 mm), dark, more grouped cavities than those found in cystic atrophy. Benign glands [Fig. 6(3)] have even smaller, mostly grouped cavities (). The cavities could be dark or opaque. Fat [Fig. 6(4)] had a honeycomb structure, consisting of a unique pattern of bright stripes, alternated by dark dots.
The remaining histopathological classifications were more difficult to distinguish by OCT. Malignancy Gleason patterns 3 and 4 [Figs. 6(5) and 6(6)] present as homogeneous tissue structures. High signal surrounding the probe and low signal in depth represented a low signal penetration (), probably due to a high cell density. Stroma [Fig. 6(7)], similarly to malignant tissue, appeared homogeneous but with an increased signal depth (). Parallel ellipsoidal-shaped lines on both sides of the probe were unique for stroma and probably caused by fibrotic tissue, yet this was not seen in all stromal B-scans. Inflammation [Fig. 6(8)] had a homogeneous pattern with high cell density and signal penetration of , with vague outer boundaries.
Qualitative Blind Assessment of B-Scans
Blind assessment of OCT B-scans was done by two reviewers (see Table 2). The percentages represent the number of B-scans correctly scored by the reviewer. Grouping scores for categories with similar histology and intended treatment improved the results. Combining benign glands, cystic atrophy, and regular atrophy resulted in scores of 100% and 88.9%. The combination of stroma and inflammation resulted in scores of 56% and 59%. Malignancy Gleason patterns 3, 4, and 5 resulted in correct scores of 79% and 88%. In particular, stroma and inflammation were most often misidentified, totaling and out of 41 for reviewers 1 and 2, respectively. Grouping malignant and benign scores gave an overall sensitivity and specificity for malignancy detection of 79% and 88% for reviewer 1, and 88% and 81% for reviewer 2. Negative predictive values were high (94% and 96% for reviewers 1 and 2, respectively). Positive predictive values were relatively low (66% and 57% for reviewers 1 and 2, respectively).
Results of visual OCT B-scan scoring.
|Reviewer 1||Reviewer 2|
|Gleason pattern 3||2/5||40%||3/5||60%|
|Gleason pattern 4||9/19||47%||6/19||32%|
|Benign, cystic, regular||36/36||100%||32/36||89%|
|Malignancy detection (Gl 3/4/5)|
The interobserver agreement on OCT images was calculated and resulted in moderate agreement between observers (weighted kappa of 0.50 for overall test values).
When test results were combined into similar groups with similar histology and intended treatment, benign cystic structures (benign/cystic and regular glands), benign stromal structures (stroma/inflammation), and malignant structures (Gleason pattern 3/4/5), there was substantial agreement between observers (kappa of 0.64). When the results were grouped in two groups: benign structures (stroma and inflammation) and malignant structures (Gleason pattern 3/4/5), there was moderate agreement between observers (kappa of 0.57).
Quantitative Analysis: Attenuation Coefficient
Automated attenuation coefficient calculations were performed for additional differentiation between stroma and malignancy. Twelve patients were included; eight patients could not be included in the analysis, as the tumor or inflammation was not scanned by OCT.
We noticed that the attenuation coefficient (malignant and benign) differed per patient.
The attenuation coefficient was higher for malignancy than for stroma ( versus ). This yields for most patients except four (Table 3, Fig. 7). Data were tested for a normal distribution with a Kolmogorov–Smirnov test, and normality was rejected in both benign and malignant data (). A Mann–Whitney U test showed a significant difference in attenuation coefficient (). Additional analyses with an ROC curve as accuracy measurement of the attenuation coefficient in discrimination of malignancy and stroma showed an AUC of 0.62. For a threshold value of , the test has a sensitivity of 68% and a specificity of 49% (Fig. 8). Separate categories were analyzed with a Kruskal Wallis test (Fig. 9), which showed a significant difference in median attenuation coefficient between stroma, inflammation, Gleason 3, and Gleason 4 (4.6, 4.1, 5.9, and , respectively) (). Since attenuation coefficient strongly differed per patient, we also decided to perform a paired samples Wilcoxon signed-rank test. The dot-and-line diagram is shown in Fig. 10. The paired samples -test does not show a significant difference in optical attenuation coefficient between benign and malignant tissues in the prostate per individual patient .
Overview of the mean attenuation (mm−1) coefficient per histological category per patient. Only patients with malignant OCT data were included in the analysis. (standard deviation) [number of b-scan excluded/total number of b-scans (percentage of b-scans excluded)].
|AC stroma (mm−1)||AC inflammation (mm−1)||AC Gleason 3 (mm−1)||AC Gleason 4 (mm−1)|
|1||4.9 (1.0)||4.5 (0.8)||—||6.5 (1.1)|
|[30/95 (32%)]||[26/80 (33%)]||[1/32 (3%)]|
|2||4.5 (0.8)||—||4.8 (0.8)||—|
|[256/672 (38%)]||[35/80 (44%)]|
|3||5.00 (0.9)||4.8 (0.6)||—||4.9 (0.7)|
|[307/672 (46%)]||[13/32 (41%)]||[184/1136 (16%)]|
|6||4.6 (1.5)||—||—||6.1 (1.4)|
|[57/128 (45%)]||[2/48 (4%)]|
|7||4.6 (0.9)||—||—||4.6 (0.6)|
|[78/208 (38%)]||[89/384 (23%)]|
|8||4.5 (0.9)||4.5 (1.4)||—||5.6 (1.0)|
|[256/672 (38%)]||[8/16 (50%)]||[2/16 (13%)]|
|10||4.8 (0.8)||—||—||4.6 (0.6)|
|[98/400 (25%)]||[15/48 (31%)]|
|11||4.8 (1.1)||—||—||5.2 (0.8)|
|[55/96 (57%)]||[37/128 (29%)]|
|13||4.3 (0.8)||—||6.1 (0.7)||6.1 (0.4)|
|[93/144 (65%)]||[1/48 (2%)]||[4/16 (25%)]|
|17||4.6 (0.7)||—||—||4.9 (0.9)|
|18||6.4 (2.1)||—||—||5.4 (1.0)|
|20||4.4 (0.8)||—||—||4.3 (0.8)|
|Mean||4.4 (0.9)||4.7 (1.0)||5.5 (1.0)||5.1 (0.9)|
This study demonstrates that needle-based OCT can identify unique tissue patterns in benign and malignant prostatic tissue and is even able to distinguish benign tissue from malignant prostate (). Visually, all types of prostatic tissue were identifiable on OCT, although malignancy, stroma, and inflammation present as similar patterns and therefore it was more challenging to differentiate between those categories. Additional computed OCT analysis by means of the attenuation coefficient did contribute to discrimination of benign from malignant tissue in the prostate. However, when analyzed on a per patient basis, there was no significant difference in optical attenuation coefficient. The unique aspect of this study is the one to one correlation of OCT and histopathology in prostate cancer. In previous work, the correlation of histopathology and OCT scans was a serious drawback.10,14 For this reason, a customized tool for OCT measurement and prostate slicing was developed. The feasibility of precise correlation of histopathology and OCT of the device was described in earlier work.25 Using this method, we performed a single blind qualitative and quantitative accuracy study of OCT in a larger cohort of 20 patients ex vivo.
Limitations on Data Matching
The correlation of 3-D OCT scans and 2-D histopathology was based on two assumptions. First, OCT probe rotation in the prostate tissue was approximated by corresponding structures seen in OCT and histology. The estimated rotational inaccuracy is 10 deg, which we consider a small inaccuracy. Second, when both sides of the trajectory contained identical histological structures, it was assumed that all tissue surrounding the OCT probe was from the same histological category. Although this seems plausible, histopathology showed a large heterogeneity in prostatic tissue; therefore, it might be possible that some areas contain tissue from another category. In order to minimize measurement errors, we identified five B-scans on both sides of the ROIs as the margin of error and excluded them from data analysis.
Because the OCT-probe only samples part of the tissue, we missed the tumor in eight patients. For this reason, we did not include this OCT data in the visual assessment and we were not able to analyze the attenuation coefficient on malignancy in these patients.
Limitations in Visual Assessment
For visual assessment, the observers were affiliated with OCT but were not formally trained for detection of prostate cancer. Therefore, it is plausible that their assessment scores can be improved by OCT assessment training. This training could consist of an explanation of the unique tissue characteristics followed by an extended test of prostatic B-scans. This method has been successfully executed before in other imaging studies, e.g., using MRI to visualize prostate cancer, training significantly increased diagnostic accuracy.27,28
Limitations in Quantitative Assessment Using the Optical Attenuation Coefficient
It was remarkable that the optical attenuation coefficient for one category (e.g., Gleason 4) could differ within a patient with values ranging 1 to . This intrapatient variation might be because prostate tissue is heterogonous and is therefore challenging for the pathologist to delineate an area as benign or malignant. Moreover, the pathologist analyzes the whole mount prostatic histology slide using a microscope and delineates this using software on a digital lower resolution image of the same slide, a process prone to errors. We partially corrected for this by removing 5 B scans from each side of each ROI as described earlier. Yet, it is also known from histological studies that healthy prostate tissue differs between patients and therefore tumor tissue might differ between patients as well.29 Furthermore, a fair amount of fits () was excluded from the analysis because the calculated fit did not exactly match the measured data. Most likely, these specific regions of the prostate are too heterogeneous for a correct fit, e.g., cysts, but they might still have a clinical significance. Technology such as automated texture analysis and pattern recognition could contribute to solve this problem and could increase the performance of OCT for prostate cancer detection.30 Furthermore, when more data are available, convolutional neural networks could be trained to perform the analyses on the OCT data, as it now starting in conventional histopathology.31
Perspective to Other Work in the Field
Ex vivo application of OCT on the human prostate was first described in the year 2000. 2-D cross sectional OCT scans were generated of specimens after radical prostatectomy.32 It was claimed that OCT could distinguish malignant from benign prostate tissue on the basis of architectural differences in the tissue. OCT has also been described for prostatic nerve identification in order to spare them during radical prostatectomy experiments.33,34 Dangle et al.35 evaluated surgical margins, seminal vesicle invasion, and capsular invasion with OCT ex vivo in prostates after radical prostatectomy. They found that OCT overestimated the amount of margin involvement, but the NPV was high. For surgical margins on OCT sensitivity, specificity, PPV, and NPV was 70%, 84%, 33%, and 96%, respectively. Recently, a study by Lopater et al.36 described the initial application of full-field OCT for the detection of prostate cancer in prostate biopsies. This ex vivo OCT imaging method approaches the resolution of traditional histological slides with a resolution of about . The images were scored solely by architectural structures and solely in the 2-D setting, sensitivity, specificity, PPV, and NPV values were 63%, 74%, 55.5%, and 80% respectively. In our study, sensitivity (79% and 88%) and specificity (88% and 81%) for malignancy detection were higher than in these previous studies. PPVs (66% and 57%) and NPVs (94% and 96%) were higher or comparable to these previous studies. Although PPVs were higher than previous studies, they were still considerably lower than the NPVs. This indicates that OCT, as it is investigated here, can be a good test to exclude disease in the prostate. The advantage of the St. Jude Intravascular Imaging System is that it is commercially available and applicable in a sterile environment. Consequently, the application of this OCT system in the outpatient clinic or even in a surgery setting is relatively easy. The small probe diameter even allows for in vivo insertion in the prostate. Images are acquired in 3-D, which increases the amount of information substantially when compared to 2-D histopathology. Needle-based OCT lacks the spatial resolution to function as a stand-alone diagnostic modality in prostate cancer. However, it can make the process of conventional histopathology faster and possibly more accurate. Since OCT results are quantifiable, it can reduce interrater variability, which is high for conventional histopathology.37
This study validates one to one correlation of histopathology and OCT using the uniquely designed prostate slicing device in a cohort of 20 patients. Correlation showed that most histological tissues have a unique pattern and therefore could be visually identified on OCT, such as cysts, lines, view in depth, or signal intensity. The abilities of OCT for prostate cancer identification were explored and qualitative visual analysis confirmed the hypothesized high cell density of malignant tissue. In addition, the optical attenuation coefficient contributes to the differentiation between stroma and malignancy, although a per-patient analysis did not show a significant difference. These findings may serve as a basis for an in vivo study combining OCT and prostate biopsy histopathology; allowing comparison of both tests and exploring the clinical potential of OCT in digital pathology of the prostate.
Fitting the Optical Attenuation Coefficient
Quantification of the optical attenuation coefficient was performed using custom-written code (Matlab 7.11.0 R2010b, The Mathworks Inc., Natick, Massachusetts). Methods were derived from our previous work.14,16,26 First, OCT amplitude data were loaded into the software. Isolation of tissue-related data from the original data set was achieved by a succession of image processing steps applied to all B-Scans. First, low SNR regions were excluded by applying a pixel value threshold on the amplitude data. Next, to remove the catheter from the B-Scans, an algorithm utilizing region connectivity as well as a priori knowledge about the catheter-tissue geometry was implemented. Following this, a Gaussian low pass filter was applied to the catheter-free images; this step served to improve performance of the ensuing edge detection using the Sobel method. The detected tissue edge was then smoothed further by carrying out local regression using weighted linear least-squares (assigning lower weight to outliers and zero weight to data outside six mean absolute deviations) and a second-degree polynomial model. OCT amplitude data of the prostate tissue were isolated from the original B-Scan by selecting all data below this smoothed edge. Finally, these data were straightened and then divided into FRs with the following dimensions: in the fast axis, slow axis, and depth (located 20 pixels below the tissue edge), respectively. These individual FRs were then laterally averaged, yielding a single average A-line for each FR. Subsequently, an FR-specific attenuation coefficient () was determined by nonlinear least squares fitting the following equation to this average A-line: . is the averaged OCT amplitude in depth, is the position in depth, is the position of the tissue boundary, and are free running parameters (amplitude and attenuation coefficient, respectively). The fit always starts 20 pixels below the tissue surface to ensure the absence of tissue edge reflection. The noise is defined as the average of the last 100 pixels of the data. The system dependent parameter (describing the confocal point spread function) is defined as: , where is the position of the focus in depth, is the Rayleigh length, and is the refractive index of the medium. The sensitivity roll-off is defined as: . Here, is the maximal imaging depth of the OCT system and is the ratio between the spectral resolution to the sampling interval. The contribution of the confocal point spread function and the sensitivity roll-off were determined by fitting the equation for to the OCT amplitude of a sample with a low concentration of Intralipid (0.003%) using four different catheters, for which negligible scattering is assumed, and and were the free running parameters.38,39 The median values for and are 0.73 mm and 1.4, respectively, which corresponds to reported values.40,41 After fitting, all fits were judged by an experienced observer and nonmatching fits (e.g., because of a cyst in the data) were excluded from the analysis. In this way, only accurate representations of the optical attenuation coefficient were included in the analysis. Finally, all obtained attenuation coefficient values were multiplied by 1.4 (refractive index for tissue) yielding values per mm.42
This study was funded by an unrestricted research grant from the Cure for Cancer Foundation.