Description: We present a novel reconstruction algorithm based on a general cone-beam CT forward model, which is capable of incorporating the blur and noise correlations that are exhibited in flat-panel CBCT measurement data. Specifically, the proposed model may include scintillator blur, focal-spot blur, and noise correlations due to light spread in the scintillator. The proposed algorithm (GPL-BC) uses a Gaussian Penalized-Likelihood objective function, which incorporates models of blur and correlated noise. In a simulation study, GPL-BC was able to achieve lower bias as compared with deblurring followed by FDK as well as a model-based reconstruction method without integration of measurement blur. In the same study, GPL-BC was able to achieve better line-pair reconstructions (in terms of segmented-image accuracy) as compared with deblurring followed by FDK, a model-based method without blur, and a model-based method with blur but not noise correlations. A prototype extremities quantitative cone-beam CT test-bench was used to image a physical sample of human trabecular bone. These data were used to compare reconstructions using the proposed method and model-based methods without blur and/or correlation to a registered $mu $ CT image of the same bone sample. The GPL-BC reconstructions resulted in more accurate trabecular bone segmentation. Multiple trabecular bone metrics, including trabecular thickness (Tb.Th.) were computed for each reconstruction approach as well as the $mu $ CT volume. The GPL-BC reconstruction provided the most accurate Tb.Th. measurement, 0.255 mm, as compared with the $mu $ CT derived value of 0.193 mm, followed by the GPL-B reconstruction, the GPL-I reconstruction, and then the FDK reconstruction (0.271 mm, 0.309 mm, and 0.335 mm, respect- vely).

Description: Computer-aided diagnosis (CAD) of retinal pathologies is a current active area in medical image analysis. Due to the increasing use of retinal optical coherence tomography (OCT) imaging technique, a CAD system in retinal OCT is essential to assist ophthalmologist in the early detection of ocular diseases and treatment monitoring. This paper presents a novel CAD system based on a multi-scale convolutional mixture of expert (MCME) ensemble model to identify normal retina, and two common types of macular pathologies, namely, dry age-related macular degeneration, and diabetic macular edema. The proposed MCME modular model is a data-driven neural structure, which employs a new cost function for discriminative and fast learning of image features by applying convolutional neural networks on multiple-scale sub-images. MCME maximizes the likelihood function of the training data set and ground truth by considering a mixture model, which tries also to model the joint interaction between individual experts by using a correlated multivariate component for each expert module instead of only modeling the marginal distributions by independent Gaussian components. Two different macular OCT data sets from Heidelberg devices were considered for the evaluation of the method, i.e., a local data set of OCT images of 148 subjects and a public data set of 45 OCT acquisitions. For comparison purpose, we performed a wide range of classification measures to compare the results with the best configurations of the MCME method. With the MCME model of four scale-dependent experts, the precision rate of 98.86%, and the area under the receiver operating characteristic curve (AUC) of 0.9985 were obtained on average.

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 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: We consider the challenges in estimating the state-related changes in brain connectivity networks with a large number of nodes. Existing studies use the sliding-window analysis or time-varying coefficient models, which are unable to capture both smooth and abrupt changes simultaneously, and rely on ad-hoc approaches to the high-dimensional estimation. To overcome these limitations, we propose a Markov-switching dynamic factor model, which allows the dynamic connectivity states in functional magnetic resonance imaging (fMRI) data to be driven by lower-dimensional latent factors. We specify a regime-switching vector autoregressive (SVAR) factor process to quantity the time-varying directed connectivity. The model enables a reliable, data-adaptive estimation of change-points of connectivity regimes and the massive dependencies associated with each regime. We develop a three-step estimation procedure: 1) extracting the factors using principal component analysis, 2) identifying connectivity regimes in a low-dimensional subspace based on the factor-based SVAR model, and 3) constructing high-dimensional state connectivity metrics based on the subspace estimates. Simulation results show that our estimator outperforms $K$ -means clustering of time-windowed coefficients, providing more accurate estimate of time-evolving connectivity. It achieves percentage of reduction in mean squared error by 60% when the network dimension is comparable to the sample size. When applied to the resting-state fMRI data, our method successfully identifies modular organization in the resting-statenetworksin consistencywith other studies. It further reveals changes in brain states with variations across subjects and distinct large-scale directed connectivity patterns across states.

Description: Inhomogeneities in the transmit radio frequency magnetic field ( ${text{B}}_{1}^{+}$ ) reduce the quality of magnetic resonance (MR) images. This quality can be improved by using high-permittivity pads that tailor the ${text{B}}_{1}^{+}$ fields. The design of an optimal pad is application-specific and not straightforward and would therefore benefit from a systematic optimization approach. In this paper, we propose such a method to efficiently design dielectric pads. To this end, a projection-based model order reduction technique is used that significantly decreases the dimension of the design problem. Subsequently, the resulting reduced-order model is incorporated in an optimization method in which a desired field in a region of interest can be set. The method is validated by designing a pad for imaging the cerebellum at 7 T. The optimal pad that is found is used in an MR measurement to demonstrate its effectiveness in improving the image quality.

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.