ALBERT

Description: Accurately recovering the hippocampal shapes against rough and noisy segmentations is as challenging as achieving good anatomical correspondence between the individual shapes. To address these issues, we propose a mesh-to-volume registration approach, characterized by a progressive model deformation. Our model implements flexible weighting scheme for model rigidity under a multi-level neighborhood for vertex connectivity. This method induces a large-to-small scale deformation of a template surface to build the pairwise correspondence by minimizing geometric distortion while robustly restoring the individuals' shape characteristics. We evaluated the proposed method's 1) accuracy and robustness in smooth surface reconstruction, 2) sensitivity in detecting significant shape differences between healthy control and disease groups (mild cognitive impairment and Alzheimer's disease), 3) robustness in constructing the anatomical correspondence between individual shape models, and 4) applicability in identifying subtle shape changes in relation to cognitive abilities in a healthy population. We compared the performance of the proposed method with other well-known methods—SPHARM-PDM, ShapeWorks and LDDMM volume registration with template injection—using various metrics of shape similarity, surface roughness, volume, and shape deformity. The experimental results showed that the proposed method generated smooth surfaces with less volume differences and better shape similarity to input volumes than others. The statistical analyses with clinical variables also showed that it was sensitive in detecting subtle shape changes of hippocampus.

Description: We propose a conditional random field (CRF) based classifier for segmentation of small enhanced pathologies. Specifically, we develop a temporal hierarchical adaptive texture CRF (THAT-CRF) and apply it to the challenging problem of gad enhancing lesion segmentation in brain MRI of patients with multiple sclerosis. In this context, the presence of many nonlesion enhancements (such as blood vessels) renders the problem more difficult. In addition to voxel-wise features, the framework exploits multiple higher order textures to discriminate the true lesional enhancements from the pool of other enhancements. Since lesional enhancements show more variation over time as compared to the nonlesional ones, we incorporate temporal texture analysis in order to study the textures of enhanced candidates over time. The parameters of the THAT-CRF model are learned based on 2380 scans from a multi-center clinical trial. The effect of different components of the model is extensively evaluated on 120 scans from a separate multi-center clinical trial. The incorporation of the temporal textures results in a general decrease of the false discovery rate. Specifically, THAT-CRF achieves overall sensitivity of 95% along with false discovery rate of 20% and average false positive count of 0.5 lesions per scan. The sensitivity of the temporal method to the trained time interval is further investigated on five different intervals of 69 patients. Moreover, superior performance is achieved by the reviewed labelings of our model compared to the fully manual labeling when applied to the context of separating different treatment arms in a real clinical trial.

Description: We present a novel general-purpose compression method for tomographic images, termed 3D adaptive sparse representation based compression (3D-ASRC). In this paper, we focus on applications of 3D-ASRC for the compression of ophthalmic 3D optical coherence tomography (OCT) images. The 3D-ASRC algorithm exploits correlations among adjacent OCT images to improve compression performance, yet is sensitive to preserving their differences. Due to the inherent denoising mechanism of the sparsity based 3D-ASRC, the quality of the compressed images are often better than the raw images they are based on. Experiments on clinical-grade retinal OCT images demonstrate the superiority of the proposed 3D-ASRC over other well-known compression methods.

Description: Accurate segmentation is usually crucial in transrectal ultrasound (TRUS) image based prostate diagnosis; however, it is always hampered by heavy speckles. Contrary to the traditional view that speckles are adverse to segmentation, we exploit intrinsic properties induced by speckles to facilitate the task, based on the observations that sizes and orientations of speckles provide salient cues to determine the prostate boundary. Since the speckle orientation changes in accordance with a statistical prior rule, rotation-invariant texture feature is extracted along the orientations revealed by the rule. To address the problem of feature changes due to different speckle sizes, TRUS images are split into several arc-like strips. In each strip, every individual feature vector is sparsely represented, and representation residuals are obtained. The residuals, along with the spatial coherence inherited from biological tissues, are combined to segment the prostate preliminarily via graph cuts. After that, the segmentation is fine-tuned by a novel level sets model, which integrates 1) the prostate shape prior, 2) dark-to-light intensity transition near the prostate boundary, and 3) the texture feature just obtained. The proposed method is validated on two 2-D image datasets obtained from two different sonographic imaging systems, with the mean absolute distance on the mid gland images only $1.06pm 0.53~{hbox {mm}}$ and $1.25pm 0.77~{hbox {mm}}$ , respectively. The method is also extended to segment apex and base images, producing competitive results over the state of the art.

Description: Alzheimer's disease (AD) is a progressive and irreversible neurodegenerative disorder that has recently seen serious increase in the number of affected subjects. In the last decade, neuroimaging has been shown to be a useful tool to understand AD and its prodromal stage, amnestic mild cognitive impairment (MCI). The majority of AD/MCI studies have focused on disease diagnosis, by formulating the problem as classification with a binary outcome of AD/MCI or healthy controls. There have recently emerged studies that associate image scans with continuous clinical scores that are expected to contain richer information than a binary outcome. However, very few studies aim at modeling multiple clinical scores simultaneously, even though it is commonly conceived that multivariate outcomes provide correlated and complementary information about the disease pathology. In this article, we propose a sparse multi-response tensor regression method to model multiple outcomes jointly as well as to model multiple voxels of an image jointly. The proposed method is particularly useful to both infer clinical scores and thus disease diagnosis, and to identify brain subregions that are highly relevant to the disease outcomes. We conducted experiments on the Alzheimer's Disease Neuroimaging Initiative (ADNI) dataset, and showed that the proposed method enhances the performance and clearly outperforms the competing solutions.

Description: A plethora of techniques for cardiac deformation imaging with 3D ultrasound, typically referred to as 3D speckle tracking techniques, are available from academia and industry. Although the benefits of single methods over alternative ones have been reported in separate publications, the intrinsic differences in the data and definitions used makes it hard to compare the relative performance of different solutions. To address this issue, we have recently proposed a framework to simulate realistic 3D echocardiographic recordings and used it to generate a common set of ground-truth data for 3D speckle tracking algorithms, which was made available online. The aim of this study was therefore to use the newly developed database to contrast non-commercial speckle tracking solutions from research groups with leading expertise in the field. The five techniques involved cover the most representative families of existing approaches, namely block-matching, radio-frequency tracking, optical flow and elastic image registration. The techniques were contrasted in terms of tracking and strain accuracy. The feasibility of the obtained strain measurements to diagnose pathology was also tested for ischemia and dyssynchrony.

Description: Ultra-high field MRI is an area of great interest for clinical research and basic science due to the increased signal-to-noise, spatial resolution and magnetic-susceptibility-based contrast. However, the fact that the electromagnetic wavelength in tissue is comparable to the relevant body dimensions means that the uniformity of the excitation field is much poorer than at lower field strengths. In addition to techniques such as transmit arrays, one simple but effective method to counteract this effect is to use high permittivity “pads”. Very high permittivities enable thinner, flexible pads to be used, but the limiting factor is wavelength effects within the pads themselves, which can lead to image artifacts. So far, all studies have used simple continuous rectangular/circular pad geometries. In this work we investigate how the wavelength effects can be partially mitigated utilizing shaped pad with holes. Several arrangements have been simulated, including low order pre-fractal geometries, which maintain the overall coverage of the pad, but can provide better image homogeneity in the region of interest or higher sensitivity depending on the setup. Experimental data in the form of in vivo human images at 7T were acquired to validate the simulation results.

Description: In this paper, we consider combined TV denoising and diffusion tensor fitting in DTI using the affine-invariant Riemannian metric on the space of diffusion tensors. Instead of first fitting the diffusion tensors, and then denoising them, we define a suitable TV type energy functional which incorporates the measured DWIs (using an inverse problem setup) and which measures the nearness of neighboring tensors in the manifold. To approach this functional, we propose generalized forward- backward splitting algorithms which combine an explicit and several implicit steps performed on a decomposition of the functional. We validate the performance of the derived algorithms on synthetic and real DTI data. In particular, we work on real 3D data. To our knowledge, the present paper describes the first approach to TV regularization in a combined manifold and inverse problem setup.

Description: We present a generic method for automatic detection of abnormal regions in medical images as deviations from a normative data base. The algorithm decomposes an image, or more broadly a function defined on the image grid, into the superposition of a normal part and a residual term. A statistical model is constructed with regional sparse learning to represent normative anatomical variations among a reference population (e.g., healthy controls), in conjunction with a Markov random field regularization that ensures mutual consistency of the regional learning among partially overlapping image blocks. The decomposition is performed in a principled way so that the normal part fits well with the learned normative model, while the residual term absorbs pathological patterns, which may then be detected through a statistical significance test. The decomposition is applied to multiple image features from an individual scan, detecting abnormalities using both intensity and shape information. We form an iterative scheme that interleaves abnormality detection with deformable registration, gradually improving robustness of the spatial normalization and precision of the detection. The algorithm is evaluated with simulated images and clinical data of brain lesions, and is shown to achieve robust deformable registration and localize pathological regions simultaneously. The algorithm is also applied on images from Alzheimer’s disease patients to demonstrate the generality of the method.

Description: Staining and scanning of tissue samples for microscopic examination is fraught with undesirable color variations arising from differences in raw materials and manufacturing techniques of stain vendors, staining protocols of labs, and color responses of digital scanners. When comparing tissue samples, color normalization and stain separation of the tissue images can be helpful for both pathologists and software. Techniques that are used for natural images fail to utilize structural properties of stained tissue samples and produce undesirable color distortions. The stain concentration cannot be negative. Tissue samples are stained with only a few stains and most tissue regions are characterized by at most one effective stain. We model these physical phenomena that define the tissue structure by first decomposing images in an unsupervised manner into stain density maps that are sparse and non-negative. For a given image, we combine its stain density maps with stain color basis of a pathologist-preferred target image, thus altering only its color while preserving its structure described by the maps. Stain density correlation with ground truth and preference by pathologists were higher for images normalized using our method when compared to other alternatives. We also propose a computationally faster extension of this technique for large whole-slide images that selects an appropriate patch sample instead of using the entire image to compute the stain color basis.

Description: Prospective authors are requested to submit new, unpublished manuscripts for inclusion in the upcoming event described in this call for papers.

