Search NASA⌕ Search

SEARCH · Search NASA

Results for “Estimation method”

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.

At least 1,495 records · Page 83

Anomaly Detection in Power System State Estimation: Review and New Directions

Foundational and state-of-the-art anomaly-detection methods through power system state estimation are reviewed. Traditional components for bad data detection, such as chi-square testing, residual-based methods, and hypothesis testing, are discussed to explain the motivations for recent anomaly-detection methods given the increasing complexity of power grids, energy management systems, and cyber-threats. In particular, state estimation anomaly detection based on data-driven quickest-change detection and artificial intelligence are discussed, and directions for research are suggested with particular emphasis on considerations of the future smart grid.

42 ENGINEERING↗

Overcoming Communications Outages in Inverter Downtime Analysis

This paper presents two methods of detecting inverter downtime and estimating lost production from downtime events using timeseries system production measurements. The methods focus on distinguishing communications interruptions from true production outages and are successful in most cases. To enable fleet-scale analysis of inverter availability, the methods are designed to be semi-autonomous. The first method uses only inverter and system meter AC power measurements while the second uses cumulative production and expected energy data to classify outages and estimate lost production. The methods are demonstrated using real field data and the results are discussed.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Enhancing statistical moment calculations for stochastic Galerkin solutions with Monte Carlo techniques

Here, we provide a method for enhancing stochastic Galerkin moment calculations to the linear elliptic equation with random diffusivity using an ensemble of Monte Carlo solutions. This hybrid approach combines the accuracy of low-order stochastic Galerkin and the computational efficiency of Monte Carlo methods to provide statistical moment estimates which are significantly more accurate than performing each method individually. The hybrid approach involves computing a low-order stochastic Galerkin solution, after which Monte Carlo techniques are used to estimate the residual. We show that the combined stochastic Galerkin solution and residual is superior in both time and accuracy for a one-dimensional test problem and a more computational intensive two-dimensional linear elliptic problem for both the mean and variance quantities.

97 MATHEMATICS AND COMPUTING↗

Optimization performance, fidelity, and cost: SIAM VQE

This dataset contains files storing results from classically-simulated quantum subroutines within a dynamical mean-field theory workflow, and jupyter notebooks processing the data in these files to generate plots. The files store: (1) Results from variational quantum eigensolver (VQE) simulations searching for optimal parameters allowing parametrized quantum circuits to prepare approximations to ground states of different Anderson impurity models (AIMs) (2) Results from simulations of a quantum Lanczos algorithm (QLA) estimating the Lanczos coefficients defining the continued-fraction representation of an (AIM) Green’s function Description: Any file named vqe_gs_results* stores approximations to the ground state and energy of a given AIM estimated using three different methods: (1) Numerical diagonalization (2) Ideal VQE simulation (3) VQE simulation with sampling noise For each VQE simulations metadata about the optimization (optimization results plus number of quantum circuits that would have been executed on real hardware) is also stored. Any file named qla_dos_results* estimations for the Lanczos coefficients defining the Green’s function of an AIM. The stored estimations are achieved using different methods: (1) Numerical Lanczos algorithm from initial states obtained from numerical diagonalization (2) Simulated quantum Lanczos algorithm from initial states prepared from parametrized quantum circuits yielded by corresponding ideal and noisy VQE subroutines. The dataset is used and described in M. Karabin et al., "Quantum solver for single-impurity Anderson models with particle-hole symmetry", Phys. Rev. Research 8, 033066 (2026). DOI: https://doi.org/10.1103/7ys3-tl4l

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Large Scale Study of Ligand–Protein Relative Binding Free Energy Calculations: Actionable Predictions from Statistically Robust Protocols

