Affine and polynomial mutual information coregistration for artifact elimination in diffusion tensor imaging of newborns

We have investigated the use of two different image coregistration algorithms for identifying local regions of erroneously high fractional anisotropy (FA) as derived from diffusion tensor imaging (DTI) data sets in newborns. The first algorithm uses conventional affine registration of each of the diffusion-weighted images to the unweighted (b = 0) image for each slice, while the second algorithm uses second-order polynomial warping. Similarity between images was determined using the mutual information (MI) criterion, which is the preferred 'cost' criterion for coregistration of images with significantly different image intensity distributions. We have found that subtle differences exist in the FA values resulting from affine and second-order polynomial coregistration and demonstrate that nonlinear distortions introduce artifacts of spatial extent similar to real white matter structures in the newborn subcortex. We show that polynomial coregistration systematically reduces the presence of erroneous regions of high FA and that such artifacts can be identified by visual inspection of FA maps resulting from affine and polynomial coregistrations. Furthermore, we show that nonlinear distortions may be particularly pronounced when acquiring image slices of axial orientation at the height of the nasal cavity. Finally, we show that third-order polynomial MI coregistration (using the images resulting from second-order coregistration as input) has no observable effect on the resulting FA maps.     

Image registration by maximization of mutual information based on edge width matching using particle swarm optimization

Mutual information (MI) based image registration has been found to be quite effective in many medical image applications. However, standard MI hampers the convergence of registration transformation parameters since it contains local maxima. In this paper, a novel registration method is proposed. At first,MI based on edge width matching is computed to avoid great change of joint probability distribution and get less local maxima. Particle swarm optimization (PSO), which combines local search methods with global ones balancing exploration and exploitation, is done to search the optimal registration parameter.PSO has less computational complexity as its complex behavior follows only a few simple rules. It could avoid local maxima and reach global optimal results. This method is applicable to a variety of multimodal images, and suitable to different interpolation methods. Theoretical analysis and experiments show that this method is effective and accurate to register multimodal medical images.     

Continuous EEG-fMRI in patients with partial epilepsy and focal interictal slow-wave discharges on EEG

Purpose

To verify whether in patients with partial epilepsy and routine electroenecephalogram (EEG) showing focal interictal slow-wave discharges without spikes combined EEG–functional magnetic resonance imaging (fMRI) would localize the corresponding epileptogenic focus, thus providing reliable information on the epileptic source.

Methods

Eight patients with partial epileptic seizures whose routine scalp EEG recordings on presentation showed focal interictal slow-wave activity underwent EEG–fMRI. EEG data were continuously recorded for 24 min (four concatenated sessions) from 18 scalp electrodes, while fMRI scans were simultaneously acquired with a 1.5-Tesla magnetic resonance imaging (MRI) scanner. After recording sessions and MRI artefact removal, EEG data were analyzed offline. We compared blood oxygen level-dependent (BOLD) signal changes on fMRI with EEG recordings obtained at rest and during activation (with and without focal interictal slow-wave discharges).

Results

In all patients, when the EEG tracing showed the onset of focal slow-wave discharges on a few lateralized electrodes, BOLD-fMRI activation in the corresponding brain area significantly increased. We detected significant concordance between focal EEG interictal slow-wave discharges and focal BOLD activation on fMRI. In patients with lesional epilepsy, the epileptogenic area corresponded to the sites of increased focal BOLD signal.

Conclusions

Even in patients with partial epilepsy whose standard EEGs show focal interictal slow-wave discharges without spikes, EEG–fMRI can visualize related focal BOLD activation thus providing useful information for pre-surgical planning.     

Quantitative evaluation of optimal imaging parameters for single-cell detection in MRI using simulation

