Research Papers: Imaging

Analysis and compensation for the effect of the catheter position on image intensities in intravascular optical coherence tomography

[+] Author Affiliations
Shengnan Liu, Jeroen Eggermont, Alexander Broersen, Jouke Dijkstra

Leiden University Medical Center, Department of Radiology C2-S, Division of Image Processing, P.O. Box 9600, Leiden 2300 RC, The Netherlands

Ron Wolterbeek

Leiden University Medical Center, Department of Medical Statistics and Bioinformatics, P.O. Box 9600, Leiden 2300 RC, The Netherlands

Carol A. G. R. Busk

University of Southern Denmark, Institute of Forensic Medicine, Odense C 5000, Denmark

Helle Precht

University College Lillebaelt, Conrad Research Center, Odense SØ 5220, Denmark

Boudewijn P. F. Lelieveldt

Leiden University Medical Center, Department of Radiology C2-S, Division of Image Processing, P.O. Box 9600, Leiden 2300 RC, The Netherlands

Delft University of Technology, Intelligent Systems Department, P.O. Box 5031, Delft 2600 GA, The Netherlands

J. Biomed. Opt. 21(12), 126005 (Dec 07, 2016). doi:10.1117/1.JBO.21.12.126005
History: Received August 22, 2016; Accepted November 9, 2016
Text Size: A A A

Open Access Open Access

Abstract.  Intravascular optical coherence tomography (IVOCT) is an imaging technique that is used to analyze the underlying cause of cardiovascular disease. Because a catheter is used during imaging, the intensities can be affected by the catheter position. This work aims to analyze the effect of the catheter position on IVOCT image intensities and to propose a compensation method to minimize this effect in order to improve the visualization and the automatic analysis of IVOCT images. The effect of catheter position is modeled with respect to the distance between the catheter and the arterial wall (distance-dependent factor) and the incident angle onto the arterial wall (angle-dependent factor). A light transmission model incorporating both factors is introduced. On the basis of this model, the interaction effect of both factors is estimated with a hierarchical multivariant linear regression model. Statistical analysis shows that IVOCT intensities are significantly affected by both factors with p<0.001, as either aspect increases the intensity decreases. This effect differs for different pullbacks. The regression results were used to compensate for this effect. Experiments show that the proposed compensation method can improve the performance of the automatic bioresorbable vascular scaffold strut detection.

Figures in this Article

Cardiovascular disease (CVD) is a major cause of death worldwide.1 One of the underlying processes that cause CVD is atherosclerosis, which is the long-term accumulation of plaque in the vessel wall. The extent and composition of atherosclerosis can be visualized in vivo with intravascular optical coherence tomography (IVOCT) at a higher resolution of 10 to 20  μm24 compared to other in-vivo imaging modalities, such as intravascular ultrasound (IVUS), computed tomography angiography or magnetic resonance imaging.

IVOCT is an optical imaging modality using near-infrared (NIR) light as the imaging source. The images are acquired using a catheter, which is inserted into the coronary artery. Images of arterial cross-sections are reconstructed from the echo time delay and the intensity of backscattered light. Due to the high scattering of NIR light in blood caused by red blood cells, the artery is flushed with saline or a contrast medium to clear the blood inside the artery. The image intensity is assumed to be only tissue dependent, thus different types of tissue appear different.2

In practice, however, the signal magnitude may not be dependent only on the tissue type but also on the position of the catheter with respect to the vessel wall, which causes nontissue-related effects on the IVOCT image intensities.5,6 An example is given in Fig. 1. The average intensities are calculated within the thin superficial uniform tissue layer of a nonpathological artery segment. Nevertheless, there is a clear variation in the profile of the average intensities [Fig. 1(b)].

Graphic Jump Location
Fig. 1
F1 :

(a): An IVOCT image of a nonpathological artery wall: the artery has a regular and almost circular shape; three arterial layers, intima (I), media (M), and adventitia (A), are clearly visible (as shown in zoomed-in top-left corner). (b): Polar representation of the image in (a) sampled clockwise along radial A-lines from the catheter center shown as a bright line on the top of the image. For each A-scan, the average intensity within a superficial thin layer (50  μm) is calculated and shown as the green profile. The white curve is the smoothed green profile.

The importance of analyzing the effect on intensities caused by the position of the catheter has been well depicted in the field of IVUS. Courtney et al.7 showed that the IVUS image intensities are significantly related to the catheter position. Their study concluded that when the distance or the angle toward the luminal wall increases, the intensity will decrease for both intima-media tissue and adventitia tissue. Earlier work810 shows that the reflected ultrasound signal is critically dependent on the angle of incidence and varies for different types of arterial plaques.

In the literature of IVOCT image analysis, statistical values of the intensities are commonly used as key features for both automated detection algorithms and the quantitative studies. For example, mean intensity has been applied as one of the textural features for automated tissue characterization.11 Recently, a stent strut detection algorithm has been proposed by training a supervised artificial neural network classifier with statistical features including the maximum, mean, median intensities, and so on.12 Furthermore, percentile values of intensity distributions are often used as thresholds. For example, the fifth percentile has been used as a threshold for noise removal.1316 Other percentiles were used as cutoff values to determine the trailing shadows15,17 for metal strut detection, or the black core regions18 in the detection of the bioresorbable vascular scaffold (BVS) strut. With the assistance of the BVS strut detection, median and peak intensity values within the black core regions were quantitatively analyzed to track the resorption of the BVS struts over time in IVOCT images at 6, 12, 24, and 36 month poststenting, respectively.19

