Snapshot imaging Mueller matrix polarimeter using modified Savart polariscopes.

We present a snapshot imaging Mueller matrix polarimeter using modified Savart polariscopes (MSP-SIMMP). The MSP-SIMMP contains both the polarizing optics and the analyzing optics encoding all Mueller matrix components of the sample into the interferogram by the spatial modulation technique. An interference model and the methods of reconstruction and calibration are discussed. To demonstrate the feasibility of the proposed MSP-SIMMP, the numerical simulation and the laboratory experiment of a design example are presented. The remarkable advantage of the MSP-SIMMP is easy to calibrate. Moreover, compared with conventional imaging Mueller matrix polarimeters with rotating parts, the advantage of the proposed instrument is simple, compact, snapshot-enabled, and stationary (no moving parts).

Calibration of spatially modulated snapshot imaging polarimeter based on phase-shift interferometry.

The snapshot imaging polarimeters (SIPs) using spatial modulation have gained increasing popularity due to their capability of obtaining all four Stokes parameters in a single measurement. However, the existing reference beam calibration techniques cannot extract the modulation phase factors of the spatially modulated system. In this paper, a calibration technique based on a phase-shift interference (PSI) theory is proposed to address this issue. The proposed technique can accurately extract and demodulate the modulation phase factors through measuring the reference object at different polarization analyzer orientations and performing a PSI algorithm. Using the snapshot imaging polarimeter with modified Savart polariscopes as an example, the basic principle of the proposed technique is analyzed in detail. Subsequently, the feasibility of this calibration technique was demonstrated by a numerical simulation and a laboratory experiment. This work provides a different perspective for the calibration of a spatially modulated snapshot imaging polarimeter.

On the generalizability of diffusion MRI signal representations across acquisition parameters, sequences and tissue types: Chronicles of the MEMENTO challenge.

Diffusion MRI (dMRI) has become an invaluable tool to assess the microstructural organization of brain tissue. Depending on the specific acquisition settings, the dMRI signal encodes specific properties of the underlying diffusion process. In the last two decades, several signal representations have been proposed to fit the dMRI signal and decode such properties. Most methods, however, are tested and developed on a limited amount of data, and their applicability to other acquisition schemes remains unknown. With this work, we aimed to shed light on the generalizability of existing dMRI signal representations to different diffusion encoding parameters and brain tissue types. To this end, we organized a community challenge - named MEMENTO, making available the same datasets for fair comparisons across algorithms and techniques. We considered two state-of-the-art diffusion datasets, including single-diffusion-encoding (SDE) spin-echo data from a human brain with over 3820 unique diffusion weightings (the MASSIVE dataset), and double (oscillating) diffusion encoding data (DDE/DODE) of a mouse brain including over 2520 unique data points. A subset of the data sampled in 5 different voxels was openly distributed, and the challenge participants were asked to predict the remaining part of the data. After one year, eight participant teams submitted a total of 80 signal fits. For each submission, we evaluated the mean squared error, the variance of the prediction error and the Bayesian information criteria. The received submissions predicted either multi-shell SDE data (37%) or DODE data (22%), followed by cartesian SDE data (19%) and DDE (18%). Most submissions predicted the signals measured with SDE remarkably well, with the exception of low and very strong diffusion weightings. The prediction of DDE and DODE data seemed more challenging, likely because none of the submissions explicitly accounted for diffusion time and frequency. Next to the choice of the model, decisions on fit procedure and hyperparameters play a major role in the prediction performance, highlighting the importance of optimizing and reporting such choices. This work is a community effort to highlight strength and limitations of the field at representing dMRI acquired with trending encoding schemes, gaining insights into how different models generalize to different tissue types and fiber configurations over a large range of diffusion encodings.

Acquisition of a full-resolution image and aliasing reduction for a spatially modulated imaging polarimeter with two snapshots.

