Entropy-based particle correspondence for shape populations

Enabling supra-aortic vessels inclusion in statistical shape models of the aorta: a novel non-rigid registration method

Statistical Shape Models (SSMs) are well-established tools for assessing the variability of 3D geometry and for broadening a limited set of shapes. They are widely used in medical imaging due to their ability to model complex geometries and their high efficiency as generative models. The principal step behind these techniques is a registration phase, which, in the case of complex geometries, can be a critical issue due to the correspondence problem, as it necessitates the development of correspondence mapping between shapes. The thoracic aorta, with its high level of morphological complexity, poses a multi-scale deformation problem due to the presence of several branch vessels with varying diameters. Moreover, branch vessels exhibit significant variability in shape, making the correspondence optimization even more challenging. Consequently, existing studies have focused on developing SSMs based only on the main body of the aorta, excluding the supra-aortic vessels from the analysis. In this work, we present a novel non-rigid registration algorithm based on optimizing a differentiable distance function through a modified gradient descent approach. This strategy enables the inclusion of custom, domain-specific constraints in the objective function, which act as landmarks during the registration phase. The algorithm's registration performance was tested and compared to an alternative Statistical Shape modeling framework, and subsequently used for the development of a comprehensive SSM of the thoracic aorta, including the supra-aortic vessels. The developed SSM was further evaluated against the alternative framework in terms of generalisation, specificity, and compactness to assess its effectiveness.

Feasibility study for the automatic surgical planning method based on statistical model

PURPOSE

In this study, we proposed establishing an automatic computer-assisted surgical planning approach based on average population models.

METHODS

We built the average population models from humerus datasets using the Advanced Normalization Toolkits (ANTs) and Shapeworks. Experiments include (1) evaluation of the average population models before surgical planning and (2) validation of the average population models in the context of predicting clinical landmarks on the humerus from the new dataset that was not involved in the process of building the average population model. The evaluation experiment consists of explained variation and distance model. The validation experiment calculated the root-mean-square error (RMSE) between the expert-determined clinical ground truths and the landmarks transferred from the average population model to the new dataset. The evaluation results and validation results when using the templates built from ANTs were compared to when using the mean shape generated from Shapeworks.

RESULTS

The average population models predicted clinical locations on the new dataset with acceptable errors when compared to the ground truth determined by an expert. However, the templates built from ANTs present better accuracy in landmark prediction when compared to the mean shape built from the Shapeworks.

CONCLUSION

The average population model could be utilized to assist anatomical landmarks checking automatically and following surgical decisions for new patients who are not involved in the dataset used to generate the average population model.

Establishing Surface Correspondence for Post-surgical Cortical Thickness Changes in Temporal Lobe Epilepsy

In pre- and post-surgical surface shape analysis, establishing shape correspondence is necessary to investigate the postoperative surface changes. However, structural absence after the operation accompanies focal non-rigid changes, which leads to challenges in existing surface registration methods. In this paper, we present a fully automatic particle-based method to establish surface correspondence that can handle partial structural abnormality in the temporal lobe resection. Our method optimizes the coordinates of points which are modeled as particles on surfaces in a hierarchical way to reduce a chance of being trapped in a local minimum during the optimization. In the experiments, we evaluate the effectiveness of our method in comparison with conventional spherical registration (FreeSurfer) on two scenarios: cortical thickness changes in healthy controls within a short scan-rescan time window and patients with temporal lobe resection. The post-surgical scan is acquired at least 1 year after the presurgical scan. In region of interest-wise (ROI-wise) analysis, no changes on cortical thickness are found in both methods on the healthy control group. In patients, since there is no ground truth available, we instead investigated the disagreement between our method and FreeSurfer. We see poorly matched ROIs and large cortical thickness changes using FreeSurfer. On the contrary, our method shows well-matched ROIs and subtle cortical thickness changes. This suggests that the proposed method can establish a stable shape correspondence, which is not fully captured in a conventional spherical registration.

Comparative analysis of squamate brains unveils multi-level variation in cerebellar architecture associated with locomotor specialization

Ecomorphological studies evaluating the impact of environmental and biological factors on the brain have so far focused on morphology or size measurements, and the ecological relevance of potential multi-level variations in brain architecture remains unclear in vertebrates. Here, we exploit the extraordinary ecomorphological diversity of squamates to assess brain phenotypic diversification with respect to locomotor specialization, by integrating single-cell distribution and transcriptomic data along with geometric morphometric, phylogenetic, and volumetric analysis of high-definition 3D models. We reveal significant changes in cerebellar shape and size as well as alternative spatial layouts of cortical neurons and dynamic gene expression that all correlate with locomotor behaviours. These findings show that locomotor mode is a strong predictor of cerebellar structure and pattern, suggesting that major behavioural transitions in squamates are evolutionarily correlated with mosaic brain changes. Furthermore, our study amplifies the concept of ‘cerebrotype’, initially proposed for vertebrate brain proportions, towards additional shape characters.

The cerebellum is critical in sensory-motor control and is structurally diverse across vertebrates. Here, the authors investigate the evolutionary relationship between locomotory mode and cerebellum architecture across squamates by integrating study of gene expression, cell distribution, and 3D morphology.

Surface-Based Spatial Pyramid Matching of Cortical Regions for Analysis of Cognitive Performance