The distributions of the intensities can appear different due to the catheter position, which may increase the variation of those statistical numbers. To the best of our knowledge, only one study with respect to this issue on OCT image intensities has been reported which was about nonperpendicular incident light causing a significant variance in the measurement of the articular cartilage.6 In the followup studies of the bioresorption progress of the BVS strut, the bias in light intensity caused by the eccentric catheter was claimed to be minimized with normalization by manual selecting reference regions.20 Results from IVUS cannot be applied directly to IVOCT, due to the differences in physical properties between both modalities.

The aim of this work is to analyze the effect of the catheter position, with regards to both the distance to the vessel wall and the incident angle of light, on IVOCT image intensities. Based on this analysis, a compensation algorithm is proposed to reduce this effect. As an application of compensated images, images with foam cells have been enhanced and compared with histological cross sections. Furthermore, the compensation algorithm is used in combination with an existing BVS detection algorithm.

The general structure of the paper is as follows. In Sec. 2.1, a light transmission model incorporating both distance-dependent and angle-dependent factors is introduced. In Sec. 2.2, a hierarchical multivariant linear regression model is proposed to further investigate the relationship and estimate both factors. The regression result is further used in Sec. 2.3 to propose a method to compensate images. Results are presented in Sec. 3. Furthermore, a BVS struts detection experiment with the compensated images was carried out in Sec. 4. All the experiments and results are discussed in Sec. 5 with limitations and future works given as well. Conclusions are drawn in Sec. 6.

Images of nonpathological segments from nine IVOCT pullbacks recorded with a C7XR swept-source OCT system and a C7 Dragonfly Imaging Catheter (St. Jude Medical, Minnesota) were used. The technical details are listed in Tables 1 and 2.

Table Grahic Jump Location
Table 1Technical details of the IVOCT system.
Table Grahic Jump Location
Table 2Number of selected frames in each pullback.
Distance and Incident Angle Extended Light Transmission Model

A schematic overview of light propagation for IVOCT imaging is shown in Fig. 2. The light emitted from the catheter first travels through the flush medium before reaching the arterial wall with a distance denoted as xt. At the interface between the flush medium and the arterial wall, both reflection and refraction occur. θ represents the incident angle of the light entering the arterial wall. Δx represents the light transmitting distance of the refracted light beam inside the arterial tissue. For the convenience of explanation, we introduce x=xt+Δx.

Graphic Jump Location
Fig. 2
F2 :

Light transmission. xt denotes the distance between the light source and the arterial wall, θ is the incident angle of the light beam.

Light transmission model

As the light propagates inside the arterial wall, the intensity of an OCT signal is typically modeled as the first-order scattering function of x and Δx as21Display Formula

Ib(x)12IinμbT(x)·e2μtΔx,(1)
where Iin denotes the light intensity upon entering the arterial wall, Ib(x) denotes the backscattered light intensity from the distance x, μb represents the backscattering coefficient, and μt is the total attenuation coefficient (summation of scatter and absorption). T(x) is the confocal function, which is defined as follows:22Display Formula
T(x)=[(xz0zR)2+1]1/2,(2)
where z0 and zR are the beam waist and the Rayleigh length, respectively.

The intensity entering the luminal wall is affected by two factors: (1) the attenuation in the flush medium region (FMR) and (2) the reflection and refraction at the interface of flush medium and the arterial wall. In a well-flushed artery, the FMR region can be regarded as homogeneous, nonscattering and weakly attenuating, obeying Lambert–Beer law.23 With a constant attenuation, μf, the light decay is determined by the distance from the catheter to the lumen wall, xt.

The interaction of the light is more complex at the interface between FMR and the lumen tissue due to the relative irregular surface of the arterial wall. To analyze the effect of the incident angle on image intensities, the total effect of the incident angle is normalized into 0 to 1 by using a term similar to the Fresnel transmission ratio. Thus, Iin is modeled as Display Formula

IinI0·Tr(θ,ni,nt)β2·eμfxt,(3)
where β2 is the parameter to be estimated. Tr(θ,ni,nt) is the Fresnel like function, which is calculated with the incident angle θ, index of refraction of the incident medium ni and transmission medium nt, respectively.23 With Eq. (3) substituted in Eq. (1) and taking the logarithm from both sides results in Display Formula
lnIb(x)μfxt+βlnTr(θ,ni,nt)+lnT(xt+Δx)2μtΔx+C(I0,μb),(4)
where C(I0,μb)ln(I0·μb) is a constant term.

Parameter Estimation of the Linear Model with Hierarchical Linear Regression
Hierarchical linear regression

Hierarchical linear models are specifically utilized for data with hierarchical structures.24 Here, a hierarchical linear model is designed to analyze the potential relationship between OCT image intensities and three factors: distance, (x); angle, lnTr(θ,ni,nt); and the constant term C(I0,μb). The linear model for regression is Display Formula

lnIb(x)=β0+β1·x+β2·lnTr(θ,ni,nt).(5)
In order to keep the consistency of the notations, the parameters were denoted as β0, β1, and β2. The A-lines can be hierarchized into different frames, which in turn can be hierarchized into different pullbacks. Based on this observation, a three-level linear model is considered for this study (see Fig. 3).

Graphic Jump Location
Fig. 3
F3 :