A snapshot imaging polarimeter using spatial modulation can encode four Stokes parameters allowing instantaneous polarization measurement from a single interferogram. However, the reconstructed polarization images could suffer a severe aliasing signal if the high-frequency component of the intensity image is prominent and occurs in the polarization channels, and the reconstructed intensity image also suffers reduction of spatial resolution due to low-pass filtering. In this work, a method using two anti-phase snapshots is proposed to address the two problems simultaneously. The full-resolution target image and the pure interference fringes can be obtained from the sum and the difference of the two anti-phase interferograms, respectively. The polarization information reconstructed from the pure interference fringes does not contain the aliasing signal from the high-frequency component of the object intensity image. The principles of the method are derived and its feasibility is tested by both computer simulation and a verification experiment. This work provides a novel method for spatially modulated imaging polarization technology with two snapshots to simultaneously reconstruct a full-resolution object intensity image and high-quality polarization components.

Spatiotemporal denoising of low-dose cardiac CT image sequences using RecycleGAN.

Electrocardiogram (ECG)-gated multi-phase computed tomography angiography (MP-CTA) is frequently used for diagnosis of coronary artery disease. Radiation dose may become a potential concern as the scan needs to cover a wide range of cardiac phases during a heart cycle. A common method to reduce radiation is to limit the full-dose acquisition to a predefined range of phases while reducing the radiation dose for the rest. Our goal in this study is to develop a spatiotemporal deep learning method to enhance the quality of low-dose CTA images at phases acquired at reduced radiation dose. Recently, we demonstrated that a deep learning method, Cycle-Consistent generative adversarial networks (CycleGAN), could effectively denoise low-dose CT images through spatial image translation without labeled image pairs in both low-dose and full-dose image domains. As CycleGAN does not utilize the temporal information in its denoising mechanism, we propose to use RecycleGAN, which could translate a series of images ordered in time from the low-dose domain to the full-dose domain through an additional recurrent network. To evaluate RecycleGAN, we use the XCAT phantom program, a highly realistic simulation tool based on real patient data, to generate MP-CTA image sequences for 18 patients (14 for training, 2 for validation and 2 for test). Our simulation results show that RecycleGAN can achieve better denoising performance than CycleGAN based on both visual inspection and quantitative metrics. We further demonstrate the superior denoising performance of RecycleGAN using clinical MP-CTA images from 50 patients.

Optimizing the performance of local canonical correlation analysis in fMRI using spatial constraints.

The benefits of locally adaptive statistical methods for fMRI research have been shown in recent years, as these methods are more proficient in detecting brain activations in a noisy environment. One such method is local canonical correlation analysis (CCA), which investigates a group of neighboring voxels instead of looking at the single voxel time course. The value of a suitable test statistic is used as a measure of activation. It is customary to assign the value to the center voxel for convenience. The method without constraints is prone to artifacts, especially in a region of localized strong activation. To compensate for these deficiencies, the impact of different spatial constraints in CCA on sensitivity and specificity are investigated. The ability of constrained CCA (cCCA) to detect activation patterns in an episodic memory task has been studied. This research shows how any arbitrary contrast of interest can be analyzed by cCCA and how accurate P-values optimized for the contrast of interest can be computed using nonparametric methods. Results indicate an increase of up to 20% in detecting activation patterns for some of the advanced cCCA methods, as measured by ROC curves derived from simulated and real fMRI data.

Prediction of dMRI signals with neural architecture search.

BACKGROUND: There is growing interest in the neuroscience community in estimating and mapping microscopic properties of brain tissue non-invasively using magnetic resonance measurements. Machine learning methods are actively investigated to predict the signals measured in diffusion magnetic resonance imaging (dMRI). NEW METHOD: We applied the neural architecture search (NAS) to train a recurrent neural network to generate a multilayer perceptron to predict the dMRI data of unknown signals based on the different acquisition parameters and training data. The search space of NAS is the number of neurons in each layer of the multilayer perceptron network. To our best knowledge, this is the first time to apply NAS to solve the dMRI signal prediction problem. RESULTS: The experimental results demonstrate that the proposed NAS method can achieve fast training and predict dMRI signals accurately. For dMRI signals with four acquisition strategies of double diffusion encoding (DDE), double oscillating diffusion encoding (DODE), multi-shell and DSI-like pulsed gradient spin-echo (PGSE), the mean squared errors of the multilayer perceptron network designed by NAS are 0.0043, 0.0034, 0.0147 and 0.0199, respectively. COMPARISON WITH EXISTING METHOD(S): We also compared NAS with other machine learning prediction methods, such as support vector regression (SVR), decision tree (DT) and random forest (RF), k-nearest neighbors (KNN), adaboost regressor (AR), gradient boosting regressor (GBR) and extra-trees regressor (ET). NAS achieved the better prediction performance in most cases. CONCLUSION: In this study, NAS was developed for the prediction of dMRI signals and could become an effective prediction tool.

