Optical coherence tomography (OCT) is an imaging modality with a high resolution.1 Since the light can be transmitted via an optical fiber, OCT has been introduced into intravascular imaging by embedding an optical fiber into an imaging catheter.2 The acquired intravascular OCT (IVOCT) images with a high resolution (10 to ) boosted the understanding of the composition of atherosclerosis in the coronary artery, which is one of the fatal diseases that is responsible for most of the death in the world.3 In IVOCT images, different tissue types can be recognized,4 e.g., the fibrous tissue is bright, both lipid-rich and calcified plaque appear as low intensities, while the border of the former tends to be diffused and the latter is sharper.
To achieve automated analysis of IVOCT images for diagnosis and risk assessment, much attention has been spent on the quantitative characterization of the tissue composition of the coronary vessel wall.188.8.131.52.184.108.40.206.13.14.–15 Inspired by the literature on tissue optics in general OCT,16220.127.116.11.–21 one of the important features that has been analyzed is the OCT-derived attenuation and backscatter coefficient. For accurate estimation, the modeling of the transmitting light through biological tissue usually requires taking into consideration the complex interactions18,19 such as multiple scatter and concentration dependency. On the other hand, estimations have been reported at an adequate accuracy for tissue quantification based on a simplified model involving a single exponential decaying term.20,21
Based on the simplified model, one of the commonly used approaches for the estimation of the attenuation is the least square fitting either an exponential curve to a one-dimensional (1-D) OCT signal or a linear function to the logarithm of 1-D OCT signal. The curve-fitting (CF) attenuation has been widely used to quantify biological tissue such as bladder urothelial carcinoma,22 human skin,23 nasopharyngeal carcinoma cells,24 collagen content in ovarian tissue,25 breast tissue,26 and human atrial tissue.27 Specifically for the analysis of tissue in the coronary arteries, an earlier work of Xu et al.8 has shown that the CF approach can be used to estimate both the attenuation and backscatter coefficients, which can improve the differentiation of atherosclerotic tissue. In order to improve the reproducibility and the robustness of the CF approach, Gargesha et al.28 estimated the attenuation coefficient within volumes of interest. The method has been validated with different combinations of noise reduction filtering and outlier rejection to optimize a goodness-of-fit. Furthermore, attenuation and backscatter coefficients were analyzed to characterize the neointima formation after the stent implemented from different manufacturers.29,30 van Soest et al.9 proposed an automated optimization scheme to estimate the attenuation coefficient by expanding fitting windows within regions assumed to be homogeneous. This approach has been validated with a later study from the same group.31 Since the automated algorithm allows for a high robustness and reproducibility, the CF-attenuation is increasingly used for region-based tissue characterization combined with machine learning tools, such as random forest,32 Gaussian likelihood map,26 and supported vector machine.27
Another approach for the estimation of attenuation coefficient has been recently introduced by Vermeer et al.10 It is based on a depth-resolved (DR) model and allows for local pixel-wise estimation with two fundamental assumptions, the light should be attenuated completely and the backscatter and attenuation are linearly related to each other. Despite strict assumptions that have been made, their measurements did show a good consistency in both homogeneous and layered phantom experiments. This DR attenuation has been used for the detection of the cerebral edema in OCT images.33 Taking the confocal function and the fall-off of the system into consideration, the DR model has been reported to estimate the attenuation more accurately and potentially further improve the pixel-wise tissue characterization in combination with a noise removal approach.34 In addition to estimating the attenuation, it should be noted that the DR model has been widely used for the enhancement of OCT images.3536.37.38.–39 Due to its recent introduction, no quantitative analysis with the DR model for tissue characterization has been reported on IVOCT images.
The aim of this work is to explore the potential of DR attenuation coefficient, backscatter term, and image intensities in distinguishing different atherosclerotic tissue types with a robust implementation of the DR approach. In order to achieve this goal, the DR attenuation is introduced to IVOCT, the model is further extended to estimate a backscatter term, and a cut-off border is determined to remove noisy regions. Six types of tissue were labeled by two pathologists referring to their histological counterparts. Combining with the image intensities, the estimated attenuation coefficient and the backscatter term are analyzed for the differentiation of these six types of tissues. To the best of our knowledge, no such study has been published so far for intravascular tissue analysis.
The structure of this paper is organized as follows. In Sec. 2.1, the light propagation model and the DR-attenuation are introduced, the estimation of a backscatter term is derived in Sec. 2.2, and the DR model is combined with the CF approach to determine the cut-off depth border in Sec. 2.3. The experimental materials and the implementation of the experiment are described in Sec. 3. The results are reported in Sec. 4 and further discussed in Sec. 5, where the estimated values are extensively compared with that in the literature. The conclusion is drawn in Sec. 6.
Depth-Resolved Model Analysis
We continue to use the model of light transmission as it has been introduced in the work of Vermeer et al.1019 and the converting factor of the digitizer, is the source light power (), the attenuation is denoted as , and the backscatter coefficient is denoted as .
With an assumption that the backscattered light is a fixed fraction to the attenuated light the integration from to infinite of Eq. (1) will be
The attenuation can be estimated with the OCT image intensities as
Backscatter Term Estimation
Suppose an A-line from a homogeneous region starts at depth and ends at depth . Within this region, both the attenuation and the backscatter coefficients at each depth are the same and , thus Eq. (1) can be simplified as
Within the homogeneous region, and can be estimated as the slope and the offset of the linear function fitted with , as it has been proposed that in the work of Schmitt et al.16 Noting that the is a constant relative to , linear fitting gives the same estimation of the attenuation coefficient. However, the estimation of is necessary for the estimation of backscatter coefficient.
With the DR model, we consider and are the start and the end of a pixel, then approximates to 1 since is then very small. Dividing Eq. (5) by a factor of and taking the logarithm, we get a similar intersection to be our backscatter term as
DR-CF Cut-Off Algorithm
One of the main concerns regarding the pixel-wise DR algorithm is the assumption that the integral of intensity should be taken to infinity. However, infinity in practice is usually substituted with a large value, the maximum depth of the signal in the case of OCT images. This causes little effect on the estimated attenuation of the superficial tissue region, where the signal is much stronger than the noise. However, as it approaches the end of an A-line, the estimated attenuation coefficient with the noise-dominant signal will be exaggerated and thus is unreliable.
In order to avoid this, we applied the DR-CF algorithm to determine a cut-off depth. Noting that
The second term is a constant, the logarithm of the integral intensities follows a linear function with respect to the depth within this region, and the attenuation can be calculated from the slope. With the DR-CF algorithm, the left term is always decreasing, thus it guarantees that the estimated values are positive and the high noise tolerance of the CF-approach is maintained. Since the deeper signals are influenced strongly by noise, it is difficult to perceive the local variations, but the general decreasing trend can still be estimated with the DR-CF algorithm.
The cut-off algorithm is described in Algorithm 1. The linear fitting is implemented in MATLAB (R2016a, The MathWorks, Inc., Natick, Massachusetts) with the least square approximation. In Algorithm 1 step 2, the size of the sliding window size is chosen to be 21 pixels for each A-line. After fitting, an error image is calculated, and then a threshold of the squared error was experimentally chosen to create a binary image in step 3. In order to determine a continuous border with the regional maximum got in step 3, a morphological closing of is used in step 4, where the octagon operator is used with a radius of 6 pixels.
Determine the DR-CF cut off border.
|1. Calculate the ;|
|2. From the backside of each A-line, estimate the attenuation with linear fitting according to Eq. (8) with a sliding window calculate the image of the fitting error (the total norm of fitting residuals) is found.|
|3. Threshold the fitting error and extract a binary image marking the regional maximum.|
|4. Morphologically close with a octagon of radius , and the backborder defined binary image is the defined cut-off depth border.|
Experiment and Implementation
In order to analyze the potential of the IVOCT-estimated optical values for distinguishing the coronary artery tissue, six types of tissues were delineated according to their pathological counterparts harvested from two cadaver hearts, resulting in 135 of IVOCT images for the analysis. The data and the implementation for the analysis are described in more detail as follows.
Coronary segments were dissected from the cadaver hearts and imaged ex vivo by optical frequency domain imaging (OFDI). The samples were pinned on a rubber sheet and immersed in a water bath with 0.9% saline solution. A 0.014-in. guidewire was introduced into the vessel followed by an OFDI catheter (LUNAWAVE, Terumo). Sequential images were acquired at a pullback rate of (, 512 A-lines per frame). The source light has a center wavelength of 1300 nm and a bandwidth of 124 nm. The depth resolution for each A-line used in this work is , which has been calculated with the refractive index stored in the raw data as 1.447.
Tissue preparation and plaque assessment
After the intracoronary imaging, samples were x-rayed and then decalcified using ethylene diamine tetraacetic acid as required. Coronary segments were sectioned at 3- to 4-mm intervals and embedded in paraffin to obtain 4- to -thick sections followed by hematoxylin and eosin (H&E) and Movat pentachrome stain.
Matching of IVOCT and histology, and differentiation and delineation of specific tissue type on histology and IVOCT
In each pathological cross section, two pathologists (CN and ST) identified the following tissue categories in histology: lipid, calcification, mixed calcified plaque, necrotic core, and macrophage. The remaining tissue was identified as fibrous plaque. The identified tissues in histology were independently delineated by the pathologists such that each region of interest (ROIs) included only one type of tissue [Figs. 1(a)–1(d)].
The histological cross section and the corresponding IVOCT cross section were matched using the lumen shape and tissue appearance in both images by two cardiologists (Y.S. and Y.O.). IVOCT cross sections with insufficient image quality and significant artifacts such as tangential artifacts were not used for the analysis. Based on the contours delineated by the pathologists in histology, the two cardiologists delineated the same tissue in matched IVOCT images in QCU-CMS (Quantitative Coronary Ultrasound-Clinical Measurement Systems, LKEB, Leiden University Medical Center, Leiden, The Netherlands) [Fig. 1(a’)–1(d’)]. When the border of the tissue could clearly be detected on IVOCT, the delineation followed the borderlines. When the border could not be clearly perceived, the contours in histological images were superimposed on the IVOCT images in order to yield a high reproducibility. If there were more than two layers of different ROIs in one cross section, only the most adluminal ROI was employed for the analysis considering the light property.
The two pathologists (G.N. and S.T.) were blinded to the IVOCT delineation and IVOCT light property analysis. At the stage of IVOCT delineation, to enable accurate colocalization of the histology with IVOCT images, two cardiologists (Y.S. and Y.O.) were unblinded to the histological identification. However, these cardiologists were blind to the following IVOCT light property analysis. The IVOCT light property analysis was performed independently from the either histological identification or the IVOCT delineation process.
The experimental flow of analyzing the IVOCT values was designed as it is shown in Fig. 2.
In the Terumo Lunawave imaging system, images are saved at 8-bit depth to allow for a high system sensitivity, resulting into a lower dynamic range than those saved at 12- to 14-bit depth.40 The algorithm has been proposed based on a transition model that is linearly related to the detected light signal. Once the image is compressed nonlinearly, the estimated values can be in different ranges. As optical properties, it is expected that the estimated attenuation coefficient and backscatter term are in the same range for the same type of tissue rather than depending on the imaging systems. In order to overcome this issue, we used a matching look-up table to match 8-bit images generated with Lunawave from Terumo (Tokyo, Japan) to 16-bit range of images generated with Illumien Optis from St. Jude Medical (St. Paul, Minnesota).41
A median filter with small window size is applied to the pixel-wise DR estimated values to remove noise. Pixel-wise DR estimation allows statistical analysis within local regions. Although it might blur the edges between different types of tissues, the median filter will hardly affect the statistical distribution within the local region.
The median filter window size was chosen as [3,11] for calculation of and the same window size was applied for noise removal. With a pixel size of , this length accordingly covers a range about . The window size 11 is chosen to be the maximum odd number that covers approximately the thickness of small identifiable structures, e.g., a thin cap defined in the literature42 ().
The data were processed and analyzed with MATLAB. For each delineated region, the distribution of intensities, attenuation, and the backscatter term are analyzed. The mean, standard deviation (std), median, maximum (max), 95th percentile (95th per.), minimum (min), skewness, and kurtosis values were calculated within each delineated region. Values between different tissue types were compared with the two-sample -test. Both three-dimensional (3-D) scatter plots of the intensity, attenuation coefficient, and backscatter term and their principal component analysis (PCA) projected two-dimensional (2-D) scatter plot are generated.
In MATLAB, the current code has been implemented with 3-D matrices calculation. Median filtering with a window size of 11 pixels on a PC with a 2.83-GHz dual-core processor currently takes 7 s for every 100 IVOCT polar images of 512 A-lines with 900 depth points. The total time of calculating the attenuation coefficient and the backscatter term with cut-off border determination is 2.7 min. Meanwhile as implemented in QCU-CMS, median filtering with a window size of 11 pixel points currently takes 24 s on a PC with a 2.66-GHz quad-core processor with hyperthreading for every 553 IVOCT polar images of 512 A-lines with 1024 depth points. The total time of calculating the attenuation coefficient and the backscatter term without cut-off border determination is 23 and 22 s, respectively. Currently, the processing time is reasonably fast and it is not a major concern. For the future clinical orientated on-the-fly application, the code and the implementation can be further optimized.
The statistical results are presented in Table 1. The average and the standard deviation values were reported for all the regions from each tissue type, and the number of analyzed regions is given under the name of the tissue. The reported IVOCT intensity is the Terumo image intensity. Using the known calibration factor of the IVOCT images, the unit of the estimated attenuation coefficient is .
|Intensity (no unit) average±std||Attenuation (mm−1) average±std||Backscatter (no unit) average±std|
The comparison results of median values, mean values, standard deviation, maximum values, 95th percentile, minimum values, skewness, and kurtosis are shown in Figs. 3Fig. 4Fig. 5Fig. 6Fig. 7Fig. 8Fig. 9–10, respectively. The whisker box plot is given with six tissue groups. When the result of the two-sample -test indicates the numbers are not significantly different, the two groups are marked with an X-mark. In the graph in the second row of each figure, the comparison of six tissues is shown as a relationship grid graph. For each paired tissue types, the square is colored as green if they are significantly different, red otherwise. On the right bottom panel, three grid graphs are combined with an “OR” operation; green indicates at least one square in the previous three graphs are marked as green, otherwise the square is marked as red. The 3-D scatter plots can be seen in Figs. 11 and 12. The PCA projected 2-D scatter plots are shown on the right-hand side of the top row for each figure.
In general, we found that the values within the delineated regions had a strong correlation with tissue types at the 5% significance level by comparing the median, mean, std, 95th percentile, and minimum values of IVOCT intensity, estimated attenuation coefficient, and the backscatter term. This can be observed from the combined comparative graphs in Figs. 3 and 4 where the red color is not presented. Also from the PCA projected 2-D plot, the other values did not show clusters as distinguishable as they are in the mean and median values. The PCA projected 2-D plot of std values shows clear clusters of fibrous and the calcified plaques. Furthermore, it is worth noting that clear cluster patterns of the calcification, fibrous, and macrophages can be observed from the 2-D plot of the skewness values.
Most of the time, the three IVOCT numbers are complementary to each other for the differentiation of the six tissue types. For example, the median values (Fig. 3) of the attenuation are not significantly different between tissue types of mixed calcified (mixed) and the fibrous but the median values of intensities and backscatter values are significantly different, and vice versa when comparing necrotic core to lipid-rich regions. This complement can be observed in all the compared figures (Figs. 3–9) but is most notable in median, mean, std, and maximum values (Figs. 3–6). In order to give an overview of the estimated values for the entire pullback, a “carpet view” image is generated. For a given stack of polar images, the maximum value for each A-line is determined in the first 1-mm depth below the lumen surface. In order to maintain the relative size of the plaques, the images are displayed with the lumen contour circumference normalized by the largest circumference over the entire pullback. The maximum projected attenuation and backscatter carpet views are shown in Figs. 13 and 14, where the largest circumference is 7.71 mm. In order to show the values in a proper window level, key points of generating the color map are appointed both experimentally and referring to the work by Gnanadesigan et al.31 As it can be seen in Fig. 15, multiple pathological images show macrophage infiltration. The corresponding locations are marked as “M” in blue, and they are much brighter than other tissues in both carpet-views. This is consistent with our quantification analysis showing that both attenuation coefficient and the backscatter term are high. In order to show the estimated values for all six types of tissues, the attenuation coefficient and the backscatter term of IVOCT images in Fig. 1 are generated as well.
This work aims at exploring the potential of DR attenuation coefficient, backscatter term, and image intensity in distinguishing different atherosclerotic tissue types with a robust implementation of the DR approach. We describe the implementation of the DR approach for the estimation of the attenuation coefficient in IVOCT images and extend it to estimate a backscatter term. To exclude the noisy regions, we describe an algorithm to automatically determine the cut-off border. Combining the estimated attenuation coefficient and the backscatter term with the image intensity, we report that it is possible to differentiate six types of clinically important tissues2,4,5 including the mixed calcified plaque, calcified plaque, fibrous plaque, lipid-rich plaque, plaque with macrophage, and plaque with necrotic core.
The light transmission model of OCT signals was used as it has been introduced in the work of Vermeer et al. The attenuating term is modeled as , while it has been modeled as in some CF approaches in the work of van Soest et al. The direct result of using this factor 2 in the exponential term is that the estimated attenuation coefficients are expected to be smaller than those previously reported in some of the studies. This difference affects only the absolute values rather than the comparative trend; therefore, we focus more on the comparative trend of the estimated values when comparing to the literature.
The attenuation coefficients of neointima have been reported as without lipid-laden and with lipid-laden in the work of Yonetsu et al.29 For the in-stent restenosis lesions of drug-eluting stents (DES) and bare-metal stents (BMS), the attenuation coefficients have been estimated as and , respectively, in the study of Nagoshi et al.30 These values have been estimated with the LightLab software. Despite the estimation algorithm being unknown, these reported values are within the similar range of estimated attenuation coefficients of the fibrous tissue (mean: , median: ). Taking the factor 2 into consideration, the estimated attenuation coefficient of the fibrous tissue (mean: , median: ) is comparable to that reported as 2 to in the work of van Soest et al.9 and relatively lower than that reported as in the work of Xu et al.8
For regions with macrophages infiltration, the range reported by van Soest et al. is , and even can be as high as in a given example, which tends to be much higher than the other tissue types. Although we found attenuation coefficients of regions with calcification and macrophages much lower than what has been reported in these two studies, the comparative relationship of the values is the same, where the value of macrophages (mean: , median: ) is the highest, and the value of the calcification (mean: , median: ) is lower than that of fibrous tissue (mean: , median: ).
For the lipid-rich regions, we found a mean attenuation value of and median . This is lower than what has been reported in the literature; 8 and 31, but higher than fibrotic tissue which fits in the reported trends. A reason for this is that the cardiologists did draw the similar shape in IVOCT as was seen in the pathological slides (Fig. 1). Due to the high attenuation, the back border of the lipidic region is not seen in IVOCT and the deeper regions are dominated by noise. To avoid the inclusion of noise in the statistical parameters, we introduced the cut-off algorithm to remove the noisy regions. In order to compare relatively high values within the superficial limited regions of the lipid plaque, a higher percentile (95th) was considered for comparison as well. It should be noted that the 95th percentile of the estimated values within the lipid-rich region can be as high as and that of macrophages can be , which is much higher than that for other tissue types. The comparative relationship in the mean, median, and 95th percentile values of the estimated attenuation coefficient is Cal < Fib < Lip < Mac.
With respect to the necrotic tissue cores, since the plaques with macrophages and lipid were labeled separately, the necrotic cores chosen in this work are expected to have fewer macrophages infiltrations. This is shown in the estimated attenuation that the necrotic core (mean , median ) is lower than lipid (mean , median ) and macrophages (mean , median ), and is just above the fibrotic tissue (mean , median ). In the work of van Soest,9 the necrotic core () has been “read” as lipid plaque when compared to the study of Xu et al. We carefully inspected the regions in the IVOCT images, which were annotated as necrotic core by pathologists, but could not identify bright regions or spots as were shown in the work of van Soest thus resulting in a lower mean and median attenuation and backscatter. Therefore, the discrepancy compared to the work of van Soest is caused by a difference in the data.
It should be noted that two clusters of values are observed in the median and mean values of the calcified plaque in the PCA projected 2-D plot (Fig. 3 and 4). This can be caused by the large variation in the median/mean values of the intensities (mean: , median: ) and the backscatter term (mean: , median: ). However, it is worth noting that the distribution of the median/mean values in the attenuation is relatively condensed (mean: , median: ), and the range of values is significantly different from any other tissue type.
As it has been stated, the DR attenuation was estimated with the precondition that the attenuation and the backscatter coefficients are linearly related to each other. Rather than estimating the backscatter coefficient, the proposed calculation of the backscatter term in our work is based on the principle that the “offset” of the linear fitting is the backscatter-related term, the exponential of which has been shown to be linearly related to the actual backscatter coefficient in the work of Schmitt et al.16 Applying this principle, the DR model was extended to estimate this term pixel wise. Our analysis of the optical model shows that this term is related to both the backscatter coefficient and the initial power of the incident light. Results show that this term provides different comparative graphs and additional information for tissue differentiation. For example, fibrotic tissue was not significantly different from macrophages for the mean intensities, but they are statistically different in mean values of the estimated backscatter term (Fig. 4).
To the best of our knowledge, only three studies have reported the backscatter coefficient in the literature of tissue analysis within vessel walls.8,29,30 The estimated values with the LightLab software have been reported for neointima as without lipid-laden and with lipid-laden,29 and have been reported as and for the in-stent restenosis lesions of DES and BMS, respectively.30 This range is comparable to that of our estimated values of the fibrous region (mean: , median: ).
The backscatter coefficient has been measured by Xu et al.8 as for fibrous tissue, for lipid-rich tissue, and for the calcified plaque. Although the absolute values are not comparable to our values due to different experimental settings and a fundamentally different approach for the estimation as mentioned before, we still observe similar patterns. The fibrous plaque has been reported to have “high backscatter and low attenuation,” while our results do show that backscatter for fibrotic tissue is relatively high with relatively low attenuation. When comparing fibrous regions to macrophage regions, it is worth noting that the mean, median, and the 95th percentile values for the former are lower, suggesting a more compact distribution. This may relate to the speckle pattern caused by the inhomogeneous distribution of the foam cells. The estimated value in calcified regions (mean: , median: ) is significantly lower than that in fibrous regions (mean: , median: ). The lipid is not significantly different from fibrous regions in the mean and median of the estimated backscatter term, but the prior is higher than the latter in the 95th percentile values (95th percentile: versus ).
Most studies in the literature focused on the comparisons of only the attenuation coefficient,9,2021.22.23.–24,28,32 or the IVOCT intensities,6,7 while a few studies combined the attenuation and the backscatter coefficients.8,29,30 The two-sample -test results suggest that the IVOCT image intensity does additionally contribute to the differentiation of the tissue types. For example, for the fibrotic versus necrotic tissue, the mean values of neither the attenuation nor the backscatter term are significantly different (Fig. 7), while the mean values of the intensities are significantly different.
On the other hand, it has been found in our previous work that the IVOCT intensities can be affected by the eccentric position of the catheter when light travels a long distance through the flush medium or enters the arterial wall with a large incident angle.43 In this analysis, this bias on intensities has been minimized by using only the lesions with front edges that are almost perpendicular to the incident light from the catheter. However, it is usually difficult to avoid this bias in clinical data. It has been shown that the intensities can be normalized with the summation of the total amount of incident light, yet additional statistical analysis on healthy tissues is required. From the principle equations for estimating attenuation coefficient and backscatter term, one may note that these two optical parameters are normalized with the summation of the total amount of intensities, and thus are expected to be less affected by this particular bias. Therefore, practical analysis with in vivo data is suggested to be focused more on these two estimated optical properties.
The catheter-related terms were not considered in this work. As far as we know, even in the most recently published studies in the literature of the experimental optical measurements, the catheter-related term was reported to be calibrated either with OCT measurements in a weakly scattering media,19 or a particular phantom.44 Since we had access only to the image data and not the used catheter, the catheter specific calibration could not be done. However, as it has been discussed in the work of Vermeer et al., the maximum bias caused by the confocal function to the estimated DR attenuation is at depth where is the focus position and the is the Rayleigh length. Calculating with used in the work of Ughi et al.,32 the maximum bias is expected to be around , which is relatively small compared to the scale for the estimated DR attenuation coefficient with IVOCT images.
Another potential issue could be that different tissue types have different refractive indices. This is not taken into account during the reconstruction process. This could lead to a difference in image resolution in different tissue types, while a fixed resolution is used. Ideally, different refractive indices should be applied for different types of tissues during the reconstruction. However, it requires prior knowledge of tissue compositions, which is almost impossible in practice.
When the resolution term is modeled as a function related to depth and tissue type, the impact on the accuracy of estimated optical coefficients can vary. Equation (4) will become8 the refractive indices of fibrous tissue, lipid pool, and calcification tend to be around 1.35, 1.43, and 1.50 with a center wavelength of 1320 nm, thus the pixel size in the previous two types of tissue should be larger than , while calcified lesions yield a smaller pixel size. Then, it is expected that estimated values can be maximally overestimated by a factor of 1.447/1.35 in normal/fibrotic tissue or underestimated by a factor of 1.447/1.50 when a large calcified lesion is present. Taking macrophages in this study as an example, the median value of estimated attenuation is . When the pixel size is considered to be the only factor that affects the estimation, the real value is expected to be within an interval of [3.00, 3.34]. Given the large difference in attenuation to the other types of tissues, the impact caused by the imaging resolution in detecting the presence of macrophages is expected to be small.
For the sake of real-time computing, a median filter is applied rather than a more complex speckle noise removal filter such as the iwhTV denoising published recently34 or the entropy filter used by Jimenez et al.45 A fixed window size has been selected as around . As it was discussed, the median filter blurs only the attenuation at the border of different tissues, and it affects the local statistical numbers minimally.
The statistical numbers for different tissue types have been calculated within the manually delineated regions. Other tissue types such as loose connective tissue, proteoglycans, calcified nodules, and smooth muscle cell-rich fibrous tissue were not assessed in this study because of the limited study sample size and the absence of specific staining.42 Artifacts would inevitably affect the light property analysis. These points can be limitations in the future for tissue classification and subjects for a further investigation.
Generating a distribution map indicating the probability of each type of plaque is one of the interesting directions in future work. Toward this goal as it has been proposed in the literature,26 the following works need to be done. Enough features should be extracted, including the first- and second-order statistical features27 of intensities, attenuation coefficient, backscatter term, the power spectrum parameters, etc. Feature selection should be used to remove strongly linear-related features, which can result in a noninvertible covariance matrix. Then, a likelihood density map can be generated for each type of tissue. For a fair comparison between tissue types to determine the most likely tissue type for each pixel, the prior probabilities are necessary to be obtained in advance. Usually, manually reviewing a pullback series is time-consuming and points of interest easily cannot be overlooked due to the huge amount of data (up to over 1000 images). Integrated with the automated lumen detection in QCU-CMS,46 carpet views as shown in Figs. 13 and 14 can be generated robust and fast without the need of predefining the segmentation of the lumen tissue. This overview of the examined coronary vessel wall can potentially improve the abnormality determination and further help with the treatment decision47,48 for coronary diseases, e.g., for stent placement and risk assessment of the patient in an automated way which is lacking in the current clinical practice.
In order to differentiate atherosclerotic tissue with IVOCT images, the DR approach has been implemented in the IVOCT images to robustly estimate the attenuation coefficient and the backscatter term. A cut-off algorithm has been proposed to efficiently exclude noisy regions. The reported values for each type of tissue are consistent with the findings in the literature. Despite assumptions made in the DR approach, statistical analysis shows that the estimated attenuation coefficient and the backscatter term can differentiate six tissue types in the lumen vessel; mixed calcified plaque, fibrous, calcification, lipid, macrophages, and necrotic core in combination with the image intensity.
Dr. Onuma reports to be a member of the advisory board of Abbott Vascular; Professor Serruys reports personal fees from Abbott, AstraZeneca, Biotronik, Cardialysis, from GLG Research, Medtronic, Sinomedical Sciences Technology, Socit Europa Digital Publishing, Stentys France, Svelte Medical Systems, Volcano, and St. Jude Medical, outside the submitted work.
Shengnan Liu received her BS degree in mathematics and applied mathematics in 2009 and her MS degree in controlling science and engineering in 2012. She is a PhD candidate at Leiden University Medical Center. Her current research interests include intensity analysis, optical properties, and tissue characterization with optical coherence tomography images.
Yohei Sotomi graduated from Osaka University and completed his cardiology training at Sakurabashi Watanabe Hospital, Japan. He is a clinical research fellow in cardiovascular medicine at the Academic Medical Center, University of Amsterdam, The Netherlands, and additionally works as a clinical researcher at Cardialysis, The Netherlands. His research interests include bioresorbable vascular scaffold, advanced intracoronary imaging (optical coherence tomography and intravascular ultrasound), and its clinical application.
Jeroen Eggermont received his MSc degree in computer science from Leiden University in 1998 and his PhD from Leiden University in 2005. He is a scientific researcher in the Division of Image Processing. Currently, his work is mainly focused on the automatic analysis of intravascular ultrasound and optical coherence tomography images.
Gaku Nakazawa is an assistant professor of the Department of Cardiology at Tokai University. He has worked as an interventional cardiologist since 2003. During his career, he has worked in CVPath Institute as a research fellow and studied vascular pathology from 2006 to 2009, which resulted in publishing many manuscripts regarding pathology of coronary arteries and drug-eluting stents.
Sho Torii is a research associate at Tokai University. He has been working at the Division of Cardiology since 2011 and his main work is cardiovascular pathology.
Takeshi Ijichi graduated from Tokai University and completed his cardiology training at the Tokai University School of Medicine. He is a clinical physician in cardiovascular medicine at Tokai University School of Medicine, Japan. His research interests include coronary imaging (optical coherence tomography and intravascular ultrasound), coronary pathology, cardiac regenerative therapy, and its clinical application.
Yoshinobu Onuma graduated from Tohoku University, Japan, and completed his cardiology fellowship at Mitsui Memorial Hospital. He is a supervising cardiologist at Cardialysis and works in the Throaxcenter at Erasmus Medical Center, The Netherlands. He was granted his doctorate in interventional cardiology at the Erasmus University in 2014. He is a leading expert in the field of bioresorbable scaffolds, metallic stents, quantitative coronary angiography, intravascular imaging, and multislice computer tomography.
Patrick W. Serruys is an emeritus professor of medicine at Erasmus Medical Center, Rotterdam, and a professor of cardiology at the Imperial College, London. He is also a Doctor Honoris Causa in medicine from the University of Athens; medicine at the Complutense University of Madrid; engineering of the Melbourne University. He has been working as a principal investigator of many clinical trials, resulting in an enormous number of publications and awards.
Boudewijn P. F. Lelieveldt is the head of the Division of Image Processing (LKEB), Department of Radiology, Leiden University Medical Center, Leiden, The Netherlands. He is also appointed as a medical delta professor at the Department of Intelligent Systems, Delft University of Technology, Delft, The Netherlands. He serves as a member of the editorial board of Medical Image Analysis and the International Journal of Cardiovascular Imaging, and is an associate editor of IEEE Transactions on Medical Imaging.
Jouke Dijkstra is an associate professor at LUMC-LKEB, where he is the leader of the vascular and molecular imaging section. He has been working at the Division of Image Processing since 1996 and is responsible for the research in vessel imaging, in particular intravascular ultrasound, optical coherence tomography, and computed tomographic angiography. The research resulted in several scientific publications and analytical software packages, which have been transferred to industry and are used world-wide.