Robust adaptive 3-D segmentation of vessel laminae from fluorescence confocal microscope images and parallel GPU implementation

Ultraliser: a framework for creating multiscale, high-fidelity and geometrically realistic 3D models for in silico neuroscience

Ultraliser is a neuroscience-specific software framework capable of creating accurate and biologically realistic 3D models of complex neuroscientific structures at intracellular (e.g. mitochondria and endoplasmic reticula), cellular (e.g. neurons and glia) and even multicellular scales of resolution (e.g. cerebral vasculature and minicolumns). Resulting models are exported as triangulated surface meshes and annotated volumes for multiple applications in in silico neuroscience, allowing scalable supercomputer simulations that can unravel intricate cellular structure-function relationships. Ultraliser implements a high-performance and unconditionally robust voxelization engine adapted to create optimized watertight surface meshes and annotated voxel grids from arbitrary non-watertight triangular soups, digitized morphological skeletons or binary volumetric masks. The framework represents a major leap forward in simulation-based neuroscience, making it possible to employ high-resolution 3D structural models for quantification of surface areas and volumes, which are of the utmost importance for cellular and system simulations. The power of Ultraliser is demonstrated with several use cases in which hundreds of models are created for potential application in diverse types of simulations. Ultraliser is publicly released under the GNU GPL3 license on GitHub (BlueBrain/Ultraliser).

SIGNIFICANCE

There is crystal clear evidence on the impact of cell shape on its signaling mechanisms. Structural models can therefore be insightful to realize the function; the more realistic the structure can be, the further we get insights into the function. Creating realistic structural models from existing ones is challenging, particularly when needed for detailed subcellular simulations. We present Ultraliser, a neuroscience-dedicated framework capable of building these structural models with realistic and detailed cellular geometries that can be used for simulations.

From whole-organ imaging to in-silico blood flow modeling: A new multi-scale network analysis for revisiting tissue functional anatomy

We present a multi-disciplinary image-based blood flow perfusion modeling of a whole organ vascular network for analyzing both its structural and functional properties. We show how the use of Light-Sheet Fluorescence Microscopy (LSFM) permits whole-organ micro-vascular imaging, analysis and modelling. By using adapted image post-treatment workflow, we could segment, vectorize and reconstruct the entire micro-vascular network composed of 1.7 million vessels, from the tissue-scale, inside a ∼ 25 × 5 × 1 = 125mm3 volume of the mouse fat pad, hundreds of times larger than previous studies, down to the cellular scale at micron resolution, with the entire blood perfusion modeled. Adapted network analysis revealed the structural and functional organization of meso-scale tissue as strongly connected communities of vessels. These communities share a distinct heterogeneous core region and a more homogeneous peripheral region, consistently with known biological functions of fat tissue. Graph clustering analysis also revealed two distinct robust meso-scale typical sizes (from 10 to several hundred times the cellular size), revealing, for the first time, strongly connected functional vascular communities. These community networks support heterogeneous micro-environments. This work provides the proof of concept that in-silico all-tissue perfusion modeling can reveal new structural and functional exchanges between micro-regions in tissues, found from community clusters in the vascular graph.

Three-Dimensional Segmentation and Reconstruction of Neuronal Nuclei in Confocal Microscopic Images

The detailed architectural examination of the neuronal nuclei in any brain region, using confocal microscopy, requires quantification of fluorescent signals in three-dimensional stacks of confocal images. An essential prerequisite to any quantification is the segmentation of the nuclei which are typically tightly packed in the tissue, the extreme being the hippocampal dentate gyrus (DG), in which nuclei frequently appear to overlap due to limitations in microscope resolution. Segmentation in DG is a challenging task due to the presence of a significant amount of image artifacts and densely packed nuclei. Accordingly, we established an algorithm based on continuous boundary tracing criterion aiming to reconstruct the nucleus surface and to separate the adjacent nuclei. The presented algorithm neither uses a pre-built nucleus model, nor performs image thresholding, which makes it robust against variations in image intensity and poor contrast. Further, the reconstructed surface is used to study morphology and spatial arrangement of the nuclear interior. The presented method is generally dedicated to segmentation of crowded, overlapping objects in 3D space. In particular, it allows us to study quantitatively the architecture of the neuronal nucleus using confocal-microscopic approach.

Preprocessing of Medical Image Data for Three-Dimensional Bioprinted Customized-Neural-Scaffolds

