Slic-Seg: A minimally interactive segmentation of the placenta from sparse and motion-corrupted fetal MRI in multiple views

Segmentation of the placenta from fetal MRI is challenging due to sparse acquisition, inter-slice motion, and the widely varying position and shape of the placenta between pregnant women. We propose a minimally interactive framework that combines multiple volumes acquired in different views to obtain accurate segmentation of the placenta. In the first phase, a minimally interactive slice-by-slice propagation method called Slic-Seg is used to obtain an initial segmentation from a single motion-corrupted sparse volume image. It combines high-level features, online Random Forests and Conditional Random Fields, and only needs user interactions in a single slice. In the second phase, to take advantage of the complementary resolution in multiple volumes acquired in different views, we further propose a probability-based 4D Graph Cuts method to refine the initial segmentations using inter-slice and inter-image consistency. We used our minimally interactive framework to examine the placentas of 16 mid-gestation patients from MRI acquired in axial and sagittal views respectively. The results show the proposed method has 1) a good performance even in cases where sparse scribbles provided by the user lead to poor results with the competitive propagation approaches; 2) a good interactivity with low intra- and inter-operator variability; 3) higher accuracy than state-of-the-art interactive segmentation methods; and 4) an improved accuracy due to the co-segmentation based refinement, which outperforms single volume or intensity-based Graph Cuts.

Fully-Automated μMRI Morphometric Phenotyping of the Tc1 Mouse Model of Down Syndrome

We describe a fully automated pipeline for the morphometric phenotyping of mouse brains from μMRI data, and show its application to the Tc1 mouse model of Down syndrome, to identify new morphological phenotypes in the brain of this first transchromosomic animal carrying human chromosome 21. We incorporate an accessible approach for simultaneously scanning multiple ex vivo brains, requiring only a 3D-printed brain holder, and novel image processing steps for their separation and orientation. We employ clinically established multi-atlas techniques–superior to single-atlas methods–together with publicly-available atlas databases for automatic skull-stripping and tissue segmentation, providing high-quality, subject-specific tissue maps. We follow these steps with group-wise registration, structural parcellation and both Voxel- and Tensor-Based Morphometry–advantageous for their ability to highlight morphological differences without the laborious delineation of regions of interest. We show the application of freely available open-source software developed for clinical MRI analysis to mouse brain data: NiftySeg for segmentation and NiftyReg for registration, and discuss atlases and parameters suitable for the preclinical paradigm. We used this pipeline to compare 29 Tc1 brains with 26 wild-type littermate controls, imaged ex vivo at 9.4T. We show an unexpected increase in Tc1 total intracranial volume and, controlling for this, local volume and grey matter density reductions in the Tc1 brain compared to the wild-types, most prominently in the cerebellum, in agreement with human DS and previous histological findings.

Rapid processing of PET list-mode data for efficient uncertainty estimation and data analysis

In this technical note we propose a rapid and scalable software solution for the processing of PET list-mode data, which allows the efficient integration of list mode data processing into the workflow of image reconstruction and analysis. All processing is performed on the graphics processing unit (GPU), making use of streamed and concurrent kernel execution together with data transfers between disk and CPU memory as well as CPU and GPU memory. This approach leads to fast generation of multiple bootstrap realisations, and when combined with fast image reconstruction and analysis, it enables assessment of uncertainties of any image statistic and of any component of the image generation process (e.g. random correction, image processing) within reasonable time frames (e.g. within five minutes per realisation). This is of particular value when handling complex chains of image generation and processing. The software outputs the following: (1) estimate of expected random event data for noise reduction; (2) dynamic prompt and random sinograms of span-1 and span-11 and (3) variance estimates based on multiple bootstrap realisations of (1) and (2) assuming reasonable count levels for acceptable accuracy. In addition, the software produces statistics and visualisations for immediate quality control and crude motion detection, such as: (1) count rate curves; (2) centre of mass plots of the radiodistribution for motion detection; (3) video of dynamic projection views for fast visual list-mode skimming and inspection; (4) full normalisation factor sinograms. To demonstrate the software, we present an example of the above processing for fast uncertainty estimation of regional SUVR (standard uptake value ratio) calculation for a single PET scan of (18)F-florbetapir using the Siemens Biograph mMR scanner.

Music-of-light stethoscope: a demonstration of the photoacoustic effect

In this paper we present a system aimed at demonstrating the photoacoustic (PA) effect for educational purposes. PA imaging is a hybrid imaging modality that requires no contrast agent and has a great potential for spine and brain lesion characterisation, breast cancer and blood flow monitoring notably in the context of fetal surgery. It relies on combining light excitation with ultrasound reception. Our brief was to present and explain PA imaging in a public-friendly way suitable for a variety of ages and backgrounds.

We developed a simple, accessible demonstration unit using readily available materials. We used a modulated light emitting diode (LED) torch and an electronic stethoscope. The output of a music player was used for light modulation and the chest piece of the stethoscope covered by a black tape was used as an absorbing target and an enclosed chamber.

This demonstration unit was presented to the public at the Bloomsbury Festival On Light in October 2015. Our stall was visited by over 100 people of varying ages. Twenty families returned in-depth evaluation questionnaires, which show that our explanations of the photoacoustic effect were well understood. Their interest in biomedical engineering was increased.

Hand-eye calibration for rigid laparoscopes using an invariant point

PURPOSE

Laparoscopic liver resection has significant advantages over open surgery due to less patient trauma and faster recovery times, yet it can be difficult due to the restricted field of view and lack of haptic feedback. Image guidance provides a potential solution but one current challenge is in accurate "hand-eye" calibration, which determines the position and orientation of the laparoscope camera relative to the tracking markers.

METHODS

In this paper, we propose a simple and clinically feasible calibration method based on a single invariant point. The method requires no additional hardware, can be constructed by theatre staff during surgical setup, requires minimal image processing and can be visualised in real time. Real-time visualisation allows the surgical team to assess the calibration accuracy before use in surgery. In addition, in the laboratory, we have developed a laparoscope with an electromagnetic tracking sensor attached to the camera end and an optical tracking marker attached to the distal end. This enables a comparison of tracking performance.

RESULTS

We have evaluated our method in the laboratory and compared it to two widely used methods, "Tsai's method" and "direct" calibration. The new method is of comparable accuracy to existing methods, and we show RMS projected error due to calibration of 1.95 mm for optical tracking and 0.85 mm for EM tracking, versus 4.13 and 1.00 mm respectively, using existing methods. The new method has also been shown to be workable under sterile conditions in the operating room.

CONCLUSION