Texture transformer super-resolution for low-dose computed tomography.

Computed tomography (CT) is widely used to diagnose many diseases. Low-dose CT has been actively pursued to lower the ionizing radiation risk. A relatively smoother kernel is typically used in low-dose CT to suppress image noise, which may sacrifice spatial resolution. In this work, we propose a texture transformer network to simultaneously reduce image noise and improve spatial resolution in CT images. This network, referred to as Texture Transformer for Super Resolution (TTSR), is a reference-based deep-learning image super-resolution method built upon a generative adversarial network (GAN). The noisy low-resolution CT (LRCT) image and the routine-dose high-resolution (HRCT) image are severed as the query and key in a transformer, respectively. Image translation is optimized through deep neural network (DNN) texture extraction, correlation embedding, and attention-based texture transfer and synthesis to achieve joint feature learning between LRCT and HRCT images for super-resolution CT (SRCT) images. To evaluate SRCT performance, we use the data from both simulations of the XCAT phantom program and the real patient data. Peak signal-to-noise ratio (PSNR), structural similarity index measure (SSIM), and feature similarity (FSIM) index are used as quantitative metrics. For comparison of SRCT performance, the cubic spline interpolation, SRGAN (a GAN super-resolution with an additional content loss), and GAN-CIRCLE (a GAN super-resolution with cycle consistency) were used. Compared to the other two methods, TTSR can restore more details in SRCT images and achieve better PSNR, SSIM, and FSIM for both simulation and real-patient data. In addition, we show that TTSR can yield better image quality and demand much less computation time than high-resolution low-dose CT images denoised by block-matching and 3D filtering (BM3D) and GAN-CIRCLE. In summary, the proposed TTSR method based on texture transformer and attention mechanism provides an effective and efficient tool to improve spatial resolution and suppress noise of low-dose CT images.

Effects of motion, attenuation, and scatter corrections on gated cardiac SPECT reconstruction.

PURPOSE: In gated cardiac single photon emission computed tomography (SPECT), image reconstruction is often hampered by various degrading factors including depth-dependent spatial blurring, attenuation, scatter, motion blurring, and low data counts. Consequently, there has been significant development in image reconstruction methods for improving the quality of reconstructed images. The goal of this work is to investigate how these degrading factors will impact the reconstructed myocardium when different reconstruction methods are used. METHODS: The authors conduct a comparative study of the effects of these degrading factors on the accuracy of myocardium by several reconstruction algorithms, including (1) a clinical spatiotemporal processing method, (2) maximum likelihood (ML) estimation, (3) 3D maximum a posteriori (MAP) estimation, (4) 3D MAP with posttemporal filtering, and (5) motion-compensated spatiotemporal (4D) reconstruction. To quantify the reconstruction results, the authors use the following measures on different aspects of the myocardium: (1) overall error level in the myocardium, (2) regional accuracy of the left ventricle (LV) wall, (3) uniformity of the LV, (4) accuracy of regional time activity curves by normalized cross-correlation coefficient, and (5) perfusion defect detectability. The authors also assess the effectiveness of degrading corrections in reconstruction by considering an upper bound for each reconstruction method, which represents what would be achieved by each method if the acquired data were free from attenuation and scatter degradations. In the experiments the authors use Monte Carlo simulated cardiac gated SPECT imaging based on the 4D NURBS-based cardiac-torso (NCAT) phantom with different patient geometry and lesion settings, in which the simulated ground truth is known for the purpose of quantitative evaluation. RESULTS: The results demonstrate that use of temporal processing in reconstruction (Methods 1, 4, and 5 above) can greatly improve the reconstructed myocardium in terms of both error level and perfusion defect detection. In low-count gated studies, it can have even greater impact than other degrading factors. Both attenuation and scatter corrections can lead to reduced error levels in the myocardium in all methods; in particular, with 4D the bias can be reduced by as much as four-fold compared to no correction. There is a slight increase in noise level observed with scatter correction. A significant improvement in heart wall appearance is demonstrated in reconstruction results from three sets of clinical acquisitions as correction for degradations is combined with refinement of temporal filtering. CONCLUSIONS: Correction for degrading factors such as resolution, attenuation, scatter, and motion blur can all lead to improved image quality in cardiac gated SPECT reconstruction. However, their effectiveness could also vary with the reconstruction algorithms used. Both attenuation and scatter corrections can effectively reduce the bias level of the reconstructed LV wall, though scatter correction is also observed to increase the variance level. Use of temporal processing in reconstruction can have greater impact on the accuracy of the myocardium than correction of other degrading factors. Overall, use of degrading corrections in 4D reconstruction is shown to be most effective for improving both reconstruction accuracy of the myocardium and detectability of perfusion defects in gated images.

