Brain activation difference evoked by different binocular disparities of stereograms: An fMRI study.

The binocular disparity of two retina images is a main cue of stereoscopic vision. However, the global dependency between brain response and binocular disparity still remains unclear. Here, we used functional Magnetic Resonance Imaging (fMRI) to identify stereopsis-related brain regions with a modified Random Dot Stereogram (RDS) and plotted the activation variation curves under different disparity size. In order to eliminate the confounding shape difference between the stereogram and the plane, commonly seen in RDS, we modified the RDS to a checkerboard version. We found that V3A, V7 and MT+/V5 in dorsal visual stream were activated in stereoscopic experiment, while little activation was found in ventral visual regions. According to the activation trends, 13 subjects were divided into three groups: 5 subjects with turning points (a shift from increased to decreased activation), 5 subjects without turning points and 3 subjects with activation unrelated to disparity. We inferred that the dorsal visual stream primarily processes spatial depth information, rather than shape information.

Concomitant Imaging Dose and Cancer Risk in Image Guided Thoracic Radiation Therapy.

PURPOSE: Kilovoltage cone beam computed tomography (CT) (kVCBCT) imaging guidance improves the accuracy of radiation therapy but imposes an extra radiation dose to cancer patients. This study aimed to investigate concomitant imaging dose and associated cancer risk in image guided thoracic radiation therapy. METHODS AND MATERIALS: The planning CT images and structure sets of 72 patients were converted to CT phantoms whose chest circumferences (Cchest) were calculated retrospectively. A low-dose thorax protocol on a Varian kVCBCT scanner was simulated by a validated Monte Carlo code. Computed doses to organs and cardiac substructures (for 5 selected patients of various dimensions) were regressed as empirical functions of Cchest, and associated cancer risk was calculated using the published models. The exposures to nonthoracic organs in children were also investigated. RESULTS: The structural mean doses decreased monotonically with increasing Cchest. For all 72 patients, the median doses to the heart, spinal cord, breasts, lungs, and involved chest were 1.68, 1.33, 1.64, 1.62, and 1.58 cGy/scan, respectively. Nonthoracic organs in children received 0.6 to 2.8 cGy/scan if they were directly irradiated. The mean doses to the descending aorta (1.43 ± 0.68 cGy), left atrium (1.55 ± 0.75 cGy), left ventricle (1.68 ± 0.81 cGy), and right ventricle (1.85 ± 0.84 cGy) were significantly different (P<.05) from the heart mean dose (1.73 ± 0.82 cGy). The blade shielding alleviated the exposure to nonthoracic organs in children by an order of magnitude. CONCLUSIONS: As functions of patient size, a series of models for personalized estimation of kVCBCT doses to thoracic organs and cardiac substructures have been proposed. Pediatric patients received much higher doses than did the adults, and some nonthoracic organs could be irradiated unexpectedly by the default scanning protocol. Increased cancer risks and disease adverse events in the thorax were strongly related to higher imaging doses and smaller chest dimensions.

Free-Breathing 3D Imaging of Right Ventricular Structure and Function Using Respiratory and Cardiac Self-Gated Cine MRI.

Providing a movie of the beating heart in a single prescribed plane, cine MRI has been widely used in clinical cardiac diagnosis, especially in the left ventricle (LV). Right ventricular (RV) morphology and function are also important for the diagnosis of cardiopulmonary diseases and serve as predictors for the long term outcome. The purpose of this study is to develop a self-gated free-breathing 3D imaging method for RV quantification and to evaluate its performance by comparing it with breath-hold 2D cine imaging in 7 healthy volunteers. Compared with 2D, the 3D RV functional measurements show a reduction of RV end-diastole volume (RVEDV) by 10%, increase of RV end-systole volume (RVESV) by 1.8%, reduction of RV systole volume (RVSV) by 21%, and reduction of RV ejection fraction (RVEF) by 12%. High correlations between the two techniques were found (RVEDV: 0.94; RVESV: 0.85; RVSV: 0.95; and RVEF: 0.89). Compared with 2D, the 3D image quality measurements show a small reduction in blood SNR, myocardium-blood CNR, myocardium contrast, and image sharpness. In conclusion, the proposed self-gated free-breathing 3D cardiac cine imaging technique provides comparable image quality and correlated functional measurements to those acquired with the multiple breath-hold 2D technique in RV.