We have proposed a new method of hand-eye calibration, based on a single invariant point. Initial experience has shown that the method provides visual feedback, satisfactory accuracy and can be performed during surgery. We also show that an EM sensor placed near the camera would provide significantly improved image overlay accuracy.

Maximum-likelihood joint image reconstruction and motion estimation with misaligned attenuation in TOF-PET/CT

This work is an extension of our recent work on joint activity reconstruction/motion estimation (JRM) from positron emission tomography (PET) data. We performed JRM by maximization of the penalized log-likelihood in which the probabilistic model assumes that the same motion field affects both the activity distribution and the attenuation map. Our previous results showed that JRM can successfully reconstruct the activity distribution when the attenuation map is misaligned with the PET data, but converges slowly due to the significant cross-talk in the likelihood. In this paper, we utilize time-of-flight PET for JRM and demonstrate that the convergence speed is significantly improved compared to JRM with conventional PET data.

Maximum-Likelihood Joint Image Reconstruction/Motion Estimation in Attenuation-Corrected Respiratory Gated PET/CT Using a Single Attenuation Map

This work provides an insight into positron emission tomography (PET) joint image reconstruction/motion estimation (JRM) by maximization of the likelihood, where the probabilistic model accounts for warped attenuation. Our analysis shows that maximum-likelihood (ML) JRM returns the same reconstructed gates for any attenuation map (μ-map) that is a deformation of a given μ-map, regardless of its alignment with the PET gates. We derived a joint optimization algorithm accordingly, and applied it to simulated and patient gated PET data. We first evaluated the proposed algorithm on simulations of respiratory gated PET/CT data based on the XCAT phantom. Our results show that independently of which μ-map is used as input to JRM: (i) the warped μ-maps correspond to the gated μ-maps, (ii) JRM outperforms the traditional post-registration reconstruction and consolidation (PRRC) for hot lesion quantification and (iii) reconstructed gated PET images are similar to those obtained with gated μ-maps. This suggests that a breath-held μ-map can be used. We then applied JRM on patient data with a μ-map derived from a breath-held high resolution CT (HRCT), and compared the results with PRRC, where each reconstructed PET image was obtained with a corresponding cine-CT gated μ-map. Results show that JRM with breath-held HRCT achieves similar reconstruction to that using PRRC with cine-CT. This suggests a practical low-dose solution for implementation of motion-corrected respiratory gated PET/CT.

Voxelwise atlas rating for computer assisted diagnosis: Application to congenital heart diseases of the great arteries

Atlas-based analysis methods rely on the morphological similarity between the atlas and target images, and on the availability of labelled images. Problems can arise when the deformations introduced by pathologies affect the similarity between the atlas and a patient's image. The aim of this work is to exploit the morphological dissimilarities between atlas databases and pathological images to diagnose the underlying clinical condition, while avoiding the dependence on labelled images. We propose a voxelwise atlas rating approach (VoxAR) relying on multiple atlas databases, each representing a particular condition. Using a local image similarity measure to assess the morphological similarity between the atlas and target images, a rating map displaying for each voxel the condition of the atlases most similar to the target is defined. The final diagnosis is established by assigning the condition of the database the most represented in the rating map. We applied the method to diagnose three different conditions associated with dextro-transposition of the great arteries, a congenital heart disease. The proposed approach outperforms other state-of-the-art methods using annotated images, with an accuracy of 97.3% when evaluated on a set of 60 whole heart MR images containing healthy and pathological subjects using cross validation.

Locally rigid, vessel-based registration for laparoscopic liver surgery

PURPOSE

Laparoscopic liver resection has significant advantages over open surgery due to less patient trauma and faster recovery times, yet is difficult for most lesions due to the restricted field of view and lack of haptic feedback. Image guidance provides a potential solution but is challenging in a soft deforming organ such as the liver. In this paper, we therefore propose a laparoscopic ultrasound (LUS) image guidance system and study the feasibility of a locally rigid registration for laparoscopic liver surgery.

METHODS

We developed a real-time segmentation method to extract vessel centre points from calibrated, freehand, electromagnetically tracked, 2D LUS images. Using landmark-based initial registration and an optional iterative closest point (ICP) point-to-line registration, a vessel centre-line model extracted from preoperative computed tomography (CT) is registered to the ultrasound data during surgery.

RESULTS

Using the locally rigid ICP method, the RMS residual error when registering to a phantom was 0.7 mm, and the mean target registration error (TRE) for two in vivo porcine studies was 3.58 and 2.99 mm, respectively. Using the locally rigid landmark-based registration method gave a mean TRE of 4.23 mm using vessel centre lines derived from CT scans taken with pneumoperitoneum and 6.57 mm without pneumoperitoneum.

CONCLUSION

In this paper we propose a practical image-guided surgery system based on locally rigid registration of a CT-derived model to vascular structures located with LUS. In a physical phantom and during porcine laparoscopic liver resection, we demonstrate accuracy of target location commensurate with surgical requirements. We conclude that locally rigid registration could be sufficient for practically useful image guidance in the near future.

Estimating anatomical trajectories with Bayesian mixed-effects modeling

We introduce a mass-univariate framework for the analysis of whole-brain structural trajectories using longitudinal Voxel-Based Morphometry data and Bayesian inference. Our approach to developmental and aging longitudinal studies characterizes heterogeneous structural growth/decline between and within groups. In particular, we propose a probabilistic generative model that parameterizes individual and ensemble average changes in brain structure using linear mixed-effects models of age and subject-specific covariates. Model inversion uses Expectation Maximization (EM), while voxelwise (empirical) priors on the size of individual differences are estimated from the data. Bayesian inference on individual and group trajectories is realized using Posterior Probability Maps (PPM). In addition to parameter inference, the framework affords comparisons of models with varying combinations of model order for fixed and random effects using model evidence. We validate the model in simulations and real MRI data from the Alzheimer's Disease Neuroimaging Initiative (ADNI) project. We further demonstrate how subject specific characteristics contribute to individual differences in longitudinal volume changes in healthy subjects, Mild Cognitive Impairment (MCI), and Alzheimer's Disease (AD).

USING FLORBETAPIR PET TO INCREASE DIAGNOSTIC CERTAINTY IN ATYPICAL DEMENTIAS

