In vivo acoustic super-resolution and super-resolved velocity mapping using microbubbles

Overcoming the acoustic diffraction limit in photoacoustic imaging by the localization of flowing absorbers

The resolution of photoacoustic imaging deep inside scattering media is limited by the acoustic diffraction limit. In this Letter, taking inspiration from super-resolution imaging techniques developed to beat the optical diffraction limit, we demonstrate that the localization of individual optical absorbers can provide super-resolution photoacoustic imaging well beyond the acoustic diffraction limit. As a proof-of-principle experiment, photoacoustic cross-sectional images of microfluidic channels were obtained with a 15 MHz linear capacitive micromachined ultrasonic transducer array, while absorbing beads were flown through the channels. The localization of individual absorbers allowed us to obtain a super-resolved cross-sectional image of the channels by reconstructing both the channel width and position with an accuracy better than λ/10. Given the discrete nature of endogenous absorbers such as red blood cells, or that of exogenous particular contrast agents, localization is a promising approach to push the current resolution limits of photoacoustic imaging.

Microbubble Axial Localization Errors in Ultrasound Super-Resolution Imaging

Acoustic super-resolution imaging has allowed the visualization of microvascular structure and flow beyond the diffraction limit using standard clinical ultrasound systems through the localization of many spatially isolated microbubble signals. The determination of each microbubble position is typically performed by calculating the centroid, finding a local maximum, or finding the peak of a 2-D Gaussian function fit to the signal. However, the backscattered signal from a microbubble depends not only on diffraction characteristics of the waveform, but also on the microbubble behavior in the acoustic field. Here, we propose a new axial localization method by identifying the onset of the backscattered signal. We compare the accuracy of localization methods using in vitro experiments performed at 7-cm depth and 2.3-MHz center frequency. We corroborate these findings with simulation results based on the Marmottant model. We show experimentally and in simulations that detecting the onset of the returning signal provides considerably increased accuracy for super-resolution. Resulting experimental cross-sectional profiles in super-resolution images demonstrate at least 5.8 times improvement in contrast ratio and more than 1.8 times reduction in spatial spread (provided by 90% of the localizations) for the onset method over centroiding, peak detection, and 2-D Gaussian fitting methods. Simulations estimate that these latter methods could create errors in relative bubble positions as high as at these experimental settings, while the onset method reduced the interquartile range of these errors by a factor of over 2.2. Detecting the signal onset is, therefore, expected to considerably improve the accuracy of super-resolution.

Toward optimization of in vivo super-resolution ultrasound imaging using size-selected microbubble contrast agents

PURPOSE

Microvascular processes play key roles in many diseases including diabetes. Improved understanding of the microvascular changes involved in disease development could offer crucial insight into the relationship of these changes to disease pathogenesis. Super-resolution ultrasound (SR-US) imaging has showed the potential to visualize microvascular detail down to the capillary level (i.e., subwavelength resolution), but optimization is still necessary. The purpose of this study was to investigate in vivo SR-US imaging of skeletal muscle microvascularity using microbubble (MB) contrast agents of various size and concentration while evaluating different ultrasound (US) system level parameters such as imaging frame rate and image acquisition length.

METHODS

An US system equipped with a linear array transducer was used in a harmonic imaging mode at low transmit power. C57BL/6J mice fed a normal diet were used in this study. An assortment of size-selected MB contrast agents (1-2 μm, 3-4 μm, and 5-8 μm in diameter) were slowly infused in the tail vein at various doses (1.25 × 10⁷ , 2.5 × 10⁷ , or 5 × 10⁷  MBs). US image data were collected before MB injection and thereafter for 10 min at 30 frames per s (fps). The US transducer was fixed throughout and between each imaging period to help capture microvascular patterns along the same image plane. An adaptive SR-US image processing technique was implemented using custom Matlab software.

RESULTS

Experimental findings illustrate the use of larger MB results in better SR-US images in terms of skeletal muscle microvascular detail. A dose of 2.5 × 10⁷  MBs resulted in SR-US images with optimal spatial resolution. An US imaging rate of at least 20 fps and image acquisition length of at least 8 min also resulted in SR-US images with pronounced microvascular detail.

CONCLUSIONS

This study indicates that MB size and dose and US system imaging rate and data acquisition length have significant impact on the quality of in vivo SR-US images of skeletal muscle microvascularity.