Multilevel linear model.

Implementation

The lumen border in the Cartesian images was used to estimate the incident angle. To compensate for the polar to Cartesian transformation, the lumen border points were resampled with respect to the depth. The angle was calculated in a window of nine points.

The index of refraction of the flush solution is 1.449 (read from the stored data). The refraction index of intima is about 1.358.25 Therefore, the incident angle is the only variable during the calculation of the transmission ratio for each point.

Intensities of only a thin inner layer of the arterial wall are used for the statistical analysis, then Δx0, and thus x=xt+Δxxt. The general trend of the signal regards to the distance xt is decreasing due to both the attenuation of the flush medium and the confocal function. Approximating this term as linear, the object model for hierarchical linear regression can be written as Display Formula

lnIb(x)=C(I0,μb)+β1·xt+β2lnTr(θ).(6)
This can be equalized to the hierarchical linear regression model if we denote β0=C(I0,μb); thus, the linear regression can be used to investigate the linear relationship regarding the distance and the incident angle.

Compensation

The linear model that describes the effect of the catheter position can also be used for the compensation of this effect.

Based on the linear regression model, the primary goal for the compensation is to normalize the IVOCT image intensities within the superficial layer of the nonpathological artery. This can be achieved with the following equation involving the regression result: Display Formula

Icompensated(y)=Ioriginal(y)·eβ^0Ib(x).(7)
As defined, Ioriginal(y) and Icompensated(y) denote the original and the compensated IVOCT signals at the depth y. β0 is the estimated constant term in the regression model. With a thin layer with thickness Δx selected, Ib(x) is the average intensity within the superficial thin layer: Display Formula
Ib(x)=y=xtxt+ΔxIoriginal(y)/Δx.(8)

Noting that the following mathematical equation holds Display Formula

Ib(x)β^1·xIb(t)dt,(9)
where β^1 is the estimated parameter. The intensities can be normalized as follows: Display Formula
Icompensated(y)=Ioriginal(y)·eβ^0β^1·xIb(t)dtIoriginal(y)·eβ^0β^1·xMIb(t),(10)
where M is a large depth selected far enough away from the lumen border. Noting that Display Formula
y=xMIcompensated=eβ^0β^1=constant,(11)
the principle of this compensation method is to normalize the summation of intensities behind the lumen border. Since arterial tissues are strong scattering and weak absorbing, the summation of the IVOCT intensities should be approximately constant (linearly related to the total emitting energy from the catheter) for most arterial tissue types; thus, this method is not limited to compensate only the IVOCT images of nonpathological arteries.

Hierarchical Linear Regression

The hierarchical linear regression considers three fixed effects and two random effects. The F-tests result for each of the fixed effects specified in the model indicate that all three effects contribute to the model statistically significantly with p<0.001.

Table 3 shows the results of the fixed model. It was found that the constant related to light source (β0=6.5121, SE=0.1615, p<0.001), the distance between catheter and artery wall (β1=0.0023, SE=0.0000, p<0.001) and the logarithm of “Fresnel” transmission ratio (β2=2.8178, SE=0.0725, p<0.001) were significant predictors.

Table Grahic Jump Location
Table 3Estimates of fixed effects.003
Table Footer NoteaDependent variable: natural logarithm of intensities.

Table 4 shows the results for the two random factors and residual covariance matrices. Results indicate that both the defined random effects—the frame number (Nu=0.0060, SE=0.0007, p<0.001) and pull-back number (Nu=0.2342, SE=0.1173, p<0.05)—contribute to the covariance statistically significantly with almost two thirds of the total variance. However, the influence of the frame number is relatively very small (1.6%) compared to the other contributors. Based on this observation, this random effect can be ignored during modeling.

Table Grahic Jump Location
Table 4Estimates of covariance parameters.004
Table Footer NoteaDependent variable: natural logarithm of intensities.

The histogram of the residual (ϵ=0.1276, SE=0.0006, p<0.0001) is distributed symmetrically around zero with a mean value of 2.01E11 and a standard deviation of 0.3568, thus indicating the model can fit the data well.

Compensation

With the results of the linear regression, images can be compensated using Eq. (10) proposed in Sec. 2.3. Figure 4 demonstrates the compensation of IVOCT images of nonpathological arteries. The nonuniform image intensities behind the lumen border and even a small shadow artifact on the lumen wall were compensated. The nonuniform image intensities and even a small unusual shadow on the lumen wall were compensated.

Graphic Jump Location
Fig. 4
F4 :

(a and c) The images before compensation and (b and d) the compensated images. The image in (a) is a clear example of the effects caused by the eccentric position of the catheter, the compensated image (b) shows more evenly distributed intensities around the artery wall. In the second example (c), the light intensities are more or less homogeneous but there is information missing at 1 o’clock. (d) After compensation, the shadow region is removed.

Regarding the application of the compensation algorithm, two experiments were carried out. In the first experiment, compensated IVOCT images with foam cells are compared with the original images assisted with the histological images. In the second experiment, the compensation algorithm was applied to IVOCT images with BVS struts, followed with the strut detection and the evaluation.

Foam Cells Visualization

The compensated images are compared with histological cross sections. Selected frames from two ex-vivo OCT pullbacks on explanted hearts26 were used. OCT imaging was performed with the Ilumien PCI Optimisation system and a C7 Dragonfly Imaging Catheter of LightLab Imaging, St. Jude Medical, Minnesota. The proximal 5 cm of the vessels were cut out and standard paraffin embedding was performed. For every 200  μm, 3-μm-thick sections were cut and stained with haematoxylin–eosin. These slices were annotated by a pathologist and matched with the corresponding OCT frames based on anatomical landmarks.