IMPACT STATEMENT

Nerve damage, which can be devastating, triggers several biological cascades, which result in the insufficiencies of the human nervous system to provide complete nerve repair and regain of function. Since no therapeutic strategy exists to provide immediate attention and intervention to patients with newly acquired nerve damage, we propose a strategy in which accelerated medical image processing through graphical processing unit implementation and three-dimensional printing are combined to produce a time-efficient, patient-specific (custom-neural-scaffold) solution to nerve damage. This work aims to beneficially shorten the time required for medical decision-making so that improved patient outcomes are achieved.

Segmentation of Vasculature From Fluorescently Labeled Endothelial Cells in Multi-Photon Microscopy Images

Vasculature is known to be of key biological significance, especially in the study of tumors. As such, considerable effort has been focused on the automated segmentation of vasculature in medical and pre-clinical images. The majority of vascular segmentation methods focus on bloodpool labeling methods; however, particularly, in the study of tumors, it is of particular interest to be able to visualize both the perfused and the non-perfused vasculature. Imaging vasculature by highlighting the endothelium provides a way to separate the morphology of vasculature from the potentially confounding factor of perfusion. Here, we present a method for the segmentation of tumor vasculature in 3D fluorescence microscopic images using signals from the endothelial and surrounding cells. We show that our method can provide complete and semantically meaningful segmentations of complex vasculature using a supervoxel-Markov random field approach. We show that in terms of extracting meaningful segmentations of the vasculature, our method outperforms both state-of-the-art method, specific to these data, as well as more classical vasculature segmentation methods.

Computational 3D imaging to quantify structural components and assembly of protein networks

Traditionally, protein structures have been described by the secondary structure architecture and fold arrangement. However, the relatively novel method of 3D confocal microscopy of fluorescent-protein-tagged networks in living cells allows resolving the detailed spatial organization of these networks. This provides new possibilities to predict network functionality, as structure and function seem to be linked at various scales. Here, we propose a quantitative approach using 3D confocal microscopy image data to describe protein networks based on their nano-structural characteristics. This analysis is constructed in four steps: (i) Segmentation of the microscopic raw data into a volume model and extraction of a spatial graph representing the protein network. (ii) Quantifying protein network gross morphology using the volume model. (iii) Quantifying protein network components using the spatial graph. (iv) Linking these two scales to obtain insights into network assembly. Here, we quantitatively describe the filamentous temperature sensitive Z protein network of the moss Physcomitrella patens and elucidate relations between network size and assembly details. Future applications will link network structure and functionality by tracking dynamic structural changes over time and comparing different states or types of networks, possibly allowing more precise identification of (mal) functions or the design of protein-engineered biomaterials for applications in regenerative medicine.

STATEMENT OF SIGNIFICANCE

Protein networks are highly complex and dynamic structures that play various roles in biological environments. Analyzing the detailed spatial structure of these networks may lead to new insight into biological functions and malfunctions. Here, we propose a tool set that extracts structural information at two scales of the protein network and allows therefore to address questions such as "how is the network built?" or "how networks grow?".

Cerebral vascular structure in the motor cortex of adult mice is stable and is not altered by voluntary exercise

The cerebral vasculature provides blood flow throughout the brain, and local changes in blood flow are regulated to match the metabolic demands of the active brain regions. This neurovascular coupling is mediated by real-time changes in vessel diameter and depends on the underlying vascular network structure. Neurovascular structure is configured during development by genetic and activity-dependent factors. In adulthood, it can be altered by experiences such as prolonged hypoxia, sensory deprivation and seizure. Here, we have sought to determine whether exercise could alter cerebral vascular structure in the adult mouse. We performed repeated in vivo two-photon imaging in the motor cortex of adult transgenic mice expressing membrane-anchored green fluorescent protein in endothelial cells (tyrosine endothelial kinase 2 receptor (Tie2)-Cre:mTmG). This strategy allows for high-resolution imaging of the vessel walls throughout the lifespan. Vascular structure, as measured by capillary branch point number and position, segment diameter and length remained stable over a time scale of months as did pericyte number and position. Furthermore, we compared the vascular structure before, during, and after periods of voluntary wheel running and found no alterations in these same parameters. In both running and control mice, we observed a low rate of capillary segment subtraction. Interestingly, these rare subtraction events preferentially remove short vascular loops.