The accurate and reliable prediction of protein–ligand binding affinities can play a central role in the drug discovery process as well as in personalized medicine. Of considerable importance during lead optimization are the alchemical free energy methods that furnish an estimation of relative binding free energies (RBFE) of similar molecules. Recent advances in these methods have increased their speed, accuracy, and precision. This is evident from the increasing number of retrospective as well as prospective studies employing them. However, such methods still have limited applicability in real-world scenarios due to a number of important yet unresolved issues. Here, we report the findings from a large data set comprising over 500 ligand transformations spanning over 300 ligands binding to a diverse set of 14 different protein targets which furnish statistically robust results on the accuracy, precision, and reproducibility of RBFE calculations. We use ensemble-based methods which are the only way to provide reliable uncertainty quantification given that the underlying molecular dynamics is chaotic. These are implemented using TIES (Thermodynamic Integration with Enhanced Sampling). Results achieve chemical accuracy in all cases. Ensemble simulations also furnish information on the statistical distributions of the free energy calculations which exhibit non-normal behavior. We find that the “enhanced sampling” method known as replica exchange with solute tempering degrades RBFE predictions. We also report definitively on numerous associated alchemical factors including the choice of ligand charge method, flexibility in ligand structure, and the size of the alchemical region including the number of atoms involved in transforming one ligand into another. Our findings provide a key set of recommendations that should be adopted for the reliable application of RBFE methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reconstruction of neuronal input through modeling single-neuron dynamics and computations

Mathematical models provide a mathematical description of neuron activity, which can better understand and quantify neural computations and corresponding biophysical mechanisms evoked by stimulus. In this paper, based on the output spike train evoked by the acupuncture mechanical stimulus, we present two different levels of models to describe the input-output system to achieve the reconstruction of neuronal input. The reconstruction process is divided into two steps: First, considering the neuronal spiking event as a Gamma stochastic process. The scale parameter and the shape parameter of Gamma process are, respectively, defined as two spiking characteristics, which are estimated by a state-space method. Then, leaky integrate-and-fire (LIF) model is used to mimic the response system and the estimated spiking characteristics are transformed into two temporal input parameters of LIF model, through two conversion formulas. We test this reconstruction method by three different groups of simulation data. All three groups of estimates reconstruct input parameters with fairly high accuracy. We then use this reconstruction method to estimate the non-measurable acupuncture input parameters. Results show that under three different frequencies of acupuncture stimulus conditions, estimated input parameters have an obvious difference. The higher the frequency of the acupuncture stimulus is, the higher the accuracy of reconstruction is.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Milky Way mass with K giants and BHB stars using LAMOST, SDSS/SEGUE, and Gaia : 3D spherical Jeans equation and tracer mass estimator

ABSTRACT We measure the enclosed Milky Way mass profile to Galactocentric distances of ∼70 and ∼50 kpc using the smooth, diffuse stellar halo samples of Bird et al. The samples are Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and Sloan Digital Sky Survey/Sloan Extension for Galactic Understanding and Exploration (SDSS/SEGUE) K giants (KG) and SDSS/SEGUE blue horizontal branch (BHB) stars with accurate metallicities. The 3D kinematics are available through LAMOST and SDSS/SEGUE distances and radial velocities and Gaia DR2 proper motions. Two methods are used to estimate the enclosed mass: 3D spherical Jeans equation and Evans et al. tracer mass estimator (TME). We remove substructure via the Xue et al. method based on integrals of motion. We evaluate the uncertainties on our estimates due to random sampling noise, systematic distance errors, the adopted density profile, and non-virialization and non-spherical effects of the halo. The tracer density profile remains a limiting systematic in our mass estimates, although within these limits we find reasonable agreement across the different samples and the methods applied. Out to ∼70 and ∼50 kpc, the Jeans method yields total enclosed masses of 4.3 ± 0.95 (random) ±0.6 (systematic) × 1011 M⊙ and 4.1 ± 1.2 (random) ±0.6 (systematic) × 1011 M⊙ for the KG and BHB stars, respectively. For the KG and BHB samples, we find a dark matter virial mass of $M_{200}=0.55^{+0.15}_{-0.11}$ (random) ±0.083 (systematic) × 1012 M⊙ and $M_{200}=1.00^{+0.67}_{-0.33}$ (random) ±0.15 (systematic) × 1012 M⊙, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

Accelerated Sensitivity Analysis in High-Dimensional Stochastic Reaction Networks

