J. T. Ratnanather

Population of anatomically variable 4D XCAT adult phantoms for imaging research and optimization.

PURPOSE The authors previously developed the 4D extended cardiac-torso (XCAT) phantom for multimodality imaging research. The XCAT consisted of highly detailed whole-body models for the standard male and female adult, including the cardiac and respiratory motions. In this work, the authors extend the XCAT beyond these reference anatomies by developing a series of anatomically variable 4D XCAT adult phantoms for imaging research, the first library of 4D computational phantoms. METHODS The initial anatomy of each phantom was based on chest-abdomen-pelvis computed tomography data from normal patients obtained from the Duke University database. The major organs and structures for each phantom were segmented from the corresponding data and defined using nonuniform rational B-spline surfaces. To complete the body, the authors manually added on the head, arms, and legs using the original XCAT adult male and female anatomies. The structures were scaled to best match the age and anatomy of the patient. A multichannel large deformation diffeomorphic metric mapping algorithm was then used to calculate the transform from the template XCAT phantom (male or female) to the target patient model. The transform was applied to the template XCAT to fill in any unsegmented structures within the target phantom and to implement the 4D cardiac and respiratory models in the new anatomy. Each new phantom was refined by checking for anatomical accuracy via inspection of the models. RESULTS Using these methods, the authors created a series of computerized phantoms with thousands of anatomical structures and modeling cardiac and respiratory motions. The database consists of 58 (35 male and 23 female) anatomically variable phantoms in total. Like the original XCAT, these phantoms can be combined with existing simulation packages to simulate realistic imaging data. Each new phantom contains parameterized models for the anatomy and the cardiac and respiratory motions and can, therefore, serve as a jumping point from which to create an unlimited number of 3D and 4D variations for imaging research. CONCLUSIONS A population of phantoms that includes a range of anatomical variations representative of the public at large is needed to more closely mimic a clinical study or trial. The series of anatomically variable phantoms developed in this work provide a valuable resource for investigating 3D and 4D imaging devices and the effects of anatomy and motion in imaging. Combined with Monte Carlo simulation programs, the phantoms also provide a valuable tool to investigate patient-specific dose and image quality, and optimization for adults undergoing imaging procedures.

Dynamic programming generation of boundaries of local coordinatized submanifolds in the neocortex: application to the planum temporale.

Dynamic programming is used to define boundaries of cortical submanifolds with focus on the planum temporale (PT) of the superior temporal gyrus (STG), which has been implicated in a variety of neuropsychiatric disorders. To this end, automated methods are used to generate the PT manifold from 10 high-resolution MRI subvolumes ROI masks encompassing the STG. A procedure to define the subvolume ROI masks from original MRI brain scans is developed. Bayesian segmentation is then used to segment the subvolumes into cerebrospinal fluid, gray matter (GM), and white matter (WM). 3D isocontouring using the intensity value at which there is equal probability of GM and WM is used to reconstruct the triangulated graph representing the STG cortical surface, enabling principal curvature at each point on the graph to be computed. Dynamic programming is used to delineate the PT manifold by tracking principal curves from the retro-insular end of the Heschls gyrus (HG) to the STG, along the posterior STG up to the start of the ramus and back to the retro-insular end of the HG. A coordinate system is then defined on the PT manifold. The origin is defined by the retro-insular end of the HG and the y-axis passes through the point on the posterior STG where the ramus begins. Automated labeling of GM in the STG is robust with L(1) distances between Bayesian and manual segmentation in the range 0.001-0.12 (n = 20). PT reconstruction is also robust with 90% of the vertices of the reconstructed PT within about 1 voxel (n = 20) from semiautomated contours. Finally, the reliability index (based on interrater intraclass correlation) for the surface area derived from repeated reconstructions is 0.96 for the left PT and 0.94 for the right PT, thus demonstrating the robustness of dynamic programming in defining a coordinate system on the PT. It provides a method with potential significance in the study of neuropsychiatric disorders.