Ultrasound localization microscopy to image and assess microvasculature in a rat kidney

The recent development of ultrasound localization microscopy, where individual microbubbles (contrast agents) are detected and tracked within the vasculature, provides new opportunities for imaging the vasculature of entire organs with a spatial resolution below the diffraction limit. In stationary tissue, recent studies have demonstrated a theoretical resolution on the order of microns. In this work, single microbubbles were localized in vivo in a rat kidney using a dedicated high frame rate imaging sequence. Organ motion was tracked by assuming rigid motion (translation and rotation) and appropriate correction was applied. In contrast to previous work, coherence-based non-linear phase inversion processing was used to reject tissue echoes while maintaining echoes from very slowly moving microbubbles. Blood velocity in the small vessels was estimated by tracking microbubbles, demonstrating the potential of this technique to improve vascular characterization. Previous optical studies of microbubbles in vessels of approximately 20 microns have shown that expansion is constrained, suggesting that microbubble echoes would be difficult to detect in such regions. We therefore utilized the echoes from individual MBs as microscopic sensors of slow flow associated with such vessels and demonstrate that highly correlated, wideband echoes are detected from individual microbubbles in vessels with flow rates below 2 mm/s.

Neutrophil-inspired propulsion in a combined acoustic and magnetic field

Systems capable of precise motion in the vasculature can offer exciting possibilities for applications in targeted therapeutics and non-invasive surgery. So far, the majority of the work analysed propulsion in a two-dimensional setting with limited controllability near boundaries. Here we show bio-inspired rolling motion by introducing superparamagnetic particles in magnetic and acoustic fields, inspired by a neutrophil rolling on a wall. The particles self-assemble due to dipole–dipole interaction in the presence of a rotating magnetic field. The aggregate migrates towards the wall of the channel due to the radiation force of an acoustic field. By combining both fields, we achieved a rolling-type motion along the boundaries. The use of both acoustic and magnetic fields has matured in clinical settings. The combination of both fields is capable of overcoming the limitations encountered by single actuation techniques. We believe our method will have far-reaching implications in targeted therapeutics.

Devising effective swimming and propulsion strategies in microenvironments is attractive for drug delivery applications. Here Ahmed et al. demonstrate a micropropulsion strategy in which a combination of magnetic and acoustic fields is used to assemble and propel colloidal particles along channel walls.

3-D In Vitro Acoustic Super-Resolution and Super-Resolved Velocity Mapping Using Microbubbles

Standard clinical ultrasound (US) imaging frequencies are unable to resolve microvascular structures due to the fundamental diffraction limit of US waves. Recent demonstrations of 2-D super-resolution both in vitro and in vivo have demonstrated that fine vascular structures can be visualized using acoustic single bubble localization. Visualization of more complex and disordered 3-D vasculature, such as that of a tumor, requires an acquisition strategy which can additionally localize bubbles in the elevational plane with high precision in order to generate super-resolution in all three dimensions. Furthermore, a particular challenge lies in the need to provide this level of visualization with minimal acquisition time. In this paper, we develop a fast, coherent US imaging tool for microbubble localization in 3-D using a pair of US transducers positioned at 90°. This allowed detection of point scatterer signals in 3-D with average precisions equal to [Formula: see text] in axial and elevational planes, and [Formula: see text] in the lateral plane, compared to the diffraction limited point spread function full-widths at half-maximum of 488, 1188, and [Formula: see text] of the original imaging system with a single transducer. Visualization and velocity mapping of 3-D in vitro structures was demonstrated far beyond the diffraction limit. The capability to measure the complete flow pattern of blood vessels associated with disease at depth would ultimately enable analysis of in vivo microvascular morphology, blood flow dynamics, and occlusions resulting from disease states.

Optimizing Sensitivity of Ultrasound Contrast-Enhanced Super-Resolution Imaging by Tailoring Size Distribution of Microbubble Contrast Agent

