ALBERT

Description: Low-dose computed tomography (LDCT) images are often severely degraded by amplified mottle noise and streak artifacts. These artifacts are often hard to suppress without introducing tissue blurring effects. In this paper, we propose to process LDCT images using a novel image-domain algorithm called “artifact suppressed dictionary learning (ASDL).” In this ASDL method, orientation and scale information on artifacts is exploited to train artifact atoms, which are then combined with tissue feature atoms to build three discriminative dictionaries. The streak artifacts are cancelled via a discriminative sparse representation operation based on these dictionaries. Then, a general dictionary learning processing is applied to further reduce the noise and residual artifacts. Qualitative and quantitative evaluations on a large set of abdominal and mediastinum CT images are carried out and the results show that the proposed method can be efficiently applied in most current CT systems.

Description: When segmenting intraretinal layers from multiple optical coherence tomography (OCT) images forming a mosaic or a set of repeated scans, it is attractive to exploit the additional information from the overlapping areas rather than discarding it as redundant, especially in low contrast and noisy images. However, it is currently not clear how to effectively combine the multiple information sources available in the areas of overlap. In this paper, we propose a novel graph-theoretic method for multi-surface multi-field co-segmentation of intraretinal layers, assuring consistent segmentation of the fields across the overlapped areas. After 2-D en-face alignment, all the fields are segmented simultaneously, imposing a priori soft interfield-intrasurface constraints for each pair of overlapping fields. The constraints penalize deviations from the expected surface height differences, taken to be the depth-axis shifts that produce the maximum cross-correlation of pairwise-overlapped areas. The method's accuracy and reproducibility are evaluated qualitatively and quantitatively on 212 OCT images (20 nine-field, 32 single-field acquisitions) from 26 patients with glaucoma. Qualitatively, the obtained thickness maps show no stitching artifacts, compared to pronounced stitches when the fields are segmented independently. Quantitatively, two ophthalmologists manually traced four intraretinal layers on 10 patients, and the average error ( $4.58 pm 1.46 ~ {rm mu m}$ ) was comparable to the average difference between the observers ( $5.86pm 1.72 ~ {rm mu m}$ ). Furthermore, we show the benefit of the proposed approach in co-segmenting longitudinal scans. As opposed to segmenting layers in each of the fields independently, the proposed co-segmentation method obtains consistent segmentations across the overlapped areas, producing accurate, reproducibl- , and artifact-free results.

Description: Diffuse optical tomography (DOT) is a noninvasive technique which measures hemodynamic changes in the tissue with near infrared light, which has been increasingly used to study brain functions. Due to the nature of light propagation in the tissue, the reconstruction problem is severely ill-posed. For linearized DOT problems, sparsity regularization has achieved promising results over conventional Tikhonov regularization in recent experimental research. As extensions to standard sparsity, it is widely known that structured sparsity based methods are often superior in terms of reconstruction accuracy, when the data follows some structures. In this paper, we exploit the structured sparsity of diffuse optical images. Based on the functional specialization of the brain, it is observed that the in vivo absorption changes caused by a specific brain function would be clustered in certain region(s) and not randomly distributed. Thus, a new algorithm is proposed for this clustered sparsity reconstruction (CSR). Results of numerical simulations and phantom experiments have demonstrated the superiority of the proposed method over the state-of-the-art methods. An example from human in vivo measurements further confirmed the advantages of the proposed CSR method.

Description: Provides instructions and guidelines to prospective authors who wish to submit manuscripts.

Description: In ophthalmology, it is now common practice to record every surgical procedure and to archive the resulting videos for documentation purposes. In this paper, we present a solution to automatically segment and categorize surgical tasks in real-time during the surgery, using the video recording. The goal would be to communicate information to the surgeon in due time, such as recommendations to the less experienced surgeons. The proposed solution relies on the content-based video retrieval paradigm: it reuses previously archived videos to automatically analyze the current surgery, by analogy reasoning. Each video is segmented, in real-time, into an alternating sequence of idle phases, during which no clinically-relevant motions are visible, and action phases. As soon as an idle phase is detected, the previous action phase is categorized and the next action phase is predicted. A conditional random field is used for categorization and prediction. The proposed system was applied to the automatic segmentation and categorization of cataract surgery tasks. A dataset of 186 surgeries, performed by ten different surgeons, was manually annotated: ten possibly overlapping surgical tasks were delimited in each surgery. Using the content of action phases and the duration of idle phases as sources of evidence, an average recognition performance of $A_{z} = 0.832 pm 0.070$ was achieved.

Description: In this study, we propose a novel pathological lung segmentation method that takes into account neighbor prior constraints and a novel pathology recognition system. Our proposed framework has two stages; during stage one, we adapted the fuzzy connectedness (FC) image segmentation algorithm to perform initial lung parenchyma extraction. In parallel, we estimate the lung volume using rib-cage information without explicitly delineating lungs. This rudimentary, but intelligent lung volume estimation system allows comparison of volume differences between rib cage and FC based lung volume measurements. Significant volume difference indicates the presence of pathology, which invokes the second stage of the proposed framework for the refinement of segmented lung. In stage two, texture-based features are utilized to detect abnormal imaging patterns (consolidations, ground glass, interstitial thickening, tree-inbud, honeycombing, nodules, and micro-nodules) that might have been missed during the first stage of the algorithm. This refinement stage is further completed by a novel neighboring anatomy-guided segmentation approach to include abnormalities with weak textures, and pleura regions. We evaluated the accuracy and efficiency of the proposed method on more than 400 CT scans with the presence of a wide spectrum of abnormalities. To our best of knowledge, this is the first study to evaluate all abnormal imaging patterns in a single segmentation framework. The quantitative results show that our pathological lung segmentation method improves on current standards because of its high sensitivity and specificity and may have considerable potential to enhance the performance of routine clinical tasks.

Description: In pinhole single photon emission computed tomography (SPECT), multi-pinhole collimators can increase sensitivity but may lead to projection overlap, or multiplexing, which can cause image artifacts. In this work, we explore whether a stacked-detector configuration with a germanium and a silicon detector, used with $^{123}{rm I}$ (27–32, 159 keV), where little multiplexing occurs in the Si projections, can reduce image artifacts caused by highly-multiplexed Ge projections. Simulations are first used to determine a reconstruction method that combines the Si and Ge projections to maximize image quality. Next, simulations of different pinhole configurations (varying projection multiplexing) in conjunction with digital phantoms are used to examine whether additional Si projections mitigate artifacts from the multiplexing in the Ge projections. Reconstructed images using both Si and Ge data are compared to those using Ge data alone. Normalized mean-square error and normalized standard deviation provide a quantitative evaluation of reconstructed images' error and noise, respectively, and are used to evaluate the impact of the additional nonmultiplexed data on image quality. For a qualitative comparison, the differential point response function is used to examine multiplexing artifacts. Results show that in cases of highly-multiplexed Ge projections, the addition of low-multiplexed Si projections helps to reduce image artifacts both quantitatively and qualitatively.

Description: An easily implementable tissue cancellation method for dual energy mammography is proposed to reduce anatomical noise and enhance lesion visibility. For dual energy calibration, the images of an imaging object are directly mapped onto the images of a customized calibration phantom. Each pixel pair of the low and high energy images of the imaging object was compared to pixel pairs of the low and high energy images of the calibration phantom. The correspondence was measured by absolute difference between the pixel values of imaged object and those of the calibration phantom. Then the closest pixel pair of the calibration phantom images is marked and selected. After the calibration using direct mapping, the regions with lesion yielded different thickness from the background tissues. Taking advantage of the different thickness, the visibility of cancerous lesions was enhanced with increased contrast-to-noise ratio, depending on the size of lesion and breast thickness. However, some tissues near the edge of imaged object still remained after tissue cancellation. These remaining residuals seem to occur due to the heel effect, scattering, nonparallel X-ray beam geometry and Poisson distribution of photons. To improve its performance further, scattering and the heel effect should be compensated.

Description: Provides a listing of current staff, committee members and society officers.

Description: Diffusion tensor imaging is widely used in brain connectivity research. As more and more studies recruit large numbers of subjects, it is important to design registration methods which are not only theoretically rigorous, but also computationally efficient. However, the requirement of reorienting diffusion tensors complicates and considerably slows down registration procedures, due to the correlated impacts of registration forces at adjacent voxel locations. Based on the diffeomorphic Demons algorithm (Vercauteren , 2009), we propose a fast local trust region algorithm for handling inseparable registration forces for quadratic energy functions. The method guarantees that, at any time and at any voxel location, the velocity is always within its local trust region. This local regularization allows efficient calculation of the transformation update with numeric integration instead of completely solving a large linear system at every iteration. It is able to incorporate exact reorientation and regularization into the velocity optimization, and preserve the linear complexity of the diffeomorphic Demons algorithm. In an experiment with 84 diffusion tensor images involving both pair-wise and group-wise registrations, the proposed algorithm achieves better registration in comparison with other methods solving large linear systems (Yeo , 2009). At the same time, this algorithm reduces the computation time and memory demand tenfold.

Description: Routine ultrasound exam in the second and third trimesters of pregnancy involves manually measuring fetal head and brain structures in 2-D scans. The procedure requires a sonographer to find the standardized visualization planes with a probe and manually place measurement calipers on the structures of interest. The process is tedious, time consuming, and introduces user variability into the measurements. This paper proposes an automatic fetal head and brain (AFHB) system for automatically measuring anatomical structures from 3-D ultrasound volumes. The system searches the 3-D volume in a hierarchy of resolutions and by focusing on regions that are likely to be the measured anatomy. The output is a standardized visualization of the plane with correct orientation and centering as well as the biometric measurement of the anatomy. The system is based on a novel framework for detecting multiple structures in 3-D volumes. Since a joint model is difficult to obtain in most practical situations, the structures are detected in a sequence, one-by-one. The detection relies on Sequential Estimation techniques, frequently applied to visual tracking. The interdependence of structure poses and strong prior information embedded in our domain yields faster and more accurate results than detecting the objects individually. The posterior distribution of the structure pose is approximated at each step by sequential Monte Carlo. The samples are propagated within the sequence across multiple structures and hierarchical levels. The probabilistic model helps solve many challenges present in the ultrasound images of the fetus such as speckle noise, signal drop-out, shadows caused by bones, and appearance variations caused by the differences in the fetus gestational age. This is possible by discriminative learning on an extensive database of scans comprising more than two thousand volumes and more than thirteen thousand annotations. The average difference between ground truth and automatic - easurements is below 2 mm with a running time of 6.9 s (GPU) or 14.7 s (CPU). The accuracy of the AFHB system is within inter-user variability and the running time is fast, which meets the requirements for clinical use.

