Challenges in measuring angles between craniofacial structures.

OBJECTIVE: Three-dimensional (3D) angular measurements between craniofacial planes pose challenges to quantify maxillary and mandibular skeletal discrepancies in surgical treatment planning. This study aims to compare the reproducibility and reliability of two modules to measure angles between planes or lines in 3D virtual surface models. METHODOLOGY: Twenty oriented 3D virtual surface models de-identified and constructed from CBCT scans were randomly selected. Three observers placed landmarks and oriented planes to determine angular measurements of pitch, roll and yaw using (1) 3D pre-existing planes, (2) 3D planes created from landmarks and (3) lines created from landmarks. Inter- and intra-observer reproducibility and repeatability were examined using the Intra-Class Correlation (ICC) test. One observer repeated the measurements with an interval of 15 days. ANOVA was applied to compare the 3 methods. RESULTS: The three methods tested provided statistically similar, reproducible and reliable angular measurements of the facial structures. A strong ICC varying from 0.92 to 1.00 was found for the intra-observer agreement. The inter-observer ICC varied from 0.84 to 1.00. CONCLUSION: Measurements of 3D angles between facial planes in a common coordinate system are reproducible and repeatable either using 3D pre-existing planes, created based on landmarks or angles between lines created from landmarks.

3D superimposition of craniofacial imaging-The utility of multicentre collaborations.

Clinical applications of 3D image registration and superimposition have contributed to better understanding growth changes and clinical outcomes. The use of 3D dental and craniofacial imaging in dentistry requires validate image analysis methods for improved diagnosis, treatment planning, navigation and assessment of treatment response. Volumetric 3D images, such as cone-beam computed tomography, can now be superimposed by voxels, surfaces or landmarks. Regardless of the image modality or the software tools, the concepts of regions or points of reference affect all quantitative of qualitative assessments. This study reviews current state of the art in 3D image analysis including 3D superimpositions relative to the cranial base and different regional superimpositions, the development of open source and commercial tools for 3D analysis, how this technology has increased clinical research collaborations from centres all around the globe, some insight on how to incorporate artificial intelligence for big data analysis and progress towards personalized orthodontics.

Challenges in measuring angles between craniofacial structures

Abstract Objective: Three-dimensional (3D) angular measurements between craniofacial planes pose challenges to quantify maxillary and mandibular skeletal discrepancies in surgical treatment planning. This study aims to compare the reproducibility and reliability of two modules to measure angles between planes or lines in 3D virtual surface models. Methodology: Twenty oriented 3D virtual surface models de-identified and constructed from CBCT scans were randomly selected. Three observers placed landmarks and oriented planes to determine angular measurements of pitch, roll and yaw using (1) 3D pre-existing planes, (2) 3D planes created from landmarks and (3) lines created from landmarks. Inter- and intra-observer reproducibility and repeatability were examined using the Intra-Class Correlation (ICC) test. One observer repeated the measurements with an interval of 15 days. ANOVA was applied to compare the 3 methods. Results: The three methods tested provided statistically similar, reproducible and reliable angular measurements of the facial structures. A strong ICC varying from 0.92 to 1.00 was found for the intra-observer agreement. The inter-observer ICC varied from 0.84 to 1.00. Conclusion: Measurements of 3D angles between facial planes in a common coordinate system are reproducible and repeatable either using 3D pre-existing planes, created based on landmarks or angles between lines created from landmarks.

Automatic quantification framework to detect cracks in teeth.

Studies show that cracked teeth are the third most common cause for tooth loss in industrialized countries. If detected early and accurately, patients can retain their teeth for a longer time. Most cracks are not detected early because of the discontinuous symptoms and lack of good diagnostic tools. Currently used imaging modalities like Cone Beam Computed Tomography (CBCT) and intraoral radiography often have low sensitivity and do not show cracks clearly. This paper introduces a novel method that can detect, quantify, and localize cracks automatically in high resolution CBCT (hr-CBCT) scans of teeth using steerable wavelets and learning methods. These initial results were created using hr-CBCT scans of a set of healthy teeth and of teeth with simulated longitudinal cracks. The cracks were simulated using multiple orientations. The crack detection was trained on the most significant wavelet coefficients at each scale using a bagged classifier of Support Vector Machines. Our results show high discriminative specificity and sensitivity of this method. The framework aims to be automatic, reproducible, and open-source. Future work will focus on the clinical validation of the proposed techniques on different types of cracks ex-vivo. We believe that this work will ultimately lead to improved tracking and detection of cracks allowing for longer lasting healthy teeth.