Super-paramagnetic iron oxide (SPIO) nanoparticles are actively investigated to enhance disease detection through molecular imaging using magnetic resonance imaging (MRI). Detection of the cells labeled by SPIO depends on the MRI protocols and pulse sequence parameters that can be optimized. To evaluate the sensitivity and specificity of the image acquisition methods and to obtain optimal imaging parameters for single-cell detection, we further developed an MRI simulator. The simulator models an object (tissue) at a microscopic level to evaluate effects of spatial distribution and concentration of nanoparticles on the resulting image. In this study, the simulator was used to evaluate and compare imaging of the labeled cells by the gradient-echo (GE), true-FISP [fast imaging employing steady-state acquisition (FIESTA)] and echo-planar imaging (EPI) pulse sequences. Effects of the imaging and object parameters, such as field strength, imaging protocol and pulse sequence parameters, imaging resolution, cell iron load, position of SPIO within the voxel and cell division within the voxel, were investigated in the work. The results suggest that true-FISP has the highest sensitivity for single-cell detection by MRI.     

Quantitative magnetic resonance imaging in consecutive patients evaluated for surgical treatment of temporal lobe epilepsy

We present the results of quantitative Magnetic Resonance Imaging (MRI) in 55 consecutively referred patients with clinical evidence of temporal lobe epilepsy (TLE). The Cavalieri method was used in combination with point counting to provide unbiased estimates of the volume of the left and right hippocampus, amygdala, temporal lobe, lateral ventricles and cerebral hemisphere, and pixel by pixel maps of the T2 relaxation time were computed for both central and anterior sections of the hippocampus. The 99th centiles of hippocampal volume, hippocampal volume asymmetry and T2 relaxation times in 20 control subjects provided limits which identified the presence of MTS. The results of the quantitative MRI were compared with the results of conventional diagnostic MRI, foramen ovale (FO) recording and the WADA test. Thirty-one patients were found to have unilateral MTS (17 left and 14 right) and 7 bilateral MTS. No evidence of MTS was detected in 16 patients. Of the 31 patients diagnosed with unilateral MTS on the basis of hippocampal volume and T2 measurement, 74% and 77% would respectively have received the same diagnosis on the basis of hippocampal volume and T2 measurements alone. In comparison to FO recording, quantitative MRI has a sensitivity of 55% and a specificity of 86%, while conventional diagnostic MRI has a sensitivity of 42% and a specificity of 80% for detection of MTS. Unilateral abnormalities were detected by FO recording in 30% cent of patients who appeared normal on quantitative MRI. WADA test results were available for 40 patients. The findings were consistent with quantitative MRI showing reduced memory function ipsilateral to unilateral MTS in 18 patients, but reduced memory function contralateral to unilateral MTS in two patients, and reduced memory function without MR abnormality in seven patients. WADA testing revealed unilateral memory impairments where MRI found bilateral pathology in 4 patients and in 4 patients in whom quantitative MRI detected unilateral MTS there was no evidence of reduced memory during WADA testing of the corresponding cerebral hemisphere. In the patients with unilateral right MTS a highly significant negative correlation (p = 0.0003) was observed between age of onset and the volume of the contralateral temporal lobe.

Quantitative MR imaging of the hippocampus (i.e. volume and T2 measurement) is preferable to conventional radiological reporting for providing objective evidence of the presence of MTS on which to base the referral of patients for surgery, and since it has associated morbidity FO recording is now only being used in selected patients. Furthermore, stereology provides a convenient method for estimating the volume of other brain structures, which is relevant to obtaining a better understanding of the effects of laterality and age of onset of TLE.     

基于互信息方法的燃油品质光谱分析模型构建与简化

近红外漫反射光谱和紫外吸收光谱分别用于燃油的辛烷值和单芳香族化合物含量的测定,偏最小二乘回归(partial least squares regression,PLSR)用于光谱多元校正模型的构建。基于互信息(mutual infor-mation,MI)理论的变量筛选方法用于模型优化以提高模型的预测精度,降低模型的复杂度。结果表明,MI-PLSR可以有效的提高燃油品质模型的预测精度,简化分析模型。辛烷值的预测均方根误差(root meansquare error of prediction,RMSEP)由0.288减小为0.111,预测相关系数R从0.985提高到0.998,建模变量由401减小为112;单芳香族化合物含量的RMSEP从0.753减小为0.478,R由0.996提高为0.998,建模变量由572缩减为37。说明振动光谱结合MI-PLSR方法可用于燃油品质检测,具有高效率低成本的特点。     

Quantitative measurement of leakage volume and permeability in gliomas, meningiomas and brain metastases with dynamic contrast-enhanced MRI