Amyloid PET or CSF can be used to determine Alzheimer pathology in vivo. Few studies have assessed the additional value of amyloid imaging where CSF results are equivocal. We recruited 20 cognitive patients (65.5+/–7.6 y) with MRI, neuropsychology, and CSF Aβ1–42 and tau measured during their diagnostic assessment. Individuals were selected to have a range of CSF results; ten had amnestic and ten non-amnestic presentations. Following the investigations, the treating neurologist gave a diagnosis (AD or non-AD). Four controls (63+/–7.0y) also had CSF examination. All subjects had Florbetapir PET imaging, reported as positive/negative. The clinicians were given the PET results and asked to review their diagnoses. Eighteen patients had positive Florbetapir scans; two patients and all controls were Florbetapir negative. Following initial investigations, thirteen patients were diagnosed with AD, and seven with non-AD pathology. Providing the Florbetapir result led to a change in diagnosis in seven patients, five of whom had atypical phenotypes. For all seven the CSF results were close to or in a “grey” area, where results overlapped for positive and negative PET scans. Even in individuals with CSF measures of Aβ1–42, and tau, Florbetapir PET imaging may have diagnostic utility, particularly in atypical cases and/or equivocal CSF results.

Interventional Photoacoustic Imaging of the Human Placenta with Ultrasonic Tracking for Minimally Invasive Fetal Surgeries

Image guidance plays a central role in minimally invasive fetal surgery such as photocoagulation of inter-twin placental anastomosing vessels to treat twin-to-twin transfusion syndrome (TTTS). Fetoscopic guidance provides insufficient sensitivity for imaging the vasculature that lies beneath the fetal placental surface due to strong light scattering in biological tissues. Incomplete photocoagulation of anastamoses is associated with postoperative complications and higher perinatal mortality. In this study, we investigated the use of multi-spectral photoacoustic (PA) imaging for better visualization of the placental vasculature. Excitation light was delivered with an optical fiber with dimensions that are compatible with the working channel of a fetoscope. Imaging was performed on an ex vivo normal term human placenta collected at Caesarean section birth. The photoacoustically-generated ultrasound signals were received by an external clinical linear array ultrasound imaging probe. A vein under illumination on the fetal placenta surface was visualized with PA imaging, and good correspondence was obtained between the measured PA spectrum and the optical absorption spectrum of deoxygenated blood. The delivery fiber had an attached fiber optic ultrasound sensor positioned directly adjacent to it, so that its spatial position could be tracked by receiving transmissions from the ultrasound imaging probe. This study provides strong indications that PA imaging in combination with ultrasonic tracking could be useful for detecting the human placental vasculature during minimally invasive fetal surgery.

Interventional Photoacoustic Imaging of the Human Placenta with Ultrasonic Tracking for Minimally Invasive Fetal Surgeries

Image guidance plays a central role in minimally invasive fetal surgery such as photocoagulation of inter-twin placental anastomosing vessels to treat twin-to-twin transfusion syndrome (TTTS). Fetoscopic guidance provides insufficient sensitivity for imaging the vasculature that lies beneath the fetal placental surface due to strong light scattering in biological tissues. Incomplete photocoagulation of anastamoses is associated with postoperative complications and higher perinatal mortality. In this study, we investigated the use of multi-spectral photoacoustic (PA) imaging for better visualization of the placental vasculature. Excitation light was delivered with an optical fiber with dimensions that are compatible with the working channel of a fetoscope. Imaging was performed on an ex vivo normal term human placenta collected at Caesarean section birth. The photoacoustically-generated ultrasound signals were received by an external clinical linear array ultrasound imaging probe. A vein under illumination on the fetal placenta surface was visualized with PA imaging, and good correspondence was obtained between the measured PA spectrum and the optical absorption spectrum of deoxygenated blood. The delivery fiber had an attached fiber optic ultrasound sensor positioned directly adjacent to it, so that its spatial position could be tracked by receiving transmissions from the ultrasound imaging probe. This study provides strong indications that PA imaging in combination with ultrasonic tracking could be useful for detecting the human placental vasculature during minimally invasive fetal surgery.

Toward adaptive radiotherapy for head and neck patients: Uncertainties in dose warping due to the choice of deformable registration algorithm

PURPOSE

The aims of this work were to evaluate the performance of several deformable image registration (DIR) algorithms implemented in our in-house software (NiftyReg) and the uncertainties inherent to using different algorithms for dose warping.

METHODS

The authors describe a DIR based adaptive radiotherapy workflow, using CT and cone-beam CT (CBCT) imaging. The transformations that mapped the anatomy between the two time points were obtained using four different DIR approaches available in NiftyReg. These included a standard unidirectional algorithm and more sophisticated bidirectional ones that encourage or ensure inverse consistency. The forward (CT-to-CBCT) deformation vector fields (DVFs) were used to propagate the CT Hounsfield units and structures to the daily geometry for "dose of the day" calculations, while the backward (CBCT-to-CT) DVFs were used to remap the dose of the day onto the planning CT (pCT). Data from five head and neck patients were used to evaluate the performance of each implementation based on geometrical matching, physical properties of the DVFs, and similarity between warped dose distributions. Geometrical matching was verified in terms of dice similarity coefficient (DSC), distance transform, false positives, and false negatives. The physical properties of the DVFs were assessed calculating the harmonic energy, determinant of the Jacobian, and inverse consistency error of the transformations. Dose distributions were displayed on the pCT dose space and compared using dose difference (DD), distance to dose difference, and dose volume histograms.

RESULTS

All the DIR algorithms gave similar results in terms of geometrical matching, with an average DSC of 0.85 ± 0.08, but the underlying properties of the DVFs varied in terms of smoothness and inverse consistency. When comparing the doses warped by different algorithms, we found a root mean square DD of 1.9% ± 0.8% of the prescribed dose (pD) and that an average of 9% ± 4% of voxels within the treated volume failed a 2%pD DD-test (DD2%-pp). Larger DD2%-pp was found within the high dose gradient (21% ± 6%) and regions where the CBCT quality was poorer (28% ± 9%). The differences when estimating the mean and maximum dose delivered to organs-at-risk were up to 2.0%pD and 2.8%pD, respectively.

CONCLUSIONS

The authors evaluated several DIR algorithms for CT-to-CBCT registrations. In spite of all methods resulting in comparable geometrical matching, the choice of DIR implementation leads to uncertainties in dose warped, particularly in regions of high gradient and/or poor imaging quality.

Simulated field maps for susceptibility artefact correction in interventional MRI

PURPOSE

Intraoperative MRI (iMRI) is a powerful modality for acquiring images of the brain to facilitate precise image-guided neurosurgery. Diffusion-weighted MRI (DW-MRI) provides critical information about location, orientation and structure of nerve fibre tracts, but suffers from the "susceptibility artefact" stemming from magnetic field perturbations due to the step change in magnetic susceptibility at air-tissue boundaries in the head. An existing approach to correcting the artefact is to acquire a field map by means of an additional MRI scan. However, to recover true field maps from the acquired field maps near air-tissue boundaries is challenging, and acquired field maps are unavailable in historical MRI data sets. This paper reports a detailed account of our method to simulate field maps from structural MRI scans that was first presented at IPCAI 2014 and provides a thorough experimental and analysis section to quantitatively validate our technique.

