Cerebral blood flow estimation from perfusion-weighted MRI using FT-based MMSE filtering method

INTRODUCTION: Perfusion-weighted MRI can be used for estimating blood flow parameters using bolus tracking technique based on dynamic susceptibility contrast MRI. In order to extract flow parameters, several deconvolution techniques have been proposed, of which the singular value decomposition (SVD) and Fourier transform (FT)-based techniques are more popular and widely used. In this work, an FT-based method has been proposed that involves derivation of an optimal shaped filter (defined as a filter function) estimated using minimum mean-squared error (MMSE) technique in the frequency domain. The proposed technique has been compared with the well-established SVD technique using simulation experiments. SIMULATION METHODS: Simulation was performed in multiple steps. An arterial input function (AIF) was first defined based on a certain blood flow value. The T2* signal change was then derived from this AIF, and noise was added to the signal. Then, a unique and optimal shaped filter function Phi(f) was derived in order to obtain the best estimate of scaled residue function. One way is by minimizing the mean-squared error between the noiseless and noisy scaled residue function, i.e., using an MMSE method. The effect of low and moderate noise and distorted AIF on cerebral blood flow (CBF) estimates was obtained by using FT-based MMSE method. Results were compared with the SVD technique. In this work, SVD technique was assumed to be the standard reference deconvolution technique. RESULTS AND DISCUSSION: For low-noise condition, the FT-based technique was more stable than the SVD technique, while for moderate noise, both techniques consistently underestimated CBF. SVD technique was found to be more stable in presence of AIF distortions. However, SVD technique was found to be unstable due to AIF delay compared to the FT-based MMSE method. The shaped filter function was found to be sensitive to effect of AIF distortions.     

Effects of echo time variation on perfusion assessment using dynamic susceptibility contrast MR imaging at 3 tesla

The implications of changing the echo time of a gradient-echo echo planar imaging sequence applied to dynamic susceptibility contrast magnetic resonance imaging (DSC-MRI) for perfusion imaging at 3T were investigated. Four echo times in the range of 21 to 45 ms were examined in a total of 17 patients who received a dose of 0.1 mmol/kg bodyweight Gadobutrol (Gadovist, 1.0 mmol/ml). As the primary optimization parameter, the concentration-to-noise ratio (SNRc) was selected as it takes effects of variations in baseline as well as in signal drop into account. In an analysis of gray matter, white matter and arterial regions of interest, SNRc showed the highest values for the shortest applied echo time in all cases. Maps of regional cerebral blood volume (rCBV) and blood flow (rCBF) were calculated using deconvolution based on singular value decomposition. The quality of rCBF and rCBV images was judged to be good or excellent in all cases, independent of the echo time. Calculated gray matter/white matter ratios of rCBF and rCBV displayed no significant dependence on the applied echo time. Considering the better SNRc and arterial signal saturation aspects, we found that the shortest investigated echo time was the superior one. We thus suggest that short echo times should be applied, taking technical limitations and clinical demands into consideration.     

Application of FT-based MMSE deconvolution method for cerebral blood flow measurement in patients with leukoaraiosis

Introduction

The bolus-tracking (BT) technique is the most popular perfusion-weighted (PW) dynamic susceptibility contrast MRI method used for estimating cerebral blood flow (CBF), cerebral blood volume and mean transit time. The BT technique uses a convolution model that establishes the input–output relationship between blood flow and the vascular tracer concentration. Singular value decomposition (SVD)- and Fourier transform (FT)-based deconvolution methods are popular and widely used for estimating PW MRI parameters. However, from the published literature, it appears that SVD is more widely accepted than other methods. In a previous article, an FT-based minimum mean-squared error (MMSE) technique was proposed and simulation experiments were performed to compare it with the well-established circular SVD (oSVD) method. In this study, the FT-based MMSE method has been used to estimate relative CBF (rCBF) in 13 patients with white matter lesions (WMLs) (leukoaraiosis), and results are compared with the widely used oSVD method.

Materials and Methods