The spatial properties and function of the tumor vasculature differ with the tumor type and grade. T1-weighted dynamic contrast-enhanced imaging technique enables the simultaneous quantification of some functional parameters of the vasculature. These are the fractional contrast-enhancing volumes of the tissue compartments (blood volume and leakage/extravascular extracellular volume) and the exchange parameters (perfusion and permeability). The relatively long monitoring duration of 12 min used here made it necessary to divide the extravascular extracellular compartment into two subcompartments, a slowly and a fast enhancing one with different permeabilities. Forty-one gliomas (WHO grades II-IV), six meningiomas and eight distant metastases were investigated. It was shown that the technique noninvasively provides information for separating different tumor types and characterizing their microenvironment. Fast permeability describes vessel permeability and was significantly increased in meningiomas as compared with intra-axial tumors. The corresponding volume of the fast enhancing compartment was significantly increased in meningiomas compared to all gliomas taken together. Slow permeability describes diffusion within the extravascular extracellular space and was significantly reduced in low-grade gliomas, indicating short diffusion distances. The slowly enhancing extravascular extracellular space was found to be increased in high-grade gliomas and distant metastases. Blood volume differed significantly among some tumor entities and glioma grades. Perfusion was shown to increase linearly with blood volume for volumes of up to 20%, flattening out thereafter. The scatter plots of extravascular extracellular volume and blood volume were shown to differ among the tumor entities.     

Quantitative human in vivo evaluation of high resolution MRI for vessel wall morphometry after percutaneous transluminal angioplasty

Visualization of the vessel wall after transluminal angioplasty is important to monitor the restenosis progress. Intravascular ultrasound proved its capabilities as an invasive procedure in many studies. The aim of this study was to evaluate the feasibility of high-resolution MRI as a non-invasive tool for follow-up after PTA. High-resolution magnetic resonance images (pixel size: 0.49 * 0.49 mm(2)) were acquired on a 1.0 T clinical scanner. Morphometry was conducted after conversion of DICOM images into TIFF format using ScionImage on a PC. In-vitro studies using a polyvinylchloride tube were evaluated by two independent investigators. Goldstandard was a caliper rule and direct radiography. Five patients were monitored before and 24 h, six weeks, three months and six months after PTA. In vivo measurements promised a good concordance for both investigators for area as well as for diameter measurements. Area measurements showed correlations up to r = 0.86 (p < 0.001) whereas the correlations of diameters were slightly inferior (r between 0.58 and 0.84; p < 0.005). Relocation of the same slice position in the follow up studies could be guaranteed using anatomic landmarks in the images. As a non-invasive tool to assess restenosis after PTA high-resolution MRI promises to be a reproducible technique. It is easy to identify the same vascular region in different studies due to neighboring anatomic landmarks. Progression of disease as well as success of pharmacologic treatment to prevent restenosis may be monitored.     

How does brain MRI lesion volume change on serial scans in patients with multiple sclerosis?

Although lesion load changes on conventional T₂-weighted brain magnetic resonance imaging (MRI) scans from patients with multiple sclerosis (MS) are used to monitor the effect of treatment, there is no clear definition of how lesion load changes over years according to the lesion load present at a baseline evaluation. In the present study, we evaluated the relationship between lesion load changes over time and lesion load at a baseline evaluation in a group of untreated patients with MS. We scanned nineteen patients on two separate occasions with a mean interval 16.4 months between the two examinations. In each scanning session, a scan with forty contiguous 3-mm-thick axial slices was acquired. We assessed MRI lesion loads using a semi-automated local thresholding technique. Both a linear (p < 0.0001) and a quadratic component (p = 0.0008) of the baseline volume were significant in describing the follow-up volume. The equation to model this finding was as follows: V_f = β₀ V_b + β₁ (V_b)², where V_f is the lesion volume at follow-up, V_b is the lesion volume at baseline, β₀ = 0.834 (SE = 0.098), and β₁ = 0.014 (SE = 0.003) (mL)⁻¹. Our data indicate that lesion volume changes detectable on serial brain MRI studies from patients with MS are dependent on the extent of lesion burden present on the baseline MRI scans. This finding has to be considered when planning phase III trials.     