METHODS

We perform automatic air-tissue segmentation of intraoperative MRI scans, feed the segmentation into a field map simulation step and apply the acquired and the simulated field maps to correct DW-MRI data sets.

RESULTS

We report results for 12 patient data sets acquired during anterior temporal lobe resection surgery for the surgical management of focal epilepsy. We find a close agreement between acquired and simulated field maps and observe a statistically significant reduction in the susceptibility artefact in DW-MRI data sets corrected using simulated field maps in the vicinity of the resection. The artefact reduction obtained using acquired field maps remains better than that using the simulated field maps in all evaluated regions of the brain.

CONCLUSIONS

The proposed simulated field maps facilitate susceptibility artefact reduction near the resection. Accurate air-tissue segmentation is key to achieving accurate simulation. The proposed simulation approach is adaptable to different iMRI and neurosurgical applications.

Practical PET Respiratory Motion Correction in Clinical PET/MR

Respiratory motion during PET acquisition may lead to blurring in resulting images and underestimation of uptake parameters. The advent of integrated PET/MR scanners allows us to exploit the integration of modalities, using high spatial resolution and high-contrast MR images to monitor and correct PET images degraded by motion. We proposed a practical, anatomy-independent MR-based correction strategy for PET data affected by respiratory motion and showed that it can improve image quality both for PET acquired simultaneously to the motion-capturing MR and for PET acquired up to 1 h earlier during a clinical scan.

METHODS

To estimate the respiratory motion, our method needs only an extra 1-min dynamic MR scan, acquired at the end of the clinical PET/MR protocol. A respiratory signal was extracted directly from the PET list-mode data. This signal was used to gate the PET data and to construct a motion model built from the dynamic MR data. The estimated motion was then incorporated into the PET image reconstruction to obtain a single motion-corrected PET image. We evaluated our method in 2 steps. The PET-derived respiratory signal was compared with an MR measure of diaphragmatic displacement via a pencil-beam navigator. The motion-corrected images were compared with uncorrected images with visual inspection, line profiles, and standardized uptake value (SUV) in focally avid lesions.

RESULTS

We showed a strong correlation between the PET-derived and MR-derived respiratory signals for 9 patients, with a mean correlation of 0.89. We then showed 4 clinical case study examples ((18)F-FDG and (68)Ga-DOTATATE) using the motion-correction technique, demonstrating improvements in image sharpness and reduction of respiratory artifacts in scans containing pancreatic, liver, and lung lesions as well as cardiac scans. The mean increase in peak SUV (SUV(peak)) and maximum SUV (SUV(max)) in a patient with 4 pancreatic lesions was 23.1% and 34.5% in PET acquired simultaneously with motion-capturing MR, and 17.6% and 24.7% in PET acquired 50 min before as part of the clinical scan.

CONCLUSION

We showed that a respiratory signal can be obtained from raw PET data and that the clinical PET image quality can be improved using only a short additional PET/MR acquisition. Our method does not need external respiratory hardware or modification of the normal clinical MR sequences.

In vivo imaging of tau pathology using multi-parametric quantitative MRI

As the number of people diagnosed with Alzheimer's disease (AD) reaches epidemic proportions, there is an urgent need to develop effective treatment strategies to tackle the social and economic costs of this fatal condition. Dozens of candidate therapeutics are currently being tested in clinical trials, and compounds targeting the aberrant accumulation of tau proteins into neurofibrillary tangles (NFTs) are the focus of substantial current interest. Reliable, translatable biomarkers sensitive to both tau pathology and its modulation by treatment along with animal models that faithfully reflect aspects of the human disease are urgently required. Magnetic resonance imaging (MRI) is well established as a valuable tool for monitoring the structural brain changes that accompany AD progression. However the descent into dementia is not defined by macroscopic brain matter loss alone: non-invasive imaging measurements sensitive to protein accumulation, white matter integrity and cerebral haemodynamics probe distinct aspects of AD pathophysiology and may serve as superior biomarkers for assessing drug efficacy. Here we employ a multi-parametric array of five translatable MRI techniques to characterise the in vivo pathophysiological phenotype of the rTg4510 mouse model of tauopathy (structural imaging, diffusion tensor imaging (DTI), arterial spin labelling (ASL), chemical exchange saturation transfer (CEST) and glucose CEST). Tau-induced pathological changes included grey matter atrophy, increased radial diffusivity in the white matter, decreased amide proton transfer and hyperperfusion. We demonstrate that the above markers unambiguously discriminate between the transgenic group and age-matched controls and provide a comprehensive profile of the multifaceted neuropathological processes underlying the rTg4510 model. Furthermore, we show that ASL and DTI techniques offer heightened sensitivity to processes believed to precede detectable structural changes and, as such, provides a platform for the study of disease mechanisms and therapeutic intervention.

A modality-agnostic patch-based technique for lesion filling in multiple sclerosis

Multiple Sclerosis lesions influence the process of image analysis, leading to tissue segmentation problems and biased morphometric estimates. With the aim of reducing this bias, existing techniques fill segmented lesions as normal appearing white matter. However, due to lesion segmentation errors or the presence of neighbouring structures, such as the ventricles and deep grey matter structures, filling all lesions as white matter like intensities is prone to introduce errors and artefacts. In this paper, we present a novel lesion filling strategy based on in-painting techniques for image completion. This technique makes use of a patch-based Non-Local Means algorithm that fills the lesions with the most plausible texture, rather than normal appearing white matter. We demonstrate that this strategy introduces less bias and fewer artefacts and spurious edges than previous techniques. The advantages of the proposed methodology are that it preserves both anatomical structure and signal-to-noise characteristics even when the lesions are neighbouring grey matter and cerebrospinal fluid, and avoids excess blurring or rasterisation due to the choice of segmentation plane, and lesion shape, size and/or position.

The NifTK software platform for image-guided interventions: platform overview and NiftyLink messaging

PURPOSE

To perform research in image-guided interventions, researchers need a wide variety of software components, and assembling these components into a flexible and reliable system can be a challenging task. In this paper, the NifTK software platform is presented. A key focus has been high-performance streaming of stereo laparoscopic video data, ultrasound data and tracking data simultaneously.

METHODS