A web-based system for neural network based classification in temporomandibular joint osteoarthritis.

OBJECTIVE: The purpose of this study is to describe the methodological innovations of a web-based system for storage, integration and computation of biomedical data, using a training imaging dataset to remotely compute a deep neural network classifier of temporomandibular joint osteoarthritis (TMJOA). METHODS: This study imaging dataset consisted of three-dimensional (3D) surface meshes of mandibular condyles constructed from cone beam computed tomography (CBCT) scans. The training dataset consisted of 259 condyles, 105 from control subjects and 154 from patients with diagnosis of TMJ OA. For the image analysis classification, 34 right and left condyles from 17 patients (39.9 ±â€¯11.7 years), who experienced signs and symptoms of the disease for less than 5 years, were included as the testing dataset. For the integrative statistical model of clinical, biological and imaging markers, the sample consisted of the same 17 test OA subjects and 17 age and sex matched control subjects (39.4 ±â€¯15.4 years), who did not show any sign or symptom of OA. For these 34 subjects, a standardized clinical questionnaire, blood and saliva samples were also collected. The technological methodologies in this study include a deep neural network classifier of 3D condylar morphology (ShapeVariationAnalyzer, SVA), and a flexible web-based system for data storage, computation and integration (DSCI) of high dimensional imaging, clinical, and biological data. RESULTS: The DSCI system trained and tested the neural network, indicating 5 stages of structural degenerative changes in condylar morphology in the TMJ with 91% close agreement between the clinician consensus and the SVA classifier. The DSCI remotely ran with a novel application of a statistical analysis, the Multivariate Functional Shape Data Analysis, that computed high dimensional correlations between shape 3D coordinates, clinical pain levels and levels of biological markers, and then graphically displayed the computation results. CONCLUSIONS: The findings of this study demonstrate a comprehensive phenotypic characterization of TMJ health and disease at clinical, imaging and biological levels, using novel flexible and versatile open-source tools for a web-based system that provides advanced shape statistical analysis and a neural network based classification of temporomandibular joint osteoarthritis.

Diagnostic Index: An open-source tool to classify TMJ OA condyles.

Osteoarthritis (OA) of temporomandibular joints (TMJ) occurs in about 40% of the patients who present TMJ disorders. Despite its prevalence, OA diagnosis and treatment remain controversial since there are no clear symptoms of the disease, especially in early stages. Quantitative tools based on 3D imaging of the TMJ condyle have the potential to help characterize TMJ OA changes. The goals of the tools proposed in this study are to ultimately develop robust imaging markers for diagnosis and assessment of treatment efficacy. This work proposes to identify differences among asymptomatic controls and different clinical phenotypes of TMJ OA by means of Statistical Shape Modeling (SSM), obtained via clinical expert consensus. From three different grouping schemes (with 3, 5 and 7 groups), our best results reveal that that the majority (74.5%) of the classifications occur in agreement with the groups assigned by consensus between our clinical experts. Our findings suggest the existence of different disease-based phenotypic morphologies in TMJ OA. Our preliminary findings with statistical shape modeling based biomarkers may provide a quantitative staging of the disease. The methodology used in this study is included in an open source image analysis toolbox, to ensure reproducibility and appropriate distribution and dissemination of the solution proposed.

The UNC-Wisconsin Rhesus Macaque Neurodevelopment Database: A Structural MRI and DTI Database of Early Postnatal Development.

Rhesus macaques are commonly used as a translational animal model in neuroimaging and neurodevelopmental research. In this report, we present longitudinal data from both structural and diffusion MRI images generated on a cohort of 34 typically developing monkeys from 2 weeks to 36 months of age. All images have been manually skull stripped and are being made freely available via an online repository for use by the research community.

White Matter Fiber-based Analysis of T1w/T2w Ratio Map.

