3D reconstruction of a carotid bifurcation from 2D transversal ultrasound images

Visualizing and analyzing the morphological structure of carotid bifurcations are important for understanding the etiology of carotid atherosclerosis, which is a major cause of stroke and transient ischemic attack. For delineation of vasculatures in the carotid artery, ultrasound examinations have been widely employed because of a noninvasive procedure without ionizing radiation. However, conventional 2D ultrasound imaging has technical limitations in observing the complicated 3D shapes and asymmetric vasodilation of bifurcations. This study aims to propose image-processing techniques for better 3D reconstruction of a carotid bifurcation in a rat by using 2D cross-sectional ultrasound images. A high-resolution ultrasound imaging system with a probe centered at 40 MHz was employed to obtain 2D transversal images. The lumen boundaries in each transverse ultrasound image were detected by using three different techniques; an ellipse-fitting, a correlation mapping to visualize the decorrelation of blood flow, and the ellipse-fitting on the correlation map. When the results are compared, the third technique provides relatively good boundary extraction. The incomplete boundaries of arterial lumen caused by acoustic artifacts are somewhat resolved by adopting the correlation mapping and the distortion in the boundary detection near the bifurcation apex was largely reduced by using the ellipse-fitting technique. The 3D lumen geometry of a carotid artery was obtained by volumetric rendering of several 2D slices. For the 3D vasodilatation of the carotid bifurcation, lumen geometries at the contraction and expansion states were simultaneously depicted at various view angles. The present 3D reconstruction methods would be useful for efficient extraction and construction of the 3D lumen geometries of carotid bifurcations from 2D ultrasound images.     

Intracranial vascular feature changes in time of flight MR angiography in patients undergoing carotid revascularization surgery

PurposeTo characterize the intracranial vascular features extracted from time of flight (TOF) images and their changes from baseline to follow-up in patients undergoing carotid revascularization, using arterial spin labeling (ASL) cerebral blood flow (CBF) measurement as a reference.MethodsIn this retrospective study, brain TOF and ASL images of 99 subjects, acquired before, within 48 h, and/or 6 months after, carotid revascularization surgery were analyzed. TOF images were analyzed using a custom software (iCafe) to quantify intracranial vascular features, including total vessel length, total vessel volume, and number of branches. Mean whole-brain CBF was calculated from ASL images. ASL scans showing low ASL signal in the entire flow territory of an internal carotid artery (ICA), which may be caused by labeling failure, were excluded. Changes and correlations between time points were analyzed separately for TOF intracranial vascular features and ASL CBF.ResultsSimilar to ASL CBF, TOF vascular features (i.e. total vessel length, total vessel volume and number of branches) increased dramatically from baseline to post-surgery, then returned to a level slightly higher than the baseline in long-term follow-up (All P < 0.05). Correlation between time points was observed for all three TOF vascular features but not for ASL CBF.ConclusionIntracranial vascular features, including total vessel length, total vessel volume and number of branches, extracted from TOF images are useful in detecting brain blood flow changes induced by carotid revascularization surgery.     

Doubly tuned solenoidal resonators for small animal imaging and spectroscopy at 1.5 Tesla

The design and construction of solenoidal resonators for use with small animals in a 1.5-Tesla clinical imaging system are described. The coils have been designed to exploit the B1 distributions of two resonant modes of a four-turn solenoid whose windings are in parallel. Both singly and doubly tuned versions have been constructed. 1H images of normal and pathologic anatomy in mice and rats as well as a 31P spectrum of a Walker 256 rat sarcoma are presented. A primary advantage of this design is that the coils are easy to build and implement while providing the necessary sensitivity to allow high quality images to be obtained with no changes to the hardware or software of the clinical unit.     

小波分析在灰度遥感图像去云中的应用

针对薄云覆盖的灰度遥感图像其清晰度低、难于判读地物信息等问题, 提出了小波分析方法用于灰度遥感图像去云。以256 pixel×256 pixel的灰度遥感图像为例, 进行1～8级分解, 重构时对细节系数放大1～10倍, 先选定小波函数确定图像信息熵取得最高时的分解级数为7和相应细节系数的放大倍数, 再确定双正交2.2型小波和离散近似Mery小波, 以0.1为步长的细节放大倍数作进一步分析。结果表明, 分解的级数较低时, 图像信息熵对细节系数放大倍数的依赖性很强; 放大倍数为1.1时得到的图像其信息熵取得极大值; 采用离散近似Mery小波, 当放大倍数为4.8和5.7时还出现2个峰值。     

