EarthRef.org Reference Database (ERR)
Development and Maintenance by the EarthRef.org Database Team

Detailed Reference Information
Goff et al. 2006
Goff, J.A., Jenkins, C. and Calder, B. (2006). Maximum a posteriori resampling of noisy, spatially correlated data. Geochemistry Geophysics Geosystems 7: doi: 10.1029/2006GC001297. issn: 1525-2027.

In any geologic application, noisy data are sources of consternation for researchers, inhibiting interpretability and marring images with unsightly and unrealistic artifacts. Filtering is the typical solution to dealing with noisy data. However, filtering commonly suffers from ad hoc (i.e., uncalibrated, ungoverned) application. We present here an alternative to filtering: a newly developed method for correcting noise in data by finding the "best" value given available information. The motivating rationale is that data points that are close to each other in space cannot differ by "too much," where "too much" is governed by the field covariance. Data with large uncertainties will frequently violate this condition and therefore ought to be corrected, or "resampled." Our solution for resampling is determined by the maximum of the a posteriori density function defined by the intersection of (1) the data error probability density function (pdf) and (2) the conditional pdf, determined by the geostatistical kriging algorithm applied to proximal data values. A maximum a posteriori solution can be computed sequentially going through all the data, but the solution depends on the order in which the data are examined. We approximate the global a posteriori solution by randomizing this order and taking the average. A test with a synthetic data set sampled from a known field demonstrates quantitatively and qualitatively the improvement provided by the maximum a posteriori resampling algorithm. The method is also applied to three marine geology/geophysics data examples, demonstrating the viability of the method for diverse applications: (1) three generations of bathymetric data on the New Jersey shelf with disparate data uncertainties; (2) mean grain size data from the Adriatic Sea, which is a combination of both analytic (low uncertainty) and word-based (higher uncertainty) sources; and (3) side-scan backscatter data from the Martha's Vineyard Coastal Observatory which are, as is typical for such data, affected by speckle noise. Compared to filtering, maximum a posteriori resampling provides an objective and optimal method for reducing noise, and better preservation of the statistical properties of the sampled field. The primary disadvantage is that maximum a posteriori resampling is a computationally expensive procedure.

BACKGROUND DATA FILES

Abstract

Keywords
Mathematical Geophysics, Spatial analysis, Mathematical Geophysics, Stochastic processes (3235, 4468, 4475, 7857), Marine Geology and Geophysics, Seafloor morphology, geology, and geophysics
Journal
Geochemistry Geophysics Geosystems
Publisher
American Geophysical Union
2000 Florida Avenue N.W.
Washington, D.C. 20009-1277
USA
1-202-462-6900
1-202-328-0566
service@agu.org
Click to clear formClick to return to previous pageClick to submit