Description: In clinical practice, echocardiographers are often unkeen to make the significant time investment to make additional multiple measurements of Doppler velocity. Main hurdle to obtaining multiple measurements is the time required to manually trace a series of Doppler traces. To make it easier to analyze more beats, we present the description of an application system for automated aortic Doppler envelope quantification, compatible with a range of hardware platforms. It analyses long Doppler strips, spanning many heartbeats, and does not require electrocardiogram to separate individual beats. We tested its measurement of velocity-time-integral and peak-velocity against the reference standard defined as the average of three experts who each made three separate measurements. The automated measurements of velocity-time-integral showed strong correspondence $({rm R}^{2} =0.94)$ and good Bland–Altman agreement $({rm SD} = 1.39~{hbox {cm}})$ with the reference consensus expert values, and indeed performed as well as the individual experts ( ${rm R}^{2} = 0.90$ to 0.96, ${rm SD} = 1.05$ to 1.53 cm). The same performance was observed for peak-velocities; ( ${rm R}^{2} =0.98$ , ${rm SD} = 3.07~{rm cm/s}$ ) and ( ${rm R}^{2} = 0.93$ to 0.98, ${rm SD} = 2.96$ to 5.18 cm/s). This automated technology allows $〉 10$ times as many beats to be analyzed compared to the conventional manua- approach. This would make clinical and research protocols more precise for the same operator effort.

Description: We study the problem of joint registration and deformation analysis of endometrial tissue using 3D magnetic resonance imaging (MRI) and 2D trans-vaginal ultrasound (TVUS) measurements. In addition to the different imaging techniques involved in the two modalities, this problem is complicated due to: 1) different patient pose during MRI and TVUS observations, 2) the 3D nature of MRI and 2D nature of TVUS measurements, 3) the unknown intersecting plane for TVUS in MRI volume, and 4) the potential deformation of endometrial tissue during TVUS measurement process. Focusing on the shape of the tissue, we use expert manual segmentation of its boundaries in the two modalities and apply, with modification, recent developments in shape analysis of parametric surfaces to this problem. First, we extend the 2D TVUS curves to generalized cylindrical surfaces through replication, and then we compare them with MRI surfaces using elastic shape analysis. This shape analysis provides a simultaneous registration (optimal reparameterization) and deformation (geodesic) between any two parametrized surfaces. Specifically, it provides optimal curves on MRI surfaces that match with the original TVUS curves. This framework results in an accurate quantification and localization of the deformable endometrial cells for radiologists, and growth characterization for gynecologists and obstetricians. We present experimental results using semi-synthetic data and real data from patients to illustrate these ideas.

Description: Repetitive and alternating lower limb movements are a specific component of human gait. Due to technical challenges, the neural mechanisms underlying such movements have not been previously studied with functional magnetic resonance imaging. In this study, we present a novel treadmill device employed to investigate the kinematics and the brain activation patterns involved in alternating and repetitive movements of the lower limbs. Once inside the scanner, 19 healthy subjects were guided by two visual cues and instructed to perform a motor task which involved repetitive and alternating movements of both lower limbs while selecting their individual comfortable amplitude on the treadmill. The device facilitated the performance of coordinated stepping while registering the concurrent lower-limb displacements, which allowed us to quantify some movement primary kinematic features such as amplitude and frequency. During stepping, significant blood oxygen level dependent signal increases were observed bilaterally in primary and secondary sensorimotor cortex, the supplementary motor area, premotor cortex, prefrontal cortex, superior and inferior parietal lobules, putamen and cerebellum, regions that are known to be involved in lower limb motor control. Brain activations related to individual adjustments during motor performance were identified in a right lateralized network including striatal, extrastriatal, and fronto-parietal areas.

Description: 2D/3D registration of patient vasculature from preinterventional computed tomography angiography (CTA) to interventional X-ray angiography is of interest to improve guidance in percutaneous coronary interventions. In this paper we present a novel feature based 2D/3D registration framework, that is based on probabilistic point correspondences, and show its usefulness on aligning 3D coronary artery centerlines derived from CTA images with their 2D projection derived from interventional X-ray angiography. The registration framework is an extension of the Gaussian mixture model (GMM) based point-set registration to the 2D/3D setting, with a modified distance metric. We also propose a way to incorporate orientation in the registration, and show its added value for artery registration on patient datasets as well as in simulation experiments. The oriented GMM registration achieved a median accuracy of 1.06 mm, with a convergence rate of 81% for nonrigid vessel centerline registration on 12 patient datasets, using a statistical shape model. The method thereby outperformed the iterative closest point algorithm, the GMM registration without orientation, and two recently published methods on 2D/3D coronary artery registration.

Description: In a positron emission tomography (PET) study, the local uptake of the tracer is dependent on vascular delivery and retention. For dynamic studies the measured uptake time-course information can be best interpreted when knowledge of the time-course of tracer in the blood is available. This is certainly true for the most established tracers such as $^{18}F$ -Fluorodeoxyglucose (FDG) and $^{15}O$ -Water (H $_2$ O). Since direct sampling of blood as part of PET studies is increasingly impractical, there is ongoing interest in image-extraction of blood time-course information. But analysis of PET-measured blood pool signals is complicated because they will typically involve a combination of arterial, venous and tissue information. Thus, a careful appreciation of these components is needed to interpret the available data. To facilitate this process, we propose a novel Markov chain model for representation of the circulation of a tracer atom in the body. The model represents both arterial and venous time-course patterns. Under reasonable conditions equilibration of tracer activity in arterial and venous blood is achieved by the end of the PET study—consistent with empirical measurement. Statistical inference for Markov model parameters is a challenge. A penalized nonlinear least squares process, incorporating a generalized cross-validation score, is proposed. Random effects analysis is used to adaptively specify the structure of the penalty function based on historical samples of directly measured blood data. A collection of arterially sampled data from PET studies with FDG and H $_2$ O is used to illustrate the methodology. These data analyses are highly supportive of the overall modeling approach. An adapt- tion of the model to the problem of extraction of arterial blood signals from imaging data is also developed and promising preliminary results for cerebral and thoracic imaging studies with FDG and H $_2$ O are obtained.

Description: Digital image-based elasto-tomography (DIET) is a prototype system for breast cancer screening. A breast is imaged while being vibrated, and the observed surface motion is used to infer the internal stiffness of the breast, hence identifying tumors. This paper describes a computer vision system for accurately measuring 3-D surface motion. A model-based segmentation is used to identify the profile of the breast in each image, and the 3-D surface is reconstructed by fitting a model to the profiles. The surface motion is measured using a modern optical flow implementation customized to the application, then trajectories of points on the 3-D surface are given by fusing the optical flow with the reconstructed surfaces. On data from human trials, the system is shown to exceed the performance of an earlier marker-based system at tracking skin surface motion. We demonstrate that the system can detect a 10 mm tumor in a silicone phantom breast.

Description: Gold immunochromatographic strip assay provides a rapid, simple, single-copy and on-site way to detect the presence or absence of the target analyte. This paper aims to develop a method for accurately segmenting the test line and control line of the gold immunochromatographic strip (GICS) image for quantitatively determining the trace concentrations in the specimen, which can lead to more functional information than the traditional qualitative or semi-quantitative strip assay. The canny operator as well as the mathematical morphology method is used to detect and extract the GICS reading-window. Then, the test line and control line of the GICS reading-window are segmented by the cellular neural network (CNN) algorithm, where the template parameters of the CNN are designed by the switching particle swarm optimization (SPSO) algorithm for improving the performance of the CNN. It is shown that the SPSO-based CNN offers a robust method for accurately segmenting the test and control lines, and therefore serves as a novel image methodology for the interpretation of GICS. Furthermore, quantitative comparison is carried out among four algorithms in terms of the peak signal-to-noise ratio. It is concluded that the proposed CNN algorithm gives higher accuracy and the CNN is capable of parallelism and analog very-large-scale integration implementation within a remarkably efficient time.

Description: Prostate cancer is one of the major causes of cancer death for men in the western world. Magnetic resonance imaging (MRI) is being increasingly used as a modality to detect prostate cancer. Therefore, computer-aided detection of prostate cancer in MRI images has become an active area of research. In this paper we investigate a fully automated computer-aided detection system which consists of two stages. In the first stage, we detect initial candidates using multi-atlas-based prostate segmentation, voxel feature extraction, classification and local maxima detection. The second stage segments the candidate regions and using classification we obtain cancer likelihoods for each candidate. Features represent pharmacokinetic behavior, symmetry and appearance, among others. The system is evaluated on a large consecutive cohort of 347 patients with MR-guided biopsy as the reference standard. This set contained 165 patients with cancer and 182 patients without prostate cancer. Performance evaluation is based on lesion-based free-response receiver operating characteristic curve and patient-based receiver operating characteristic analysis. The system is also compared to the prospective clinical performance of radiologists. Results show a sensitivity of 0.42, 0.75, and 0.89 at 0.1, 1, and 10 false positives per normal case. In clinical workflow the system could potentially be used to improve the sensitivity of the radiologist. At the high specificity reading setting, which is typical in screening situations, the system does not perform significantly different from the radiologist and could be used as an independent second reader instead of a second radiologist. Furthermore, the system has potential in a first-reader setting.

Description: Quantification of regional myocardial motion and deformation from cardiac ultrasound is fostering considerable research efforts. Despite the tremendous improvements done in the field, all existing approaches still face a common limitation which is intrinsically connected with the formation of the ultrasound images. Specifically, the reduced lateral resolution and the absence of phase information in the lateral direction highly limit the accuracy in the computation of lateral displacements. In this context, this paper introduces a novel setup for the estimation of cardiac motion with ultrasound. The framework includes an unconventional beamforming technique and a dedicated motion estimation algorithm. The beamformer aims at introducing phase information in the lateral direction by producing transverse oscillations. The estimator directly exploits the phase information in the two directions by decomposing the image into two 2-D single-orthant analytic signals. An in silico evaluation of the proposed framework is presented on five ultra-realistic simulated echocardiographic sequences, where the proposed motion estimator is contrasted against other two phase-based solutions exploiting the presence of transverse oscillations and against block-matching on standard images. An implementation of the new beamforming strategy on a research ultrasound platform is also shown along with a preliminary in vivo evaluation on one healthy subject.

Description: In this paper different model-based methods of classification of global patterns in dermoscopic images are proposed. Global patterns identification is included in the pattern analysis framework, the melanoma diagnosis method most used among dermatologists. The modeling is performed in two senses: first a dermoscopic image is modeled by a finite symmetric conditional Markov model applied to $L^{ast}a^{ast}b^{ast}$ color space and the estimated parameters of this model are treated as features. In turn, the distribution of these features are supposed that follow different models along a lesion: a Gaussian model, a Gaussian mixture model, and a bag-of-features histogram model. For each case, the classification is carried out by an image retrieval approach with different distance metrics. The main objective is to classify a whole pigmented lesion into three possible patterns: globular, homogeneous, and reticular. An extensive evaluation of the performance of each method has been carried out on an image database extracted from a public Atlas of Dermoscopy. The best classification success rate is achieved by the Gaussian mixture model-based method with a 78.44% success rate in average. In a further evaluation the multicomponent pattern is analyzed obtaining a 72.91% success rate.

