Implicit Neural Representations for Simultaneous Reduction and Continuous Reconstruction of Multi-Altitude Climate Data
Explore the source record for details and available documents.
97 MATHEMATICS AND COMPUTING↗
SEARCH · Search NASA
Results for “continuous reconstruction”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
A continuous surface reconstruction method on point cloud captured from a 3D surface photogrammetry system
Purpose: To accurately and efficiently reconstruct a continuous surface from noisy point clouds captured by a surface photogrammetry system (VisionRT). Methods: The authors have developed a level-set based surface reconstruction method on point clouds captured by a surface photogrammetry system (VisionRT). The proposed method reconstructs an implicit and continuous representation of the underlying patient surface by optimizing a regularized fitting energy, offering extra robustness to noise and missing measurements. By contrast to explicit/discrete meshing-type schemes, their continuous representation is particularly advantageous for subsequent surface registration and motion tracking by eliminating the need for maintaining explicit point correspondences as in discrete models. The authors solve the proposed method with an efficient narrowband evolving scheme. The authors evaluated the proposed method on both phantom and human subject data with two sets of complementary experiments. In the first set of experiment, the authors generated a series of surfaces each with different black patches placed on one chest phantom. The resulting VisionRT measurements from the patched area had different degree of noise and missing levels, since VisionRT has difficulties in detecting dark surfaces. The authors applied the proposed method to point clouds acquired under these different configurations, and quantitatively evaluated reconstructed surfaces by comparing against a high-quality reference surface with respect to root mean squared error (RMSE). In the second set of experiment, the authors applied their method to 100 clinical point clouds acquired from one human subject. In the absence of ground-truth, the authors qualitatively validated reconstructed surfaces by comparing the local geometry, specifically mean curvature distributions, against that of the surface extracted from a high-quality CT obtained from the same patient. Results: On phantom point clouds, their method achieved submillimeter reconstruction RMSE under different configurations, demonstrating quantitatively the faith of the proposed method in preserving local structural properties of the underlying surface in the presence of noise and missing measurements, and its robustness toward variations of such characteristics. On point clouds from the human subject, the proposed method successfully reconstructed all patient surfaces, filling regions where raw point coordinate readings were missing. Within two comparable regions of interest in the chest area, similar mean curvature distributions were acquired from both their reconstructed surface and CT surface, with mean and standard deviation of (μ{sub recon} = − 2.7 × 10{sup −3} mm{sup −1}, σ{sub recon} = 7.0 × 10{sup −3} mm{sup −1}) and (μ{sub CT} = − 2.5 × 10{sup −3} mm{sup −1}, σ{sub CT} = 5.3 × 10{sup −3} mm{sup −1}), respectively. The agreement of local geometry properties between the reconstructed surfaces and the CT surface demonstrated the ability of the proposed method in faithfully representing the underlying patient surface. Conclusions: The authors have integrated and developed an accurate level-set based continuous surface reconstruction method on point clouds acquired by a 3D surface photogrammetry system. The proposed method has generated a continuous representation of the underlying phantom and patient surfaces with good robustness against noise and missing measurements. It serves as an important first step for further development of motion tracking methods during radiotherapy.
60 APPLIED LIFE SCIENCES↗
CPU Isosurface Ray Tracing of Adaptive Mesh Refinement Data
Adaptive mesh refinement (AMR) is a key technology for large-scale simulations that allows for adaptively changing the simulation mesh resolution, resulting in significant computational and storage savings. However, visualizing such AMR data poses a significant challenge due to the difficulties introduced by the hierarchical representation when reconstructing continuous field values. Here, we detail a comprehensive solution for interactive isosurface rendering of block-structured AMR data. We contribute a novel reconstruction strategy—the octant method—which is continuous, adaptive and simple to implement. Furthermore, we present a generally applicable hybrid implicit isosurface ray-tracing method, which provides better rendering quality and performance than the built-in sampling-based approach in OSPRay. Finally, we integrate our octant method and hybrid isosurface geometry into OSPRay as a module, providing the ability to create high-quality interactive visualizations combining volume and isosurface representations of BS-AMR data. Lastly, we evaluate the rendering performance, memory consumption and quality of our method on two gigascale block-structured AMR datasets.
97 MATHEMATICS AND COMPUTING↗
TH-AB-202-08: A Robust Real-Time Surface Reconstruction Method On Point Clouds Captured From a 3D Surface Photogrammetry System
Purpose: Surface photogrammetry (e.g. VisionRT, C-Rad) provides a noninvasive way to obtain high-frequency measurement for patient motion monitoring in radiotherapy. This work aims to develop a real-time surface reconstruction method on the acquired point clouds, whose acquisitions are subject to noise and missing measurements. In contrast to existing surface reconstruction methods that are usually computationally expensive, the proposed method reconstructs continuous surfaces with comparable accuracy in real-time. Methods: The key idea in our method is to solve and propagate a sparse linear relationship from the point cloud (measurement) manifold to the surface (reconstruction) manifold, taking advantage of the similarity in local geometric topology in both manifolds. With consistent point cloud acquisition, we propose a sparse regression (SR) model to directly approximate the target point cloud as a sparse linear combination from the training set, building the point correspondences by the iterative closest point (ICP) method. To accommodate changing noise levels and/or presence of inconsistent occlusions, we further propose a modified sparse regression (MSR) model to account for the large and sparse error built by ICP, with a Laplacian prior. We evaluated our method on both clinical acquired point clouds under consistent conditions and simulated point clouds with inconsistent occlusions. The reconstruction accuracy was evaluated w.r.t. root-mean-squared-error, by comparing the reconstructed surfaces against those from the variational reconstruction method. Results: On clinical point clouds, both the SR and MSR models achieved sub-millimeter accuracy, with mean reconstruction time reduced from 82.23 seconds to 0.52 seconds and 0.94 seconds, respectively. On simulated point cloud with inconsistent occlusions, the MSR model has demonstrated its advantage in achieving consistent performance despite the introduced occlusions. Conclusion: We have developed a real-time and robust surface reconstruction method on point clouds acquired by photogrammetry systems. It serves an important enabling step for real-time motion tracking in radiotherapy. This work is supported in part by NIH grant R01 CA169102-02.
60 APPLIED LIFE SCIENCES↗
Real-Time event reconstruction for Nuclear Physics Experiments using Artificial Intelligence
Charged track reconstruction is a critical task in nuclear physics experiments, enabling the identification and analysis of particles produced in high-energy collisions. Machine learning (ML) has emerged as a powerful tool for this purpose, addressing the challenges posed by complex detector geometries, high event multiplicities, and noisy data. Traditional methods rely on pattern recognition algorithms like the Kalman filter, but ML techniques, such as neural networks, graph neural networks (GNNs), and recurrent neural networks (RNNs), offer improved accuracy and scalability. By learning from simulated and real detector data, ML models can identify and classify tracks, predict trajectories, and handle ambiguities caused by overlapping or missing hits. Moreover, ML-based approaches can process data in near-real-time, enhancing the efficiency of experiments at large-scale facilities like the Large Hadron Collider (LHC) and Jefferson Lab (JLAB). As detector technologies and computational resources evolve, ML-driven charged track reconstruction continues to push the boundaries of precision and discovery in nuclear physics. In these proceedings, we highlight advancements in charged track identification leveraging Artificial Intelligence within the CLAS12 detector, achieving a notable enhancement in experimental statistics compared to traditional methods. Additionally, we showcase real-time event reconstruction capabilities, including the inference of charged particle properties, such as momentum, direction, and species identification, at speeds matching data acquisition rates. These innovations enable the extraction of physics observables directly from the experiment in real-time.
Gavalian, Gagik (ORCID:0000000267385457)↗
Numerical Modeling of Atmospheric Rime Ice Accretion on an Airfoil Using an Eulerian Approach
This paper presents an approach to numerically simulate the inherently unsteady rime ice accretion problem on a two-dimensional airfoil and elucidate the associated variations under different icing conditions. The airflow field and the water impingement on the airfoil are obtained based on an Eulerian two-phase model. A dynamic mesh strategy is employed to unsteadily account for the changes in the ice profile and its impact on the air and droplet flow by continuously reconstructing the computational grid at each time-step through smoothing and layering mechanisms. All main icing modules including the airflow field, droplet trajectory, icing thickness profile, and mesh management are fully coupled within the same computational framework without resorting to any external tools. Classical icing theory is employed to model the rime ice roughness, and it is assumed that the ice accretes in a direction normal to the airfoil surface. The governing Reynolds-averaged Navier–Stokes (RANS) conservation equations along with the energy and continuity equations are solved to produce the velocity and temperature fields. A convective film heat transfer coefficient is computed based on the surface heat flux and a recovery temperature which takes into account the dissipative heat release in the boundary layer in the vicinity of the airfoil surface. With the implemented strategy and calculating the convective heat transfer coefficient, the water film thickness is also calculated along with the ice shape. Furthermore, the model is validated by comparing the local collection efficiency distribution and ice shape with experimental data, and the results show that the implemented approach provides acceptable predictions of ice accretion profiles and rates.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
Multiple independent records of local glacier variability on Nuussuaq, West Greenland, during the Holocene
The sensitivity of mountain glaciers to small-scale climate fluctuations makes records of their past extent among the best proxies for identifying spatio-temporal climate variability. Here we build on the few existing Holocene records of local glacier change in Greenland by using three independent geochronological methods—proglacial lake sediment analysis, cosmogenic 10 Be surface-exposure dating, and in situ 14 C burial modeling—to reconstruct continuous records of Holocene glacier variability on Nuussuaq, West Greenland. 10 Be ages of perched boulders indicate eastern Nuussuaq was deglaciated between ~11.0 and 10.5 ka. Radiocarbon-dated sediments from two lakes on Nuussuaq contain mineral-rich layers between ~9.6 and 9.0 and ~8.7–8.0 cal ka BP that may be correlative with nearby ice sheet moraines deposited in the early Holocene. Multiple proxies for glacier size indicate frequent, high-amplitude glacier fluctuations superimposed on net glacier growth during the late Holocene, with significant ice expansion phases at ~3.7 ka, 2.8 ka, and throughout the past ~2 ka. Mean 10 Be ages from five nested moraine crests confirm that local glacier extents on Nuussuaq culminated during both the Little Ice Age [~1470 C.E. (n=3) and 1750 C.E. (n=3)] and the preceding centuries (~520–1320 C.E.; n=11). Results reveal that local glaciers on Nuussuaq episodically advanced and retreated at centennial timescales throughout the Holocene, most likely in response to regional climate changes in West Greenland superimposed on the progressive insolation-driven cooling trend in the Northern Hemisphere. Our new 10 Be moraine chronologies coupled with other glacier-size proxies corroborate an emerging pattern of significant summer cooling and glacier expansion in the centuries prior to the Little Ice Age in the Arctic.
54 ENVIRONMENTAL SCIENCES↗
Quantitative probing of the fast particle motion during the solidification of battery electrodes
Synchrotron x-ray tomography represents one of the most powerful imaging techniques that can visualize the three-dimensional microstructural complexity with a high spatial resolution. To date, however, the temporal resolution of this technique has been handicapped by the intrinsic limitations in the data acquisition speed and the lack of a suitable data processing method. The present study aims to fundamentally resolve this challenge by implementing a data analysis method that is based on differentiating discrete tomography reconstruction from a continuous scan reconstruction. With this methodological development, we capture and quantify the particle motion during the solidification of battery electrodes and reveal the statistics of the dynamically evolving motion in the drying process, which has been challenging to resolve. We discover that the particle motion exhibits a strong dependence on its geometric location within a drying electrode. Our results also imply that the final electrode quality can be controlled by balancing the solvent evaporation rate and the particle mobility in the region close to the drying surface.
25 ENERGY STORAGE↗
Asymptotic, multigroup flux reconstruction and consistent discontinuity factors
Recent theoretical work has led to an asymptotically derived expression for reconstructing the neutron flux from lattice functions and multigroup diffusion solutions. The leading-order asymptotic term is the standard expression for flux reconstruction, i.e., it is the product of a shape function, obtained through a lattice calculation, and the multigroup diffusion solution. The first-order asymptotic correction term is significant only where the gradient of the diffusion solution is not small. Inclusion of this first-order correction term can significantly improve the accuracy of the reconstructed flux. One may define discontinuity factors (DFs) to make certain angular moments of the reconstructed flux continuous across interfaces between assemblies in 1-D. Indeed, the standard assembly discontinuity factors make the zeroth moment (scalar flux) of the reconstructed flux continuous. The inclusion of the correction term in the flux reconstruction provides an additional degree of freedom that can be used to make two angular moments of the reconstructed flux continuous across interfaces by using current DFs in addition to flux DFs. Thus, numerical results demonstrate that using flux and current DFs together can be more accurate than using only flux DFs, and that making the second angular moment continuous can be more accurate than making the zeroth moment continuous.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
A hybrid reconstructed discontinuous Galerkin and continuous Galerkin finite element method for incompressible flows on unstructured grids
Not Available
Pandare, Aditya K. (ORCID:0000000301637321)↗
AIVT: Inference of turbulent thermal convection from measured 3D velocity data by physics-informed Kolmogorov-Arnold networks
We propose the artificial intelligence velocimetry-thermometry (AIVT) method to reconstruct a continuous and differentiable representation of the temperature and velocity in turbulent convection from measured three-dimensional (3D) velocity data. AIVT is based on physics-informed Kolmogorov-Arnold networks and trained by optimizing a loss function that minimizes residuals of the velocity data, boundary conditions, and governing equations. We apply AIVT to a set of simultaneously measured 3D temperature and velocity data of Rayleigh-Bénard convection, obtained by combining particle image thermometry and Lagrangian particle tracking. This enables us to directly compare machine learning results to true volumetric, simultaneous temperature and velocity measurements. We demonstrate that AIVT can reconstruct and infer continuous, instantaneous velocity and temperature fields and their gradients from sparse experimental data at a high resolution, providing an additional approach for understanding thermal turbulence.
Science & Technology - Other Topics↗
Maximum a posteriori estimation of crystallographic phases in X-ray diffraction tomography
A maximum a posteriori approach is proposed for X-ray diffraction tomography for reconstructing three-dimensional spatial distribution of crystallographic phases and orientations of polycrystalline materials. The approach maximizes the a posteriori density which includes a Poisson log-likelihood and an a priori term that reinforces expected solution properties such as smoothness or local continuity. The reconstruction method is validated with experimental data acquired from a section of the spinous process of a porcine vertebra collected at the 1-ID-C beamline of the Advanced Photon Source, at Argonne National Laboratory. The reconstruction results show significant improvement in the reduction of aliasing and streaking artefacts, and improved robustness to noise and undersampling compared to conventional analytical inversion approaches. The approach has the potential to reduce data acquisition times, and significantly improve beamtime efficiency.
Biomineralization↗
A multidimensional approach to quantum state tomography of photoelectron wavepackets
There is a growing interest in reconstructing the density matrix of photoelectron wavepackets, in particular in complex systems where decoherence can be introduced either by a partial measurement of the system or through coupling with a stochastic environment. To this end, several methods to reconstruct the density matrix, quantum state tomography protocols, have been developed and tested on photoelectrons ejected from noble gases following absorption of extreme ultraviolet (XUV) photons from attosecond pulses. It remains a challenge to obtain model-free, single scan protocols that can reconstruct the density matrix with high fidelities. Current methods require extensive measurements or involve complex fitting of the signal. Efficient single-scan reconstructions would be of great help to increase the number of systems that can be studied. We propose a new and more efficient protocol that is able to reconstruct the continuous variable density matrix of a photoelectron in a single time delay scan. It is based on measuring the coherences of a photoelectron created by absorption of an XUV pulse using a broadband infrared (IR) probe that is scanned in time and a narrowband IR reference that is temporally fixed to the XUV pulse. We illustrate its performance for a Fano resonance in He as well as mixed states in Ar arising from spin-orbit splitting. We show that the protocol results in excellent fidelities and near-perfect estimation of the purity.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Adaptive Multilinear Tensor Product Wavelets
Many foundational visualization techniques including isosurfacing, direct volume rendering and texture mapping rely on piecewise multilinear interpolation over the cells of a mesh. However, there has not been much focus within the visualization community on techniques that efficiently generate and encode globally continuous functions defined by the union of multilinear cells. Wavelets provide a rich context for analyzing and processing complicated datasets. In this paper, we exploit adaptive regular refinement as a means of representing and evaluating functions described by a subset of their nonzero wavelet coefficients. We analyze the dependencies involved in the wavelet transform and describe how to generate and represent the coarsest adaptive mesh with nodal function values such that the inverse wavelet transform is exactly reproduced via simple interpolation (subdivision) over the mesh elements. This allows for an adaptive, sparse representation of the function with on-demand evaluation at any point in the domain. In conclusion, we focus on the popular wavelets formed by tensor products of linear B-splines, resulting in an adaptive, nonconforming but crack-free quadtree (2D) or octree (3D) mesh that allows reproducing globally continuous functions via multilinear interpolation over its cells.
97 MATHEMATICS AND COMPUTING↗
Velocity reconstruction in the era of DESI and Rubin/LSST. I. Exploring spectroscopic, photometric, and hybrid samples
Peculiar velocities of galaxies and halos can be reconstructed from their spatial distribution alone. This technique is analogous to the baryon acoustic oscillations reconstruction, using the continuity equation to connect density and velocity fields. The resulting reconstructed velocities can be used to measure imprints of galaxy velocities on the cosmic microwave background like the kinematic Sunyaev-Zel’dovich effect or the moving lens effect. As the precision of these measurements increases, characterizing the performance of the velocity reconstruction becomes crucial to allow unbiased and statistically optimal inference. In this paper, we quantify the relevant performance metrics: the variance of the reconstructed velocities and their correlation coefficient with the true velocities. We show that the relevant velocities to reconstruct for kSZ and moving lens are actually the halo—rather than galaxy—velocities. We quantify the impact of redshift-space distortions, photometric redshift errors, satellite galaxy fraction, incorrect cosmological parameter assumptions and smoothing scale on the reconstruction performance. Here, we also investigate hybrid reconstruction methods, where velocities inferred from spectroscopic samples are evaluated at the positions of denser photometric samples. We find that using exclusively the photometric sample is better than performing a hybrid analysis. The 2 Gpc/ℎ length simulations from abacussummit with realistic galaxy samples for DESI and Rubin LSST allow us to perform this analysis in a controlled setting. In the companion paper [B. Hadzhiyska, S. Ferraro, B. Ried Guachalla, and E. Schaan, companion paper, Phys. Rev. D 109, 103534 (2024).], we further include the effects of evolution along the light cone and give realistic performance estimates for DESI luminous red galaxies, emission line galaxies, and Rubin LSST-like samples.
79 ASTRONOMY AND ASTROPHYSICS↗
Determination of activity of {sup 51}Cr on gamma radiation measurements
A method of determining the activity of intensive distributed -sources on the measurement of the continuous spectrum of radiation, for example the internal bremsstrahlung, is developed. The recurrent formula for reconstructing of a continuous spectrum, registered in a Ge detector, at distorting it in the detector. The method of precise measurements of the spectrum of {sup 51}Cr internal bremsstrahlung using two point sources of low activity is described.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Entanglement wedge reconstruction and the information paradox
When absorbing boundary conditions are used to evaporate a black hole in AdS/CFT, we show that there is a phase transition in the location of the quantum Ryu-Takayanagi surface, at precisely the Page time. The new RT surface lies slightly inside the event horizon, at an infalling time approximately the scrambling time β/2 π logS BH into the past. We can immediately derive the Page curve, using the Ryu-Takayanagi formula, and the Hayden-Preskill decoding criterion, using entanglement wedge reconstruction. Because part of the interior is now encoded in the early Hawking radiation, the decreasing entanglement entropy of the black hole is exactly consistent with the semiclassical bulk entanglement of the late-time Hawking modes, despite the absence of a firewall.By studying the entanglement wedge of highly mixed states, we can understand the state dependence of the interior reconstructions. A crucial role is played by the existence of tiny, non-perturbative errors in entanglement wedge reconstruction. Directly after the Page time, interior operators can only be reconstructed from the Hawking radiation if the initial state of the black hole is known. As the black hole continues to evaporate, reconstructions become possible that simultaneously work for a large class of initial states. Using similar techniques, we generalise Hayden-Preskill to show how the amount of Hawking radiation required to reconstruct a large diary, thrown into the black hole, depends on both the energy and the entropy of the diary. Finally we argue that, before the evaporation begins, a single, state-independent interior reconstruction exists for any code space of microstates with entropy strictly less than the Bekenstein-Hawking entropy, and show that this is sufficient state dependence to avoid the AMPSS typical-state firewall paradox.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
SU-E-T-18: A Comparison of Planning Techniques for Bilateral Reconstructed Chest Wall Patients Undergoing Whole Breast Irradiation
Purpose: As we continuously see more bilateral reconstructed chest wall cases, new challenges are being presented to deliver left-sided breast irradiation. We herein compare three Deep Inspiration Breath Hold (DIBH) planning techniques (tangents, VMAT, and IMRT) and two free breathing techniques (VMAT and IMRT). Methods: Three left-sided chest wall patients with bilateral implants were studied. Tangents, VMAT, and IMRT plans were created for DIBH scans. VMAT and IMRT plans were created for free breathing scans. All plans were normalized so that 95% of the prescription dose was delivered to 95% of the planning target volume (PTV). The maximum point dose was constrained to less than 120% of the prescription dose. Since the success of DIBH delivery largely depends on patient’s ability to perform consistent breath hold during beam on time, smaller number of Monitor Units (MU) is in general desired. For each patient, the following information was collected to compare the planning techniques: heart mean dose, left and right lung V20 Gy, contra-lateral (right) breast mean dose, cord max dose, and MU. Results: The average heart mean dose over all patients are 1561, 692, 985, 1245, and 1121 cGy, for DIBH tangents, VMAT, IMRT, free breathing VMAT and IMRT, respectively. For left lung V20 are 60%, 28%, 26%, 30%, and 29%. For contra-lateral breast mean dose are 244, 687, 616, 783, 438 cGy. MU are 253, 853, 2048, 1035, and 1874 MUs. Conclusion: In the setting of bilateral chest wall reconstruction, opposed tangent beams cannot consistently achieve desired heart and left lung sparing. DIBH consistently achieves better healthy tissue sparing. VMAT appears to be preferential to IMRT for planning and delivering radiation to patients with bilaterally reconstructed chest walls being treated with DIBH.
60 APPLIED LIFE SCIENCES↗