Thirteen patients with leukoaraiosis were imaged on a 1.5-T Siemens whole-body scanner. After acquiring the localizer and structural scans consisting of FLAIR (fluid attenuated with inversion recovery), T₁-weighted and T₂-weighted images, perfusion study was implemented as part of the MRI protocol. For each patient and method, two values were calculated: (a) rCBF for normal white matter (NWM) ROI, obtained by dividing the average CBF value in NWM ROI with average CBF in gray matter (GM) ROI, and (b) rCBF for WML ROI, obtained by dividing the average CBF value in WML ROI with average CBF in GM ROI. Results for the two deconvolution methods were computed.

Results and Discussion

A significant (P<.05) decrease in estimated rCBF was observed in the WML in all the patients using the MMSE method, while for the oSVD method, the decrease was observed in all but one patient. Initial results suggest that the MMSE method is comparable to the oSVD method for estimating rCBF in NMW while it may be better than oSVD for estimating rCBF in lesions of low flow. Studies involving a larger patient population may be required to further validate the findings of this work.     

Model-free arterial spin labelling for cerebral blood flow quantification: introduction of regional arterial input functions identified by factor analysis

PURPOSE: To identify regional arterial input functions (AIFs) using factor analysis of dynamic studies (FADS) when quantification of perfusion is performed using model-free arterial spin labelling. MATERIAL AND METHODS: Five healthy volunteers and one patient were examined on a 3-T Philips unit using quantitative STAR labelling of arterial regions (QUASAR). Two sets of images were retrieved, one where the arterial signal had been crushed and another where it was retained. FADS was applied to the arterial signal curves to acquire the AIFs. Perfusion maps were obtained using block-circulant SVD deconvolution and regional AIFs obtained by FADS. In the volunteers, the ASL experiment was repeated within 24 h. The patient was also examined using dynamic susceptibility contrast MRI. RESULTS: In the healthy volunteers, CBF was 64+/-10 ml/[min 100 g] (mean+/-S.D.) in GM and 24+/-4 ml/[min 100 g] in WM, while the mean aBV was 0.94% in GM and 0.25% in WM. DISCUSSION: Good CBF image quality and reasonable quantitative CBF values were obtained using the combined QUASAR/FADS technique. We conclude that FADS may be a useful supplement in the evaluation of ASL data using QUASAR.     

Study of cerebrovascular reserve capacity by magnetic resonance perfusion weighted imaging and photoacoustic imaging

The aim of this work was to assess the feasibility of photoacoustic imaging (PAI) and MR imaging for evaluating the cerebrovascular reserve capacity (CVRC) in animal models. Wistar-Kyoto (WKY) rats and spontaneous hypertensive rats (SHR) were used for MRI. BALB/c mice were used for PAI. MR perfusion weighted imaging (PWI) was performed on a 1.5-T whole-body MR system before and after oral administration of acetazolamide (ACZ). The region of interest (ROI) was chosen in the bilateral frontal lobe for measuring regional cerebral blood flow (rCBF), regional cerebral blood volume (rCBV) and mean transit time (MTT). The vessel diameters of the superficial layer of the cortex were measured by PAI in the resting and ACZ-activated mice. The results showed that there was a statistical difference between the resting and ACZ-activated animals in vessel diameter, rCBV and rCBF values. The increments in rCBV and rCBF of WKY rats between resting and ACZ test states were significantly higher than that of SHR. The pathological findings of small arterial walls and lumen of the brain were also different between WKY and SHR rats. The diameters of blood vessels in the superficial layer of the brain measured by PAI were enlarged after the ACZ tolerance test. This result was also observed in the MRI CBV map, where the signal of the vessel in the superficial layer of the cortex became redder after the ACZ stimulation, suggesting the increase of blood flow. It can be concluded that MR PWI and PAI combined with the ACZ test might be useful in evaluating the CVRC and revealing the pathologic changes in cerebral vessels.     

Sevoflurane and nitrous oxide increase regional cerebral blood flow (rCBF) and regional cerebral blood volume (rCBV) in a drug-specific manner in human volunteers