Functional and anatomical properties of human visual cortical fields.

Human visual cortical fields (VCFs) vary in size and anatomical location across individual subjects. Here, we used functional magnetic resonance imaging (fMRI) with retinotopic stimulation to identify VCFs on the cortical surface. We found that aligning and averaging VCF activations across the two hemispheres provided clear delineation of multiple retinotopic fields in visual cortex. The results show that VCFs have consistent locations and extents in different subjects that provide stable and accurate landmarks for functional and anatomical mapping. Interhemispheric comparisons revealed minor differences in polar angle and eccentricity tuning in comparable VCFs in the left and right hemisphere, and somewhat greater intersubject variability in the right than left hemisphere. We then used the functional boundaries to characterize the anatomical properties of VCFs, including fractional anisotropy (FA), magnetization transfer ratio (MTR) and the ratio of T1W and T2W images and found significant anatomical differences between VCFs and between hemispheres.

Distributed CT image reconstruction algorithm based on the alternating direction method.

With the development of compressive sensing theory, image reconstruction from few-view projections has been paid considerable research attention in the field of computed tomography (CT). Total variation (TV)-based CT image reconstruction has been shown experimentally to be capable of producing accurate reconstructions from sparse-view data. Motivated by the need of solving few-view reconstruction problem with large scale data, a general block distribution reconstruction algorithm based on TV minimization and the alternating direction method (ADM) has been developed in this study. By utilizing the inexact ADM, which involves linearization and proximal point techniques, the algorithm is relatively simple and hence convenient for the derivation and distributed implementation. And because the data as well as the computation are distributed to individual nodes, an outstanding acceleration factor is achieved. Experimental results demonstrate that the proposed method can accelerate the alternating direction total variation minimization (ADTVM) algorithm with nearly no loss of accuracy, which means compared with ADTVM, the proposed algorithm has a better accuracy with same running time.

Interior reconstruction method based on rotation-translation scanning model.

In various applications of computed tomography (CT), it is common that the reconstructed object is over the field of view (FOV) or we may intend to sue a FOV which only covers the region of interest (ROI) for the sake of reducing radiation dose. These kinds of imaging situations often lead to interior reconstruction problems which are difficult cases in the reconstruction field of CT, due to the truncated projection data at every view angle. In this paper, an interior reconstruction method is developed based on a rotation-translation (RT) scanning model. The method is implemented by first scanning the reconstructed region, and then scanning a small region outside the support of the reconstructed object after translating the rotation centre. The differentiated backprojection (DBP) images of the reconstruction region and the small region outside the object can be respectively obtained from the two-time scanning data without data rebinning process. At last, the projection onto convex sets (POCS) algorithm is applied to reconstruct the interior region. Numerical simulations are conducted to validate the proposed reconstruction method.

Ultrashort echo time bi-component analysis of cortical bone--a field dependence study.

PURPOSE: The purpose of this study is to investigate the effect of differing field strength on the T2* of cortical bone at 1.5 T and 3 T. METHODS: Ultrashort echo time pulse sequences were used to study six bovine and nine human bone samples at 1.5 T and 3 T using single- and bi-component T2* analysis. RESULTS: On average, the bound water T2* of bovine bone decreased by 16% (from 0.32 ms at 1.5 T to 0.27 ms at 3 T, P < 0.01) and the bound water T2* of human bone decreased by 21% (from 0.42 ms at 1.5 T to 0.33 ms at 3 T, P < 0.01) at the higher field strength. The free water T2* of bovine bone decreased by 50% (from 4.23 ms at 1.5 T to 2.12 ms at 3 T, P < 0.001) and the free water T2* of human bone decreased by 68% (from 7.65 ms at 1.5 T to 2.46 ms at 3 T, P < 0.001) at the higher field strength. Bound and free water fractions showed only minor change with field strength in bovine (< 2%, P > 0.05) and human bone (< 4%, P > 0.05). CONCLUSION: Ultrashort echo time bi-component analysis provides consistent bound and free water fractions at 1.5 T and 3 T, thereby allowing field-independent comparisons.

