Real-Time Joint Coupling of the Spine for Inverse Kinematics

In this paper we propose a simple model for the coupling behavior of the human spine for an inverse kinematics framework. Our spine model exhibits anatomically correct motions of the vertebrae of virtual mannequins by coupling standard swing and revolute joint models. The adjustement of the joints is made with several simple (in)equality constraints, resulting in a reduction of the solution space dimensionality for the inverse kinematics solver. By reducing the solution space dimensionality to feasible spine shapes, we prevent the inverse kinematics algorithm from providing infeasible postures for the spine.In this paper, we exploit how to apply these simple constraints to the human spine by a strict decoupling of the swing and torsion motion of the vertebrae. We demonstrate the validity of our approach on various experiments.

Analytical form of the Probability Density Function of travel times in rectangular zone with Tchebyshev’s metric

Geometrical dependencies are being researched for analytical representation of the probability density function (pdf) for the travel time between a random, and a known or another random point in Tchebyshev’s metric. In the most popular case - a rectangular area of service - the pdf of this random variable depends directly on the position of the server. Two approaches have been introduced for the exact analytical calculation of the pdf: Ad-hoc approach – useful for a ‘manual’ solving of a specific case; by superposition – an algorithmic approach for the general case. The main concept of each approach is explained, and a short comparison is done to prove the faithfulness.

Distance-based kriging relying on proxy simulations for inverse conditioning

Let us consider a large set of candidate parameter fields, such as hydraulic conductivity maps, on which we can run an accurate forward flow and transport simulation. We address the issue of rapidly identifying a subset of candidates whose response best match a reference response curve. In order to keep the number of calls to the accurate flow simulator computationally tractable, a recent distance-based approach relying on fast proxy simulations is revisited, and turned into a non-stationary kriging method where the covariance kernel is obtained by combining a classical kernel with the proxy. Once the accurate simulator has been run for an initial subset of parameter fields and a kriging metamodel has been inferred, the predictive distributions of misfits for the remaining parameter fields can be used as a guide to select candidate parameter fields in a sequential way. The proposed algorithm, Proxy-based Kriging for Sequential Inversion (ProKSI), relies on a variant of the Expected Improvement, a popular criterion for kriging-based global optimization. A statistical benchmark of ProKSI’s performances illustrates the efficiency and the robustness of the approach when using different kinds of proxies.

Inverse relationship between fractionated electrograms and atrial fibrosis in persistent atrial fibrillation: combined magnetic resonance imaging and high-density mapping

OBJECTIVES This study sought to evaluate the relationship between fibrosis imaged by delayed-enhancement (DE) magnetic resonance imaging (MRI) and atrial electrograms (Egms) in persistent atrial fibrillation (AF). BACKGROUND Atrial fractionated Egms are strongly related to slow anisotropic conduction. Their relationship to atrial fibrosis has not yet been investigated. METHODS Atrial high-resolution MRI of 18 patients with persistent AF (11 long-lasting persistent AF) was registered with mapping geometry (NavX electro-anatomical system (version 8.0, St. Jude Medical, St. Paul, Minnesota)). DE areas were categorized as dense or patchy, depending on their DE content. Left atrial Egms during AF were acquired using a high-density, 20-pole catheter (514 ± 77 sites/map). Fractionation, organization/regularity, local mean cycle length (CL), and voltage were analyzed with regard to DE. RESULTS Patients with long-lasting persistent versus persistent AF had larger left atrial (LA) surface area (134 ± 38 cm(2) vs. 98 ± 9 cm(2), p = 0.02), a higher amount of atrial DE (70 ± 16 cm(2) vs. 49 ± 10 cm(2), p = 0.01), more complex fractionated atrial Egm (CFAE) extent (54 ± 16 cm(2) vs. 28 ± 15 cm(2), p = 0.02), and a shorter baseline AF CL (147 ± 10 ms vs. 182 ± 14 ms, p = 0.01). Continuous CFAE (CFEmean [NavX algorithm that quantifies Egm fractionation] <80 ms) occupied 38 ± 19% of total LA surface area. Dense DE was detected at the left posterior left atrium. In contrast, the right posterior left atrium contained predominantly patchy DE. Most CFAE (48 ± 14%) occurred at non-DE LA sites, followed by 41 ± 12% CFAE at patchy DE and 11 ± 6% at dense DE regions (p = 0.005 and p = 0.008, respectively); 19 ± 6% CFAE sites occurred at border zones of dense DE. Egms were less fractionated, with longer CL and lower voltage at dense DE versus non-DE regions: CFEmean: 97 ms versus 76 ms, p < 0.0001; local CL: 153 ms versus 143 ms, p < 0.0001; mean voltage: 0.63 mV versus 0.86 mV, p < 0.0001. CONCLUSIONS Atrial fibrosis as defined by DE MRI is associated with slower and more organized electrical activity but with lower voltage than healthy atrial areas. Ninety percent of continuous CFAE sites occur at non-DE and patchy DE LA sites. These findings are important when choosing the ablation strategy in persistent AF.

Prediction of Cranio-Maxilofacial Surgical Planning using an Inverse Soft Tissue Modelling Approach

In cranio-maxillofacial surgery, the determination of a proper surgical plan is an important step to attain a desired aesthetic facial profile and a complete denture closure. In the present paper, we propose an efficient modeling approach to predict the surgical planning on the basis of the desired facial appearance and optimal occlusion. To evaluate the proposed planning approach, the predicted osteotomy plan of six clinical cases that underwent CMF surgery were compared to the real clinical plan. Thereafter, simulated soft-tissue outcomes were compared using the predicted and real clinical plan. This preliminary retrospective comparison of both osteotomy planning and facial outlook shows a good agreement and thereby demonstrates the potential application of the proposed approach in cranio-maxillofacial surgical planning prediction.

