Lumen segmentation using a Mask R-CNN in carotid arteries with stenotic atherosclerotic plaque.

In patients at high risk for ischemic stroke, clinical carotid ultrasound is often used to grade stenosis, determine plaque burden and assess stroke risk. Analysis currently requires a trained sonographer to manually identify vessel and plaque regions, which is time and labor intensive. We present a method for automatically determining bounding boxes and lumen segmentation using a Mask R-CNN network trained on sonographer assisted ground-truth carotid lumen segmentations. Automatic lumen segmentation also lays the groundwork for developing methods for accurate plaque segmentation, and wall thickness measurements in cases with no plaque. Different training schemes are used to identify the Mask R-CNN model with the highest accuracy. Utilizing a single-channel B-mode training input, our model produces a mean bounding box intersection over union (IoU) of 0.81 and a mean lumen segmentation IoU of 0.75. However, we encountered errors in prediction when the jugular vein is the most prominently visualized vessel in the B-mode image. This was due to the fact that our dataset has limited instances of B-mode images with both the jugular vein and carotid artery where the vein is dominantly visualized. Additional training datasets are anticipated to mitigate this issue.

Changes in carotid artery texture by ultrasound and elastin features in a murine model.

Objective: In humans, arterial grayscale ultrasound texture features independently predict adverse cardiovascular disease (CVD) events and change with medical interventions. We performed this study to examine how grayscale ultrasound texture features and elastin fibers change in plaque-free segments of the arterial wall in a murine model prone to atherosclerosis. Methods: A total of 10 Apoetm1Unc/J mice (n = 5 male, n = 5 female) were imaged at 6, 16, and 24 weeks of age. Two mice were euthanized at 6 and 16 weeks and the remaining mice at 24 weeks. Texture features were extracted from the ultrasound images of the distal 1.0 mm of the common carotid artery wall, and elastin measures were extracted from histology images. Two-way analysis of variance was used to evaluate associations between week, sex, and grayscale texture features. Texture feature and elastin number comparisons between weeks were conducted using the sex-by-week two-way interaction contrasts. Sex-specific correlations between the number of elastin fibers and grayscale texture features were analyzed by conducting non-parametric Spearman's rank correlation analyses. Results: Arterial wall homogeneity changed significantly in male mice from 6 to 24 weeks, with a mean (SD) of 0.14 (0.03) units at 6 weeks and 0.18 (0.03) units at 24 weeks (p = 0.026). Spatial gray level dependence matrices-homogeneity (SGLD-HOM) also correlated with carotid artery plaque score (rs = 0.707, p = 0.033). Elastin fibers in the region of interest decreased from 6 to 24 weeks for both male and female mice, although only significantly in male mice. The mean (SD) number of elastin fibers for male mice was 5.32 (1.50) at 6 weeks and 3.59 (0.38) at 24 weeks (p = 0.023). For female mice, the mean (SD) number of elastin fibers was 3.98 (0.38) at 6 weeks and 3.46 (0.19) at 24 weeks (p = 0.051). Conclusion: Grayscale ultrasound texture features that are associated with increased risk for CVD events in humans were used in a murine model, and the grayscale texture feature SGLD-HOM was shown to change in male mice from 6 weeks to 24 weeks. Structural alterations of the arterial wall (change in elastin fiber number) were observed during this time and may differ by sex.

Ultrasonic array imaging of highly attenuative materials with spatio-temporal singular value decomposition.

With the increasing use of high density polyethylene (HDPE) pipes in nuclear industry which puts safety at the forefront, ultrasonic array imaging methods play a vital role in the structural integrity of HDPE pipe materials. However, the viscoelastic attenuation of HDPE pipe materials significantly decreases the level of signals, leading to a low signal-to-noise ratio caused by electronic noise. In this work, a domain-adapted spatio-temporal singular value decomposition (STSVD) processing algorithm combined with the total focusing method is proposed to improve the ultrasonic array image quality. First, the real-valued radio frequency (RF) data or A-scan signals are demodulated into the complex analytic signals containing in-phase and quadrature (I/Q) components. Then, the STSVD processing algorithm is used to filter the I/Q data, and the filtered I/Q data is converted into RF signals. Finally, the total focusing method is applied to the processed RF signals to produce the image of the region under detection as a stage of post-processing. Experiments are carried out with an ultrasonic linear phased array in contact with the HDPE pipe materials containing multiple side-drilled holes and through-wall notches. Results show that the proposed method can produce images with high quality to provide good inspection and characterization of defects in highly attenuative materials, especially the deeper defects.