Labeled cortical mantle distance maps of the cingulate quantify differences between dementia of the Alzheimer type and healthy aging

The cingulate gyri in 37 subjects with and without early dementia of the Alzheimer type (DAT) were studied by using MRI at 1.0 mm3 isotropic resolution. Groups were segregated into young controls (n = 10), age-matched normal controls (n = 10), very mild DAT (n = 8), and mild DAT (n = 9). By using automated Bayesian segmentation of the cortex and gray matter/white matter (GM/WM) isosurface generation, tissue compartments were labeled into gray, white, and cerebrospinal fluid as a function of distance from the GM/WM isosurface. Cortical mantle distance maps are generated profiling the GM volume and cortical mantle distribution as a function of distance from the cortical surface. Probabilistic tests based on generalizations of Wilcoxon–Mann–Whitney tests were applied to quantify cortical mantle distribution changes with normal and abnormal aging. We find no significant change between young controls and healthy aging as measured by the GM volume and cortical mantle distribution as a function of distance in both anterior and posterior regions of the cingulate. Significant progression of GM loss is seen in the very mild DAT and mild DAT groups in all areas of the cingulate. Posterior regions show both GM volume loss as well as significant cortical mantle distribution decrease with the onset of mild DAT. The “shape of the cortical mantle” as measured by the cortical mantle distance profiles manifests a pronounced increase in variability with mild DAT.

Validating Cortical Surface Analysis of Medial Prefrontal Cortex

This paper describes cortical analysis of 19 high resolution MRI subvolumes of medial prefrontal cortex (MPFC), a region that has been implicated in major depressive disorder. An automated Bayesian segmentation is used to delineate the MRI subvolumes into cerebrospinal fluid (CSF), gray matter (GM), white matter (WM), and partial volumes of either CSF/GM or GM/WM. The intensity value at which there is equal probability of GM and GM/WM partial volume is used to reconstruct MPFC cortical surfaces based on a 3-D isocontouring algorithm. The segmented data and the generated surfaces are validated by comparison with hand segmented data and semiautomated contours, respectively. The L(1) distances between Bayesian and hand segmented data are 0.05-0.10 (n = 5). Fifty percent of the voxels of the reconstructed surface lie within 0.12-0.28 mm (n = 14) from the semiautomated contours. Cortical thickness metrics are generated in the form of frequency of occurrence histograms for GM and WM labelled voxels as a function of their position from the cortical surface. An algorithm to compute the surface area of the GM/WM interface of the MPFC subvolume is described. These methods represent a novel approach to morphometric chacterization of regional cortex features which may be important in the study of psychiatric disorders such as major depression.

A set of 4D pediatric XCAT reference phantoms for multimodality research.

PURPOSE The authors previously developed an adult population of 4D extended cardiac-torso (XCAT) phantoms for multimodality imaging research. In this work, the authors develop a reference set of 4D pediatric XCAT phantoms consisting of male and female anatomies at ages of newborn, 1, 5, 10, and 15 years. These models will serve as the foundation from which the authors will create a vast population of pediatric phantoms for optimizing pediatric CT imaging protocols. METHODS Each phantom was based on a unique set of CT data from a normal patient obtained from the Duke University database. The datasets were selected to best match the reference values for height and weight for the different ages and genders according to ICRP Publication 89. The major organs and structures were segmented from the CT data and used to create an initial pediatric model defined using nonuniform rational B-spline surfaces. The CT data covered the entire torso and part of the head. To complete the body, the authors manually added on the top of the head and the arms and legs using scaled versions of the XCAT adult models or additional models created from cadaver data. A multichannel large deformation diffeomorphic metric mapping algorithm was then used to calculate the transform from a template XCAT phantom (male or female 50th percentile adult) to the target pediatric model. The transform was applied to the template XCAT to fill in any unsegmented structures within the target phantom and to implement the 4D cardiac and respiratory models in the new anatomy. The masses of the organs in each phantom were matched to the reference values given in ICRP Publication 89. The new reference models were checked for anatomical accuracy via visual inspection. RESULTS The authors created a set of ten pediatric reference phantoms that have the same level of detail and functionality as the original XCAT phantom adults. Each consists of thousands of anatomical structures and includes parameterized models for the cardiac and respiratory motions. Based on patient data, the phantoms capture the anatomic variations of childhood, such as the development of bone in the skull, pelvis, and long bones, and the growth of the vertebrae and organs. The phantoms can be combined with existing simulation packages to generate realistic pediatric imaging data from different modalities. CONCLUSIONS The development of patient-derived pediatric computational phantoms is useful in providing variable anatomies for simulation. Future work will expand this ten-phantom base to a host of pediatric phantoms representative of the public at large. This can provide a means to evaluate and improve pediatric imaging devices and to optimize CT protocols in terms of image quality and radiation dose.

