Multilevel oblique corpectomy for cervical spondylotic myelopathy preserves segmental motion.

PURPOSE: To document the neurological outcome, spinal alignment and segmental range of movement after oblique cervical corpectomy (OCC) for cervical compressive myelopathy. METHODS: This retrospective study included 109 patients--93 with cervical spondylotic myelopathy and 16 with ossified posterior longitudinal ligament in whom spinal curvature and range of segmental movements were assessed on neutral and dynamic cervical radiographs. Neurological function was measured by Nurick's grade and modified Japanese Orthopedic Association (JOA) scores. Eighty-eight patients (81%) underwent either a single- or two-level corpectomy; the remaining (19%) undergoing three- or four-level corpectomies. The average duration of follow-up was 30.52 months. RESULTS: The Nurick's grade and the JOA scores showed statistically significant improvements after surgery (p < 0.001). The mean postoperative segmental angle in the neutral position straightened by 4.7 ± 6.5°. The residual segmental range of movement for a single-level corpectomy was 16.7° (59.7% of the preoperative value), for two-level corpectomy it was 20.0° (67.2%) and for three-level corpectomies it was 22.9° (74.3%). 63% of patients with lordotic spines continued to have lordosis postoperatively while only one became kyphotic without clinical worsening. Four patients with preoperative kyphotic spines showed no change in spine curvature. None developed spinal instability. CONCLUSIONS: The OCC preserves segmental motion in the short-term, however, the tendency towards straightening of the spine, albeit without clinical worsening, warrants serial follow-up imaging to determine whether this motion preservation is long lasting.

Dependence of brain intravoxel incoherent motion perfusion parameters on the cardiac cycle.

Measurement of microvascular perfusion with Intravoxel Incoherent Motion (IVIM) MRI is gaining interest. Yet, the physiological influences on the IVIM perfusion parameters ("pseudo-diffusion" coefficient D*, perfusion fraction f, and flow related parameter fD*) remain insufficiently characterized. In this article, we hypothesize that D* and fD*, which depend on blood speed, should vary during the cardiac cycle. We extended the IVIM model to include time dependence of D* = D*(t), and demonstrate in the healthy human brain that both parameters D* and fD* are significantly larger during systole than diastole, while the diffusion coefficient D and f do not vary significantly. The results non-invasively demonstrate the pulsatility of the brain's microvasculature.

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.

Respiratory motion artifact suppression in diffusion-weighted MR imaging of the spine.

Diffusion-weighted spin-echo imaging of the spine has been successfully implemented for differentiation of benign fracture edema and tumor infiltration of the vertebral body. Nevertheless, this technique still suffers from insufficient image quality in numerous patients due to motion artifacts. The aim of this study was to investigate the impact of variable respiratory motion artifact suppression techniques on image quality in diffusion-weighted spin-echo imaging of the spine. In addition to phase-encoding reordering, a newly implemented right hemi-diaphragmaitc navigator for respiratory gating was used. Subjective and objective image quality parameters were compared. Respiratory motion artifact suppression has a major impact on image quality in diffusion-weighted imaging of the spine. Phase-encoding reordering does not enhance image quality while right hemi-diaphragmatic respiratory navigator gating significantly improves image quality at the cost of data acquisition time. Navigator gating should be used if standard spin-echo diffusion-weighted imaging demonstrates insufficient image quality.

Assessment of left coronary artery beat-to-beat repositioning in order to propose an optimal acquisition window for coronary MRA