Anesthesia for diagnostic procedures, e.g., MRI measurements, has increasingly used sevoflurane and nitrous oxide in recent years. Sevoflurane and nitrous oxide are known cerebrovasodilatators, however, which potentially interferes with MRI examination of cerebral hemodynamics. To compare the effects of relevant equianesthetic concentrations (0.4 MAC) of both drugs on regional cerebral blood flow (rCBF) and regional cerebral blood volume (rCBV) we used contrast-enhanced magnetic resonance imaging (MRI) perfusion measurement, which has the advantage of providing regional anatomic resolution.

Sevoflurane increased rCBF more than did nitrous oxide in all regions except in parietal and frontal gray matter. Nitrous oxide, by contrast, increased rCBV in most of the gray matter regions more than did sevoflurane. In summary we show that, in contrast to nitrous oxide, sevoflurane supratentorially reversed the anterior-posterior gradient in rCBF and typically redistributed rCBF to infratentorial gray matter. In contrast, nitrous oxide increased rCBV more than did sevoflurane. Both inhalational anesthetics had a drug-specific influence on cerebral hemodynamics, which is of importance when interpreting MRI studies of cerebral hemodynamics in anesthetized patients.     

Quantitative assessment of intracranial tumor response to dexamethasone using diffusion, perfusion and permeability magnetic resonance imaging

There is increasing interest in obtaining quantitative imaging parameters to aid in the assessment of tumor responses to treatment. In this study, the feasibility of performing integrated diffusion, perfusion and permeability magnetic resonance imaging (MRI) for characterizing responses to dexamethasone in intracranial tumors was assessed. Eight patients with glioblastoma, five with meningioma and three with metastatic carcinoma underwent MRI prior to and 48-72 h following dexamethasone administration. The MRI protocol enabled quantification of the volume transfer constant (K(trans)), extracellular space volume fraction (nu(e)), plasma volume fraction (nu(p)), regional cerebral blood flow (rCBF), regional cerebral blood volume (rCBV), longitudinal relaxation time (T(1)) and mean diffusivity (D(av)). All subjects successfully completed the imaging protocol for the presteroid and poststeroid scans. Significant reductions were observed after the treatment for K(trans), nu(e) and nu(p) in enhancing tumor as well as for T(1) and D(av) in the edematous brain in glioblastoma; on the other hand, for meningioma, significant differences were seen only in edematous brain T(1) and D(av). No significant difference was observed for any parameter in metastatic carcinoma, most likely due to the small sample size. In addition, no significant difference was observed for enhancing tumor rCBF and rCBV in any of the tumor types, although the general trend was for rCBV to be reduced and for rCBF to be more variable. The yielded parameters provide a wealth of physiologic information and contribute to the understanding of dexamethasone actions on different types of intracranial tumors.     

Comparison of FAIR technique with different inversion times and post contrast dynamic perfusion MRI in chronic occlusive cerebrovascular disease

The purpose of this study was to examine the signal change occurring with different inversion times (TIs) of the flow-sensitive alternating inversion recovery (FAIR) technique and to compare with the perfusion image obtained with Gd-DTPA injection. The subjects were 11 patients with unilateral occlusive cerebrovascular disease. Two FAIR images with different TIs (800 ms and 1600 ms) were measured for each patient and dynamic perfusion MRI was performed to produce four kinds of parameter maps: mean transit time (MTT), time to peak (TTP), relative cerebral blood flow (rCBF) and relative cerebral blood volume (rCBV) maps. Asymmetry ratios (ARs) between the affected and contra-lateral vascular sides were measured in both FAIR images and the four dynamic parameter maps. The AR of the MTT map of the four parameters showed the highest correlation with that of the FAIR images, especially with that of TI = 1600 ms (r = 0.829), and the AR of the rCBV map revealed the worst correlation with the FAIR images. The AR of the FAIR image with TI = 800 ms was less correlated with that of MTT than that with TI = 1600 ms. These results suggested that the signal intensity of the FAIR image was influenced by flow transition time and the change in TI could be used to select the flow with a different transition time. Our study suggested that a longer TI in the FAIR technique might be more useful than a shorter TI for evaluating chronic occlusive cerebrovascular disease in the clinical setting.     

Correcting saturation effects of the arterial input function in dynamic susceptibility contrast-enhanced MRI: a Monte Carlo simulation