Patient Specific Computerized Phantoms to Estimate Dose in Pediatric CT

We create a series of detailed computerized phantoms to estimate patient organ and effective dose in pediatric CT and investigate techniques for efficiently creating patient-specific phantoms based on imaging data. The initial anatomy of each phantom was previously developed based on manual segmentation of pediatric CT data. Each phantom was extended to include a more detailed anatomy based on morphing an existing adult phantom in our laboratory to match the framework (based on segmentation) defined for the target pediatric model. By morphing a template anatomy to match the patient data in the LDDMM framework, it was possible to create a patient specific phantom with many anatomical structures, some not visible in the CT data. The adult models contain thousands of defined structures that were transformed to define them in each pediatric anatomy. The accuracy of this method, under different conditions, was tested using a known voxelized phantom as the target. Errors were measured in terms of a distance map between the predicted organ surfaces and the known ones. We also compared calculated dose measurements to see the effect of different magnitudes of errors in morphing. Despite some variations in organ geometry, dose measurements from morphing predictions were found to agree with those calculated from the voxelized phantom thus demonstrating the feasibility of our methods.

Statistical Analysis of Cortical Morphometrics Using Pooled Distances Based on Labeled Cortical Distance Maps

Neuropsychiatric disorders have been demonstrated to manifest shape differences in cortical structures. Labeled Cortical Distance Mapping (LCDM) is a powerful tool in quantifying such morphometric differences and characterizes the morphometry of the laminar cortical mantle of cortical structures. Specifically, LCDM data are distances of labeled gray matter (GM) voxels with respect to the gray/white matter cortical surface. Volumes and descriptive measures (such as means and variances for each subject) based on LCDM distances provide descriptive summary information on some of the shape characteristics. However, additional morphometrics are contained in the data and their analysis may provide additional clues to underlying differences in cortical characteristics. To use more of this information, we pool (merge) LCDM distances from subjects in the same group. These pooled distances can help detect morphometric differences between groups, but do not provide information about the locations of such differences in the tissue in question. In this article, we check for the influence of the assumption violations on the analysis of pooled LCDM distances. We demonstrate that the classical parametric tests are robust to the non-normality and within sample dependence of LCDM distances and nonparametric tests are robust to within sample dependence of LCDM distances. We specify the types of alternatives for which the tests are more sensitive. We also show that the pooled LCDM distances provide powerful results for group differences in distribution of LCDM distances. As an illustrative example, we use GM in the ventral medial prefrontal cortex (VMPFC) in subjects with major depressive disorder (MDD), subjects at high risk (HR) of MDD, and healthy subjects. Significant morphometric differences were found in VMPFC due to MDD or being at HR. In particular, the analysis indicated that distances in left and right VMPFCs tend to decrease due to MDD or being at HR, possibly as a result of thinning. The methodology can also be applied to other cortical structures.

Ventromedial Prefrontal Cortex Thinning in Preschool-Onset Depression

