DSA image registration using non-uniform MRF model and pivotal control points

In order to reduce the motion artifacts in DSA, non-rigid image registration is commonly used before subtracting the mask from the contrast image. Since DSA registration requires a set of spatially non-uniform control points, a conventional MRF model is not very efficient. In this paper, we introduce the concept of pivotal and non-pivotal control points to address this, and propose a non-uniform MRF for DSA registration. We use quad-trees in a novel way to generate the non-uniform grid of control points. Our MRF formulation produces a smooth displacement field and therefore results in better artifact reduction than that of registering the control points independently. We achieve improved computational performance using pivotal control points without compromising on the artifact reduction. We have tested our approach using several clinical data sets, and have presented the results of quantitative analysis, clinical assessment and performance improvement on a GPU. (C) 2013 Elsevier Ltd. All rights reserved.

GLOBAL REGISTRATION OF MULTIPLE POINT CLOUDS USING SEMIDEFINITE PROGRAMMING

Consider N points in R-d and M local coordinate systems that are related through unknown rigid transforms. For each point, we are given (possibly noisy) measurements of its local coordinates in some of the coordinate systems. Alternatively, for each coordinate system, we observe the coordinates of a subset of the points. The problem of estimating the global coordinates of the N points (up to a rigid transform) from such measurements comes up in distributed approaches to molecular conformation and sensor network localization, and also in computer vision and graphics. The least-squares formulation of this problem, although nonconvex, has a well-known closed-form solution when M = 2 (based on the singular value decomposition (SVD)). However, no closed-form solution is known for M >= 3. In this paper, we demonstrate how the least-squares formulation can be relaxed into a convex program, namely, a semidefinite program (SDP). By setting up connections between the uniqueness of this SDP and results from rigidity theory, we prove conditions for exact and stable recovery for the SDP relaxation. In particular, we prove that the SDP relaxation can guarantee recovery under more adversarial conditions compared to earlier proposed spectral relaxations, and we derive error bounds for the registration error incurred by the SDP relaxation. We also present results of numerical experiments on simulated data to confirm the theoretical findings. We empirically demonstrate that (a) unlike the spectral relaxation, the relaxation gap is mostly zero for the SDP (i.e., we are able to solve the original nonconvex least-squares problem) up to a certain noise threshold, and (b) the SDP performs significantly better than spectral and manifold-optimization methods, particularly at large noise levels.

Diffraction Of Elastic Waves By Two Parallel Rigid Strips Embedded In An Infinite Orthotropic Medium

The elastodynamic response of a pair of parallel rigid strips embedded in an infinite orthotropic medium due to elastic waves incident normally on the strips has been investigated. The mixed boundary value problem has been solved by the Integral Equation method. The normal stress and the vertical displacement have been derived in closed form. Numerical values of stress intensity factors at inner and outer edges of the strips and vertical displacement at points in the plane of the strips for several orthotropic materials have been calculated and plotted graphically to show the effect of material orthotropy.

Analysis of a finite composite plate with smooth rigid pin

A continuum method of analysis is presented in this paper for the problem of a smooth rigid pin in a finite composite plate subjected to uniaxial loading. The pin could be of interference, push or clearance fit. The plate is idealized to an orthotropic sheet. As the load on the plate is progressively increased, the contact along the pin-hole interface is partial above certain load levels in all three types of fit. In misfit pins (interference or clearance), such situations result in mixed boundary value problems with moving boundaries and in all of them the arc of contact and the stress and displacement fields vary nonlinearly with the applied load. In infinite domains similar problems were analysed earlier by ‘inverse formulation’ and, now, the same approach is selected for finite plates. Finite outer domains introduce analytical complexities in the satisfaction of boundary conditions. These problems are circumvented by adopting a method in which the successive integrals of boundary error functions are equated to zero. Numerical results are presented which bring out the effects of the rectangular geometry and the orthotropic property of the plate. The present solutions are the first step towards the development of special finite elements for fastener joints.

Structured model-following neuro-adaptive design for attitude maneuver of rigid bodies

A new structured model-following adaptive approach is presented in this paper to achieve large attitude maneuvers of rigid bodies. First, a nominal controller is designed using the dynamic inversion philosophy. Next, a neuro- adaptive design is proposed to augment the nominal design in order to assure robust performance in the presence of parameter inaccuracies as well as unknown constant external disturbances. The structured approach proposed in this paper (where kinematic and dynamic equations are handled separately), reduces the complexity of the controller structure. From simulation studies, this adaptive controller is found to be very effective in assuring robust performance.

Coalescence of rigid droplets in a stirred dispersion-II. Band-limited force fluctuations