Optimizing brain MRI protocols in the follow-up of patients with multiple sclerosis T2-weighted MRI of the brain after the administration of gadopentetate dimeglumine

The aim of our study was to determine whether T2-weighted (T2w) MRI of the brain could be performed immediately after the administration of gadopentetate dimeglumine (gadolinium DTPA) in patients with multiple sclerosis (MS) without a loss in image quality or diagnostic reliability. Sixteen patients with clinically diagnosed MS were included in the study. Twenty-four patients with various cerebral pathologies (14 patients with multiple lacunar lesions) were examined in order to exclude masking of T2 hyperintense lesions other than MS lesions. Images of 10 patients without pathological changes served as a control condition for the qualitative analysis. In these 50 patients, T1w and T2w MRI was performed before and after the administration of gadolinium DTPA. Signal intensities were measured within T2 hyperintense cerebral lesions, in T1-enhancing lesions and in normal appearing brain tissue on T2w turbo spin-echo (TSE) sequences. Both quantitative and qualitative analysis did not show significant differences between T2w pre- and postcontrast series. T2w MRI performed prior to and after the administration of gadolinium DTPA provides similar information in patients with MS. With a TR of 3.2 s, not a single lesion was obscured on T2w postcontrast series. Acquisition of T2w MR images immediately after the administration of gadolinium DTPA allows for shorter examination time and assures sufficient time for contrast enhancement in cerebral lesions with a disrupted blood-brain barrier.     

Fractality in the neuron axonal topography of the human brain based on 3-D diffusion MRI

We report our results on the effect of incorporation of inorganic fullerene like nanoparticles (IF) and inorganic nanotubes (INT) of WS₂ into hybrid LED device structures. To disperse into a uniform fashion, the semiconducting INT/IF WS₂ NTs were functionalized with SDS (sodium dodecylsulphate). The IF/INT WS₂ nanotubes were used in combination with PEDOT:PSS and P3HT to realize the following LED device structures: ITO/(PEDOT:PSS):(WS₂:SDS)/P3HT/LiF-Al; ITO/PEDOT:PSS/P3HT/WS₂:SDS/LiF-Al. Morphological, optical and electrochemical analysis were performed to obtain the HOMO and the LUMO energy levels to hypothesize the most efficient device structure. The spectral positions of the electroluminescent bands were found out to be device-dependent and exhibits blue shift when the proposed nanostructure is dip coated on top of P3HT. Electro-optical analysis indicate that the WS₂:SDS based P3HT/semiconductor film can improve the charge recombination probability owing to its dual functionality as hole blocking layer and electron injection moiety.     

Automatic segmentation of brain MRI in high-dimensional local and non-local feature space based on sparse representation

Automatic extraction of the varying regions of magnetic resonance images is required as a prior step in a diagnostic intelligent system. The sparsest representation and high-dimensional feature are provided based on learned dictionary. The classification is done by employing the technique that computes the reconstruction error locally and non-locally of each pixel. The acquired results from the real and simulated images are superior to the best MRI segmentation method with regard to the stability advantages. In addition, it is segmented exactly through a formula taken from the distance and sparse factors. Also, it is done automatically taking sparse factor in unsupervised clustering methods whose results have been improved.     

Comparison of analytical strategies for EEG-correlated fMRI data in patients with epilepsy

The simultaneous recording of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) can be used to localize interictal epileptiform discharges (IEDs). Previous studies have reported varying degrees of concordance of EEG-fMRI with electroclinical findings. The aim of the present study is to evaluate to what extent this variability is determined by the analytical strategy or by the properties of the EEG data. For that purpose, 42 IED sets obtained in 29 patients with epilepsy were reanalyzed using a finite impulse response approach, which estimates the hemodynamic response function (HRF) from the data and allows non-causal effects. Cardiac effects were treated as additional confounders in the model. This approach was compared to the classical approach assuming a fixed HRF for each voxel in the brain. The performance of each method was assessed by comparing the fMRI results to the EEG focus. The flexible model revealed more significantly activated voxels, which resulted in more activated brain regions concordant with the EEG focus (26 vs. 16). Correction for cardiac effects improved the results in 7 out of the 42 data sets. Furthermore, design theory for event-related experiments was applied in order to determine the influence of the number of IEDs and their temporal distribution on the success of an experiment. It appeared that this success is highly dependent upon the number of IEDs present during the recording and less on their temporal spacing. We conclude that the outcome of EEG-fMRI can be improved by using an optimized analytical strategy, but also depends on the number of IEDs occurring during the recording.     

