CUDA accelerated method for motion correction in MR PROPELLER imaging

In PROPELLER, raw data are collected in N strips, each locating at the center of k-space and consisting of M_x sampling points in frequency encoding direction and L lines in phase encoding direction. Phase correction, rotation correction, and translation correction are used to remove artifacts caused by physiological motion and physical movement, but their time complexities reach O(M_x × M_x × L × N), O(N × R_A × M_x × L × (M_x × L + R_N × R_N)), and O(N × (R_N × R_N + M_x × L)) where R_N × R_N is the coordinate space each strip gridded onto and R_A denotes the rotation range. A CUDA accelerated method is proposed in this paper to improve their performances. Although our method is implemented on a general PC with Geforce 8800GT and Intel Core(TM)2 E6550 2.33 GHz, it can directly run on more modern GPUs and achieve a greater speedup ratio without being changed. Experiments demonstrate that (1) our CUDA accelerated phase correction achieves exactly the same result with the non-accelerated implementation, (2) the results of our CUDA accelerated rotation correction and translation correction have only slight differences with those of their non-accelerated implementation, (3) images reconstructed from the motion correction results of CUDA accelerated methods proposed in this paper satisfy the clinical requirements, and (4) the speed up ratio is close to 6.5.     

FPGA implementation of real-time SENSE reconstruction using pre-scan and Emaps sensitivities

Sensitivity Encoding (SENSE) is a widely used technique in Parallel Magnetic Resonance Imaging (MRI) to reduce scan time. Reconfigurable hardware based architecture for SENSE can potentially provide image reconstruction with much less computation time. Application specific hardware platform for SENSE may dramatically increase the power efficiency of the system and can decrease the execution time to obtain MR images. A new implementation of SENSE on Field Programmable Gate Array (FPGA) is presented in this study, which provides real-time SENSE reconstruction right on the receiver coil data acquisition system with no need to transfer the raw data to the MRI server, thereby minimizing the transmission noise and memory usage. The proposed SENSE architecture can reconstruct MR images using receiver coil sensitivity maps obtained using pre-scan and eigenvector (E-maps) methods. The results show that the proposed system consumes remarkably less computation time for SENSE reconstruction, i.e., 0.164 ms @ 200 MHz, while maintaining the quality of the reconstructed images with good mean SNR (29 + dB), less RMSE (< 5 × 10^− 2) and comparable artefact power (< 9 × 10^− 4) to conventional SENSE reconstruction. A comparison of the center line profiles of the reconstructed and reference images also indicates a good quality of the reconstructed images. Furthermore, the results indicate that the proposed architectural design can prove to be a significant tool for SENSE reconstruction in modern MRI scanners and its low power consumption feature can be remarkable for portable MRI scanners.     

Chromatic aberration free liquid crystal adaptive optics for flood illuminated retinal camera

A high-resolution, flood-illumination retinal camera using liquid crystal (LC) adaptive optics (AO) is presented. The retinal camera uses light at 780 nm for ocular aberration measurement while light at 655 nm and 593 nm for retinal imaging. In order to avoid chromatic aberrations due to wavelength dependence of LC, we adopt an open-loop technique, in which dynamic correction of aberrations is applied only to the imaging light. A compensation pattern projected on the LC wavefront corrector is adjusted to provide phase wrapping of 2 π for illumination light. We confirmed feasibility of this technique by performing in vivo retinal imaging experiments. Photoreceptors were clearly revealed at both imaging light at 655 nm and 593 nm. Feasibility of the technique was also supported by comparison of the retinal images taken by the present open-loop technique with those taken by the conventional closed-loop one and by analysis of the spatial distribution of the photoreceptors.     

Real-time automatic registration in optical surgical navigation

