A Mechanical Approach for Smooth Surface Fitting to Delineate Vessel Walls in Optical Coherence Tomography Images

Fully Automated Construction of Three-dimensional Finite Element Simulations from Optical Coherence Tomography

Despite recent advances in diagnosis and treatment, atherosclerotic coronary artery diseases remain a leading cause of death worldwide. Various imaging modalities and metrics can detect lesions and predict patients at risk; however, identifying unstable lesions is still difficult. Current techniques cannot fully capture the complex morphology-modulated mechanical responses that affect plaque stability, leading to catastrophic failure and mute the benefit of device and drug interventions. Finite Element (FE) simulations utilizing intravascular imaging OCT (Optical Coherence Tomography) are effective in defining physiological stress distributions. However, creating 3D FE simulations of coronary arteries from OCT images is challenging to fully automate given OCT frame sparsity, limited material contrast, and restricted penetration depth. To address such limitations, we developed an algorithmic approach to automatically produce 3D FE-ready digital twins from labeled OCT images. The 3D models are anatomically faithful and recapitulate mechanically relevant tissue lesion components, automatically producing morphologies structurally similar to manually constructed models whilst including more minute details. A mesh convergence study highlighted the ability to reach stress and strain convergence with average errors of just 5.9% and 1.6% respectively in comparison to FE models with approximately twice the number of elements in areas of refinement. Such an automated procedure will enable analysis of large clinical cohorts at a previously unattainable scale and opens the possibility for in-silico methods for patient specific diagnoses and treatment planning for coronary artery disease.

Segmentation of anatomical layers and imaging artifacts in intravascular polarization sensitive optical coherence tomography using attending physician and boundary cardinality losses

Intravascular ultrasound and optical coherence tomography are widely available for assessing coronary stenoses and provide critical information to optimize percutaneous coronary intervention. Intravascular polarization-sensitive optical coherence tomography (PS-OCT) measures the polarization state of the light scattered by the vessel wall in addition to conventional cross-sectional images of subsurface microstructure. This affords reconstruction of tissue polarization properties and reveals improved contrast between the layers of the vessel wall along with insight into collagen and smooth muscle content. Here, we propose a convolutional neural network model, optimized using two new loss terms (Boundary Cardinality and Attending Physician), that takes advantage of the additional polarization contrast and classifies the lumen, intima, and media layers in addition to guidewire and plaque shadows. Our model segments the media boundaries through fibrotic plaques and continues to estimate the outer media boundary behind shadows of lipid-rich plaques. We demonstrate that our multi-class classification model outperforms existing methods that exclusively use conventional OCT data, predominantly segment the lumen, and consider subsurface layers at most in regions of minimal disease. Segmentation of all anatomical layers throughout diseased vessels may facilitate stent sizing and will enable automated characterization of plaque polarization properties for investigation of the natural history and significance of coronary atheromas.

Validation of biomechanical assessment of coronary plaque vulnerability based on intravascular optical coherence tomography and digital subtraction angiography

Background

It has been suggested that biomechanical factors may influence plaque development. However, key determinants for assessing plaque vulnerability remain speculative.

Methods

In this study, a two-dimensional (2D) structural mechanical analysis and a three-dimensional (3D) fluid-structure interaction (FSI) analysis were conducted based on intravascular optical coherence tomography (IV-OCT) and digital subtraction angiography (DSA) data sets. In the 2D study, 103 IV-OCT slices were analyzed. An in-depth morpho-mechanic analysis and a weighted least absolute shrinkage and selection operator (LASSO) regression analysis were conducted to identify the crucial features related to plaque vulnerability via the tuning parameter (λ). In the 3D study, the coronary model was reconstructed by fusing the IV-OCT and DSA data, and a FSI analysis was subsequently performed. The relationship between vulnerable plaque and wall shear stress (WSS) was investigated.

Results

The influential factors were selected using the minimum criteria (λ-min) and one-standard error criteria (λ-1se). In addition to the common vulnerable factor of the minimum fibrous cap thickness (FCTmin), four biomechanical factors were selected by λ-min, including the average/maximal displacements and average/maximal stress, and two biomechanical factors were selected by λ-1se, including the average/maximal displacements. Additionally, the positions of the vulnerable plaques were consistent with the sites of high WSS.

Conclusions

Functional indices are crucial for plaque status assessment. An evaluation based on biomechanical simulations might provide insights into risk identification and guide therapeutic decisions.

Fully automated construction of three-dimensional finite element simulations from Optical Coherence Tomography

