Background Removal and Vessel Filtering of Noncontrast Ultrasound Images of Microvasculature

Quantitative Biomarkers Derived from a Novel, Contrast-Free Ultrasound, High-Definition Microvessel Imaging for Differentiating Choroidal Tumors

Angiogenesis has an essential role in the de novo evolution of choroidal melanoma as well as choroidal nevus transformation into melanoma. Differentiating early-stage melanoma from nevus is of high clinical importance; thus, imaging techniques that provide objective information regarding tumor microvasculature structures could aid accurate early detection. Herein, we investigated the feasibility of quantitative high-definition microvessel imaging (qHDMI) for differentiation of choroidal tumors in humans. This new ultrasound-based technique encompasses a series of morphological filtering and vessel enhancement techniques, enabling the visualization of tumor microvessels as small as 150 microns and extracting vessel morphological features as new tumor biomarkers. Distributional differences between the malignant melanomas and benign nevi were tested on 37 patients with choroidal tumors using a non-parametric Wilcoxon rank-sum test, and statistical significance was declared for biomarkers with p-values < 0.05. The ocular oncology diagnosis was choroidal melanoma (malignant) in 21 and choroidal nevus (benign) in 15 patients. The mean thickness of benign and malignant masses was 1.70 ± 0.40 mm and 3.81 ± 2.63 mm, respectively. Six HDMI biomarkers, including number of vessel segments (p = 0.003), number of branch points (p = 0.003), vessel density (p = 0.03), maximum tortuosity (p = 0.001), microvessel fractal dimension (p = 0.002), and maximum diameter (p = 0.003) exhibited significant distributional differences between the two groups. Contrast-free HDMI provided noninvasive imaging and quantification of microvessels of choroidal tumors. The results of this pilot study indicate the potential use of qHDMI as a complementary tool for characterization of small ocular tumors and early detection of choroidal melanoma.

Reduced osteoprotegerin expression by osteocytes may contribute to rebound resorption after denosumab discontinuation

Denosumab is an anti-RANKL Ab that potently suppresses bone resorption, increases bone mass, and reduces fracture risk. Discontinuation of denosumab causes rapid rebound bone resorption and bone loss, but the molecular mechanisms are unclear. We generated humanized RANKL mice and treated them with denosumab to examine the cellular and molecular conditions associated with rebound resorption. Denosumab potently suppressed both osteoclast and osteoblast numbers in cancellous bone in humanized RANKL mice. The decrease in osteoclast number was not associated with changes in osteoclast progenitors in bone marrow. Long-term, but not short-term, denosumab administration reduced osteoprotegerin (OPG) mRNA in bone. Localization of OPG expression revealed that OPG mRNA is produced by a subpopulation of osteocytes. Long-term denosumab administration reduced osteocyte OPG mRNA, suggesting that OPG expression declines as osteocytes age. Consistent with this, osteocyte expression of OPG was more prevalent near the surface of cortical bone in humans and mice. These results suggest that new osteocytes are an important source of OPG in remodeling bone and that suppression of remodeling reduces OPG abundance by reducing new osteocyte formation. The lack of new osteocytes and the OPG they produce may contribute to rebound resorption after denosumab discontinuation.

Multi-layered adaptive neoangiogenesis Intra-Operative quantification (MANIOQ)

Quantification of vascularization volume can provide valuable information for diagnosis and prognosis in vascular pathologies. It can be adapted to inform the surgical management of gliomas, aggressive brain tumors characterized by exuberant sprouting of new blood vessels (neoangiogenesis). Filtered ultrafast Doppler data can provide two main parameters: vascularization index (VI) and fractional moving blood volume (FMBV) that clinically reflect tumor micro vascularization. Current protocols lack robust, automatic, and repeatable filtering methods. We present a filtrating method called Multi-layered Adaptive Neoangiogenesis Intra-Operative Quantification (MANIOQ). First, an adaptive clutter filtering is implemented, based on singular value decomposition (SVD) and hierarchical clustering. Second a method for noise equalization is applied, based on the subtraction of a weighted noise profile. Finally, an in vivo analysis of the periphery of the B-mode hyper signal area allows to measure the vascular infiltration extent of the brain tumors. Ninety ultrasound acquisitions were processed from 23 patients. Compared to reference methods in the literature, MANIOQ provides a more robust tissue filtering, and noise equalization allows for the first time to keep axial and lateral gain compensation (TGC and LGC). MANIOQ opens the way to an intra-operative clinical analysis of gliomas micro vascularization.

Custom-made flow phantoms for quantitative ultrasound microvessel imaging

Morphologically realistic flow phantoms are essential experimental tools for quantitative ultrasound-based microvessel imaging. As new quantitative flow imaging tools are developed, the need for more complex vessel-mimicking phantoms is indisputable. In this article, we propose a method for fabricating phantoms with sub-millimeter channels consisting of branches and curvatures in various shapes and sizes suitable for quantifying vessel morphological features. We used different tissue-mimicking materials (TMMs) compatible with ultrasound imaging as the base and metal wires of different diameters (0.15-1.25 mm) to create wall-less channels. The TMMs used are silicone rubber, plastisol, conventional gelatin, and medical gelatin. Mother channels in these phantoms were made in diameters of 1.25 mm or 0.3 mm and the daughter channels in diameters 0.3 mm or 0.15 mm. Bifurcations were created by soldering wires together at branch points. Quantitative parameters were assessed, and accuracy of measurements from the ground truth were determined. Channel diameters were seen to have increased (76-270%) compared to the initial state in the power Doppler images, partly due to blood mimicking fluid pressure. Amongst the microflow phantoms made from the different TMMs, the medical gelatin phantom was selected as the best option for microflow imaging, fulfilling the objective of being easy to fabricate with high transmittance while having a speed of sound and acoustic attenuation close to human tissue. A flow velocity of 0.85 ± 0.01 mm/s, comparable to physiological flow velocity was observed in the smallest diameter phantom (medical gelatin branch) presented here. We successfully constructed more complex geometries, including tortuous and multibranch channels using the medical gelatin as the TMM. We anticipate this will create new avenues for validating quantitative ultrasound microvessel imaging techniques.

Improved Ultrafast Power Doppler Imaging Using United Spatial-Angular Adaptive Scaling Wiener Postfilter

