Abstract
Remote sensing and geophysical imaging techniques are often limited in terms of spatial resolution. This prevents the characterization of physical properties and processes at scales finer than the spatial resolution provided by the imaging sensor. In the last decade, multiple-point statistics simulation has been successfully used for downscaling problems. In this approach, the missing fine-scale structures are imported from a training image which describes the correspondence between coarse and equivalent fine-scale structures. However, in many cases, large variations in the amplitude of the imaged physical attribute, known as trends, pose a challenge for the detection and simulation of these fine-scale features. Here, we develop a novel multiple-point statistics simulation method for downscaling coarse-resolution images with trends. The proposed algorithm relies on a multi-scale sequential simulation framework. Trends in the data are handled by an inbuilt decomposition of the target variable into a deterministic trend component and a stochastic residual component at multiple scales. We also introduce the application of kernel weighting for computing distances between data events and probability aggregation operations for integrating different support data based on a distance-to-probability transformation function. The algorithm is benchmarked against two-point and multiple-point statistics simulation methods, and a deterministic interpolation technique. Results show that the approach is able to cope with non-stationary data sets and scenarios in which the statistics of the training image differ from the conditioning data statistics. Two case studies using digital elevation models of mountain ranges in Switzerland illustrate the method.
This is a preview of subscription content, log in via an institution to check access.
Similar content being viewed by others
Change history
10 October 2019
The original version of this article unfortunately contained mistakes in equations 9, 10, 12 and 13.
References
Allard D, Comunian A, Renard P (2012) Probability aggregation methods in geoscience. Math Geosci 44(5):545–581
Atkinson PM (2009) Issues of uncertainty in super-resolution mapping and their implications for the design of an inter-comparison study. Int J Remote Sens 30(20):5293–5308
Atkinson PM (2013) Downscaling in remote sensing. Int J Appl Earth Obs Geoinf 22:106–114
Bertero M, Boccacci P (1998) Introduction to inverse problems in imaging. IOP Publishing Ltd, Bristol
Boucher A (2009) Sub-pixel mapping of coarse satellite remote sensing images with stochastic simulations from training images. Math Geosci 41(3):265–290
Boucher A, Kyriakidis PC (2007) Integrating fine scale information in super-resolution land-cover mapping. Photogramm Eng Remote Sens 73(8):913–921
Burt P, Adelson E (1983) The Laplacian pyramid as a compact image code. IEEE Trans Commun 31(4):532–540
Dimitrakopoulos R, Mustapha H, Gloaguen E (2010) High-order statistics of spatial random fields: exploring spatial cumulants for modeling complex non-Gaussian and non-linear phenomena. Math Geosci 42(1):65–99
Farr TG, Rosen PA, Caro E, Crippen R, Duren R, Hensley S, Kobrick M, Paller M, Rodriguez E, Roth L, Seal D, Shaffer S, Shimada J, Umland J, Werner M, Oskin M, Burbank D, Alsdorf D (2007) The Shuttle Radar Topography Mission. Rev Geophys 45(2):1–13
Goovaerts P (1997) Geostatistics for natural resources evaluation. Oxford University Press, New York
Guardiano FB, Srivastava RM (1993) Multivariate geostatistics: beyond bivariate moments. Quantitative geology and geostatistics. Springer, Berlin, pp 133–144
Hengl T, Bajat B, Blagojević D, Reuter HI (2008) Geostatistical modeling of topography using auxiliary maps. Comput Geosci 34(12):1886–1899
Hoffimann J, Scheidt C, Barfod A, Caers J (2017) Stochastic simulation by image quilting of process-based geological models. Comput Geosci 106:18–32
Journel AG (1999) Markov models for cross-covariances. Math Geol 31(8):955–964
Journel AG, Deutsch CV (1993) Entropy and spatial disorder. Math Geol 25(3):329–355
Journel AG, Huijbregts CJ (1978) Mining geostatistics. Academic Press, NewYork
Journel AG, Rossi ME (1989) When do we need a trend model in kriging? Math Geol 21(7):715–739
Keys R (1981) Cubic convolution interpolation for digital image processing. IEEE Trans Acoust Speech Signal Process 29(6):1153–1160
Kwatra V, Schödl A, Essa I, Turk G, Bobick A (2003) Graphcut textures: Image and video synthesis using graph cuts. ACM Trans Graph 22(3):277–286
Kyriakidis PC, Shortridge AM, Goodchild MF (1999) Geostatistics for conflation and accuracy assessment of digital elevation models. Int J Geogr Inf Sci 13(7):677–707
Kyriakidis PC, Yoo EH (2005) Geostatistical prediction and simulation of point values from areal data. Geogr Anal 37(2):124–151
Liu Y, Journel AG (2009) A package for geostatistical integration of coarse and fine scale data. Comput Geosci 35(3):527–547
Mariethoz G, Renard P, Straubhaar J (2010) The direct sampling method to perform multiple-point geostatistical simulations. Water Resour Res 46(11):W11536
Mariethoz G, Renard P, Straubhaar J (2011) Extrapolating the fractal characteristics of an image using scale-invariant multiple-point statistics. Math Geosci 43(7):783–797
Meerschman E, Pirot G, Mariethoz G, Straubhaar J, Meirvenne MV, Renard P (2013) A practical guide to performing multiple-point statistical simulations with the direct sampling algorithm. Comput Geosci 52:307–324
Oriani F, Ohana-Levi N, Marra F, Straubhaar J, Mariethoz G, Renard P, Karnieli A, Morin E (2017) Simulating small-scale rainfall fields conditioned by weather state and elevation: a data-driven approach based on rainfall radar images. Water Resour Res 53(10):8512–8532
Ravalec ML, Noetinger B, Hu LY (2000) The FFT moving average (FFT-MA) generator: an efficient numerical method for generating and conditioning Gaussian simulations. Math Geol 32(6):701–723
Remy N, Boucher A, Wu J (2009) Applied Geostatistics with SGeMS: A User’s Guide. Cambridge University Press, Cambridge
Renard P, Allard D (2013) Connectivity metrics for subsurface flow and transport. Adv Water Resour 51:168–196
Straubhaar J, Renard P, Mariethoz G (2016) Conditioning multiple-point statistics simulations to block data. Spat Stat 16:53–71
Strebelle S (2002) Conditional simulation of complex geological structures using multiple-point statistics. Math Geol 34(1):1–21
Tang Y, Zhang J, Jing L, Li H (2015) Digital elevation data fusion using multiple-point geostatistical simulation. IEEE J Sel Top Appl Earth Observ Remote Sens 8(10):4922–4934
Tran TT (1994) Improving variogram reproduction on dense simulation grids. Comput Geosci 20(7–8):1161–1168
Wang Z, Bovik A, Sheikh H, Simoncelli E (2004) Image quality assessment: from error visibility to structural similarity. IEEE Trans Image Process 13(4):600–612
Wiederkehr M, Möri A (2013) SwissALTI3D: a new tool for geological mapping. Bulletin fuer Angewandte Geologie 18(1):61–69
Zagayevskiy Y, Deutsch CV (2015) Multivariate grid-free geostatistical simulation with point or block scale secondary data. Stoch Env Res Risk Assess 30(6):1613–1633
Zhang T, Switzer P, Journel A (2006) Filter-based classification of training image patterns for spatial simulation. Math Geol 38(1):63–80
Acknowledgements
The authors would like to thank Raphael Nussbaumer for providing the area-to-point simulation code and for the insightful discussions on adapting the approach for DEM downscaling. We are also grateful for the constructive comments of two anonymous reviewers which helped to improve the overall quality of the manuscript. This research was funded by the Swiss National Science Foundation, Grant 200021_159756.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Rasera, L.G., Gravey, M., Lane, S.N. et al. Downscaling Images with Trends Using Multiple-Point Statistics Simulation: An Application to Digital Elevation Models. Math Geosci 52, 145–187 (2020). https://doi.org/10.1007/s11004-019-09818-4
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11004-019-09818-4