Despite recent advances in diagnosis and treatment, atherosclerotic coronary artery diseases remain a leading cause of death worldwide. Various imaging modalities and metrics can detect lesions and predict patients at risk; however, identifying unstable lesions is still difficult. Current techniques cannot fully capture the complex morphology-modulated mechanical responses that affect plaque stability, leading to catastrophic failure and mute the benefit of device and drug interventions. Finite Element (FE) simulations utilizing intravascular imaging OCT (Optical Coherence Tomography) are effective in defining physiological stress distributions. However, creating 3D FE simulations of coronary arteries from OCT images is challenging to fully automate given OCT frame sparsity, limited material contrast, and restricted penetration depth. To address such limitations, we developed an algorithmic approach to automatically produce 3D FE-ready digital twins from labeled OCT images. The 3D models are anatomically faithful and recapitulate mechanically relevant tissue lesion components, automatically producing morphologies structurally similar to manually constructed models whilst including more minute details. A mesh convergence study highlighted the ability to reach stress and strain convergence with average errors of just 5.9% and 1.6% respectively in comparison to FE models with approximately twice the number of elements in areas of refinement. Such an automated procedure will enable analysis of large clinical cohorts at a previously unattainable scale and opens the possibility for in-silico methods for patient specific diagnoses and treatment planning for coronary artery disease.

Comparison of multilayer and single-layer coronary plaque models on stress/strain calculations based on optical coherence tomography images

Mechanical stress and strain conditions are closely related to atherosclerotic plaque progression and rupture and have been under intensive investigations in recent years. It is well known that arteries have a three-layer structure: intima, media and adventitia. However, in vivo image-based multilayer plaque models are not available in the current literature due to lack of multilayer image segmentation data. A multilayer segmentation and repairing technique was introduced to segment coronary plaque optical coherence tomography (OCT) image to obtain its three-layer vessel structure. A total of 200 OCT slices from 20 patients (13 male; 7 female) were used to construct multilayer and single-layer 3D thin-slice models to calculate plaque stress and strain and compare model differences. Our results indicated that the average maximum plaque stress values of 20 patients from multilayer and single-layer models were 385.13 ± 110.09 kPa and 270.91 ± 95.86 kPa, respectively. The relative difference was 42.2%, with single-layer stress serving as the base value. The average mean plaque stress values from multilayer and single-layer models were 129.59 ± 32.77 kPa and 93.27 ± 18.20 kPa, respectively, with a relative difference of 38.9%. The maximum and mean plaque strain values obtained from the multilayer models were 11.6% and 19.0% higher than those from the single-layer models. Similarly, the maximum and mean cap strains showed increases of 9.6% and 12.9% over those from the single-layer models. These findings suggest that use of multilayer models could improve plaque stress and strain calculation accuracy and may have large impact on plaque progression and vulnerability investigation and potential clinical applications. Further large-scale studies are needed to validate our findings.

Human Coronary Plaque Optical Coherence Tomography Image Repairing, Multilayer Segmentation and Impact on Plaque Stress/Strain Calculations

Coronary vessel layer structure may have a considerable impact on plaque stress/strain calculations. Most current plaque models use single-layer vessel structures due to the lack of available multilayer segmentation techniques. In this paper, an automatic multilayer segmentation and repair method was developed to segment coronary optical coherence tomography (OCT) images to obtain multilayer vessel geometries for biomechanical model construction. Intravascular OCT data were acquired from six patients (one male; mean age: 70.0) using a protocol approved by the local institutional review board with informed consent obtained. A total of 436 OCT slices were selected in this study. Manually segmented data were used as the gold standard for method development and validation. The edge detection method and cubic spline surface fitting were applied to detect and repair the internal elastic membrane (IEM), external elastic membrane (EEM) and adventitia-periadventitia interface (ADV). The mean errors of automatic contours compared to manually segmented contours were 1.40%, 4.34% and 6.97%, respectively. The single-layer mean plaque stress value from lumen was 117.91 kPa, 10.79% lower than that from three-layer models (132.33 kPa). On the adventitia, the single-layer mean plaque stress value was 50.46 kPa, 156.28% higher than that from three-layer models (19.74 kPa). The proposed segmentation technique may have wide applications in vulnerable plaque research.

In Vivo Intravascular Optical Coherence Tomography (IVOCT) Structural and Blood Flow Imaging Based Mechanical Simulation Analysis of a Blood Vessel

INTRODUCTION