Ultrafast power Doppler imaging (uPDI) using high-frame-rate plane-wave transmission is a new microvascular imaging modality that offers high Doppler sensitivity. However, due to the unfocused transmission of plane waves, the echo signal is subject to interference from noise and clutter, resulting in a low signal-to-noise ratio (SNR) and poor image quality. Adaptive beamforming techniques are effective in suppressing noise and clutter for improved image quality. In this study, an adaptive beamformer based on a united spatial-angular adaptive scaling Wiener (uSA-ASW) postfilter is proposed to improve the resolution and contrast of uPDI. In the proposed method, the signal power and noise power of the Wiener postfilter are estimated by uniting spatial and angular signals, and a united generalized coherence factor (uGCF) is introduced to dynamically adjust the noise power estimation and enhance the robustness of the method. Simulation and in vivo data were used to verify the effectiveness of the proposed method. The results show that the uSA-ASW can achieve higher resolution and significant improvements in image contrast and background noise suppression compared with conventional delay-and-sum (DAS), coherence factor (CF), spatial-angular CF (SACF), and adaptive scaling Wiener (ASW) postfilter methods. In the simulations, uSA-ASW improves contrast-to-noise ratio (CNR) by 34.7 dB (117.3%) compared with DAS, while reducing background noise power (BNP) by 52 dB (221.4%). The uSA-ASW method provides full-width at half-maximum (FWHM) reductions of [Formula: see text] (59.5%) and [Formula: see text] (56.9%), CNR improvements of 25.6 dB (199.9%) and 42 dB (253%), and BNP reductions of 46.1 dB (319.3%) and 12.9 dB (289.1%) over DAS in the experiments of contrast-free human neonatal brain and contrast-free human liver, respectively. In the contrast-free experiments, uSA-ASW effectively balances the performance of noise and clutter suppression and enhanced microvascular visualization. Overall, the proposed method has the potential to become a reliable microvascular imaging technique for aiding in more accurate diagnosis and detection of vascular-related diseases in clinical contexts.

Mapping Pharmacologically Evoked Neurovascular Activation and Its Suppression in a Rat Model of Tremor Using Functional Ultrasound: A Feasibility Study

Functional ultrasound (fUS), an emerging hemodynamic-based functional neuroimaging technique, is especially suited to probe brain activity and primarily used in animal models. Increasing use of pharmacological models for essential tremor extends new research to the utilization of fUS imaging in such models. Harmaline-induced tremor is an easily provoked model for the development of new therapies for essential tremor (ET). Furthermore, harmaline-induced tremor can be suppressed by the same classic medications used for essential tremor, which leads to the utilization of this model for preclinical testing. However, changes in local cerebral activities under the effect of tremorgenic doses of harmaline have not been completely investigated. In this study, we explored the feasibility of fUS imaging for visualization of cerebral activation and deactivation associated with harmaline-induced tremor and tremor-suppressing effects of propranolol. The spatial resolution of fUS using a high frame rate imaging enabled us to visualize time-locked and site-specific changes in cerebral blood flow associated with harmaline-evoked tremor. Intraperitoneal administration of harmaline generated significant neural activity changes in the primary motor cortex and ventrolateral thalamus (VL Thal) regions during tremor and then gradually returned to baseline level as tremor subsided with time. To the best of our knowledge, this is the first functional ultrasound study to show the neurovascular activation of harmaline-induced tremor and the therapeutic suppression in a rat model. Thus, fUS can be considered a noninvasive imaging method for studying neuronal activities involved in the ET model and its treatment.

Contrast-free ultrasound imaging for blood flow assessment of the lower limb in patients with peripheral arterial disease: a feasibility study

While being a relatively prevalent condition particularly among aging patients, peripheral arterial disease (PAD) of lower extremities commonly goes undetected or misdiagnosed due to its symptoms being nonspecific. Additionally, progression of PAD in the absence of timely intervention can lead to dire consequences. Therefore, development of non-invasive and affordable diagnostic approaches can be highly beneficial in detection and treatment planning for PAD patients. In this study, we present a contrast-free ultrasound-based quantitative blood flow imaging technique for PAD diagnosis. The method involves monitoring the variations of blood flow in the calf muscle in response to thigh-pressure-cuff-induced occlusion. Four quantitative metrics are introduced for analysis of these variations. These metrics include post-occlusion to baseline flow intensity variation (PBFIV), total response region (TRR), Lag0 response region (L0RR), and Lag4 (and more) response region (L4 + RR). We examine the feasibility of this method through an in vivo study consisting of 14 PAD patients with abnormal ankle-brachial index (ABI) and 8 healthy volunteers. Ultrasound data acquired from 13 legs in the patient group and 13 legs in the healthy group are analyzed. Out of the four utilized metrics, three exhibited significantly different distributions between the two groups (p-value < 0.05). More specifically, p-values of 0.0015 for PBFIV, 0.0183 for TRR, and 0.0048 for L0RR were obtained. The results of this feasibility study indicate the diagnostic potential of the proposed method for the detection of PAD.

Noninvasive prediction of axillary lymph node breast cancer metastasis using morphometric analysis of nodal tumor microvessels in a contrast-free ultrasound approach

PURPOSE

Changes in microcirculation of axillary lymph nodes (ALNs) may indicate metastasis. Reliable noninvasive imaging technique to quantify such variations is lacking. We aim to develop and investigate a contrast-free ultrasound quantitative microvasculature imaging technique for detection of metastatic ALN in vivo.

EXPERIMENTAL DESIGN

The proposed ultrasound-based technique, high-definition microvasculature imaging (HDMI) provides superb images of tumor microvasculature at sub-millimeter size scales and enables quantitative analysis of microvessels structures. We evaluated the new HDMI technique on 68 breast cancer patients with ultrasound-identified suspicious ipsilateral axillary lymph nodes recommended for fine needle aspiration biopsy (FNAB). HDMI was conducted before the FNAB and vessel morphological features were extracted, analyzed, and the results were correlated with the histopathology.

RESULTS

Out of 15 evaluated quantitative HDMI biomarkers, 11 were significantly different in metastatic and reactive ALNs (10 with P << 0.01 and one with 0.01 < P < 0.05). We further showed that through analysis of these biomarkers, a predictive model trained on HDMI biomarkers combined with clinical information (i.e., age, node size, cortical thickness, and BI-RADS score) could identify metastatic lymph nodes with an area under the curve of 0.9 (95% CI [0.82,0.98]), sensitivity of 90%, and specificity of 88%.

CONCLUSIONS

The promising results of our morphometric analysis of HDMI on ALNs offer a new means of detecting lymph node metastasis when used as a complementary imaging tool to conventional ultrasound. The fact that it does not require injection of contrast agents simplifies its use in routine clinical practice.

High-Quality Ultrafast Power Doppler Imaging Based on Spatial Angular Coherence Factor

The morphological and hemodynamic changes of microvessels are demonstrated to be related to the diseased conditions in tissues. Ultrafast power Doppler imaging (uPDI) is a novel modality with a significantly increased Doppler sensitivity, benefiting from the ultrahigh frame rate plane-wave imaging (PWI) and advanced clutter filtering. However, unfocused plane-wave transmission often leads to a low imaging quality, which degrades the subsequent microvascular visualization in power Doppler imaging. Coherence factor (CF)-based adaptive beamformers have been widely studied in conventional B-mode imaging. In this study, we propose a spatial and angular coherence factor (SACF) beamformer for improved uPDI (SACF-uPDI) by calculating the spatial CF across apertures and the angular CF across transmit angles, respectively. To identify the superiority of SACF-uPDI, simulations, in vivo contrast-enhanced rat kidney, and in vivo contrast-free human neonatal brain studies were conducted. Results demonstrate that SACF-uPDI can effectively enhance contrast and resolution and suppress background noise simultaneously, compared with conventional uPDI methods based on delay-and-sum (DAS) (DAS-uPDI) and CF (CF-uPDI). In the simulations, SACF-uPDI can improve the lateral and axial resolutions compared with those of DAS-uPDI, from 176 to [Formula: see text] of lateral resolution, and from 111 to [Formula: see text] of axial resolution. In the in vivo contrast-enhanced experiments, SACF achieves 15.14- and 5.6-dB higher contrast-to-noise ratio (CNR), 15.25- and 3.68-dB lower noise power, and 240- and 15- [Formula: see text] narrower full-width at half-maximum (FWHM) than DAS-uPDI and CF-uPDI, respectively. In the in vivo contrast-free experiments, SACF achieves 6.11- and 1.09-dB higher CNR, 11.93- and 4.01-dB lower noise power, and 528- and 160- [Formula: see text] narrower FWHM than DAS-uPDI and CF-uPDI, respectively. In conclusion, the proposed SACF-uPDI method can efficiently improve the microvascular imaging quality and has the potential to facilitate clinical applications.