Small hyperintense hepatic lesions on T1-weighted images in patients with cirrhosis: evaluation with serial MRI and imaging features for clinical benignity

PURPOSE: The aim of this study was to evaluate the frequency and magnetic resonance imaging (MRI) features of clinically benign, small (<2 cm) hyperintense hepatic lesions in the cirrhotic liver on T1-weighted MR images seen at serial MRI. MATERIALS AND METHODS: This study included 189 patients with cirrhosis, who underwent hepatic MRI more than twice with an interval of at least 12 months. The initial MR images were reviewed for the presence of small hyperintense lesions on T1-weighted images. The size, location and signal intensity on T2-weighted images as well as enhancement patterns of the corresponding lesions were recorded. RESULTS: On the initial T1-weighted MR images, 43 small hyperintense hepatic lesions were detected in 23 (12%) of 189 patients. Twelve (28%) of 43 lesions showed early enhancement and were pathologically diagnosed as hepatocellular carcinoma (HCC) during the follow-up period. Thirty-one (72%) of 43 lesions showed no early enhancement with various signal intensity on T2-weighted images (hyperintensity=4, isointensity=20, hypointensity=7). Among these 31 lesions, 12 showed no interval change, while 11 disappeared (n=10) or decreased in size (n=1). In the remaining eight lesions, seven were diagnosed as HCC on the basis of pathologic confirmation or the interval growth. CONCLUSION: Small hyperintense hepatic lesions on T1-weighted magnetic resonance (MR) images without early enhancement on the arterial-phase contrast-enhanced dynamic studies in patients with cirrhosis usually showed no interval growth or disappeared during the serial MRI. These lesions with additional findings of iso- or hypointensity on the T2-weighted MR images without "washout effect" on the contrast-enhanced equilibrium-phase images may more frequently be clinically benign or hyperplastic nodules than HCCs.     

Technical aspects and evaluation methodology for the application of two automated brain MRI tumor segmentation methods in radiation therapy planning

The purpose of this study was to design the steps necessary to create a tumor volume outline from the results of two automated multispectral magnetic resonance imaging segmentation methods and integrate these contours into radiation therapy treatment planning. Algorithms were developed to create a closed, smooth contour that encompassed the tumor pixels resulting from two automated segmentation methods: k-nearest neighbors and knowledge guided. These included an automatic three-dimensional (3D) expansion of the results to compensate for their undersegmentation and match the extended contouring technique used in practice by radiation oncologists. Each resulting radiation treatment plan generated from the automated segmentation and from the outlining by two radiation oncologists for 11 brain tumor patients was compared against the volume and treatment plan from an expert radiation oncologist who served as the control. As part of this analysis, a quantitative and qualitative evaluation mechanism was developed to aid in this comparison. It was found that the expert physician reference volume was irradiated within the same level of conformity when using the plans generated from the contours of the segmentation methods. In addition, any uncertainty in the identification of the actual gross tumor volume by the segmentation methods, as identified by previous research into this area, had small effects when used to generate 3D radiation therapy treatment planning due to the averaging process in the generation of margins used in defining a planning target volume.     

Estimation of mutual information objective function based on Fourier shift theorem: an application to eddy current distortion correction in diffusion tensor imaging

Diffusion tensor imaging requires correction of eddy current distortion in diffusion-weighted images. An effective retrospective correction approach is to transform a diffusion-weighted image to maximize the mutual information (MI) between the transformed diffusion-weighted image and the corresponding T2-weighted image. In the literature, either linear interpolation or partial volume interpolation is applied to estimate the MI objective function. However, these interpolation methods induce artifacts to the MI objective function, thus compromising correction results. In this work, the MI objective function is estimated based on interpolation using Fourier shift theorem. This method eliminates the artifacts incurred with the aforementioned interpolation methods. The algorithm is further improved by approximating pixel values using their nearest neighbors in the up-sampled spatial domain, resulting in dramatically increased computational efficiency without compromising the correction results. The effects of varying the number of quantization levels and using Parzen window filtering to smooth the MI objective function are also investigated to obtain optimized algorithm parameters. The diffusion tensor image quality after applying the proposed distortion correction method is significantly improved visually.     