PURPOSE: To develop, test, evaluate and apply a novel tool for the white matter fiber-based analysis of T1w/T2w ratio maps quantifying myelin content. BACKGROUND: The cerebral white matter in the human brain develops from a mostly non-myelinated state to a nearly fully mature white matter myelination within the first few years of life. High resolution T1w/T2w ratio maps are believed to be effective in quantitatively estimating myelin content on a voxel-wise basis. We propose the use of a fiber-tract-based analysis of such T1w/T2w ratio data, as it allows us to separate fiber bundles that a common regional analysis imprecisely groups together, and to associate effects to specific tracts rather than large, broad regions. METHODS: We developed an intuitive, open source tool to facilitate such fiber-based studies of T1w/T2w ratio maps. Via its Graphical User Interface (GUI) the tool is accessible to non-technical users. The framework uses calibrated T1w/T2w ratio maps and a prior fiber atlas as an input to generate profiles of T1w/T2w values. The resulting fiber profiles are used in a statistical analysis that performs along-tract functional statistical analysis. We applied this approach to a preliminary study of early brain development in neonates. RESULTS: We developed an open-source tool for the fiber based analysis of T1w/T2w ratio maps and tested it in a study of brain development.

Autotract: Automatic cleaning and tracking of fibers.

We propose a new tool named Autotract to automate fiber tracking in diffusion tensor imaging (DTI). Autotract uses prior knowledge from a source DTI and a set of corresponding fiber bundles to extract new fibers for a target DTI. Autotract starts by aligning both DTIs and uses the source fibers as seed points to initialize a tractography algorithm. We enforce similarity between the propagated source fibers and automatically traced fibers by computing metrics such as fiber length and fiber distance between the bundles. By analyzing these metrics, individual fiber tracts can be pruned. As a result, we show that both bundles have similar characteristics. Additionally, we compare the automatically traced fibers against bundles previously generated and validated in the target DTI by an expert. This work is motivated by medical applications in which known bundles of fiber tracts in the human brain need to be analyzed for multiple datasets.

UNC-Emory Infant Atlases for Macaque Brain Image Analysis: Postnatal Brain Development through 12 Months.

Computational anatomical atlases have shown to be of immense value in neuroimaging as they provide age appropriate reference spaces alongside ancillary anatomical information for automated analysis such as subcortical structural definitions, cortical parcellations or white fiber tract regions. Standard workflows in neuroimaging necessitate such atlases to be appropriately selected for the subject population of interest. This is especially of importance in early postnatal brain development, where rapid changes in brain shape and appearance render neuroimaging workflows sensitive to the appropriate atlas choice. We present here a set of novel computation atlases for structural MRI and Diffusion Tensor Imaging as crucial resource for the analysis of MRI data from non-human primate rhesus monkey (Macaca mulatta) data in early postnatal brain development. Forty socially-housed infant macaques were scanned longitudinally at ages 2 weeks, 3, 6, and 12 months in order to create cross-sectional structural and DTI atlases via unbiased atlas building at each of these ages. Probabilistic spatial prior definitions for the major tissue classes were trained on each atlas with expert manual segmentations. In this article we present the development and use of these atlases with publicly available tools, as well as the atlases themselves, which are publicly disseminated to the scientific community.

Early Postnatal Myelin Content Estimate of White Matter via T1w/T2w Ratio.

PURPOSE: To develop and evaluate a novel processing framework for the relative quantification of myelin content in cerebral white matter (WM) regions from brain MRI data via a computed ratio of T1 to T2 weighted intensity values. DATA: We employed high resolution (1mm3 isotropic) T1 and T2 weighted MRI from 46 (28 male, 18 female) neonate subjects (typically developing controls) scanned on a Siemens Tim Trio 3T at UC Irvine. METHODS: We developed a novel, yet relatively straightforward image processing framework for WM myelin content estimation based on earlier work by Glasser et al. We first co-register the structural MRI data to correct for motion. Then, background areas are masked out via a joint T1w and T2 foreground mask computed. Raw T1w/T2w-ratios images are computed next. For purpose of calibration across subjects, we first coarsely segment the fat-rich facial regions via an atlas co-registration. Linear intensity rescaling based on median T1w/T2w-ratio values in those facial regions yields calibrated T1w/T2w-ratio images. Mean values in lobar regions are evaluated using standard statistical analysis to investigate their interaction with age at scan. RESULTS: Several lobes have strongly positive significant interactions of age at scan with the computed T1w/T2w-ratio. Most regions do not show sex effects. A few regions show no measurable effects of change in myelin content change within the first few weeks of postnatal development, such as cingulate and CC areas, which we attribute to sample size and measurement variability. CONCLUSIONS: We developed and evaluated a novel way to estimate white matter myelin content for use in studies of brain white matter development.