A new messaging library called NiftyLink is introduced that uses the OpenIGTLink protocol and provides the user with easy-to-use asynchronous two-way messaging, high reliability and comprehensive error reporting. A small suite of applications called NiftyGuide has been developed, containing lightweight applications for grabbing data, currently from position trackers and ultrasound scanners. These applications use NiftyLink to stream data into NiftyIGI, which is a workstation-based application, built on top of MITK, for visualisation and user interaction. Design decisions, performance characteristics and initial applications are described in detail. NiftyLink was tested for latency when transmitting images, tracking data, and interleaved imaging and tracking data.

RESULTS

NiftyLink can transmit tracking data at 1,024 frames per second (fps) with latency of 0.31 milliseconds, and 512 KB images with latency of 6.06 milliseconds at 32 fps. NiftyIGI was tested, receiving stereo high-definition laparoscopic video at 30 fps, tracking data from 4 rigid bodies at 20-30 fps and ultrasound data at 20 fps with rendering refresh rates between 2 and 20 Hz with no loss of user interaction.

CONCLUSION

These packages form part of the NifTK platform and have proven to be successful in a variety of image-guided surgery projects. Code and documentation for the NifTK platform are available from http://www.niftk.org . NiftyLink is provided open-source under a BSD license and available from http://github.com/NifTK/NiftyLink . The code for this paper is tagged IJCARS-2014.

Toward adaptive radiotherapy for head and neck patients: Feasibility study on using CT-to-CBCT deformable registration for "dose of the day" calculations

PURPOSE

The aim of this study was to evaluate the appropriateness of using computed tomography (CT) to cone-beam CT (CBCT) deformable image registration (DIR) for the application of calculating the "dose of the day" received by a head and neck patient.

METHODS

NiftyReg is an open-source registration package implemented in our institution. The affine registration uses a Block Matching-based approach, while the deformable registration is a GPU implementation of the popular B-spline Free Form Deformation algorithm. Two independent tests were performed to assess the suitability of our registrations methodology for "dose of the day" calculations in a deformed CT. A geometric evaluation was performed to assess the ability of the DIR method to map identical structures between the CT and CBCT datasets. Features delineated in the planning CT were warped and compared with features manually drawn on the CBCT. The authors computed the dice similarity coefficient (DSC), distance transformation, and centre of mass distance between features. A dosimetric evaluation was performed to evaluate the clinical significance of the registrations errors in the application proposed and to identify the limitations of the approximations used. Dose calculations for the same intensity-modulated radiation therapy plan on the deformed CT and replan CT were compared. Dose distributions were compared in terms of dose differences (DD), gamma analysis, target coverage, and dose volume histograms (DVHs). Doses calculated in a rigidly aligned CT and directly in an extended CBCT were also evaluated.

RESULTS

A mean value of 0.850 in DSC was achieved in overlap between manually delineated and warped features, with the distance between surfaces being less than 2 mm on over 90% of the pixels. Deformable registration was clearly superior to rigid registration in mapping identical structures between the two datasets. The dose recalculated in the deformed CT is a good match to the dose calculated on a replan CT. The DD is smaller than 2% of the prescribed dose on 90% of the body's voxels and it passes a 2% and 2 mm gamma-test on over 95% of the voxels. Target coverage similarity was assessed in terms of the 95%-isodose volumes. A mean value of 0.962 was obtained for the DSC, while the distance between surfaces is less than 2 mm in 95.4% of the pixels. The method proposed provided adequate dose estimation, closer to the gold standard than the other two approaches. Differences in DVH curves were mainly due to differences in the OARs definition (manual vs warped) and not due to differences in dose estimation (dose calculated in replan CT vs dose calculated in deformed CT).

CONCLUSIONS

Deforming a planning CT to match a daily CBCT provides the tools needed for the calculation of the "dose of the day" without the need to acquire a new CT. The initial clinical application of our method will be weekly offline calculations of the "dose of the day," and use this information to inform adaptive radiotherapy (ART). The work here presented is a first step into a full implementation of a "dose-driven" online ART.

Global image registration using a symmetric block-matching approach

Most medical image registration algorithms suffer from a directionality bias that has been shown to largely impact subsequent analyses. Several approaches have been proposed in the literature to address this bias in the context of nonlinear registration, but little work has been done for global registration. We propose a symmetric approach based on a block-matching technique and least-trimmed square regression. The proposed method is suitable for multimodal registration and is robust to outliers in the input images. The symmetric framework is compared with the original asymmetric block-matching technique and is shown to outperform it in terms of accuracy and robustness. The methodology presented in this article has been made available to the community as part of the NiftyReg open-source package.

Automatic structural parcellation of mouse brain MRI using multi-atlas label fusion

Multi-atlas segmentation propagation has evolved quickly in recent years, becoming a state-of-the-art methodology for automatic parcellation of structural images. However, few studies have applied these methods to preclinical research. In this study, we present a fully automatic framework for mouse brain MRI structural parcellation using multi-atlas segmentation propagation. The framework adopts the similarity and truth estimation for propagated segmentations (STEPS) algorithm, which utilises a locally normalised cross correlation similarity metric for atlas selection and an extended simultaneous truth and performance level estimation (STAPLE) framework for multi-label fusion. The segmentation accuracy of the multi-atlas framework was evaluated using publicly available mouse brain atlas databases with pre-segmented manually labelled anatomical structures as the gold standard, and optimised parameters were obtained for the STEPS algorithm in the label fusion to achieve the best segmentation accuracy. We showed that our multi-atlas framework resulted in significantly higher segmentation accuracy compared to single-atlas based segmentation, as well as to the original STAPLE framework.

Multi-modal measurement of the myelin-to-axon diameter g-ratio in preterm-born neonates and adult controls

Infants born prematurely are at increased risk of adverse functional outcome. The measurement of white matter tissue composition and structure can help predict functional performance and this motivates the search for new multi-modal imaging biomarkers. In this work we develop a novel combined biomarker from diffusion MRI and multi-component T2 relaxation measurements in a group of infants born very preterm and scanned between 30 and 40 weeks equivalent gestational age. We also investigate this biomarker on a group of seven adult controls, using a multi-modal joint model-fitting strategy. The proposed emergent biomarker is tentatively related to axonal energetic efficiency (in terms of axonal membrane charge storage) and conduction velocity and is thus linked to the tissue electrical properties, giving it a good theoretical justification as a predictive measurement of functional outcome.

The pattern of atrophy in familial Alzheimer disease: volumetric MRI results from the DIAN study

OBJECTIVE

To assess regional patterns of gray and white matter atrophy in familial Alzheimer disease (FAD) mutation carriers.

METHODS