Description: Describes the above-named upcoming conference event. May include topics to be covered or calls for papers.

Description: Multispectral optoacoustic tomography (MSOT) utilizes broadband ultrasound detection for imaging biologically-relevant optical absorption features at a range of scales. Due to the multiscale and multispectral features of the technology, MSOT comes with distinct requirements in implementation and data analysis. In this work, we investigate the interplay between scale, which depends on ultrasonic detection frequency, and optical multispectral spectral analysis, two dimensions that are unique to MSOT and represent a previously unexplored challenge. We show that ultrasound frequency-dependent artifacts suppress multispectral features and complicate spectral analysis. In response, we employ a wavelet decomposition to perform spectral unmixing on a per-scale basis (or per ultrasound frequency band) and showcase imaging of fine-scale features otherwise hidden by low frequency components. We explain the proposed algorithm by means of simple simulations and demonstrate improved performance in imaging data of blood vessels in human subjects.

Description: Provides instructions and guidelines to prospective authors who wish to submit manuscripts.

Description: The purpose of this study was to increase the sensitivity of XFCT imaging by optimizing the data acquisition geometry for reduced scatter X-rays. The placement of detectors and detector energy window were chosen to minimize scatter X-rays. We performed both theoretical calculations and Monte Carlo simulations of this optimized detector configuration on a mouse-sized phantom containing various gold concentrations. The sensitivity limits were determined for three different X-ray spectra: a monoenergetic source, a Gaussian source, and a conventional X-ray tube source. Scatter X-rays were minimized using a backscatter detector orientation (scatter direction $〉 110^{circ }$ to the primary X-ray beam). The optimized configuration simultaneously reduced the number of detectors and improved the image signal-to-noise ratio. The sensitivity of the optimized configuration was 10 $mu {rm g}/{rm mL}$ (10 pM) at 2 mGy dose with the mono-energetic source, which is an order of magnitude improvement over the unoptimized configuration (102 pM without the optimization). Similar improvements were seen with the Gaussian spectrum source and conventional X-ray tube source. The optimization improvements were predicted in the theoretical model and also demonstrated in simulations. The sensitivity of XFCT imaging can be enhanced by an order of magnitude with the data acquisition optimization, greatly enhancing the potential of this modality for future use in clinical molecular imaging.

Description: Optoacoustic tomography (OAT), also known as photoacoustic tomography, is an emerging computed biomedical imaging modality that exploits optical contrast and ultrasonic detection principles. Iterative image reconstruction algorithms that are based on discrete imaging models are actively being developed for OAT due to their ability to improve image quality by incorporating accurate models of the imaging physics, instrument response, and measurement noise. In this work, we investigate the use of discrete imaging models based on Kaiser–Bessel window functions for iterative image reconstruction in OAT. A closed-form expression for the pressure produced by a Kaiser–Bessel function is calculated, which facilitates accurate computation of the system matrix. Computer-simulation and experimental studies are employed to demonstrate the potential advantages of Kaiser–Bessel function-based iterative image reconstruction in OAT.

Description: The multi-channel nature of digital histopathological images presents an opportunity to exploit the correlated color channel information for better image modeling. Inspired by recent work in sparsity for single channel image classification, we propose a new simultaneous sparsity model for multi-channel histopathological image representation and classification (SHIRC). Essentially, we represent a histopathological image as a sparse linear combination of training examples under suitable channel-wise constraints. Classification is performed by solving a newly formulated simultaneous sparsity-based optimization problem. A practical challenge is the correspondence of image objects (cellular and nuclear structures) at different spatial locations in the image. We propose a robust locally adaptive variant of SHIRC (LA-SHIRC) to tackle this issue. Experiments on two challenging real-world image data sets: 1) mammalian tissue images acquired by pathologists of the animal diagnostics lab (ADL) at Pennsylvania State University, and 2) human intraductal breast lesions, reveal the merits of our proposal over state-of-the-art alternatives. Further, we demonstrate that LA-SHIRC exhibits a more graceful decay in classification accuracy against the number of training images which is highly desirable in practice where generous training per class is often not available.

Description: Advertisement: This publication offers open access options for authors. IEEE open access publishing.

Description: Statistical analysis of motion patterns of body joints is potentially useful for detecting and quantifying pathologies. However, building a statistical motion model across different subjects remains a challenging task, especially for a complex joint like the wrist. We present a novel framework for simultaneous registration and segmentation of multiple 3-D (CT or MR) volumes of different subjects at various articulated positions. The framework starts with a pose model generated from 3-D volumes captured at different articulated positions of a single subject (template). This initial pose model is used to register the template volume to image volumes from new subjects. During this process, the Grow-Cut algorithm is used in an iterative refinement of the segmentation of the bone along with the pose parameters. As each new subject is registered and segmented, the pose model is updated, improving the accuracy of successive registrations. We applied the algorithm to CT images of the wrist from 25 subjects, each at five different wrist positions and demonstrated that it performed robustly and accurately. More importantly, the resulting segmentations allowed a statistical pose model of the carpal bones to be generated automatically without interaction. The evaluation results show that our proposed framework achieved accurate registration with an average mean target registration error of $0.34 pm 0.27$ mm. The automatic segmentation results also show high consistency with the ground truth obtained semi-automatically. Furthermore, we demonstrated the capability of the resulting statistical pose and shape models by using them to generate a measurement tool for scaphoid-lunate dissociation diagnosis, which achieved 90% sensitivity and specificity.

Description: The aim of this work is the evaluation of the design for a nonconventional PET scanner, the voxel imaging PET (VIP), based on pixelated room-temperature CdTe detectors yielding a true 3-D impact point with a density of 450 ${hbox {channels}}/{hbox {cm}}^{3}$ , for a total 6 336 000 channels in a seamless ring shaped volume. The system is simulated and evaluated following the prescriptions of the NEMA NU 2-2001 and the NEMA NU 4-2008 standards. Results show that the excellent energy resolution of the CdTe detectors (1.6% for 511 keV photons), together with the small voxel pitch (1 $, times ,$ 1 $,times ,$ 2 ${hbox {mm}}^{3}$ ), and the crack-free ring geometry, give the design the potential to overcome the current limitations of PET scanners and to approach the intrinsic image resolution limits set by physics. The VIP is expected to reach a competitive sensitivity and a superior signal purity with respect to values commonly quoted for state-of-the-art scintillating crystal PETs. The system can provide 14 cps/kBq with a scatter fraction of 3.95% and 21 cps/kBq with a scatter fraction of 0.73% according to NEMA NU 2-2001 and NEMA NU 4-2008, respectively. The calculated NEC curve has a peak value of 122 kcps at 5.3 kBq/mL for NEMA NU 2-2001 and 908 kcps at 1.6 MBq/mL for NEMA NU 4-2008. The proposed scanner can achieve an image resolution of $sim 1~{hbox {mm}}$ full-width at half-maximum in all directions. The virtually noise-free data sample leads to direct positive impact on the quality of the reconstructed images. As a consequence, high-quality high-resolution images can be obtained with significantly lower number of events compared to - onventional scanners. Overall, simulation results suggest the VIP scanner can be operated either at normal dose for fast scanning and high patient throughput, or at low dose to decrease the patient radioactivity exposure. The design evaluation presented in this work is driving the development and the optimization of a fully operative prototype to prove the feasibility of the VIP concept.

Description: Diffusion tensor imaging (DTI) is unable to represent the diffusion signal arising from multiple crossing fascicles and freely diffusing water molecules. Generative models of the diffusion signal, such as multi-fascicle models, overcome this limitation by providing a parametric representation for the signal contribution of each population of water molecules. These models are of great interest in population studies to characterize and compare the brain microstructural properties. Central to population studies is the construction of an atlas and the registration of all subjects to it. However, the appropriate definition of registration and atlasing methods for multi-fascicle models have proven challenging. This paper proposes a mathematical framework to register and analyze multi-fascicle models. Specifically, we define novel operators to achieve interpolation, smoothing and averaging of multi-fascicle models. We also define a novel similarity metric to spatially align multi-fascicle models. Our framework enables simultaneous comparisons of different microstructural properties that are confounded in conventional DTI. The framework is validated on multi-fascicle models from 24 healthy subjects and 38 patients with tuberous sclerosis complex, 10 of whom have autism. We demonstrate the use of the multi-fascicle models registration and analysis framework in a population study of autism spectrum disorder.

Description: Here, we present an innovative imaging technology for breast cancer using gamma-ray stimulated spectroscopy based on the nuclear resonance fluorescence (NRF) technique. In NRF, a nucleus of a given isotope selectively absorbs gamma rays with energy exactly equal to one of its quantized energy states, emitting an outgoing gamma ray with energy nearly identical to that of the incident gamma ray. Due to its application of NRF, gamma-ray stimulated spectroscopy is sensitive to trace element concentration changes, which are suspected to occur at early stages of breast cancer, and therefore can be potentially used to noninvasively detect and diagnose cancer in its early stages. Using Monte-Carlo simulations, we have designed and demonstrated an imaging system that uses gamma-ray stimulated spectroscopy for visualizing breast cancer. We show that gamma-ray stimulated spectroscopy is able to visualize breast cancer lesions based primarily on the differences in the concentrations of trace elements between diseased and healthy tissue, rather than differences in density that are crucial for X-ray mammography. The technique shows potential for early breast cancer detection; however, improvements are needed in gamma-ray laser technology for the technique to become a clinically feasible method of detecting and diagnosing cancer at early stages.

Description: The National Library of Medicine (NLM) is developing a digital chest X-ray (CXR) screening system for deployment in resource constrained communities and developing countries worldwide with a focus on early detection of tuberculosis. A critical component in the computer-aided diagnosis of digital CXRs is the automatic detection of the lung regions. In this paper, we present a nonrigid registration-driven robust lung segmentation method using image retrieval-based patient specific adaptive lung models that detects lung boundaries, surpassing state-of-the-art performance. The method consists of three main stages: 1) a content-based image retrieval approach for identifying training images (with masks) most similar to the patient CXR using a partial Radon transform and Bhattacharyya shape similarity measure, 2) creating the initial patient-specific anatomical model of lung shape using SIFT-flow for deformable registration of training masks to the patient CXR, and 3) extracting refined lung boundaries using a graph cuts optimization approach with a customized energy function. Our average accuracy of 95.4% on the public JSRT database is the highest among published results. A similar degree of accuracy of 94.1% and 91.7% on two new CXR datasets from Montgomery County, MD, USA, and India, respectively, demonstrates the robustness of our lung segmentation approach.