Quantitative Biomarkers Derived from a Novel Contrast-Free Ultrasound High-Definition Microvessel Imaging for Distinguishing Thyroid Nodules

Low specificity in current ultrasound modalities for thyroid cancer detection necessitates the development of new imaging modalities for optimal characterization of thyroid nodules. Herein, the quantitative biomarkers of a new high-definition microvessel imaging (HDMI) were evaluated for discrimination of benign from malignant thyroid nodules. Without the help of contrast agents, this new ultrasound-based quantitative technique utilizes processing methods including clutter filtering, denoising, vessel enhancement filtering, morphological filtering, and vessel segmentation to resolve tumor microvessels at size scales of a few hundred microns and enables the extraction of vessel morphological features as new tumor biomarkers. We evaluated quantitative HDMI on 92 patients with 92 thyroid nodules identified in ultrasound. A total of 12 biomarkers derived from vessel morphological parameters were associated with pathology results. Using the Wilcoxon rank-sum test, six of the twelve biomarkers were significantly different in distribution between the malignant and benign nodules (all p < 0.01). A support vector machine (SVM)-based classification model was trained on these six biomarkers, and the receiver operating characteristic curve (ROC) showed an area under the curve (AUC) of 0.9005 (95% CI: [0.8279,0.9732]) with sensitivity, specificity, and accuracy of 0.7778, 0.9474, and 0.8929, respectively. When additional clinical data, namely TI-RADS, age, and nodule size were added to the features, model performance reached an AUC of 0.9044 (95% CI: [0.8331,0.9757]) with sensitivity, specificity, and accuracy of 0.8750, 0.8235, and 0.8400, respectively. Our findings suggest that tumor vessel morphological features may improve the characterization of thyroid nodules.

Region-based SVD processing of high-frequency ultrafast ultrasound to visualize cutaneous vascular networks

Observing alterations in cutaneous vasculature in response to any disease or pathology is considered as a potential diagnostic marker in the progression and cure of a disease. To observe skin morphologies and tissue conditions, high-frequency ultrasound (HFUS) has been used in dermatology, although its ability to selectively visualize micro-vessels is limited due to insufficient Doppler sensitivity to peripheral slow-speed blood flow. In recent studies, this issue has been improved by increasing the sensitivity of Doppler imaging to slow flow, leveraging advanced cutter filtering approaches based on singular value decomposition (SVD) techniques that aid to effectively extract flow signals overlapped with tissue echo signals. Nevertheless, in skin imaging, variations in flow speed, diameter, and depth of the blood vessels at different skin layers can make clutter filtering challenging because these variations are problematic in selecting the optimal cut-off value for the SVD filtering. In this study, we aimed to devise a novel region-based SVD filtering approach for ultrafast HFUS data to visualize cutaneous vascular networks. The proposed method divides the acquired high-framerate data into two regions based on B-mode cutaneous morphological identification (dermis layer and subcutaneous tissue). Singular value decomposition processing was performed on each region to effectively extract the desired flow signal, and the processed regions were merged to generate a single power Doppler image, thereby highlighting the appearance of a complete cutaneous vascular network. Finally, top-hat transform was applied to the power Doppler image to further suppress the background noises and enhances the visibility of the micro-vessels. Experimental observations of the human cutaneous circulation showed that the image quality (contrast-to-noise ratio) through the region-based SVD filtering was measured to be 4.1 dB (before top-hat filtering) and 5.2 dB (after top-hat filtering), which were improved from 3.4 dB and 4.0 dB obtained using the global SVD approach with and without top-hat filtering, respectively. We envisioned that this approach would provide diverse applications in the diagnosis of cutaneous disorders.

Resistivity index mapping in Kidney based on ultrasensitive Pulsed-Wave Doppler and automatic spectrogram envelope detection

In recent years, ultrasensitive Pulsed-Wave Doppler (uPWD) ultrasound (US) has emerged as an alternative imaging approach for microcirculation imaging and as a complementary tool to other imaging modalities, such as positron emission tomography (PET). uPWD is based on the acquisition of a large set of highly spatiotemporally coherent frames, which allows high-quality images of a wide field of view to be obtained. In addition, these acquired frames allow calculation of the resistivity index (RI) of the pulsatile flow detected over the entire field of view, which is of great interest to clinicians, for example, in monitoring the transplanted kidney course. This work aims to develop and evaluate a method to automatically obtain an RI map of the kidney based on the uPWD approach. The effect of time gain compensation (TGC) on the visualization of vascularization and aliasing on the blood flow frequency response, was also assessed. A pilot study conducted in patients referred for renal transplant Doppler examination showed that the proposed method provided relative errors of about 15% for RI measurements with respect to conventional pulsed-wave (PW) Doppler.

Morphometric analysis of tumor microvessels for detection of hepatocellular carcinoma using contrast-free ultrasound imaging: A feasibility study

Introduction

A contrast-free ultrasound microvasculature imaging technique was evaluated in this study to determine whether extracting morphological features of the vascular networks in hepatic lesions can be beneficial in differentiating benign and malignant tumors (hepatocellular carcinoma (HCC) in particular).

Methods

A total of 29 lesions from 22 patients were included in this work. A post-processing algorithm consisting of clutter filtering, denoising, and vessel enhancement steps was implemented on ultrasound data to visualize microvessel structures. These structures were then further characterized and quantified through additional image processing. A total of nine morphological metrics were examined to compare different groups of lesions. A two-sided Wilcoxon rank sum test was used for statistical analysis.

Results

In the malignant versus benign comparison, six of the metrics manifested statistical significance. Comparing only HCC cases with the benign, only three of the metrics were significantly different. No statistically significant distinction was observed between different malignancies (HCC versus cholangiocarcinoma and metastatic adenocarcinoma) for any of the metrics.

Discussion

Obtained results suggest that designing predictive models based on such morphological characteristics on a larger sample size may prove helpful in differentiating benign from malignant liver masses.

Combining Spatial Registration With Clutter Filtering for Power-Doppler Imaging in Peripheral Perfusion Applications