Ultrasound contrast-enhanced super-resolution imaging has recently attracted attention because of its extraordinary ability to image vascular features much smaller than the ultrasound diffraction limit. This method requires sensitive detection of separable microbubble events despite a noisy tissue background to indicate the microvasculature, and any approach that could improve the sensitivity of the ultrasound system to individual microbubbles would be highly beneficial. In this study, we evaluated the effect of varying microbubble size on super-resolution imaging sensitivity. Microbubble preparations were size sorted into different mean diameters and then were imaged at equal concentrations. Commercially manufactured Definity and Optison were also imaged for comparison. Both in vitro experiments in phantom vessels and in vivo experiments imaging rat tumors revealed that the sensitivity of contrast-enhanced super-resolution imaging can be improved by using microbubbles with a larger diameter.

A Temporal and Spatial Analysis Approach to Automated Segmentation of Microbubble Signals in Contrast-Enhanced Ultrasound Images: Application to Quantification of Active Vascular Density in Human Lower Limbs

Contrast-enhanced ultrasound (CEUS) using microbubble contrast agents has shown great promise in visualising and quantifying active vascular density. Most existing approaches for vascular density quantification using CEUS are calculated based on image intensity and are susceptible to confounding factors and imaging artefact. Poor reproducibility is a key challenge to clinical translation. In this study, a new automated temporal and spatial signal analysis approach is developed for reproducible microbubble segmentation and quantification of contrast enhancement in human lower limbs. The approach is evaluated in vitro on phantoms and in vivo in lower limbs of healthy volunteers before and after physical exercise. In this approach, vascular density is quantified based on the relative areas microbubbles occupy instead of their image intensity. Temporal features of the CEUS image sequences are used to identify pixels that contain microbubble signals. A microbubble track density (MTD) measure, the ratio of the segmented microbubble area to the whole tissue area, is calculated as a surrogate for active capillary density. In vitro results reveal a good correlation (r² = 0.89) between the calculated MTD measure and the known bubble concentration. For in vivo results, a significant increase (129% in average) in the MTD measure is found in lower limbs of healthy volunteers after exercise, with excellent repeatability over a series of days (intra-class correlation coefficient = 0.96). This compares to the existing state-of-the-art approach of destruction and replenishment analysis on the same patients (intra-class correlation coefficient ≤0.78). The proposed new approach shows great potential as an accurate and highly reproducible clinical tool for quantification of active vascular density.

Experimental performance assessment of the sub-band minimum variance beamformer for ultrasound imaging

Recent progress in adaptive beamforming techniques for medical ultrasound has shown that current resolution limits can be surpassed. One method of obtaining improved lateral resolution is the Minimum Variance (MV) beamformer. The frequency domain implementation of this method effectively divides the broadband ultrasound signals into sub-bands (MVS) to conform with the narrow-band assumption of the original MV theory. This approach is investigated here using experimental Synthetic Aperture (SA) data from wire and cyst phantoms. A 7MHz linear array transducer is used with the SARUS experimental ultrasound scanner for the data acquisition. The lateral resolution and the contrast obtained, are evaluated and compared with those from the conventional Delay-and-Sum (DAS) beamformer and the MV temporal implementation (MVT). From the wire phantom the Full-Width-at-Half-Maximum (FWHM) measured at a depth of 52mm, is 16.7μm (0.08λ) for both MV methods, while the corresponding values for the DAS case are at least 24 times higher. The measured Peak-Side-lobe-Level (PSL) may reach -41dB using the MVS approach, while the values from the DAS and MVT beamforming are above -24dB and -33dB, respectively. From the cyst phantom, the power ratio (PR), the contrast-to-noise ratio (CNR), and the speckle signal-to-noise ratio (sSNR) measured at a depth of 30mm are at best similar for MVS and DAS, with values ranging between -29dB and -30dB, 1.94 and 2.05, and 2.16 and 2.27 respectively. In conclusion the MVS beamformer is not suitable for imaging continuous targets, and significant resolution gains were obtained only for isolated targets.

Subwavelength motion-correction for ultrafast ultrasound localization microscopy

Ultrafast Ultrasound Localization Microscopy uses microbubbles that are individually localized with a resolution below 10μm. Positions of the microbubbles are accumulated to create a super resolution image, which bypass the diffraction-limit of spatial resolution. However, microbubbles localization is affected by physiological motions at the micrometric scale. Here, we demonstrate a phase correlation method for rigid motion correction. Spatiotemporal filters extract tissue dominated images, which are tracked to correct linear motions and improve the precision of microbubbles' localization, improving the quality of the image. It is the first proof of concept towards a full motion correction strategy and super-resolution imaging in moving tissues.