High signal-to-noise FLASH imaging at 8 Tesla.

A radio frequency (RF) and gradient spoiled fast low angle shot technique was used to acquire images from the human brain at 8 Tesla. The resulting FLASH images, obtained with a 17 degrees nutation, a 70 ms repetition time, and a 17 ms echo time, displayed an average signal-to-noise ratio (SNR) of 220:1 (slice thickness 2.2 mm, field-of-view 24 cm, matrix 256 x 128). These images were compared with images obtained at 1.5 Tesla using identical parameters yielding a signal-to-noise of less than 10:1. As such, the 8 Tesla images display a remarkable improvement in SNR with increasing field strength. The images also show little evidence of susceptibility distortion, chemical shift, or RF penetration limitations.     

Computerized quantification of carotid artery stenosis using MRA axial images

Currently, the North American symptomatic carotid endarterectomy trial, European carotid surgery trial, and common carotid method are used to measure the carotid stenosis for determining candidates for carotid endarterectomy using the projection angiography from different modalities such as digital subtraction angiography, rotational angiography, computed tomography angiography and magnetic resonance angiography. A new computerized carotid stenosis measuring system was developed using MR angiography axial image to overcome the drawbacks of conventional carotid stenosis measuring methods, to reduce the variability of inter-observer and intra-observer. The gray-level thresholding is one of the most popular and efficient methods for image segmentation. We segmented the carotid artery and lumen from three-dimensional time-of-flight MRA axial images using gray-level thresholding technique. Using the measured intima-media thickness value of common carotid artery for each case, we separated carotid artery wall from the segmented carotid artery region. After that, the regions of segmented carotid without artery wall were divided into region of blood flow and plaque. The calculation of carotid stenosis degree was performed as follows: carotid stenosis grading = (area measure of plaque/area measure of blood flow region and plaque) * 100%. No previous study has developed the carotid stenosis measuring method using MRA axial image. The new computerized stenosis measuring system has advantage over conventional caliper measuring methods; it will not only greatly increase the speed of stenosis measuring but also reduce the variability between readers. It should also reduce the variability between different institutions.     

Causal dynamic MRI reconstruction via nuclear norm minimization

This work addresses the problem of online reconstruction of dynamic magnetic resonance images (MRI). The proposed method reconstructs the difference between the images of previous and current time frames. This difference image is modeled as a rank deficient matrix and is solved from the partially sampled k-space data via nuclear norm minimization. Our proposed method has been compared against state-of-the-art offline and online reconstruction methods. Our method has similar reconstruction accuracy as the offline method and significantly higher accuracy compared to the online technique. It is about an order of magnitude faster than the online technique compared against. Our experimental data consisted of dynamic MRI data that were collected at 6 to 7 frames per second and having resolutions of 128×128 and 256×256 pixels per frame. Experimental evaluation indicates that our proposed method is capable of reconstructing 128×128 images at the rate of 4 frames per second and 256×256 images at the rate of 2 frames per second. The previous online method requires about 3.75 s for reconstructing each image. The improvement in reconstruction speed is clearly discernible. Moreover, our method has a reconstruction error that is about half that of the previous online method.     

Fast inversion recovery T1 contrast and chemical shift contrast in high-resolution snapshot FLASH MR images.

Fast MR imaging attracts the interest of both clinicians and physicists because new diagnostic information arises with reduced artifacts due to short investigational times. With the acceleration of the Snapshot FLASH MR sequence, the measurement of high-resolution images with 256 x 256 matrix is reported, together with contrasting prepulses that are applied to attain contrast in combination with higher in-plane resolution. Measuring times are in the range of a second. For whole-body imaging, a TR = 5.2 msec and a TE = 2.6 msec could be attained measuring omit 256 x 256 matrix images. Artifact-free images demonstrating T1 contrast and contrast from chemical shift are performed on moving organs (heart, intestine) in different experiments. These applications can easily be performed in a couple of minutes for clinical use. Especially in the lung, short TE and high resolution result in a new imaging quality of pulmonary and mediastinal vessels.     