Two examples with bright spots are given in Fig. 5. In Figs 5(a)5(c), the image intensities marked with red arrows were darkened due to the residual of blood within lumen. In the compensated image, the darkened regions were compensated and the edges between the calcified region and the fibrous region are more clear. In Figs. 5(d) and 5(e), the image regions near side branches were darkened due to the eccentric catheter position. The foam cells, marked with the red arrows, are more accentuated in the compensated image than in the original image.

Graphic Jump Location
Fig. 5
F5 :

Demonstration of the effect of compensation. (a–f) From right to left are histology images, original IVOCT images, and compensated images. In the histology images, red arrow indicates the location with form cells, and stars mark the calcified lesions as landmarks. The original and the compensated images are displayed at the same contrast and brightness level, and the red arrows mark noticeable regions for comparison.

Bioresorbable Vascular Scaffold Strut Detection

In order to examine the impact of compensating images on automated image segmentation, it is tested on the BVS strut detection, as proposed by Wang et al.18 For this purpose, eight pullbacks were used, which were acquired with a C7XR swept-source OCT system and a C7 Dragonfly Imaging Catheter (St. Jude Medical, Minnesota) at 6 to 12 months poststenting. All the stents are the ABSORB 1.1 BVS (Abbott Vascular, Santa Clara, California). The manual drawn ground truth (GT) data contains 7933 black cores in total. The experiment is carried out with the in-house developed software QCU-CMS (LUMC, Leiden, The Netherlands).

The compensated pullbacks were rescaled linearly by aligning the 99.5 percentile value of the histogram to that of the original pullbacks, thus the same detection scheme can be applied, as it is described in the work of Wang et al.18 The results were evaluated by counting the true positive (TP), the false positive (FP), and false negative (FN), and the F-score was calculated as the measurement of the detecting performance.27

With the described data and experimental settings, the compensated images were used for automatic BVS strut detection. The detection results can be seen in Table 5, where the outperformed F-score is marked with bold font. The compensation algorithm improves the BVS struts detection for 6 out of 8 pullbacks by 1.0 to 6.6 percent in F-score. An example of an improved BVS detection image can be seen in Fig. 6. In QCU-CMS, all the delineations are displayed in the original image. In order to make a clear distinction, the delineations of the GT, the detection results with the original images, and the detection results with the compensated images are shown separated in A, B, and D. The box-and-whisker plot in C shows the intensity percentile within 10 BVS black core regions.

Table Grahic Jump Location
Table 5The stent struts detection results.
Graphic Jump Location
Fig. 6
F6 :

An example of image with improved BVS detection. (a) The cross section with GT BVS delineated with cyan color. (b) Detection results with the original image delineated with white color. (d) Detection results with the compensated image delineated with white color. (c) Within each black core of the GT, the percentile of both original and compensated image can be seen.

The aim of this work is to investigate the influence of the position of the catheter on IVOCT intensities and to use this knowledge to compensate for it. Two aspects of the catheter’s position were analyzed: (1) the distance between the catheter and the arterial wall and (2) the angle of incidence of the light entering the arterial wall.

Hierarchical Linear Regression

The statistic analysis focuses on intima-media regions on the artery wall, thus the tissue-dependent effects were minimized in the study. Similar region selection criteria have been used by Courtney on IVUS images.7 The statistical results show that the amount of light that enters the artery wall is significantly related to the catheter position. In the linear model, an angle-related transmission ratio has been used to model the trend. The trend of this transmission ratio is conforming to empirical observations.

As the angle of incidence increases, the IVOCT intensities decrease accordingly. The more the angle of incidence approaches to a critical angle, the faster the IVOCT intensities decrease. When the angle of incidence approaches this critical angle, light propagation into the tissue in the artery wall is limited. This can explain the “signal dropout” reported by van Soest et al.5 When the incident angle becomes equal or even larger than a critical angle, there is hardly light entering the artery wall, and the transmission ratio will approximate to zero. This results in the appearance of disconnecting tissue along the arterial wall, which has been reported as dissection artifacts in an IVUS study by Mario et al.10 Due to the lack of valid signal, these artifacts cannot be compensated.

As random variables, the pullback number and frame number contribute to the hierarchical linear model significantly in terms of covariance. The covariance contribution of frame number is relatively small enough to be ignored. The covariance of the pullback number occupies almost two thirds of the total covariance. Since the same flush medium was used, this suggests that the distance-dependency can differ between IVOCT imaging catheters. This confirms the statement of van Soest et al. that the parameters of catheters differ from each other.22

The range of the estimated incident angles in the experiment is relatively small due to the elliptical shape of the artery wall, which is an inevitable issue of in-vivo IVOCT data. Since the angle related term is between 0 and 1, the logarithm operation can enlarge the range of the transmission term in the linear model. Another potential issue related to the angle estimation is that the angle of incidence has been estimated with 2-D IVOCT images: the best estimation that can be achieved at present. The estimation can approximate the spatial angle well because the imaging catheter has elasticity to resist over-bending and can thus be assumed parallel to the longitudinal direction of the artery. Nevertheless, it would be interesting to measure the angle of incidence in 3-D in a future study.