Description: Image quality in 3-D optoacoustic (photoacoustic) tomography is greatly influenced by both the measurement system, in particular the number and spatial arrangement of ultrasound sensors, and the ability to account for the spatio-temporal response of the sensor element(s) in the reconstruction algorithm. Herein we present a reconstruction procedure based on the inversion of a time-domain forward model incorporating the spatial impulse response due to the shape of the transducer, which is subsequently applied in a tomographic system based on a translation-rotation scan of a linear detector array. The proposed method was also adapted to cope with the data-intensive requirements of high-resolution volumetric optoacoustic imaging. The processing of 2 $cdot$ 10 $^4$ individual signals resulted in well-resolved images of both $sim 200 mu{rm m}$ absorbers in phantoms and complex vascular structures in biological tissue. The results reported herein demonstrate that the introduced model-based methodology exhibits a better contrast and resolution than standard back-projection and model-based algorithms that assume point detectors. Moreover, the capability of handling large datasets anticipates that model-based methods incorporating the sensor properties can become standard practice in volumetric opto acoustic image formation.

Description: It is known that a reduction of the field-of-view in 3-D X-ray imaging is proportional to a reduction in radiation dose. The resulting truncation, however, is incompatible with conventional reconstruction algorithms. Recently, a novel method for region of interest reconstruction that uses neither prior knowledge nor extrapolation has been published, named approximated truncation robust algorithm for computed tomography (ATRACT). It is based on a decomposition of the standard ramp filter into a 2-D Laplace filtering and a 2-D Radon-based residual filtering step. In this paper, we present two variants of the original ATRACT. One is based on expressing the residual filter as an efficient 2-D convolution with an analytically derived kernel. The second variant is to apply ATRACT in 1-D to further reduce computational complexity. The proposed algorithms were evaluated by using a reconstruction benchmark, as well as two clinical data sets. The results are encouraging since the proposed algorithms achieve a speed-up factor of up to 245 compared to the 2-D Radon-based ATRACT. Reconstructions of high accuracy are obtained, e.g., even real-data reconstruction in the presence of severe truncation achieve a relative root mean square error of as little as 0.92% with respect to nontruncated data.

Description: A new step-by-step comprehensive MR physics simulator (MRISIMUL) of the Bloch equations is presented. The aim was to develop a magnetic resonance imaging (MRI) simulator that makes no assumptions with respect to the underlying pulse sequence and also allows for complex large-scale analysis on a single computer without requiring simplifications of the MRI model. We hypothesized that such a simulation platform could be developed with parallel acceleration of the executable core within the graphic processing unit (GPU) environment. MRISIMUL integrates realistic aspects of the MRI experiment from signal generation to image formation and solves the entire complex problem for densely spaced isochromats and for a densely spaced time axis. The simulation platform was developed in MATLAB whereas the computationally demanding core services were developed in CUDA-C. The MRISIMUL simulator imaged three different computer models: a user-defined phantom, a human brain model and a human heart model. The high computational power of GPU-based simulations was compared against other computer configurations. A speedup of about 228 times was achieved when compared to serially executed C-code on the CPU whereas a speedup between 31 to 115 times was achieved when compared to the OpenMP parallel executed C-code on the CPU, depending on the number of threads used in multithreading (2–8 threads). The high performance of MRISIMUL allows its application in large-scale analysis and can bring the computational power of a supercomputer or a large computer cluster to a single GPU personal computer.

Description: Positron emission tomography (PET) images usually suffer from poor signal-to-noise ratio (SNR) due to the high level of noise and low spatial resolution, which adversely affect its performance for lesion detection and quantification. The complementary information present in high-resolution anatomical images from multi-modality imaging systems could potentially be used to improve the ability to detect and/or quantify lesions. However, previous methods that use anatomical priors usually require matched organ/lesion boundaries. In this study, we investigated the use of anatomical information to suppress noise in PET images while preserving both quantitative accuracy and the amplitude of prominent signals that do not have corresponding boundaries on computerized tomography (CT). The proposed approach was realized through a postreconstruction filter based on the nonlocal means (NLM) filter, which reduces noise by computing the weighted average of voxels based on the similarity measurement between patches of voxels within the image. Anatomical knowledge obtained from CT was incorporated to constrain the similarity measurement within a subset of voxels. In contrast to other methods that use anatomical priors, the actual number of neighboring voxels and weights used for smoothing were determined from a robust measurement on PET images within the subset. Thus, the proposed approach can be robust to signal mismatches between PET and CT. A 3-D search scheme was also investigated for the volumetric PET/CT data. The proposed anatomically guided median nonlocal means filter (AMNLM) was first evaluated using a computer phantom and a physical phantom to simulate realistic but challenging situations where small lesions are located in homogeneous regions, which can be detected on PET but not on CT. The proposed method was further assessed with a clinical study of a patient with lung lesions. The performance of the proposed method was compared to Gaussian, edge-preserving - ilateral and NLM filters, as well as median nonlocal means (MNLM) filtering without an anatomical prior. The proposed AMNLM method yielded improved lesion contrast and SNR compared with other methods even with imperfect anatomical knowledge, such as missing lesion boundaries and mismatched organ boundaries.

Description: Preoperative diagnostic magnetic resonance (MR) breast images can provide good contrast between different tissues and 3-D information about suspicious tissues. Aligning preoperative diagnostic MR images with a patient in the theatre during breast conserving surgery could assist surgeons in achieving the complete excision of cancer with sufficient margins. Typically, preoperative diagnostic MR breast images of a patient are obtained in the prone position, while surgery is performed in the supine position. The significant shape change of breasts between these two positions due to gravity loading, external forces and related constraints makes the alignment task extremely difficult. Our previous studies have shown that either nonrigid intensity-based image registration or biomechanical modelling alone are limited in their ability to capture such a large deformation. To tackle this problem, we proposed in this paper a nonlinear biomechanical model-based image registration method with a simultaneous optimization procedure for both the material parameters of breast tissues and the direction of the gravitational force. First, finite element (FE) based biomechanical modelling is used to estimate a physically plausible deformation of the pectoral muscle and the major deformation of breast tissues due to gravity loading. Then, nonrigid intensity-based image registration is employed to recover the remaining deformation that FE analyses do not capture due to the simplifications and approximations of biomechanical models and the uncertainties of external forces and constraints. We assess the registration performance of the proposed method using the target registration error of skin fiducial markers and the Dice similarity coefficient (DSC) of fibroglandular tissues. The registration results on prone and supine MR image pairs are compared with those from two alternative nonrigid registration methods for five breasts. Overall, the proposed algorithm achieved the best registration - erformance on fiducial markers (target registration error, $8.44 pm 5.5$ mm for 45 fiducial markers) and higher overlap rates on segmentation propagation of fibroglandular tissues (DSC value $〉 82%$ ).

Description: Recent theoretical results on low-rank matrix reconstruction have inspired significant interest in low-rank modeling of MRI images. Existing approaches have focused on higher-dimensional scenarios with data available from multiple channels, timepoints, or image contrasts. The present work demonstrates that single-channel, single-contrast, single-timepoint $k$ -space data can also be mapped to low-rank matrices when the image has limited spatial support or slowly varying phase. Based on this, we develop a novel and flexible framework for constrained image reconstruction that uses low-rank matrix modeling of local $k$ -space neighborhoods (LORAKS). A new regularization penalty and corresponding algorithm for promoting low-rank are also introduced. The potential of LORAKS is demonstrated with simulated and experimental data for a range of denoising and sparse-sampling applications. LORAKS is also compared against state-of-the-art methods like homodyne reconstruction, $ell_{1}$ -norm minimization, and total variation minimization, and is demonstrated to have distinct features and advantages. In addition, while calibration-based support and phase constraints are commonly used in existing methods, the LORAKS framework enables calibrationless use of these constraints.

Description: In this paper, we propose a novel method to detect glandular structures in microscopic images of human tissue. We first convert the image from Cartesian space to polar space and then introduce a novel random field model to locate the possible boundary of a gland. Next, we develop a visual feature-based support vector regressor to verify if the detected contour corresponds to a true gland. And finally, we combine the outputs of the random field and the regressor to form the GlandVision algorithm for the detection of glandular structures. Our approach can not only detect the existence of the gland, but also can accurately locate it with pixel accuracy. In the experiments, we treat the task of detecting glandular structures as object (gland) detection and segmentation problems respectively. The results indicate that our new technique outperforms state-of-the-art computer vision algorithms in respective fields.

Description: System designs in single photon emission tomography (SPECT) can be evaluated based on the fundamental trade-off between bias and variance that can be achieved in the reconstruction of emission tomograms. This trade off can be derived analytically using the Cramer-Rao type bounds, which imply the calculation and the inversion of the Fisher information matrix (FIM). The inverse of the FIM expresses the uncertainty associated to the tomogram, enabling the comparison of system designs. However, computing, storing and inverting the FIM is not practical with 3-D imaging systems. In order to tackle the problem of the computational load in calculating the inverse of the FIM, a method based on the calculation of the local impulse response and the variance, in a single point, from a single row of the FIM, has been previously proposed for system design. However this approximation (circulant approximation) does not capture the global interdependence between the variables in shift-variant systems such as SPECT, and cannot account e.g., for data truncation or missing data. Our new formulation relies on subsampling the FIM. The FIM is calculated over a subset of voxels arranged in a grid that covers the whole volume. Every element of the FIM at the grid points is calculated exactly, accounting for the acquisition geometry and for the object. This new formulation reduces the computational complexity in estimating the uncertainty, but nevertheless accounts for the global interdependence between the variables, enabling the exploration of design spaces hindered by the circulant approximation. The graphics processing unit accelerated implementation of the algorithm reduces further the computation times, making the algorithm a good candidate for real-time optimization of adaptive imaging systems. This paper describes the subsampled FIM formulation and implementation details. The advantages and limitations of the new approximation are explored, in comparison with the circulant approxima- ion, in the context of design optimization of a parallel-hole collimator SPECT system and of an adaptive imaging system (similar to the commercially available D-SPECT).

Description: This paper presents the evaluation results of the methods submitted to Challenge US: Biometric Measurements from Fetal Ultrasound Images, a segmentation challenge held at the IEEE International Symposium on Biomedical Imaging 2012. The challenge was set to compare and evaluate current fetal ultrasound image segmentation methods. It consisted of automatically segmenting fetal anatomical structures to measure standard obstetric biometric parameters, from 2D fetal ultrasound images taken on fetuses at different gestational ages (21 weeks, 28 weeks, and 33 weeks) and with varying image quality to reflect data encountered in real clinical environments. Four independent sub-challenges were proposed, according to the objects of interest measured in clinical practice: abdomen, head, femur, and whole fetus. Five teams participated in the head sub-challenge and two teams in the femur sub-challenge, including one team who tackled both. Nobody attempted the abdomen and whole fetus sub-challenges. The challenge goals were two-fold and the participants were asked to submit the segmentation results as well as the measurements derived from the segmented objects. Extensive quantitative (region-based, distance-based, and Bland–Altman measurements) and qualitative evaluation was performed to compare the results from a representative selection of current methods submitted to the challenge. Several experts (three for the head sub-challenge and two for the femur sub-challenge), with different degrees of expertise, manually delineated the objects of interest to define the ground truth used within the evaluation framework. For the head sub-challenge, several groups produced results that could be potentially used in clinical settings, with comparable performance to manual delineations. The femur sub-challenge had inferior performance to the head sub-challenge due to the fact that it is a harder segmentation problem and that the techniques presented relied more on the femur's appea- ance.