General simultaneous motion estimation and image reconstruction (G-SMEIR).

To achieve better performance for 4D multi-frame reconstruction with the parametric motion model (MF-PMM), a general simultaneous motion estimation and image reconstruction (G-SMEIR) method is proposed. In G-SMEIR, projection domain motion estimation and image domain motion estimation are performed alternatively to achieve better 4D reconstruction. This method can mitigate the local optimum trapping problem in either domain. To improve computational efficiency, the image domain motion estimation is accelerated by adapting fast convergent algorithms and graphics processing unit (GPU) computing. The proposed G-SMEIR method is tested using a cone-beam computed tomography (CBCT) simulation study of 4D XCAT phantom at different dose levels and compared with 3D total variation-based reconstruction (3D TV), 4D reconstruction with image domain motion estimation (IM4D), and SMEIR. G-SMEIR shows strong denoising capability and achieves similar performance at regular dose and half dose. The root mean squared error (RMSE) of G-SMEIR is the best among the four methods and improved about 12% over SMEIR for all respiratory phase images at full dose. G-SMEIR also achieved the best structural similarity index (SSIM) values among all methods. More importantly, G-SMEIR leads to more than 40% improvement of the mean deviation from the phantom tumor motion over SMEIR. A preliminary patient CBCT image reconstruction also shows better image quality of G-SMEIR than that of the frame-by-frame reconstruction (3D TV) and MF-PMM either using image domain motion estimation (IM4D) or using projection domain motion estimation (SMEIR) alone. G-SMEIR with a flexible combination of image domain and projection domain motion estimation provides an effective tool for 4D tomographic reconstruction.

Investigation of Low-Dose CT Image Denoising Using Unpaired Deep Learning Methods.

Low-dose computed tomography (LDCT) is desired due to prevalence and ionizing radiation of CT, but suffers elevated noise. To improve LDCT image quality, an image-domain denoising method based on cycle-consistent generative adversarial network ("CycleGAN") is developed and compared with two other variants, IdentityGAN and GAN-CIRCLE. Different from supervised deep learning methods, these unpaired methods can effectively learn image translation from the low-dose domain to the full-dose (FD) domain without the need of aligning FDCT and LDCT images. The results on real and synthetic patient CT data show that these methods can achieve peak signal-to-noise ratio (PSNR) and structural similarity index (SSIM) comparable to, if not better than, the other state-of-the-art denoising methods. Among CycleGAN, IdentityGAN, and GAN-CIRCLE, the later achieves the best denoising performance with the shortest computation time. Subsequently, GAN-CIRCLE is used to demonstrate that the increasing number of training patches and of training patients can improve denoising performance. Finally, two non-overlapping experiments, i.e. no counterparts of FDCT and LDCT images in the training data, further demonstrate the effectiveness of unpaired learning methods. This work paves the way for applying unpaired deep learning methods to enhance LDCT images without requiring aligned full-dose and low-dose images from the same patient.

H2RSPET: a 0.5 mm resolution high-sensitivity small-animal PET scanner, a simulation study.