Compensation

The compensation of the intensities in the OCT images can enhance the visualization of arterial tissues in IVOCT images. An eccentric position of the catheter can result in inhomogeneous intensity in homogeneous tissue, which requires the contrast and brightness levels to be constantly adjusted during visual inspection. Our proposed compensation algorithm improves the visualization by balancing the signal levels within each pullback, which can be seen in both in-vivo (Fig. 4) and ex-vivo images (Fig. 5).

The principle of the algorithm is that for each pixel at the lumen boundary, the same amount of light enters the tissue, such that further analysis can be carried out without the bias caused by the catheter position. It is designed to overall enhance the absolute intensities for each A-line rather than changing the relative trend. Therefore, it will not affect the parameters like attenuation and backscatter. This is why the proposed algorithm preserves the dark trend within regions with weak backscattering (calcified lesion, dark square inside BVS stent struts, etc) or behind tissue with high attenuation (foam cells, e.g.).

Noting Figs. 5(b)5(d), shadow artifacts caused by the residual blood in the protective sheath were compensated as well. This is because our compensation algorithm compensates the total energy behind the lumen border for each A-line. The shadow artifacts caused by factors within the lumen, e.g., residual blood and thrombosis, may be the result of a sudden drop of the total energy compared to A-lines in the neighborhood. Therefore, normalizing the total energy can compensate such local shadow artifacts as well.

In the BVS strut detection experiment, the overall result shows that the compensation algorithm can improve the performance of the BVS struts detection. As it is shown in Fig. 6(c), the box-and-whisker plot of the percentile values within the black cores, we can observe that the percentiles in the compensated image are lower than those in the original images. Furthermore, it is worth noting that the lower percentiles over all the black cores are more condensly distributed. This can be the main reason that the proposed algorithm improves detection, since the lower percentiles are used as thresholds in the BVS strut detection.18

Meanwhile, there are two pullbacks with lower F-score, pullback 1 and 8. The first pullback was from a patient with a previously implanted metal stent strut followed by a BVS stent treatment, which results in multiple artificial layer structures. This also explains the large amount of the FP in the detection results on the original images in Table 5. The FP ratio on the compensated images is even higher, and this may be because the deeper layer structures are enhanced. Despite a higher FP ratio, the detection with the compensated images gives a higher TP rate. In the eighth pullback, the F-score with the compensated image is slightly lower than that with the original images, but it should be noted that both F-scores are already very high.

Limitations

There are two catheter-related terms for modeling the IVOCT signal. One is the confocal function and the other is the spectral coherence term. Both terms can be different for different catheters. This can be the reason for the variation that lies in the estimated intercept for the hierarchical regression results. As a result, the compensation factor can be different as well. Noting that the compensation factor serves mainly to normalize the total energy of the incident light onto the arterial tissue, this difference causes only a shift of the histogram of the compensated pullback rather than the shape of the histogram. In these situations, it would be sufficient to apply the compensation algorithm with a fixed “average” factor to remove the bias caused by the catheter position. For the quantification analysis on the data from different pullbacks,19,20 the compensating images with fixed factor can be still useful to remove the bias caused by the catheter position, but for further comparison of the statistical numbers, the differences caused by different catheter parameters should be taken into account using statistical tools, such as hierarchical models and so on. It can be useful to characterize the variance caused by different catheters and model it into the compensation method. However, in practice, this will be difficult since most hospitals will not do additional calibration measurements of the catheter before or after the procedure.

Another limitation regarding the catheter-related parameters is that the overall distance-dependent decreasing trend is implemented as a linear regression model after taking the logarithm. So, the algorithm is designed to compensate the general descending trend rather than the Gaussian-shape variation in the confocal function. Taking the numbers reported in the work of Ughi et al.11 as an example, z01.5  mm (focal point) and zR2  mm (Rayleigh length), the confocal function T(x) at depth x=4  mm equals to 0.625. After the overall linear decreasing trend is compensated, the variation of this term should be around 1 with a variation of ±18%. That is to say, the total incident light is ether overestimated or underestimated by about 18%. For a full depth of 2 mm, the variation of intensities of each pixel can be 0.045%, which is relatively small.

Future Directions

In this paper, the compensation algorithm shows its potential to be a preprocessing step for automatic BVS stent strut detection. For further incorporating the algorithm into the standard BVS detection workflow, more validation studies are needed, thus the parameters can be further optimized with the compensated images. Meanwhile, it is interesting to see if the compensation algorithm can improve the automatic detection of metal stent struts15 as well. Further experiments will be carried out for the metal stent strut detection in the future.

Both aspects regarding the catheter position, the distance from the catheter to the artery wall and the angle of light incident upon the artery wall, significantly affect IVOCT image intensities. The hierarchical linear regression result shows that as either aspect increases the intensity decreases. Using the hierarchical linear regression result, a compensation algorithm is proposed to enhance the IVOCT images and improve the BVS strut detection. Using the pathological images as a gold standard, it is illustrated that the foam cells appear more accentuated in the compensated image than in the original image. Quantitative evaluation results of the BVS strut detection with the compensated images show that the proposed method improves the performance for the BVS struts detection.

There are no competing interests to disclose.

Support to Shengnan Liu provided by the Chinese Scholarship Council is gratefully acknowledged.