Murine cardiac fibrosis localization using adaptive Bayesian cardiac strain imaging in vivo.

An adaptive Bayesian regularized cardiac strain imaging (ABR-CSI) algorithm for in vivo murine myocardial function assessment is presented. We report on 31 BALB/CJ mice (n = 17 females, n = 14 males), randomly stratified into three surgical groups: myocardial infarction (MI, n = 10), ischemia-reperfusion (IR, n = 13) and control (sham, n = 8) imaged pre-surgery (baseline- BL), and 1, 2, 7 and 14 days post-surgery using a high frequency ultrasound imaging system (Vevo 2100). End-systole (ES) radial and longitudinal strain images were used to generate cardiac fibrosis maps using binary thresholding. Percentage fibrotic myocardium (PFM) computed from regional fibrosis maps demonstrated statistically significant differences post-surgery in scar regions. For example, the MI group had significantly higher PFMRadial (%) values in the anterior mid region (p = 0.006) at Day 14 (n = 8, 42.30 ± 14.57) compared to BL (n = 12, 1.32 ± 0.85). A random forest classifier automatically detected fibrotic regions from ground truth Masson's trichrome stained histopathology whole slide images. Both PFMRadial (r = 0.70) and PFMLongitudinal (r = 0.60) results demonstrated strong, positive correlation with PFMHistopathology (p < 0.001).

Improving Minimum Variance Beamforming with Sub-Aperture Processing for Photoacoustic Imaging.

Minimum variance (MV) beamforming improves resolution and reduces sidelobes when compared to delay-and-sum (DAS) beamforming for photoacoustic imaging (PAI). However, some level of sidelobe signal and incoherent clutter persist degrading MV PAI quality. Here, an adaptive beamforming algorithm (PSAPMV) combining MV formulation and sub-aperture processing is proposed. In PSAPMV, the received channel data are split into two complementary nonoverlapping sub-apertures and beamformed using MV. A weighting matrix based on similarity between sub-aperture beamformed images was derived and multiplied with the full aperture MV image resulting in suppression of sidelobe and incoherent clutter in the PA image. Numerical simulation experiments with point targets, diffuse inclusions and microvasculature networks are used to validate PSAPMV. Quantitative evaluation was done in terms of main-lobe-to-side-lobe ratio, full width at half maximum (FWHM), contrast ratio (CR) and generalized contrast-to-noise ratio (gCNR). PSAPMV demonstrated improved beamforming performance both qualitatively and quantitatively. PSAPMV had higher resolution (FWHM =0.19 mm) than MV (0.21 mm) and DAS (0.22mm) in point target simulations, better target detectability (gCNR =0.99) than MV (0.89) and DAS (0.84) for diffuse inclusions and improved contrast (CR in microvasculature simulation, DAS = 15.38, MV = 22.42, PSAPMV = 51.74 dB).

Bayesian Regularized Strain Imaging for Assessment of Murine Cardiac Function In vivo.

A cardiac strain imaging framework with adaptive Bayesian regularization (ABR) is proposed for in vivo assessment of murine cardiac function. The framework uses ultrasound (US) radio-frequency data collected with a high frequency (fc = 30MHz) imaging system and a multi-level block matching algorithm with ABR to derive inter-frame cardiac displacements. Lagrangian cardiac strain (radial, er and longitudinal, el) tensors were derived by segmenting the myocardial wall starting at the ECG R-wave and accumulating interframe deformations over a cardiac cycle. In vivo feasibility was investigated through a longitudinal study with two mice (one ischemia-perfusion (IR) injury and one sham) imaged at five sessions (pre-surgery (BL) and 1,2,7 and 14 days post-surgery). End-systole (ES) strain images and segmental strain curves were derived for quantitative evaluation. Both mice showed periodic variation of er and el strain at BL with segmental synchroneity. Infarcted regions of IR mouse at Day 14 were associated with reduced or sign reversed ES er and el values while the sham mouse had similar or higher strain than at BL. Infarcted regions identified in vivo were associated with increased collagen content confirmed with Masson's Trichrome stained ex vivo heart sections.Clinical Relevance-Higher quality cardiac strain images derived with RF data and Bayesian regularization can potentially improve the sensitivity and accuracy of non-invasive assessment of cardiovascular disease models.

Enhancement of in vivo cardiac photoacoustic signal specificity using spatiotemporal singular value decomposition.