With the goal of developing a total-body small-animal PET system with a high spatial resolution of ∼0.5 mm and a high sensitivity >10% for mouse/rat studies, we simulated four scanners using the graphical processing unit-based Monte Carlo simulation package (gPET) and compared their performance in terms of spatial resolution and sensitivity. We also investigated the effect of depth-of-interaction (DOI) resolution on the spatial resolution. All the scanners are built upon 128 DOI encoding dual-ended readout detectors with lutetium yttrium oxyorthosilicate (LYSO) arrays arranged in 8 detector rings. The solid angle coverages of the four scanners are all ∼0.85 steradians. Each LYSO element has a cross-section of 0.44 × 0.44 mm2 and the pitch size of the LYSO arrays are all 0.5 mm. The four scanners can be divided into two groups: (1) H2RS110-C10 and H2RS110-C20 with 40 × 40 LYSO arrays, a ring diameter of 110 mm and axial length of 167 mm, and (2) H2RS160-C10 and H2RS160-C20 with 60 × 60 LYSO arrays, a diameter of 160 mm and axial length of 254 mm. C10 and C20 denote the crystal thickness of 10 and 20 mm, respectively. The simulation results show that all scanners have a spatial resolution better than 0.5 mm at the center of the field-of-view (FOV). The radial resolution strongly depends on the DOI resolution and radial offset, but not the axial resolution and tangential resolution. Comparing the C10 and C20 designs, the former provides better resolution, especially at positions away from the center of the FOV, whereas the latter has 2× higher sensitivity (∼10% versus ∼20%). This simulation study provides evidence that the 110 mm systems are a good choice for total-body mouse studies at a lower cost, whereas the 160 mm systems are suited for both total-body mouse and rat studies.

Single source precursor synthesized CuS nanoparticles for NIR phototherapy of cancer and photodegradation of organic carcinogen.

Herein, we report cost effective and body compatible CuS nanoparticles (NPs) derived from a single source precursor as photothermal agent for healing deep cancer and photocatalytic remediation of organic carcinogens. These NPs efficiently kill MCF7 cells (both in vivo and in vitro) under NIR irradiation by raising the temperature of tumor cells. Such materials can be used for the treatment of deep cancer as they can produce a heating effect using high wavelength and deeply penetrating NIR radiation. Furthermore, CuS NPs under solar light irradiation efficiently convert p-nitrophenol (PNP), an environmental carcinogen, to p-aminophenol (PAP) of pharmaceutical implication. In a nutshell, CuS can be used for the treatment of deep cancer and for the remediation of carcinogenic pollutants. There seems an intrinsic connection between the two functions of CuS NPs that need to be explored in length.

Regularized Fully 5D Reconstruction of Cardiac Gated Dynamic SPECT Images.

In our recent work, we proposed an image reconstruction procedure aimed to unify gated imaging and dynamic imaging in nuclear cardiac imaging. With this procedure the goal is to obtain an image sequence from a single acquisition which shows simultaneously both cardiac motion and tracer distribution change over the course of imaging. In this work, we further develop and demonstrate this procedure for fully 5D (3D space plus time plus gate) reconstruction in gated, dynamic cardiac SPECT imaging, where the challenge is even greater without the use of multiple fast camera rotations. For 5D reconstruction, we develop and compare two iterative algorithms: one is based on the modified block sequential regularized EM (BSREM-II) algorithm, and the other is based on the one-step late (OSL) algorithm. In our experiments, we simulated gated cardiac imaging with the NURBS-based cardiac-torso (NCAT) phantom and Tc99m-Teboroxime as the imaging agent, where acquisition with the equivalent of only three full camera rotations was used during the course of a 12-minute postinjection period. We conducted a thorough evaluation of the reconstruction results using a number of quantitative measures. Our results demonstrate that the 5D reconstruction procedure can yield gated dynamic images which show quantitative information for both perfusion defect detection and cardiac motion.

Carbon-11 and Carbon-12 beam range verifications through prompt gamma and annihilation gamma measurements: Monte Carlo simulations.