Mozaffarian  D.  et al., “Heart disease and stroke statistics—2016 update,” Circulation. 133, (15 ), e599  (2016).CrossRef
Tearney  G. J.  et al., “Consensus standards for acquisition, measurement, and reporting of intravascular optical coherence tomography studies: a report from the international working group for intravascular optical coherence tomography standardization and validation,” J. Am. Coll. Cardiol.. 59, , 1058 –1072 (2012). 0735-1097 CrossRef
Prati  F.  et al., “OCT-based diagnosis and management of STEMI associated with intact fibrous cap,” JACC Cardiovasc. Imaging. 6, , 283 –287 (2013).CrossRef
Evelyn  R., , van Leeuwen  A., and Serruys  P. W., Eds., Optical Coherence Tomography in Cardiovascular Research. , 1st ed.,  Informa Healthcare ,  Oxon  (2007).
van Soest  G.  et al., “Pitfalls in plaque characterization by OCT: imageartifacts in native coronary arteries,” JACC Cardiovasc. Imaging. 4, (7 ), 810 –813 (2011).CrossRef
Huang  Y.-P.  et al., “Effects of optical beam angle on quantitative optical coherence tomography (OCT) in normal and surface degenerated bovine articular cartilage,” Phys. Med. Biol.. 56, (2 ), 491 –509 (2011). 0031-9155 CrossRef
Courtney  B. K.  et al., “Effects of transducer position on backscattered intensity in coronary arteries,” Ultrasound Med. Biol.. 28, (1 ), 81 –91 (2002).CrossRef
Picano  E.  et al., “Angle dependence of ultrasonic backscatter in arterial tissues: a study in vitro,” Circulation. 72, (3 ), 572 –576 (1985). 0009-7322 CrossRef
de Kroon  M.  et al., “Angle-dependent backscatter from the arterial wall,” Ultrasound Med. Biol.. 17, (2 ), 121 –126 (1991).CrossRef
Mario  C. D.  et al., “The angle of incidence of the ultrasonic beam: a critical factor for the image quality in intravascular ultrasonography,” Am. Heart J.. 125, (2 ), 442 –448 (1993).CrossRef
Ughi  G. J.  et al., “Automated tissue characterization of in vivo atherosclerotic plaques by intravascular optical coherence tomography images,” Biomed. Opt. Express. 4, , 1014 –1030 (2013). 2156-7085 CrossRef
Nam  H. S.  et al., “Automated detection of vessel lumen and stent struts in intravascular optical coherence tomography to evaluate stent apposition and neointimal coverage,” Med. Phys.. 43, , 1662 –1675 (2016). 0094-2405 CrossRef
Menguy  P.-Y.  et al., “Preliminary results for the supervised detection of lumen and stent from OCT pullbacks,” IRBM. 37, (5–6 ), 271 –275 (2016).CrossRef
Ughi  G. J.  et al., “Automatic segmentation of in-vivo intra-coronary optical coherence tomography images to assess stent strut apposition and coverage,” Int. J. Cardiovasc. Imaging. 28, , 229 –241 (2012).CrossRef
Wang  A.  et al., “Automatic stent strut detection in intravascular optical coherence tomographic pullback runs,” Int. J. Cardiovasc. Imaging. 29, , 29 –38 (2013).CrossRef
Wang  A., and Tang  X., “Automatic segmentation of vessel lumen in intravascular optical coherence tomography images,” in  IEEE Int. Conf. on Mechatronics and Automation , pp. 948 –953 (2016).CrossRef
Dubuisson  F.  et al., “Automated peroperative assessment of stents apposition from OCT pullbacks,” Comput. Biol. Med.. 59, , 98 –105 (2015). 0010-4825 CrossRef
Wang  A.  et al., “Automatic detection of bioresorbable vascular scaffold struts in intravascular optical coherence tomography pullback runs,” Biomed. Opt. Express. 5, , 3589 –3602 (2014). 2156-7085 CrossRef
Nakatani  S.  et al., “Temporal evolution of strut light intensity after implantation of bioresorbable polymeric intracoronary scaffolds in the absorb cohort b trial,” Circ. J.. 78, (8 ), 1873 –1881 (2014).CrossRef
Nakatani  S.  et al., “Bioresorption and vessel wall integration of a fully bioresorbable polymeric everolimus-eluting scaffold: optical coherence tomography, intravascular ultrasound, and histological study in a porcine model with 4-year follow-up,” JACC: Cardiovasc. Interventions. 9, , 838 –851 (2016).CrossRef
Valery  T., Tissue Optics: Light Scattering Methods and Instruments for Medical Diagnosis. , 2nd ed., Vol. 23, ,  SPIE ,  Bellingham, Washington  (2007).
van Soest  G.  et al., “Atherosclerotic tissue characterization in vivo by optical coherence tomography attenuation imaging,” J. Biomed. Opt.. 15, (1 ), 011105  (2010). 1083-3668 CrossRef
Pedrotti  F. L., and Pedrotti  L. S., Introduction to Optics. , 2nd ed.,  Prentice-Hall International ,  Englewood Cliffs  (1993).
Sullivan  L. M., , Dukes  K. A., and Losina  E., “Tutorial in biostatistics. an introduction to hierarchical linear modelling,” Stat. Med.. 18, , 855 –888 (1999). 1097-0258 CrossRef
van der Meer  F. J., “Vascular applications of quantitative optical coherence tomography,” Dissertation, Universiteit van Amsterdam (2005).
Precht  H.  et al., “Optimisation of post mortem cardiac computed tomography compared to optical coherence tomography and histopathology,” J. Forensic Radiol. Imaging. 2, (2 ), 85 –90 (2014).CrossRef
Goutte  C., and Gaussier  E., “A probabilistic interpretation of precision, recall and f-score, with implication for evaluation,” in  European Conf. on Information Retrieval , pp. 345 –359,  Springer ,  Berlin Heidelberg  (2005).CrossRef