Personalized estimation of dose to red bone marrow and the associated leukaemia risk attributable to pelvic kilo-voltage cone beam computed tomography scans in image-guided radiotherapy.

The aim of this study is to investigate the imaging dose to red bone marrow (RBM) and the associated leukaemia risks attributable to pelvic kilo-voltage cone beam computed tomography (kVCBCT) scans in image-guided radiation therapy (IGRT). The RBM doses of 42 patients (age 2.7-86.4 years) were calculated using Monte Carlo simulations. The trabecular spongiosa was segmented to substitute RBM rather than the whole bone. Quantitative correlations between anthropometric variables such as age, physical bone density (PBD) and RBM dose were established. Personalized leukaemia risk was evaluated using an improved Boice model which included the age-associated RBM involvement. An incremental leukaemia risk of 29%-82% (mean = 45%) was found to be associated with 40 pelvic kVCBCT scans in the subject group used in a typical external beam radiation therapy course. Higher risks were observed in children. Due to the enhanced photoelectric effect in high atomic number materials, PBD was observed to strongly affect the RBM dose. Considerable overestimations (9%-42%, mean = 28%) were observed if the whole bone doses were used as surrogates of RBM doses. The personalized estimation of RBM dose and associated leukaemia risk caused by pelvic kVCBCT scans is clinically feasible with the proposed empirical models. Higher radiogenic cancer risks are associated with repeated kVCBCT scans in IGRT of cancer patients, especially children.

Fast reconstruction of flat region in a super-short scan based on MD-FBP algorithm.

In circular cone-beam computed tomography (CT), although the minimum data filtered-backprojection (MD-FBP) algorithm has many significant applications, such as handling super-short scan problem, its reconstruction efficiency is limited by the heavy calculation of backprojection. In this paper, aiming at the image reconstruction of flat region in a super-short scan, an improved method based on MD-FBP algorithm is developed using an integral operation with fixed integral interval during the implementation of backprojection, which has an improvement in reconstruction efficiency and parallel performance compared with the original MD-FBP algorithm. It is found that if the thickness of the flat region is less than 0.0349 R (R is the scanning radius), the uncertainty of the method can be ignored. When the thickness of reconstructed region is a little fat, it can also be reconstructed by increasing the scanning radius befittingly. The results of numerical simulation and real data experiments have demonstrated the correctness and merits of the proposed method.

Personalized assessment of kV cone beam computed tomography doses in image-guided radiotherapy of pediatric cancer patients.

PURPOSE: To develop a quantitative method for the estimation of kV cone beam computed tomography (kVCBCT) doses in pediatric patients undergoing image-guided radiotherapy. METHODS AND MATERIALS: Forty-two children were retrospectively analyzed in subgroups of different scanned regions: one group in the head-and-neck and the other group in the pelvis. Critical structures in planning CT images were delineated on an Eclipse treatment planning system before being converted into CT phantoms for Monte Carlo simulations. A benchmarked EGS4 Monte Carlo code was used to calculate three-dimensional dose distributions of kVCBCT scans with full-fan high-quality head or half-fan pelvis protocols predefined by the manufacturer. Based on planning CT images and structures exported in DICOM RT format, occipital-frontal circumferences (OFC) were calculated for head-and-neck patients using DICOMan software. Similarly, hip circumferences (HIP) were acquired for the pelvic group. Correlations between mean organ doses and age, weight, OFC, and HIP values were analyzed with SigmaPlot software suite, where regression performances were analyzed with relative dose differences (RDD) and coefficients of determination (R(2)). RESULTS: kVCBCT-contributed mean doses to all critical structures decreased monotonically with studied parameters, with a steeper decrease in the pelvis than in the head. Empirical functions have been developed for a dose estimation of the major organs at risk in the head and pelvis, respectively. If evaluated with physical parameters other than age, a mean RDD of up to 7.9% was observed for all the structures in our population of 42 patients. CONCLUSIONS: kVCBCT doses are highly correlated with patient size. According to this study, weight can be used as a primary index for dose assessment in both head and pelvis scans, while OFC and HIP may serve as secondary indices for dose estimation in corresponding regions. With the proposed empirical functions, it is possible to perform an individualized quantitative dose assessment of kVCBCT scans.