The coalescence of nearly rigid liquid droplets in a turbulent flow field is viewed as the drainage of a thin film of liquid under the action of a stochastic force representing the effect of turbulence. The force squeezing the drop pair is modelled as a correlated random function of time. The drops are assumed to coalesce once the film thickness becomes smaller than a critical thickness while they are regarded as separated if their distance of separation is larger than a prescribed distance. A semi-analytical solution is derived to determine the coalescence efficiency. The veracity of the solution procedure is established via a Monte-Carlo solution scheme. The model predicts a reversing trend of the dependence of the coalescence efficiency on the drop radii, the film liquid viscosity and the turbulence energy dissipation per unit mass, as the relative fluctuation increases. However, the dependence on physical parameters is weak (especially at high relative fluctuation) so that for the smallest droplets (which are nearly rigid) the coalescence efficiency may be treated as an empirical constant. The predictions of this model are compared with those of a white-noise force model. The results of this paper and those in Muralidhar and Ramkrishna (1986, Ind. Engng Chem. Fundam. 25, 554-56) suggest that dynamic drop deformation is the key factor that influences the coalescence efficiency.

Pulsatile blood flow in a rigid pulmonary alveolar sheet with porous walls

This paper deals with the pulsatile blood flow in the lung alveolar sheets by idealizing each of them as a channel covered by porous media. As the blood flow in the lung is of low Reynolds number, a creeping flow is assumed in the channel. The analytical and numerical results for the velocity and pressure distribution in the porous medium are presented. The effect of an imposed slip condition is also studied. Comparisons with the corresponding results for the steady-state case are made at the end.

Is 3'-nucleotide rigid?