Description: We present a methodology to recover the geometrical calibration of conventional X-ray settings with the help of an ordinary video camera and visible fiducials that are present in the scene. After calibration, equivalent points of interest can be easily identifiable with the help of the epipolar geometry. The same procedure also allows the measurement of real anatomic lengths and angles and obtains accurate 3D locations from image points. Our approach completely eliminates the need for X-ray-opaque reference marks (and necessary supporting frames) which can sometimes be invasive for the patient, occlude the radiographic picture, and end up projected outside the imaging sensor area in oblique protocols. Two possible frameworks are envisioned: a spatially shifting X-ray anode around the patient/object and a moving patient that moves/rotates while the imaging system remains fixed. As a proof of concept, experiences with a device under test (DUT), an anthropomorphic phantom and a real brachytherapy session have been carried out. The results show that it is possible to identify common points with a proper level of accuracy and retrieve three-dimensional locations, lengths and shapes with a millimetric level of precision. The presented approach is simple and compatible with both current and legacy widespread diagnostic X-ray imaging deployments and it can represent a good and inexpensive alternative to other radiological modalities like CT.

Description: Fiber based fluorescence lifetime imaging has shown great potential for intraoperative diagnosis and guidance of surgical procedures. Here we describe a novel method addressing a significant challenge for the practical implementation of this technique, i.e., the real-time display of the quantified biochemical or functional tissue properties superimposed on the interrogated area. Specifically, an aiming beam (450 nm) generated by a continuous-wave laser beam was merged with the pulsed fluorescence excitation light in a single delivery/collection fiber and then imaged and segmented using a color-based algorithm. We demonstrate that this approach enables continuous delineation of the interrogated location and dynamic augmentation of the acquired frames with the corresponding fluorescence decay parameters. The method was evaluated on a fluorescence phantom and fresh tissue samples. Current results demonstrate that 34 frames per second can be achieved for augmenting videos of 640 $times $ 512 pixels resolution. Also we show that the spatial resolution of the fluorescence lifetime map depends on the tissue optical properties, the scanning speed, and the frame rate. The dice similarity coefficient between the fluorescence phantom and the reconstructed maps was estimated to be as high as 93%. The reported method could become a valuable tool for augmenting the surgeon's field of view with diagnostic information derived from the analysis of fluorescence lifetime data in real-time using handheld, automated, or endoscopic scanning systems. Current method provides also a means for maintaining the tissue light exposure within safety limits. This study provides a framework for using an aiming beam with other point spectroscopy applications.

Description: Simultaneous Multi-Slice (SMS) magnetic resonance imaging (MRI) is a rapidly evolving technique for increasing imaging speed. Controlled aliasing techniques utilize periodic undersampling patterns to help mitigate the loss in signal-to-noise ratio (SNR) in SMS MRI. To evaluate the performance of different undersampling patterns, a quantitative description of the image SNR loss is needed. Additionally, eddy current effects in echo planar imaging (EPI) lead to slice-specific Nyquist ghosting artifacts. These artifacts cannot be accurately corrected for each individual slice before or after slice-unaliasing. In this work, we propose a hybrid-space sensitivity encoding (SENSE) reconstruction framework for SMS MRI by adopting a three-dimensional representation of the SMS acquisition. Analytical SNR loss maps are derived for SMS acquisitions with arbitrary phase encoding undersampling patterns. Moreover, we propose a matrix-decoding correction method that corrects the slice-specific Nyquist ghosting artifacts in SMS EPI acquisitions. Brain images demonstrate that the proposed hybrid-space SENSE reconstruction generates images with comparable quality to commonly used split-slice-generalized autocalibrating partially parallel acquisition reconstruction. The analytical SNR loss maps agree with those calculated by a Monte Carlo based method, but require less computation time for high quality maps. The analytical maps enable a fair comparison between the performances of coherent and incoherent SMS undersampling patterns. Phantom and brain SMS EPI images show that the matrix-decoding method performs better than the single-slice and slice-averaged Nyquist ghosting correction methods under the hybrid-space SENSE reconstruction framework.

Description: Photoacoustic (PA) images utilize pulsed lasers and ultrasound transducers to visualize targets with higher optical absorption than the surrounding medium. However, they are susceptible to acoustic clutter and background noise artifacts that obfuscate biomedical structures of interest. We investigated three spatial-angular compounding methods to improve PA image quality for biomedical applications, implemented by combining multiple images acquired as an ultrasound probe was rotated about the elevational axis with the laser beam and target fixed. Compounding with conventional averaging was based on the pose information of each PA image, while compounding with weighted and selective averaging utilized both the pose and image content information. Weighted-average compounding enhanced PA images with the least distortion of signal size, particularly when there were large (i.e., 2.5 mm and 7 $^{circ}$ ) perturbations from the initial probe position. Selective-average compounding offered the best improvement in image quality with up 181, 1665, and 1568 times higher contrast, CNR, and SNR, respectively, compared to the mean values of individual PA images. The three presented spatial compounding methods have promising potential to enhance image quality in multiple photoacoustic applications.

Description: Analysis of cranial nerve systems, such as the anterior visual pathway (AVP), from MRI sequences is challenging due to their thin long architecture, structural variations along the path, and low contrast with adjacent anatomic structures. Segmentation of a pathologic AVP (e.g., with low-grade gliomas) poses additional challenges. In this work, we propose a fully automated partitioned shape model segmentation mechanism for AVP steered by multiple MRI sequences and deep learning features. Employing deep learning feature representation, this framework presents a joint partitioned statistical shape model able to deal with healthy and pathological AVP. The deep learning assistance is particularly useful in the poor contrast regions, such as optic tracts and pathological areas. Our main contributions are: 1) a fast and robust shape localization method using conditional space deep learning, 2) a volumetric multiscale curvelet transform-based intensity normalization method for robust statistical model, and 3) optimally partitioned statistical shape and appearance models based on regional shape variations for greater local flexibility. Our method was evaluated on MRI sequences obtained from 165 pediatric subjects. A mean Dice similarity coefficient of 0.779 was obtained for the segmentation of the entire AVP (optic nerve only $=0.791$ ) using the leave-one-out validation. Results demonstrated that the proposed localized shape and sparse appearance-based learning approach significantly outperforms current state-of-the-art segmentation approaches and is as robust as the manual segmentation.

Description: Analytical (closed-form) inversion schemes have been the standard approach for image reconstruction in optoacoustic tomography due to their fast reconstruction abilities and low memory requirements. Yet, the need for quantitative imaging and artifact reduction has led to the development of more accurate inversion approaches, which rely on accurate forward modeling of the optoacoustic wave generation and propagation. In this way, multiple experimental factors can be incorporated, such as the exact detection geometry, spatio-temporal response of the transducers, and acoustic heterogeneities. The model-based inversion commonly results in very large sparse matrix formulations that require computationally extensive and memory demanding regularization schemes for image reconstruction, hindering their effective implementation in real-time imaging applications. Herein, we introduce a new discretization procedure for efficient model-based reconstructions in two-dimensional optoacoustic tomography that allows for parallel implementation on a graphics processing unit (GPU) with a relatively low numerical complexity. By on-the-fly calculation of the model matrix in each iteration of the inversion procedure, the new approach results in imaging frame rates exceeding 10 Hz, thus enabling real-time image rendering using the model-based approach.

Description: This paper introduces a statistical estimation framework for magnetic resonance (MR) fingerprinting, a recently proposed quantitative imaging paradigm. Within this framework, we present a maximum likelihood (ML) formalism to estimate multiple MR tissue parameter maps directly from highly undersampled, noisy ${bf k}$ -space data. A novel algorithm, based on variable splitting, the alternating direction method of multipliers, and the variable projection method, is developed to solve the resulting optimization problem. Representative results from both simulations and in vivo experiments demonstrate that the proposed approach yields significantly improved accuracy in parameter estimation, compared to the conventional MR fingerprinting reconstruction. Moreover, the proposed framework provides new theoretical insights into the conventional approach. We show analytically that the conventional approach is an approximation to the ML reconstruction; more precisely, it is exactly equivalent to the first iteration of the proposed algorithm for the ML reconstruction, provided that a gridding reconstruction is used as an initialization.

Description: We propose a unified Bayesian framework to detect both hyper- and hypo-active communities within whole-brain fMRI data. Specifically, our model identifies dense subgraphs that exhibit population-level differences in functional synchrony between a control and clinical group. We derive a variational EM algorithm to solve for the latent posterior distributions and parameter estimates, which subsequently inform us about the afflicted network topology. We demonstrate that our method provides valuable insights into the neural mechanisms underlying social dysfunction in autism, as verified by the Neurosynth meta-analytic database. In contrast, both univariate testing and community detection via recursive edge elimination fail to identify stable functional communities associated with the disorder.

Description: A recent entry into the rapidly evolving field of integrated PET/MR scanners is presented in this paper: a whole body hybrid PET/MR system (SIGNA PET/MR, GE Healthcare) capable of simultaneous acquisition of both time-of-flight (TOF) PET and high resolution MR data. The PET ring was integrated into an existing 3T MR system resulting in a (patient) bore opening of 60 cm diameter, with a 25 cm axial FOV. PET performance was evaluated both on the standalone PET ring and on the same detector integrated into the MR system, to assess the level of mutual interference between both subsystems. In both configurations we obtained detector performance data. PET detector performance was not significantly affected by integration into the MR system. The global energy resolution was within 2% (10.3% versus 10.5%), and the system coincidence time resolution showed a maximum change of 〈 3% (385 ps versus 394 ps) when measured outside MR and during simultaneous PET/MRI acquisitions, respectively. To evaluate PET image quality and resolution, the NEMA IQ phantom was acquired with MR idle and with MR active. Impact of PET on MR IQ was assessed by comparing SNR with PET acquisition on and off. B0 and B1 homogeneities were acquired before and after the integration of the PET ring inside the magnet. In vivo brain and whole body head-to-thighs data were acquired to demonstrate clinical image quality.

Description: The performance of photoacoustic microscopy (PAM) degrades due to the turbidity of the skull that introduces attenuation and distortion of both laser and stimulated ultrasound. In this manuscript, we demonstrated that a newly developed skull optical clearing solution (SOCS) could enhance not only the transmittance of light, but also that of ultrasound in the skull in vitro. Thus the photoacoustic signal was effectively elevated, and the relative strength of the artifacts induced by the skull could be suppressed. Furthermore in vivo studies demonstrated that SOCS could drastically enhance the performance of photoacoustic microscopy for cerebral microvasculature imaging.