Synthesis, radiolabeling, biodistribution and fluorescent imaging of histidine-coupled hematoporphyrin.

INTRODUCTION: Hematoporphyrin (Hp) and hematoporphyrin derivatives (HpDs) have been widely used as photosensitizers in photodynamic therapy (PDT). Radiolabeling of HpDs is helpful for preclinical and clinical studies of PDT. METHODS: The histidine-coupled hematoporphyrin (His-Hp) was synthesized and radiolabeled with [(99m)Tc(CO)(3)(H(2)O)(3)](+). Biodistribution of the radioligand and fluorescent imaging of His-Hp in mice bearing S180 tumor were investigated. RESULTS: [(99m)Tc(CO)(3)](+)-labeled His-Hp was electrically neutral, hydrophilic and stable. The biodistribution of the radioligand in S180 tumor-bearing mice was similar with that of nonlabeled HpD in the literature. The uptake of His-Hp in tumors and livers was confirmed by fluorescent imaging. CONCLUSIONS: The complex [(99m)Tc(CO)(3)](+)-His-Hp might be suitable for in vivo dose evaluation of HpD in PDT.

Investigation of motion artifacts associated with fat saturation technique in 3D flash imaging.

PURPOSE: Fast low-angle shot (FLASH) imaging is widely used in dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) because it permits fast and accurate T1 measurement in vivo. Suppression of the fat signal is necessary for most FLASH applications; otherwise, fat will appear hyperintense. The fat saturation technique is one popular method to reduce fat images on clinical scanners. However, fat saturation combined with the 3D FLASH sequence in breast DCE-MRI scans results in heavy ghosting artifacts caused by heartbeat. We used simulation and experimental scans to determine the cause of these artifact-enhancement phenomena. METHODS: We simulated imaging of motion in the x, y, and z directions, with and without fat saturation, to investigate the origin of artifacts. Fourier transform (FT) of the whole field of view was used in the simulation, and we assumed that the uniform phantom was static during one TR. The amplitude of each echo was considered a factor in the FT data. Images were reconstructed using FT data from different phantom positions multiplied by the amplitude factor. Phantom experiments and volunteer studies were implemented to verify the conclusion. RESULTS: Both phantom and volunteer results showed artifacts similar to those in simulation images. We found that FLASH sequence without fat saturation is insensitive to motion. Fat saturation radiofrequency pulses placed before each group of echoes disrupted the steady state of the signal amplitude and produced a low-pass filter effect that enhanced the motion artifacts. CONCLUSIONS: We conclude that the low-pass filter effect associated with the fat saturation technique is responsible for dramatically increased motion artifacts.

Flow compensation for the fast spin echo triple-echo Dixon sequence.

The fast spin echo (FSE) triple-echo Dixon (FTED) sequence uses bipolar triple-echo readout during each echo-spacing period of FSE to collect all the images necessary for Dixon water and fat separation in a single scan. In comparison to other FSE implementations of the Dixon technique, the triple echo readout used in FTED incurs minimal deadtime in the pulse sequence design and thus greatly enhances the overall scan efficiency. A potential drawback of FTED is that the time dependence of the gradient moment along the frequency encode direction becomes more complicated than in FSE and flow compensation based on the gradient moment (GM) nulling is difficult to achieve. In this work, the first order GM along the frequency encode direction of FTED was examined and two different methods to minimize the GM were proposed. The first method nulls the GM at all the locations of the refocusing radiofrequency pulses so that the Carr-Purcell Meiboom-Gill condition is always maintained. The second method minimizes the GM of the spin echo component of the FSE signal at the echo locations. The efficacy of both methods in reducing the first order GM and flow-related artefacts was demonstrated both in phantom and in images in vivo.