Real-time GPU-accelerated processing and volumetric display for wide-field laser-scanning optical-resolution photoacoustic microscopy

Fast signal processing and real-time displays are essential for practical imaging modality in various fields of applications. However, the imaging speed in optical-resolution photoacoustic microscopy (OR-PAM), in particular, depends on factors such as the pulse repetition rate of the laser, scanning method, field of view (FOV), and signal processing time. In the past, efforts to increase acquisition speed either focused on developing new scanning methods or using lasers with higher pulse repetition rates. However, high-speed signal processing is also important for real-time volumetric display in OR-PAM. In this study, we carried out parallel signal processing using a graphics processing unit (GPU) to enable fast signal processing and wide-field real-time displays in laser-scanning OR-PAM. The average total GPU processing time for a B-mode PAM image was approximately 1.35 ms at a display speed of 480 fps when the data samples were acquired with 736 (axial) × 500 (lateral) points/B-mode-frame at a pulse repetition rate of 300 kHz. In addition, we successfully displayed maximum amplitude projection images of a mouse's ear as volumetric images with an FOV of 3 mm × 3 mm (500 × 500 pixels) at 1.02 s, corresponding to 0.98 fps.

Computer Vision Techniques for Transcatheter Intervention

Minimally invasive transcatheter technologies have demonstrated substantial promise for the diagnosis and the treatment of cardiovascular diseases. For example, transcatheter aortic valve implantation is an alternative to aortic valve replacement for the treatment of severe aortic stenosis, and transcatheter atrial fibrillation ablation is widely used for the treatment and the cure of atrial fibrillation. In addition, catheter-based intravascular ultrasound and optical coherence tomography imaging of coronary arteries provides important information about the coronary lumen, wall, and plaque characteristics. Qualitative and quantitative analysis of these cross-sectional image data will be beneficial to the evaluation and the treatment of coronary artery diseases such as atherosclerosis. In all the phases (preoperative, intraoperative, and postoperative) during the transcatheter intervention procedure, computer vision techniques (e.g., image segmentation and motion tracking) have been largely applied in the field to accomplish tasks like annulus measurement, valve selection, catheter placement control, and vessel centerline extraction. This provides beneficial guidance for the clinicians in surgical planning, disease diagnosis, and treatment assessment. In this paper, we present a systematical review on these state-of-the-art methods. We aim to give a comprehensive overview for researchers in the area of computer vision on the subject of transcatheter intervention. Research in medical computing is multi-disciplinary due to its nature, and hence, it is important to understand the application domain, clinical background, and imaging modality, so that methods and quantitative measurements derived from analyzing the imaging data are appropriate and meaningful. We thus provide an overview on the background information of the transcatheter intervention procedures, as well as a review of the computer vision techniques and methodologies applied in this area.

Algorithms for improved 3-D reconstruction of live mammalian embryo vasculature from optical coherence tomography data

BACKGROUND

Robust reconstructions of the three-dimensional network of blood vessels in developing embryos imaged by optical coherence tomography (OCT) are needed for quantifying the longitudinal development of vascular networks in live mammalian embryos, in support of developmental cardiovascular research. Past computational methods [such as speckle variance (SV)] have demonstrated the feasibility of vascular reconstruction, but multiple challenges remain including: the presence of vessel structures at multiple spatial scales, thin blood vessels with weak flow, and artifacts resulting from bulk tissue motion (BTM).

METHODS

In order to overcome these challenges, this paper introduces a robust and scalable reconstruction algorithm based on a combination of anomaly detection algorithms and a parametric dictionary based sparse representation of blood vessels from structural OCT data.

RESULTS

Validation results using confocal data as the baseline demonstrate that the proposed method enables the detection of vessel segments that are either partially missed or weakly reconstructed using the SV method. Finally, quantitative measurements of vessel reconstruction quality indicate an overall higher quality of vessel reconstruction with the proposed method.

CONCLUSIONS

Results suggest that sparsity-integrated speckle anomaly detection (SSAD) is potentially a valuable tool for performing accurate quantification of the progression of vascular development in the mammalian embryonic yolk sac as imaged using OCT.

Medical image segmentation on GPUs--a comprehensive review