To prevent systematic errors in quantitative brain perfusion studies using dynamic susceptibility contrast-enhanced magnetic resonance imaging (DSC-MRI), a reliable determination of the arterial input function (AIF) is essential. We propose a novel algorithm for correcting distortions of the AIF caused by saturation of the peak amplitude and discuss its relevance for longitudinal studies. The algorithm is based on the assumption that the AIF can be separated into a reliable part at low contrast agent concentrations and an unreliable part at high concentrations. This unreliable part is reconstructed, applying a theoretical framework based on a transport-diffusion theory and using the bolus-shape in the tissue. A validation of the correction scheme is tested by a Monte Carlo simulation. The input of the simulation was a wide range of perfusion, and the main aim was to compare this input to the determined perfusion parameters. Another input of the simulation was an AIF template derived from in vivo measurements. The distortions of this template was modeled via a Rician distribution for image intensities. As for a real DSC-MRI experiment, the simulation returned the AIF and the tracer concentration-dependent signal in the tissue. The novel correction scheme was tested by deriving perfusion parameters from the simulated data for the corrected and the uncorrected case. For this analysis, a common truncated singular value decomposition approach was applied. We find that the saturation effect caused by Rician-distributed noise leads to an overestimation of regional cerebral blood flow and regional cerebral blood volume, as compared to the input parameter. The aberration can be amplified by a decreasing signal-to-noise ratio (SNR) or an increasing tracer concentration. We also find that the overestimation can be successfully eliminated by the proposed saturation-correction scheme. In summary, the correction scheme will allow DSC-MRI to be expanded towards higher tracer concentrations and lower SNR and will help to increase the measurement to measurement reproducibility for longitudinal studies.     

Analysis of input functions from different arterial branches with gamma variate functions and cluster analysis for quantitative blood volume measurements

Regional cerebral blood volume (rCBV) provides valuable information about the nature and progress of diseases of the central nervous system. While relative rCBV maps can be derived directly from dynamic susceptibility contrast data, the arterial input function (AIF) has to be measured for absolute rCBV quantification. For determination of the AIF pixels located completely within a feeding artery must be selected. However, by using a region-of-interest (ROI) based selection some confounding effects can occur, especially if single shot echo planar imaging (EPI) with low spatial resolution is used. In this study we analyzed the influence of partial volume effects and spatial misregistration due to frequency shifts induced by paramagnetic contrast agents. We analyzed AIFs from the internal carotid artery (ICA), the vertebral artery (VA) and the middle cerebral artery (MCA) using gamma variate function based parameterization. The concentration time curves (CTC) of several pixels which were selected on the basis of strong signal drop appeared distorted during the bolus passage. Moreover, the amplitudes of input functions derived from the MCA were smaller by a factor of three as compared to those of the ICA and VA. Simulations revealed that these effects can be attributed to a spatial shift of the vessel along phase-encoding direction during the passage of the bolus. We therefore developed a procedure for a pixel selection based on cluster analysis which classifies pixels according to the parameters of the fitted gamma variate functions. This approach accounted for misregistration of the vessel and yielded very consistent results for a group of normal subjects.     

Age dependency of the regional cerebral blood volume (rCBV) measured with dynamic susceptibility contrast MR imaging (DSC)

The changes of the regional cerebral blood volume (rCBV) with age were studied using dynamic susceptibility contrast MRI (DSC). We examined an unselected, random sample of 71 consecutive patients referred for work-up of suspected intracranial tumors (35 normal examinations, 36 tumors) with a standard 1.5 T clinical MR system. Determination of the rCBV was performed with a T₂¹-weighted simultaneous dual (SD) FLASH sequence (TR/TE1/TE2/α = 32/25/16/10°, 55 images) after bolus injection of Gd-DTPA. Absolute quantification of the rCBV was achieved by normalizing the measured tissue concentration-time curves with the integrated arterial input function (AIF), which was simultaneously measured in the brain feeding arteries. The rCBV (mean ± SD) was 8.4 ± 2.9 ml/100 g and 4.2 ± 1.7 ml/100 g in gray and white matter, respectively, with a decline of about 3% and 6% per decade for white and gray matter, respectively. We conclude that DSC using a SD FLASH sequence allows the simultaneous measurement of the AIF and the tissue concentration-time curve and thus an absolute quantification of the rCBV, which is the basis for interperson comparisons and follow-up studies.     