An image-guided surgical navigation system requires the improvement of the patient-to-image registration time to enhance the convenience of the registration procedure. A critical step in achieving this aim is performing a fully automatic patient-to-image registration. This study reports on a design of custom fiducial markers and the performance of a real-time automatic patient-to-image registration method using these markers on the basis of an optical tracking system for rigid anatomy. The custom fiducial markers are designed to be automatically localized in both patient and image spaces. An automatic localization method is performed by registering a point cloud sampled from the three dimensional (3D) pedestal model surface of a fiducial marker to each pedestal of fiducial markers searched in image space. A head phantom is constructed to estimate the performance of the real-time automatic registration method under four fiducial configurations. The head phantom experimental results demonstrate that the real-time automatic registration method is more convenient, rapid, and accurate than the manual method. The time required for each registration is approximately 0.1 s. The automatic localization method precisely localizes the fiducial markers in image space. The averaged target registration error for the four configurations is approximately 0.7 mm. The automatic registration performance is independent of the positions relative to the tracking system and the movement of the patient during the operation.     

Assessment of cerebral blood perfusion reserve with acetazolamide using 3D spiral ASL MRI: Preliminary experience in pediatric patients

PurposeTo demonstrate the clinical feasibility of a new non-Cartesian cylindrically-distributed spiral 3D pseudo-continuous arterial spin labeling (pCASL) magnetic resonance imaging (MRI) pulse sequence in pediatric patients in quantifying cerebral blood flow (CBF) response to an acetazolamide (ACZ) vasodilator challenge.Materials and methodsMRI exams were performed on two 3 Tesla Philips Ingenia systems using 32 channel head coil arrays. After local institutional review board approval, the 3D spiral-based pCASL technique was added to a standard brain MRI exam and evaluated in 13 pediatric patients (average age: 11.7 ± 6.4 years, range: 1.4–22.2 years). All patients were administered ACZ for clinically indicated reasons. Quantitative whole-brain CBF measurements were computed pre- and post-ACZ to assess cerebrovascular reserve.Results3D spiral pCASL data were successfully reconstructed in all 13 cases. In 11 patients, CBF increased 2.8% to 93.2% after administration of ACZ. In the two remaining patients, CBF decreased by 2.4 to 6.0% after ACZ. The group average change in CBF due to ACZ was approximately 25.0% and individual changes were statistically significant (p < 0.01) in all patients using a paired t-test analysis. CBF perfusion data were diagnostically useful in supporting conventional MR angiography and clinical findings.Conclusion3D cylindrically-distributed spiral pCASL MRI provides a robust approach to assess cerebral blood flow and reserve in pediatric patients.     

Alterations of the optic pathway between unilateral and bilateral optic nerve damage in multiple sclerosis as revealed by the combined use of advanced diffusion kurtosis imaging and visual evoked potentials