Description: Cryo-balloon catheters have attracted an increasing amount of interest in the medical community as they can reduce patient risk during left atrial pulmonary vein ablation procedures. As cryo-balloon catheters are not equipped with electrodes, they cannot be localized automatically by electro-anatomical mapping systems. As a consequence, X-ray fluoroscopy has remained an important means for guidance during the procedure. Most recently, image guidance methods for fluoroscopy-based procedures have been proposed, but they provide only limited support for cryo-balloon catheters and require significant user interaction. To improve this situation, we propose a novel method for automatic cryo-balloon catheter detection in fluoroscopic images by detecting the cryo-balloon catheter's built-in X-ray marker. Our approach is based on a blob detection algorithm to find possible X-ray marker candidates. Several of these candidates are then excluded using prior knowledge. For the remaining candidates, several catheter specific features are introduced. They are processed using a machine learning approach to arrive at the final X-ray marker position. Our method was evaluated on 75 biplane fluoroscopy images from 40 patients, from two sites, acquired with a biplane angiography system. The method yielded a success rate of 99.0% in plane A and 90.6% in plane B, respectively. The detection achieved an accuracy of $1.00~{rm mm}pm 0.82~{rm mm}$ in plane A and $1.13~{rm mm}pm 0.24~{rm mm}$ in plane B. The localization in 3-D was associated with an average error of $0.36~{rm mm}pm 0.86~{rm mm}$ .

Description: Segmentation of the wrist bones in CT images has been frequently used in different clinical applications including arthritis evaluation, bone age assessment and image-guided interventions. The major challenges include non-uniformity and spongy textures of the bone tissue as well as narrow inter-bone spaces. In this work, we propose an automatic wrist bone segmentation technique for CT images based on a statistical model that captures the shape and pose variations of the wrist joint across 60 example wrists at nine different wrist positions. To establish the correspondences across the training shapes at neutral positions, the wrist bone surfaces are jointly aligned using a group-wise registration framework based on a Gaussian Mixture Model. Principal component analysis is then used to determine the major modes of shape variations. The variations in poses not only across the population but also across different wrist positions are incorporated in two pose models. An intra-subject pose model is developed by utilizing the similarity transforms at all wrist positions across the population. Further, an inter-subject pose model is used to model the pose variations across different wrist positions. For segmentation of the wrist bones in CT images, the developed model is registered to the edge point cloud extracted from the CT volume through an expectation maximization based probabilistic approach. Residual registration errors are corrected by application of a non-rigid registration technique. We validate the proposed segmentation method by registering the wrist model to a total of 66 unseen CT volumes of average voxel size of 0.38 mm. We report a mean surface distance error of 0.33 mm and a mean Jaccard index of 0.86.

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

Description: The paper presents a new endoscope motion tracking method that is based on a novel external endoscope tracking device and our modified stochastic optimization method for boosting endoscopy navigation. We designed a novel tracking prototype where a 2-D motion sensor was introduced to directly measure the insertion–retreat linear motion and also the rotation of the endoscope. With our developed stochastic optimization method, which embeds traceable particle swarm optimization in the Condensation algorithm, a full six degrees-of-freedom endoscope pose (position and orientation) can be recovered from 2-D motion sensor measurements. Experiments were performed on a dynamic bronchial phantom with maximal simulated respiratory motion around 24.0 mm. The experimental results demonstrate that our proposed method provides a promising endoscope motion tracking approach with more effective and robust performance than several current available tracking techniques. The average tracking accuracy of the position improved from 6.5 to 3.3 mm, which further approaches the clinical requirement of 2.0 mm in practice.

Description: Radical prostatectomy is performed on approximately 40% of men with organ-confined prostate cancer. Pathologic information obtained from the prostatectomy specimen provides important prognostic information and guides recommendations for adjuvant treatment. The current pathology protocol in most centers involves primarily qualitative assessment. In this paper, we describe and evaluate our system for automatic prostate cancer detection and grading on hematoxylin & eosin-stained tissue images. Our approach is intended to address the dual challenges of large data size and the need for high-level tissue information about the locations and grades of tumors. Our system uses two stages of AdaBoost-based classification. The first provides high-level tissue component labeling of a superpixel image partitioning. The second uses the tissue component labeling to provide a classification of cancer versus noncancer, and low-grade versus high-grade cancer. We evaluated our system using 991 sub-images extracted from digital pathology images of 50 whole-mount tissue sections from 15 prostatectomy patients. We measured accuracies of 90% and 85% for the cancer versus noncancer and high-grade versus low-grade classification tasks, respectively. This system represents a first step toward automated cancer quantification on prostate digital histopathology imaging, which could pave the way for more accurately informed postprostatectomy patient care.

Description: We present DevStaR, an automated computer vision and machine learning system that provides rapid, accurate, and quantitative measurements of C. elegans embryonic viability in high-throughput (HTP) applications. A leading genetic model organism for the study of animal development and behavior, C. elegans is particularly amenable to HTP functional genomic analysis due to its small size and ease of cultivation, but the lack of efficient and quantitative methods to score phenotypes has become a major bottleneck. DevStaR addresses this challenge using a novel hierarchical object recognition machine that rapidly segments, classifies, and counts animals at each developmental stage in images of mixed-stage populations of C. elegans. Here, we describe the algorithmic design of the DevStaR system and demonstrate its performance in scoring image data acquired in HTP screens.

Description: Atrial fibrillation is an increasingly prevalent cardiovascular disease; changes in atrial structure and function induced by atrial fibrillation and its treatments are often spatially heterogeneous. However, spatial heterogeneity of function is difficult to assess with standard imaging techniques. This paper describes a method to assess global and regional mechanical function by combining cardiac magnetic resonance imaging and finite-element surface fitting. We used this fitted surface to derive measures of left atrial volume, regional motion, and spatial heterogeneity of motion in 23 subjects, including healthy volunteers and atrial fibrillation patients. We fit the surfaces using a Newton optimization scheme in under 1 min on a standard laptop, with a root mean square error of $2.3pm 0.5$ mm, less than 9% of the mean fitted radius, and an inter-operator variability of less than 10%. Fitted surfaces showed clear definition of the phases of left atrial motion (filling, passive emptying, active contraction) in both volume-time and regional radius-time curves. Averaged surfaces of healthy volunteers and atrial fibrillation patients provided evidence of substantial regional variation in both amount and timing of regional motion, indicating spatial heterogeneity of function, even in healthy adults.

Description: Magnetic resonance spectroscopic imaging (MRSI) is a powerful tool capable of providing spatially localized maps of metabolite concentrations. Its utility, however, is often depreciated by spectral leakage artifacts resulting from low spatial resolution measurements through an effort to reduce acquisition times. Though model-based techniques can help circumvent these drawbacks, they require strong prior knowledge, and can introduce additional artifacts when the underlying models are inaccurate. We introduce a novel scheme in which a generative model is estimated from the raw MRSI data via a regularized variational framework that minimizes the model approximation error within a measurement-prescribed subspace. As additional a priori information, our approach relies only upon a measured field inhomogeneity map at high spatial resolution. We demonstrate the feasibility of our approach on both synthetic and experimental data.

Description: The absorbed energy distribution (AED) in X-ray imaging detectors is an important factor that affects both energy resolution and image quality through the Swank factor and detective quantum efficiency. In the diagnostic energy range (20–140 keV), escape of characteristic photons following photoelectric absorption and Compton scatter photons are primary sources of absorbed-energy dispersion in X-ray detectors. In this paper, we describe the development of an analytic model of the AED in compound X-ray detector materials, based on the cascaded-systems approach, that includes the effects of escape and reabsorption of characteristic and Compton-scatter photons. We derive analytic expressions for both semi-infinite slab and pixel geometries and validate our approach by Monte Carlo simulations. The analytic model provides the energy-dependent X-ray response function of arbitrary compound materials without time-consuming Monte Carlo simulations. We believe this model will be useful for correcting spectral distortion artifacts commonly observed in photon-counting applications and optimal design and development of novel X-ray detectors.

Description: We investigate an approach to evaluation of emission-tomography (ET) imaging systems used for region-of-interest (ROI) estimation tasks. In the evaluation we employ the concept of “emission counts” (EC), which are the number of events per voxel emitted during a scan. We use the reduction in posterior variance of ROI EC, compared to the prior ROI EC variance, as the metric of primary interest, which we call the “posterior variance reduction index” (PVRI). Systems that achieve a higher PVRI are considered superior to systems with lower PVRI. The approach is independent of the reconstruction method and is applicable to all photon-limited data types including list-mode data. We analyzed this approach using a model of 2-D tomography, and compared our results to the classical theory of tomographic sampling. We found that performance evaluations using the PVRI index were consistent with the classical theory. System evaluation based on EC posterior variance is an intuitively appealing and physically meaningful method that is useful for evaluation of system performance in ROI quantitation tasks.

Description: This paper presents a novel probabilistic voxel selection strategy for medical image registration in time-sensitive contexts, where the goal is aggressive voxel sampling (e.g., using less than 1% of the total number) while maintaining registration accuracy and low failure rate. We develop a Bayesian framework whereby, first, a voxel sampling probability field (VSPF) is built based on the uncertainty on the transformation parameters. We then describe a practical, multi-scale registration algorithm, where, at each optimization iteration, different voxel subsets are sampled based on the VSPF. The approach maximizes accuracy without committing to a particular fixed subset of voxels. The probabilistic sampling scheme developed is shown to manage the tradeoff between the robustness of traditional random voxel selection (by permitting more exploration) and the accuracy of fixed voxel selection (by permitting a greater proportion of informative voxels).