Acoustic window planning for ultrasound acquisition

Autonomous robotic ultrasound has recently gained considerable interest, especially for collaborative applications. Existing methods for acquisition trajectory planning are solely based on geometrical considerations, such as the pose of the transducer with respect to the patient surface.

PURPOSE

This work aims at establishing acoustic window planning to enable autonomous ultrasound acquisitions of anatomies with restricted acoustic windows, such as the liver or the heart.

METHODS

We propose a fully automatic approach for the planning of acquisition trajectories, which only requires information about the target region as well as existing tomographic imaging data, such as X-ray computed tomography. The framework integrates both geometrical and physics-based constraints to estimate the best ultrasound acquisition trajectories with respect to the available acoustic windows. We evaluate the developed method using virtual planning scenarios based on real patient data as well as for real robotic ultrasound acquisitions on a tissue-mimicking phantom.

RESULTS

The proposed method yields superior image quality in comparison with a naive planning approach, while maintaining the necessary coverage of the target.

CONCLUSION

We demonstrate that by taking image formation properties into account acquisition planning methods can outperform naive plannings. Furthermore, we show the need for such planning techniques, since naive approaches are not sufficient as they do not take the expected image quality into account.

Fast Vascular Ultrasound Imaging With Enhanced Spatial Resolution and Background Rejection

Ultrasound super-localization microscopy techniques presented in the last few years enable non-invasive imaging of vascular structures at the capillary level by tracking the flow of ultrasound contrast agents (gas microbubbles). However, these techniques are currently limited by low temporal resolution and long acquisition times. Super-resolution optical fluctuation imaging (SOFI) is a fluorescence microscopy technique enabling sub-diffraction limit imaging with high temporal resolution by calculating high order statistics of the fluctuating optical signal. The aim of this work is to achieve fast acoustic imaging with enhanced resolution by applying the tools used in SOFI to contrast-enhance ultrasound (CEUS) plane-wave scans. The proposed method was tested using numerical simulations and evaluated using two in-vivo rabbit models: scans of healthy kidneys and VX-2 tumor xenografts. Improved spatial resolution was observed with a reduction of up to 50% in the full width half max of the point spread function. In addition, substantial reduction in the background level was achieved compared to standard mean amplitude persistence images, revealing small vascular structures within tumors. The scan duration of the proposed method is less than a second while current super-localization techniques require acquisition duration of several minutes. As a result, the proposed technique may be used to obtain scans with enhanced spatial resolution and high temporal resolution, facilitating flow-dynamics monitoring. Our method can also be applied during a breath-hold, reducing the sensitivity to motion artifacts.

Characterization of Contrast Agent Microbubbles for Ultrasound Imaging and Therapy Research

The high efficiency with which gas microbubbles can scatter ultrasound compared with the surrounding blood pool or tissues has led to their widespread employment as contrast agents in ultrasound imaging. In recent years, their applications have been extended to include super-resolution imaging and the stimulation of localized bio-effects for therapy. The growing exploitation of contrast agents in ultrasound and in particular these recent developments have amplified the need to characterize and fully understand microbubble behavior. The aim in doing so is to more fully exploit their utility for both diagnostic imaging and potential future therapeutic applications. This paper presents the key characteristics of microbubbles that determine their efficacy in diagnostic and therapeutic applications and the corresponding techniques for their measurement. In each case, we have presented information regarding the methods available and their respective strengths and limitations, with the aim of presenting information relevant to the selection of appropriate characterization methods. First, we examine methods for determining the physical properties of microbubble suspensions and then techniques for acoustic characterization of both suspensions and single microbubbles. The next section covers characterization of microbubbles as therapeutic agents, including as drug carriers for which detailed understanding of their surface characteristics and drug loading capacity is required. Finally, we discuss the attempts that have been made to allow comparison across the methods employed by various groups to characterize and describe their microbubble suspensions and promote wider discussion and comparison of microbubble behavior.

3-D Ultrasound Localization Microscopy for Identifying Microvascular Morphology Features of Tumor Angiogenesis at a Resolution Beyond the Diffraction Limit of Conventional Ultrasound