Range uncertainty remains a big concern in particle therapy, as it may cause target dose degradation and normal tissue overdosing. Positron emission tomography (PET) and prompt gamma imaging (PGI) are two promising modalities for range verification. However, the relatively long acquisition time of PET and the relatively low yield of PGI pose challenges for real-time range verification. In this paper, we explore using the primary Carbon-11 (C-11) ion beams to enhance the gamma yield compared to the primary C-12 ion beams to improve PET and PGI by using Monte Carlo simulations of water and PMMA phantoms at four incident energies (95, 200, 300, and 430 MeV u-1). Prompt gammas (PGs) and annihilation gammas (AGs) were recorded for post-processing to mimic PGI and PET imaging, respectively. We used both time-of-flight (TOF) and energy selections for PGI, which boosted the ratio of PGs to background neutrons to 2.44, up from 0.87 without the selections. At the lowest incident energy (100 MeVu-1), PG yield from C-11 was 0.82 times of that from C-12, while AG yield from C-11 was 6 ∼ 11 folds higher than from C-12 in PMMA. At higher energies, PG differences between C-11 and C-12 were much smaller, while AG yield from C-11 was 30%∼90% higher than from C-12 using minute-acquisition. With minute-acquisition, the AG depth distribution of C-11 showed a sharp peak coincident with the Bragg peak due to the decay of the primary C-11 ions, but that of C-12 had no such one. The high AG yield and distinct peaks could lead to more precise range verification of C-11 than C-12. These results demonstrate that using C-11 ion beams for potentially combined PGI and PET has great potential to improve online single-spot range verification accuracy and precision.

A quantitative evaluation study of four-dimensional gated cardiac SPECT reconstruction.

In practice, gated cardiac SPECT images suffer from a number of degrading factors, including distance-dependent blur, attenuation, scatter and increased noise due to gating. Recently, we proposed a motion-compensated approach for four-dimensional (4D) reconstruction for gated cardiac SPECT and demonstrated that use of motion-compensated temporal smoothing could be effective for suppressing the increased noise due to lowered counts in individual gates. In this work, we further develop this motion-compensated 4D approach by also taking into account attenuation and scatter in the reconstruction process, which are two major degrading factors in SPECT data. In our experiments, we conducted a thorough quantitative evaluation of the proposed 4D method using Monte Carlo simulated SPECT imaging based on the 4D NURBS-based cardiac-torso (NCAT) phantom. In particular, we evaluated the accuracy of the reconstructed left ventricular myocardium using a number of quantitative measures including regional bias-variance analyses and wall intensity uniformity. The quantitative results demonstrate that use of motion-compensated 4D reconstruction can improve the accuracy of the reconstructed myocardium, which in turn can improve the detectability of perfusion defects. Moreover, our results reveal that while traditional spatial smoothing could be beneficial, its merit would become diminished with the use of motion-compensated temporal regularization. As a preliminary demonstration, we also tested our 4D approach on patient data. The reconstructed images from both simulated and patient data demonstrated that our 4D method can improve the definition of the LV wall.

gPET: a GPU-based, accurate and efficient Monte Carlo simulation tool for PET.

Monte Carlo (MC) simulation method plays an essential role in the refinement and development of positron emission tomography (PET) systems. However, most existing MC simulation packages suffer from long execution time for practical PET simulations. To fully address this issue, we developed and validated gPET, a graphics processing unit (GPU)-based MC simulation tool for PET. gPET was built on the NVidia CUDA platform. The simulation process was modularized into three functional parts and carried out by the GPU parallel threads: (1) source management, including positron decay, transport and annihilation; (2) gamma transport inside the phantom; and (3) signal detection and processing inside the detector. A hybrid of voxelized (for patient phantoms) and parametrized (for detectors) geometries were employed to sufficiently support particle navigations. Multiple inputs and outputs were available. Hence, a user can flexibly examine different aspects of a PET simulation. We evaluated the performance of gPET in three test cases with benchmark work from GATE8.0, in terms of the testing of the functional modules, the physics models used for gamma transport inside the detector, and the geometric configuration of an irregularly shaped PET detector. Both accuracy and efficiency were quantified. In all test cases, the differences between gPET and GATE for the coincidences with respect to the energy and crystal index distributions are below 3.18% and 2.54%, respectively. The speedup factor is 500 for gPET on a single Titan Xp GPU (1.58 GHz) over GATE8.0 on a single core of Intel i7-6850K CPU (3.6 GHz) for all test cases. In summary, gPET is an accurate and efficient MC simulation tool for PET.

Predication of different stages of Alzheimer's disease using neighborhood component analysis and ensemble decision tree.