Description: Segmentation of the spinal column from computed tomography (CT) images is a preprocessing step for a range of image-guided interventions. One intervention that would benefit from accurate segmentation is spinal needle injection. Previous spinal segmentation techniques have primarily focused on identification and separate segmentation of each vertebra. Recently, statistical multi-object shape models have been introduced to extract common statistical characteristics between several anatomies. These models can be used for segmentation purposes because they are robust, accurate, and computationally tractable. In this paper, we develop a statistical multi-vertebrae shape+pose model and propose a novel registration-based technique to segment the CT images of spine. The multi-vertebrae statistical model captures the variations in shape and pose simultaneously, which reduces the number of registration parameters. We validate our technique in terms of accuracy and robustness of multi-vertebrae segmentation of CT images acquired from lumbar vertebrae of 32 subjects. The mean error of the proposed technique is below 2 mm, which is sufficient for many spinal needle injection procedures, such as facet joint injections.

Description: In this paper, we report on the first edition of the HEp-2 Cells Classification contest, held at the 2012 edition of the International Conference on Pattern Recognition, and focused on indirect immunofluorescence (IIF) image analysis. The IIF methodology is used to detect autoimmune diseases by searching for antibodies in the patient serum but, unfortunately, it is still a subjective method that depends too heavily on the experience and expertise of the physician. This has been the motivation behind the recent initial developments of computer aided diagnosis systems in this field. The contest aimed to bring together researchers interested in the performance evaluation of algorithms for IIF image analysis: 28 different recognition systems able to automatically recognize the staining pattern of cells within IIF images were tested on the same undisclosed dataset. In particular, the dataset takes into account the six staining patterns that occur most frequently in the daily diagnostic practice: centromere, nucleolar, homogeneous, fine speckled, coarse speckled, and cytoplasmic. In the paper, we briefly describe all the submitted methods, analyze the obtained results, and discuss the design choices conditioning the performance of each method.

Description: Tissue displacements required for mechanical property reconstruction in magnetic resonance elastography (MRE) are acquired in a magnetic resonance imaging (MRI) scanner, therefore, anatomical information is available from other imaging sequences. Despite its availability, few attempts to incorporate prior spatial information in the MRE reconstruction process have been reported. This paper implements and evaluates soft prior regularization (SPR), through which homogeneity in predefined spatial regions is enforced by a penalty term in a nonlinear inversion strategy. Phantom experiments and simulations show that when predefined regions are spatially accurate, recovered property values are stable for SPR weighting factors spanning several orders of magnitude, whereas inaccurate segmentation results in bias in the reconstructed properties that can be mitigated through proper choice of regularization weighting. The method was evaluated in vivo by estimating viscoelastic mechanical properties of frontal lobe gray and white matter for five repeated scans of a healthy volunteer. Segmentations of each tissue type were generated using automated software, and statistically significant differences between frontal lobe gray and white matter were found for both the storage modulus and loss modulus $(p 〈 10^{-6})$ . Provided homogeneous property assumptions are reasonable, SPR produces accurate quantitative property estimates for tissue structures which are finer than the resolution currently achievable with fully distributed MRE.

Description: Recent research has demonstrated that improved image segmentation can be achieved by multiple template fusion utilizing both label and intensity information. However, intensity weighted fusion approaches use local intensity similarity as a surrogate measure of local template quality for predicting target segmentation and do not seek to characterize template performance. This limits both the usefulness and accuracy of these techniques. Our work here was motivated by the observation that the local intensity similarity is a poor surrogate measure for direct comparison of the template image with the true image target segmentation. Although the true image target segmentation is not available, a high quality estimate can be inferred, and this in turn allows a principled estimate to be made of the local quality of each template at contributing to the target segmentation. We developed a fusion algorithm that uses probabilistic segmentations of the target image to simultaneously infer a reference standard segmentation of the target image and the local quality of each probabilistic segmentation. The concept of comparing templates to a hidden reference standard segmentation enables accurate assessments of the contribution of each template to inferring the target image segmentation to be made, and in practice leads to excellent target image segmentation. We have used the new algorithm for the multiple-template-based segmentation and parcellation of magnetic resonance images of the brain. Intensity and label map images of each one of the aligned templates are used to train a local Gaussian mixture model based classifier. Then, each classifier is used to compute the probabilistic segmentations of the target image. Finally, the generated probabilistic segmentations are fused together using the new fusion algorithm to obtain the segmentation of the target image. We evaluated our method in comparison to other state-of-the-art segmentation methods. We demonstrated that our new fusio- algorithm has higher segmentation performance than these methods.

Description: The electroencephalography source estimation problem consists of inferring cortical activation from measurements of electrical potential taken on the scalp surface. This inverse problem is intrinsically ill-posed. In particular the dimensionality of cortical sources greatly exceeds the number of electrode measurements, and source estimation requires regularization to obtain a unique solution. In this work, we introduce a novel regularization function called cortical graph smoothing, which exploits knowledge of anatomical connectivity available from diffusion-weighted imaging. Given a weighted graph description of the anatomical connectivity of the brain, cortical graph smoothing penalizes the weighted sum of squares of differences of cortical activity across the graph edges, thus encouraging solutions with consistent activation across anatomically connected regions. We explore the performance of the cortical graph smoothing source estimates for analysis of the event related potential for simple motor tasks, and compare against the commonly used minimum norm, weighted minimum norm, LORETA and sLORETA source estimation methods. Evaluated over a series of 18 subjects, the proposed cortical graph smoothing method shows superior localization accuracy compared to the minimum norm method, and greater relative peak intensity than the other comparison methods.

Description: A novel framework for the classification of acute rejection versus nonrejection status of renal transplants from 2-D dynamic contrast-enhanced magnetic resonance imaging is proposed. The framework consists of four steps. First, kidney objects are segmented from adjacent structures with a level set deformable boundary guided by a stochastic speed function that accounts for a fourth-order Markov–Gibbs random field model of the kidney/background shape and appearance. Second, a Laplace-based nonrigid registration approach is used to account for local deformations caused by physiological effects. Namely, the target kidney object is deformed over closed, equispaced contours (iso-contours) to closely match the reference object. Next, the cortex is segmented as it is the functional kidney unit that is most affected by rejection. To characterize rejection, perfusion is estimated from contrast agent kinetics using empirical indexes, namely, the transient phase indexes (peak signal intensity, time-to-peak, and initial up-slope), and a steady-phase index defined as the average signal change during the slowly varying tissue phase of agent transit. We used a $k_n$ -nearest neighbor classifier to distinguish between acute rejection and nonrejection. Performance of our method was evaluated using the receiver operating characteristics (ROC). Experimental results in 50 subjects, using a combinatoric $k_n$ -classifier, correctly classified 92% of training subjects, 100% of the test subjects, and yielded an area under the ROC curve that approached the ideal value. Our proposed framework thus holds promise as a reliable noninvasive diagnostic tool.

Description: Longitudinal imaging studies are frequently used to investigate temporal changes in brain morphology and often require spatial correspondence between images achieved through image registration. Beside morphological changes, image intensity may also change over time, for example when studying brain maturation. However, such intensity changes are not accounted for in image similarity measures for standard image registration methods. Hence, 1) local similarity measures, 2) methods estimating intensity transformations between images, and 3) metamorphosis approaches have been developed to either achieve robustness with respect to intensity changes or to simultaneously capture spatial and intensity changes. For these methods, longitudinal intensity changes are not explicitly modeled and images are treated as independent static samples. Here, we propose a model-based image similarity measure for longitudinal image registration that estimates a temporal model of intensity change using all available images simultaneously.

Description: We present a new approach for the automated segmentation of synapses in image stacks acquired by electron microscopy (EM) that relies on image features specifically designed to take spatial context into account. These features are used to train a classifier that can effectively learn cues such as the presence of a nearby post-synaptic region. As a result, our algorithm successfully distinguishes synapses from the numerous other organelles that appear within an EM volume, including those whose local textural properties are relatively similar. Furthermore, as a by-product of the segmentation, our method flawlessly determines synaptic orientation, a crucial element in the interpretation of brain circuits. We evaluate our approach on three different datasets, compare it against the state-of-the-art in synapse segmentation and demonstrate our ability to reliably collect shape, density, and orientation statistics over hundreds of synapses.

Description: Describes the above-named upcoming special issue or section. May include topics to be covered or calls for papers.

Description: The use of diffusion tensor imaging (DTI) for studying the human heart in vivo is very challenging due to cardiac motion. This paper assesses the effects of cardiac motion on the human myocardial fiber architecture. To this end, a model for analyzing the effects of cardiac motion on signal intensity is presented. A Monte-Carlo simulation based on polarized light imaging data is then performed to calculate the diffusion signals obtained by the displacement of water molecules, which generate diffusion weighted (DW) images. Rician noise and in vivo motion data obtained from DENSE acquisition are added to the simulated cardiac DW images to produce motion-induced datasets. An algorithm based on principal components analysis filtering and temporal maximum intensity projection (PCATMIP) is used to compensate for motion-induced signal loss. Diffusion tensor parameters derived from motion-reduced DW images are compared to those derived from the original simulated DW images. Finally, to assess cardiac motion effects on in vivo fiber architecture, in vivo cardiac DTI data processed by PCATMIP are compared to those obtained from one trigger delay (TD) or one single phase acquisition. The results showed that cardiac motion produced overestimated fractional anisotropy and mean diffusivity as well as a narrower range of fiber angles. The combined use of shifted TD acquisitions and postprocessing based on image registration and PCATMIP effectively improved the quality of in vivo DW images and subsequently, the measurement accuracy of fiber architecture properties. This suggests new solutions to the problems associated with obtaining in vivo human myocardial fiber architecture properties in clinical conditions.

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

Description: Electrical impedance tomography (EIT) is a medical imaging technique in which current is applied on electrodes on the surface of the body, the resulting voltage is measured, and an inverse problem is solved to recover the conductivity and/or permittivity in the interior. Images are then formed from the reconstructed conductivity and permittivity distributions. In the 2-D geometry, EIT is clinically useful for chest imaging. In this work, an implementation of a D-bar method for complex admittivities on a general 2-D domain is presented. In particular, reconstructions are computed on a chest-shaped domain for several realistic phantoms including a simulated pneumothorax, hyperinflation, and pleural effusion. The method demonstrates robustness in the presence of noise. Reconstructions from trigonometric and pairwise current injection patterns are included.