Description: A model-based inversion method was used to obtain quantitative estimates of histological parameters from multispectral images of the colon and to examine their potential for discriminating between normal and pathological tissues. Pixel-wise estimates of the mucosal blood volume fraction, density of the scattering particles and thickness were derived using a two-stage method. In the first (forward) stage reflectance spectra corresponding to given instances of the parameter values were computed using Monte Carlo simulation of photon propagation through a multi-layered tissue. In the second (inversion) stage the parameter values were obtained via optimization using an iterated conditional modes algorithm based on Discrete Markov Random Fields. The method was validated on computer generated data contaminated with noise giving a mean normalized root mean square deviation (NRMSD) of 2.04. Validation on ex vivo images demonstrated that parametric maps show gross correspondence with histological features of mucosa characteristic of cancerous, precancerous and noncancerous colon lesions. The key signs of abnormality were shown to be the increase in the blood volume fraction and decrease in the density of scattering particles.

Description: Optoacoustic (photoacoustic) imaging based on cylindrically focused 1-D transducer arrays comes with powerful characteristics in visualizing optical contrast. Parallel reading of multiple detectors arranged around a tissue cross section enables capturing data for generating images of this plane within micro-seconds. Dedicated small animals scanners and handheld systems using 1-D cylindrically focused ultrasound transducer arrays have demonstrated real-time cross-sectional imaging and high in-plane resolution. Yet, the resolution achieved along the axis perpendicular to the focal plane, i.e., the elevation resolution, is determined by the focusing capacities of the detector and is typically lower than the in-plane resolution. Herein, we investigated whether deconvolution of the sensitivity field of the transducer could lead to tangible image improvements. We showcase the findings on experimental measurements from phantoms and animals and discuss the features and the limitations of the approach in improving resolution along the elevation dimension.

Description: One distinguishing property of life is its temporal dynamics, and it is hence only natural that time lapse experiments play a crucial role in modern biomedical research areas such as signaling pathways, drug discovery or developmental biology. Such experiments yield a very large number of images that encode complex cellular activities, and reliable automated cell tracking emerges naturally as a prerequisite for further quantitative analysis. However, many existing cell tracking methods are restricted to using only a small number of features to allow for manual tweaking. In this paper, we propose a novel cell tracking approach that embraces a powerful machine learning technique to optimize the tracking parameters based on user annotated tracks. Our approach replaces the tedious parameter tuning with parameter learning and allows for the use of a much richer set of complex tracking features, which in turn affords superior prediction accuracy. Furthermore, we developed an active learning approach for efficient training data retrieval, which reduces the annotation effort to only 17%. In practical terms, our approach allows life science researchers to inject their expertise in a more intuitive and direct manner. This process is further facilitated by using a glyph visualization technique for ground truth annotation and validation. Evaluation and comparison on several publicly available benchmark sequences show significant performance improvement over recently reported approaches. Code and software tools are provided to the public.

Description: In this paper, we propose a novel approach to landmark-based shape representation that is based on transportation theory, where landmarks are considered as sources and destinations, all possible landmark connections as roads, and established landmark connections as goods transported via these roads. Landmark connections, which are selectively established, are identified through their statistical properties describing the shape of the object of interest, and indicate the least costly roads for transporting goods from sources to destinations. From such a perspective, we introduce three novel shape representations that are combined with an existing landmark detection algorithm based on game theory. To reduce computational complexity, which results from the extension from 2-D to 3-D segmentation, landmark detection is augmented by a concept known in game theory as strategy dominance. The novel shape representations, game-theoretic landmark detection and strategy dominance are combined into a segmentation framework that was evaluated on 3-D computed tomography images of lumbar vertebrae and femoral heads. The best shape representation yielded symmetric surface distance of 0.75 mm and 1.11 mm, and Dice coefficient of 93.6% and 96.2% for lumbar vertebrae and femoral heads, respectively. By applying strategy dominance, the computational costs were further reduced for up to three times.

Description: Proton radiography has applications in patient alignment and verification procedures for proton beam radiation therapy. In this paper, we report an experiment which used 200 MeV protons to generate proton energy-loss and scattering radiographs of a hand phantom. The experiment used the first-generation proton computed tomography (CT) scanner prototype, which was installed on the research beam line of the clinical proton synchrotron at Loma Linda University Medical Center. It was found that while both radiographs displayed anatomical details of the hand phantom, the energy-loss radiograph had a noticeably higher resolution. Nonetheless, scattering radiography may yield more contrast between soft and bone tissue than energy-loss radiography, however, this requires further study. This study contributes to the optimization of the performance of the next-generation of clinical proton CT scanners. Furthermore, it demonstrates the potential of proton imaging (proton radiography and CT), which is now within reach of becoming available as a new, potentially low-dose medical imaging modality.

Description: This paper presents a predictive framework for the statistical personalization of ventricular fibers. To this end, the relationship between subject-specific geometry of the left (LV) and right ventricles (RV) and fiber orientation is learned statistically from a training sample of ex vivo diffusion tensor imaging datasets. More specifically, the axes in the shape space which correlate most with the myocardial fiber orientations are extracted and used for prediction in new subjects. With this approach and unlike existing fiber models, inter-subject variability is taken into account to generate latent shape predictors that are statistically optimal to estimate fiber orientation at each individual myocardial location. The proposed predictive model was applied to the task of personalizing fibers in 10 canine subjects. The results indicate that the ventricular shapes are good predictors of fiber orientation, with an improvement of 11.4% in accuracy over the average fiber model.

Description: The image reconstruction problem encountered in diffuse optical tomographic imaging is ill-posed in nature, necessitating the usage of regularization to result in stable solutions. This regularization also results in loss of resolution in the reconstructed images. A frame work, that is attributed by model-resolution, to improve the reconstructed image characteristics using the basis pursuit deconvolution method is proposed here. The proposed method performs this deconvolution as an additional step in the image reconstruction scheme. It is shown, both in numerical and experimental gelatin phantom cases, that the proposed method yields better recovery of the target shapes compared to traditional method, without the loss of quantitativeness of the results.

Description: Noninvasive imaging of cardiac electrical function has begun to move towards clinical adoption. Here, we consider one common formulation of the problem, in which the goal is to estimate the spatial distribution of electrical activation times during a cardiac cycle. We address the challenge of understanding the robustness and uncertainty of solutions to this formulation. This formulation poses a nonconvex, nonlinear least squares optimization problem. We show that it can be relaxed to be convex, at the cost of some degree of physiological realism of the solution set, and that this relaxation can be used as a framework to study model inaccuracy and solution uncertainty. We present two examples, one using data from a healthy human subject and the other synthesized with the ECGSIM software package. In the first case, we consider uncertainty in the initial guess and regularization parameter. In the second case, we mimic the presence of an ischemic zone in the heart in a way which violates a model assumption. We show that the convex relaxation allows understanding of spatial distribution of parameter sensitivity in the first case, and identification of model violation in the second.

Description: Multi-voxel pattern analysis (MVPA) of functional magnetic resonance imaging (fMRI) data is an emerging approach for probing the neural correlates of cognition. MVPA allows cognitive states to be modeled as distributed patterns of neural activity and classified according to stimulus conditions. In practice, building a robust, generalizable classification model can be challenging because the number of voxels (features) far exceeds the number of stimulus instances/data observations. To avoid model overfitting, there is a need to select informative voxels before building a classification model. In this paper, we propose a robust feature (voxel) selection framework using mutual information (MI) and partial least square regression (PLS) to establish an informativeness index for prioritizing selection of voxels based on the degree of their association to the experimental conditions. We evaluated the robustness of our proposed framework by assessing performance of standard classification algorithms, when combined with our feature selection approach, in a publicly-available fMRI dataset of object-level representation widely used to benchmark MVPA performance (Haxby, 2001). The computational results suggest that our feature selection framework based on MI and PLS drastically improves the classification accuracy relative to those previously reported in the literature. Our results also suggest that highly informative voxels may provide meaningful insight into the functional-anatomic relationship of brain activity and stimulus conditions.

Description: Dynamic imaging data are currently analyzed with a tracer-kinetic theory developed for individual time curves measured over whole organs. The assumption is that voxels represent isolated systems which all receive indicator through the same arterial inlet. This leads to well-known systematic errors, but also fails to exploit the spatial structure of the data. In this study, a more general theoretical framework is developed which makes full use of the specific structure of image data. The theory encodes the fact that voxels receive indicator from their immediate neighbors rather than from an upstream arterial input. This results in a tracer-kinetic field theory where the tissue parameters are functions of space which can be measured by analyzing the temporal and spatial patterns in the concentrations. The implications are evaluated through a number of field models for common tissue types. The key benefits of a tracer-kinetic field theory are that: 1) long-standing systematic errors can be corrected, specifically the issue of bolus dispersion and the contamination of large-vessel blood flow on tissue perfusion measurements; 2) additional tissue parameters can be measured that characterize convective or diffusive exchange between voxels; 3) the need to measure a separate arterial input function can be eliminated.

Description: We propose a novel global optimization-based approach to segmentation of 3-D prostate transrectal ultrasound (TRUS) and T2 weighted magnetic resonance (MR) images, enforcing inherent axial symmetry of prostate shapes to simultaneously adjust a series of 2-D slice-wise segmentations in a “global” 3-D sense. We show that the introduced challenging combinatorial optimization problem can be solved globally and exactly by means of convex relaxation. In this regard, we propose a novel coherent continuous max-flow model (CCMFM), which derives a new and efficient duality-based algorithm, leading to a GPU-based implementation to achieve high computational speeds. Experiments with 25 3-D TRUS images and 30 3-D T2w MR images from our dataset, and 50 3-D T2w MR images from a public dataset, demonstrate that the proposed approach can segment a 3-D prostate TRUS/MR image within 5–6 s including 4–5 s for initialization, yielding a mean Dice similarity coefficient of $93.2%pm 2.0%$ for 3-D TRUS images and $88.5%pm 3.5%$ for 3-D MR images. The proposed method also yields relatively low intra- and inter-observer variability introduced by user manual initialization, suggesting a high reproducibility, independent of observers.