Stable periodic vortex shedding studied using computational fluid dynamics, laser sheet flow visualization, and MR imaging

Recirculating and detached flow patterns close to the carotid bifurcation are thought to play an important role in the development of carotid stenoses by promoting atherosclerosis. The aim of this study was to investigate a flow regime with strong transient characteristics, including vortex shedding and transport to develop methodologies appropriate to the analysis of carotid stenoses. The existence of a regular periodic vortex street behind a cylindrical flow obstruction was predicted and analysed in detail by Theodore van Karman in the early 20th century. This model was chosen in our study for both ease of phantom construction and of theoretical modelling using finite element computational fluid dynamics (CFD). The results of the theoretical calculations have been compared with two methods of flow visualization-laser sheet imaging and real-time echo planar magnitude MR imaging. Flow was investigated over a range of Reynold's number from 40 through 400 through which vortex shedding is predicted. Good overall agreement was obtained between the theoretical (16 mm-CFD) and experimental (16+/-2 mm-Laser, 17+/-2 mm-MRI) estimates of the Karman Vortex street wavelength for a Reynolds number of 200.     

Spin-echo fluorine magnetic resonance imaging at 2 T: in vivo spatial distribution of halothane in the rabbit head

Spin-echo 19F magnetic resonance imaging was performed at 2.0 T to explore the in vivo spatial distribution of halothane in the rabbit head. Because the halothane concentration is low in vivo, and because the measured relaxation times of the 19F resonance peak for halothane were T1 approximately equal to 1.0 sec and T2 approximately equal to 3.5-65 msec, 1-3-h imaging times were required (TR = 1 sec, TE = 9 msec) in order to obtain adequate images with a 64 X 256 raw data matrix and a 20-mm slice thickness. With this technique, halothane was primarily detected in lipophilic regions of the rabbit head, but little or no halothane was observed in brain tissue. Because T2 was shorter in brain tissue than in surrounding fat, a shorter TE than we could obtain is needed for optimal spin-echo imaging of brain halothane.     

3D dynamical ultrasonic model of pulsating vessel walls

The aim of this work is to introduce a novel 3-D model of pulsating vessels, through which the dynamic acoustic response of arterial regions can be predicted. Blood flow is numerically simulated by considering the fluid-dynamic displacements of the scatterers (erythrocytes), while a mechanical model calculates the wall displacement due to fluid pressure. The acoustic characteristics of each region are simulated through the FIELD software. Two numerical phantoms of a carotid artery surrounded by elastic tissue have been developed to illustrate the model. One of them includes a plaque involving a 50% stenosis. B-mode and M-mode images are produced and segmented to obtain the wall displacement profile. A cylindrical holed phantom made of cryogel mimicking material has been constructed for the model validation. In pulsatile flow conditions, fluid and wall displacements have been measured by Doppler ultrasound methods and quantitatively compared to simulated M-mode images, showing a fairly good agreement.     

Colours in a complex fluid flow

The addition of colours to flow visualization allows a third dimension to be added to the intrinsically two-dimensional information rendered by an image.In the present work, the velocity domain and the wall heat transfer in the near field of a jet in cross-flow is experimentally investigated by making extensive use of coloured images. Tests are performed in a low turbulence wind tunnel at a jet Reynolds number equal to 8000 and for velocity ratios ranging from 1 to 5. Data are obtained with Particle Image Velocimetry (PIV) and with infrared thermography applied to the steady state heated-thin-foil technique.     

A comparison of three SPRITE techniques for the quantitative 3D imaging of the 23Na spin density on a 4T whole-body machine