Description: This paper describes a task-based, information-theoretic approach to the assessment of image quality in diagnostic sonography. We expand the Kullback–Leibler divergence metric $J$ , which quantifies the diagnostic information contained within recorded radio-frequency echo signals, into a spatial-frequency integral comprised of two spectral components: one describes patient features for low-contrast diagnostic tasks and the other describes instrumentation properties. The latter quantity is the acquisition information spectrum (AIS), which measures the density of object information that an imaging system is able to transfer to the echo data at each spatial frequency. AIS is derived based on unique properties of acoustic scattering in tissues that generate object contrast. Predictions made by the $J$ integral expression were validated through Monte Carlo studies using echo-signal data from simulated lesions. Our analysis predicts the diagnostic performance of any sonographic system at specific diagnostic tasks based on engineering properties of the instrument that constitute image quality.

Description: Provides a listing of the editors, board members, and current staff for this issue of the publication.

Description: Electron and light microscopy imaging can now deliver high-quality image stacks of neural structures. However, the amount of human annotation effort required to analyze them remains a major bottleneck. While machine learning algorithms can be used to help automate this process, they require training data, which is time-consuming to obtain manually, especially in image stacks. Furthermore, due to changing experimental conditions, successive stacks often exhibit differences that are severe enough to make it difficult to use a classifier trained for a specific one on another. This means that this tedious annotation process has to be repeated for each new stack. In this paper, we present a domain adaptation algorithm that addresses this issue by effectively leveraging labeled examples across different acquisitions and significantly reducing the annotation requirements. Our approach can handle complex, nonlinear image feature transformations and scales to large microscopy datasets that often involve high-dimensional feature spaces and large 3D data volumes. We evaluate our approach on four challenging electron and light microscopy applications that exhibit very different image modalities and where annotation is very costly. Across all applications we achieve a significant improvement over the state-of-the-art machine learning methods and demonstrate our ability to greatly reduce human annotation effort.

Description: Efficient and accurate segmentation of cellular structures in microscopic data is an essential task in medical imaging. Many state-of-the-art approaches to image segmentation use structured models whose parameters must be carefully chosen for optimal performance. A popular choice is to learn them using a large-margin framework and more specifically structured support vector machines (SSVM). Although SSVMs are appealing, they suffer from certain limitations. First, they are restricted in practice to linear kernels because the more powerful nonlinear kernels cause the learning to become prohibitively expensive. Second, they require iteratively finding the most violated constraints, which is often intractable for the loopy graphical models used in image segmentation. This requires approximation that can lead to reduced quality of learning. In this paper, we propose three novel techniques to overcome these limitations. We first introduce a method to “kernelize” the features so that a linear SSVM framework can leverage the power of nonlinear kernels without incurring much additional computational cost. Moreover, we employ a working set of constraints to increase the reliability of approximate subgradient methods and introduce a new way to select a suitable step size at each iteration. We demonstrate the strength of our approach on both 2-D and 3-D electron microscopic (EM) image data and show consistent performance improvement over state-of-the-art approaches.

Description: In this study, we proposed an efficient nonrigid magnetic resonance (MR) to transrectal ultrasound (TRUS) deformable registration method in order to improve the accuracy of targeting suspicious regions during a three dimensional (3-D) TRUS guided prostate biopsy. The proposed deformable registration approach employs the multi-channel modality independent neighborhood descriptor (MIND) as the local similarity feature across the two modalities of MR and TRUS, and a novel and efficient duality-based convex optimization-based algorithmic scheme was introduced to extract the deformations and align the two MIND descriptors. The registration accuracy was evaluated using 20 patient images by calculating the TRE using manually identified corresponding intrinsic fiducials in the whole gland and peripheral zone. Additional performance metrics [Dice similarity coefficient (DSC), mean absolute surface distance (MAD), and maximum absolute surface distance (MAXD)] were also calculated by comparing the MR and TRUS manually segmented prostate surfaces in the registered images. Experimental results showed that the proposed method yielded an overall median TRE of 1.76 mm. The results obtained in terms of DSC showed an average of $80.8pm7.8hbox{%}$ for the apex of the prostate, $92.0pm3.4hbox{%}$ for the mid-gland, $81.7pm6.4hbox{%}$ for the base and $85.7pm4.7hbox{%}$ for the whole gland. The surface distance calculations showed an overall average of $1.84pm0.52$ mm for MAD and $6.90pm2.07$ mm for MAXD.

Description: In this work we present a framework for reliably detecting significant differences in quantitative magnetic resonance imaging and evaluate it with diffusion tensor imaging (DTI) experiments. As part of this framework we propose a new spatially regularized maximum likelihood estimator that simultaneously estimates the quantitative parameters and the spatially-smoothly-varying noise level from the acquisitions. The noise level estimation method does not require repeated acquisitions. We show that the amount of regularization in this method can be set a priori to achieve a desired coefficient of variation of the estimated noise level. The noise level estimate allows the construction of a Cramér-Rao-lower-bound based test statistic that reliably assesses the significance of differences between voxels within a scan or across different scans. We show that the regularized noise level estimate improves upon existing methods and results in a substantially increased precision of the uncertainty estimates of the DTI parameters. It enables correct specification of the null distribution of the test statistic and with it the test statistic obtains the highest sensitivity and specificity. The source code of the estimation framework, test statistic and experiment scripts are made available to the community.

Description: Magnetic particle imaging (MPI) shows promise for medical imaging, particularly in angiography of patients with chronic kidney disease. As the first biomedical imaging technique that truly depends on nanoscale materials properties, MPI requires highly optimized magnetic nanoparticle tracers to generate quality images. Until now, researchers have relied on tracers optimized for MRI ${rm T}2^ {ast} $ -weighted imaging that are sub-optimal for MPI. Here, we describe new tracers tailored to MPI's unique physics, synthesized using an organic-phase process and functionalized to ensure biocompatibility and adequate in vivo circulation time. Tailored tracers showed up to 3 $,times$ greater signal-to-noise ratio and better spatial resolution than existing commercial tracers in MPI images of phantoms.

Description: Prospective authors are requested to submit new, unpublished manuscripts for inclusion in the upcoming event described in this call for papers.

Description: Prospective authors are requested to submit new, unpublished manuscripts for inclusion in the upcoming event described in this call for papers.

Description: These instructions give guidelines for preparing papers for this publication. Presents information for authors publishing in this journal.

Description: This page or pages intentionally left blank.

Description: A promising, recently explored, alternative to ultra-high $b$ -value diffusion weighted imaging (UHB-DWI) is apparent ultra-high $b$ -value diffusion-weighted image reconstruction (AUHB-DWR), where a computational model is used to assist in the reconstruction of apparent DW images at ultra-high $b$ -values. Firstly, we present a novel approach to AUHB-DWR that aims to improve image quality. We formulate the reconstruction of an apparent DW image as a hidden conditional random field (HCRF) in which tissue model diffusion parameters act as hidden states in this random field. The second contribution of this paper is a new generation of fully connected conditional random fields, called the hidden stochastically fully connected conditional random fields (HSFCRF) that allows for efficient inference with significantly reduced computational complexity via stochastic clique structures. The proposed AUHB-DWR algorithms, HCRF and HSFCRF, are evaluated quantitatively in nine different patient cases using Fisher's criteria, probability of error, and coefficient of variation metrics to validate its effectiveness for the purpose of improving intensity delineation between expert identified suspected cancerous and healthy tissue within the prostate gland. The proposed methods are also examined using a prostate phantom, where the apparent ultra-high $b$ -value DW images reconstructed using the tested AUHB-DWR methods are compared with real captured UHB-DWI. The results illustrate that the proposed AUHB-DWR methods has improved reconstruction quality and improved intensity delineation compared with existing AUHB-DWR approaches.

Description: A numerical method is shown for calculating the noise correlation coefficient in arrays of magnetic resonance imaging (MRI) coils loaded with capacitively-loaded ring metamaterial lenses, and in the presence of a conducting half-space resembling a sample. This numerical method is validated by comparison with experimental results obtained in two different experimental procedures for double check: noise resistance measurements with a network analyzer and noise correlation measurements in an MRI system. It is found that, for practical array configurations such as overlapping coils or capacitively-decoupled coils, the noise correlation coefficient turns negative for coils loaded with metamaterial lenses. In particular, the analysis is carried out with metamaterial structures known as magnetoinductive lenses, which have been demonstrated in previous works to improve the signal-to-noise ratio of MRI coils. Results are also shown to demonstrate that negative noise correlations have as an effect the improvement of the g-factor in coil arrays for parallel MRI.

Description: Functional magnetic resonance imaging (fMRI) is an indirect measure of neural activity which is modeled as a convolution of the latent neuronal response and the hemodynamic response function (HRF). Since the sources of HRF variability can be nonneural in nature, the measured fMRI signal does not faithfully represent underlying neural activity. Therefore, it is advantageous to deconvolve the HRF from the fMRI signal. However, since both latent neural activity and the voxel-specific HRF is unknown, the deconvolution must be blind. Existing blind deconvolution approaches employ highly parameterized models, and it is unclear whether these models have an over fitting problem. In order to address these issues, we 1) present a nonparametric deconvolution method based on homomorphic filtering to obtain the latent neuronal response from the fMRI signal and, 2) compare our approach to the best performing existing parametric model based on the estimation of the biophysical hemodynamic model using the Cubature Kalman Filter/Smoother. We hypothesized that if the results from nonparametric deconvolution closely resembled that obtained from parametric deconvolution, then the problem of over fitting during estimation in highly parameterized deconvolution models of fMRI could possibly be over stated. Both simulations and experimental results demonstrate support for our hypothesis since the estimated latent neural response from both parametric and nonparametric methods were highly correlated in the visual cortex. Further, simulations showed that both methods were effective in recovering the simulated ground truth of the latent neural response.

Description: Prospective authors are requested to submit new, unpublished manuscripts for inclusion in the upcoming event described in this call for papers.