Hemodynamic investigation of peritumoral impaired blood oxygenation-level dependent cerebrovascular reactivity in patients with diffuse glioma

IntroductionThe presence of peritumorally impaired blood oxygenation-level dependent cerebrovascular reactivity (BOLD-CVR) has been unequivocally demonstrated in patients with diffuse glioma, and may have value to better identify tumor infiltration zone. Since BOLD-CVR does not measure hemodynamic changes directly, we performed additional MR perfusion studies to better characterize the peritumoral hemodynamic environment.MethodsSeventeen patients with WHO grade III and IV diffuse glioma underwent high resolution advanced hemodynamic MR imaging including BOLD-CVR and MR perfusion. The obtained multiparametric hemodynamic factors (i.e., regional cerebral blood flow (rCBF), relative cerebral blood volume (rCBV), mean transit time (MTT), time-to-peak (TTP) and BOLD-CVR, were analyzed within 10 concentric expanding 3 mm volumes of interest (VOIs) up to 30 mm from the tumor tissue mask.ResultsBOLD-CVR impairment was found within the tumor tissue mask and the peritumoral VOIs up to 21 mm as compared to the contralateral flipped CVR analysis (p<0.05). In the affected hemisphere, we observed positive spatial correlations including all VOIs between BOLD-CVR and rCBV values (r=0.27; p<0.001), rCBF (r=0.42; p<0.001) and a negative correlation between BOLD-CVR and TTP (r=-0.47; p<0.001).ConclusionsPeritumorally impaired BOLD-CVR is associated with concomitant hemodynamic alterations with severity correlating to tumor volume. The distribution of these multiparametric hemodynamic MRI patterns may be considered for future studies characterizing the hemodynamic peritumoral environment, thereby better identifying the extent of tumor infiltration.     

The influence of hyperoxia on regional cerebral blood flow (rCBF), regional cerebral blood volume (rCBV) and cerebral blood flow velocity in the middle cerebral artery (CBFVMCA) in human volunteers

Conflicting results reported on the effects of hyperoxia on cerebral hemodynamics have been attributed mainly to methodical and species differences. In the present study contrast-enhanced magnetic resonance imaging (MRI) perfusion measurement was used to analyze the influence of hyperoxia (fraction of inspired oxygen (FiO2) = 1.0) on regional cerebral blood flow (rCBF) and regional cerebral blood volume (rCBV) in awake, normoventilating volunteers (n = 19). Furthermore, the experiment was repeated in 20 volunteers for transcranial Doppler sonography (TCD) measurement of cerebral blood flow velocity in the middle cerebral artery (CBFV(MCA)). When compared to normoxia (FiO2 = 0.21), hyperoxia heterogeneously influenced rCBV (4.95 +/- 0.02 to 12.87 +/- 0.08 mL/100g (FiO2 = 0.21) vs. 4.50 +/- 0.02 to 13.09 +/- 0.09 mL/100g (FiO2 = 1.0). In contrast, hyperoxia diminished rCBF in all regions (68.08 +/- 0.38 to 199.58 +/- 1.58 mL/100g/min (FiO2 = 0.21) vs. 58.63 +/- 0.32 to 175.16 +/- 1.51 mL/100g/min (FiO2 = 1.0)) except in parietal and left frontal gray matter. CBFV(MCA) remained unchanged regardless of the inspired oxygen fraction (62 +/- 9 cm/s (FiO2 = 0.21) vs. 64 +/- 8 cm/s (FiO2 = 1.0)). Finding CBFV(MCA) unchanged during hyperoxia is consistent with the present study's unchanged rCBF in parietal and left frontal gray matter. In these fronto-parietal regions predominantly fed by the middle cerebral artery, the vasoconstrictor effect of oxygen was probably counteracted by increased perfusion of foci of neuronal activity controlling general behavior and arousal.     

Improved T(1)-weighted dynamic contrast-enhanced MRI to probe microvascularity and heterogeneity of human glioma

Dynamic contrast-enhanced (DCE) T(1)-weighted magnetic resonance imaging (MRI) is a powerful tool capable of providing quantitative assessment of contrast uptake and characterization of microvascular structure in human gliomas. The kinetics of the bolus injection doped with increasing concentrations of gadopentate dimeglumine (Gd-DTPA) depends on tissue as well as pulse sequence parameters. A simple method is described that overcomes the limitation of relative signal increase measurement and may lead to improved accuracy in quantification of perfusion indices of glioma. Based on an analysis of the contrast behavior of spoiled gradient-recalled echo sequence; a parameter K with arbitrary unit 5.0 is introduced, which provides a better approximation to the differential T(1) relaxation rate. DCE-MRI measurements of relative cerebral blood volume (rCBV) and cerebral blood flow (rCBF) were calculated in 25 patients with brain tumors (15=high-grade glioma, 10=low-grade glioma). The mean rCBV was 6.46 +/- 2.45 in high-grade glioma and 2.89 +/- 1.47 in the low-grade glioma. The rCBF was 3.94 +/- 1.47 in high-grade glioma while 2.25 +/- 0.87 in low-grade glioma. A significant difference in rCBF and rCBV was found between high- and low-grade gliomas. This simple and robust technique reveals the complexity of tumor vasculature and heterogeneity that may aid in therapeutic management especially in nonenhancing high-grade gliomas. We conclude that the precontrast medium steady-state residue parameter K may be useful in improved quantification of perfusion indices in human glioma using T(1)-weighted DCE-MRI.     

Biexponential analysis of diffusion-related signal decay in normal human cortical and deep gray matter

Diffusion imaging with high-b factors, high spatial resolution and cerebrospinal fluid signal suppression was performed in order to characterize the biexponential nature of the diffusion-related signal decay with b-factor in normal cortical gray and deep gray matter (GM). Integration of inversion pulses with a line scan diffusion imaging sequence resulted in 91% cerebrospinal fluid signal suppression, permitting accurate measurement of the fast diffusion coefficient in cortical GM (1.142+/-0.106 microm2/ms) and revealing a marked similarity with that found in frontal white matter (WM) (1.155+/-0.046 microm2/ms). The reversal of contrast between GM and WM at low vs high b-factors is shown to be due to a significantly faster slow diffusion coefficient in cortical GM (0.338+/-0.027 microm2/ms) than in frontal WM (0.125+/-0.014 microm2/ms). The same characteristic diffusion differences between GM and WM are observed in other brain tissue structures. The relative component size showed nonsignificant differences among all tissues investigated. Cellular architecture in GM and WM are fundamentally different and may explain the two- to threefold higher slow diffusion coefficient in GM.     

Visualisation of changes in regional cerebral blood flow (rCBF) produced by ketamine using long TE gradient-echo sequences: Preliminary results

Autoradiographic studies have shown that low dose ketamine produces increases in regional glucose utilisation and blood flow in the hippocampus, cerebral cortex, and olfactory lobe in the rat brain, probably due to antagonism at the NMDA receptor. Functional MRI using deoxyhaemoglobin contrast can be used to study changes in regional cerebral blood flow (rCBF). Long TE gradient-echo sequences were used to study changes in rCBF produced by low dose ketamine in rats anaesthetised with nitrous oxide, supplemented with either halothane (HAL) or fentanyl/fluanisone/midazolam (FFM) combination. Images from rats in the FFM group showed a 10–14% increase in signal intensity in the hippocampus, cerebral cortex, and olfactory lobe following either a single bolus or a low dose infusion of ketamine (p < .05). These changes were significantly reduced in the HAL group (p < .005). Halothane is known to attenuate the changes in regional glucose utilisation produced by the noncompetitive NMDA antagonist dizocilpine (MK-801), and its effects on ketamine-induced changes in rCBF seen in this study may be due to a similar effect. The potential use of functional MRI in studying the effect of pharmacological interventions on rCBF is discussed.     

Graded image segmentation of brain tissue in the presence of inhomogeneous radio frequency fields

Image segmentation is used increasingly to interpret MR spectroscopic data of the brain, using image contrast to identify gray matter (GM), white matter (WM), and cerebral spinal fluid (CSF). T(1)- or T(2)-weighted images are typically used, but poor shimming, susceptibility effects, and small variations in B(1) and receptivity cause difficulties in tissue identification. Quantitative imaging of T(1) can reduce many of these difficulties but is still subject to complications when B(1) has large variations like those observed with the surface coils often used for spectroscopy. In this study, B(1) imaging was implemented to support quantitative imaging of T(1) with either a surface coil or a volume coil. The T(1) observed by this method is a continuous function across mixtures of WM/GM and GM/CSF, and this function was measured and used to convert the images of T(1) to maps of percent GM, WM, and CSF.     

基于峰值信噪比和小波方向特性的图像奇异值去噪技术

提出一种利用小波变换子图像不同的方向特性和峰值信噪比进行奇异值分解的图像去噪算法。由于图像经过小波变换后,低频子图像集中了原图像的大部分能量噪声,故仅作简单维纳滤波;而噪声则主要集中在小波域中的三个不同方向的高频子图中,且系数较小,因此可以利用奇异值分解进行去噪处理,即用较大的奇异值和对应的特征向量重构出去噪图像,然而由于奇异值分解固有的行列方向性,对于高频对角线子图重构出的图像去噪效果不理想,故采取旋转至行列方向后再进行常用的奇异值滤波;最后将去噪后的低频和高频子图进行小波反变换重构出最终的去噪图像,其中重构所需的奇异值个数由图像的峰值信噪比确定。 实验结果表明,该方法在有效去噪的同时较好的保留了原有的高频细节信息。     

Comparison of MR perfusion imaging and microsphere measurements of regional cerebral blood flow in a rat model of middle cerebral artery occlusion

The purpose of this investigation was to correlate magnetic resonance (MR) perfusion measurements with absolute regional cerebral blood flow (rCBF) in a rat model of focal ischemia. The MR perfusion measurements were made using dynamic first-pass bolus tracking of a susceptibility contrast agent, whereas rCBF was measured using radioactive microspheres. Two simple MR perfusion parameters, the maximum change in ( ) and time delay to ( ), were derived from the signal intensity versus time curves on a pixel-to-pixel basis, without applying curve-fitting procedures or tracer kinetic theory. In each hemisphere, and were compared with the rCBF measurements in four selected regions of interest. Sixteen MR bolus tracking series were performed in 12 rats with occlusion of the middle cerebral artery. In all of the individual series there was a significant correlation (.0001 ≤ p ≤ .02) between and the microsphere rCBF measurements, with correlation coefficients ranging from .784 to .983. Pooling the data resulted in a correlation coefficient of .809 (p = .0001). There was a nonlinear correlation between the and rCBF. For both parameters there was considerable variation between different measurements regarding both the slope of the regression line and its intercept with the y-axis. Our results justify the use of as a relative measure of perfusion during acute cerebral ischemia. Because of the interindividual variation, calibration of MR perfusion measurements for the estimation of absolute flow values must be considered unreliable. The may have physiological relevance as a marker of collateral flow.     

Optimization of 3D MP-RAGE for neonatal brain imaging at 3.0 T

Three-dimensional (3D) magnetic resonance imaging (MRI) has shown great potential for studying the impact of prematurity and pathology on brain development. We have investigated the potential of optimized T1-weighted 3D magnetization-prepared rapid gradient-echo imaging (MP-RAGE) for obtaining contrast between white matter (WM) and gray matter (GM) in neonates at 3 T. Using numerical simulations, we predicted that the inversion time (TI) for obtaining strongest contrast at 3 T is approximately 2 s for neonates, whereas for adults, this value is approximately 1.3 s. The optimal neonatal TI value was found to be insensitive to reasonable variations of the assumed T1 relaxation times. The maximum theoretical contrast for neonates was found to be approximately one third of that for adults. Using the optimized TI values, MP-RAGE images were obtained from seven neonates and seven adults at 3 T, and the contrast-to-noise ratio (CNR) was measured for WM versus five GM regions. Compared to adults, neonates exhibited lower CNR between cortical GM and WM and showed a different pattern of regional variation in CNR. These results emphasize the importance of sequence optimization specifically for neonates and demonstrate the challenge in obtaining strong contrast in neonatal brain with T1-weighted 3D imaging.