ALBERT

Description: The simultaneous removal of noise and preservation of the integrity of 3-D magnetic resonance (MR) images is a difficult and important task. In this paper, we consider characterizing MR images with 3-D operators, and present a novel 4-D transform-domain method termed ‘modified nonlocal tensor-SVD (MNL-tSVD)’ for MR image denoising. The proposed method is based on the grouping, hard-thresholding and aggregation paradigms, and can be viewed as a generalized nonlocal extension of tensor-SVD (t-SVD). By keeping MR images in its natural three-dimensional form, and collaboratively filtering similar patches, MNL-tSVD utilizes both the self-similarity property and 3-D structure of MR images to preserve more actual details and minimize the introduction of new artifacts. We show the adaptability of MNL-tSVD by incorporating it into a two-stage denoising strategy with a few adjustments. In addition, analysis of the relationship between MNL-tSVD and current the state-of-the-art 4-D transforms is given. Experimental comparisons over simulated and real brain data sets at different Rician noise levels show that MNL-tSVD can produce competitive performance compared with related approaches.

Description: Recent advances in imaging genetics produce large amounts of data including functional MRI images, single nucleotide polymorphisms (SNPs), and cognitive assessments. Understanding the complex interactions among these heterogeneous and complementary data has the potential to help with diagnosis and prevention of mental disorders. However, limited efforts have been made due to the high dimensionality, group structure, and mixed type of these data. In this paper, we present a novel method to detect conditional associations between imaging genetics data. We use projected distance correlation to build a conditional dependency graph among high-dimensional mixed data, and then use multiple testing to detect significant group level associations (e.g., regions of interest-gene). In addition, we introduce a scalable algorithm based on orthogonal greedy algorithm, yielding the greedy projected distance correlation (G-PDC). This can reduce the computational cost, which is critical for analyzing large volume of imaging genomics data. The results from our simulations demonstrate a higher degree of accuracy with G-PDC than distance correlation, Pearson’s correlation, and partial correlation, especially when the correlation is nonlinear. Finally, we apply our method to the Philadelphia Neurodevelopmental data cohort with 866 samples including fMRI images and SNP profiles. The results uncover several statistically significant and biologically interesting interactions, which are further validated with many existing studies. The MATLAB code is available at https://sites.google.com/site/jianfang86/gPDC .

Description: We present a multi-scale approach of tumor modeling in order to predict its evolution during radiotherapy. Within this context we focus on three different scales of tumor modeling: microscopic (individual cells in a voxel), mesoscopic (population of cells in a voxel) and macroscopic (whole tumor), with transition interfaces between these three scales. At the cellular level, the description is based on phase transfer probabilities in the cellular cycle. At the mesoscopic scale we represent populations of cells according to different stages in a cell cycle. Finally, at the macroscopic scale, the tumor description is based on the use of FDG PET image voxels. These three scales exist naturally: biological data are collected at the macroscopic scale, but the pathological behavior of the tumor is based on an abnormal cell-cycle at the microscopic scale. On the other hand, the introduction of a mesoscopic scale is essential in order to reduce the gap between the two extreme, in terms of resolution, description levels. It also reduces the computational burden of simulating a large number of individual cells. As an application of the proposed multi-scale model, we simulate the effect of oxygen on tumor evolution during radiotherapy. Two consecutive FDG PET images of 17 rectal cancer patients undergoing radiotherapy are used to simulate the tumor evolution during treatment. The simulated results are compared with those obtained on a third FDG PET image acquired two weeks after the beginning of the treatment.

Description: In this paper, we investigate the estimation of the maximum target registration error (TRE) magnitude of the target location while using point-based rigid registration in the image guided system. Under the uniform restriction of fiducial localization error (FLE) magnitude, we explicitly formulate the estimation as an optimization problem. Through analyzing the approximated problem which assumes the rigidity of the fiducial set holds with the perturbation of FLE, we present a strict lower bound for the maximum TRE magnitude. The simulations show that the lower bound is close to the actual maximum TRE magnitude for the target locations lying far away from the fiducial points. Unlike the expected TRE magnitude in which all fiducial points contribute, the lower bound is only related to the fiducial points serving as the vertices of the convex hull of the fiducial set. Our analysis provides a new perspective of investigating the problem of TRE estimation and is helpful for the surgeons to learn about the worst situation during using the image guided system.