Description: Zero-echo Time (ZTE) imaging is a promising technique for magnetic resonance imaging (MRI) of short-T2 tissue nuclei in tissues. A problem inherent to the method currently hindering its translation to the clinic is the presence of a spatial encoding gradient during excitation, which causes the hard pulse to become spatially selective, resulting in blurring and shadow artifacts in the image. While shortening radio-frequency (RF) pulse duration alleviates this problem the resulting elevated RF peak power and specific absorption rate (SAR) in practice impede such a solution. In this work, an approach is described to correct the artifacts by applying quadratic phase-modulated RF excitation and iteratively solving an inverse problem formulated from the signal model of ZTE imaging. A simple pulse sequence is also developed to measure the excitation profile of the RF pulse. Results from simulations, phantom and in vivo studies, demonstrate the effectiveness of the method in correcting image artifacts caused by inhomogeneous excitation. The proposed method may contribute toward establishing ZTE MRI as a routine 3D pulse sequence for imaging protons and other nuclei with quasi solid-state behavior on clinical scanners.

Description: This paper proposes the application of multiple importance sampling in fully 3-D positron emission tomography to speed up the iterative reconstruction process. The proposed method combines the results of lines of responses (LOR) driven and voxel driven projections keeping their advantages, like importance sampling, performance and parallel execution on graphics processing units. Voxel driven methods can focus on point like features while LOR driven approaches are efficient in reconstructing homogeneous regions. The theoretical basis of the combination is the application of the mixture of the samples generated by the individual importance sampling methods, emphasizing a particular method where it is better than others. The proposed algorithms are built into the Tera-tomo system.

Description: The reconstruction of dynamic magnetic resonance data from an undersampled ${rm k}$ -space has been shown to have a huge potential in accelerating the acquisition process of this imaging modality. With the introduction of compressed sensing (CS) theory, solutions for undersampled data have arisen which reconstruct images consistent with the acquired samples and compliant with a sparsity model in some transform domain. Fixed basis transforms have been extensively used as sparsifying transforms in the past, but recent developments in dictionary learning (DL) have been shown to outperform them by training an overcomplete basis that is optimal for a particular dataset. We present here an iterative algorithm that enables the application of DL for the reconstruction of cardiac cine data with Cartesian undersampling. This is achieved with local processing of spatio-temporal 3D patches and by independent treatment of the real and imaginary parts of the dataset. The enforcement of temporal gradients is also proposed as an additional constraint that can greatly accelerate the convergence rate and improve the reconstruction for high acceleration rates. The method is compared to and shown to systematically outperform ${rm k}$ - ${rm t}$ FOCUSS, a successful CS method that uses a fixed basis transform.

Description: We propose a patient-specific model based on partial differential equation to predict the evolution of lung tumors during radiotherapy. The evolution of tumor cell density is formulated by three terms: 1) advection describing the advective flux transport of tumor cells, 2) proliferation representing the tumor cell proliferation modeled as Gompertz differential equation, and 3) treatment quantifying the radiotherapeutic efficacy from linear quadratic formulation. We consider that tumor cell density variation can be derived from positron emission tomography images, the novel idea is to model the advection term by calculating 3D optical flow field from sequential images. To estimate patient-specific parameters, we propose an optimization between the predicted and observed images, under a global constraint that the tumor volume decreases exponentially as radiation dose increases. A thresholding on the predicted tumor cell densities is then used to define tumor contours, tumor volumes and maximum standardized uptake values (SUVmax). Results obtained on seven patients show a satisfying agreement between the predicted tumor contours and those drawn by an expert.

Description: The goal of this work is to reliably and accurately localize anatomical landmarks in 3-D computed tomography scans, particularly for the deformable registration of whole-body scans, which show huge variation in posture, and the spatial distribution of anatomical features. Parts-based graphical models (GM) have shown attractive properties for this task because they capture naturally anatomical relationships between landmarks. Unfortunately, standard GMs are learned from manually annotated training images and the quantity of landmarks is limited by the high cost of expert annotation. We propose a novel method that automatically learns new corresponding landmarks from a database of 3-D whole-body CT scans, using a limited initial set of expert-labeled ground-truth landmarks. The newly learned landmarks, called B-landmarks, are used to build enriched GMs. We compare our method of deformable registration based on such GM landmarks to a conventional deformable registration method and to a “baseline” state-of-the-art GM. The results show our method finds new relevant anatomical correspondences and improves by up to 35% the matching accuracy of highly variable skeletal and soft-tissue landmarks of clinical interest.

Description: We investigate using dual time-point PET data to perform Patlak modeling. This approach can be used for whole body dynamic PET studies in which we compute voxel-wise estimates of Patlak parameters using two frames of data for each bed position. Our approach directly uses list-mode arrival times for each event to estimate the Patlak parametric image. We use a penalized likelihood method in which the penalty function uses spatially variant weighting to ensure a count independent local impulse response. We evaluate performance of the method in comparison to fractional changes in SUV values (%DSUV) between the two frames using Cramer Rao analysis and Monte Carlo simulation. Receiver operating characteristic (ROC) curves are used to compare performance in differentiating tumors relative to background based on the dynamic data sets. Using area under the ROC curve as a performance metric, we show superior performance of Patlak relative to %DSUV over a range of dynamic data sets and parameters. These results suggest that Patlak analysis may be appropriate for analysis of dual time-point whole body PET data and could lead to superior detection of tumors relative to %DSUV metrics.

Description: Harmonic motion imaging for focused ultrasound (HMIFU) utilizes an amplitude-modulated HIFU beam to induce a localized focal oscillatory motion simultaneously estimated. The objective of this study is to develop and show the feasibility of a novel fast beamforming algorithm for image reconstruction using GPU-based sparse-matrix operation with real-time feedback. In this study, the algorithm was implemented onto a fully integrated, clinically relevant HMIFU system. A single divergent transmit beam was used while fast beamforming was implemented using a GPU-based delay-and-sum method and a sparse-matrix operation. Axial HMI displacements were then estimated from the RF signals using a 1-D normalized cross-correlation method and streamed to a graphic user interface with frame rates up to 15 Hz, a 100-fold increase compared to conventional CPU-based processing. The real-time feedback rate does not require interrupting the HIFU treatment. Results in phantom experiments showed reproducible HMI images and monitoring of 22 in vitro HIFU treatments using the new 2-D system demonstrated reproducible displacement imaging, and monitoring of 22 in vitro HIFU treatments using the new 2-D system showed a consistent average focal displacement decrease of 46.7 $pm$ 14.6% during lesion formation. Complementary focal temperature monitoring also indicated an average rate of displacement increase and decrease with focal temperature at $0.84pm 1.15%/^{circ}{rm C}$ , and $2.03pm 0.93%/^{circ}{rm C}$ , respectively. These results reinforce the HMIFU capability of estimating and monitoring stiffness related changes in real time. Current ongoing studies include clinical translation of the presented system for monitoring of HIFU treatment for breast and pancreatic tumor applications.

Description: Shear wave speed can be used to assess tissue elasticity, which is associated with tissue health. Ultrasound shear wave elastography techniques based on measuring the propagation speed of the shear waves induced by acoustic radiation force are becoming promising alternatives to biopsy in liver fibrosis staging. However, shear waves generated by such methods are typically very weak. Therefore, the penetration may become problematic, especially for overweight or obese patients. In this study, we developed a new method called external vibration multi-directional ultrasound shearwave elastography (EVMUSE), in which external vibration from a loudspeaker was used to generate a multi-directional shear wave field. A directional filter was then applied to separate the complex shear wave field into several shear wave fields propagating in different directions. A 2-D shear wave speed map was reconstructed from each individual shear wave field, and a final 2-D shear wave speed map was constructed by compounding these individual wave speed maps. The method was validated using two homogeneous phantoms and one multi-purpose tissue-mimicking phantom. Ten patients undergoing liver magnetic resonance elastography (MRE) were also studied with EVMUSE to compare results between the two methods. Phantom results showed EVMUSE was able to quantify tissue elasticity accurately with good penetration. In vivo EVMUSE results were well correlated with MRE results, indicating the promise of using EVMUSE for liver fibrosis staging.

Description: This study demonstrates real-time scan plane control dependent on three-dimensional needle bending, as measured from magnetic resonance imaging (MRI)-compatible optical strain sensors. A biopsy needle with embedded fiber Bragg grating (FBG) sensors to measure surface strains is used to estimate its full 3-D shape and control the imaging plane of an MR scanner in real-time, based on the needle's estimated profile. The needle and scanner coordinate frames are registered to each other via miniature radio-frequency (RF) tracking coils, and the scan planes autonomously track the needle as it is deflected, keeping its tip in view. A 3-D needle annotation is superimposed over MR-images presented in a 3-D environment with the scanner's frame of reference. Scan planes calculated based on the FBG sensors successfully follow the tip of the needle. Experiments using the FBG sensors and RF coils to track the needle shape and location in real-time had an average root mean square error of 4.2 mm when comparing the estimated shape to the needle profile as seen in high resolution MR images. This positional variance is less than the image artifact caused by the needle in high resolution SPGR (spoiled gradient recalled) images. Optical fiber strain sensors can estimate a needle's profile in real-time and be used for MRI scan plane control to potentially enable faster and more accurate physician response.

Description: One of the technical challenges in cine magnetic resonance imaging (MRI) is to reduce the acquisition time to enable the high spatio-temporal resolution imaging of a cardiac volume within a short scan time. Recently, compressed sensing approaches have been investigated extensively for highly accelerated cine MRI by exploiting transform domain sparsity using linear transforms such as wavelets, and Fourier. However, in cardiac cine imaging, the cardiac volume changes significantly between frames, and there often exist abrupt pixel value changes along time. In order to effectively sparsify such temporal variations, it is necessary to exploit temporal redundancy along motion trajectories. This paper introduces a novel patch-based reconstruction method to exploit geometric similarities in the spatio-temporal domain. In particular, we use a low rank constraint for similar patches along motion, based on the observation that rank structures are relatively less sensitive to global intensity changes, but make it easier to capture moving edges. A Nash equilibrium formulation with relaxation is employed to guarantee convergence. Experimental results show that the proposed algorithm clearly reconstructs important anatomical structures in cardiac cine image and provides improved image quality compared to existing state-of-the-art methods such as k-t FOCUSS, k-t SLR, and MASTeR.

Description: Recently inexpensive graphical processing units (GPUs) have become established as a viable alternative to traditional CPUs for many medical image processing applications. GPUs offer the potential of very significant improvements in performance at low cost and with low power consumption. One way in which GPU programs differ from traditional CPU programs is that increasingly elaborate calculations per voxel may not impact of the overall processing time because memory accesses can dominate execution time. This paper presents a new GPU based elastic image registration program named Ezys. The Ezys image registration algorithm belongs to the wide class of diffeomorphic demons but uses surface preserving image smoothing and regularization filters designed for a GPU that would be computationally expensive on a CPU. We describe the methods used in Ezys and present results from two important neuroscience applications. Firstly inter-subject registration for transfer of anatomical labels and secondly longitudinal intra-subject registration to quantify atrophy in individual subjects. Both experiments showed that Ezys registration compares favorably with other popular elastic image registration programs. We believe Ezys is a useful tool for neuroscience and other applications, and also demonstrates the value of developing of novel image processing filters specifically designed for GPUs.