Sodium density maps acquired with three SPRITE-based methods have been compared in terms of the resulting quantitative information as well as image quality and acquisition times. Consideration of factors relevant for the clinical implementation of SPRITE shows that the Conical-SPRITE variant is preferred because of a 20-fold reduction in acquisition time, slightly improved image quality, and no loss of quantitative information. The acquisition of a 3D data set (32x32x16; FOV=256x256x160 mm) for the quantitative determination of sodium density is demonstrated. In vivo Conical-SPRITE 23Na images of the brain of a healthy volunteer were acquired in 30 min with a resolution of 7.5x7.5x7.5 mm and a signal-to-noise ratio of 23 in cerebrospinal fluid and 17 in brain tissue.     

Quantitative 'real-time' imaging of multi-phase flow in ceramic monoliths

An extension of the RARE technique has been developed which acquires multiple images from a single radio-frequency excitation. This pulse sequence has been used to image, in real-time, gas flow through stagnant liquid within parallel-channel ceramic monoliths. From these images, gas-phase volume fractions, and distributions of gas bubble length and velocity as a function of gas flow rate (50-300 cm3 min(-1)) and channel size (300 and 400 channels per square inch, cpsi) are obtained directly. Increasing the gas flow rate increased the number of large bubbles and the average bubble velocity. A bimodal distribution in the bubble velocities was observed for flow within the larger channel size (300 cpsi) in contrast to a broad unimodal distribution characterizing two-phase flow within the smaller channel size (400 cpsi).     

Rock porosity quantification by dual-energy X-ray computed microtomography

Porous media investigation by X-ray microtomography allows obtaining valuable quantitative and qualitative information, while preserving sample integrity. Modern X-ray nanotomography or Synchrotron radiation systems may distinguish structures sized only hundreds of nanometers. However, pores sized less than a few microns (microporosity) may be undetectable due to the system’s spatial resolution and noise in microfocus sources, compromising the quality of the measurement. In this study a dual-energy methodology was developed to generate density-based images from two scans made at two different voltages (80 kV and 130 KV) with a microfocus bench-top microtomography system. The images obtained were quantized in 256 gray levels, where the lowest value (zero) corresponded to voids and the highest value (255) corresponded to the densest regions mapped. From density images and single energy images, porosity was evaluated and compared. Results indicate that density images present better results than single energy images when both are compared with porosity obtained by the helium injection method. In addition, images acquired in dual-energy show good agreement with the sample’s real density values.     

Reliability of mean transit time obtained using perfusion-weighted MR imaging; comparison with positron emission tomography

The purpose of this project was to assess the reliability of the cerebral mean transit time (MTT) obtained using perfusion-weighted MR imaging by comparing it with the MTT obtained when performing positron emission tomography (PET). Ten patients with chronic occlusive cerebrovascular disease were investigated. They had either unilateral internal carotid artery occlusion or middle cerebral artery occlusion. The regions-of-interest were placed in non-infarcted areas within the territory of the middle cerebral artery on the affected side. Control regions-of-interest were placed in mirrored regions of the contralateral side. Linear regression analyses were performed using the parameters of the MTT obtained with perfusion-weighted MR imaging and the MTT, cerebral blood flow, vascular reactivity, and oxygen extraction fraction obtained with PET. The respective MTTs of the affected and non-affected sides obtained with perfusion-weighted MR imaging versus those with PET were 7.3 +/- 2.2 s and 6.0 +/- 1.2 s versus 8.2 +/- 3.0 s and 6.4 +/- 1.7 s. The MTT obtained using perfusion-weighted MR imaging and PET demonstrated statistically significant correlation (r = 0.87, p < 0.0001). The MTT obtained with perfusion-weighted MR imaging correlated statistically with cerebral blood flow (r = -0.74, p < 0.001), vascular reactivity (r = -0.73, p < 0.001) and oxygen extraction fraction (r = 0.61, p < 0.01). Similarly, the MTT obtained using PET statistically correlated with cerebral blood flow (r = -0.78, p < 0.0001), vascular reactivity (r = -0.51, p < 0.05) and oxygen extraction fraction (r = 0.68, p < 0.01). The reliability of the MTT obtained using perfusion-weighted MR imaging appears to be approximately equal to that obtained with positron emission tomography.     

Comparison of carotid vessel wall area measurements using three different contrast-weighted black blood MR imaging techniques