Power-Doppler ultrasonic (PD-US) imaging is sensitive to echoes from blood cell motion in the microvasculature but generally nonspecific because of difficulties with filtering nonblood-echo sources. We are studying the potential for using PD-US imaging for routine assessments of peripheral blood perfusion without contrast media. The strategy developed is based on an experimentally verified computational model of tissue perfusion that simulates typical in vivo conditions. The model considers directed and diffuse blood perfusion states in a field of moving clutter and noise. A spatial registration method is applied to minimize tissue motion prior to clutter and noise filtering. The results show that in-plane clutter motion is effectively minimized. While out-of-plane motion remains a strong source of clutter-filter leakage, those registration errors are readily minimized by straightforward modification of scanning techniques and spatial averaging.

Flow imaging using differential Golay encoded ultrasound

In our research we present a new method of differential compression of the Golay encoded ultrasound (DCGEU) in the standard beamforming mode to visualize the slow (<1cm/s) blood mimicking fluid flow in small diameter tubes. The proposed DCGEU method is based on synthesis of several subsequent B-mode frames acquired with certain time intervals (30 ms in this study) followed by the visualization of differential beamformed radio frequency (RF) echoes, which yielded the images of the scatterers moving slowly in the vessel and suppressing the static echoes outside the vessel. In order to extract small backscattered echoes from the vessel area we took an advantage of improved sensitivity of the complementary Golay coded sequences (CGCS). The validation of the proposed DCGEU method was carried out in two stages. In the first one, we compared the flow images in small tubes with a diameter of 1 mm and 2.5 mm, reconstructed from numerically simulated acoustic data for the standard transmission of short pulses and 16-bits long CGCS signals. In the second stage of the research, the experimental data were acquired in a flow phantom with silicone tubes with an internal diameter of 1.5 mm and 4.5 mm and a fluid flow velocity of 0.9 cm/s. The experiments were carried out using preprogrammed Verasonics Vantage™ research ultrasound system equipped with ALT L12-5/50 mm MHz linear array transducer with 7.8 MHz center frequency. It was evidenced both in simulations and experiments that the DCGEU provided a good flow image along the entire length of tubing with virtually angle independent detection in comparison with the conventional short pulse interrogation.

Ultrasound high-definition microvasculature imaging with novel quantitative biomarkers improves breast cancer detection accuracy

OBJECTIVES

To overcome the limitations of power Doppler in imaging angiogenesis, we sought to develop and investigate new quantitative biomarkers of a contrast-free ultrasound microvasculature imaging technique for differentiation of benign from malignant pathologies of breast lesion.

METHODS

In this prospective study, a new high-definition microvasculature imaging (HDMI) was tested on 521 patients with 527 ultrasound-identified suspicious breast masses indicated for biopsy. Four new morphological features of tumor microvessels, microvessel fractal dimension (mvFD), Murray's deviation (MD), bifurcation angle (BA), and spatial vascularity pattern (SVP) as well as initial biomarkers were extracted and analyzed, and the results correlated with pathology. Multivariable logistic regression analysis was used to study the performance of different prediction models, initial biomarkers, new biomarkers, and combined new and initial biomarkers in differentiating benign from malignant lesions.

RESULTS

The new HDMI biomarkers, mvFD, BA, MD, and SVP, were statistically significantly different in malignant and benign lesions, regardless of tumor size. Sensitivity and specificity of the new biomarkers in lesions > 20 mm were 95.6% and 100%, respectively. Combining the new and initial biomarkers together showed an AUC, sensitivity, and specificity of 97% (95% CI: 95-98%), 93.8%, and 89.2%, respectively, for all lesions regardless of mass size. The classification was further improved by adding the Breast Imaging Reporting and Data System (BI-RADS) score to the prediction model, showing an AUC, sensitivity, and specificity of 97% (95% CI: 95-98%), 93.8%, and 89.2%, respectively.

CONCLUSION

The addition of new quantitative HDMI biomarkers significantly improved the accuracy in breast lesion characterization when used as a complementary imaging tool to the conventional ultrasound.

KEY POINTS

• Novel quantitative biomarkers extracted from tumor microvessel images increase the sensitivity and specificity in discriminating malignant from benign breast masses. • New HDMI biomarkers Murray's deviation, bifurcation angles, microvessel fractal dimension, and spatial vascularity pattern outperformed the initial biomarkers. • The addition of BI-RADS scores based on US descriptors to the multivariable analysis using all biomarkers remarkably increased the sensitivity, specificity, and AUC in all size groups.

Improved Ultrafast Power Doppler Imaging by Using Spatiotemporal Non-Local Means Filtering

The change of microvasculature is associated with the occurrence and development of many diseases. Ultrafast power Doppler imaging (uPDI) is an emerging technology for the visualization of microvessels due to the development of ultrafast plane wave (PW) imaging and advanced clutter filters. However, the low signal-to-noise ratio (SNR) caused by unfocused transmit of PW imaging deteriorates the subsequent imaging of microvasculature. Nonlocal means (NLM) filtering has been demonstrated to be effective in the denoising of both natural and medical images, including ultrasound power Doppler images. However, the feasibility and performance of applying an NLM filter on the ultrasound radio frequency (RF) data have not been investigated so far. In this study, we propose to apply an NLM filter on the spatiotemporal domain of clutter filtered blood flow RF data (St-NLM) to improve the quality of uPDI. Experiments were conducted to compare the proposed method with three different methods (under various similarity window sizes), including conventional uPDI without NLM filtering (Non-NLM), NLM filtering on the obtained power Doppler images (PD-NLM), and NLM filtering on the spatial domain of clutter filtered blood flow RF data (S-NLM). Phantom experiments, in vivo contrast-enhanced human spinal cord tumor experiments, and in vivo contrast-free human liver experiments were performed to demonstrate the superiority of the proposed St-NLM method over the other three methods. Qualitative and quantitative results show that the proposed St-NLM method can effectively suppress the background noise, improve the contrast between vessels and background, and preserve the details of small vessels at the same time. In the human liver study, the proposed St-NLM method achieves 31.05-, 24.49-, and 11.15-dB higher contrast-to-noise ratios (CNRs) and 36.86-, 36.86-, and 15.22-dB lower noise powers than Non-NLM, PD-NLM, and S-NLM, respectively. In the human spinal cord tumor, the full-width at half-maximums (FWHMs) of vessel cross Section are 76, 201, and [Formula: see text] for St-NLM, Non-NLM, and S-NLM, respectively. The proposed St-NLM method can enhance the microvascular visualization in uPDI and has the potential for the diagnosis of many microvessel-change-related diseases.

Estimating the neovascularity of human finger tendon through high frequency ultrasound micro-Doppler imaging

OBJECTIVE

Neovascularization of injured tendons prolongs the proliferative phase of healing, but prolonged neovascularization may cause improper healing and pain. Currently, ultrasound Doppler imaging is used for measuring the neovascularization of injured tendons (e.g., Achilles tendon). However, the resolution of state-of-the-art clinical ultrasound machines is insufficient for visualizing the neovascularization in finger tendons. In this study, a high-resolution micro-Doppler imaging (HFDI) based on 40-MHz ultrafast ultrasound imaging was proposed for visualizing the neovascularization in injured finger tendons during multiple rehabilitation phases.