Angiogenesis has been known as a hallmark of solid tumor cancers for decades, yet ultrasound has been limited in its ability to detect the microvascular changes associated with malignancy. Here, we demonstrate the potential of 'ultrasound localization microscopy' applied volumetrically in combination with quantitative analysis of microvascular morphology, as an approach to overcome this limitation. This pilot study demonstrates our ability to image complex microvascular patterns associated with tumor angiogenesis in-vivo at a resolution of tens of microns - substantially better than the diffraction limit of traditional clinical ultrasound, yet using an 8 MHz clinical ultrasound probe. Furthermore, it is observed that data from healthy and tumor-bearing tissue exhibit significant differences in microvascular pattern and density. Results suggests that with continued development of these novel technologies, ultrasound has the potential to detect biomarkers of cancer based on the microvascular 'fingerprint' of malignant angiogenesis rather than through imaging of blood flow dynamics or the tumor mass itself.

High Resolution Ultrasound Superharmonic Perfusion Imaging: In Vivo Feasibility and Quantification of Dynamic Contrast-Enhanced Acoustic Angiography

Mapping blood perfusion quantitatively allows localization of abnormal physiology and can improve understanding of disease progression. Dynamic contrast-enhanced ultrasound is a low-cost, real-time technique for imaging perfusion dynamics with microbubble contrast agents. Previously, we have demonstrated another contrast agent-specific ultrasound imaging technique, acoustic angiography, which forms static anatomical images of the superharmonic signal produced by microbubbles. In this work, we seek to determine whether acoustic angiography can be utilized for high resolution perfusion imaging in vivo by examining the effect of acquisition rate on superharmonic imaging at low flow rates and demonstrating the feasibility of dynamic contrast-enhanced superharmonic perfusion imaging for the first time. Results in the chorioallantoic membrane model indicate that frame rate and frame averaging do not affect the measured diameter of individual vessels observed, but that frame rate does influence the detection of vessels near and below the resolution limit. The highest number of resolvable vessels was observed at an intermediate frame rate of 3 Hz using a mechanically-steered prototype transducer. We also demonstrate the feasibility of quantitatively mapping perfusion rate in 2D in a mouse model with spatial resolution of ~100 μm. This type of imaging could provide non-invasive, high resolution quantification of microvascular function at penetration depths of several centimeters.

Towards Dynamic Contrast Specific Ultrasound Tomography

We report on the first study demonstrating the ability of a recently-developed, contrast-enhanced, ultrasound imaging method, referred to as cumulative phase delay imaging (CPDI), to image and quantify ultrasound contrast agent (UCA) kinetics. Unlike standard ultrasound tomography, which exploits changes in speed of sound and attenuation, CPDI is based on a marker specific to UCAs, thus enabling dynamic contrast-specific ultrasound tomography (DCS-UST). For breast imaging, DCS-UST will lead to a more practical, faster, and less operator-dependent imaging procedure compared to standard echo-contrast, while preserving accurate imaging of contrast kinetics. Moreover, a linear relation between CPD values and ultrasound second-harmonic intensity was measured (coefficient of determination = 0.87). DCS-UST can find clinical applications as a diagnostic method for breast cancer localization, adding important features to multi-parametric ultrasound tomography of the breast.

Modeling of Cerebral Oxygen Transport Based on In vivo Microscopic Imaging of Microvascular Network Structure, Blood Flow, and Oxygenation

Oxygen is delivered to brain tissue by a dense network of microvessels, which actively control cerebral blood flow (CBF) through vasodilation and contraction in response to changing levels of neural activity. Understanding these network-level processes is immediately relevant for (1) interpretation of functional Magnetic Resonance Imaging (fMRI) signals, and (2) investigation of neurological diseases in which a deterioration of neurovascular and neuro-metabolic physiology contributes to motor and cognitive decline. Experimental data on the structure, flow and oxygen levels of microvascular networks are needed, together with theoretical methods to integrate this information and predict physiologically relevant properties that are not directly measurable. Recent progress in optical imaging technologies for high-resolution in vivo measurement of the cerebral microvascular architecture, blood flow, and oxygenation enables construction of detailed computational models of cerebral hemodynamics and oxygen transport based on realistic three-dimensional microvascular networks. In this article, we review state-of-the-art optical microscopy technologies for quantitative in vivo imaging of cerebral microvascular structure, blood flow and oxygenation, and theoretical methods that utilize such data to generate spatially resolved models for blood flow and oxygen transport. These “bottom-up” models are essential for the understanding of the processes governing brain oxygenation in normal and disease states and for eventual translation of the lessons learned from animal studies to humans.