Description: In this work, we demonstrated that an optimized detector angular configuration based on the anisotropic energy distribution of background scattered X-rays improves X-ray fluorescence computed tomography (XFCT) detection sensitivity. We built an XFCT imaging system composed of a bench-top fluoroscopy X-ray source, a CdTe X-ray detector, and a phantom motion stage. We imaged a 6.4-cm-diameter phantom containing different concentrations of gold solution and investigated the effect of detector angular configuration on XFCT image quality. Based on our previous theoretical study, three detector angles were considered. The X-ray fluorescence detector was first placed at 145 $^{circ}$ (approximating back-scatter) to minimize scatter X-rays. XFCT image quality was compared to images acquired with the detector at 60 $^{circ}$ (forward-scatter) and 90 $^{circ}$ (side-scatter). The datasets for the three different detector positions were also combined to approximate an isotropically arranged detector. The sensitivity was optimized with detector in the 145 $^{circ}$ back-scatter configuration counting the 78-keV gold ${rm K}beta 1$ X-rays. The improvement arose from the reduced energy of scattered X-ray at the 145 $^{circ}$ position and the large energy separation from gold ${rm K} beta 1$ X-rays. The lowest detected concentration in this configuration was 2.5 mgAu/mL (or 0.25% Au with SNR $ = 4.3$ ). This concentration could not be detected with the 60- formula formulatype="inline"〉 $^{circ}$ , 90 $^{circ}$ , or isotropic configurations (SNRs $ = 1.3$ , 0, 2.3, respectively). XFCT imaging dose of 14 mGy was in the range of typical clinical X-ray CT imaging doses. To our knowledge, the sensitivity achieved in this experiment is the highest in any XFCT experiment using an ordinary bench-top X-ray source in a phantom larger than a mouse ( ${〉} 3$ cm).

Description: We have developed a method for the simultaneous estimation of local diffusion and the global fiber tracts based upon the information entropy flow that computes the maximum entropy trajectories between locations and depends upon the global structure of the multi-dimensional and multi-modal diffusion field. Computation of the entropy spectrum pathways requires only solving a simple eigenvector problem for the probability distribution for which efficient numerical routines exist, and a straight forward integration of the probability conservation through ray tracing of the convective modes guided by a global structure of the entropy spectrum coupled with a small scale local diffusion. The intervoxel diffusion is sampled by multi b-shell multi q-angle diffusion weighted imaging data expanded in spherical waves. This novel approach to fiber tracking incorporates global information about multiple fiber crossings in every individual voxel and ranks it in the most scientifically rigorous way. This method has potential significance for a wide range of applications, including studies of brain connectivity.

Description: The variation between images obtained with different scanners or different imaging protocols presents a major challenge in automatic segmentation of biomedical images. This variation especially hampers the application of otherwise successful supervised-learning techniques which, in order to perform well, often require a large amount of labeled training data that is exactly representative of the target data. We therefore propose to use transfer learning for image segmentation. Transfer-learning techniques can cope with differences in distributions between training and target data, and therefore may improve performance over supervised learning for segmentation across scanners and scan protocols. We present four transfer classifiers that can train a classification scheme with only a small amount of representative training data, in addition to a larger amount of other training data with slightly different characteristics. The performance of the four transfer classifiers was compared to that of standard supervised classification on two magnetic resonance imaging brain-segmentation tasks with multi-site data: white matter, gray matter, and cerebrospinal fluid segmentation; and white-matter-/MS-lesion segmentation. The experiments showed that when there is only a small amount of representative training data available, transfer learning can greatly outperform common supervised-learning approaches, minimizing classification errors by up to 60%.

Description: We present a graph-based framework for airway tree reconstruction from computerized tomography (CT) scans and evaluate the performance of different feature categories and their combinations on five lung cohorts. The approach consists of two main processing steps. First, potential airway branch and connection candidates are identified and represented by a graph structure with weighted nodes and edges, respectively. Second, an optimization algorithm is utilized for generating an airway detection result by selecting a subset of airway branches and connections based on graph weights derived from image features. The performance of the algorithm with different feature categories and their combinations was assessed on a set of 50 lung CT scans from five different cohorts, including normal and diseased lungs. Results show trade-offs between feature categories/combinations in terms of correctly (true positive) and incorrectly (false positive) identified airways. Also, the performance of features in dependence of lung cohort was analyzed. Across all cohorts, a good trade-off with high true positive rate (TPR) and low false positive rate (FPR) was achieved by a combination of gray-value, local shape, and structural features. This combination enabled extracting 91.80% of reference airways (TPR) in combination with a low FPR of 1.00%. In addition, this variant was evaluated on the public EXACT'09 test set, and a comparison with other airway detection approaches is provided. One of the main advantages of the presented method is that it is robust against local disturbances/artifacts or other ambiguities that are frequently occurring in lung CT scans.

Description: In this paper, we discuss about applications of different methods for decomposing a signal over elementary waveforms chosen in a family called a dictionary (atomic representations) in optical coherence tomography (OCT). If the representation is learned from the data, a nonparametric dictionary is defined with three fundamental properties of being data-driven, applicability on 3D, and working in multi-scale, which make it appropriate for processing of OCT images. We discuss about application of such representations including complex wavelet based K-SVD, and diffusion wavelets on OCT data. We introduce complex wavelet based K-SVD to take advantage of adaptability in dictionary learning methods to improve the performance of simple dual tree complex wavelets in speckle reduction of OCT datasets in 2D and 3D. The algorithm is evaluated on 144 randomly selected slices from twelve 3D OCTs taken by Topcon 3D OCT-1000 and Cirrus Zeiss Meditec. Improvement of contrast to noise ratio (CNR) (from 0.9 to 11.91 and from 3.09 to 88.9, respectively) is achieved. Furthermore, two approaches are proposed for image segmentation using diffusion. The first method is designing a competition between extended basis functions at each level and the second approach is defining a new distance for each level and clustering based on such distances. A combined algorithm, based on these two methods is then proposed for segmentation of retinal OCTs, which is able to localize 12 boundaries with unsigned border positioning error of $9.22 pm 3.05 mu {hbox {m}}$ , on a test set of 20 slices selected from 13 3D OCTs.

Description: The spectral graph wavelet transform (SGWT) has recently been developed to compute wavelet transforms of functions defined on non-Euclidean spaces such as graphs. By capitalizing on the established framework of the SGWT, we adopt a fast and efficient computation of a discretized Laplace-Beltrami (LB) operator that allows its extension from arbitrary graphs to differentiable and closed 2-D manifolds (smooth surfaces embedded in the 3-D Euclidean space). This particular class of manifolds are widely used in bioimaging to characterize the morphology of cells, tissues, and organs. They are often discretized into triangular meshes, providing additional geometric information apart from simple nodes and weighted connections in graphs. In comparison with the SGWT, the wavelet bases constructed with the LB operator are spatially localized with a more uniform “spread” with respect to underlying curvature of the surface. In our experiments, we first use synthetic data to show that traditional applications of wavelets in smoothing and edge detectio can be done using the wavelet bases constructed with the LB operator. Second, we show that multi-resolutional capabilities of the proposed framework are applicable in the classification of Alzheimer's patients with normal subjects using hippocampal shapes. Wavelet transforms of the hippocampal shape deformations at finer resolutions registered higher sensitivity (96%) and specificity (90%) than the classification results obtained from the direct usage of hippocampal shape deformations. In addition, the Laplace-Beltrami method requires consistently a smaller number of principal components (to retain a fixed variance) at higher resolution as compared to the binary and weighted graph Laplacians, demonstrating the potential of the wavelet bases in adapting to the geometry of the underlying manifold.

Description: Linear predictive models are applied to functional MRI (fMRI) data to estimate boundaries that predict experimental task states for scans. These boundaries are visualized as statistical parametric maps (SPMs) and range from low to high spatial reproducibility across subjects (e.g., Strother , 2004; LaConte , 2003). Such inter-subject pattern reproducibility is an essential characteristic of interpretable SPMs that generalize across subjects. Therefore, we introduce a flexible hybrid model that optimizes reproducibility by simultaneously enhancing the prediction power and reproducibility. This hybrid model is formed by a weighted summation of the optimization functions of a linear discriminate analysis (LDA) model and a generalized canonical correlation (gCCA) model (Afshin-Pour , 2012). LDA preserves the model’s ability to discriminate the fMRI scans of multiple brain states while gCCA finds a linear combination for each subject’s scans such that the estimated boundary map is reproducible. The hybrid model is implemented in a split-half resampling framework (Strother , 2010) which provides reproducibility (r) and prediction (p) quality metrics. Then the model was compared with LDA, and Gaussian Naive Bayes (GNB). For simulated fMRI data, the hybrid model outperforms the other two techniques in terms of receiver operating characteristic (ROC) curves, particularly for detecting less predictable but spatially reproducible networks. These techniques were applied to real fMRI data to estimate the maps for two task contrasts. Our results indicate that compared to LDA and GNB, the hybrid model can provide maps with large increases in reproducibility for small reductions in prediction, which are jointly closer to the ideal performance point of (p=1, r=1).

Description: This paper introduces a new algorithm for recognizing surgical tasks in real-time in a video stream. The goal is to communicate information to the surgeon in due time during a video-monitored surgery. The proposed algorithm is applied to cataract surgery, which is the most common eye surgery. To compensate for eye motion and zoom level variations, cataract surgery videos are first normalized. Then, the motion content of short video subsequences is characterized with spatiotemporal polynomials: a multiscale motion characterization based on adaptive spatiotemporal polynomials is presented. The proposed solution is particularly suited to characterize deformable moving objects with fuzzy borders, which are typically found in surgical videos. Given a target surgical task, the system is trained to identify which spatiotemporal polynomials are usually extracted from videos when and only when this task is being performed. These key spatiotemporal polynomials are then searched in new videos to recognize the target surgical task. For improved performances, the system jointly adapts the spatiotemporal polynomial basis and identifies the key spatiotemporal polynomials using the multiple-instance learning paradigm. The proposed system runs in real-time and outperforms the previous solution from our group, both for surgical task recognition ( $A_z = 0.851$ on average, as opposed to $A_z = 0.794$ previously) and for the joint segmentation and recognition of surgical tasks ( $A_z = 0.856$ on average, as opposed to $A_z = 0.832$ previously).