Existing sensitivity analysis approaches are not able to handle efficiently stochastic reaction networks with a large number of parameters and species, which are typical in the modeling and simulation of complex biochemical phenomena. In this paper, a two-step strategy for parametric sensitivity analysis for such systems is proposed, exploiting advantages and synergies between two recently proposed sensitivity analysis methodologies for stochastic dynamics. The first method performs sensitivity analysis of the stochastic dynamics by means of the Fisher Information Matrix on the underlying distribution of the trajectories; the second method is a reduced-variance, finite-difference, gradient-type sensitivity approach relying on stochastic coupling techniques for variance reduction. Here in this paper we demonstrate that these two methods can be combined and deployed together by means of a new sensitivity bound which incorporates the variance of the quantity of interest as well as the Fisher Information Matrix estimated from the first method. The first step of the proposed strategy labels sensitivities using the bound and screens out the insensitive parameters in a controlled manner. In the second step of the proposed strategy, a finite-difference method is applied only for the sensitivity estimation of the (potentially) sensitive parameters that have not been screened out in the first step. Results on an epidermal growth factor network with fifty parameters and on a protein homeostasis with eighty parameters demonstrate that the proposed strategy is able to quickly discover and discard the insensitive parameters and in the remaining potentially sensitive parameters it accurately estimates the sensitivities. The new sensitivity strategy can be several times faster than current state-of-the-art approaches that test all parameters, especially in “sloppy” systems. In particular, the computational acceleration is quantified by the ratio between the total number of parameters over the number of the sensitive parameters.

97 MATHEMATICS AND COMPUTING↗

Estimation of constituent properties of concrete materials with an artificial neural network based method

Multi-scale models are developed for heterogeneous concrete materials to estimate their macroscopic mechanical properties in terms of micro-structural data. One crucial challenge of those models is the identification of local properties of constituent phases. In this paper, we present an efficient method based on Artificial Neural Networks (ANN). Typical concrete materials are taken as example. A macroscopic analytical strength criterion is established from three steps of nonlinear homogenization procedure. The macroscopic strength of materials is determined as a function of the frictional coefficient and cohesion of solid cement particles at nanometer scale, intra-particle pores, inter-particle pores and aggregates (inclusions). The objective is to identify the nanoscopic frictional coefficient and cohesion of cement particle from measured macroscopic values of uniaxial compression and tensile strengths. For this purpose, a numerical method based on the ANN is developed. With the analytical macroscopic strength criterion, sensitivity studies are first realized to identify the most important micro-structural parameters influencing the macroscopic strength of concrete. A simplified analytical macroscopic strength criterion is then proposed. A large dataset is further constructed through the inversion of the analytical strength criterion by using the aggregates volume fraction, porosity, macroscopic uniaxial tensile and compressive strengths as input variables and the frictional coefficient and cohesion of cement particles as output unknowns. An ANN model containing four hidden layers and 100 neurons in each layer is constructed and trained by using this dataset. Various types of validation of the ANN model are performed. It is found that the proposed ANN based model can effectively predict the frictional coefficient and cohesion of porous cement paste at the microscopic scale with a very good accuracy.

36 MATERIALS SCIENCE↗

Time-temperature-superposition analysis of diverse datasets by the minimum-arclength method: long-term prediction with uncertainty margins

In a recent publication, we carried out an extensive analysis of an unsupervised method of determining optimum shift factors in time-temperature-superposition of accelerated-aging data that involves minimizing the vertical arclength to obtain the master curve. For synthetic Arrhenius data with a variety of noise distributions, the work showed that, in conjunction with bootstrap-resampling, the method can produce reliable estimates of the mean activation energy along with uncertainty quantification. Here, we apply the above method to six different datasets taken from the published literature and demonstrate accurate prediction of mean activation energy from the data as-is without the need for any pre-processing or fitting. We also compare uncertainty margins computed by second-order bootstrap with that by linear regression theory and show that the former appears to provide consistent margins in the presence of common noise types in real data, including intra-isotherm measurement-errors, sample-to-sample variations, and intrinsic deviation from perfect Arrhenius behavior.

36 MATERIALS SCIENCE↗

Sources of Propane Consumed in California

