Head motion detection using FID navigators.

This work explores a concept for motion detection in brain MR examinations using high channel-count RF coil arrays. It applies ultrashort (<100 μsec) free induction decay signals, making use of the knowledge that motion induces variations in these signals when compared to a reference free induction decay signal. As a proof-of-concept, the method was implemented in a standard structural MRI sequence. The stability of the free induction decay-signal was verified in phantom experiments. Human experiments demonstrated that the observed variations in the navigator data provide a sensitive measure for detection of relevant and common subject motion patterns. The proposed methodology provides a means to monitor subject motion throughout a MRI scan while causing little or no impact on the sequence timing and image contrast. It could hence complement available motion detection and correction methods, thus further reducing motion sensitivity in MR applications.

Calibration of surface contamination monitors for the detection of iodine incorporation in the thyroid gland.

In Switzerland, individuals exposed to the risk of activity intake are required to perform regular monitoring. Monitoring consists in a screening measurement and is meant to be performed using commonly available laboratory instruments. More particularly, iodine intake is measured using a surface contamination monitor. The goal of the present paper is to report the calibration method developed for thyroid screening instruments. It consists of measuring the instrument response to a known activity located in the thyroid gland of a standard neck phantom. One issue of this procedure remains that the iodine radioisotopes have a short half-life. Therefore, the adequacy and limitations to simulate the short-lived radionuclides with so-called mock radionuclides of longer half-life were also evaluated. In light of the results, it has been decided to use only the appropriate iodine sources to perform the calibration.

Volumetric quantification of atherosclerotic plaque in CT considering partial volume effect.

Coronary artery calcification (CAC) is quantified based on a computed tomography (CT) scan image. A calcified region is identified. Modified expectation maximization (MEM) of a statistical model for the calcified and background material is used to estimate the partial calcium content of the voxels. The algorithm limits the region over which MEM is performed. By using MEM, the statistical properties of the model are iteratively updated based on the calculated resultant calcium distribution from the previous iteration. The estimated statistical properties are used to generate a map of the partial calcium content in the calcified region. The volume of calcium in the calcified region is determined based on the map. The experimental results on a cardiac phantom, scanned 90 times using 15 different protocols, demonstrate that the proposed method is less sensitive to partial volume effect and noise, with average error of 9.5% (standard deviation (SD) of 5-7mm(3)) compared with 67% (SD of 3-20mm(3)) for conventional techniques. The high reproducibility of the proposed method for 35 patients, scanned twice using the same protocol at a minimum interval of 10 min, shows that the method provides 2-3 times lower interscan variation than conventional techniques.

Copper filtration in pediatric digital X-ray imaging: its impact on image quality and dose.

The effect of copper (Cu) filtration on image quality and dose in different digital X-ray systems was investigated. Two computed radiography systems and one digital radiography detector were used. Three different polymethylmethacrylate blocks simulated the pediatric body. The effect of Cu filters of 0.1, 0.2, and 0.3 mm thickness on the entrance surface dose (ESD) and the corresponding effective doses (EDs) were measured at tube voltages of 60, 66, and 73 kV. Image quality was evaluated in a contrast-detail phantom with an automated analyzer software. Cu filters of 0.1, 0.2, and 0.3 mm thickness decreased the ESD by 25-32%, 32-39%, and 40-44%, respectively, the ranges depending on the respective tube voltages. There was no consistent decline in image quality due to increasing Cu filtration. The estimated ED of anterior-posterior (AP) chest projections was reduced by up to 23%. No relevant reduction in the ED was noted in AP radiographs of the abdomen and pelvis or in posterior-anterior radiographs of the chest. Cu filtration reduces the ESD, but generally does not reduce the effective dose. Cu filters can help protect radiosensitive superficial organs, such as the mammary glands in AP chest projections.

Scaling behavior and equilibrium lengths of knotted polymers