Introduction: Coronary magnetic resonance angiography (MRA) is a medical imaging technique that involves collecting data from consecutive heartbeats, always at the same time in the cardiac cycle, in order to minimize heart motion artifacts. This technique relies on the assumption that coronary arteries always follow the same trajectory from heartbeat to heartbeat. Until now, choosing the acquisition window in the cardiac cycle was based exclusively on the position of minimal coronary motion. The goal of this study was to test the hypothesis that there are time intervals during the cardiac cycle when coronary beat-to-beat repositioning is optimal. The repositioning uncertainty values in these time intervals were then compared with the intervals of low coronary motion in order to propose an optimal acquisition window for coronary MRA. Methods: Cine breath-hold x-ray angiograms with synchronous ECG were collected from 11 patients who underwent elective routine diagnostic coronarography. Twenty-three bifurcations of the left coronary artery were selected as markers to evaluate repositioning uncertainty and velocity during cardiac cycle. Each bifurcation was tracked by two observers, with the help of a user-assisted algorithm implemented in Matlab (The Mathworks, Natick, MA, USA) that compared the trajectories of the markers coming from consecutive heartbeats and computed the coronary repositioning uncertainty with steps of 50ms until 650ms after the R-wave. Repositioning uncertainty was defined as the diameter of the smallest circle encompassing the points to be compared at the same time after the R-wave. Student's t-tests with a false discovery rate (FDR, q=0.1) correction for multiple comparison were applied to see whether coronary repositioning and velocity vary statistically during cardiac cycle. Bland-Altman plots and linear regression were used to assess intra- and inter-observer agreement. Results: The analysis of left coronary artery beat-to-beat repositioning uncertainty shows a tendency to have better repositioning in mid systole (less than 0.84±0.58mm) and mid diastole (less than 0.89±0.6mm) than in the rest of the cardiac cycle (highest value at 50ms=1.35±0.64mm). According to Student's t-tests with FDR correction for multiple comparison (q=0.1), two intervals, in mid systole (150-200ms) and mid diastole (550-600ms), provide statistically better repositioning in comparison with the early systole and the early diastole. Coronary velocity analysis reveals that left coronary artery moves more slowly in end systole (14.35±11.35mm/s at 225ms) and mid diastole (11.78±11.62mm/s at 625ms) than in the rest of the cardiac cycle (highest value at 25ms: 55.96±22.34mm/s). This was confirmed by Student's t-tests with FDR correction for multiple comparison (q=0.1, FDR-corrected p-value=0.054): coronary velocity values at 225, 575 and 625ms are not much different between them but they are statistically inferior to all others. Bland-Altman plots and linear regression show that intra-observer agreement (y=0.97x+0.02 with R²=0.93 at 150ms) is better than inter-observer (y=0.8x+0.11 with R²=0.67 at 150ms). Discussion: The present study has demonstrated that there are two time intervals in the cardiac cycle, one in mid systole and one in mid diastole, where left coronary artery repositioning uncertainty reaches points of local minima. It has also been calculated that the velocity is the lowest in end systole and mid diastole. Since systole is less influenced by heart rate variability than diastole, it was finally proposed to test an acquisition window between 150 and 200ms after the R-wave.

Perfusion measurement in brain gliomas with intravoxel incoherent motion MRI.

BACKGROUND AND PURPOSE: Intravoxel incoherent motion MRI has been proposed as an alternative method to measure brain perfusion. Our aim was to evaluate the utility of intravoxel incoherent motion perfusion parameters (the perfusion fraction, the pseudodiffusion coefficient, and the flow-related parameter) to differentiate high- and low-grade brain gliomas. MATERIALS AND METHODS: The intravoxel incoherent motion perfusion parameters were assessed in 21 brain gliomas (16 high-grade, 5 low-grade). Images were acquired by using a Stejskal-Tanner diffusion pulse sequence, with 16 values of b (0-900 s/mm(2)) in 3 orthogonal directions on 3T systems equipped with 32 multichannel receiver head coils. The intravoxel incoherent motion perfusion parameters were derived by fitting the intravoxel incoherent motion biexponential model. Regions of interest were drawn in regions of maximum intravoxel incoherent motion perfusion fraction and contralateral control regions. Statistical significance was assessed by using the Student t test. In addition, regions of interest were drawn around all whole tumors and were evaluated with the help of histograms. RESULTS: In the regions of maximum perfusion fraction, perfusion fraction was significantly higher in the high-grade group (0.127 ± 0.031) than in the low-grade group (0.084 ± 0.016, P < .001) and in the contralateral control region (0.061 ± 0.011, P < .001). No statistically significant difference was observed for the pseudodiffusion coefficient. The perfusion fraction correlated moderately with dynamic susceptibility contrast relative CBV (r = 0.59). The histograms of the perfusion fraction showed a "heavy-tailed" distribution for high-grade but not low-grade gliomas. CONCLUSIONS: The intravoxel incoherent motion perfusion fraction is helpful for differentiating high- from low-grade brain gliomas.