Validation of Osteoarthritis synthetic defect database via non-rigid registration.

Temporomandibular joint (TMJ) disorders are a group of conditions that cause pain and dysfunction in the jaw joint and the muscles controlling jaw movement. However, diagnosis and treatment of these conditions remain controversial. To date, there is no single sign, symptom, or test that can clearly diagnose early stages of osteoarthritis (OA). Instead, the diagnosis is based on a consideration of several factors, including radiological evaluation. The current radiological diagnosis scores of TMJ pathology are subject to misdiagnosis. We believe these scores are limited by the acquisition procedures, such as oblique cuts of the CT and head positioning errors, and can lead to incorrect diagnoses of flattening of the head of the condyle, formation of osteophytes, or condylar pitting. This study consists of creating and validating a methodological framework to simulate defects in CBCT scans of known location and size, in order to create synthetic TMJ OA database. User-generated defects were created using a non-rigid deformation protocol in CBCT. All segmentation evaluation, surface distances and linear distances from the user-generated to the simulated defects showed our methodological framework to be very precise and within a voxel (0.5 mm) of magnitude. A TMJ OA synthetic database will be created next, and evaluated by expert radiologists, and this will serve to evaluate how sensitive the current radiological diagnosis tools are.

Automatic Tissue Segmentation of Neonate Brain MR Images with Subject-specific Atlases.

Automatic tissue segmentation of the neonate brain using Magnetic Resonance Images (MRI) is extremely important to study brain development and perform early diagnostics but is challenging due to high variability and inhomogeneity in contrast throughout the image due to incomplete myelination of the white matter tracts. For these reasons, current methods often totally fail or give unsatisfying results. Furthermore, most of the subcortical midbrain structures are misclassified due to a lack of contrast in these regions. We have developed a novel method that creates a probabilistic subject-specific atlas based on a population atlas currently containing a number of manually segmented cases. The generated subject-specific atlas is sharp and adapted to the subject that is being processed. We then segment brain tissue classes using the newly created atlas with a single-atlas expectation maximization based method. Our proposed method leads to a much lower failure rate in our experiments. The overall segmentation results are considerably improved when compared to using a non-subject-specific, population average atlas. Additionally, we have incorporated diffusion information obtained from Diffusion Tensor Images (DTI) to improve the detection of white matter that is not visible at this early age in structural MRI (sMRI) due to a lack of myelination. Although this necessitates the acquisition of an additional sequence, the diffusion information improves the white matter segmentation throughout the brain, especially for the mid-brain structures such as the corpus callosum and the internal capsule.

Incorporating 3-dimensional models in online articles.

INTRODUCTION: The aims of this article are to introduce the capability to view and interact with 3-dimensional (3D) surface models in online publications, and to describe how to prepare surface models for such online 3D visualizations. METHODS: Three-dimensional image analysis methods include image acquisition, construction of surface models, registration in a common coordinate system, visualization of overlays, and quantification of changes. Cone-beam computed tomography scans were acquired as volumetric images that can be visualized as 3D projected images or used to construct polygonal meshes or surfaces of specific anatomic structures of interest. The anatomic structures of interest in the scans can be labeled with color (3D volumetric label maps), and then the scans are registered in a common coordinate system using a target region as the reference. The registered 3D volumetric label maps can be saved in .obj, .ply, .stl, or .vtk file formats and used for overlays, quantification of differences in each of the 3 planes of space, or color-coded graphic displays of 3D surface distances. RESULTS: All registered 3D surface models in this study were saved in .vtk file format and loaded in the Elsevier 3D viewer. In this study, we describe possible ways to visualize the surface models constructed from cone-beam computed tomography images using 2D and 3D figures. The 3D surface models are available in the article's online version for viewing and downloading using the reader's software of choice. These 3D graphic displays are represented in the print version as 2D snapshots. Overlays and color-coded distance maps can be displayed using the reader's software of choice, allowing graphic assessment of the location and direction of changes or morphologic differences relative to the structure of reference. The interpretation of 3D overlays and quantitative color-coded maps requires basic knowledge of 3D image analysis. CONCLUSIONS: When submitting manuscripts, authors can now upload 3D models that will allow readers to interact with or download them. Such interaction with 3D models in online articles now will give readers and authors better understanding and visualization of the results.