We use numerical simulations to investigate how the chain length and topology of freely fluctuating knotted polymer rings affect their various spatial characteristics such as the radius of the smallest sphere enclosing momentary configurations of simulated polymer chains. We describe how the average value of a characteristic changes with the chain size and how this change depends on the topology of the modeled polymers. Although the scaling profiles of a spatial characteristic for distinct knot types do not intersect (at least, in the range of our data), the profiles for nontrivial knots intersect the corresponding profile obtained for phantom polymers, i.e., those that are free to explore all available topological states. For each knot type, this point of intersection defines its equilibrium length with respect to the spatial characteristic. At this chain length, a polymer forming a given knot type will not tend to increase or decrease. on average, the value of the spatial characteristic when the polymer is released from its topological constraint. We show interrelations between equilibrium lengths defined with respect to spatial characteristics of different character and observe that they are related to the lengths of ideal geometric configurations of the corresponding knot types.

A gamma camera count rate saturation correction method for whole-body planar imaging.

Whole-body (WB) planar imaging has long been one of the staple methods of dosimetry, and its quantification has been formalized by the MIRD Committee in pamphlet no 16. One of the issues not specifically addressed in the formalism occurs when the count rates reaching the detector are sufficiently high to result in camera count saturation. Camera dead-time effects have been extensively studied, but all of the developed correction methods assume static acquisitions. However, during WB planar (sweep) imaging, a variable amount of imaged activity exists in the detector's field of view as a function of time and therefore the camera saturation is time dependent. A new time-dependent algorithm was developed to correct for dead-time effects during WB planar acquisitions that accounts for relative motion between detector heads and imaged object. Static camera dead-time parameters were acquired by imaging decaying activity in a phantom and obtaining a saturation curve. Using these parameters, an iterative algorithm akin to Newton's method was developed, which takes into account the variable count rate seen by the detector as a function of time. The algorithm was tested on simulated data as well as on a whole-body scan of high activity Samarium-153 in an ellipsoid phantom. A complete set of parameters from unsaturated phantom data necessary for count rate to activity conversion was also obtained, including build-up and attenuation coefficients, in order to convert corrected count rate values to activity. The algorithm proved successful in accounting for motion- and time-dependent saturation effects in both the simulated and measured data and converged to any desired degree of precision. The clearance half-life calculated from the ellipsoid phantom data was calculated to be 45.1 h after dead-time correction and 51.4 h with no correction; the physical decay half-life of Samarium-153 is 46.3 h. Accurate WB planar dosimetry of high activities relies on successfully compensating for camera saturation which takes into account the variable activity in the field of view, i.e. time-dependent dead-time effects. The algorithm presented here accomplishes this task.

Optimization of spatial resolution for peripheral magnetic resonance angiography.

RATIONALE AND OBJECTIVES: To determine optimum spatial resolution when imaging peripheral arteries with magnetic resonance angiography (MRA). MATERIALS AND METHODS: Eight vessel diameters ranging from 1.0 to 8.0 mm were simulated in a vascular phantom. A total of 40 three-dimensional flash MRA sequences were acquired with incremental variations of fields of view, matrix size, and slice thickness. The accurately known eight diameters were combined pairwise to generate 22 "exact" degrees of stenosis ranging from 42% to 87%. Then, the diameters were measured in the MRA images by three independent observers and with quantitative angiography (QA) software and used to compute the degrees of stenosis corresponding to the 22 "exact" ones. The accuracy and reproducibility of vessel diameter measurements and stenosis calculations were assessed for vessel size ranging from 6 to 8 mm (iliac artery), 4 to 5 mm (femoro-popliteal arteries), and 1 to 3 mm (infrapopliteal arteries). Maximum pixel dimension and slice thickness to obtain a mean error in stenosis evaluation of less than 10% were determined by linear regression analysis. RESULTS: Mean errors on stenosis quantification were 8.8% +/- 6.3% for 6- to 8-mm vessels, 15.5% +/- 8.2% for 4- to 5-mm vessels, and 18.9% +/- 7.5% for 1- to 3-mm vessels. Mean errors on stenosis calculation were 12.3% +/- 8.2% for observers and 11.4% +/- 15.1% for QA software (P = .0342). To evaluate stenosis with a mean error of less than 10%, maximum pixel surface, the pixel size in the phase direction, and the slice thickness should be less than 1.56 mm2, 1.34 mm, 1.70 mm, respectively (voxel size 2.65 mm3) for 6- to 8-mm vessels; 1.31 mm2, 1.10 mm, 1.34 mm (voxel size 1.76 mm3), for 4- to 5-mm vessels; and 1.17 mm2, 0.90 mm, 0.9 mm (voxel size 1.05 mm3) for 1- to 3-mm vessels. CONCLUSION: Higher spatial resolution than currently used should be selected for imaging peripheral vessels.