METHOD

The vessel visibility was enhanced through a block-wise singular value decomposition filter and several curvilinear structure enhancement strategies, including the bowler-hat transform and Hessian-based vessel enhancement filtering. HFDI was verified through small animal kidney and spleen imaging because the related vessel structure patterns of mice are well studied. Five patients with finger tendon injuries underwent HFDI examination at various rehabilitation phases after surgery (weeks 1156), and finger function evaluations were performed for comparisons.

RESULTS

The results of small animal experiments revealed that the proposed HFDI provides excellent microvasculature imaging performance; the contrast-to-noise ratio of HFDI was approximately 15 dB higher than that of the conventional singular value decomposition filter, and the minimum detectable vessel size for mouse kidney was 35 m without the use of contrast agent. In the human study, neovascularization was clearly observed in injured finger tendons during the early phase of healing (weeks 1121), but it regressed from week 52 to 56. Finger rehabilitation appears to help reduce neovascularization; neovascular density decreased by approximately 1.8%8.0% in participants after 4 weeks of rehabilitation.

CONCLUSION

The experimental results verified the performance of HFDI for microvasculature imaging and its potential for injured finger tendon evaluations.

Volumetric imaging and morphometric analysis of breast tumor angiogenesis using a new contrast-free ultrasound technique: a feasibility study

BACKGROUND

There is a strong correlation between the morphological features of new tumor vessels and malignancy. However, angiogenic heterogeneity necessitates 3D microvascular data of tumor microvessels for more reliable quantification. To provide more accurate information regarding vessel morphological features and improve breast lesion characterization, we introduced a quantitative 3D high-definition microvasculature imaging (q3D-HDMI) as a new easily applicable and robust tool to morphologically characterize microvasculature networks in breast tumors using a contrast-free ultrasound-based imaging approach.

METHODS

In this prospective study, from January 2020 through December 2021, a newly developed q3D-HDMI technique was evaluated on participants with ultrasound-identified suspicious breast lesions recommended for core needle biopsy. The morphological features of breast tumor microvessels were extracted from the q3D-HDMI. Leave-one-out cross-validation (LOOCV) was applied to test the combined diagnostic performance of multiple morphological parameters of breast tumor microvessels. Receiver operating characteristic (ROC) curves were used to evaluate the prediction performance of the generated pooled model.

RESULTS

Ninety-three participants (mean age 52 ± 17 years, 91 women) with 93 breast lesions were studied. The area under the ROC curve (AUC) generated with q3D-HDMI was 95.8% (95% CI 0.901-1.000), yielding a sensitivity of 91.7% and a specificity of 98.2%, that was significantly higher than the AUC generated with the q2D-HDMI (p = 0.02). When compared to q2D-HDMI, the tumor microvessel morphological parameters obtained from q3D-HDMI provides distinctive information that increases accuracy in differentiating breast tumors.

CONCLUSIONS

The proposed quantitative volumetric imaging technique augments conventional breast ultrasound evaluation by increasing specificity in differentiating malignant from benign breast masses.

Hybrid high-definition microvessel imaging/shear wave elastography improves breast lesion characterization

BACKGROUND

Low specificity in current breast imaging modalities leads to increased unnecessary follow-ups and biopsies. The purpose of this study is to evaluate the efficacy of combining the quantitative parameters of high-definition microvasculature imaging (HDMI) and 2D shear wave elastography (SWE) with clinical factors (lesion depth and age) for improving breast lesion differentiation.

METHODS

In this prospective study, from June 2016 through April 2021, patients with breast lesions identified on diagnostic ultrasound and recommended for core needle biopsy were recruited. HDMI and SWE were conducted prior to biopsies. Two new HDMI parameters, Murray's deviation and bifurcation angle, and a new SWE parameter, mass characteristic frequency, were included for quantitative analysis. Lesion malignancy prediction models based on HDMI only, SWE only, the combination of HDMI and SWE, and the combination of HDMI, SWE and clinical factors were trained via elastic net logistic regression with 70% (360/514) randomly selected data and validated with the remaining 30% (154/514) data. Prediction performances in the validation test set were compared across models with respect to area under the ROC curve as well as sensitivity and specificity based on optimized threshold selection.

RESULTS

A total of 508 participants (mean age, 54 years ± 15), including 507 female participants and 1 male participant, with 514 suspicious breast lesions (range, 4-72 mm, median size, 13 mm) were included. Of the lesions, 204 were malignant. The SWE-HDMI prediction model, combining quantitative parameters from SWE and HDMI, with AUC of 0.973 (95% CI 0.95-0.99), was significantly higher than the result predicted with the SWE model or HDMI model alone. With an optimal cutoff of 0.25 for the malignancy probability, the sensitivity and specificity were 95.5% and 89.7%, respectively. The specificity was further improved with the addition of clinical factors. The corresponding model defined as the SWE-HDMI-C prediction model had an AUC of 0.981 (95% CI 0.96-1.00).

CONCLUSIONS

The SWE-HDMI-C detection model, a combination of SWE estimates, HDMI quantitative biomarkers and clinical factors, greatly improved the accuracy in breast lesion characterization.

Quantitative Biomarkers for Cancer Detection Using Contrast-Free Ultrasound High-Definition Microvessel Imaging: Fractal Dimension, Murray's Deviation, Bifurcation Angle & Spatial Vascularity Pattern

A growing body of evidence indicates that there is a strong correlation between microvascular morphological features and malignant tumors. Therefore, quantification of these features might allow more accurate differentiation of benign and malignant tumors. The main objective of this research project is to improve the quantification of microvascular networks depicted in contrast-free ultrasound microvessel images. To achieve this goal, a new series of quantitative microvessel morphological parameters are introduced for differentiation of breast masses using contrast-free ultrasound-based high-definition microvessel imaging (HDMI). Using HDMI, we quantified and analyzed four new parameters: 1) microvessel fractal dimension (mvFD), a marker of tumor microvascular complexity; 2) Murray's deviation (MD), the diameter mismatch, defined as the deviation from Murray's law; 3) bifurcation angle (BA), abnormally decreased angle; and 4) spatial vascular pattern (SVP), indicating tumor vascular distribution pattern, either intratumoral or peritumoral. The new biomarkers have been tested on 60 patients with breast masses. Validation of the feature's extraction algorithm was performed using a synthetic data set. All the proposed parameters had the power to discriminate the breast lesion malignancy (p < 0.05), displaying BA as the most sensitive test, with a sensitivity of 90.6%, and mvFD as the most specific test, with a specificity of 92%. The results of all four new biomarkers showed an AUC = 0.889, sensitivity of 80% and specificity of 91.4% In conclusion, the added value of the proposed quantitative morphological parameters, as new biomarkers of angiogenesis within breast masses, paves the way for more accurate breast cancer detection with higher specificity.

Dynamic Filtering of Adherent and Non-adherent Microbubble Signals Using Singular Value Thresholding and Normalized Singular Spectrum Area Techniques