Description: In ophthalmic ultrasonography the crystalline lens is known to be the main source of phase aberration, causing a significant decrease in resolution and distortion effects on axial B-scans. This paper proposes a computationally efficient method to correct the phase aberration arising from the crystalline lens, including refraction effects using a bending ray tracing approach based on Fermat's principle. This method is used as a basis to perform eye-adapted beamforming (BF), with appropriate focusing delays for a 128-element 20-MHz linear array in both emission and reception. Implementation was achieved on an in-house developed experimental ultrasound scanning device, the ECODERM. The proposed BF was tested in vitro by imaging a wire phantom through an eye phantom consisting of a synthetic gelatin lens anatomically set up in an appropriate liquid (turpentine) to approach the in vivo velocity ratio. Both extremes of accommodation shapes of the human crystalline lens were investigated. The performance of the developed BF was evaluated in relation to that in homogeneous medium and compared to a conventional delay-and-sum (DAS) BF and a second adapted BF which was simplified to ignore the lens refraction. Global expectations provided by our method with the transducer array are reviewed by an analysis quantifying both image quality and spatial fidelity, as well as the detrimental effects of a crystalline lens in conventional reconstruction. Compared to conventional array imaging, the results indicated a two-fold improvement in the lateral resolution, greater sensitivity and a considerable reduction of spatial distortions that were sufficient to envisage reliable biometry directly in B-mode, especially phakometry.

Description: Nonparametric hemodynamic response function (HRF) estimation in functional near-infrared spectroscopy (fNIRS) data plays an important role when investigating the temporal dynamics of a brain region response during activations. Assuming the drift arising from both physical and physiological effects in fNIRS data is Lipschitz continuous; a novel algorithm for joint HRF and drift estimation is derived in this paper. The proposed algorithm estimates the HRF by applying a first-order differencing to the fNIRS time series samples in order to remove the drift effect. An estimate of the drift is then obtained using a wavelet thresholding technique applied to the residuals generated by removing the estimated induced activation response from the fNIRS time-series. It is shown that the proposed HRF estimator is $sqrt {N}$ consistent whereas the estimator of the drift is asymptotically optimal. The de-drifted fNIRS oxygenated (HbO) and deoxygenated (HbR) hemoglobin responses are then obtained by removing the corresponding estimated drifts from the fNIRS time-series. Its performance is assessed using both simulated and real fNIRS data sets. The application results reveal that the proposed joint HRF and drift estimation method is efficient both computationally and in terms of accuracy. In comparison to traditional model based methods used for HRF estimation, the proposed novel method avoids the selection of a model to remove the drift component. As a result, the proposed method finds an optimal estimate of the fNIRS drift and offers a model-free approach to de-drift the HbO/HbR responses.

Description: Most conventional spine interventions are performed under X-ray fluoroscopy guidance. In recent years, there has been a growing interest to develop nonionizing imaging alternatives to guide these procedures. Ultrasound guidance has emerged as a leading alternative. However, a challenging problem is automatic identification of the spinal anatomy in ultrasound data. In this paper, we propose a local phase-based bone feature enhancement technique that can robustly identify the spine surface in ultrasound images. The local phase information is obtained using a gradient energy tensor filter. This information is used to construct local phase tensors in ultrasound images, which highlight the spine surface. We show that our proposed approach results in a more distinct enhancement of the bone surfaces compared to recently proposed techniques based on monogenic scale-space filters and logarithmic Gabor filters. We also demonstrate that registration accuracy of a statistical shape+pose model of the spine to 3-D ultrasound images can be significantly improved, using the proposed method, compared to those obtained using monogenic scale-space filters and logarithmic Gabor filters.

Description: We propose a technique for myocardial motion estimation based on image registration using both B-mode echocardiographic images and tissue Doppler sequences acquired interleaved. The velocity field is modeled continuously using B-splines and the spatiotemporal transform is constrained to be diffeomorphic. Images before scan conversion are used to improve the accuracy of the estimation. The similarity measure includes a model of the speckle pattern distribution of B-mode images. It also penalizes the disagreement between tissue Doppler velocities and the estimated velocity field. Registration accuracy is evaluated and compared to other alternatives using a realistic synthetic dataset, obtaining mean displacement errors of about 1 mm. Finally, the method is demonstrated on data acquired from six volunteers, both at rest and during exercise. Robustness is tested against low image quality and fast heart rates during exercise. Results show that our method provides a robust motion estimate in these situations.

Description: Noninvasive functional imaging of awake, unrestrained small animals using motion-compensation removes the need for anesthetics and enables an animal's behavioral response to stimuli or administered drugs to be studied concurrently with imaging. While the feasibility of motion-compensated radiotracer imaging of awake rodents using marker-based optical motion tracking has been shown, markerless motion tracking would avoid the risk of marker detachment, streamline the experimental workflow, and potentially provide more accurate pose estimates over a greater range of motion. We have developed a stereoscopic tracking system which relies on native features on the head to estimate motion. Features are detected and matched across multiple camera views to accumulate a database of head landmarks and pose is estimated based on 3D-2D registration of the landmarks to features in each image. Pose estimates of a taxidermal rat head phantom undergoing realistic rat head motion via robot control had a root mean square error of 0.15 and 1.8 mm using markerless and marker-based motion tracking, respectively. Markerless motion tracking also led to an appreciable reduction in motion artifacts in motion-compensated positron emission tomography imaging of a live, unanesthetized rat. The results suggest that further improvements in live subjects are likely if nonrigid features are discriminated robustly and excluded from the pose estimation process.

Description: The paper presents a computational model of magnetic resonance (MR) flow imaging. The model consists of three components. The first component is used to generate complex vascular structures, while the second one provides blood flow characteristics in the generated vascular structures by the lattice Boltzmann method. The third component makes use of the generated vascular structures and flow characteristics to simulate MR flow imaging. To meet computational demands, parallel algorithms are applied in all the components. The proposed approach is verified in three stages. In the first stage, experimental validation is performed by an in vitro phantom. Then, the simulation possibilities of the model are shown. Flow and MR flow imaging in complex vascular structures are presented and evaluated. Finally, the computational performance is tested. Results show that the model is able to reproduce flow behavior in large vascular networks in a relatively short time. Moreover, simulated MR flow images are in accordance with the theoretical considerations and experimental images. The proposed approach is the first such an integrative solution in literature. Moreover, compared to previous works on flow and MR flow imaging, this approach distinguishes itself by its computational efficiency. Such a connection of anatomy, physiology and image formation in a single computer tool could provide an in silico solution to improving our understanding of the processes involved, either considered together or separately.

Description: Magnetic resonance parameter mapping (e.g., $T_{1}$ mapping, $T_{2}$ mapping, $T^{ast }_{2}$ mapping) is a valuable tool for tissue characterization. However, its practical utility has been limited due to long data acquisition time. This paper addresses this problem with a new model-based parameter mapping method. The proposed method utilizes a formulation that integrates the explicit signal model with sparsity constraints on the model parameters, enabling direct estimation of the parameters of interest from highly undersampled, noisy $ {bf k}$ -space data. An efficient greedy-pursuit algorithm is described to solve the resulting constrained parameter estimation problem. Estimation-theoretic bounds are also derived to analyze the benefits of incorporating sparsity constraints and benchmark the performance of the proposed method. The theoretical properties and empirical performance of the proposed method are illustrated in a $T_{2}$ mapping application example using computer simulations.

Description: —Incorporating prior knowledge into image segmentation algorithms has proven useful for obtaining more accurate and plausible results. Two important constraints, containment and exclusion of regions, have gained attention in recent years mainly due to their descriptive power. In this paper, we augment the level set framework with the ability to handle these two intuitive geometric relationships, containment and exclusion, along with a distance constraint between boundaries of multi-region objects. Level set’s important property of automatically handling topological changes of evolving contours/surfaces enables us to segment spatially-recurring objects (e.g., multiple instances of multi-region cells in a large microscopy image) while satisfying the two aforementioned constraints. In addition, the level set approach gives us a very simple and natural way to compute the distance between contours/surfaces and impose constraints on it. The downside, however, is a local optimization framework in which the final segmentation solution depends on the initialization. In fact, here, we sacrifice the optimizability (local instead of global solution) in exchange for lower space complexity (less memory usage) and faster runtime (especially for large microscopic images) as well as no grid artifacts. Nevertheless, the result from validating our method on several biomedical applications showed the utility and advantages of this augmented level set framework (even with rough initialization that is distant from the desired boundaries). We also compared our framework with its counterpart methods in the discrete domain and reported the pros and cons of each of these methods in terms of metrication error and efficiency in memory usage and runtime.

Description: While tomographic imaging of cardiac structure and kinetics has improved substantially, electrophysiological mapping of the heart is still restricted to the surface with little or no depth information beneath. The progress in reconstructing 3-D action potential from surface voltage data has been hindered by the intrinsic ill-posedness of the problem and the lack of a unique solution in the absence of prior assumptions. In this work, we propose a novel adaption of the total-variation (TV) prior to exploit the unique spatial property of transmural action potential of being piecewise smooth with a steep boundary (gradient) separating depolarized and repolarized regions. We present a variational TV-prior instead of a common discrete TV-prior for improved robustness to mesh resolution, and solve the TV-minimization by a sequence of weighted, first-order L2-norm minimization. In a large set of phantom experiments, the proposed method is shown to outperform existing quadratic methods in preserving the steep gradient of action potential along the border of infarcts, as well as in capturing the disruption to the normal path of electrical wavefronts. Real-data experiments also further demonstrate the potential of the proposed method in revealing the location and shape of infarcts when quadratic methods fail to do so.

Description: Flow diversion is an emerging endovascular treatment option for cerebral aneurysms. Quantitative assessment of hemodynamic changes induced by flow diversion can aid clinical decision making in the treatment of cerebral aneurysms. In this article, besides summarizing past key research efforts, we propose a novel metric for the angiographic assessment of flow diverter deployments in the treatment of cerebral aneurysms. By analyzing the frequency spectra of signals derived from digital subtraction angiography (DSA) series, the metric aims to quantify the prevalence of frequency components that correspond to the patient-specific heart rate. Indicating the decoupling of aneurysms from healthy blood circulation, our proposed metric could advance clinical guidelines for treatment success prediction. The very promising results of a retrospective feasibility study on 26 DSA series warrant future efforts to study the validity of the proposed metric within a clinical setting.