MP2RAGE, a self bias-field corrected sequence for improved segmentation and T-1-mapping at high field

The large spatial inhomogeneity in transmit B, field (B-1(+)) observable in human MR images at hi h static magnetic fields (B-0) severely impairs image quality. To overcome this effect in brain T-1-weighted images the, MPRAGE sequence was modified to generate two different images at different inversion times MP2RAGE By combining the two images in a novel fashion, it was possible to create T-1-weigthed images where the result image was free of proton density contrast, T-2* contrast, reception bias field, and, to first order transmit field inhomogeneity. MP2RAGE sequence parameters were optimized using Bloch equations to maximize contrast-to-noise ratio per unit of time between brain tissues and minimize the effect of B-1(+) variations through space. Images of high anatomical quality and excellent brain tissue differentiation suitable for applications such as segmentation and voxel-based morphometry were obtained at 3 and 7 T. From such T-1-weighted images, acquired within 12 min, high-resolution 3D T-1 maps were routinely calculated at 7 T with sub-millimeter voxel resolution (0.65-0.85 mm isotropic). T-1 maps were validated in phantom experiments. In humans, the T, values obtained at 7 T were 1.15 +/- 0.06 s for white matter (WM) and 1.92 +/- 0.16 s for grey matter (GM), in good agreement with literature values obtained at lower spatial resolution. At 3 T, where whole-brain acquisitions with 1 mm isotropic voxels were acquired in 8 min the T-1 values obtained (0.81 +/- 0.03 S for WM and 1.35 +/- 0.05 for GM) were once again found to be in very good agreement with values in the literature. (C) 2009 Elsevier Inc. All rights reserved.

COMMIT: Convex Optimization Modeling for Microstructure Informed Tractography.

Tractography is a class of algorithms aiming at in vivo mapping the major neuronal pathways in the white matter from diffusion magnetic resonance imaging (MRI) data. These techniques offer a powerful tool to noninvasively investigate at the macroscopic scale the architecture of the neuronal connections of the brain. However, unfortunately, the reconstructions recovered with existing tractography algorithms are not really quantitative even though diffusion MRI is a quantitative modality by nature. As a matter of fact, several techniques have been proposed in recent years to estimate, at the voxel level, intrinsic microstructural features of the tissue, such as axonal density and diameter, by using multicompartment models. In this paper, we present a novel framework to reestablish the link between tractography and tissue microstructure. Starting from an input set of candidate fiber-tracts, which are estimated from the data using standard fiber-tracking techniques, we model the diffusion MRI signal in each voxel of the image as a linear combination of the restricted and hindered contributions generated in every location of the brain by these candidate tracts. Then, we seek for the global weight of each of them, i.e., the effective contribution or volume, such that they globally fit the measured signal at best. We demonstrate that these weights can be easily recovered by solving a global convex optimization problem and using efficient algorithms. The effectiveness of our approach has been evaluated both on a realistic phantom with known ground-truth and in vivo brain data. Results clearly demonstrate the benefits of the proposed formulation, opening new perspectives for a more quantitative and biologically plausible assessment of the structural connectivity of the brain.

Calculation of correction factors for ionization chamber measurements with small fields in low-density media.