BACKGROUND: There is a spectrum of the progression from healthy control (HC) to mild cognitive impairment (MCI) without conversion to Alzheimer's disease (AD), to MCI with conversion to AD (cMCI), and to AD. This study aims to predict the different disease stages using brain structural information provided by magnetic resonance imaging (MRI) data. NEW METHOD: The neighborhood component analysis (NCA) is applied to select most powerful features for prediction. The ensemble decision tree classifier is built to predict which group the subject belongs to. The best features and model parameters are determined by cross validation of the training data. RESULTS: Our results show that 16 out of a total of 429 features were selected by NCA using 240 training subjects, including MMSE score and structural measures in memory-related regions. The boosting tree model with NCA features can achieve prediction accuracy of 56.25% on 160 test subjects. COMPARISON WITH EXISTING METHOD(S): Principal component analysis (PCA) and sequential feature selection (SFS) are used for feature selection, while support vector machine (SVM) is used for classification. The boosting tree model with NCA features outperforms all other combinations of feature selection and classification methods. CONCLUSIONS: The results suggest that NCA be a better feature selection strategy than PCA and SFS for the data used in this study. Ensemble tree classifier with boosting is more powerful than SVM to predict the subject group. However, more advanced feature selection and classification methods or additional measures besides structural MRI may be needed to improve the prediction performance.

Thermal Image Enhancement through the Deconvolution Methods for Low-Cost Infrared Cameras.

An infrared (IR) thermal camera may provide a tool for real-time temperature monitoring for precise disease treatment using heat generated by light-induced photosensitisers, i.e. photothermal/ablation therapies. In this work, we quantitatively demonstrated that the spatial resolution of a low-cost low-resolution IR camera could be improved via two deconvolution methods. The camera point spread function (PSF) was modeled experimentally and used to develop the deconvolution methods: 1) Richardson-Lucy blind deconvolution (BD); and 2) total variation constrained deconvolution (TD). The experimental results showed the improved spatial resolution (at 50% modulation transfer function (MTF): from the original 1.1 cycles/mm to 2.6 cycles/mm for the BD method and to 4.8 cycles/mm for the TD method) as well as contrast-to-noise ratio. With a properly chosen parameter, the TD method can resolve 1-mm size objects with the accurate temperature reading. The thermal image from the low-resolution IR camera enhanced by the TD method is comparable to that from a high-resolution IR camera. These results show that the TD method provides an effective way to improve the thermal image quality from a low-cost IR camera to monitor temperature of an object of 1-mm size, which meets the needed precision for advanced laser scanning protocols in photothermal/ablation therapies.

4D cone-beam computed tomography (CBCT) using a moving blocker for simultaneous radiation dose reduction and scatter correction.

Four-dimensional (4D) x-ray cone-beam computed tomography (CBCT) is important for a precise radiation therapy for lung cancer. Due to the repeated use and 4D acquisition over a course of radiotherapy, the radiation dose becomes a concern. Meanwhile, the scatter contamination in CBCT deteriorates image quality for treatment tasks. In this work, we propose the use of a moving blocker (MB) during the 4D CBCT acquisition ('4D MB') and to combine motion-compensated reconstruction to address these two issues simultaneously. In 4D MB CBCT, the moving blocker reduces the x-ray flux passing through the patient and collects the scatter information in the blocked region at the same time. The scatter signal is estimated from the blocked region for correction. Even though the number of projection views and projection data in each view are not complete for conventional reconstruction, 4D reconstruction with a total-variation (TV) constraint and a motion-compensated temporal constraint can utilize both spatial gradient sparsity and temporal correlations among different phases to overcome the missing data problem. The feasibility simulation studies using the 4D NCAT phantom showed that 4D MB with motion-compensated reconstruction with 1/3 imaging dose reduction could produce satisfactory images and achieve 37% improvement on structural similarity (SSIM) index and 55% improvement on root mean square error (RMSE), compared to 4D reconstruction at the regular imaging dose without scatter correction. For the same 4D MB data, 4D reconstruction outperformed 3D TV reconstruction by 28% on SSIM and 34% on RMSE. A study of synthetic patient data also demonstrated the potential of 4D MB to reduce the radiation dose by 1/3 without compromising the image quality. This work paves the way for more comprehensive studies to investigate the dose reduction limit offered by this novel 4D MB method using physical phantom experiments and real patient data based on clinical relevant metrics.