Forward Modeling Image Analysis

1. t5 t4 t3 t2 t1 time Physics-Based Constraints in the Forward Modeling Analysis of Time-Correlated Image Data James L. Carroll and Christopher D. Tomkins LANL P-21 • We tested 1 through 5 degrees of temporal freedom, and 8 different noise levels, resulting in a 5X8 matrix of experiments, each tested with 12 different numbers of optimization steps (up to 300,000), for both an independent forward modeling analysis, and a dynamic time series analysis, each averaged over 100 independent runs. Time Series Constraints Introduction We study the application of physics based constraints on the forward modeling analysis of hydrodynamic explosive data. These constraints impose a temporal regularization, restricting the space of potential reconstructions to the space of physically realizable solutions, thus improving the quality of density reconstructions produced. We study the potential applicability of this approach to experiments at Los Alamos National Laboratory’s DARHT facility with simulated data. • The forward-modeling framework makes possible a global optimization procedure Now, physics-based constraints on the evolution of the time-series data will also constrain the (global) solution Prior knowledge provides additional constraints at each time Application t1 t2 t3 t4 Results SOLUTION:Evaluated Density DATA: Transmission (experiment) • High Explosives (HE) driven experiments to study material properties under intense pressures and high velocities. • Radiographs of chosen instants during dynamic conditions. • Metals and other materials may flow like liquids under the high temperatures and pressures produced by HE. • The principle challenge is to summarize this vast body of data in order to draw meaningful conclusions. Each experiment in the 5X8 Data constrain solution at each time t1 t2 t3 t4 • array resulted in a graph like the one shown on the right. From this matrix of results, we extracted the following summary information: • These physics-based constraints will maximize information extracted from each dataset Concept: Can we learn something about the solution at time 3 (blue) from the data at surrounding times? Approach: use physics to constrain solution at each time based upon time-series of data. WHEN WILL THIS APPROACH HAVE GREATEST VALUE? When certain conditions are met: • 1) Must have the time between measurements (Dt) on the order of a relevant time scale of the flow; and 2) Must have non-perfect data (due to noise, background levels, etc). Consider an evolving interface: Static Cylinder Set-up Static Cylinder shot Static Cylinder Radiograph • Final Model Error • Advantage (the difference between the time series approach and the independent analysis approach) • The number of optimization steps before the optimization subjectively “converged” • The noise level where optimization harmed the solution instead of improving it (the “overfit point”). • This summary information is presented in the following graphs: Firing Point DARHT Phase 2: “Second Axis” Lab Space and Control Rooms Optics and Detector Bunker Data must be correlated in time Perfect data might render the information in the time series un-necessary, while noisierdata should add value to global optimization. Phase 1: “First Axis” Slide 1 DARHT is the USA’s premiere hydrodynamics testing facility, providing two intense energy radiographic views of explosives driven dynamic Final Error Methodology Physics Model: Radius Evolving through Time experiments. Axis 2, which was completed in 2009, has the capability of producing and imaging four distinct pulses. This provides a view of the dynamics of the experiments through time. The purpose of our research is to explore potential analysis techniques to make better use of this temporal time series information. Signal To Noise Ratio • Degrees of Temporal Freedom: • 1: • 2: • 3: • 4: • … • n: Forward Modeling Image Analysis Advantage Signal to Noise Ratio Conclusions Simulated Radiographs: • For all degrees of freedom: • Time series advantage first grows as the signal to noise level drops, then falls in extreme noise • Time series analysis acts as regularization, and delays the onset of overfit as noise is added • This regularization can involve the imposition of physics based constraints on the space of acceptable solutoins • Time series analysis provides improvement even when the problem is underconstrained • Time series analysis involves a more complex optimization problem than does static analysis, converging more slowly Time Step 1 Time Step 2 Time Step 3 Time Step 4 True Densities: Unobserved Starting Position Time Step 1 Time Step 2 Time Step 3 Time Step 4 Signal To Noise Levels Tested: The forward modeling analysis approach (sometimes called “analysis by synthesis”) involves optimizing the parameters of a model, in order to minimize the differences between a synthetic radiograph and the actual data. Essentially, it solves a computer vision problem through the application of computer graphics and function optimization. This approach has been shown to perform well, especially under the highly constrained conditions commonly found in scientific imaging experiments. This work is supported by The Advanced Radiography Science Campaign (C3), PM E. Mullen. LA-UR-12-26815 UNCLASSIFIED