SIGNIFICANCE: Photoacoustic imaging (PAI) can be used to infer molecular information about myocardial health non-invasively in vivo using optical excitation at ultrasonic spatial resolution. For clinical and preclinical linear array imaging systems, conventional delay-and-sum (DAS) beamforming is typically used. However, DAS cardiac PA images are prone to artifacts such as diffuse quasi-static clutter with temporally varying noise-reducing myocardial signal specificity. Typically, multiple frame averaging schemes are utilized to improve the quality of cardiac PAI, which affects the spatial and temporal resolution and reduces sensitivity to subtle PA signal variation. Furthermore, frame averaging might corrupt myocardial oxygen saturation quantification due to the presence of natural cardiac wall motion. In this paper, a spatiotemporal singular value decomposition (SVD) processing algorithm is proposed to reduce DAS PAI artifacts and subsequent enhancement of myocardial signal specificity. AIM: Demonstrate enhancement of PA signals from myocardial tissue compared to surrounding tissues and blood inside the left-ventricular (LV) chamber using spatiotemporal SVD processing with electrocardiogram (ECG) and respiratory signal (ECG-R) gated in vivo murine cardiac PAI. APPROACH: In vivo murine cardiac PAI was performed by collecting single wavelength (850 nm) photoacoustic channel data on eight healthy mice. A three-dimensional (3D) volume of complex PAI data over a cardiac cycle was reconstructed using a custom ECG-R gating algorithm and DAS beamforming. Spatiotemporal SVD was applied on a two-dimensional Casorati matrix generated using the 3D volume of PAI data. The singular value spectrum (SVS) was then filtered to remove contributions from diffuse quasi-static clutter and random noise. Finally, SVD processed beamformed images were derived using filtered SVS and inverse SVD computations. RESULTS: Qualitative comparison with DAS and minimum variance (MV) beamforming shows that SVD processed images had better myocardial signal specificity, contrast, and target detectability. DAS, MV, and SVD images were quantitatively evaluated by calculating contrast ratio (CR), generalized contrast-to-noise ratio (gCNR), and signal-to-noise ratio (SNR). Quantitative evaluations were done at three cardiac time points (during systole, at end-systole (ES), and during diastole) identified from co-registered ultrasound M-Mode image. Mean CR, gCNR, and SNR values of SVD images at ES were 245, 115.15, and 258.17 times higher than DAS images with statistical significance evaluated with one-way analysis of variance. CONCLUSIONS: Our results suggest that significantly better-quality images can be realized using spatiotemporal SVD processing for in vivo murine cardiac PAI.

Locally optimized correlation-guided Bayesian adaptive regularization for ultrasound strain imaging.

Ultrasound strain imaging utilizes radio-frequency (RF) ultrasound echo signals to estimate the relative elasticity of tissue under deformation. Due to the diagnostic value inherent in tissue elasticity, ultrasound strain imaging has found widespread clinical and preclinical applications. Accurate displacement estimation using pre and post-deformation RF signals is a crucial first step to derive high quality strain tensor images. Incorporating regularization into the displacement estimation framework is a commonly employed strategy to improve estimation accuracy and precision. In this work, we propose an adaptive variation of the iterative Bayesian regularization scheme utilizing RF similarity metric signal-to-noise ratio previously proposed by our group. The regularization scheme is incorporated into a 2D multi-level block matching (BM) algorithm for motion estimation. Adaptive nature of our algorithm is attributed to the dynamic variation of iteration number based on the normalized cross-correlation (NCC) function quality and a similarity measure between pre-deformation and motion compensated post-deformation RF signals. The proposed method is validated for either quasi-static and cardiac elastography or strain imaging applications using uniform and inclusion phantoms and canine cardiac deformation simulation models. Performance of adaptive Bayesian regularization was compared to conventional NCC and Bayesian regularization with fixed number of iterations. Results from uniform phantom simulation study show significant improvement in lateral displacement and strain estimation accuracy. For instance, at 1.5% lateral strain in a uniform phantom, Bayesian regularization with five iterations incurred a lateral strain error of 104.49%, which was significantly reduced using our adaptive approach to 27.51% (p  < 0.001). Contrast-to-noise (CNR e ) ratios obtained from inclusion phantom indicate improved lesion detectability for both axial and lateral strain images. For instance, at 1.5% lateral strain, Bayesian regularization with five iterations had lateral CNR e of -0.31 dB which was significantly increased using the adaptive approach to 7.42 dB (p  < 0.001). Similar results are seen with cardiac deformation modelling with improvement in myocardial strain images. In vivo feasibility was also demonstrated using data from a healthy murine heart. Overall, the proposed method makes Bayesian regularization robust for clinical and preclinical applications.