Feasibility of high temporal resolution breast DCE-MRI using compressed sensing theory.

PURPOSE: To investigate the feasibility of high temporal resolution breast DCE-MRI using compressed sensing theory. METHODS: Two experiments were designed to investigate the feasibility of using reference image based compressed sensing (RICS) technique in DCE-MRI of the breast. The first experiment examined the capability of RICS to faithfully reconstruct uptake curves using undersampled data sets extracted from fully sampled clinical breast DCE-MRI data. An average approach and an approach using motion estimation and motion compensation (ME/MC) were implemented to obtain reference images and to evaluate their efficacy in reducing motion related effects. The second experiment, an in vitro phantom study, tested the feasibility of RICS for improving temporal resolution without degrading the spatial resolution. RESULTS: For the uptake-curve reconstruction experiment, there was a high correlation between uptake curves reconstructed from fully sampled data by Fourier transform and from undersampled data by RICS, indicating high similarity between them. The mean Pearson correlation coefficients for RICS with the ME/MC approach and RICS with the average approach were 0.977 +/- 0.023 and 0.953 +/- 0.031, respectively. The comparisons of final reconstruction results between RICS with the average approach and RICS with the ME/MC approach suggested that the latter was superior to the former in reducing motion related effects. For the in vitro experiment, compared to the fully sampled method, RICS improved the temporal resolution by an acceleration factor of 10 without degrading the spatial resolution. CONCLUSIONS: The preliminary study demonstrates the feasibility of RICS for faithfully reconstructing uptake curves and improving temporal resolution of breast DCE-MRI without degrading the spatial resolution.

Quantitative analysis of clinical dynamic contrast-enhanced MR imaging for evaluating treatment response in human breast cancer.

PURPOSE: To develop a method that combines a fixed-T1, fuzzy c-means (FCM) technique with a reference region (RR) model (T1-FCM method) to estimate pharmacokinetic parameters without measuring the arterial input function or baseline T1, or T1(0), and to demonstrate its feasibility in the assessment of treatment response to neoadjuvant chemotherapy (NAC) in patients with breast cancer by using data from dynamic contrast material-enhanced magnetic resonance (MR) imaging. MATERIALS AND METHODS: This study was approved by the human investigation committees of the two participating institutions. All patients gave written informed consent. A conventional dual-flip-angle gradient-echo method was used to evaluate the effects of noise and the T1 in the tissue itself on the accuracy of T1 estimation. Both conventional RR and fixed-T1 methods were used to evaluate the effects of noise and preselected T1(0) on the estimation of pharmacokinetic parameters by means of a simulation study. Thirty-three women (age range, 32-66 years; mean age, 45 years) with pathologically proved breast tumors were examined to evaluate the feasibility of using the T1-FCM method as a means of assessing treatment response to NAC. A nonparametric Mann-Whitney U test was used to assess the difference in each of the MR imaging parameters between patients with a major histologic response to treatment and those with a nonmajor histologic response. RESULTS: With use of the dual-flip-angle method, the accuracy and distribution of T1 estimation are dependent on the T1 in the tissue itself. The T1-FCM method is more accurate than other methods and is relatively insensitive to the effects of noise and incorrect T1(0) selection. Preliminary clinical data revealed a significant difference (P < .01) in the change of the volume transfer constant after two cycles of NAC between the major and nonmajor histologic response groups. CONCLUSION: Results of the simulation study demonstrate that the T1-FCM method appears to be relatively insensitive to noisy dynamic contrast-enhanced MR imaging data. This method could prove useful in the evaluation of breast cancer therapy.

Fast finite Hilbert transform via DE quadrature scheme.