Shengnan Liu is a PhD candidate at Leiden University Medical Center. She received her BSc degrees in mathematics and applied mathematics in 2009 and MSc degrees in controlling science and engineering in 2012. Her current research interests include intensity analysis, optical properties, and tissue characterization with optical coherence tomography images.

Jeroen Eggermont is a scientific researcher in the Division of Image Processing. He received his MSc degree in computer science from Leiden University in 1998. In 2005, he received his PhD from Leiden University. Currently, his work is mainly focused on the automatic analysis of intravascular ultrasound and optical coherence tomography images.

Ron Wolterbeek is a senior researcher at Department of Medical Statistic and Bioinformatics, Leiden University. His research of interest lies in medical statistics.

Alexander Broersen is a scientific researcher in the Division of Image Processing. He received his MSc degree in computer science at the University of Twente in 2003. He defended his PhD at Eindhoven University of Technology. In 2008, he started as a postdoctoral researcher at Leiden University Medical Center. His work involves the automatic diagnostic vascular analysis by comparing computed tomography angiography with corresponding intravascular ultrasound and IVOCT datasets.

Carol A. G. R. Busk is an undergraduate student in the Forensic Institute at the University of Southern Denmark since February 2015.

Helle Precht is an associate professor in radiography at the University College Lillebaelt. She is particularly interested in digital radiography, radiation protection and imaging quality optimization. She participates in national and international research at the European level.

Boudewijn P.F. Lelieveldt is heading the Divison of Image Processing (LKEB) at the Department of Radiology, Leiden University Medical Center, Leiden, The Netherlands. He is also appointed as 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 worldwide.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation

Shengnan Liu ; Jeroen Eggermont ; Ron Wolterbeek ; Alexander Broersen ; Carol A. G. R. Busk, et al.
"Analysis and compensation for the effect of the catheter position on image intensities in intravascular optical coherence tomography", J. Biomed. Opt. 21(12), 126005 (Dec 07, 2016). ; http://dx.doi.org/10.1117/1.JBO.21.12.126005


Figures

Graphic Jump Location
Fig. 1
F1 :

(a): An IVOCT image of a nonpathological artery wall: the artery has a regular and almost circular shape; three arterial layers, intima (I), media (M), and adventitia (A), are clearly visible (as shown in zoomed-in top-left corner). (b): Polar representation of the image in (a) sampled clockwise along radial A-lines from the catheter center shown as a bright line on the top of the image. For each A-scan, the average intensity within a superficial thin layer (50  μm) is calculated and shown as the green profile. The white curve is the smoothed green profile.

Graphic Jump Location
Fig. 2
F2 :

Light transmission. xt denotes the distance between the light source and the arterial wall, θ is the incident angle of the light beam.

Graphic Jump Location
Fig. 3
F3 :

Multilevel linear model.

Graphic Jump Location
Fig. 4
F4 :

(a and c) The images before compensation and (b and d) the compensated images. The image in (a) is a clear example of the effects caused by the eccentric position of the catheter, the compensated image (b) shows more evenly distributed intensities around the artery wall. In the second example (c), the light intensities are more or less homogeneous but there is information missing at 1 o’clock. (d) After compensation, the shadow region is removed.

Graphic Jump Location
Fig. 5
F5 :

Demonstration of the effect of compensation. (a–f) From right to left are histology images, original IVOCT images, and compensated images. In the histology images, red arrow indicates the location with form cells, and stars mark the calcified lesions as landmarks. The original and the compensated images are displayed at the same contrast and brightness level, and the red arrows mark noticeable regions for comparison.

Graphic Jump Location
Fig. 6
F6 :

An example of image with improved BVS detection. (a) The cross section with GT BVS delineated with cyan color. (b) Detection results with the original image delineated with white color. (d) Detection results with the compensated image delineated with white color. (c) Within each black core of the GT, the percentile of both original and compensated image can be seen.

Tables

Table Grahic Jump Location
Table 1Technical details of the IVOCT system.
Table Grahic Jump Location
Table 2Number of selected frames in each pullback.
Table Grahic Jump Location
Table 3Estimates of fixed effects.003
Table Footer NoteaDependent variable: natural logarithm of intensities.
Table Grahic Jump Location
Table 4Estimates of covariance parameters.004
Table Footer NoteaDependent variable: natural logarithm of intensities.
Table Grahic Jump Location
Table 5The stent struts detection results.

References