The quantity of interest for high-energy photon beam therapy recommended by most dosimetric protocols is the absorbed dose to water. Thus, ionization chambers are calibrated in absorbed dose to water, which is the same quantity as what is calculated by most treatment planning systems (TPS). However, when measurements are performed in a low-density medium, the presence of the ionization chamber generates a perturbation at the level of the secondary particle range. Therefore, the measured quantity is close to the absorbed dose to a volume of water equivalent to the chamber volume. This quantity is not equivalent to the dose calculated by a TPS, which is the absorbed dose to an infinitesimally small volume of water. This phenomenon can lead to an overestimation of the absorbed dose measured with an ionization chamber of up to 40% in extreme cases. In this paper, we propose a method to calculate correction factors based on the Monte Carlo simulations. These correction factors are obtained by the ratio of the absorbed dose to water in a low-density medium □D(w,Q,V1)(low) averaged over a scoring volume V₁ for a geometry where V₁ is filled with the low-density medium and the absorbed dose to water □D(w,QV2)(low) averaged over a volume V₂ for a geometry where V₂ is filled with water. In the Monte Carlo simulations, □D(w,QV2)(low) is obtained by replacing the volume of the ionization chamber by an equivalent volume of water, according to the definition of the absorbed dose to water. The method is validated in two different configurations which allowed us to study the behavior of this correction factor as a function of depth in phantom, photon beam energy, phantom density and field size.

Free-breathing 3 T magnetic resonance t(2)-mapping of the heart.

OBJECTIVES: This study sought to establish an accurate and reproducible T(2)-mapping cardiac magnetic resonance (CMR) methodology at 3 T and to evaluate it in healthy volunteers and patients with myocardial infarct. BACKGROUND: Myocardial edema affects the T(2) relaxation time on CMR. Therefore, T(2)-mapping has been established to characterize edema at 1.5 T. A 3 T implementation designed for longitudinal studies and aimed at guiding and monitoring therapy remains to be implemented, thoroughly characterized, and evaluated in vivo. METHODS: A free-breathing navigator-gated radial CMR pulse sequence with an adiabatic T(2) preparation module and an empirical fitting equation for T(2) quantification was optimized using numerical simulations and was validated at 3 T in a phantom study. Its reproducibility for myocardial T(2) quantification was then ascertained in healthy volunteers and improved using an external reference phantom with known T(2). In a small cohort of patients with established myocardial infarction, the local T(2) value and extent of the edematous region were determined and compared with conventional T(2)-weighted CMR and x-ray coronary angiography, where available. RESULTS: The numerical simulations and phantom study demonstrated that the empirical fitting equation is significantly more accurate for T(2) quantification than that for the more conventional exponential decay. The volunteer study consistently demonstrated a reproducibility error as low as 2 ± 1% using the external reference phantom and an average myocardial T(2) of 38.5 ± 4.5 ms. Intraobserver and interobserver variability in the volunteers were -0.04 ± 0.89 ms (p = 0.86) and -0.23 ± 0.91 ms (p = 0.87), respectively. In the infarction patients, the T(2) in edema was 62.4 ± 9.2 ms and was consistent with the x-ray angiographic findings. Simultaneously, the extent of the edematous region by T(2)-mapping correlated well with that from the T(2)-weighted images (r = 0.91). CONCLUSIONS: The new, well-characterized 3 T methodology enables robust and accurate cardiac T(2)-mapping at 3 T with high spatial resolution, while the addition of a reference phantom improves reproducibility. This technique may be well suited for longitudinal studies in patients with suspected or established heart disease.

Evaluation of monoenergetic imaging to reduce metallic instrumentation artifacts in computed tomography of the cervical spine.

OBJECT Monoenergetic imaging with dual-energy CT has been proposed to reduce metallic artifacts in comparison with conventional polychromatic CT. The purpose of this study is to systematically evaluate and define the optimal dual-energy CT imaging parameters for specific cervical spinal implant alloy compositions. METHODS Spinal fixation rods of cobalt-chromium or titanium alloy inserted into the cervical spine section of an Alderson Rando anthropomorphic phantom were imaged ex vivo with fast-kilovoltage switching CT at 80 and 140 peak kV. The collimation width and field of view were varied between 20 and 40 mm and medium to large, respectively. Extrapolated monoenergetic images were generated at 70, 90, 110, and 130 kiloelectron volts (keV). The standard deviation of voxel intensities along a circular line profile around the spine was used as an index of the magnitude of metallic artifact. RESULTS The metallic artifact was more conspicuous around the fixation rods made of cobalt-chromium than those of titanium alloy. The magnitude of metallic artifact seen with titanium fixation rods was minimized at monoenergies of 90 keV and higher, using a collimation width of 20 mm and large field of view. The magnitude of metallic artifact with cobalt-chromium fixation rods was minimized at monoenergies of 110 keV and higher; collimation width or field of view had no effect. CONCLUSIONS Optimization of acquisition settings used with monoenergetic CT studies might yield reduced metallic artifacts.