Project Scope: The objective of this study is to specify the sources of propane consumed in California. It answers the questions, where does the propane used in California come from and how was it produced? The results of this study provide comprehensive, transparent, and verifiable estimates, based on the 2018 market. The information provided in this report is suitable for use to assess the life cycle carbon intensity of propane used as a transportation fuel in California. As the 2009 Low Carbon Fuel Standard (LCFS) aims to reduce California’s greenhouse gas (GHG) emissions and other smog-forming and toxic air pollutants, the appropriate designation of carbon intensity for propane as a transportation fuel is important for evaluating propane’s potential to contribute to GHG goals and understandings in the context of various actions. This study focuses on estimating the shares of total propane consumed in the state of California produced from petroleum refineries, natural gas plants, and bituminous sands sources inside California and elsewhere. Results: An estimated 590 million gallons of propane were consumed in California in 2018, of which, 59.5% originated from refinery production and 40.5% originated from natural gas plants. The majority of this was sourced from refinery production in California, 334 million gallons. Most of the propane imported to California for consumption was sourced from natural gas plants, 113 million gallons, with over half of the imported volume sourced from Canada. The volume sourced from bituminous sand upgrader and fractionator operations was negligible. Details from this analysis are presented in the table below which provides an overview of the propane flows estimated in this study by region and production method. The shares and volumes presented here represent a snapshot for 2018. A significant increase in propane demand, such as could be caused by increased use of propane as a transportation fuel in the state, would affect California’s propane production, imports, and exports. The method and data sources used for the estimates provided in this report also provide the framework which could be used for future updates. Key Method Considerations: The values presented here are based on a two-step approach where the first step was to determine the flows of propane into and out of California from different regions and the second step was to estimate the propane production methods in each region. A volume balance approach is used as the primary method for tracking the volume of propane in and out of California as propane production and import volumes are available by Petroleum Administration of Defense District (PADD) from EIA and neither inter-PADD propane transfers nor state-specific non-prime supplier consumption are available from a public data source. The volume balance performed for this study covered PADD 5 (the West Coast), which includes Arizona, California, Nevada, Oregon, and Washington. The volume balance used all available public datasets to determine propane production, imports, exports, and consumption. Volumes unaccounted for by these datasets were estimated using the resulting volume balance by assuming market equilibrium. Consumption within each state in PADD 5 was estimated based on known import, export, and production volumes and this amount was used to develop the volume balance. EIA only tracks consumption at the state level by prime supplier sales. The volume balance approach provides the basis to correct for additional propane consumed in-states where propane is transferred to California. To determine the California propane sources and trade in 2018. volume of propane consumed in California, the volume balance approach is again used where it was estimated all imported volumes not specifically flagged for re-export were consumed, and the remaining consumption was produced in-state. The California Energy Commission (CEC) provided the total volume of propane imported and exported from California in 2018; this volume data set along with commodity tracking from the Canada Energy Regulator (CER) and the International Trade Commission (ITC) which tracks port of entry and final destination was used to determine where propane originated from and where it was ultimately consumed. For example, the CER tracks propane leaving Canada and entering each state within the U.S. Imported propane from Canada to California – marked for California – is assumed to be consumed in California. When no further data were available, import volumes were assumed to be consumed in California without pass-through (i.e., no propane imported to California was directly sold and exported). In most cases, the production method for each propane source region was applied to the volume of propane transferred to California. In other words, the shares of propane sourced from natural gas and refineries for each production region was assigned to California imports based on their contribution to the total volume flows into California to determine the production method for propane consumed in-state. For volumes imported into California from PADD 4, Washington State, Canada, and the rest of the world (Argentina, Chile, Norway, Peru, South Korea, and Trinidad and Tobago), the volumes sourced from petroleum refineries and natural gas plants reflect the either production ratio for the region or, in cases where the sources specific to the amounts exported to California could be determined, the sources specific to the volumes transferred to California.

03 NATURAL GAS↗

SU-F-I-80: Correction for Bias in a Channelized Hotelling Model Observer Caused by Temporally Variable Non-Stationary Noise