Ultrasound molecular imaging techniques rely on the separation and identification of three types of signals: static tissue, adherent microbubbles and non-adherent microbubbles. In this study, the image filtering techniques of singular value thresholding (SVT) and normalized singular spectrum area (NSSA) were combined to isolate and identify vascular endothelial growth factor receptor 2-targeted microbubbles in a mouse hindlimb tumor model (n = 24). By use of a Verasonics Vantage 256 imaging system with an L12-5 transducer, a custom-programmed pulse inversion sequence employing synthetic aperture virtual source element imaging was used to collect contrast images of mouse tumors perfused with microbubbles. SVT was used to suppress static tissue signals by 9.6 dB while retaining adherent and non-adherent microbubble signals. NSSA was used to classify microbubble signals as adherent or non-adherent with high accuracy (receiver operating characteristic area under the curve [ROC AUC] = 0.97), matching the classification performance of differential targeted enhancement. The combined SVT + NSSA filtering method also outperformed differential targeted enhancement in differentiating MB signals from all other signals (ROC AUC = 0.89) without necessitating destruction of the contrast agent. The results from this study indicate that SVT and NSSA can be used to automatically segment and classify contrast signals. This filtering method with potential real-time capability could be used in future diagnostic settings to improve workflow and speed the clinical uptake of ultrasound molecular imaging techniques.

The Viability of 3-D Power Doppler Imaging Using Continuous Mechanical Translation: Simulation and Theoretical Analysis

Although conventional Doppler ultrasound is widely used for quantifying blood flow, it is restricted by its low sensitivity to detect slow flow. The incorporation of ultrafast ultrasound and spatial-temporal clutter filters can not only extensively boost the Doppler sensitivity to low-velocity slow flow but also facilitate the development of advanced 3-D Doppler techniques. In this work, we propose a novel 3-D Doppler method which extends 2-D imaging to 3-D through the continuous mechanical translation of a linear transducer. The viability of this method is assessed by simulations with the aids of a theoretical model. The combination of simulations and the theoretical model provides unique insights into the inherent mechanisms involved in the performance of this 3-D Doppler method and the roles of factors, such as tissue vibration characteristics, blood flow velocity, elevational point-spread-function profile, probe translating speed, and signal energy ratios.

Multidimensional Clutter Filtering of Aperture Domain Data for Improved Blood Flow Sensitivity

Singular value decomposition (SVD) is a valuable factorization technique used in clutter rejection filtering for power Doppler imaging. Conventionally, SVD is applied to a Casorati matrix of radio frequency data, which enables filtering based on spatial or temporal characteristics. In this article, we propose a clutter filtering method that uses a higher order SVD (HOSVD) applied to a tensor of aperture data, e.g., delayed channel data. We discuss temporal, spatial, and aperture domain features that can be leveraged in filtering and demonstrate that this multidimensional approach improves sensitivity toward blood flow. Further, we show that HOSVD remains more robust to short ensemble lengths than conventional SVD filtering. Validation of this technique is shown using Field II simulations and in vivo data.

In vivo 3D Power Doppler Imaging Using Continuous Translation and Ultrafast Ultrasound

The introduction of ultrafast ultrasound and spatiotemporal filtering has significantly improved the sensitivity of Doppler ultrasound imaging. This work describes the development of a novel 3D power Doppler imaging technique which uses a 1D-array ultrasound probe that mechanically translates at a constant speed. The continuous translation allows for a fast scan of a large 3D volume without requiring complex hardware. The technique was realized in a prototype and its feasibility illustrated using phantom and in-vivo kidney and breast lesion experiments. Although this 3D Doppler imaging technique is limited in some aspects, it enables power Doppler imaging of a large volume in a short acquisition time with less computational costs.

Quantitative assessment of ensemble coherency in contrast-free ultrasound microvasculature imaging

Purpose

Contrast‐free visualization of microvascular blood flow (MBF) using ultrasound can play a valuable role in diagnosis and detection of diseases. In this study, we demonstrate the importance of quantifying ensemble coherence for robust MBF imaging. We propose a novel approach to quantify ensemble coherence by estimating the local spatiotemporal correlation (LSTC) image, and evaluate its efficacy through simulation and in vivo studies.

Methods

The in vivo patient studies included three volunteers with a suspicious breast tumor, 15 volunteers with a suspicious thyroid tumor, and two healthy volunteers for renal MBF imaging. The breast data displayed negligible prior motion and were used for simulation analysis involving synthetically induced motion, to assess its impact on ensemble coherency and motion artifacts in MBF images. The in vivo thyroid data involved complex physiological motion due to its proximity to the pulsating carotid artery, which was used to assess the in vivo efficacy of the proposed technique. Further, in vivo renal MBF images demonstrated the feasibility of using the proposed ensemble coherence metric for curved array‐based MBF imaging involving phase conversion. All ultrasound data were acquired at high imaging frame rates and the tissue signal was suppressed using spatiotemporal clutter filtering. Thyroid tissue motion was estimated using two‐dimensional normalized cross correlation‐based speckle tracking, which was subsequently used for ensemble motion correction. The coherence of the MBF image was quantified based on Casorati correlation of the Doppler ensemble.

Results

The simulation results demonstrated that an increase in ensemble motion corresponded with a decrease in ensemble coherency, which reciprocally degraded the MBF images. Further the data acquired from breast tumors demonstrated higher ensemble coherency than that from thyroid tumors. Motion correction improved the coherence of the thyroid MBF images, which substantially improved its visualization. The proposed coherence metrics were also useful in assessing the ensemble coherence for renal MBF imaging. The results also demonstrated that the proposed coherence metric can be reliably estimated from downsampled ensembles (by up to 90%), thus allowing improved computational efficiency for potential applications in real‐time MBF imaging.

Conclusions

This pilot study demonstrates the importance of assessing ensemble coherency in contrast‐free MBF imaging. The proposed LSTC image quantified coherence of the Doppler ensemble for robust MBF imaging. The results obtained from this pilot study are promising, and warrant further development and in vivo validation.

Blood Flow Imaging in the Neonatal Brain Using Angular Coherence Power Doppler

Using ultrasound to image small vessels in the neonatal brain can be difficult in the presence of strong clutter from the surrounding tissue and with a neonate motion during the scan. We propose a coherence-based beamforming method, namely the short-lag angular coherence (SLAC) beamforming that suppresses incoherent noise and motion artifacts in Ultrafast data, and we demonstrate its applicability to improve detection of blood flow in the neonatal brain. Instead of estimating spatial coherence across the receive elements, SLAC utilizes the principle of acoustic reciprocity to estimate angular coherence from the beamsummed signals from different plane-wave transmits, which makes it computationally efficient and amenable to advanced beamforming techniques, such as f-k migration. The SLAC images of a simulated speckle phantom show similar edge resolution and texture size as the matching B-mode images, and reduced random noise in the background. We apply SLAC power Doppler (PD) to free-hand imaging of neonatal brain vasculature with long Doppler ensembles and show that: 1) it improves visualization of small vessels in the cortex compared to conventional PD and 2) it can be used for tracking of blood flow in the brain over time, meaning it could potentially improve the quality of free-hand functional ultrasound.