Simultaneous EEG-fMRI at ultra-high field: Artifact prevention and safety assessment.

The simultaneous recording of scalp electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) can provide unique insights into the dynamics of human brain function, and the increased functional sensitivity offered by ultra-high field fMRI opens exciting perspectives for the future of this multimodal approach. However, simultaneous recordings are susceptible to various types of artifacts, many of which scale with magnetic field strength and can seriously compromise both EEG and fMRI data quality in recordings above 3T. The aim of the present study was to implement and characterize an optimized setup for simultaneous EEG-fMRI in humans at 7T. The effects of EEG cable length and geometry for signal transmission between the cap and amplifiers were assessed in a phantom model, with specific attention to noise contributions from the MR scanner coldheads. Cable shortening (down to 12cm from cap to amplifiers) and bundling effectively reduced environment noise by up to 84% in average power and 91% in inter-channel power variability. Subject safety was assessed and confirmed via numerical simulations of RF power distribution and temperature measurements on a phantom model, building on the limited existing literature at ultra-high field. MRI data degradation effects due to the EEG system were characterized via B0 and B1(+) field mapping on a human volunteer, demonstrating important, although not prohibitive, B1 disruption effects. With the optimized setup, simultaneous EEG-fMRI acquisitions were performed on 5 healthy volunteers undergoing two visual paradigms: an eyes-open/eyes-closed task, and a visual evoked potential (VEP) paradigm using reversing-checkerboard stimulation. EEG data exhibited clear occipital alpha modulation and average VEPs, respectively, with concomitant BOLD signal changes. On a single-trial level, alpha power variations could be observed with relative confidence on all trials; VEP detection was more limited, although statistically significant responses could be detected in more than 50% of trials for every subject. Overall, we conclude that the proposed setup is well suited for simultaneous EEG-fMRI at 7T.

Comparison of three commercially available radio frequency coils for human brain imaging at 3 Tesla.

OBJECTIVE: To evaluate a transverse electromagnetic (TEM), a circularly polarized (CP) (birdcage), and a 12-channel phased array head coil at the clinical field strength of B0 = 3T in terms of signal-to-noise ratio (SNR), signal homogeneity, and maps of the effective flip angle alpha. MATERIALS AND METHODS: SNR measurements were performed on low flip angle gradient echo images. In addition, flip angle maps were generated for alpha(nominal) = 30 degrees using the double angle method. These evaluation steps were performed on phantom and human brain data acquired with each coil. Moreover, the signal intensity variation was computed for phantom data using five different regions of interest. RESULTS: In terms of SNR, the TEM coil performs slightly better than the CP coil, but is second to the smaller 12-channel coil for human data. As expected, both the TEM and the CP coils show superior image intensity homogeneity than the 12-channel coil, and achieve larger mean effective flip angles than the combination of body and 12-channel coil with reduced radio frequency power deposition. CONCLUSION: At 3T the benefits of TEM coil design over conventional lumped element(s) coil design start to emerge, though the phased array coil retains an advantage with respect to SNR performance.

3-d residual eddy current field characterisation: applied to diffusion weighted magnetic resonance imaging.

Clinical use of the Stejskal-Tanner diffusion weighted images is hampered by the geometric distortions that result from the large residual 3-D eddy current field induced. In this work, we aimed to predict, using linear response theory, the residual 3-D eddy current field required for geometric distortion correction based on phantom eddy current field measurements. The predicted 3-D eddy current field induced by the diffusion-weighting gradients was able to reduce the root mean square error of the residual eddy current field to ~1 Hz. The model's performance was tested on diffusion weighted images of four normal volunteers, following distortion correction, the quality of the Stejskal-Tanner diffusion-weighted images was found to have comparable quality to image registration based corrections (FSL) at low b-values. Unlike registration techniques the correction was not hindered by low SNR at high b-values, and results in improved image quality relative to FSL. Characterization of the 3-D eddy current field with linear response theory enables the prediction of the 3-D eddy current field required to correct eddy current induced geometric distortions for a wide range of clinical and high b-value protocols.