Segmentation of anatomical structures, from modalities like computed tomography (CT), magnetic resonance imaging (MRI) and ultrasound, is a key enabling technology for medical applications such as diagnostics, planning and guidance. More efficient implementations are necessary, as most segmentation methods are computationally expensive, and the amount of medical imaging data is growing. The increased programmability of graphic processing units (GPUs) in recent years have enabled their use in several areas. GPUs can solve large data parallel problems at a higher speed than the traditional CPU, while being more affordable and energy efficient than distributed systems. Furthermore, using a GPU enables concurrent visualization and interactive segmentation, where the user can help the algorithm to achieve a satisfactory result. This review investigates the use of GPUs to accelerate medical image segmentation methods. A set of criteria for efficient use of GPUs are defined and each segmentation method is rated accordingly. In addition, references to relevant GPU implementations and insight into GPU optimization are provided and discussed. The review concludes that most segmentation methods may benefit from GPU processing due to the methods' data parallel structure and high thread count. However, factors such as synchronization, branch divergence and memory usage can limit the speedup.

Visualization and correction of automated segmentation, tracking and lineaging from 5-D stem cell image sequences

BACKGROUND

Neural stem cells are motile and proliferative cells that undergo mitosis, dividing to produce daughter cells and ultimately generating differentiated neurons and glia. Understanding the mechanisms controlling neural stem cell proliferation and differentiation will play a key role in the emerging fields of regenerative medicine and cancer therapeutics. Stem cell studies in vitro from 2-D image data are well established. Visualizing and analyzing large three dimensional images of intact tissue is a challenging task. It becomes more difficult as the dimensionality of the image data increases to include time and additional fluorescence channels. There is a pressing need for 5-D image analysis and visualization tools to study cellular dynamics in the intact niche and to quantify the role that environmental factors play in determining cell fate.

RESULTS

We present an application that integrates visualization and quantitative analysis of 5-D (x,y,z,t,channel) and large montage confocal fluorescence microscopy images. The image sequences show stem cells together with blood vessels, enabling quantification of the dynamic behaviors of stem cells in relation to their vascular niche, with applications in developmental and cancer biology. Our application automatically segments, tracks, and lineages the image sequence data and then allows the user to view and edit the results of automated algorithms in a stereoscopic 3-D window while simultaneously viewing the stem cell lineage tree in a 2-D window. Using the GPU to store and render the image sequence data enables a hybrid computational approach. An inference-based approach utilizing user-provided edits to automatically correct related mistakes executes interactively on the system CPU while the GPU handles 3-D visualization tasks.

CONCLUSIONS

By exploiting commodity computer gaming hardware, we have developed an application that can be run in the laboratory to facilitate rapid iteration through biological experiments. We combine unsupervised image analysis algorithms with an interactive visualization of the results. Our validation interface allows for each data set to be corrected to 100% accuracy, ensuring that downstream data analysis is accurate and verifiable. Our tool is the first to combine all of these aspects, leveraging the synergies obtained by utilizing validation information from stereo visualization to improve the low level image processing tasks.

Multi-focus image fusion using epifluorescence microscopy for robust vascular segmentation

Automatic segmentation of three-dimensional mi-crovascular structures is needed for quantifying morphological changes to blood vessels during development, disease and treatment processes. Single focus two-dimensional epifluorescent imagery lead to unsatisfactory segmentations due to multiple out of focus vessel regions that have blurred edge structures and lack of detail. Additional segmentation challenges include varying contrast levels due to diffusivity of the lectin stain, leakage out of vessels and fine morphological vessel structure. We propose an approach for vessel segmentation that combines multi-focus image fusion with robust adaptive filtering. The robust adaptive filtering scheme handles noise without destroying small structures, while multi-focus image fusion considerably improves segmentation quality by deblurring out-of-focus regions through incorporating 3D structure information from multiple focus steps. Experiments using epifluorescence images of mice dura mater show an average of 30.4% improvement compared to single focus microvasculature segmentation.

3D actin network centerline extraction with multiple active contours

Fluorescence microscopy is frequently used to study two and three dimensional network structures formed by cytoskeletal polymer fibers such as actin filaments and actin cables. While these cytoskeletal structures are often dilute enough to allow imaging of individual filaments or bundles of them, quantitative analysis of these images is challenging. To facilitate quantitative, reproducible and objective analysis of the image data, we propose a semi-automated method to extract actin networks and retrieve their topology in 3D. Our method uses multiple Stretching Open Active Contours (SOACs) that are automatically initialized at image intensity ridges and then evolve along the centerlines of filaments in the network. SOACs can merge, stop at junctions, and reconfigure with others to allow smooth crossing at junctions of filaments. The proposed approach is generally applicable to images of curvilinear networks with low SNR. We demonstrate its potential by extracting the centerlines of synthetic meshwork images, actin networks in 2D Total Internal Reflection Fluorescence Microscopy images, and 3D actin cable meshworks of live fission yeast cells imaged by spinning disk confocal microscopy. Quantitative evaluation of the method using synthetic images shows that for images with SNR above 5.0, the average vertex error measured by the distance between our result and ground truth is 1 voxel, and the average Hausdorff distance is below 10 voxels.