Computer modelling of blood vessels based on biomedical imaging provides important hemodynamic and biomechanical information for vascular disease studies and diagnosis. However due to lacking well-defined physiological blood flow profiles, the accuracy of the simulation results is often not guaranteed. Flow velocity profiles of a specific cross section of a blood vessel were obtained by in vivo Doppler intravascular optical coherence tomography (IVOCT) lately. However due to the influence of the catheter, the velocity profile imaged by IVOCT can't be applied to simulation directly.

METHODS

A simulation-experiment combined method to determine the inlet flow boundary based on in vivo porcine carotid Doppler IVOCT imaging is proposed. A single conduit carotid model was created from the 3D IVOCT structural images using an image processing-computer aided design combined method.

RESULTS

With both high- resolution arterial model and near physiological blood flow profile, stress analysis by fluid-structure interaction and computational fluid dynamics were performed. The influence of the catheter to the wall shear stress, the hemodynamics and the stresses of the carotid wall under the measured inlet flow and various outlet pressure boundary conditions, are analyzed.

CONCLUSION

This study provides a solution to the difficulty of getting inlet flow boundary for numerical simulation of arteries. It paves the way for developing IVOCT based vascular stress analysis and imaging methods for the studies and diagnosis of vascular diseases.

Framework for lumen-based nonrigid tomographic coregistration of intravascular images

Purpose: Modern medical imaging enables clinicians to effectively diagnose, monitor, and treat diseases. However, clinical decision-making often relies on combined evaluation of either longitudinal or disparate image sets, necessitating coregistration of multiple acquisitions. Promising coregistration techniques have been proposed; however, available methods predominantly rely on time-consuming manual alignments or nontrivial feature extraction with limited clinical applicability. Addressing these issues, we present a fully automated, robust, nonrigid registration method, allowing for coregistering of multimodal tomographic vascular image datasets using luminal annotation as the sole alignment feature. Approach: Registration is carried out by the use of the registration metrics defined exclusively for lumens shapes. The framework is primarily broken down into two sequential parts: longitudinal and rotational registration. Both techniques are inherently nonrigid in nature to compensate for motion and acquisition artifacts in tomographic images. Results: Performance was evaluated across multimodal intravascular datasets, as well as in longitudinal cases assessing pre-/postinterventional coronary images. Low registration error in both datasets highlights method utility, with longitudinal registration errors-evaluated throughout the paired tomographic sequences-of 0.29 ± 0.14    mm ( < 2 longitudinal image frames) and 0.18 ± 0.16    mm ( < 1 frame) for multimodal and interventional datasets, respectively. Angular registration for the interventional dataset rendered errors of 7.7 ° ± 6.7 ° , and 29.1 ° ± 23.2 ° for the multimodal set. Conclusions: Satisfactory results across datasets, along with additional attributes such as the ability to avoid longitudinal over-fitting and correct nonlinear catheter rotation during nonrigid rotational registration, highlight the potential wide-ranging applicability of our presented coregistration method.

Automated Coronary Optical Coherence Tomography Feature Extraction with Application to Three-Dimensional Reconstruction

Coronary optical coherence tomography (OCT) is an intravascular, near-infrared light-based imaging modality capable of reaching axial resolutions of 10-20 µm. This resolution allows for accurate determination of high-risk plaque features, such as thin cap fibroatheroma; however, visualization of morphological features alone still provides unreliable positive predictive capability for plaque progression or future major adverse cardiovascular events (MACE). Biomechanical simulation could assist in this prediction, but this requires extracting morphological features from intravascular imaging to construct accurate three-dimensional (3D) simulations of patients' arteries. Extracting these features is a laborious process, often carried out manually by trained experts. To address this challenge, numerous techniques have emerged to automate these processes while simultaneously overcoming difficulties associated with OCT imaging, such as its limited penetration depth. This systematic review summarizes advances in automated segmentation techniques from the past five years (2016-2021) with a focus on their application to the 3D reconstruction of vessels and their subsequent simulation. We discuss four categories based on the feature being processed, namely: coronary lumen; artery layers; plaque characteristics and subtypes; and stents. Areas for future innovation are also discussed as well as their potential for future translation.

An inverse method for mechanical characterization of heterogeneous diseased arteries using intravascular imaging