A total of 192 participants with volumetric T1-weighted MRI, genotyping, and clinical diagnosis were available from the Dominantly Inherited Alzheimer Network. Of these, 69 were presymptomatic mutation carriers, 50 were symptomatic carriers (31 with Clinical Dementia Rating [CDR] = 0.5, 19 with CDR > 0.5), and 73 were noncarriers from the same families. Voxel-based morphometry was used to identify cross-sectional group differences in gray matter and white matter volume.

RESULTS

Significant differences in gray matter (p < 0.05, family-wise error-corrected) were observed between noncarriers and mildly symptomatic (CDR = 0.5) carriers in the thalamus and putamen, as well as in the temporal lobe, precuneus, and cingulate gyrus; the same pattern, but with more extensive changes, was seen in those with CDR > 0.5. Significant white matter differences between noncarriers and symptomatic carriers were observed in the cingulum and fornix; these form input and output connections to the medial temporal lobe, cingulate, and precuneus. No differences between noncarriers and presymptomatic carriers survived correction for multiple comparisons, but there was a trend for decreased gray matter in the thalamus for carriers closer to their estimated age at onset. There were no significant increases of gray or white matter in asymptomatic or symptomatic carriers compared to noncarriers.

CONCLUSIONS

Atrophy in FAD is observed early, both in areas commonly associated with sporadic Alzheimer disease and also in the putamen and thalamus, 2 regions associated with early amyloid deposition in FAD mutation carriers.

Using manifold learning for atlas selection in multi-atlas segmentation

Multi-atlas segmentation has been widely used to segment various anatomical structures. The success of this technique partly relies on the selection of atlases that are best mapped to a new target image after registration. Recently, manifold learning has been proposed as a method for atlas selection. Each manifold learning technique seeks to optimize a unique objective function. Therefore, different techniques produce different embeddings even when applied to the same data set. Previous studies used a single technique in their method and gave no reason for the choice of the manifold learning technique employed nor the theoretical grounds for the choice of the manifold parameters. In this study, we compare side-by-side the results given by 3 manifold learning techniques (Isomap, Laplacian Eigenmaps and Locally Linear Embedding) on the same data set. We assess the ability of those 3 different techniques to select the best atlases to combine in the framework of multi-atlas segmentation. First, a leave-one-out experiment is used to optimize our method on a set of 110 manually segmented atlases of hippocampi and find the manifold learning technique and associated manifold parameters that give the best segmentation accuracy. Then, the optimal parameters are used to automatically segment 30 subjects from the Alzheimer's Disease Neuroimaging Initiative (ADNI). For our dataset, the selection of atlases with Locally Linear Embedding gives the best results. Our findings show that selection of atlases with manifold learning leads to segmentation accuracy close to or significantly higher than the state-of-the-art method and that accuracy can be increased by fine tuning the manifold learning process.

Benchmarking framework for myocardial tracking and deformation algorithms: an open access database

In this paper we present a benchmarking framework for the validation of cardiac motion analysis algorithms. The reported methods are the response to an open challenge that was issued to the medical imaging community through a MICCAI workshop. The database included magnetic resonance (MR) and 3D ultrasound (3DUS) datasets from a dynamic phantom and 15 healthy volunteers. Participants processed 3D tagged MR datasets (3DTAG), cine steady state free precession MR datasets (SSFP) and 3DUS datasets, amounting to 1158 image volumes. Ground-truth for motion tracking was based on 12 landmarks (4 walls at 3 ventricular levels). They were manually tracked by two observers in the 3DTAG data over the whole cardiac cycle, using an in-house application with 4D visualization capabilities. The median of the inter-observer variability was computed for the phantom dataset (0.77 mm) and for the volunteer datasets (0.84 mm). The ground-truth was registered to 3DUS coordinates using a point based similarity transform. Four institutions responded to the challenge by providing motion estimates for the data: Fraunhofer MEVIS (MEVIS), Bremen, Germany; Imperial College London - University College London (IUCL), UK; Universitat Pompeu Fabra (UPF), Barcelona, Spain; Inria-Asclepios project (INRIA), France. Details on the implementation and evaluation of the four methodologies are presented in this manuscript. The manually tracked landmarks were used to evaluate tracking accuracy of all methodologies. For 3DTAG, median values were computed over all time frames for the phantom dataset (MEVIS=1.20mm, IUCL=0.73 mm, UPF=1.10mm, INRIA=1.09 mm) and for the volunteer datasets (MEVIS=1.33 mm, IUCL=1.52 mm, UPF=1.09 mm, INRIA=1.32 mm). For 3DUS, median values were computed at end diastole and end systole for the phantom dataset (MEVIS=4.40 mm, UPF=3.48 mm, INRIA=4.78 mm) and for the volunteer datasets (MEVIS=3.51 mm, UPF=3.71 mm, INRIA=4.07 mm). For SSFP, median values were computed at end diastole and end systole for the phantom dataset(UPF=6.18 mm, INRIA=3.93 mm) and for the volunteer datasets (UPF=3.09 mm, INRIA=4.78 mm). Finally, strain curves were generated and qualitatively compared. Good agreement was found between the different modalities and methodologies, except for radial strain that showed a high variability in cases of lower image quality.

Toward detection and localization of instruments in minimally invasive surgery

Methods for detecting and localizing surgical instruments in laparoscopic images are an important element of advanced robotic and computer-assisted interventions. Robotic joint encoders and sensors integrated or mounted on the instrument can provide information about the tool's position, but this often has inaccuracy when transferred to the surgeon's point of view. Vision sensors are currently a promising approach for determining the position of instruments in the coordinate frame of the surgical camera. In this study, we propose a vision algorithm for localizing the instrument's pose in 3-D leaving only rotation in the axis of the tool's shaft as an ambiguity. We propose a probabilistic supervised classification method to detect pixels in laparoscopic images that belong to surgical tools. We then use the classifier output to initialize an energy minimization algorithm for estimating the pose of a prior 3-D model of the instrument within a level set framework. We show that the proposed method is robust against noise using simulated data and we perform quantitative validation of the algorithm compared to ground truth obtained using an optical tracker. Finally, we demonstrate the practical application of the technique on in vivo data from minimally invasive surgery with traditional laparoscopic and robotic instruments.

WE-E-213CD-03: Inverse-Consistent Symmetric Registration of Inner Colon Surfaces Derived from Prone and Supine CT Colonography

PURPOSE

Robust registration of prone and supine colonie surfaces acquired during CT colonography may lead to faster and more accurate detection of colorectal cancer and polyps. Any directional bias when registering one surface to the other could precipitate incorrect anatomical correspondence and engender reader error. Despite this, non-rigid registration methods are often implemented asymmetrically, which could negatively influence the registration. We aimed to reduce directional bias and so increase robustness by adapting a cylindrical registration algorithm to be both symmetric and inverse-consistent.