GPU accelerated segmentation and centerline extraction of tubular structures from medical images

PURPOSE

To create a fast and generic method with sufficient quality for extracting tubular structures such as blood vessels and airways from different modalities (CT, MR and US) and organs (brain, lungs and liver) by utilizing the computational power of graphic processing units (GPUs).

METHODS

A cropping algorithm is used to remove unnecessary data from the datasets on the GPU. A model-based tube detection filter combined with a new parallel centerline extraction algorithm and a parallelized region growing segmentation algorithm is used to extract the tubular structures completely on the GPU. Accuracy of the proposed GPU method and centerline algorithm is compared with the ridge traversal and skeletonization/thinning methods using synthetic vascular datasets.

RESULTS

The implementation is tested on several datasets from three different modalities: airways from CT, blood vessels from MR, and 3D Doppler Ultrasound. The results show that the method is able to extract airways and vessels in 3-5 s on a modern GPU and is less sensitive to noise than other centerline extraction methods.

CONCLUSIONS

Tubular structures such as blood vessels and airways can be extracted from various organs imaged by different modalities in a matter of seconds, even for large datasets.

Three-dimensional vasculature reconstruction of tumour microenvironment via local clustering and classification

The vasculature inside breast cancers is one important component of the tumour microenvironment. The investigation of its spatial morphology, distribution and interactions with cancer cells, including cancer stem cells, is essential for elucidating mechanisms of tumour development and treatment response. Using confocal microscopy and fluorescent markers, we have acquired three-dimensional images of vasculature within mammary tumours and normal mammary gland of mouse models. However, it is difficult to segment and reconstruct complex vasculature accurately from the in vivo three-dimensional images owing to the existence of uneven intensity and regions with low signal-to-noise ratios (SNR). To overcome these challenges, we have developed a novel three-dimensional vasculature segmentation method based on local clustering and classification. First, images of vasculature are clustered into local regions, whose boundaries well delineate vasculature even in low SNR and uneven intensity regions. Then local regions belonging to vasculature are identified by applying a semi-supervised classification method based on three informative features of the local regions. Comparison of results using simulated and real vasculature images, from mouse mammary tumours and normal mammary gland, shows that the new method outperforms existing methods, and can be used for three-dimensional images with uneven background and low SNR to achieve accurate vasculature reconstruction.

Computational tissue volume reconstruction of a peripheral nerve using high-resolution light-microscopy and reconstruct

The development of neural cuff-electrodes requires several in vivo studies and revisions of the electrode design before the electrode is completely adapted to its target nerve. It is therefore favorable to simulate many of the steps involved in this process to reduce costs and animal testing. As the restoration of motor function is one of the most interesting applications of cuff-electrodes, the position and trajectories of myelinated fibers in the simulated nerve are important. In this paper, we investigate a method for building a precise neuroanatomical model of myelinated fibers in a peripheral nerve based on images obtained using high-resolution light microscopy. This anatomical model describes the first aim of our "Virtual workbench" project to establish a method for creating realistic neural simulation models based on image datasets. The imaging, processing, segmentation and technical limitations are described, and the steps involved in the transition into a simulation model are presented. The results showed that the position and trajectories of the myelinated axons were traced and virtualized using our technique, and small nerves could be reliably modeled based on of light microscopy images using low-cost OpenSource software and standard hardware. The anatomical model will be released to the scientific community.

Medical image processing on the GPU - past, present and future

Graphics processing units (GPUs) are used today in a wide range of applications, mainly because they can dramatically accelerate parallel computing, are affordable and energy efficient. In the field of medical imaging, GPUs are in some cases crucial for enabling practical use of computationally demanding algorithms. This review presents the past and present work on GPU accelerated medical image processing, and is meant to serve as an overview and introduction to existing GPU implementations. The review covers GPU acceleration of basic image processing operations (filtering, interpolation, histogram estimation and distance transforms), the most commonly used algorithms in medical imaging (image registration, image segmentation and image denoising) and algorithms that are specific to individual modalities (CT, PET, SPECT, MRI, fMRI, DTI, ultrasound, optical imaging and microscopy). The review ends by highlighting some future possibilities and challenges.