Purpose: Application of a channelized Hotelling model observer (CHO) over a wide range of x-ray angiography detector target dose (DTD) levels demonstrated substantial bias for conditions yielding low detectability indices (d’), including low DTD and small test objects. The purpose of this work was to develop theory and methods to correct this bias. Methods: A hypothesis was developed wherein the measured detectability index (d’b) for a known test object is positively biased by temporally variable non-stationary noise in the images. Hotelling’s T2 test statistic provided the foundation for a mathematical theory which accounts for independent contributions to the measured d’b value from both the test object (d’o) and non-stationary noise (d’ns). Experimental methods were developed to directly estimate d’o by determining d’ns and subtracting it from d’b, in accordance with the theory. Specifically, d’ns was determined from two sets of images from which the traditional test object was withheld. This method was applied to angiography images with DTD levels in the range 0 to 240 nGy and for disk-shaped iodine-based contrast targets with diameters 0.5 to 4.0 mm. Results: Bias in d’ was evidenced by d’b values which exceeded values expected from a quantum limited imaging system and decreasing object size and DTD. d’ns increased with decreasing DTD, reaching a maximum of 2.6 for DTD = 0. Bias-corrected d’o estimates demonstrated sub-quantum limited performance of the x-ray angiography for low DTD. Findings demonstrated that the source of non-stationary noise was detector electronic readout noise. Conclusion: Theory and methods to estimate and correct bias in CHO measurements from temporally variable non-stationary noise were presented. The temporal non-stationary noise was shown to be due to electronic readout noise. This method facilitates accurate estimates of d’ values over a large range of object size and detector target dose.

60 APPLIED LIFE SCIENCES↗

Bayesian forecasting and uncertainty quantifying of stream flows using Metropolis–Hastings Markov Chain Monte Carlo algorithm

This paper presents a Bayesian approach using Metropolis-Hastings Markov Chain Monte Carlo algorithm and applies this method for daily river flow rate forecast and uncertainty quantification for Zhujiachuan River using data collected from Qiaotoubao Gage Station and other 13 gage stations in Zhujiachuan watershed in China. The proposed method is also compared with the conventional maximum likelihood estimation (MLE) for parameter estimation and quantification of associated uncertainties. While the Bayesian method performs similarly in estimating the mean value of daily flow rate, it performs over the conventional MLE method on uncertainty quantification, providing relatively narrower reliable interval than the MLE confidence interval and thus more precise estimation by using the related information from regional gage stations. As a result, the Bayesian MCMC method might be more favorable in the uncertainty analysis and risk management.

58 GEOSCIENCES↗

Jet formation in cerium metal to examine material strength

Examining the evolution of material properties at extreme conditions advances our understanding of numerous high-pressure phenomena from natural events like meteorite impacts to general solid mechanics and fluid flow behavior. Recent advances in synchrotron diagnostics coupled with dynamic compression platforms have introduced new possibilities for examining in-situ, spatially resolved material response with nanosecond time resolution. In this work, we examined jet formation from a Richtmyer-Meshkov instability in cerium initially shocked into a transient, high-pressure phase, and then released to a low-pressure, higher-temperature state. Cerium's rich phase diagram allows us to study the yield stress following a shock induced solid-solid phase transition. X-ray imaging was used to obtain images of jet formation and evolution with 2–3 μm spatial resolution. From these images, an analytic method was used to estimate the post-shock yield stress, and these results were compared to continuum calculations that incorporated an experimentally validated equation-of-state (EOS) for cerium coupled with a deviatoric strength model. Reasonable agreement was observed between the calculations and the data illustrating the sensitivity of jet formation on the yield stress values. The data and analysis shown here provide insight into material strength during dynamic loading which is expected to aid in the development of strength aware multi-phase EOS required to predict the response of matter at extreme conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Investigation of nonlinear mechanical behavior of two superfine-grained graphites with DIC assisted disc splitting test

Using standardized uniaxial tensile test specimens is not practical due to the limited volumes in irradiation capsules, molten salt degradation facilities, or oxidation apparatus. The ASTM International Standard Test Method for Tensile Strength Estimate by Disc Compression of Manufactured Graphite (ASTM D8289) was developed to provide a convenient way for estimating tensile strength. Unlike the traditional uniaxial tensile test in ASTM C749 (International Standard Test Method for Tensile Stress–Strain of Carbon and Graphite), which uses dog-bone shaped specimens larger than 12.95 mm × 120.65 mm, the ASTM D8289 standard uses smaller discs with diameters of 6–12.7 mm. A digital image correlation (DIC) system, which uses a full-field noncontact surface displacement measurement technique, was applied along with the ASTM D8289 disc splitting test on IG-110 samples and compared with previous measurements from Mersen 2114 samples. Results confirmed that the DIC technique can measure the surface displacement/strain on these small (Ø6 mm × 3 mm) graphite specimens with good repeatability. However, Mersen 2114 and IG-110 samples exhibited strain discrepancies when DIC measurements were compared with analytical and finite element simulation values. The loading history and strain results also indicated different mechanical behaviors between Mersen 2114 and IG-110, particularly the nonlinear behavior of the IG-110 samples. Good agreement was observed by comparing the splitting tensile strengths of two superfine-grained grades with results from other work. The specimen size effect is discussed when comparing the splitting tensile strength with corresponding uniaxial tensile strength of these two graphite grades.