Characterization of subtle brain abnormalities in a mouse model of Hedgehog pathway antagonist-induced cleft lip and palate.

Subtle behavioral and cognitive deficits have been documented in patient cohorts with orofacial clefts (OFCs). Recent neuroimaging studies argue that these traits are associated with structural brain abnormalities but have been limited to adolescent and adult populations where brain plasticity during infancy and childhood may be a confounding factor. Here, we employed high resolution magnetic resonance microscopy to examine primary brain morphology in a mouse model of OFCs. Transient in utero exposure to the Hedgehog (Hh) signaling pathway antagonist cyclopamine resulted in a spectrum of facial dysmorphology, including unilateral and bilateral cleft lip and palate, cleft of the secondary palate only, and a non-cleft phenotype marked by midfacial hypoplasia. Relative to controls, cyclopamine-exposed fetuses exhibited volumetric differences in several brain regions, including hypoplasia of the pituitary gland and olfactory bulbs, hyperplasia of the forebrain septal region, and expansion of the third ventricle. However, in affected fetuses the corpus callosum was intact and normal division of the forebrain was observed. This argues that temporally-specific Hh signaling perturbation can result in typical appearing OFCs in the absence of holoprosencephaly--a condition classically associated with Hh pathway inhibition and frequently co-occurring with OFCs. Supporting the premise that some forms of OFCs co-occur with subtle brain malformations, these results provide a possible ontological basis for traits identified in clinical populations. They also argue in favor of future investigations into genetic and/or environmental modulation of the Hh pathway in the etiopathogenesis of orofacial clefting.

DTIPrep: quality control of diffusion-weighted images.

In the last decade, diffusion MRI (dMRI) studies of the human and animal brain have been used to investigate a multitude of pathologies and drug-related effects in neuroscience research. Study after study identifies white matter (WM) degeneration as a crucial biomarker for all these diseases. The tool of choice for studying WM is dMRI. However, dMRI has inherently low signal-to-noise ratio and its acquisition requires a relatively long scan time; in fact, the high loads required occasionally stress scanner hardware past the point of physical failure. As a result, many types of artifacts implicate the quality of diffusion imagery. Using these complex scans containing artifacts without quality control (QC) can result in considerable error and bias in the subsequent analysis, negatively affecting the results of research studies using them. However, dMRI QC remains an under-recognized issue in the dMRI community as there are no user-friendly tools commonly available to comprehensively address the issue of dMRI QC. As a result, current dMRI studies often perform a poor job at dMRI QC. Thorough QC of dMRI will reduce measurement noise and improve reproducibility, and sensitivity in neuroimaging studies; this will allow researchers to more fully exploit the power of the dMRI technique and will ultimately advance neuroscience. Therefore, in this manuscript, we present our open-source software, DTIPrep, as a unified, user friendly platform for thorough QC of dMRI data. These include artifacts caused by eddy-currents, head motion, bed vibration and pulsation, venetian blind artifacts, as well as slice-wise and gradient-wise intensity inconsistencies. This paper summarizes a basic set of features of DTIPrep described earlier and focuses on newly added capabilities related to directional artifacts and bias analysis.

UNC-Utah NA-MIC framework for DTI fiber tract analysis.

Diffusion tensor imaging has become an important modality in the field of neuroimaging to capture changes in micro-organization and to assess white matter integrity or development. While there exists a number of tractography toolsets, these usually lack tools for preprocessing or to analyze diffusion properties along the fiber tracts. Currently, the field is in critical need of a coherent end-to-end toolset for performing an along-fiber tract analysis, accessible to non-technical neuroimaging researchers. The UNC-Utah NA-MIC DTI framework represents a coherent, open source, end-to-end toolset for atlas fiber tract based DTI analysis encompassing DICOM data conversion, quality control, atlas building, fiber tractography, fiber parameterization, and statistical analysis of diffusion properties. Most steps utilize graphical user interfaces (GUI) to simplify interaction and provide an extensive DTI analysis framework for non-technical researchers/investigators. We illustrate the use of our framework on a small sample, cross sectional neuroimaging study of eight healthy 1-year-old children from the Infant Brain Imaging Study (IBIS) Network. In this limited test study, we illustrate the power of our method by quantifying the diffusion properties at 1 year of age on the genu and splenium fiber tracts.