METHODS

The registration task can be simplified by mapping both prone and supine colonie surfaces onto regular cylinders. Spatial correspondence can then be established in cylindrical space using the original surfaces' local shape indices. We implemented a symmetric formulation of the popular non-rigid B-spline image registration method in cylindrical space. A symmetric similarity measure computes the sum of squared differences between both cylindrical representations of prone-to-supine and supine-to-prone directions simultaneously. Inverse consistency of the transformation is enforced by adding an appropriately weighted penalty term to the optimisation function.

RESULTS

We selected 8 CT colonography patient cases with marked variation in luminal distension and surface morphology. We randomly allocated 4 of these for tuning an optimal set of registration parameters and 4 for validation. The mean inverse-consistency error was reduced by 32% from 4.8mm to 3.2mm by the new symmetric formulation. The mean registration error improved from 8.2mm to 7.3mm for 330 manually chosen reference points on the 4 validation sets.

CONCLUSIONS

A symmetric formulation of prone and supine surface registration improves the quality of registration. Information from both prone-to-supine and supine-to-prone directions helps enforce convergence towards a more accurate solution due to reduced directional bias. A more robust and accurate registration will facilitate interpretation of CT colonography and has the potential to improve existing computer-aided detection methods. The authors gratefully acknowledge financial support for this work from the NIHR program: “Imaging diagnosis of colorectal cancer: Interventions for efficient and acceptable diagnosis in symptomatic and screening populationsâ€.

An anatomically driven anisotropic diffusion filtering method for 3D SPECT reconstruction

In this study, we aim to reconstruct single-photon emission computed tomography images using anatomical information from magnetic resonance imaging as a priori knowledge about the activity distribution. The trade-off between anatomical and emission data is one of the main concerns for such studies. In this work, we propose an anatomically driven anisotropic diffusion filter (ADADF) as a penalized maximum likelihood expectation maximization optimization framework. The ADADF method has improved edge-preserving denoising characteristics compared to other smoothing penalty terms based on quadratic and non-quadratic functions. The proposed method has an important ability to retain information which is absent in the anatomy. To make our approach more stable to the noise-edge classification problem, robust statistics have been employed. Comparison of the ADADF method is performed with a successful anatomically driven technique, namely, the Bowsher prior (BP). Quantitative assessment using simulated and clinical neuroreceptor volumetric data show the advantage of the ADADF over the BP. For the modelled data, the overall image resolution, the contrast, the signal-to-noise ratio and the ability to preserve important features in the data are all improved by using the proposed method. For clinical data, the contrast in the region of interest is significantly improved using the ADADF compared to the BP, while successfully eliminating noise.

Abdominal Palpation Haptic Device for Colonoscopy Simulation Using Pneumatic Control

In this paper, we describe the development of a haptic device to be used in a simulator aiming to train the skills of gastroenterology assistants in abdominal palpation during colonoscopy, as well as to train team interaction skills for the colonoscopy team. To understand the haptic feedback forces to be simulated by the haptic device, we conducted an experiment with five participants of varying BMI. The applied forces and displacements were measured and hysteresis modeling was used to characterize the experimental data. These models were used to determine the haptic feedback forces required to simulate a BMI case in response to the real-time user interactions. The pneumatic haptic device consisted of a sphygmomanometer bladder as the haptic interface and a fuzzy controller to regulate the bladder pressure. The haptic device showed good steady state and dynamic response was adequate for simulating haptic interactions. Tracking accuracy averaged 94.2 percent within 300 ms of the reference input while the user was actively applying abdominal palpation and minor repositioning.

Distinct neuropsychological profiles correspond to distribution of cortical thinning in inherited prion disease caused by insertional mutation

BACKGROUND

The human prion diseases are a group of universally fatal neurodegenerative disorders associated with the auto-catalytic misfolding of the normal cell surface prion protein (PrP). Mutations causative of inherited human prion disease (IPD) include an insertion of six additional octapeptide repeats (6-OPRI) and a missense mutation (P102L) with large families segregating for each mutation residing in southern England. Here we report for the first time the neuropsychological and clinical assessments in these two groups.

METHOD

The cognitive profiles addressing all major domains were obtained for 26 patients (18 6-OPRI, 8 P102L) and the cortical thickness determined using 1.5T MRI in a subset of 10 (six 6-OPRI, four P102L).

RESULTS

The cognitive profiles were different in patients with the two mutations in the symptomatic phase of the disease. The 6-OPRI group had lower premorbid optimal levels of functioning (assessed on the NART) than the P102L group. In the symptomatic phase of the disease the 6-OPRI patients had significantly more executive dysfunction than the P102L group and were more impaired on tests of perception and nominal functions. There was anecdotal evidence of low premorbid social performance in the 6-OPRI but not P102L patients. Cortical thinning distribution correlated with the neuropsychological profile in the 6-OPRI group principally involving the parietal, occipital and posterior frontal regions. The small number of patients in the P102L group precluded statistical comparison between the groups.

CONCLUSIONS

The 6-OPRI patients had more widespread and severe cognitive dysfunction than the P102L group and this correlated with cortical thinning distribution.

The effect of motion correction on pharmacokinetic parameter estimation in dynamic-contrast-enhanced MRI

A dynamic-contrast-enhanced magnetic resonance imaging (DCE-MRI) dataset consists of many imaging frames, often acquired both before and after contrast injection. Due to the length of time spent acquiring images, patient motion is likely and image re-alignment or registration is required before further analysis such as pharmacokinetic model fitting. Non-rigid image registration procedures may be used to correct motion artefacts; however, a careful choice of registration strategy is required to reduce misregistration artefacts associated with enhancing features. This work investigates the effect of registration on the results of model-fitting algorithms for 52 DCE-MR mammography cases for 14 patients. Results are divided into two sections: a comparison of registration strategies in which a DCE-MRI-specific algorithm is preferred in 50% of cases, followed by an investigation of parameter changes with known applied deformations, inspecting the effect of magnitude and timing of motion artefacts. Increased motion magnitude correlates with increased model-fit residual and is seen to have a strong influence on the visibility of strongly enhancing features. Motion artefacts in images close to the contrast agent arrival have a disproportionate effect on discrepancies in parameter estimation. The choice of algorithm, magnitude of motion and timing of the motion are each shown to influence estimated pharmacokinetic parameters even when motion magnitude is small.

Evaluation of registration methods on thoracic CT: the EMPIRE10 challenge