ObjectivesWe investigated changes in the optic tract and optic radiation in patients with multiple sclerosis (MS) by comparing unilateral and bilateral optic nerve damage assessed based on visual evoked potentials (VEPs) using advanced diffusion MR metrics.MethodsIn 21 MS patients, diffusion MRI was performed. Maps of fractional anisotropy, apparent diffusion coefficient (ADC), and mean kurtosis (MK) were computed. On the basis of the P100 latency in VEPs, the MS patients were divided into three groups: bilateral (n = 7), unilateral (n = 7), and no abnormality (n = 7). Their optic tracts and optic radiations were analyzed with diffusion MRI-based fiber tracking. We also investigated the correlations between diffusion parameters and VEPs (n = 21).ResultsIn the optic tract, the diffusion changes in each of the three groups showed step-like changes. The diffusion changes in the optic radiations of the unilateral group were similar to those in the normal VEP group. Only the bilateral group showed significantly higher ADC and lower MK relative to the other two groups (P < 0.05, Steel–Dwass multiple-comparison test). A significant positive correlation between VEP latency and ADC and a significant negative correlation between VEP latency and MK were observed (P < 0.01, Spearman's correction).ConclusionsWe first evaluated the relationship between VEPs and DKI and concluded that the lateral geniculate nucleus may compensate for unilateral damage in the pre-geniculate optic pathway via neural plasticity.     

Accuracy and efficiency of an infrared based positioning and tracking system for patient set-up and monitoring in image guided radiotherapy

An infrared based positioning and tracking (IPT) system was introduced and its accuracy and efficiency for patient setup and monitoring were tested for daily radiotherapy treatment. The IPT system consists of a pair of floor mounted infrared stereoscopic cameras, passive infrared markers and tools used for acquiring localization information as well as a custom controlled software which can perform the positioning and tracking functions. The evaluation of IPT system characteristics was conducted based on the AAPM 147 task report. Experiments on spatial drift and reproducibility as well as static and dynamic localization accuracy were carried out to test the efficiency of the IPT system. Measurements of known translational (up to 55.0 mm) set-up errors in three dimensions have been performed on a calibration phantom. The accuracy of positioning was evaluated on an anthropomorphic phantom with five markers attached to the surface; the precision of the tracking ability was investigated through a sinusoidal motion platform. For the monitoring of the respiration, three volunteers contributed to the breathing testing in real time. The spatial drift of the IPT system was 0.65 mm within 60 min to be stable. The reproducibility of position variations were between 0.01 and 0.04 mm. The standard deviation of static marker localization was 0.26 mm. The repositioning accuracy was 0.19 mm, 0.29 mm, and 0.53 mm in the left/right (L/R), superior/inferior (S/I) and anterior/posterior (A/P) directions, respectively. The measured dynamic accuracy was 0.57 mm and discrepancies measured for the respiratory motion tracking was better than 1 mm. The overall positioning accuracy of the IPT system was within 2 mm. In conclusion, the IPT system is an accurate and effective tool for assisting patient positioning in the treatment room. The characteristics of the IPT system can successfully meet the needs for real time external marker tracking and patient positioning as well as respiration monitoring during image guided radiotherapy treatments.     

Iterative reconstruction of radially-sampled 31P bSSFP data using prior information from 1H MRI

The purpose of this study is to improve direct phosphorus (³¹P) MR imaging. Therefore, 3D density-adapted radially-sampled balanced steady-state free precession (bSSFP) sequences were developed and an iterative approach exploiting additional anatomical information from hydrogen (¹H) data was evaluated. Three healthy volunteers were examined at B₀ = 7 T in order to obtain the spatial distribution of the phosphocreatine (PCr) intensities in the human calf muscle with a nominal isotropic resolution of 10 mm in an acquisition time of 10 min. Three different bSSFP gradient schemes were investigated. The highest signal-to-noise ratio (SNR) was obtained for a scheme with two point-reflected density-adapted gradients. Furthermore, the conventional reconstruction based on a gridding algorithm was compared to an iterative method using an ¹H MRI constraint in terms of a segmented binary mask, which comprises prior knowledge. The parameters of the iterative approach were optimized and evaluated by simulations featuring ³¹P MRI parameters. Thereby, partial volume effects as well as Gibbs ringing artifacts could be reduced. In conclusion, the iterative reconstruction of ³¹P bSSFP data using an ¹H MRI constraint is appropriate for investigating regions where sharp tissue boundaries occur and leads to images that represent the real PCr distributions better than conventionally reconstructed images.     

A purpose-built neck coil for black-blood DANTE-prepared carotid artery imaging at 7 T

Atherosclerotic plaques in the bifurcation of the carotid arteries can pose a significant health risk due to possible plaque rupture and subsequent stroke. The assessment of plaques, and evaluation of the risk they pose, can be performed with Black-Blood (BB) vessel wall magnetic resonance imaging. However, resolution at standard clinical field strengths (up to 3 T) is limited, hampering reliable assessment and diagnosis. The aim of this study was to investigate the benefits of 7 T MRI using a BB application that has been successful at clinical field strengths. Therefore, for BB imaging, each sequence was preceded with ‘Delay Alternating with Nutation for Tailored Excitation’ (DANTE) preparation pulses for blood signal suppression. A coil comprising a 4-channel Tx array was designed and built to provide the required excitation coverage for the DANTE train; and a 4-channel Rx array was constructed to target the carotid bifurcation. Human and phantom results showed satisfactory blood suppression and comparable SNR and CNR to 3 T, therefore demonstrating the feasibility of the application at 7 T. However, the imposed SAR restrictions led to long scan times and subsequent motion artifacts. Thus, more accurate local SAR supervision schemes are required which could lead to a further improvement of BB DANTE vessel wall imaging at 7 T.     

Photovoltaic properties of black silicon microstructured by femtosecond pulsed laser

Prototype devices based on black silicon have been fabricated by microstructuring 250 μm thick multicrystalline n doped silicon wafers using femtosecond pulsed laser in ambient gas of SF₆ to measure its photovoltaic properties. The enhanced optical absorption of black silicon extends across the visible region and all the black silicons prepared in this work exhibit enhanced optical absorption close to 90% from 300 nm to 800 nm. The highest open-circuit voltage (V_oc) and short-circuit current (I_sc) under the illumination of He–Ne continuous laser at 632.8 nm were measured to be 53.3 mV and 0.11 mA, respectively at a maximum power conversion efficiency of 1.44%. Upon excitation with He–Ne continuous laser at 632.8 nm, external quantum efficiency (EQE) of black silicon as high as 112.9% has also been observed. Development of black silicon for photovoltaic purposes could open up a new perspective in achieving high efficient silicon-based solar cell by means of the enhanced optical absorption in the visible region. The current–voltage characteristic and photo responsivity of these prototype devices fabricated with microstructured silicon were also investigated.     

Towards intrinsic R2* imaging in the prostate at 3 and 7 tesla

PurposeHypoxia is an important marker for resistance to therapy. In this study, we quantify the macroscopic effects of R2* mapping in prostate cancer patients incorporating susceptibility matching and field strengths effects.Materials and methods91 patients were scanned without endorectal coil (ERC) at 3 T. Only when rectal gas was absent, data was included for analysis. Another group of 10 patients was scanned using a susceptibility matched ERC. To assess the residual contamination of R2 and macroscopic field non-uniformities, a group of 10 patients underwent ultra-high resolution 7 T MRI.ResultsOf the patients scanned at 3 T 60% presented rectal gas and were excluded, due to susceptibility artifacts. At 3 T the tumor was significantly different (P < 0.01) from the healthy surrounding tissue in R2* values at intrapatient level. Using the measured median R2* value of 24.9 s^− 1 at 3 T and 43.2 s^− 1 at 7 T of the peripheral zone, the minimum contribution of macroscopic susceptibility effects is 15% at 3 T.ConclusionR2* imaging might be a promising tool for hypoxia imaging, particularly when minimizing macroscopic susceptibility effects contaminating intrinsic R2* of tissue, such as rectal gas. At 3 T macroscopic effects still contribute 15% in the R2* value, compared to ultra-high resolution R2* mapping at 7 T.     

Influence and correction of temperature on optical measurement for fat and protein contents in a complex food model system

The influence of temperature on optical measurements has been studied for determining fat and protein contents in complex food systems. A model system consisting of mixtures of fat, protein, water and emulsion was developed to create an imitation of complex food systems. The changes in optical properties of the system from 25 °C to 40 °C were measured in the wavelength from 1100 nm to 1670 nm. Irregular changes as a whole were founded and therefore a statistical method was needed to correct the temperature effect. A method called global robust temperature calibration model is proposed and the correction effect was validated. The results indicated that it can significantly reduce the temperature effect on optical measurement.     

Optimization of the energy response of radiographic films by Monte Carlo method

In the present work a simple model for calculation of the energy response of radiographic films was introduced. According to the model the energy response of a radiographic film is directly proportional to the optical density on the film and thus to the number of developed grains in the emulsion. The model was simulated by Monte Carlo method using MCNP code and the relative energy response of Kodak type 2 film under a few filters of A.E.R.E./R.P.S. film badge was calculated. The simulated responses were in agreement with the experimental data in the region of 30 keV–1.5 MeV. In the next stage a multi-element filter was simulated to optimize the energy response in the above energies. The energy response varied by 25% between 40 keV and 1.5 MeV. So the dose received by the film is equivalent to the desired true dose and there would be no need to the correction factors.     

Comparison of 7 T and 3 T MRI in patients with moyamoya disease

Magnetic resonance imaging and magnetic resonance angiography (MRI/MRA) are widely used for evaluating the moyamoya disease (MMD). This study compared the diagnostic accuracy of 7 Tesla (T) and 3 T MRI/MRA in MMD. In this case control study, 12 patients [median age: 34 years; range (10–66 years)] with MMD and 12 healthy controls [median age: 25 years; range (22–59 years)] underwent both 7 T and 3 T MRI/MRA. To evaluate the accuracy of MRI/MRA in MMD, five criteria were compared between imaging systems of 7 T and 3 T: Suzuki grading system, internal carotid artery (ICA) diameter, ivy sign, flow void of the basal ganglia on T2-weighted images, and high signal intensity areas of the basal ganglia on time-of-flight (TOF) source images. No difference was observed between 7 T and 3 T MRI/MRA in Suzuki stage, ICA diameter, and ivy sign score; while, 7 T MRI/MRA showed a higher detection rate in the flow void on T2-weighted images and TOF source images (p < 0.001). Receiver operating characteristic curves of both T2 and TOF criteria showed that 7 T MRI/MRA had higher sensitivity and specificity than 3 T MRI/MRA. Our findings indicate that 7 T MRI/MRA is superior to 3 T MRI/MRA for the diagnosis of MMD in point of detecting the flow void in basal ganglia by T2-weighted and TOF images.     

Preparation of an Au-Pt alloy free from artifacts in magnetic resonance imaging

PurposeWhen magnetic resonance imaging (MRI) is performed on patients carrying metallic implants, artifacts can disturb the images around the implants, often making it difficult to interpret them appropriately. However, metallic materials are and will be indispensable as raw materials for medical devices because of their electric conductivity, visibility under X-ray fluoroscopy, and other favorable features. What is now desired is to develop a metallic material which causes no artifacts during MRI.Materials and methodsIn the present study, we prepared a single-phase and homogeneous Au-Pt alloys (Au; diamagnetic metal, and Pt; paramagnetic metal) by the processing of thermal treatment. Volume magnetic susceptibility was measured with a SQUID Flux Meter and MRI artifact was evaluated using a 1.5-T scanner.ResultsAfter final thermal treatment, an entirely recrystallized homogeneous organization was noted. The Au-35Pt alloy was shown to have a volume magnetic susceptibility of − 8.8 ppm, causing almost free from artifacts during MRI.ConclusionsWe thus prepared an Au-35Pt alloy which had a magnetic susceptibility very close to that of living tissue and caused much fewer artifacts during MRI. It is promising as a material for spinal cages, intracranial electrodes, cerebral aneurysm embolization coils, markers for MRI and so on.     

Thermal performance analysis of vacuum variable-temperature blackbody system

In this paper, the design and structure of a vacuum variable-temperature blackbody system were described, and the steady-state thermal analysis of a 3-D blackbody model was presented. Also, the thermal performance of the blackbody was evaluated using an infrared camera system. The blackbody system was constructed to operate under vacuum conditions (2.67 × 10⁻² Pa) to reduce its temperature uncertainty, which can be caused by vapor condensation at low temperatures usually below 273.15 K. A heat sink and heat shield including a cold shield were embedded around the radiator to maintain the heat balance of the blackbody. A simplified 3-D model of the blackbody including a radiator, heat sink, heat shield, cold shield, and heat source was thermophysically evaluated by performing finite elements analysis using the extended Stefan–Boltzmann’s rule, and the infrared radiating performance of the developed system was analyzed using an infrared camera system. On the basis of the results of measurements and simulations, we expect that the suggested blackbody system can serve as a highly stable reference source for the calibration and measurement of infrared optical systems within operational temperature ranges.     

Modeling of laser-diode side-pumped continuous wave Nd:YAG/KTP green laser

By reasonably assuming the distribution of the pump light in laser medium, the equations about the absorption and the gain for the end-pumped lasers are applied to the side-pumped ones, and a theoretical model for laser-diode side-pumped continuous wave intracavity-frequency-doubling lasers is given, in which the thermally induced diffraction loss and the variation of the fundamental wave radii with the pump power are considered. By using a Z-type cavity, a laser-diode side-pumped continuous wave Nd:YAG/KTP green laser is realized. The threshold pump power is 15 W, and the highest output power of the green laser is 3.75 W at the pump power of 160 W, corresponding to an optical–optical slope efficiency 2.6%. The experimental results are consistent with the numerical solutions.     

LD-end-pumped passively Q-switched Nd:YAG ceramic laser with single wall carbon nanotube saturable absorber

We report on a LD-end-pumped passively Q-switched Nd:YAG ceramic laser by using a novel single wall carbon nanotube saturable absorber (SWCNT-SA). The SWCNT wafer was fabricated by electric Arc discharge method on quartz substrate with absorption wavelength of 1064 nm. We firstly investigated the continuous wave (CW) laser performance and scattering properties of Nd:YAG ceramic sample. For the case of passively Q-switched operation, a maximum output power of 376 mW was obtained at an incident pump power of 8.68 W at 808 nm, corresponding to an optical–optical conversion efficiency of 4.3%. The repetition rate as the increase of pump power varied from 14 to 95 kHz. The minimum pulse duration of 1.2 μs and maximum pulse energy of 4.5 μJ was generated at a repetition rate of 31.8 kHz.     

Multispectral 3D phase-encoded turbo spin-echo for imaging near metal: Limitations and possibilities demonstrated by simulations and phantom experiments

To see improvements in the imaging performance near biomaterial implants we assessed a multispectral fully phase-encoded turbo spin-echo (ms3D-PE-TSE) sequence for artifact reduction capabilities and scan time efficiency in simulation and phantom experiments.For this purpose, ms3D-PE-TSE and ms3D-TSE sequences were implemented to obtain multispectral images (± 20 kHz) of a cobalt-chromium (CoCr) knee implant embedded in agarose. In addition, a knee implant computer model and the acquired ms3D-PE-TSE images were used to investigate the possibilities for scan time acceleration using field-of-view (FOV) reduction for off-resonance frequency bins and compressed sensing reconstructions of undersampled data. Both acceleration methods were combined to acquire a + 10 kHz frequency bin in a second experiment.The obtained ms3D-PE-TSE images showed no susceptibility related artifacts, while ms3D-TSE images suffered from hyper-intensity artifacts. The limitations of ms3D-TSE were apparent in the far off-resonance regions (±[10–20] kHz) located close to the implant. The scan time calculations showed that ms3D-PE-TSE can be applied in a clinically relevant timeframe (~ 12 min), when omitting the three central frequency bins. The feasibility of CS acceleration for ms3D-PE-TSE was demonstrated using retrospective reconstructions before combining CS and rFOV imaging to decrease the scan time for the + 10 kHz frequency bin from ~ 10.9 min to ~ 3.5 min, while also increasing the spatial resolution fourfold. The temporally resolved signal of ms3D-PE-TSE proved to be useful to decrease the intensity ripples after sum-of-squares reconstructions and increase the signal-to-noise ratio.The presented results suggest that the scan time limitations of ms3D-PE-TSE can be sufficiently addressed when focusing on signal acquisitions in the direct vicinity of metal implants. Because these regions cannot be measured with existing multispectral methods, the presented ms3D-PE-TSE method may enable the detection of inflammation or (pseudo-)tumors in locations close to the implant.     

Accelerating MRI fat quantification using a signal model-based dictionary to assess gastric fat volume and distribution of fat fraction

To quantify intragastric fat volume and distribution with accelerated magnetic resonance (MR) imaging using signal model-based dictionaries (DICT) in comparison to conventional parallel imaging (CG-SENSE). This study was approved by the local ethics committee and written informed consent was obtained. Seven healthy subjects were imaged after intake of a lipid emulsion and data at three different time points during the gastric emptying process was acquired in order to cover a range of fat fractions. Fully sampled and prospectively undersampled image data at a reduction factor of 4 were acquired using a multi gradient echo sequence at 1.5T. Retrospectively and prospectively undersampled data were reconstructed with DICT and CG-SENSE. Image quality of the retrospectively undersampled data was assessed relative to the fully sampled reference using the root mean square error (RMSE). In order to assess the agreement of fat volumes and intragastric fat distribution, Bland-Altman analysis and linear regression were performed on the data. The RMSE in intragastric content (ΔRMSE = 0.10 ± 0.01, P < 0.001) decreased significantly with DICT relative to CG-SENSE. CG-SENSE overestimated fat volumes (bias 2.1 ± 1.3 mL; confidence limits 5.4 and − 1.1 mL) in comparison to the prospective DICT reconstruction (bias − 0.1 ± 0.7 mL; confidence limits 1.8 and − 2.0 mL). There was a good agreement in fat distribution between the images reconstructed by retrospective DICT and the reference images (regression slope: 1.01, R² = 0.961). Accelerating gastric MRI by integrating a dictionary-based signal model allows for improved image quality and increases accuracy of fat quantification during breathholds.