Description: Electrical impedance tomography (EIT) is used to image the electrical property distribution of a tissue under test. An EIT system comprises complex hardware and software modules, which are typically designed for a specific application. Upgrading these modules is a time-consuming process, and requires rigorous testing to ensure proper functioning of new modules with the existing ones. To this end, we developed a modular and reconfigurable data acquisition (DAQ) system using National Instruments' (NI) hardware and software modules, which offer inherent compatibility over generations of hardware and software revisions. The system can be configured to use up to 32-channels. This EIT system can be used to interchangeably apply current or voltage signal, and measure the tissue response in a semi-parallel fashion. A novel signal averaging algorithm, and 512-point fast Fourier transform (FFT) computation block was implemented on the FPGA. FFT output bins were classified as signal or noise. Signal bins constitute a tissue's response to a pure or mixed tone signal. Signal bins' data can be used for traditional applications, as well as synchronous frequency-difference imaging. Noise bins were used to compute noise power on the FPGA. Noise power represents a metric of signal quality, and can be used to ensure proper tissue-electrode contact. Allocation of these computationally expensive tasks to the FPGA reduced the required bandwidth between PC, and the FPGA for high frame rate EIT. In 16-channel configuration, with a signal-averaging factor of 8, the DAQ frame rate at 100 kHz exceeded 110 frames ${rm s} ^{-1}$ , and signal-to-noise ratio exceeded 90 dB across the spectrum. Reciprocity error was found to be $〈 1%$ for frequencies up to 1 MHz. Static imaging experiments were performed on a high-conductivity inclusion - laced in a saline filled tank; the inclusion was clearly localized in the reconstructions obtained for both absolute current and voltage mode data.

Description: Iterative image reconstruction for positron emission tomography can improve image quality by using spatial regularization. The most commonly used quadratic penalty often oversmoothes sharp edges and fine features in reconstructed images, while nonquadratic penalties can preserve edges and achieve higher contrast recovery. Existing optimization algorithms such as the expectation maximization (EM) and preconditioned conjugate gradient (PCG) algorithms work well for the quadratic penalty, but are less efficient for high-curvature or nonsmooth edge-preserving regularizations. This paper proposes a new algorithm to accelerate edge-preserving image reconstruction by using two strategies: trust surrogate and optimization transfer descent. Trust surrogate approximates the original penalty by a smoother function at each iteration, but guarantees the algorithm to descend monotonically; Optimization transfer descent accelerates a conventional optimization transfer algorithm by using conjugate gradient and line search. Results of computer simulations and real 3-D data show that the proposed algorithm converges much faster than the conventional EM and PCG for smooth edge-preserving regularization and can also be more efficient than the current state-of-art algorithms for the nonsmooth $ell _{1}$ regularization.

Description: Low-dose-rate brachytherapy is a radiation treatment method for localized prostate cancer. The standard of care for this treatment procedure is to acquire transrectal ultrasound images of the prostate in order to devise a plan to deliver sufficient radiation dose to the cancerous tissue. Brachytherapy planning involves delineation of contours in these images, which closely follow the prostate boundary, i.e., clinical target volume. This process is currently performed either manually or semi-automatically, which requires user interaction for landmark initialization. In this paper, we propose a multi-atlas fusion framework to automatically delineate the clinical target volume in ultrasound images. A dataset of a priori segmented ultrasound images, i.e., atlases, is registered to a target image. We introduce a pairwise atlas agreement factor that combines an image-similarity metric and similarity between a priori segmented contours. This factor is used in an atlas selection algorithm to prune the dataset before combining the atlas contours to produce a consensus segmentation. We evaluate the proposed segmentation approach on a set of 280 transrectal prostate volume studies. The proposed method produces segmentation results that are within the range of observer variability when compared to a semi-automatic segmentation technique that is routinely used in our cancer clinic.

Description: In standard B-mode ultrasound (BMUS), segmentation of the lumen of atherosclerotic carotid arteries and studying the lumen geometry over time are difficult owing to irregular lumen shapes, noise, artifacts, and echolucent plaques. Contrast enhanced ultrasound (CEUS) improves lumen visualization, but lumen segmentation remains challenging owing to varying intensities, CEUS-specific artifacts and lack of tissue visualization. To overcome these challenges, we propose a novel method using simultaneously acquired BMUS&CEUS image sequences. Initially, the method estimates nonrigid motion (NME) from the image sequences, using intensity-based image registration. The motion-compensated image sequence is then averaged to obtain a single “epitome” image with improved signal-to-noise ratio. The lumen is segmented from the epitome image through an intensity joint-histogram classification and a graph-based segmentation. NME was validated by comparing displacements with manual annotations in 11 carotids. The average root mean square error (RMSE) was $112pm 73~mu {rm m}$ . Segmentation results were validated against manual delineations in the epitome images of two different datasets, respectively containing 11 (RMSE $191pm 43~mu {rm m}$ ) and 10 (RMSE $351pm 176~mu {rm m}$ ) carotids. From the deformation fields, we derived arterial distensibility with values comparable to the literature. The average errors in all experiments were in the inter-observer variability range. To the best of our knowledge, this is the first study exploiting combined BMUS&CEUS images for atherosclerotic carotid lumen segmentation.

Description: Prospective authors are requested to submit new, unpublished manuscripts for inclusion in the upcoming event described in this call for papers.

Description: This page or pages intentionally left blank.

Description: Multi-channel radio-frequency (RF) transmit coil arrays have been developed to mitigate many of the RF challenges associated with ultra-high field ( $geq$ 7T) magnetic resonance imaging (MRI). These arrays can be used for parallel RF transmission to produce spatially tailored RF excitation over the field of view. However, the realization of such arrays remains a challenge due to significant reactive interaction between the array coils, i.e., mutual coupling. In this paper, a novel bandstop filter (“magnetic wall”) is used in an MRI RF transmit array to decouple individual coils. The proposed decoupling method is inspired by periodic resonator designs commonly used in frequency selective surfaces and is used as a distributed RF filter to suppress the transmission of RF energy between coils in an array. The decoupling of the magnetic wall (MW) is analyzed in terms of equivalent circuits that include terms for both magnetic and electric coupling for an arbitrary number of MW resonant conductors. Both frequency—and time-domain full-wave simulations were performed to analyze a specific MW structure. The performance of the proposed method is experimentally validated for both first-order coupling and higher-order coupling with a three-coil 7T array setup. Analysis and measurements confirm that the rejection band of the MW can be tuned to provide high isolation in the presence of cross coupling between RF array coils.

Description: We develop a new approach to functional brain connectivity analysis, which deals with four fundamental aspects of connectivity not previously jointly treated. These are: temporal correlation, spurious spatial correlation, sparsity, and network construction using trajectory (as opposed to marginal) Mutual Information. We call the new method Sparse Conditional Trajectory Mutual Information (SCoTMI). We demonstrate SCoTMI on simulated and real fMRI data, showing that SCoTMI gives more accurate and more repeatable detection of network links than competing network estimation methods.

Description: The identification of small structures (blobs) from medical images to quantify clinically relevant features, such as size and shape, is important in many medical applications. One particular application explored here is the automated detection of kidney glomeruli after targeted contrast enhancement and magnetic resonance imaging. We propose a computationally efficient algorithm, termed the Hessian-based Difference of Gaussians (HDoG), to segment small blobs (e.g. glomeruli from kidney) from 3D medical images based on local convexity, intensity and shape information. The image is first smoothed and pre-segmented into small blob candidate regions based on local convexity. Two novel 3D regional features (regional blobness and regional flatness) are then extracted from the candidate regions. Together with regional intensity, the three features are used in an unsupervised learning algorithm for auto post-pruning. HDoG is first validated in a 2D form and compared with other three blob detectors from literature, which are generally for 2D images only. To test the detectability of blobs from 3D images, 240 sets of simulated images are rendered for scenarios mimicking the renal nephron distribution observed in contrast-enhanced, 3D MRI. The results show a satisfactory performance of HDoG in detecting large numbers of small blobs. Two sets of real kidney 3D MR images (6 rats, 3 human) are then used to validate the applicability of HDoG for glomeruli detection. By comparing MRI to stereological measurements, we verify that HDoG is a robust and efficient unsupervised technique for 3D blobs segmentation.

Description: We introduce a generative probabilistic model for segmentation of brain lesions in multi-dimensional images that generalizes the EM segmenter, a common approach for modelling brain images using Gaussian mixtures and a probabilistic tissue atlas that employs expectation-maximization (EM), to estimate the label map for a new image. Our model augments the probabilistic atlas of the healthy tissues with a latent atlas of the lesion. We derive an estimation algorithm with closed-form EM update equations. The method extracts a latent atlas prior distribution and the lesion posterior distributions jointly from the image data. It delineates lesion areas individually in each channel, allowing for differences in lesion appearance across modalities, an important feature of many brain tumor imaging sequences. We also propose discriminative model extensions to map the output of the generative model to arbitrary labels with semantic and biological meaning, such as “tumor core” or “fluid-filled structure”, but without a one-to-one correspondence to the hypo- or hyper-intense lesion areas identified by the generative model. We test the approach in two image sets: the publicly available BRATS set of glioma patient scans, and multimodal brain images of patients with acute and subacute ischemic stroke. We find the generative model that has been designed for tumor lesions to generalize well to stroke images, and the extended discriminative -discriminative model to be one of the top ranking methods in the BRATS evaluation.

Description: Watershed segmentation is useful for a number of image segmentation problems with a wide range of practical applications. Traditionally, the tracking of the immersion front is done by applying a fast sorting algorithm. In this work, we explore a continuous approach based on a geometric description of the immersion front which gives rise to a partial differential equation. The main advantage of using a partial differential equation to track the immersion front is that the method becomes versatile and may easily be stabilized by introducing regularization terms. Coupling the geometric approach with a proper “merging strategy” creates a robust algorithm which minimizes over- and under-segmentation even without predefined markers. Since reliable markers defined prior to segmentation can be difficult to construct automatically for various reasons, being able to treat marker-free situations is a major advantage of the proposed method over earlier watershed formulations. The motivation for the methods developed in this paper is taken from high-throughput screening of cells. A fully automated segmentation of single cells enables the extraction of cell properties from large data sets, which can provide substantial insight into a biological model system. Applying smoothing to the boundaries can improve the accuracy in many image analysis tasks requiring a precise delineation of the plasma membrane of the cell. The proposed segmentation method is applied to real images containing fluorescently labeled cells, and the experimental results show that our implementation is robust and reliable for a variety of challenging segmentation tasks.