Description: Brain tumors are the most common malignant neurologic tumors with the highest mortality and disability rate. Because of the delicate structure of the brain, the clinical use of several commonly used biopsy diagnosis is limited for brain tumors. Radiomics is an emerging technique for noninvasive diagnosis based on quantitative medical image analyses. However, current radiomics techniques are not standardized regarding feature extraction, feature selection, and decision making. In this paper, we propose a sparse representation-based radiomics (SRR) system for the diagnosis of brain tumors. First, we developed a dictionary learning- and sparse representation-based feature extraction method that exploits the statistical characteristics of the lesion area, leading to fine and more effective feature extraction compared with the traditional explicitly calculation-based methods. Then, we set up an iterative sparse representation method to solve the redundancy problem of the extracted features. Finally, we proposed a novel multi-feature collaborative sparse representation classification framework that introduces a new coefficient of regularization term to combine features from multi-modal images at the sparse representation coefficient level. Two clinical problems were used to validate the performance and usefulness of the proposed SRR system. One was the differential diagnosis between primary central nervous system lymphoma (PCNSL) and glioblastoma (GBM), and the other was isocitrate dehydrogenase 1 estimation for gliomas. The SRR system had superior PCNSL and GBM differentiation performance compared with some advanced imaging techniques and yielded 11% better performance for estimating IDH1 compared with the traditional radiomics methods.

Description: We present the first in vivo images of anisotropic conductivity distribution in the human head, measured at a frequency of approximately 10 Hz. We used magnetic resonance electrical impedance tomography techniques to encode phase changes caused by current flow within the head via two independent electrode pairs. These results were then combined with diffusion tensor imaging data to reconstruct full anisotropic conductivity distributions in 5-mm-thick slices of the brains of two participants. Conductivity values recovered in this paper were broadly consistent with literature values. We anticipate that this technique will be of use in many areas of neuroscience, most importantly in functional imaging via inverse electroencephalogram. Future studies will involve pulse sequence acceleration to maximize brain coverage and resolution.

Description: The most commonly used evaluation metrics for quality assessment of retinal vessel segmentation are sensitivity, specificity, and accuracy, which are based on pixel-to-pixel matching. However, due to the inter-observer problem that vessels annotated by different observers vary in both thickness and location, pixel-to-pixel matching is too restrictive to fairly evaluate the results of vessel segmentation. In this paper, the proposed skeletal similarity metric is constructed by comparing the skeleton maps generated from the reference and the source vessel segmentation maps. To address the inter-observer problem, instead of using a pixel-to-pixel matching strategy, each skeleton segment in the reference skeleton map is adaptively assigned with a searching range whose radius is determined based on its vessel thickness. Pixels in the source skeleton map located within the searching range are then selected for similarity calculation. The skeletal similarity consists of a curve similarity, which measures the structural similarity between the reference and the source skeleton maps and a thickness similarity, which measures the thickness consistency between the reference and the source vessel segmentation maps. In contrast to other metrics that provide a global score for the overall performance, we modify the definitions of true positive, false negative, true negative, and false positive based on the skeletal similarity, based on which sensitivity, specificity, accuracy, and other objective measurements can be constructed. More importantly, the skeletal similarity metric has better potential to be used as a pixelwise loss function for training deep learning models for retinal vessel segmentation. Through comparison of a set of examples, we demonstrate that the redefined metrics based on the skeletal similarity are more effective for quality evaluation, especially with greater tolerance to the inter-observer problem.

Description: We are developing a 1-mm 3 resolution, high-sensitivity positron emission tomography (PET) system for loco-regional cancer imaging. The completed system will comprise two $16.5times 10$ cm detector panels and contain 4 608 position sensitive avalanche photodiodes (PSAPDs) coupled to arrays of $0.9times 0.9times 1$ mm 3 LYSO crystal elements for a total of 294 912 crystal elements. For the first time, this paper summarizes the design and reports the performance of a significant portion of the final clinical PET system, comprising 1 536 PSAPDs, 98 304 crystal elements, and an active field-of-view (FOV) of $16.5times3.3$ cm. The sub-system performance parameters, such as energy, time, and spatial resolutions are predictive of the performance of the final system due to the modular design. Analysis of the multiplexed crystal flood histograms shows 84% of the crystal elements have〉99% crystal identification accuracy. The 511 keV photopeak energy resolution was 11.34±0.06% full-width half maximum (FWHM), and coincidence timing resolution was 13.92 ± 0.01 ns FWHM at 511 keV. The spatial resolution was measured using maximum likelihood expectation maximization reconstruction of a grid of point sources suspended in warm background. The averaged resolution over the central 6 cm of the FOV is 1.01 ± 0.13 mm in the X-direction, 1.84 ± 0.20 mm in the Y-direction, and 0.84 ± 0.11 mm in the Z-direction. Quantitative analysis of acquired micro-Derenzo phantom images shows better than 1.2 mm resolution at the center of the FOV, with subsequent resolution degradation in the y-direction toward the edge of the FOV caused by limited angle tomography effec- s.