An iterative Multiscale Finite-Volume algorithm converging to the exact solution

The multiscale finite volume (MsFV) method has been developed to efficiently solve large heterogeneous problems (elliptic or parabolic); it is usually employed for pressure equations and delivers conservative flux fields to be used in transport problems. The method essentially relies on the hypothesis that the (fine-scale) problem can be reasonably described by a set of local solutions coupled by a conservative global (coarse-scale) problem. In most cases, the boundary conditions assigned for the local problems are satisfactory and the approximate conservative fluxes provided by the method are accurate. In numerically challenging cases, however, a more accurate localization is required to obtain a good approximation of the fine-scale solution. In this paper we develop a procedure to iteratively improve the boundary conditions of the local problems. The algorithm relies on the data structure of the MsFV method and employs a Krylov-subspace projection method to obtain an unconditionally stable scheme and accelerate convergence. Two variants are considered: in the first, only the MsFV operator is used; in the second, the MsFV operator is combined in a two-step method with an operator derived from the problem solved to construct the conservative flux field. The resulting iterative MsFV algorithms allow arbitrary reduction of the solution error without compromising the construction of a conservative flux field, which is guaranteed at any iteration. Since it converges to the exact solution, the method can be regarded as a linear solver. In this context, the schemes proposed here can be viewed as preconditioned versions of the Generalized Minimal Residual method (GMRES), with a very peculiar characteristic that the residual on the coarse grid is zero at any iteration (thus conservative fluxes can be obtained).

Self-Navigation with Compressed Sensing for 2D Translational Motion Correction in Free-Breathing Coronary MRI: A Feasibility Study.

PURPOSE: Respiratory motion correction remains a challenge in coronary magnetic resonance imaging (MRI) and current techniques, such as navigator gating, suffer from sub-optimal scan efficiency and ease-of-use. To overcome these limitations, an image-based self-navigation technique is proposed that uses "sub-images" and compressed sensing (CS) to obtain translational motion correction in 2D. The method was preliminarily implemented as a 2D technique and tested for feasibility for targeted coronary imaging. METHODS: During a 2D segmented radial k-space data acquisition, heavily undersampled sub-images were reconstructed from the readouts collected during each cardiac cycle. These sub-images may then be used for respiratory self-navigation. Alternatively, a CS reconstruction may be used to create these sub-images, so as to partially compensate for the heavy undersampling. Both approaches were quantitatively assessed using simulations and in vivo studies, and the resulting self-navigation strategies were then compared to conventional navigator gating. RESULTS: Sub-images reconstructed using CS showed a lower artifact level than sub-images reconstructed without CS. As a result, the final image quality was significantly better when using CS-assisted self-navigation as opposed to the non-CS approach. Moreover, while both self-navigation techniques led to a 69% scan time reduction (as compared to navigator gating), there was no significant difference in image quality between the CS-assisted self-navigation technique and conventional navigator gating, despite the significant decrease in scan time. CONCLUSIONS: CS-assisted self-navigation using 2D translational motion correction demonstrated feasibility of producing coronary MRA data with image quality comparable to that obtained with conventional navigator gating, and does so without the use of additional acquisitions or motion modeling, while still allowing for 100% scan efficiency and an improved ease-of-use. In conclusion, compressed sensing may become a critical adjunct for 2D translational motion correction in free-breathing cardiac imaging with high spatial resolution. An expansion to modern 3D approaches is now warranted.