Measuring carotid artery plaque burden from MRI is a reliable method for monitoring regression and progression of atherosclerosis. However, to measure all available images would be very time consuming, and in practice the image quality (IQ) of these images may be inconsistent, which can directly impact the quality of measurement. It is hypothesized that if IQ is comparable among different contrast weighted images, then carotid artery area measurements obtained from different contrast images of the same location will produce identical results. To test this, T1, proton density and T2 weighted images were acquired from ten patients (51 +/- 7 years old). Carotid lumen and vessel wall area was measured using a custom designed software program. The results showed strong agreement evidenced with only small differences on both lumen (mean: 40.5 mm(2)) and wall (mean: 52.6 mm(2)) area measurement among different weighted images. The maximum absolute mean differences are less than 2.7 mm(2) and 4.4 mm(2), and 90(th) percentile of the absolute differences are 5.6 mm(2) and 8.2 mm(2) respectively. In conclusion, different contrast weighted images with high and comparable IQ will yield similar results in lumen and vessel wall area measurement. At each matched location, it is recommended that the image with the highest IQ be used for area measurement.     

Improvement of ultrasound speckle image velocimetry using image enhancement techniques

Ultrasound-based techniques have been developed and widely used in noninvasive measurement of blood velocity. Speckle image velocimetry (SIV), which applies a cross-correlation algorithm to consecutive B-mode images of blood flow has often been employed owing to its better spatial resolution compared with conventional Doppler-based measurement techniques. The SIV technique utilizes speckles backscattered from red blood cell (RBC) aggregates as flow tracers. Hence, the intensity and size of such speckles are highly dependent on hemodynamic conditions. The grayscale intensity of speckle images varies along the radial direction of blood vessels because of the shear rate dependence of RBC aggregation. This inhomogeneous distribution of echo speckles decreases the signal-to-noise ratio (SNR) of a cross-correlation analysis and produces spurious results. In the present study, image-enhancement techniques such as contrast-limited adaptive histogram equalization (CLAHE), min/max technique, and subtraction of background image (SB) method were applied to speckle images to achieve a more accurate SIV measurement. A mechanical sector ultrasound scanner was used to obtain ultrasound speckle images from rat blood under steady and pulsatile flows. The effects of the image-enhancement techniques on SIV analysis were evaluated by comparing image intensities, velocities, and cross-correlation maps. The velocity profiles and wall shear rate (WSR) obtained from RBC suspension images were compared with the analytical solution for validation. In addition, the image-enhancement techniques were applied to in vivo measurement of blood flow in human vein. The experimental results of both in vitro and in vivo SIV measurements show that the intensity gradient in heterogeneous speckles has substantial influence on the cross-correlation analysis. The image-enhancement techniques used in this study can minimize errors encountered in ultrasound SIV measurement in which RBCs are used as flow tracers instead of exogenous contrast agents.     

Combined two-dimensional velocity and temperature measurements using a high-speed camera and luminescent particles

This paper proposes a combined method for two-dimensional temperature and velocity measurements using temperature sensitive particles (TSParticles), a pulsed ultraviolet (UV) laser and a single high-speed camera. TSParticles were synthesized using ion-exchange particles and Eu(TTA) luminescent dye. The size and material of the particles for synthesizing TSParticles are selectable. TSParticles respond to temperature changes in a flow and can also serve as tracers for the velocity field. TSParticles were seeded into a heated water flow in a complex-shaped channel constructed of MEXFLON resin, which has a refractive index exactly equal to that of water. Particle images of flow beyond the structure can be recorded without any distortion. The TSParticles were excited by the UV pulsed laser and the luminescence from the TSParticles were recorded at 40,000 frames per second as sequential images for a lifetime-based temperature analysis. Another advantage of our approach is that high time-resolved PIV can be carried out without a high-frequency laser. The recorded images were also used for the particle image velocimetry (PIV) calculation.     

环形压气机叶棚内部分离流结构分析 第二部分:流道内分离流

1前言认识压气机叶栅内部分离流动的结构和特性,对于揭示压气机内复杂的流动机制、改善叶栅气动性能、减少流动损失和提高压气机喘振裕度具有重要意义山。本文采用激光片光层加示踪粒子流场显示技术,观察和摄取了环形压气机叶栅在约十14”大攻角下的内部流动图像;采用七孔气