PURPOSE: The finite Hilbert transform (FHT) or inverse finite Hilbert transform (IFHT) is recently found to have some important applications in computerized tomography (CT) arena, where they are used to filter the derivatives of back-projected data in the chord-line based CT reconstruction algorithms. In this paper, we implemented, improved and validated a fast numerical solution to the FHT via a double exponential (DE) integration scheme. A same strategy can be used to compute IFHT. METHODS: To overcome the underflow of floating-point numbers, we first determined the range of variable transformation from the minimum positive value of single or double precision floating point number, the integration step can be further determined by the range of variable transformation and the integration level. Two functions with their known analytical FHTs are used to validate the implementation of the FHT via DE scheme. The surface map and 2D contour of the FHT transformation error with respect to integration level and the range of the variable transformation are used to numerically determine the optimal numbers for a fast FHT. RESULTS: Given a specific precision, the lowest integration level and the optimal range of variable transformation, which are used to transform a signal with a certain degree of fluctuation, can be numerically determined by the surface map and 2D contour of the standard deviation of transformation error. These two numbers can then be taken to efficiently compute the FHT for other signals with the same or less degree of fluctuation. CONCLUSIONS: The FHT via DE scheme and the numerical method to determine the integration level and the range of transformation can be used for fast FHT in certain applications, such as data filtering in chord-line based CT reconstruction algorithms.

A clinically feasible method to estimate pharmacokinetic parameters in breast cancer.

Dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) is the MRI technique of choice for detecting breast cancer, which can be roughly classified as either quantitative or semiquantitative. The major advantage of quantitative DCE-MRI is its ability to provide pharmacokinetic parameters such as volume transfer constant (Ktrans) and extravascular extracellular volume fraction (ve). However, semiquantitative DCE-MRI is still the clinical MRI technique of choice for breast cancer diagnosis due to several major practical difficulties in the implementation of quantitative DCE-MRI in a clinical setting, including (1) long acquisition necessary to acquire 3D T1(0) map, (2) challenges in obtaining accurate artery input function (AIF), (3) long computation time required by conventional nonlinear least square (NLS) fitting, and (4) many illogical values often generated by conventional NLS method. The authors developed a new analysis method to estimate pharmacokinetic parameters Ktrans and ve from clinical DCE-MRI data, including fixed T1(0) to eliminate the long acquisition for T1(0) map and "reference region" model to remove the requirement of measuring AIF. Other techniques used in our analysis method are (1) an improved formula to calculate contrast agent (CA) concentration based on signal intensity of SPGR data, (2) FCM clustering-based techniques for automatic segmentation and generation of a clustered concentration data set (3) an empirical formula for CA time course to fit the clustered data sets, and (4) linear regression for the estimation of pharmacokinetic parameters. Preliminary results from computer simulation and clinical study of 39 patients have demonstrated (1) the feasibility of their analysis method for estimating Ktrans and ve from clinical DCE-MRI data, (2) significantly less illogical values compared to NLS method (typically less than 1% versus more than 7%), (3) relative insensitivity to the noise in DCE-MRI data; (4) reduction in computation time by a factor of more than 30 times compared to NLS method on average, (5) high statistic correlation between the method used and NLS method (correlation coefficients: 0.941 for Ktrans and 0.881 for ve), and (6) the potential clinical usefulness of the new method.

Reliable two-dimensional phase unwrapping method using region growing and local linear estimation.

In MRI, phase maps can provide useful information about parameters such as field inhomogeneity, velocity of blood flow, and the chemical shift between water and fat. As phase is defined in the (-pi,pi] range, however, phase wraps often occur, which complicates image analysis and interpretation. This work presents a two-dimensional phase unwrapping algorithm that uses quality-guided region growing and local linear estimation. The quality map employs the variance of the second-order partial derivatives of the phase as the quality criterion. Phase information from unwrapped neighboring pixels is used to predict the correct phase of the current pixel using a linear regression method. The algorithm was tested on both simulated and real data, and is shown to successfully unwrap phase images that are corrupted by noise and have rapidly changing phase.