We propose a method to analyze the relationship between the shape of functional regions of the cortex and cognitive measures, such as reading ability and vocabulary knowledge. Functional regions on the cortical surface can vary not only in size and shape but also in topology and position relative to neighboring regions. Standard diffeomorphism-based shape analysis tools do not work well here because diffeomorphisms are unable to capture these topological differences, which include region splitting and merging across subjects. State-of-the-art cortical surface shape analyses compute derived regional properties (scalars), such as regional volume, cortical thickness, curvature, and gyrification index. However, these methods cannot compare the full extent of topological or shape differences in cortical regions. We propose icosahedral spatial pyramid matching (ISPM) of region borders computed on the surface of a sphere to capture this variation in regional topology, position, and shape. We then analyze how this variation corresponds to measures of cognitive performance. We compare our method to other approaches and find that it is indeed informative to consider aspects of shape beyond the standard approaches. Analysis is performed using a subset of 27 test/retest subjects from the Human Connectome Project in order to understand both the effectiveness and reproducibility of this method.

On the Evaluation and Validation of Off-the-shelf Statistical Shape Modeling Tools: A Clinical Application

Statistical shape modeling (SSM) has proven useful in many areas of biology and medicine as a new generation of morphometric approaches for the quantitative analysis of anatomical shapes. Recently, the increased availability of high-resolution in vivo images of anatomy has led to the development and distribution of open-source computational tools to model anatomical shapes and their variability within populations with unprecedented detail and statistical power. Nonetheless, there is little work on the evaluation and validation of such tools as related to clinical applications that rely on morphometric quantifications for treatment planning. To address this lack of validation, we systematically assess the outcome of widely used off-the-shelf SSM tools, namely ShapeWorks, SPHARM-PDM, and Deformetrica, in the context of designing closure devices for left atrium appendage (LAA) in atrial fibrillation (AF) patients to prevent stroke, where an incomplete LAA closure may be worse than no closure. This study is motivated by the potential role of SSM in the geometric design of closure devices, which could be informed by population-level statistics, and patient-specific device selection, which is driven by anatomical measurements that could be automated by relating patient-level anatomy to population-level morphometrics. Hence, understanding the consequences of different SSM tools for the final analysis is critical for the careful choice of the tool to be deployed in real clinical scenarios. Results demonstrate that estimated measurements from ShapeWorks model are more consistent compared to models from Deformetrica and SPHARM-PDM. Furthermore, ShapeWorks and Deformetrica shape models capture clinically relevant population-level variability compared to SPHARM-PDM models.

Hierarchical Vertex Regression-Based Segmentation of Head and Neck CT Images for Radiotherapy Planning

Segmenting organs at risk from head and neck CT images is a prerequisite for the treatment of head and neck cancer using intensity modulated radiotherapy. However, accurate and automatic segmentation of organs at risk is a challenging task due to the low contrast of soft tissue and image artifact in CT images. Shape priors have been proved effective in addressing this challenging task. However, conventional methods incorporating shape priors often suffer from sensitivity to shape initialization and also shape variations across individuals. In this paper, we propose a novel approach to incorporate shape priors into a hierarchical learning-based model. The contributions of our proposed approach are as follows: 1) a novel mechanism for critical vertices identification is proposed to identify vertices with distinctive appearances and strong consistency across different subjects; 2) a new strategy of hierarchical vertex regression is also used to gradually locate more vertices with the guidance of previously located vertices; and 3) an innovative framework of joint shape and appearance learning is further developed to capture salient shape and appearance features simultaneously. Using these innovative strategies, our proposed approach can essentially overcome drawbacks of the conventional shape-based segmentation methods. Experimental results show that our approach can achieve much better results than state-of-the-art methods.

Genetic mapping of molar size relations identifies inhibitory locus for third molars in mice

Molar size in Mammals shows considerable disparity and exhibits variation similar to that predicted by the Inhibitory Cascade model. The importance of such developmental systems in favoring evolutionary trajectories is also underlined by the fact that this model can predict macroevolutionary patterns. Using backcross mice, we mapped QTL for molar sizes controlling for their sequential development. Genetic controls for upper and lower molars appear somewhat similar, and regions containing genes implied in dental defects drive this variation. We mapped three relationship QTLs (rQTL) modifying the control of the mesial molars on the focal third molar. These regions overlap Shh, Sostdc1, and Fst genes, which have pervasive roles in development and should be buffered against new variation. It has theoretically been shown that rQTL produces new variation channeled in the direction of adaptive changes. Our results provide evidence that evolutionary/disease patterns of tooth size variation could result from such a non-random generating process.

Automatic construction of statistical shape models using deformable simplex meshes with vector field convolution energy

BACKGROUND

In the active shape model framework, principal component analysis (PCA) based statistical shape models (SSMs) are widely employed to incorporate high-level a priori shape knowledge of the structure to be segmented to achieve robustness. A crucial component of building SSMs is to establish shape correspondence between all training shapes, which is a very challenging task, especially in three dimensions.

METHODS

We propose a novel mesh-to-volume registration based shape correspondence establishment method to improve the accuracy and reduce the computational cost. Specifically, we present a greedy algorithm based deformable simplex mesh that uses vector field convolution as the external energy. Furthermore, we develop an automatic shape initialization method by using a Gaussian mixture model based registration algorithm, to derive an initial shape that has high overlap with the object of interest, such that the deformable models can then evolve more locally. We apply the proposed deformable surface model to the application of femur statistical shape model construction to illustrate its accuracy and efficiency.

RESULTS

Extensive experiments on ten femur CT scans show that the quality of the constructed femur shape models via the proposed method is much better than that of the classical spherical harmonics (SPHARM) method. Moreover, the proposed method achieves much higher computational efficiency than the SPHARM method.

CONCLUSIONS

The experimental results suggest that our method can be employed for effective statistical shape model construction.