Description: The major advantage of multiple-instance learning (MIL) applied to a computer-aided detection (CAD) system is that it allows optimizing the latter with case-level labels instead of accurate lesion outlines as traditionally required for a supervised approach. As shown in previous work, a MIL-based CAD system can perform comparably to its supervised counterpart considering complex tasks such as chest radiograph scoring in tuberculosis (TB) detection. However, despite this remarkable achievement, the uncertainty inherent to MIL can lead to a less satisfactory outcome if analysis at lower levels (e.g., regions or pixels) is needed. This issue may seriously compromise the applicability of MIL to tasks related to quantification or grading, or detection of highly localized lesions. In this paper, we propose to reduce uncertainty by embedding a MIL classifier within an active learning (AL) framework. To minimize the labeling effort, we develop a novel instance selection mechanism that exploits the MIL problem definition through one-class classification. We adapt this mechanism to provide meaningful regions instead of individual instances for expert labeling, which is a more appropriate strategy given the application domain. In addition, and contrary to usual AL methods, a single iteration is performed. To show the effectiveness of our approach, we compare the output of a MIL-based CAD system trained with and without the proposed AL framework. The task is to detect textural abnormalities related to TB. Both quantitative and qualitative evaluations at the pixel level are carried out. Our method significantly improves the MIL-based classification.

Description: Radiation therapy is an integral part of cancer treatment, but to date it remains highly manual. Plans are created through optimization of dose volume objectives that specify intent to minimize, maximize, or achieve a prescribed dose level to clinical targets and organs. Optimization is NP-hard, requiring highly iterative and manual initialization procedures. We present a proof-of-concept for a method to automatically infer the radiation dose directly from the patient's treatment planning image based on a database of previous patients with corresponding clinical treatment plans. Our method uses regression forests augmented with density estimation over the most informative features to learn an automatic atlas-selection metric that is tailored to dose prediction. We validate our approach on 276 patients from 3 clinical treatment plan sites (whole breast, breast cavity, and prostate), with an overall dose prediction accuracies of 78.68%, 64.76%, 86.83% under the Gamma metric.

Description: The identification of tumors in the internal organs of chest, abdomen, and pelvis anatomic regions can be performed with the analysis of Dynamic Contrast Enhanced Magnetic Resonance Imaging (DCE-MRI) data. The contrast agent is accumulated differently by pathologic and healthy tissues and that results in a temporally varying contrast in an image series. The internal organs are also subject to potentially extensive movements mainly due to breathing, heart beat, and peristalsis. This contributes to making the analysis of DCE-MRI datasets challenging as well as time consuming. To address this problem we propose a novel pairwise non-rigid registration method with a Non-Parametric Bayesian Registration (NParBR) formulation. The NParBR method uses a Bayesian formulation that assumes a model for the effect of the distortion on the joint intensity statistics, a non-parametric prior for the restored statistics, and also applies a spatial regularization for the estimated registration with Gaussian filtering. A minimally biased intra-dataset atlas is computed for each dataset and used as reference for the registration of the time series. The time series registration method has been tested with 20 datasets of liver, lungs, intestines, and prostate. It has been compared to the B-Splines and to the SyN methods with results that demonstrate that the proposed method improves both accuracy and efficiency.

Description: Statistical image reconstruction (SIR) methods are studied extensively for X-ray computed tomography (CT) due to the potential of acquiring CT scans with reduced X-ray dose while maintaining image quality. However, the longer reconstruction time of SIR methods hinders their use in X-ray CT in practice. To accelerate statistical methods, many optimization techniques have been investigated. Over-relaxation is a common technique to speed up convergence of iterative algorithms. For instance, using a relaxation parameter that is close to two in alternating direction method of multipliers (ADMM) has been shown to speed up convergence significantly. This paper proposes a relaxed linearized augmented Lagrangian (AL) method that shows theoretical faster convergence rate with over-relaxation and applies the proposed relaxed linearized AL method to X-ray CT image reconstruction problems. Experimental results with both simulated and real CT scan data show that the proposed relaxed algorithm (with ordered-subsets [OS] acceleration) is about twice as fast as the existing unrelaxed fast algorithms, with negligible computation and memory overhead.

Description: In this paper, we describe a novel and generic approach to address fully-automatic segmentation of brain tumors by using multi-atlas patch-based voting techniques. In addition to avoiding the local search window assumption, the conventional patch-based framework is enhanced through several simple procedures: an improvement of the training dataset in terms of both label purity and intensity statistics, augmented features to implicitly guide the nearest-neighbor-search, multi-scale patches, invariance to cube isometries, stratification of the votes with respect to cases and labels. A probabilistic model automatically delineates regions of interest enclosing high-probability tumor volumes, which allows the algorithm to achieve highly competitive running time despite minimal processing power and resources. This method was evaluated on Multimodal Brain Tumor Image Segmentation challenge datasets. State-of-the-art results are achieved, with a limited learning stage thus restricting the risk of overfit. Moreover, segmentation smoothness does not involve any post-processing.

Description: The development of an automatic telemedicine system for computer-aided screening and grading of diabetic retinopathy depends on reliable detection of retinal lesions in fundus images. In this paper, a novel method for automatic detection of both microaneurysms and hemorrhages in color fundus images is described and validated. The main contribution is a new set of shape features, called Dynamic Shape Features, that do not require precise segmentation of the regions to be classified. These features represent the evolution of the shape during image flooding and allow to discriminate between lesions and vessel segments. The method is validated per-lesion and per-image using six databases, four of which are publicly available. It proves to be robust with respect to variability in image resolution, quality and acquisition system. On the Retinopathy Online Challenge's database, the method achieves a FROC score of 0.420 which ranks it fourth. On the Messidor database, when detecting images with diabetic retinopathy, the proposed method achieves an area under the ROC curve of 0.899, comparable to the score of human experts, and it outperforms state-of-the-art approaches.

Description: We theoretically demonstrated a new optical imaging technique based on reverse-time migration (RTM) for reconstructing optical structures in homogeneous media for the first time. RTM is a powerful wave-equation-based method to reconstruct the image of the structure by modeling the wave propagation inside the media with both forward modeling and reverse-time extrapolation. While RTM is commonly used with acoustic seismic waves, this paper represents the first effort to develop optical RTM imaging method for biomedical research. To refine the image quality, we further developed new methods to suppress the low-wavenumber artifact (LWA). When compared with the conventional means for LWA suppression such as Laplacian filtering, illumination normalization, and the ratio method, our new derivative-based and power-image methods are able to significantly reduce LWA, resulting in high-quality reconstructed images with sufficient contrasts and spatial resolutions for structure identification. The optical RTM imaging technique may provide a new platform for non-invasive optical imaging of structures in deep layers of tissues for biomedical applications.

Description: Noninvasive cardiac electrophysiological (EP) imaging aims to mathematically reconstruct the spatiotemporal dynamics of cardiac sources from body-surface electrocardiographic (ECG) data. This ill-posed problem is often regularized by a fixed constraining model. However, a fixed-model approach enforces the source distribution to follow a pre-assumed structure that does not always match the varying spatiotemporal distribution of actual sources. To understand the model-data relation and examine the impact of prior models, we present a multiple-model approach for volumetric cardiac EP imaging where multiple prior models are included and automatically picked by the available ECG data. Multiple models are incorporated as an $Lp$ -norm prior for sources, where $p$ is an unknown hyperparameter with a prior uniform distribution. To examine how different combinations of models may be favored by different measurement data, the posterior distribution of cardiac sources and hyperparameter $p$ is calculated using a Markov Chain Monte Carlo (MCMC) technique. The importance of multiple-model prior was assessed in two sets of synthetic and real-data experiments, compared to fixed-model priors (using Laplace and Gaussian priors). The results showed that the posterior combination of models (the posterior distribution of $p$ ) as determined by the ECG data differed substantially when reconstructing sources with different sizes and structures. While the use of fixed models is best suited in situations where the prior assumption fits the actual source structures, the use of an automatically adaptive set of models may have the ability to better address model-data mismatch and to provide consistent performance in reconstructing so- rces with different properties.

Description: Magnetic Resonance—Electrical Properties Tomography (MR-EPT) is an imaging modality that maps the spatial distribution of the electrical conductivity and permittivity using standard MRI systems. The presence of a body within the scanner alters the RF field, and by mapping these alterations it is possible to recover the electrical properties. The field is time-harmonic, and can be described by the Helmholtz equation. Approximations to this equation have been previously used to estimate conductivity and permittivity in terms of first or second derivatives of RF field data. Using these same approximations, an inverse approach to solving the MR-EPT problem is presented here that leverages a forward model for describing the magnitude and phase of the field within the imaging domain, and a fitting approach for estimating the electrical properties distribution. The advantages of this approach are that 1) differentiation of the measured data is not required, thus reducing noise sensitivity, and 2) different regularization schemes can be adopted, depending on prior knowledge of the distribution of conductivity or permittivity, leading to improved image quality. To demonstrate the developed approach, both Quadratic (QR) and Total Variation (TV) regularization methods were implemented and evaluated through numerical simulation and experimentally acquired data. The proposed inverse approach to MR-EPT reconstruction correctly identifies contrasts and accurately reconstructs the geometry in both simulations and experiments. The TV regularized scheme reconstructs sharp spatial transitions, which are difficult to reconstruct with other, more traditional approaches.

Description: The accurate segmentation of lung lesions from computed tomography (CT) scans is important for lung cancer research and can offer valuable information for clinical diagnosis and treatment. However, it is challenging to achieve a fully automatic lesion detection and segmentation with acceptable accuracy due to the heterogeneity of lung lesions. Here, we propose a novel toboggan based growing automatic segmentation approach (TBGA) with a three-step framework, which are automatic initial seed point selection, multi-constraints 3D lesion extraction and the final lesion refinement. The new approach does not require any human interaction or training dataset for lesion detection, yet it can provide a high lesion detection sensitivity (96.35%) and a comparable segmentation accuracy with manual segmentation $({rm P}〉0.05)$ , which was proved by a series assessments using the LIDC-IDRI dataset (850 lesions) and in-house clinical dataset (121 lesions). We also compared TBGA with commonly used level set and skeleton graph cut methods, respectively. The results indicated a significant improvement of segmentation accuracy $({rm P}〈0.05)$ . Furthermore, the average time consumption for one lesion segmentation was under 8 s using our new method. In conclusion, we believe that the novel TBGA can achieve robust, efficient and accurate lung lesion segmentation in CT images automatically.