Numerical methods for coupled reconstruction and registration in digital breast tomosynthesis.

Digital Breast Tomosynthesis (DBT) provides an insight into the fine details of normal fibroglandular tissues and abnormal lesions by reconstructing a pseudo-3D image of the breast. In this respect, DBT overcomes a major limitation of conventional X-ray mam- mography by reducing the confounding effects caused by the superposition of breast tissue. In a breast cancer screening or diagnostic context, a radiologist is interested in detecting change, which might be indicative of malignant disease. To help automate this task image registration is required to establish spatial correspondence between time points. Typically, images, such as MRI or CT, are first reconstructed and then registered. This approach can be effective if reconstructing using a complete set of data. However, for ill-posed, limited-angle problems such as DBT, estimating the deformation is com- plicated by the significant artefacts associated with the reconstruction, leading to severe inaccuracies in the registration. This paper presents a mathematical framework, which couples the two tasks and jointly estimates both image intensities and the parameters of a transformation. Under this framework, we compare an iterative method and a simultaneous method, both of which tackle the problem of comparing DBT data by combining reconstruction of a pair of temporal volumes with their registration. We evaluate our methods using various computational digital phantoms, uncom- pressed breast MR images, and in-vivo DBT simulations. Firstly, we compare both iter- ative and simultaneous methods to the conventional, sequential method using an affine transformation model. We show that jointly estimating image intensities and parametric transformations gives superior results with respect to reconstruction fidelity and regis- tration accuracy. Also, we incorporate a non-rigid B-spline transformation model into our simultaneous method. The results demonstrate a visually plausible recovery of the deformation with preservation of the reconstruction fidelity

Seismic response trends evaluation and finite element model calibration of an instrumented RC building considering soil-structure interaction and non-structural components

Peer reviewedPostprin

Earth orbital teleoperator system man-machine interface evaluation

The teleoperator system man-machine interface evaluation develops and implements a program to determine human performance requirements in teleoperator systems

Gravitational waves from black hole binary inspiral and merger: The span of third post-Newtonian effective-one-body templates

We extend the description of gravitational waves emitted by binary black holes during the final stages of inspiral and merger by introducing in the third post-Newtonian (3PN) effective-one-body (EOB) templates seven new ``flexibility'' parameters that affect the two-body dynamics and gravitational radiation emission. The plausible ranges of these flexibility parameters, notably the parameter characterising the fourth post-Newtonian effects in the dynamics, are estimated. Using these estimates, we show that the currently available standard 3PN bank of EOB templates does ``span'' the space of signals opened up by all the flexibility parameters, in that their maximized mutual overlaps are larger than 96.5%. This confirms the effectualness of 3PN EOB templates for the detection of binary black holes in gravitational-wave data from interferometric detectors. The possibility to drastically reduce the number of EOB templates using a few ``universal'' phasing functions is suggested.Comment: 23 pages, 3 figures, 4 tables, with revtex4, Minor clarifications, Final published versio

Pricing and Risk Management with High-Dimensional Quasi Monte Carlo and Global Sensitivity Analysis

We review and apply Quasi Monte Carlo (QMC) and Global Sensitivity Analysis (GSA) techniques to pricing and risk management (greeks) of representative financial instruments of increasing complexity. We compare QMC vs standard Monte Carlo (MC) results in great detail, using high-dimensional Sobol' low discrepancy sequences, different discretization methods, and specific analyses of convergence, performance, speed up, stability, and error optimization for finite differences greeks. We find that our QMC outperforms MC in most cases, including the highest-dimensional simulations and greeks calculations, showing faster and more stable convergence to exact or almost exact results. Using GSA, we are able to fully explain our findings in terms of reduced effective dimension of our QMC simulation, allowed in most cases, but not always, by Brownian bridge discretization. We conclude that, beyond pricing, QMC is a very promising technique also for computing risk figures, greeks in particular, as it allows to reduce the computational effort of high-dimensional Monte Carlo simulations typical of modern risk management.Comment: 43 pages, 21 figures, 6 table

A decentralized motion coordination strategy for dynamic target tracking

This paper presents a decentralized motion planning algorithm for the distributed sensing of a noisy dynamical process by multiple cooperating mobile sensor agents. This problem is motivated by localization and tracking tasks of dynamic targets. Our gradient-descent method is based on a cost function that measures the overall quality of sensing. We also investigate the role of imperfect communication between sensor agents in this framework, and examine the trade-offs in performance between sensing and communication. Simulations illustrate the basic characteristics of the algorithms