BACKGROUND The ventromedial prefrontal cortex (VMPFC) is a key center of affect regulation and processing, fundamental aspects of emotional competence which are disrupted in mood disorders. Structural alterations of VMPFC have consistently been observed in adult major depression and are associated with depression severity, yet it is unknown whether young children with depression demonstrate similar abnormalities. We investigated cortical thickness differences in the VMPFC of children with a history of preschool-onset depression (PO-MDD). METHODS Participants in a longitudinal study of PO-MDD underwent structural brain imaging between the ages of 7 and 12 years. Using local cortical distance metrics, cortical thickness of the VMPFC was compared in children with and without a history of PO-MDD. RESULTS Children previously diagnosed with PO-MDD (n=34) had significantly thinner right VMPFC vs. children without a history of PO-MDD [(n=95); F(1,126)=5.97, (p=.016)]. This effect was specific to children with a history of PO-MDD vs. other psychiatric conditions and was independent of comorbid anxiety or externalizing disorders. Decreases in right VMPFC thickness were predicted by preschool depressive symptoms independent of depressive symptoms in school age. LIMITATIONS Results are cross-sectional and cannot distinguish whether thinner right VMPFC represents a vulnerability marker of MDD, consequence of MDD, or marker of remitted MDD. Longitudinal imaging is needed to contextualize how this difference relates to normative VMPFC structural development. CONCLUSIONS Onset of depression at preschool age was associated with decreased cortical thickness of right VMPFC. This finding implicates the VMPFC in depression from very early stages of brain development.

The development of a population of 4D pediatric XCAT phantoms for imaging research and optimization

PURPOSE We previously developed a set of highly detailed 4D reference pediatric extended cardiac-torso (XCAT) phantoms at ages of newborn, 1, 5, 10, and 15 yr with organ and tissue masses matched to ICRP Publication 89 values. In this work, we extended this reference set to a series of 64 pediatric phantoms of varying age and height and body mass percentiles representative of the public at large. The models will provide a library of pediatric phantoms for optimizing pediatric imaging protocols. METHODS High resolution positron emission tomography-computed tomography data obtained from the Duke University database were reviewed by a practicing experienced radiologist for anatomic regularity. The CT portion of the data was then segmented with manual and semiautomatic methods to form a target model defined using nonuniform rational B-spline surfaces. A multichannel large deformation diffeomorphic metric mapping algorithm was used to calculate the transform from the best age matching pediatric XCAT reference phantom to the patient target. The transform was used to complete the target, filling in the nonsegmented structures and defining models for the cardiac and respiratory motions. The complete phantoms, consisting of thousands of structures, were then manually inspected for anatomical accuracy. The mass for each major tissue was calculated and compared to linearly interpolated ICRP values for different ages. RESULTS Sixty four new pediatric phantoms were created in this manner. Each model contains the same level of detail as the original XCAT reference phantoms and also includes parameterized models for the cardiac and respiratory motions. For the phantoms that were 10 yr old and younger, we included both sets of reproductive organs. This gave them the capability to simulate both male and female anatomy. With this, the population can be expanded to 92. Wide anatomical variation was clearly seen amongst the phantom models, both in organ shape and size, even for models of the same age and sex. The phantoms can be combined with existing simulation packages to generate realistic pediatric imaging data from different modalities. CONCLUSIONS This work provides a large cohort of highly detailed pediatric phantoms with 4D capabilities of varying age, height, and body mass. The population of phantoms will provide a vital tool with which to optimize 3D and 4D pediatric imaging devices and techniques in terms of image quality and radiation-absorbed dose.

Generating patient-specific dosimetry phantoms with whole-body diffeomorphic image registration

This work describes a protocol for creating a population of pediatric, patient specific, computational dosimetry phantoms. Pediatric CT data are mapped to a detailed adult template using multi channel (MC) - large deformation diffeomorphic metric mapping (LDDMM) applied to whole body images. Challenges are studied and overcome using 2D LDDMM applied to simulated phantoms. The protocol involves automatic placement of landmarks, landmark based affine and LDDMM registration, followed by 4 applications of MC-LDDMM, with successively decreasing smoothness constraints. Deformed adult organ surfaces typically agree with pediatric organ surfaces within 1–2 voxels (∼2.3 mm).