A survey of GPU-based medical image computing techniques

Medical imaging currently plays a crucial role throughout the entire clinical applications from medical scientific research to diagnostics and treatment planning. However, medical imaging procedures are often computationally demanding due to the large three-dimensional (3D) medical datasets to process in practical clinical applications. With the rapidly enhancing performances of graphics processors, improved programming support, and excellent price-to-performance ratio, the graphics processing unit (GPU) has emerged as a competitive parallel computing platform for computationally expensive and demanding tasks in a wide range of medical image applications. The major purpose of this survey is to provide a comprehensive reference source for the starters or researchers involved in GPU-based medical image processing. Within this survey, the continuous advancement of GPU computing is reviewed and the existing traditional applications in three areas of medical image processing, namely, segmentation, registration and visualization, are surveyed. The potential advantages and associated challenges of current GPU-based medical imaging are also discussed to inspire future applications in medicine.

Stab injury and device implantation within the brain results in inversely multiphasic neuroinflammatory and neurodegenerative responses

An estimated 25 million people in the US alone rely on implanted medical devices, ∼2.5 million implanted within the nervous system. Even though many devices perform adequately for years, the host response to medical devices often severely limits tissue integration and long-term performance. This host response is believed to be particularly limiting in the case of intracortical microelectrodes, where it has been shown that glial cell encapsulation and localized neuronal cell loss accompany intracortical microelectrode implantation. Since neuronal ensembles must be within ∼50 µm of the electrode to obtain neuronal spikes and local field potentials, developing a better understanding of the molecular and cellular environment at the device-tissue interface has been the subject of significant research. Unfortunately, immunohistochemical studies of scar maturation in correlation to device function have been inconclusive. Therefore, here we present a detailed quantitative study of the cellular events and the stability of the blood-brain barrier (BBB) following intracortical microelectrode implantation and cortical stab injury in a chronic survival model. We found two distinctly inverse multiphasic profiles for neuronal survival in device-implanted tissue compared to stab-injured animals. For chronically implanted animals, we observed a biphasic paradigm between blood-derived/trauma-induced and CNS-derived inflammatory markers driving neurodegeneration at the interface. In contrast, stab injured animals demonstrated a CNS-mediated neurodegenerative environment. Collectively these data provide valuable insight to the possibility of multiple roles of chronic neuroinflammatory events on BBB disruption and localized neurodegeneration, while also suggesting the importance to consider multiphasic neuroinflammatory kinetics in the design of therapeutic strategies for stabilizing neural interfaces.

A broadly applicable 3-D neuron tracing method based on open-curve snake

This paper presents a broadly applicable algorithm and a comprehensive open-source software implementation for automated tracing of neuronal structures in 3-D microscopy images. The core 3-D neuron tracing algorithm is based on three-dimensional (3-D) open-curve active Contour (Snake). It is initiated from a set of automatically detected seed points. Its evolution is driven by a combination of deforming forces based on the Gradient Vector Flow (GVF), stretching forces based on estimation of the fiber orientations, and a set of control rules. In this tracing model, bifurcation points are detected implicitly as points where multiple snakes collide. A boundariness measure is employed to allow local radius estimation. A suite of pre-processing algorithms enable the system to accommodate diverse neuronal image datasets by reducing them to a common image format. The above algorithms form the basis for a comprehensive, scalable, and efficient software system developed for confocal or brightfield images. It provides multiple automated tracing modes. The user can optionally interact with the tracing system using multiple view visualization, and exercise full control to ensure a high quality reconstruction. We illustrate the utility of this tracing system by presenting results from a synthetic dataset, a brightfield dataset and two confocal datasets from the DIADEM challenge.

GPU computing in medical physics: a review

The graphics processing unit (GPU) has emerged as a competitive platform for computing massively parallel problems. Many computing applications in medical physics can be formulated as data-parallel tasks that exploit the capabilities of the GPU for reducing processing times. The authors review the basic principles of GPU computing as well as the main performance optimization techniques, and survey existing applications in three areas of medical physics, namely image reconstruction, dose calculation and treatment plan optimization, and image processing.