Review of ultrasound image guidance in external beam radiotherapy part II: intra-fraction motion management and novel applications

Imaging has become an essential tool in modern radiotherapy (RT), being used to plan dose delivery prior to treatment and verify target position before and during treatment. Ultrasound (US) imaging is cost-effective in providing excellent contrast at high resolution for depicting soft tissue targets apart from those shielded by the lungs or cranium. As a result, it is increasingly used in RT setup verification for the measurement of inter-fraction motion, the subject of Part I of this review (Fontanarosa et al 2015 Phys. Med. Biol. 60 R77-114). The combination of rapid imaging and zero ionising radiation dose makes US highly suitable for estimating intra-fraction motion. The current paper (Part II of the review) covers this topic. The basic technology for US motion estimation, and its current clinical application to the prostate, is described here, along with recent developments in robust motion-estimation algorithms, and three dimensional (3D) imaging. Together, these are likely to drive an increase in the number of future clinical studies and the range of cancer sites in which US motion management is applied. Also reviewed are selections of existing and proposed novel applications of US imaging to RT. These are driven by exciting developments in structural, functional and molecular US imaging and analytical techniques such as backscatter tissue analysis, elastography, photoacoustography, contrast-specific imaging, dynamic contrast analysis, microvascular and super-resolution imaging, and targeted microbubbles. Such techniques show promise for predicting and measuring the outcome of RT, quantifying normal tissue toxicity, improving tumour definition and defining a biological target volume that describes radiation sensitive regions of the tumour. US offers easy, low cost and efficient integration of these techniques into the RT workflow. US contrast technology also has potential to be used actively to assist RT by manipulating the tumour cell environment and by improving the delivery of radiosensitising agents. Finally, US imaging offers various ways to measure dose in 3D. If technical problems can be overcome, these hold potential for wide-dissemination of cost-effective pre-treatment dose verification and in vivo dose monitoring methods. It is concluded that US imaging could eventually contribute to all aspects of the RT workflow.

Super-Resolution Ultrasound Imaging in Vivo with Transient Laser-Activated Nanodroplets

We have developed a method for super-resolution ultrasound imaging, which relies on a new class of blinking nanometer-size contrast agents: laser-activated nanodroplets (LANDs). The LANDs can be repeatedly optically triggered to undergo vaporization; the resulting spatially stationary, temporally transient microbubbles provide high ultrasound contrast for several to hundreds of milliseconds before recondensing to their native liquid nanodroplet state. By capturing high frame rate ultrasound images of blinking LANDs, we demonstrate the ability to detect individual recondensation events. Then we apply a newly developed super-resolution image processing algorithm to localize the LAND positions in vivo almost an order of magnitude better than conventional ultrasound imaging. These results pave the way for high resolution molecular imaging deep in tissue.

Experimental demonstration of passive acoustic imaging in the human skull cavity using CT-based aberration corrections

PURPOSE

Experimentally verify a previously described technique for performing passive acoustic imaging through an intact human skull using noninvasive, computed tomography (CT)-based aberration corrections Jones et al. [Phys. Med. Biol. 58, 4981-5005 (2013)].

METHODS

A sparse hemispherical receiver array (30 cm diameter) consisting of 128 piezoceramic discs (2.5 mm diameter, 612 kHz center frequency) was used to passively listen through ex vivo human skullcaps (n = 4) to acoustic emissions from a narrow-band fixed source (1 mm diameter, 516 kHz center frequency) and from ultrasound-stimulated (5 cycle bursts, 1 Hz pulse repetition frequency, estimated in situ peak negative pressure 0.11-0.33 MPa, 306 kHz driving frequency) Definity™ microbubbles flowing through a thin-walled tube phantom. Initial in vivo feasibility testing of the method was performed. The performance of the method was assessed through comparisons to images generated without skull corrections, with invasive source-based corrections, and with water-path control images.

RESULTS