Lin, Lianshan [ORNL] (ORCID:0000000203399219)↗

Spline based least squares integration for two-dimensional shape or wavefront reconstruction

In this paper, we present a novel method to handle two-dimensional shape or wavefront reconstruction from its slopes. The proposed integration method employs splines to fit the measured slope data with piecewise polynomials and uses the analytical polynomial functions to represent the height changes in a lateral spacing with the pre-determined spline coefficients. The linear least squares method is applied to estimate the height or wavefront as a final result. Numerical simulations verify that the proposed method has less algorithm errors than two other existing methods used for comparison. Especially at the boundaries, the proposed method has better performance. The noise influence is studied by adding white Gaussian noise to the slope data. Finally, experimental data from phase measuring deflectometry are tested to demonstrate the feasibility of the new method in a practical measurement.

42 ENGINEERING↗

Determination of Carbonyl Groups in Pyrolysis Bio-oils Using Potentiometric Titration: Review and Comparison of Methods

Carbonyl compounds present in bio-oils are known to be responsible for bio-oil property changes upon storage and during upgrading. As such, carbonyl content has previously been used as a method of tracking bio-oil aging and condensation reactions with less variability than viscosity measurements. Given the importance of carbonyls in bio-oils, accurate analytical methods for their quantification are very important for the bio-oil community. Potentiometric titration methods based on carbonyl oximation have long been used for the determination of carbonyl content in pyrolysis bio-oils. Here in this study, we present a modification of the traditional carbonyl oximation procedures that results in less reaction time, smaller sample size, higher precision, and more accurate carbonyl determinations. Some compounds such as carbohydrates are not measured by the traditional method (modified Nicolaides method), resulting in low estimations of the carbonyl content. Furthermore, we have shown that reaction completion for the traditional method can take up to 300 hours. The new method presented here (the modified Faix method) reduces the reaction time to 2 hours, uses triethanolamine (TEA) in the place of pyridine, and requires a smaller sample size for the analysis. Carbonyl contents determined using this new method are consistently higher than when using the traditional titration methods.

09 BIOMASS FUELS↗

Mitigating contamination in LSS surveys: a comparison of methods

ABSTRACT Future large-scale structure surveys will measure the locations and shapes of billions of galaxies. The precision of such catalogues will require meticulous treatment of systematic contamination of the observed fields. We compare several existing methods for removing such systematics from galaxy clustering measurements. We show how all the methods, including the popular pseudo-Cℓ Mode Projection and Template Subtraction methods, can be interpreted under a common regression framework and use this to suggest improved estimators. We show how methods designed to mitigate systematics in the power spectrum can be used to produce clean maps, which are necessary for cosmological analyses beyond the power spectrum, and we extend current methods to treat the next-order multiplicative contamination in observed maps and power spectra, which reduced power spectrum errors from $\Delta \chi ^2_{\rm C_\ell }\simeq 10$ to ≃ 1 in simulated analyses. Two new mitigation methods are proposed, which incorporate desirable features of current state-of-the-art methods while being simpler to implement. Investigating the performance of all the methods on a common set of simulated measurements from Year 5 of the Dark Energy Survey, we test their robustness to various analysis cases. Our proposed methods produce improved maps and power spectra when compared to current methods, while requiring almost no user tuning. We end with recommendations for systematics mitigation in future surveys, and note that the methods presented are generally applicable beyond the galaxy distribution to any field with spatial systematics.

79 ASTRONOMY AND ASTROPHYSICS↗