The increasing prevalence of finite element (FE) simulations in the study of atherosclerosis has spawned numerous inverse FE methods for the mechanical characterization of diseased tissue in vivo. Current approaches are however limited to either homogenized or simplified material representations. This paper presents a novel method to account for tissue heterogeneity and material nonlinearity in the recovery of constitutive behavior using imaging data acquired at differing intravascular pressures by incorporating interfaces between various intra-plaque tissue types into the objective function definition. Method verification was performed in silico by recovering assigned material parameters from a pair of vessel geometries: one derived from coronary optical coherence tomography (OCT); one generated from in silico-based simulation. In repeated tests, the method consistently recovered 4 linear elastic (0.1 ± 0.1% error) and 8 nonlinear hyperelastic (3.3 ± 3.0% error) material parameters. Method robustness was also highlighted in noise sensitivity analysis, where linear elastic parameters were recovered with average errors of 1.3 ± 1.6% and 8.3 ± 10.5%, at 5% and 20% noise, respectively. Reproducibility was substantiated through the recovery of 9 material parameters in two more models, with mean errors of 3.0 ± 4.7%. The results highlight the potential of this new approach, enabling high-fidelity material parameter recovery for use in complex cardiovascular computational studies.

A platform for high-fidelity patient-specific structural modelling of atherosclerotic arteries: from intravascular imaging to three-dimensional stress distributions

The pathophysiology of atherosclerotic lesions, including plaque rupture triggered by mechanical failure of the vessel wall, depends directly on the plaque morphology-modulated mechanical response. The complex interplay between lesion morphology and structural behaviour can be studied with high-fidelity computational modelling. However, construction of three-dimensional (3D) and heterogeneous models is challenging, with most previous work focusing on two-dimensional geometries or on single-material lesion compositions. Addressing these limitations, we here present a semi-automatic computational platform, leveraging clinical optical coherence tomography images to effectively reconstruct a 3D patient-specific multi-material model of atherosclerotic plaques, for which the mechanical response is obtained by structural finite-element simulations. To demonstrate the importance of including multi-material plaque components when recovering the mechanical response, a computational case study was conducted in which systematic variation of the intraplaque lipid and calcium was performed. The study demonstrated that the inclusion of various tissue components greatly affected the lesion mechanical response, illustrating the importance of multi-material formulations. This platform accordingly provides a viable foundation for studying how plaque micro-morphology affects plaque mechanical response, allowing for patient-specific assessments and extension into clinically relevant patient cohorts.

Clinical considerations and recommendations for OCT-guided carotid artery stenting

Introduction: Optical Coherence Tomography (OCT) is an intravascular imaging providing high-resolution images of morphological features of arterial wall. Nowadays, OCT is an accepted intravascular modality to study coronary arteries, stent implantation, and vessel injury. In the last decade, an increasing interest have been focused on the application of OCT in carotid arteries.Areas covered: Literature evidence in the application of OCT in carotid arteries still remains debated. So far, OCT has been used as a research tool, aiming to evaluate atherosclerotic plaques' features and stents' behavior after implantation. This paper is intended to summarize clinical evidences and practices in the use of OCT in carotid arteries district and during CAS procedures. Literature review was completed via Pubmed search using Keywords.Expert opinion: CAS is a safe and effective procedure when performed by trained physicians with a tailored approach. In this scenario, ambiguous pictures at ultrasound, angiography, and IVUS might be clarified using OCT.By providing unprecedented microstructural information on atherosclerotic plaques, OCT may identify the features of vulnerable carotid plaque and, by identifying possible defects after stent implantation as malapposition and plaque prolapse, it may help the tailoring approach to CAS.

Automated accurate lumen segmentation using L-mode interpolation for three-dimensional intravascular optical coherence tomography

Intravascular optical coherence tomography (IVOCT) lumen-based computational flow dynamics (CFD) enables physiologic evaluations such as of the fractional flow reserve (FFR) and wall sheer stress. In this study, we developed an accurate, time-efficient method for extracting lumen contours of the coronary artery. The contours of cross-sectional images containing wide intimal discontinuities due to guide wire shadowing and large bifurcations were delineated by utilizing the natural longitudinal lumen continuity of the arteries. Our algorithm was applied to 5931 pre-intervention OCT images acquired from 40 patients. For a quantitative comparison, the images were also processed through manual segmentation (the reference standard) and automated ones utilizing cross-sectional and longitudinal continuities. The results showed that the proposed algorithm outperforms other schemes, exhibiting a strong correlation (R = 0.988) and overlapping and non-overlapping area ratios of 0.931 and 0.101, respectively. To examine the accuracy of the OCT-derived FFR calculated using the proposed scheme, a CFD simulation of a three-dimensional coronary geometry was performed. The strong correlation with a manual lumen-derived FFR (R = 0.978) further demonstrated the reliability and accuracy of our algorithm with potential applications in clinical settings.