EMPIRE10 (Evaluation of Methods for Pulmonary Image REgistration 2010) is a public platform for fair and meaningful comparison of registration algorithms which are applied to a database of intrapatient thoracic CT image pairs. Evaluation of nonrigid registration techniques is a nontrivial task. This is compounded by the fact that researchers typically test only on their own data, which varies widely. For this reason, reliable assessment and comparison of different registration algorithms has been virtually impossible in the past. In this work we present the results of the launch phase of EMPIRE10, which comprised the comprehensive evaluation and comparison of 20 individual algorithms from leading academic and industrial research groups. All algorithms are applied to the same set of 30 thoracic CT pairs. Algorithm settings and parameters are chosen by researchers expert in the configuration of their own method and the evaluation is independent, using the same criteria for all participants. All results are published on the EMPIRE10 website (http://empire10.isi.uu.nl). The challenge remains ongoing and open to new participants. Full results from 24 algorithms have been published at the time of writing. This paper details the organization of the challenge, the data and evaluation methods and the outcome of the initial launch with 20 algorithms. The gain in knowledge and future work are discussed.

Consistent multi-time-point brain atrophy estimation from the boundary shift integral

Brain atrophy measurement is increasingly important in studies of neurodegenerative diseases such as Alzheimer's disease (AD), with particular relevance to trials of potential disease-modifying drugs. Automated registration-based methods such as the boundary shift integral (BSI) have been developed to provide more precise measures of change from a pair of serial MR scans. However, when a method treats one image of the pair (typically the baseline) as the reference to which the other is compared, this systematic asymmetry risks introducing bias into the measurement. Recent concern about potential biases in longitudinal studies has led to several suggestions to use symmetric image registration, though some of these methods are limited to two time-points per subject. Therapeutic trials and natural history studies increasingly involve several serial scans, it would therefore be useful to have a method that can consistently estimate brain atrophy over multiple time-points. Here, we use the log-Euclidean concept of a within-subject average to develop affine registration and differential bias correction methods suitable for any number of time-points, yielding a longitudinally consistent multi-time-point BSI technique. Baseline, 12-month and 24-month MR scans of healthy controls, subjects with mild cognitive impairment and AD patients from the Alzheimer's Disease Neuroimaging Initiative are used for testing the bias in processing scans with different amounts of atrophy. Four tests are used to assess bias in brain volume loss from BSI: (a) inverse consistency with respect to ordering of pairs of scans 12 months apart; (b) transitivity consistency over three time-points; (c) randomly ordered back-to-back scans, expected to show no consistent change over subjects; and (d) linear regression of the atrophy rates calculated from the baseline and 12-month scans and the baseline and 24-month scans, where any additive bias should be indicated by a non-zero intercept. Results indicate that the traditional BSI processing pipeline does not exhibit significant bias due to its use of windowed sinc interpolation, but with linear interpolation and asymmetric registration, bias can be pronounced. Either improved interpolation or symmetric registration alone can greatly reduce this bias, and our proposed method combining both aspects shows no significant bias in any of the four experiments.

A reduced order explicit dynamic finite element algorithm for surgical simulation

Reduced order modelling, in which a full system response is projected onto a subspace of lower dimensionality, has been used previously to accelerate finite element solution schemes by reducing the size of the involved linear systems. In the present work we take advantage of a secondary effect of such reduction for explicit analyses, namely that the stable integration time step is increased far beyond that of the full system. This phenomenon alleviates one of the principal drawbacks of explicit methods, compared with implicit schemes. We present an explicit finite element scheme in which time integration is performed in a reduced basis. Futhermore, we present a simple procedure for imposing inhomogeneous essential boundary conditions, thus overcoming one of the principal deficiencies of such approaches. The computational benefits of the procedure within a GPU-based execution framework are examined, and an assessment of the errors introduced is given. It is shown that speedups approaching an order of magnitude are feasible, without introduction of prohibitive errors, and without hardware modifications. The procedure may have applications in interactive simulation and medical image-guidance problems, in which both speed and accuracy are vital.

A comparison of voxel and surface based cortical thickness estimation methods

Cortical thickness estimation performed in-vivo via magnetic resonance imaging is an important technique for the diagnosis and understanding of the progression of neurodegenerative diseases. Currently, two different computational paradigms exist, with methods generally classified as either surface or voxel-based. This paper provides a much needed comparison of the surface-based method FreeSurfer and two voxel-based methods using clinical data. We test the effects of computing regional statistics using two different atlases and demonstrate that this makes a significant difference to the cortical thickness results. We assess reproducibility, and show that FreeSurfer has a regional standard deviation of thickness difference on same day scans that is significantly lower than either a Laplacian or Registration based method and discuss the trade off between reproducibility and segmentation accuracy caused by bending energy constraints. We demonstrate that voxel-based methods can detect similar patterns of group-wise differences as well as FreeSurfer in typical applications such as producing group-wise maps of statistically significant thickness change, but that regional statistics can vary between methods. We use a Support Vector Machine to classify patients against controls and did not find statistically significantly different results with voxel based methods compared to FreeSurfer. Finally we assessed longitudinal performance and concluded that currently FreeSurfer provides the most plausible measure of change over time, with further work required for voxel based methods.

Reduced sample sizes for atrophy outcomes in Alzheimer's disease trials: baseline adjustment

Cerebral atrophy rate is increasingly used as an outcome measure for Alzheimer's disease (AD) trials. We used the Alzheimer's disease Neuroimaging initiative (ADNI) dataset to assess if adjusting for baseline characteristics can reduce sample sizes. Controls (n = 199), patients with mild cognitive impairment (MCI) (n = 334) and AD (n = 144) had two MRI scans, 1-year apart; approximately 55% had baseline CSF tau, p-tau, and Abeta1-42. Whole brain (KN-BSI) and hippocampal (HMAPS-HBSI) atrophy rate, and ventricular expansion (VBSI) were calculated for each group; numbers required to power a placebo-controlled trial were estimated. Sample sizes per arm (80% power, 25% absolute rate reduction) for AD were (95% CI): brain atrophy = 81 (64,109), hippocampal atrophy = 88 (68,119), ventricular expansion = 118 (92,157); and for MCI: brain atrophy = 149 (122,188), hippocampal atrophy = 201 (160,262), ventricular expansion = 234 (191,295). To detect a 25% reduction relative to normal aging required increased sample sizes approximately 3-fold (AD), and approximately 5-fold (MCI). Disease severity and Abeta1-42 contributed significantly to atrophy rate variability. Adjusting for 11 predefined covariates reduced sample sizes by up to 30%. Treatment trials in AD should consider the effects of normal aging; adjusting for baseline characteristics can significantly reduce required sample sizes.