Clutter suppression in ultrasound: performance evaluation and review of low-rank and sparse matrix decomposition methods

Vessel diseases are often accompanied by abnormalities related to vascular shape and size. Therefore, a clear visualization of vasculature is of high clinical significance. Ultrasound color flow imaging (CFI) is one of the prominent techniques for flow visualization. However, clutter signals originating from slow-moving tissue are one of the main obstacles to obtain a clear view of the vascular network. Enhancement of the vasculature by suppressing the clutters is a significant and irreplaceable step for many applications of ultrasound CFI. Currently, this task is often performed by singular value decomposition (SVD) of the data matrix. This approach exhibits two well-known limitations. First, the performance of SVD is sensitive to the proper manual selection of the ranks corresponding to clutter and blood subspaces. Second, SVD is prone to failure in the presence of large random noise in the dataset. A potential solution to these issues is using decomposition into low-rank and sparse matrices (DLSM) framework. SVD is one of the algorithms for solving the minimization problem under the DLSM framework. Many other algorithms under DLSM avoid full SVD and use approximated SVD or SVD-free ideas which may have better performance with higher robustness and less computing time. In practice, these models separate blood from clutter based on the assumption that steady clutter represents a low-rank structure and that the moving blood component is sparse. In this paper, we present a comprehensive review of ultrasound clutter suppression techniques and exploit the feasibility of low-rank and sparse decomposition schemes in ultrasound clutter suppression. We conduct this review study by adapting 106 DLSM algorithms and validating them against simulation, phantom, and in vivo rat datasets. Two conventional quality metrics, signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR), are used for performance evaluation. In addition, computation times required by different algorithms for generating clutter suppressed images are reported. Our extensive analysis shows that the DLSM framework can be successfully applied to ultrasound clutter suppression.

Independent Component-Based Spatiotemporal Clutter Filtering for Slow Flow Ultrasound

Effective tissue clutter filtering is critical for non-contrast ultrasound imaging of slow blood flow in small vessels. Independent component analysis (ICA) has been considered by other groups for ultrasound clutter filtering in the past and was shown to be superior to principal component analysis (PCA)-based methods. However, it has not been considered specifically for slow flow applications or revisited since the onset of other slow flow-focused advancements in beamforming and tissue filtering, namely angled plane wave beamforming and full spatiotemporal singular value decomposition (SVD) (i.e., PCA-based) tissue filtering. In this work, we aim to develop a full spatiotemporal ICA-based tissue filtering technique facilitated by plane wave applications and compare it to SVD filtering. We compare ICA and SVD filtering in terms of optimal image quality in simulations and phantoms as well as in terms of optimal correlation to ground truth blood signal in simulations. Additionally, we propose an adaptive blood independent component sorting and selection method. We show that optimal and adaptive ICA can consistently separate blood from tissue better than principal component analysis (PCA)-based methods using simulations and phantoms. Additionally we demonstrate initial in vivo feasibility in ultrasound data of a liver tumor.

Low Rank and Sparse Decomposition of Ultrasound Color Flow Images for Suppressing Clutter in Real-Time

In this work, a novel technique for real-time clutter rejection in ultrasound Color Flow Imaging (CFI) is proposed. Suppressing undesired clutter signal is important because clutter prohibits an unambiguous view of the vascular network. Although conventional eigen-based filters are potentially efficient in suppressing clutter signal, their performance is highly dependent on proper selection of a clutter to blood boundary which is done manually. Herein, we resolve this limitation by formulating the clutter suppression problem as a foreground-background separation problem to extract the moving blood component. To that end, we adapt the fast Robust Matrix Completion (fRMC) algorithm, and utilize the in-face extended Frank-Wolfe method to minimize the rank of the matrix of ultrasound frames. Our method automates the clutter suppression process, which is critical for clinical use. We name the method RAPID (Robust mAtrix decomPosition for suppressIng clutter in ultrasounD) since the automation step can substantially streamline clutter suppression. The technique is validated with simulation, flow phantom and two sets of in-vivo data. RAPID code as well as most of the data in this paper can be downloaded from RAPID.sonography.ai.

Flexible Body-Conformal Ultrasound Patches for Image-Guided Neuromodulation

The paper presents the design and validation of body-conformal active ultrasound patches with integrated imaging and modulation modalities for image-guided neural therapy. A mechanically-flexible linear 64-element array of piezoelectric transducers with a resonance frequency of 5 MHz was designed for nerve localization. A second 8-element array using larger elements was integrated on the wearable probe for low intensity focused ultrasound neuromodulation at a resonance frequency of 1.3 MHz. Full-wave simulations were used to model the flexible arrays and estimate their generated pressure profiles. A focal depth of 10-20 mm was assumed for beamforming and focusing to support modulation of the vagus, tibial, and other nerves. A strain sensor integrated on the probe provides patient-specific feedback information on array curvature for real-time optimization of focusing and image processing. Each patch also includes high voltage (HV) multiplexers, transmit/receive switches, and pre-amplifiers that simplify connectivity and also improve the signal-to-noise ratio (SNR) of the received echo signals by  ∼ 5 dB. Experimental results from a flexible prototype show a sensitivity of 80 kPa/V with  ∼ 3 MHz bandwidth for the modulation and 20 kPa/V with  ∼ 6 MHz bandwidth for the imaging array. An algorithm for accurate and automatic localization of targeted nerves based on using nearby blood vessels (e.g., the carotid artery) as image markers is also presented.

Impact of imaging cross-section on visualization of thyroid microvessels using ultrasound: Pilot study

Non-invasive, contrast-free microvascular imaging of human thyroids can be potentially beneficial in reducing the large number of benign biopsies of suspicious nodules. However, motion incurred by thyroid due to its proximity to the pulsating carotid artery significantly impacts the visualization of blood flow in small vessels. Singular value based spatiotemporal clutter filtering (SVD-STF) improves the performance of tissue rejection in the presence of motion. However, despite effective clutter filtering, motion in thyroid imaging can impact coherent integration of the Doppler ensemble and degrade the visualization of the underlying vasculature. Recently studies have demonstrated that motion correction using 2D normalized cross-correlation based speckle tracking can address this issue, however, only in-plane motion can be tracked and corrected. Given the natural anatomical orientation of the rigid trachea, thyroid and the pulsating carotid artery, we hypothesize that imaging of thyroid microvessels may be more reliable in the longitudinal view than in the transverse. Specifically, distal presence of rigid trachea can limit out-of-plane motion in the longitudinal view. We tested this hypothesis on 48 acquisitions obtained from 24 thyroid patients having at least one suspicious nodule. In each patient, ultrasound images of the thyroid were acquired in both longitudinal and transverse views. Compounded plane-wave imaging was used to acquire the ultrasound images at high frame-rate, which is important for contrast-free small vessel blood flow imaging. Thyroid motion was tracked using 2D normalized cross-correlation based speckle tracking. Tissue clutter was rejected using singular value decomposition based spatiotemporal clutter filtering. The clutter-filtered Doppler ensemble was motion corrected prior to slow-time power Doppler integration. Signal-to-noise and contrast-to-noise ratios were computed to assess the improvement in quality of the power Doppler images. Out-of-plane motion was detected by estimating normalized ensemble cross-correlation coefficient. The results demonstrated that motion associated with the thyroid due to the carotid artery was primarily in the lateral direction, which could be estimated and corrected using 2D speckle tracking. However, the motion in the transverse view displayed increased speckle decorrelation. The average ensemble cross-correlation coefficient of the thyroid ultrasound images were significantly higher (p < 0.05) in the longitudinal view than in the transverse view. The largest improvement in SNR and CNR of the estimated PD images upon motion correction was observed in the longitudinal view (12.95 ± 3.76 dB and 16.48 ± 4.6 dB) than in the transverse view (3.72 ± 0.894 dB and 6.217 ± 1.689 dB). These preliminary results show that motion encountered by the thyroid due to carotid pulsations can be effectively tracked and corrected in the longitudinal view relative to transverse, which is important for reliably visualizing the underlying blood flow.