Conformational analysis of nucleic acids and polynucleotides is far more complex than that of proteins and polypeptides, due to five single bond rotations in addition to the sugar puckerings in the monomer. Sundaralingam1 proposed the concept of the 'rigid' nucleotides from analysis of crystal structure data, with the flexibility allowed only about the phosphodiester bonds. However, the crystal structure of deoxyguanosine-5'−phosphate2,3 indicates at gt conformation about the C-4'−C-5' bond against gg in a conformationally rigid nudeotide1. Jack et al. 4 considered the flexibility of nucleotides in tRNA about the C-4'−C-5' bond, thereby introducing the concept of 'non-rigid' ribonucleotides. Conformational flexibility of the f uranose ring in DNA and RNA and their energetics using classical and quantum chemical methods have been reported5−8. We have examined the flexibility of 3'-nucleotides. alpha, the most important of the conformational parameters defining the 3'-end of a nucleotide unit9, has a value in the range 195°−270° in all the 3'-nucleotides, dinucleoside monophosphates and higher oligomers which have been surveyed. A survey of the proposed structures of polyribonudeotides10,11 also shows the values of a to be greater than 200°. However, the structures proposed for B-DNA by Arnott and Hukins12,13 and D-DNA by Arnott et al. 14 have values of alpha of 155° and 141° respectively, much lower than the lowest observed value. The structure for B-DNA has two strong, short contacts (C-2'...OP-1 = 2.64 Å and HC-2"...OP-1 = 1.79 Å) which lead to an energetically unfavourable conformation. Hence, it is of interest to investigate whether, by allowing flexibility to the sugar moiety in the nucleotide unit, it is possible to make the structure energetically favourable. Here, conformational energy calculations were carried out to determine the range of alpha which would give rise to energetically favoured conformations with different sugar puckerings. Our analysis has shown that the theoretically obtained range is nearly the same as the preferred range in crystals, indicating the flexibility of the 3'-nucleotides.

Seismic behaviour of rigid-faced reinforced soil retaining wall models:reinforcement effect

This paper presents the results of shaking table tests on models of rigid-faced reinforced soil retaining walls in which reinforcement materials of different tensile strength were used. The construction of the model retaining walls in a laminar box mounted on a shaking table, the instrumentation and the results from the shaking table tests are described in detail and the effects of the reinforcement parameters on the acceleration response at different elevations of the retaining wall, horizontal soil pressures and face deformations are presented. It was observed from these tests that the horizontal face displacement response of the rigid-faced retaining walls was significantly affected by the inclusion of reinforcement and even low-strength polymer reinforcement was found to be efficient in significantly reducing the deformation of the face. The acceleration amplifications were, however, observed to be less influenced by the reinforcement parameters. The results obtained from this study are helpful in understanding the relative performance of reinforced soil retaining walls under the different test conditions used in the experiments.

An improved representation of vehicle incompatibility in frontal NCAP tests using a modified rigid barrier

With the objective of better understanding the significance of New Car Assessment Program (NCAP) tests conducted by the National Highway Traffic Safety Administration (NHTSA), head-on collisions between two identical cars of different sizes and between cars and a pickup truck are studied in the present paper using LS-DYNA models. Available finite element models of a compact car (Dodge Neon), midsize car (Dodge Intrepid), and pickup truck (Chevrolet C1500) are first improved and validated by comparing theanalysis-based vehicle deceleration pulses against corresponding NCAP crash test histories reported by NHTSA. In confirmation of prevalent perception, simulation-bascd results indicate that an NCAP test against a rigid barrier is a good representation of a collision between two similar cars approaching each other at a speed of 56.3 kmph (35 mph) both in terms of peak deceleration and intrusions. However, analyses carried out for collisions between two incompatible vehicles, such as an Intrepid or Neon against a C1500, point to the inability of the NCAP tests in representing the substantially higher intrusions in the front upper regions experienced by the cars, although peak decelerations in cars arc comparable to those observed in NCAP tests. In an attempt to improve the capability of a front NCAP test to better represent real-world crashes between incompatible vehicles, i.e., ones with contrasting ride height and lower body stiffness, two modified rigid barriers are studied. One of these barriers, which is of stepped geometry with a curved front face, leads to significantly improved correlation of intrusions in the upper regions of cars with respect to those yielded in the simulation of collisions between incompatible vehicles, together with the yielding of similar vehicle peak decelerations obtained in NCAP tests.

Analysis of a thick plate with a circular hole resting on a smooth rigid bed and subjected to axisymmetric normal load

A solution for the stresses and displacements in an radially infinite thick plate having a circular hole, one face of which resting on a smooth rigid bed and the other face subjected to axisymmetric normal loading is given. The solution is obtained in terms of Fourier-Bessel series and integral for the Love's stress function. Numerical results are presented for one particular ratio of thickness of plate to the hole radius and loading. It is also shown that the Poisson's ratio has a predominant effect on certain stresses and displacements. The solution would be useful in the stress analysis of bolted joints.Eine Lösung für die Spannungen und Verschiebungen in einer radial, unendlich ausgedehnten, dicken Platte mit einem kreisförmigen Loch, wobei eine Seite auf einer ebenen, starren Unterlage aufliegt, die andere Seite durch eine achsensymmetrische Vertikallast belastet ist, wird angegeben. Die Lösung wird in Form von Fourier-Bessel-Reihen und Integralen der Loveschen Spannungsfunktion angegeben. Numerische Ergebnisse werden für ein bestimmtes Verhältnis der Plattendicke zum Lochradius sowie zur Belastung angegeben. Es wird auch gezeigt, daß das Poisssonsche Verhältnis einen besonderen Einfluß auf bestimmte Spannungen und Verschiebungen hat. Die Lösung ist anwendbar für die Spannungsermittlung von Bolzenverbindungen.

Effect of interface friction on the mechanics of indentation of a finite layer resting on a rigid substrate

Copper strips of 2.5 mm thickness resting on stainless steel anvils were normally indented by wedges under nominal plane strain conditions. Inflections in the hardness-penetration characteristics were identified. Inflections separate stages where each stage has typical mechanics of deformation. These are arrived at by studying the distortion of 0.125 mm spaced grids inscribed on the deformation plane of the strip. The sensitivity of hardness and deformation mechanics to wedge angle and the interfacial friction between strip and anvil were investigated within the framework of existing slip line field models of indentation of semi-infinite and finite blocks.

Rigid and flexible region in lysozyme and the invariant features in its hydration shell

Water-mediated transformations provide a useful handle for exploring the flexibility in protein molecules and the invariant features in their hydration shells. Low-humidity monoclinic hen egg white lysozyme, resulting from such a transformation, has perhaps the lowest solvent content observed in any protein crystal so far and has a well-ordered structure. A detailed comparison involving this structure, low-humidity tetragonal lysozyme, and the other available refined crystal structures of the enzyme permits the delineation of the relatively rigid, moderately flexible and highly flexible regions of the molecule. The relatively rigid region forms a contiguous structural unit close to the molecular centroid and encompasses parts of of the main beta-structure and three alpha-helices. The hydration shell of the protein contains 30 invariant water molecules. Many of them are involved in holding different parts of the molecule together or in stabilizing local structure. Five of the six invariant water molecules attached to the substrate-binding region form part of a water cluster contiguous with the side-chains of the catalytic residues Glu-35 and Asp-52.

Rocking response of rectangular rigid blocks under random noise base excitations

The response of a rigid rectangular block resting on a rigid foundation and acted upon simultaneously by a horizontal and a vertical random white-noise excitation is considered. In the equation of motion, the energy dissipation is modeled through a viscous damping term. Under the assumption that the body does not topple, the steady-state joint probability density function of the rotation and the rotational velocity is obtained using the Fokker-Planck equation approach. Closed form solution is obtained for a specific combination of system parameters. A more general but approximate solution to the joint probability density function based on the method of equivalent non-linearization is also presented. Further, the problem of overturning of the block is approached in the framework of the diffusion methods for first passage failure studies. The overturning of the block is deemed incipient when the response trajectories in the phase plane cross the separatrix of the conservative unforced system. Expressions for the moments of first passage time are obtained via a series solution to the governing generalized Pontriagin-Vitt equations. Numerical results illustra- tive of the theoretical solutions are presented and their validity is examined through limited amount of digital simulations.