Solute transport in crystalline rocks at Äspö — II: Blind predictions, inverse modelling and lessons learnt from test STT1

Based on the results from detailed structural and petrological characterisation and on up-scaled laboratory values for sorption and diffusion, blind predictions were made for the STT1 dipole tracer test performed in the Swedish A¨ spo¨ Hard Rock Laboratory. The tracers used were nonsorbing, such as uranine and tritiated water, weakly sorbing 22Na+, 85Sr2 +, 47Ca2 +and more strongly sorbing 86Rb+, 133Ba2 +, 137Cs+. Our model consists of two parts: (1) a flow part based on a 2D-streamtube formalism accounting for the natural background flow field and with an underlying homogeneous and isotropic transmissivity field and (2) a transport part in terms of the dual porosity medium approach which is linked to the flow part by the flow porosity. The calibration of the model was done using the data from one single uranine breakthrough (PDT3). The study clearly showed that matrix diffusion into a highly porous material, fault gouge, had to be included in our model evidenced by the characteristic shape of the breakthrough curve and in line with geological observations. After the disclosure of the measurements, it turned out that, in spite of the simplicity of our model, the prediction for the nonsorbing and weakly sorbing tracers was fairly good. The blind prediction for the more strongly sorbing tracers was in general less accurate. The reason for the good predictions is deemed to be the result of the choice of a model structure strongly based on geological observation. The breakthrough curves were inversely modelled to determine in situ values for the transport parameters and to draw consequences on the model structure applied. For good fits, only one additional fracture family in contact with cataclasite had to be taken into account, but no new transport mechanisms had to be invoked. The in situ values for the effective diffusion coefficient for fault gouge are a factor of 2–15 larger than the laboratory data. For cataclasite, both data sets have values comparable to laboratory data. The extracted Kd values for the weakly sorbing tracers are larger than Swedish laboratory data by a factor of 25–60, but agree within a factor of 3–5 for the more strongly sorbing nuclides. The reason for the inconsistency concerning Kds is the use of fresh granite in the laboratory studies, whereas tracers in the field experiments interact only with fracture fault gouge and to a lesser extent with cataclasite both being mineralogically very different (e.g. clay-bearing) from the intact wall rock.

Dirichlet-Neumann inverse spectral problem for a star graph of Stieltjes strings

We solve two inverse spectral problems for star graphs of Stieltjes strings with Dirichlet and Neumann boundary conditions, respectively, at a selected vertex called root. The root is either the central vertex or, in the more challenging problem, a pendant vertex of the star graph. At all other pendant vertices Dirichlet conditions are imposed; at the central vertex, at which a mass may be placed, continuity and Kirchhoff conditions are assumed. We derive conditions on two sets of real numbers to be the spectra of the above Dirichlet and Neumann problems. Our solution for the inverse problems is constructive: we establish algorithms to recover the mass distribution on the star graph (i.e. the point masses and lengths of subintervals between them) from these two spectra and from the lengths of the separate strings. If the root is a pendant vertex, the two spectra uniquely determine the parameters on the main string (i.e. the string incident to the root) if the length of the main string is known. The mass distribution on the other edges need not be unique; the reason for this is the non-uniqueness caused by the non-strict interlacing of the given data in the case when the root is the central vertex. Finally, we relate of our results to tree-patterned matrix inverse problems.

Symmetric inverse consistent nonlinear registration driven by mutual information.

A nonlinear viscoelastic image registration algorithm based on the demons paradigm and incorporating inverse consistent constraint (ICC) is implemented. An inverse consistent and symmetric cost function using mutual information (MI) as a similarity measure is employed. The cost function also includes regularization of transformation and inverse consistent error (ICE). The uncertainties in balancing various terms in the cost function are avoided by alternatively minimizing the similarity measure, the regularization of the transformation, and the ICE terms. The diffeomorphism of registration for preventing folding and/or tearing in the deformation is achieved by the composition scheme. The quality of image registration is first demonstrated by constructing brain atlas from 20 adult brains (age range 30-60). It is shown that with this registration technique: (1) the Jacobian determinant is positive for all voxels and (2) the average ICE is around 0.004 voxels with a maximum value below 0.1 voxels. Further, the deformation-based segmentation on Internet Brain Segmentation Repository, a publicly available dataset, has yielded high Dice similarity index (DSI) of 94.7% for the cerebellum and 74.7% for the hippocampus, attesting to the quality of our registration method.

Sizing up your enemy: individual predation vulnerability predicts migratory probability

Partial migration, in which a fraction of a population migrate and the rest remain resident, occurs in an extensive range of species and can have powerful ecological consequences. The question of what drives differences in individual migratory tendency is a contentious one. It has been shown that the timing of partial migration is based upon a trade-off between seasonal fluctuations in predation risk and growth potential. Phenotypic variation in either individual predation risk or growth potential should thus mediate the strength of the trade-off and ultimately predict patterns of partial migration at the individual level (i.e. which individuals migrate and which remain resident). We provide cross-population empirical support for the importance of one component of this model—individual predation risk—in predicting partial migration in wild populations of bream Abramis brama, a freshwater fish. Smaller, high-risk individuals migrate with a higher probability than larger, low-risk individuals, and we suggest that predation risk maintains size-dependent partial migration in this system.