Simultaneous Multiscale Polyaffine Registration by Incorporating Deformation Statistics

Locally affine (polyaffine) image registration methods capture intersubject non-linear deformations with a low number of parameters, while providing an intuitive interpretation for clinicians. Considering the mandible bone, anatomical shape differences can be found at different scales, e.g. left or right side, teeth, etc. Classically, sequential coarse to fine registration are used to handle multiscale deformations, instead we propose a simultaneous optimization of all scales. To avoid local minima we incorporate a prior on the polyaffine transformations. This kind of groupwise registration approach is natural in a polyaffine context, if we assume one configuration of regions that describes an entire group of images, with varying transformations for each region. In this paper, we reformulate polyaffine deformations in a generative statistical model, which enables us to incorporate deformation statistics as a prior in a Bayesian setting. We find optimal transformations by optimizing the maximum a posteriori probability. We assume that the polyaffine transformations follow a normal distribution with mean and concentration matrix. Parameters of the prior are estimated from an initial coarse to fine registration. Knowing the region structure, we develop a blockwise pseudoinverse to obtain the concentration matrix. To our knowledge, we are the first to introduce simultaneous multiscale optimization through groupwise polyaffine registration. We show results on 42 mandible CT images.

A three-dimensional histological atlas of the human basal ganglia. II. Atlas deformation strategy and evaluation in deep brain stimulation for Parkinson disease

OBJECT: The localization of any given target in the brain has become a challenging issue because of the increased use of deep brain stimulation to treat Parkinson disease, dystonia, and nonmotor diseases (for example, Tourette syndrome, obsessive compulsive disorders, and depression). The aim of this study was to develop an automated method of adapting an atlas of the human basal ganglia to the brains of individual patients. METHODS: Magnetic resonance images of the brain specimen were obtained before extraction from the skull and histological processing. Adaptation of the atlas to individual patient anatomy was performed by reshaping the atlas MR images to the images obtained in the individual patient using a hierarchical registration applied to a region of interest centered on the basal ganglia, and then applying the reshaping matrix to the atlas surfaces. RESULTS: Results were evaluated by direct visual inspection of the structures visible on MR images and atlas anatomy, by comparison with electrophysiological intraoperative data, and with previous atlas studies in patients with Parkinson disease. The method was both robust and accurate, never failing to provide an anatomically reliable atlas to patient registration. The registration obtained did not exceed a 1-mm mismatch with the electrophysiological signatures in the region of the subthalamic nucleus. CONCLUSIONS: This registration method applied to the basal ganglia atlas forms a powerful and reliable method for determining deep brain stimulation targets within the basal ganglia of individual patients.

Quantification of Popliteal Artery Deformation During Leg Flexion in Subjects With Peripheral Artery Disease: A Pilot Study

Purpose: To quantify the in vivo deformations of the popliteal artery during leg flexion in subjects with clinically relevant peripheral artery disease (PAD). Methods: Five patients (4 men; mean age 69 years, range 56–79) with varying calcification levels of the popliteal artery undergoing endovascular revascularization underwent 3-dimensional (3D) rotational angiography. Image acquisition was performed with the leg straight and with a flexion of 70°/20° in the knee/hip joints. The arterial centerline and the corresponding branches in both positions were segmented to create 3D reconstructions of the arterial trees. Axial deformation, twisting, and curvatures were quantified. Furthermore, the relationships between the calcification levels and the deformations were investigated. Results: An average shortening of 5.9%±2.5% and twist rate of 3.8±2.2°/cm in the popliteal artery were observed. Maximal curvatures in the straight and flexed positions were 0.12±0.04 cm−1 and 0.24±0.09 cm−1, respectively. As the severity of calcification increased, the maximal curvature in the straight position increased from 0.08 to 0.17 cm−1, while an increase from 0.17 to 0.39 cm−1 was observed for the flexed position. Axial elongations and arterial twisting were not affected by the calcification levels. Conclusion: The popliteal artery of patients with symptomatic PAD is exposed to significant deformations during flexion of the knee joint. The severity of calcification directly affects curvature, but not arterial length or twisting angles. This pilot study also showed the ability of rotational angiography to quantify the 3D deformations of the popliteal artery in patients with various levels of calcification.

Characterization of granite matrix porosity and pore-space geometry by in situ and laboratory methods

Most available studies of interconnected matrix porosity of crystalline rocks are based on laboratory investigations; that is, work on samples that have undergone stress relaxation and were affected by drilling and sample preparation. The extrapolation of the results to in situ conditions is therefore associated with considerable uncertainty, and this was the motivation to conduct the ‘in situ Connected Porosity’ experiment at the Grimsel Test Site (Central Swiss Alps). An acrylic resin doped with fluorescent agents was used to impregnate the microporous granitic matrix in situ around an injection borehole, and samples were obtained by overcoring. The 3-D structure of the porespace, represented by microcracks, was studied by U-stage fluorescence microscopy. Petrophysical methods, including the determination of porosity, permeability and P -wave velocity, were also applied. Investigations were conducted both on samples that were impregnated in situ and on non-impregnated samples, so that natural features could be distinguished from artefacts. The investigated deformed granites display complex microcrack populations representing a polyphase deformation at varying conditions. The crack population is dominated by open cleavage cracks in mica and grain boundary cracks. The porosity of non-impregnated samples lies slightly above 1 per cent, which is 2–2.5 times higher than the in situ porosity obtained for impregnated samples. Measurements of seismic velocities (Vp ) on spherical rock samples as a function of confining pressure, spatial direction and water saturation for both non-impregnated and impregnated samples provide further constraints on the distinction between natural and induced crack types. The main conclusions are that (1) an interconnected network of microcracks exists in the whole granitic matrix, irrespective of the distance to ductile and brittle shear zones, and (2) conventional laboratory methods overestimate the matrix porosity. Calculations of contaminant transport through fractured media often rely on matrix diffusion as a retardation mechanism.