A non-invasive technique for 3-dimensional assessment of articular cartilage thickness based on MRI part 2: Validation using CT arthrography

Established methods for the measurement of articular cartilage thickness are invasive and cannot be sequentially applied in living subjects. In the present study, the distribution of cartilage thickness throughout entire joint surfaces was determined from MR images obtained with a fat-suppressed gradient-echo sequence at a resolution of 0.31 × 0.31 × 2.00 mm³, and compared to that derived from CT arthrography. A minimal distance algorithm was employed to produce 3D cartilage thickness maps of seven cadaveric human knee joints. The mean amount of deviation of the cartilage volumes was 5.6% (±4.6), statistical analysis showing that there was high agreement between the two methods (r = 0.995, slope = 1.037, y-intercept = -90.5 mm³). The 3D thickness maps yielded a striking agreement between the two methods, the maximum values generally yielding a deviation of none or one thickness interval of 0.5 mm. This investigation shows that accurate 3D assessment of articular cartilage thickness can be performed with MRI, this technique having the advantage that it is suitable for investigating living subjects.     

Effect of the mutual rectification of two electromagnetic waves with perpendicular polarization planes in a superlattice based on graphene

The effect of the mutual rectification of two electromagnetic waves with perpendicular polarization planes in a superlattice based on graphene was investigated. It was shown that, due to the nonadditive nature of the electronic spectrum of a superlattice based on grapheme, a constant term of current arises along its axes. The dependence of the constant of current density on the intensity of an electromagnetic wave is oscillatory in nature.     

基于1/f噪声的NPN晶体管辐照感生电荷的定量分离

本文针对NPN双极性晶体管, 在研究辐照感生的氧化层电荷及界面态对晶体管基极电流和1/f噪声的影响的基础上, 建立辐照感生氧化层电荷及界面态与基极电流和1/f噪声的定量物理模型. 根据所建立的模型, 提出一种新的分离方法, 利用1/f噪声和表面电流求出氧化层电荷密度, 利用所求得氧化层电荷密度和表面电流求出界面态密度. 利用本方法初步实现了辐照感生氧化层电荷及界面态的定量计算.     

Feasibility of similarity coefficient map in improving morphological evaluation of T2* weighted MRI for renal cancer

The purpose of this paper is to investigate the feasibility of similarity coefficient map (SCM) in improving morphological evaluation of T₂^* weighted (T₂^*W) magnatic resonance imaging (MRI) for renal cancer. Simulation studies and in vivo 12-echo T₂^*W experiments for renal cancers were performed for this purpose. The results of the first simulation study suggest that SCM can reveal small structures which are hard to be distinguished from the background tissue in T₂^*W images and the corresponding T₂^* map. The capability of improving morphological evaluation is likely due to the improvement in signal to noise ratio (SNR) and carrier to noise ratio (CNR) by SCM technique. Compared with T₂^*W images, SCM can improve SNR by a factor ranging from 1.87 to 2.47. Compared with T₂^* maps, SCM can improve SNR by a factor ranging from 3.85 to 33.31. Compared with T₂^*W images, SCM can improve CNR by a factor raging from 2.09 to 2.43. Compared with T₂^* maps SCM can improve CNR by a factor raging from 1.94 to 8.14. For a given noise level, the improvements of SNR and CNR depend mainly on the original SNRs and CNRs in T₂^*W images, respectively. In vivo experiments confirmed the results of the first simulation study. The results of the second simulation study suggest that more echoes are used to generate SCM, and higher SNR and CNR can be achieved in SCM. In conclusion, SCM can provide improved morphological evaluation of T₂^*W MR images for renal cancer by unveiling fine structures which are ambiguous or invisible in the corresponding T₂^*W MR images and T₂^* maps. What is more, in practical application, for a fixed total sampling time, one should increase the number of echoes as much as possible to achieve SCMs with better SNR and CNR.