Description: Sources of uncertainty in the boundaries of structures in medical images have motivated the use of probabilistic labels in segmentation applications. An important component in many medical image segmentation tasks is the use of a shape model, often generated by applying statistical techniques to training data. Standard statistical techniques (e.g., principal component analysis) often assume data lies in an unconstrained vector space, but probabilistic labels are constrained to the unit simplex. If these statistical techniques are used directly on probabilistic labels, relative uncertainty information can be sacrificed. A standard method for facilitating analysis of probabilistic labels is to map them to a vector space using the LogOdds transform. However, the LogOdds transform is asymmetric in one of the labels, which skews results in some applications. The isometric log-ratio (ILR) transform is a symmetrized version of the LogOdds transform, and is so named as it is an isometry between the Aitchison geometry, the inherent geometry of the simplex, and standard Euclidean geometry. We explore how to interpret the Aitchison geometry when applied to probabilistic labels in medical image segmentation applications. We demonstrate the differences when applying the LogOdds transform or the ILR transform to probabilistic labels prior to statistical analysis. Specifically, we show that statistical analysis of ILR transformed data better captures the variability of anatomical shapes in cases where multiple different foreground regions share boundaries (as opposed to foreground-background boundaries).

Description: An algorithm dedicated to automatic segmentation of breast magnetic resonance images is presented in this paper. Our approach is based on a pipeline that includes a denoising step and statistical segmentation. The noise removal preprocessing relies on an anisotropic diffusion scheme, whereas the statistical segmentation is conducted through a Markov random field model. The continuous updating of all parameters governing the diffusion process enables automatic denoising, and the partial volume effect is also addressed during the labeling step. To assess the relevance, the Jaccard similarity coefficient was computed. Experiments were conducted on synthetic data and breast magnetic resonance images extracted from a high-risk population. The relevance of the approach for the dataset is highlighted, and we demonstrate accuracy superior to that of traditional clustering algorithms. The results emphasize the benefits of both denoising guided by input data and the inclusion of spatial dependency through a Markov random field. For example, the Jaccard coefficient for the clinical data was increased by 114%, 109%, and 140% with respect to a K-means algorithm and, respectively, for the adipose, glandular and muscle and skin components. Moreover, the agreement between the manual segmentations provided by an experienced radiologist and the automatic segmentations performed with this algorithm was good, with Jaccard coefficients equal to 0.769, 0.756, and 0.694 for the above-mentioned classes.

Description: Pelvic floor dysfunction is common in women after childbirth and precise segmentation of magnetic resonance images (MRI) of the pelvic floor may facilitate diagnosis and treatment of patients. However, because of the complexity of its structures, manual segmentation of the pelvic floor is challenging and suffers from high inter and intra-rater variability of expert raters. Multiple template fusion algorithms are promising segmentation techniques for these types of applications, but they have been limited by imperfections in the alignment of templates to the target, and by template segmentation errors. A number of algorithms sought to improve segmentation performance by combining image intensities and template labels as two independent sources of information, carrying out fusion through local intensity weighted voting schemes. This class of approach is a form of linear opinion pooling, and achieves unsatisfactory performance for this application. We hypothesized that better decision fusion could be achieved by assessing the contribution of each template in comparison to a reference standard segmentation of the target image and developed a novel segmentation algorithm to enable automatic segmentation of MRI of the female pelvic floor. The algorithm achieves high performance by estimating and compensating for both imperfect registration of the templates to the target image and template segmentation inaccuracies. A local image similarity measure is used to infer a local reliability weight, which contributes to the fusion through a novel logarithmic opinion pooling. We evaluated our new algorithm in comparison to nine state-of-the-art segmentation methods and demonstrated our algorithm achieves the highest performance.

Description: Provides a listing of current staff, committee members and society officers.

Description: A model-based reconstruction technique for accelerated T2 mapping with improved accuracy is proposed using undersampled Cartesian spin-echo magnetic resonance imaging (MRI) data. The technique employs an advanced signal model for T2 relaxation that accounts for contributions from indirect echoes in a train of multiple spin echoes. An iterative solution of the nonlinear inverse reconstruction problem directly estimates spin-density and T2 maps from undersampled raw data. The algorithm is validated for simulated data as well as phantom and human brain MRI at 3T. The performance of the advanced model is compared to conventional pixel-based fitting of echo-time images from fully sampled data. The proposed method yields more accurate T2 values than the mono-exponential model and allows for retrospective undersampling factors of at least 6. Although limitations are observed for very long T2 relaxation times, respective reconstruction problems may be overcome by a gradient dampening approach. The analytical gradient of the utilized cost function is included as Appendix . The source code is made available to the community.

Description: Tissue image grading is one of the most important steps in prostate cancer diagnosis, where the pathologist relies on the gland structure to assign a Gleason grade to the tissue image. In this grading scheme, the discrimination between grade 3 and grade 4 is the most difficult, and receives the most attention from researchers. In this study, we propose a novel method (called nuclei-based method) that 1) utilizes graph theory techniques to segment glands and 2) computes a gland-score (based on the spatial arrangement of nuclei) to estimate how similar a segmented region is to a gland. Next, we create a fusion method by combining this nuclei-based method with the lumen-based method presented in our previous work to improve the performance of grade 3 versus grade 4 classification problem (the accuracy is now improved to 87.3% compared to 81.1% of the lumen-based method alone). To segment glands, we build a graph of nuclei and lumina in the image, and use the normalized cut method to partition the graph into different components, each corresponding to a gland. Unlike most state-of-the-art lumen-based gland segmentation method, the nuclei-based method is able to segment glands without lumen or glands with multiple lumina. Moreover, another important contribution in this research is the development of a set of measures to exploit the difference in nuclei spatial arrangement between grade 3 images (where nuclei form closed chain structure on the gland boundary) and grade 4 image (where nuclei distribute more randomly in the gland). These measures are combined to generate a single gland-score value, which estimates how similar a segmented region (which is a set of nuclei and lumina) is to a gland.

Description: The Sussex MK4 electrical impedance mammography system is especially designed for 3-D breast screening. It aims to diagnose breast cancer at an early stage when it is most treatable. Planar electrodes are employed in this system. The challenge with planar electrodes is the inaccuracy and poor sensitivity in the vertical direction for 3-D imaging. An enhanced image reconstruction algorithm using a duo-mesh method is proposed to improve the vertical accuracy and sensitivity. The novel part of the enhanced image reconstruction algorithm is the correction term. To evaluate the new algorithm, an image processing based error analysis method is presented, which not only can precisely assess the error of the reconstructed image but also locate the center and outline the center and outline the shape of the objects of interest. Although the enhanced image reconstruction algorithm and the image processing based error analysis method are designed for the Sussex MK4 system, they are applicable to all electrical impedance tomography systems, regardless of the hardware design. To validate the enhanced algorithm, performance results from simulations, phantoms and patients are presented.

Description: Attenuation correction is an essential requirement for quantification of positron emission tomography (PET) data. In PET/CT acquisition systems, attenuation maps are derived from computed tomography (CT) images. However, in hybrid PET/MR scanners, magnetic resonance imaging (MRI) images do not directly provide a patient-specific attenuation map. The aim of the proposed work is to improve attenuation correction for PET/MR scanners by generating synthetic CTs and attenuation maps. The synthetic images are generated through a multi-atlas information propagation scheme, locally matching the MRI-derived patient's morphology to a database of MRI/CT pairs, using a local image similarity measure. Results show significant improvements in CT synthesis and PET reconstruction accuracy when compared to a segmentation method using an ultrashort-echo-time MRI sequence and to a simplified atlas-based method.

Description: This page or pages intentionally left blank.

Description: Magnetic resonance (MR) fingerprinting is a technique for acquiring and processing MR data that simultaneously provides quantitative maps of different tissue parameters through a pattern recognition algorithm. A predefined dictionary models the possible signal evolutions simulated using the Bloch equations with different combinations of various MR parameters and pattern recognition is completed by computing the inner product between the observed signal and each of the predicted signals within the dictionary. Though this matching algorithm has been shown to accurately predict the MR parameters of interest, one desires a more efficient method to obtain the quantitative images. We propose to compress the dictionary using the singular value decomposition, which will provide a low-rank approximation. By compressing the size of the dictionary in the time domain, we are able to speed up the pattern recognition algorithm, by a factor of between 3.4–4.8, without sacrificing the high signal-to-noise ratio of the original scheme presented previously.

Description: We propose a novel, physics-based method for detecting multi-scale tubular features in ultrasound images. The detector is based on a Hessian-matrix eigenvalue method, but unlike previous work, our detector is guided by an optimal model of vessel-like structures with respect to the ultrasound-image formation process. Our method provides a voxel-wise probability map, along with estimates of the radii and orientations of the detected tubes. These results can then be used for further processing, including segmentation and enhanced volume visualization. Most Hessian-based algorithms, including the well-known Frangi filter, were developed for CTA or MRA; they implicitly assume symmetry about the vessel centerline. This is not consistent with ultrasound data. We overcome this limitation by introducing a novel filter that allows multi-scale estimation both with respect to the vessel's centerline and with respect to the vessel's border. We use manually-segmented ultrasound imagery from 35 patients to show that our method is superior to standard Hessian-based methods. We evaluate the performance of the proposed methods based on the sensitivity and specificity like measures, and finally demonstrate further applicability of our method to vascular ultrasound images of the carotid artery, as well as ultrasound data for abdominal aortic aneurysms.

Description: Identification of error in nonrigid registration is a critical problem in the medical image processing community. We recently proposed an algorithm that we call “Assessing Quality Using Image Registration Circuits” (AQUIRC) to identify nonrigid registration errors and have tested its performance using simulated cases. In this paper, we extend our previous work to assess AQUIRC's ability to detect local nonrigid registration errors and validate it quantitatively at specific clinical landmarks, namely the anterior commissure and the posterior commissure. To test our approach on a representative range of error we utilize five different registration methods and use 100 target images and nine atlas images. Our results show that AQUIRC's measure of registration quality correlates with the true target registration error (TRE) at these selected landmarks with an ${rm R}^{2}=0.542$ . To compare our method to a more conventional approach, we compute local normalized correlation coefficient (LNCC) and show that AQUIRC performs similarly. However, a multi-linear regression performed with both AQUIRC's measure and LNCC shows a higher correlation with TRE than correlations obtained with either measure alone, thus showing the complementarity of these quality measures. We conclude the paper by showing that the AQUIRC algorithm can be used to reduce registration errors for all five algorithms.

Description: Image reconstruction from low-count positron emission tomography (PET) projection data is challenging because the inverse problem is ill-posed. Prior information can be used to improve image quality. Inspired by the kernel methods in machine learning, this paper proposes a kernel based method that models PET image intensity in each pixel as a function of a set of features obtained from prior information. The kernel-based image model is incorporated into the forward model of PET projection data and the coefficients can be readily estimated by the maximum likelihood (ML) or penalized likelihood image reconstruction. A kernelized expectation-maximization algorithm is presented to obtain the ML estimate. Computer simulations show that the proposed approach can achieve better bias versus variance trade-off and higher contrast recovery for dynamic PET image reconstruction than the conventional maximum likelihood method with and without post-reconstruction denoising. Compared with other regularization-based methods, the kernel method is easier to implement and provides better image quality for low-count data. Application of the proposed kernel method to a 4-D dynamic PET patient dataset showed promising results.