For source locations at least 25 mm from the inner skull surface, the modified reconstruction algorithm successfully restored a single focus within the skull cavity at a location within 1.25 mm from the true position of the narrow-band source. The results obtained from imaging single bubbles are in good agreement with numerical simulations of point source emitters and the authors' previous experimental measurements using source-based skull corrections O'Reilly et al. [IEEE Trans. Biomed. Eng. 61, 1285-1294 (2014)]. In a rat model, microbubble activity was mapped through an intact human skull at pressure levels below and above the threshold for focused ultrasound-induced blood-brain barrier opening. During bursts that led to coherent bubble activity, the location of maximum intensity in images generated with CT-based skull corrections was found to deviate by less than 1 mm, on average, from the position obtained using source-based corrections.

CONCLUSIONS

Taken together, these results demonstrate the feasibility of using the method to guide bubble-mediated ultrasound therapies in the brain. The technique may also have application in ultrasound-based cerebral angiography.

Detection and Tracking of Multiple Microbubbles in Ultrasound B-Mode Images

The imaging of microvessels and the quantification of their blood flow is of particular interest in the characterization of tumor vasculature. The imaging resolution (50-200 μm) of high-frequency ultrasound (US) (20-50 MHz) is not sufficient to image microvessels (~10 μm) and Doppler sensitivity is not high enough to measure capillary blood flow (~1 mm/s). For imaging of blood flow in microvessels, our approach is to detect single microbubbles (MBs), track them over several frames, and to estimate their velocity. First, positions of MBs will be detected by separating B-mode frames in a moving foreground and a static background. For the crucial task of association of these positions to tracks, we implemented a modified Markov chain Monte Carlo data association (MCMCDA) algorithm, which can handle a high number of MBs. False alarms, the detection, initiation, and termination of MBs tracks are incorporated in the underlying model. To test the performance of algorithms, a US imaging simulation of a vessel tree with flowing MBs was set up (resolution 148 μm). The trajectories and flow velocity in the vessels with a lateral distance of 100 μm were reconstructed with super-resolution. In a phantom experiment, a suspension of MBs was pumped through a tube (diameter 0.4 mm) at speeds of 2.2, 4.2, 6.3, and 10.5 mm/s and was imaged with a Vevo2100 system (Visualsonics). Estimated mean speeds of the MBs were 2.1, 4.7, 7, and 10.5 mm/s. To demonstrate the applicability for in vivo measurements, a tumor xenograft-bearing mouse was imaged by this approach. The tumor vasculature was visualized with higher resolution than in a maximum intensity persistence image and the velocity values were in the expected range 0-1 mm/s.

High resolution depth-resolved imaging from multi-focal images for medical ultrasound

An ultrasound imaging technique providing sub-diffraction limit axial resolution for point sources is proposed. It is based on simultaneously acquired multi-focal images of the same object, and on the image metric of sharpness. The sharpness is extracted by image data and presents higher values for in-focus images. The technique is derived from biological microscopy and is validated here with simulated ultrasound data. A linear array probe is used to scan a point scatterer phantom that moves in depth with a controlled step. From the beamformed responses of each scatterer position the image sharpness is assessed. Values from all positions plotted together form a curve that peaks at the receive focus, which is set during the beamforming. Selection of three different receive foci for each acquired dataset will result in the generation of three overlapping sharpness curves. A set of three calibration curves combined with the use of a maximum-likelihood algorithm is then able to estimate, with high precision, the depth location of any emitter fron each single image. Estimated values are compared with the ground truth demonstrating that an accuracy of 28.6 μm (0.13λ) is achieved for a 4 mm depth range.

The polydisperse acoustic signature of rigid microbubbles

Microbubbles are used in medical ultrasound imaging as contrast agents to image the vascular bed under the mode of Ultrasound Contrast Imaging (UCI). The microbubble shell determines the acoustic response and hence the signal that is utilized to form the images in UCI. Single microbubble signals from BiSphere™ (POINT Biomedical, San Carlos, CA, USA) microbubbles were captured using a clinical ultrasound system. Three main typical responses of microbubbles were identified, a) full duration echo, b) echo with duration shorter than the incident pulse and c) echo that in part resembles that in (b) and in addition prior to that another short duration initial lower amplitude signal. These data corroborate that the shell structural and nanomechanical property provide the different responses at different microbubble sizes. These different signals present an opportunity for tracking the movement of well differentiated single microbubbles particularly with novel super-resolution imaging methods that require sparse microbubble populations.