Software-based Diffusion MR Human Brain Phantom for Evaluating Fiber-tracking Algorithms.

Fiber tracking provides insights into the brain white matter network and has become more and more popular in diffusion MR imaging. Hardware or software phantom provides an essential platform to investigate, validate and compare various tractography algorithms towards a "gold standard". Software phantoms excel due to their flexibility in varying imaging parameters, such as tissue composition, SNR, as well as potential to model various anatomies and pathologies. This paper describes a novel method in generating diffusion MR images with various imaging parameters from realistically appearing, individually varying brain anatomy based on predefined fiber tracts within a high-resolution human brain atlas. Specifically, joint, high resolution DWI and structural MRI brain atlases were constructed with images acquired from 6 healthy subjects (age 22-26) for the DWI data and 56 healthy subject (age 18-59) for the structural MRI data. Full brain fiber tracking was performed with filtered, two-tensor tractography in atlas space. A deformation field based principal component model from the structural MRI as well as unbiased atlas building was then employed to generate synthetic structural brain MR images that are individually varying. Atlas fiber tracts were accordingly warped into each synthetic brain anatomy. Diffusion MR images were finally computed from these warped tracts via a composite hindered and restricted model of diffusion with various imaging parameters for gradient directions, image resolution and SNR. Furthermore, an open-source program was developed to evaluate the fiber tracking results both qualitatively and quantitatively based on various similarity measures.

Fully automated rodent brain MR image processing pipeline on a Midas server: from acquired images to region-based statistics.

Magnetic resonance imaging (MRI) of rodent brains enables study of the development and the integrity of the brain under certain conditions (alcohol, drugs etc.). However, these images are difficult to analyze for biomedical researchers with limited image processing experience. In this paper we present an image processing pipeline running on a Midas server, a web-based data storage system. It is composed of the following steps: rigid registration, skull-stripping, average computation, average parcellation, parcellation propagation to individual subjects, and computation of region-based statistics on each image. The pipeline is easy to configure and requires very little image processing knowledge. We present results obtained by processing a data set using this pipeline and demonstrate how this pipeline can be used to find differences between populations.

Comparison of magnetic resonance imaging in live vs. post mortem rat brains.

Magnetic Resonance Imaging (MRI) is an increasingly popular technique for examining neurobiology in rodents because it is both noninvasive and nondestructive. MRI scans can be acquired from either live or post mortem specimens. In vivo scans have a key advantage in that subjects can be scanned at multiple time-points in longitudinal studies. However, repeated exposure to anesthesia and stress may confound studies. In contrast, post mortem scans offer improved image quality and increased signal-to-noise ratio (SNR) due to several key advantages: First, the images are not disrupted by motion and pulsation artifacts. Second, they allow the brain tissue to be perfused with contrast agents, enhancing tissue contrast. Third, they allow longer image acquisition times, yielding higher resolution and/or improved SNR. Fourth, they allow assessment of groups of animals at the same age without scheduling complications. Despite these advantages, researchers are often skeptical of post mortem MRI scans because of uncertainty about whether the fixation process alters the MRI measurements. To address these concerns, we present a thorough comparative study of in vivo and post mortem MRI scans in healthy male Wistar rats at three age points throughout adolescence (postnatal days 28 through 80). For each subject, an in vivo scan was acquired, followed by perfusion and two post mortem scans at two different MRI facilities. The goal was to assess robustness of measurements, to detect any changes in volumetric measurements after fixation, and to investigate any differential bias that may exist between image acquisition techniques. We present this volumetric analysis for comparison of 22 anatomical structures between in vivo and post mortem scans. No significant changes in volumetric measurements were detected; however, as hypothesized, the image quality is dramatically improved in post mortem scans. These findings illustrate the validity and utility of using post mortem scans in volumetric neurobiological studies.