Mozaffarian  D.  et al., “Heart disease and stroke statistics—2016 update,” Circulation. 133, (15 ), e599  (2016).CrossRef
Tearney  G. J.  et al., “Consensus standards for acquisition, measurement, and reporting of intravascular optical coherence tomography studies: a report from the international working group for intravascular optical coherence tomography standardization and validation,” J. Am. Coll. Cardiol.. 59, , 1058 –1072 (2012). 0735-1097 CrossRef
Prati  F.  et al., “OCT-based diagnosis and management of STEMI associated with intact fibrous cap,” JACC Cardiovasc. Imaging. 6, , 283 –287 (2013).CrossRef
Evelyn  R., , van Leeuwen  A., and Serruys  P. W., Eds., Optical Coherence Tomography in Cardiovascular Research. , 1st ed.,  Informa Healthcare ,  Oxon  (2007).
van Soest  G.  et al., “Pitfalls in plaque characterization by OCT: imageartifacts in native coronary arteries,” JACC Cardiovasc. Imaging. 4, (7 ), 810 –813 (2011).CrossRef
Huang  Y.-P.  et al., “Effects of optical beam angle on quantitative optical coherence tomography (OCT) in normal and surface degenerated bovine articular cartilage,” Phys. Med. Biol.. 56, (2 ), 491 –509 (2011). 0031-9155 CrossRef
Courtney  B. K.  et al., “Effects of transducer position on backscattered intensity in coronary arteries,” Ultrasound Med. Biol.. 28, (1 ), 81 –91 (2002).CrossRef
Picano  E.  et al., “Angle dependence of ultrasonic backscatter in arterial tissues: a study in vitro,” Circulation. 72, (3 ), 572 –576 (1985). 0009-7322 CrossRef
de Kroon  M.  et al., “Angle-dependent backscatter from the arterial wall,” Ultrasound Med. Biol.. 17, (2 ), 121 –126 (1991).CrossRef
Mario  C. D.  et al., “The angle of incidence of the ultrasonic beam: a critical factor for the image quality in intravascular ultrasonography,” Am. Heart J.. 125, (2 ), 442 –448 (1993).CrossRef
Ughi  G. J.  et al., “Automated tissue characterization of in vivo atherosclerotic plaques by intravascular optical coherence tomography images,” Biomed. Opt. Express. 4, , 1014 –1030 (2013). 2156-7085 CrossRef
Nam  H. S.  et al., “Automated detection of vessel lumen and stent struts in intravascular optical coherence tomography to evaluate stent apposition and neointimal coverage,” Med. Phys.. 43, , 1662 –1675 (2016). 0094-2405 CrossRef
Menguy  P.-Y.  et al., “Preliminary results for the supervised detection of lumen and stent from OCT pullbacks,” IRBM. 37, (5–6 ), 271 –275 (2016).CrossRef
Ughi  G. J.  et al., “Automatic segmentation of in-vivo intra-coronary optical coherence tomography images to assess stent strut apposition and coverage,” Int. J. Cardiovasc. Imaging. 28, , 229 –241 (2012).CrossRef
Wang  A.  et al., “Automatic stent strut detection in intravascular optical coherence tomographic pullback runs,” Int. J. Cardiovasc. Imaging. 29, , 29 –38 (2013).CrossRef
Wang  A., and Tang  X., “Automatic segmentation of vessel lumen in intravascular optical coherence tomography images,” in  IEEE Int. Conf. on Mechatronics and Automation , pp. 948 –953 (2016).CrossRef
Dubuisson  F.  et al., “Automated peroperative assessment of stents apposition from OCT pullbacks,” Comput. Biol. Med.. 59, , 98 –105 (2015). 0010-4825 CrossRef
Wang  A.  et al., “Automatic detection of bioresorbable vascular scaffold struts in intravascular optical coherence tomography pullback runs,” Biomed. Opt. Express. 5, , 3589 –3602 (2014). 2156-7085 CrossRef
Nakatani  S.  et al., “Temporal evolution of strut light intensity after implantation of bioresorbable polymeric intracoronary scaffolds in the absorb cohort b trial,” Circ. J.. 78, (8 ), 1873 –1881 (2014).CrossRef
Nakatani  S.  et al., “Bioresorption and vessel wall integration of a fully bioresorbable polymeric everolimus-eluting scaffold: optical coherence tomography, intravascular ultrasound, and histological study in a porcine model with 4-year follow-up,” JACC: Cardiovasc. Interventions. 9, , 838 –851 (2016).CrossRef
Valery  T., Tissue Optics: Light Scattering Methods and Instruments for Medical Diagnosis. , 2nd ed., Vol. 23, ,  SPIE ,  Bellingham, Washington  (2007).
van Soest  G.  et al., “Atherosclerotic tissue characterization in vivo by optical coherence tomography attenuation imaging,” J. Biomed. Opt.. 15, (1 ), 011105  (2010). 1083-3668 CrossRef
Pedrotti  F. L., and Pedrotti  L. S., Introduction to Optics. , 2nd ed.,  Prentice-Hall International ,  Englewood Cliffs  (1993).
Sullivan  L. M., , Dukes  K. A., and Losina  E., “Tutorial in biostatistics. an introduction to hierarchical linear modelling,” Stat. Med.. 18, , 855 –888 (1999). 1097-0258 CrossRef
van der Meer  F. J., “Vascular applications of quantitative optical coherence tomography,” Dissertation, Universiteit van Amsterdam (2005).
Precht  H.  et al., “Optimisation of post mortem cardiac computed tomography compared to optical coherence tomography and histopathology,” J. Forensic Radiol. Imaging. 2, (2 ), 85 –90 (2014).CrossRef
Goutte  C., and Gaussier  E., “A probabilistic interpretation of precision, recall and f-score, with implication for evaluation,” in  European Conf. on Information Retrieval , pp. 345 –359,  Springer ,  Berlin Heidelberg  (2005).CrossRef

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

PubMed Articles
Advertisement
  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.