Quantification of Morphological Features in Non-Contrast-Enhanced Ultrasound Microvasculature Imaging

There are significant differences in microvascular morphological features in diseased tissues, such as cancerous lesions, compared to noncancerous tissue. Quantification of microvessel morphological features could play an important role in disease diagnosis and tumor classification. However, analyzing microvessel morphology in ultrasound Doppler is a challenging task due to limitations associated with this technique. Our main objective is to provide methods for quantifying morphological features of microvasculature obtained by ultrasound Doppler imaging. To achieve this goal, we propose multiple image enhancement techniques and appropriate morphological feature extraction methods that enable quantitative analysis of microvasculature structures. Vessel segments obtained by the skeletonization of the regularized microvasculature images are further analyzed to satisfy other constraints, such as vessel segment diameter and length. Measurements of some morphological metrics, such as tortuosity, depend on preserving large vessel trunks. To address this issue, additional filtering methods are proposed. These methods are tested on in vivo images of breast lesion and thyroid nodule microvasculature, and the outcomes are discussed. Initial results show that using vessel morphological features allows for differentiation between malignant and benign breast lesions (p-value < 0.005) and thyroid nodules (p-value < 0.01). This paper provides a tool for the quantification of microvasculature images obtained by non-contrast ultrasound imaging, which may serve as potential biomarkers for the diagnosis of some diseases.

Adaptive background noise bias suppression in contrast-free ultrasound microvascular imaging

Non-invasive, contrast-free imaging of small vessel blood flow is diagnostically invaluable for detection, diagnosis and monitoring of disease. Recent advances in ultrafast imaging and tissue clutter-filtering have considerably improved the sensitivity of power Doppler (PD) imaging in detecting small vessel blood flow. However, suppression of tissue clutter exposes the depth-dependent time-gain compensated noise bias that noticeably degrades the PD image. We hypothesized that background suppression of PD images based on noise bias estimated from the entire clutter-filtered singular value spectrum can considerably improve flow signal visualization compared to currently existing techniques. To test our hypothesis, in vivo experiments were conducted on suspicious breast lesions in 10 subjects and deep-seated hepatic and renal microvasculatures in four healthy volunteers. Ultrasound PD images were acquired using a clinical ultrasound scanner, implemented with compounded plane wave imaging. The time gain compensated noise field was computed from the clutter-filtered Doppler ensemble (CFDE) based on its local spatio-temporal correlation, combined with low-rank signal estimation. Subsequently, the background bias in the PD images was suppressed by subtracting the estimated noise field. Background-suppressed PD images obtained using the proposed technique substantially improved visualization of the blood flow signal. The background bias in the noise suppressed PD images varied  <0.6 dB, independent of depth, which otherwise increased up to 13.8 dB. Further, the results demonstrated that the proposed technique efficaciously suppressed the background noise bias associated with smaller Doppler ensembles, which are challenging due to increased overlap between blood flow and noise components in the singular value spectrum. These preliminary results demonstrate the utility of the proposed technique to improve the visualization of small vessel blood flow in contrast-free PD images. The results of this feasibility study were encouraging, and warrant further development and additional in vivo validation.

3D Super-Resolution US Imaging of Rabbit Lymph Node Vasculature in Vivo by Using Microbubbles

Background Variations in lymph node (LN) microcirculation can be indicative of metastasis. The identification and quantification of metastatic LNs remains essential for prognosis and treatment planning, but a reliable noninvasive imaging technique is lacking. Three-dimensional super-resolution (SR) US has shown potential to noninvasively visualize microvascular networks in vivo. Purpose To study the feasibility of three-dimensional SR US imaging of rabbit LN microvascular structure and blood flow by using microbubbles. Materials and Methods In vivo studies were carried out to image popliteal LNs of two healthy male New Zealand white rabbits aged 6-8 weeks. Three-dimensional, high-frame-rate, contrast material-enhanced US was achieved by mechanically scanning with a linear imaging probe. Individual microbubbles were identified, localized, and tracked to form three-dimensional SR images and super-resolved velocity maps. Acoustic subaperture processing was used to improve image contrast and to generate enhanced power Doppler and color Doppler images. Vessel size and blood flow velocity distributions were evaluated and assessed by using Student paired t test. Results SR images revealed microvessels in the rabbit LN, with branches clearly resolved when separated by 30 µm, which is less than half of the acoustic wavelength and not resolvable by using power or color Doppler. The apparent size distribution of most vessels in the SR images was below 80 µm and agrees with micro-CT data, whereas most of those detected with Doppler techniques were larger than 80 µm in the images. The blood flow velocity distribution indicated that most of the blood flow in rabbit popliteal LN was at velocities lower than 5 mm/sec. Conclusion Three-dimensional super-resolution US imaging using microbubbles allows noninvasive nonionizing visualization and quantification of lymph node microvascular structures and blood flow dynamics with resolution below the wave diffraction limit. This technology has potential for studying the physiologic functions of the lymph system and for clinical detection of lymph node metastasis. Published under a CC BY 4.0 license. Online supplemental material is available for this article.

Non-Local Based Denoising Framework for In Vivo Contrast-Free Ultrasound Microvessel Imaging

Vascular networks can provide invaluable information about tumor angiogenesis. Ultrafast Doppler imaging enables ultrasound to image microvessels by applying tissue clutter filtering methods on the spatio-temporal data obtained from plane-wave imaging. However, the resultant vessel images suffer from background noise that degrades image quality and restricts vessel visibilities. In this paper, we addressed microvessel visualization and the associated noise problem in the power Doppler images with the goal of achieving enhanced vessel-background separation. We proposed a combination of patch-based non-local mean filtering and top-hat morphological filtering to improve vessel outline and background noise suppression. We tested the proposed method on a flow phantom, as well as in vivo breast lesions, thyroid nodules, and pathologic liver from human subjects. The proposed non-local-based framework provided a remarkable gain of more than 15 dB, on average, in terms of contrast-to-noise and signal-to-noise ratios. In addition to improving visualization of microvessels, the proposed method provided high quality images suitable for microvessel morphology quantification that may be used for diagnostic applications.