共查询到20条相似文献，搜索用时 31 毫秒
1.
On a controversial method for modeling a coregionalization 总被引：2，自引：0，他引：2
P. Goovaerts 《Mathematical Geology》1994,26(2):197-204
This paper reviews two alternative approaches for modeling the (cross) variograms in a coregionalization: (1) fitting the traditional linear model of coregionalization. or (2) deducing the crossvariogram model as a linear combination of prior direct (auto) variogram models while checking the Cauchy-Schwarz inequalities. We show that the second approach has no practical advantage over the traditional one and may not be valid if more than two variables are involved. In such case. Cauchy-Schwarz inequalities are necessary but not sufficient conditions for validity of a coregionalization model. 相似文献
2.
Although there are multiple methods for modeling matrix covariance functions and matrix variograms in the geostatistical literature, the linear coregionalization model is still widely used. In particular it is easy to check to ensure whether the matrix covariance function is positive definite or that the matrix variogram is conditionally negative definite. One of the difficulties in using a linear coregionalization model is in determining the number of basic structures and the corresponding covariance functions or variograms. In this paper, a new procedure is given for identifying the basic structures of the space–time linear coregionalization model and modeling the matrix variogram. This procedure is based on the near simultaneous diagonalization of the sample matrix variograms computed for a set of spatiotemporal lags. A case study using a multivariate spatiotemporal data set provided by the Environmental Protection Agency of Lombardy, Italy, illustrates how nearly simultaneous diagonalization of the empirical matrix variograms simplifies modeling of the matrix variograms. The new methodology is compared with a previous one by analyzing various indices and statistics. 相似文献
3.
In the linear model of coregionalization (LMC), when applicable to the experimental direct variograms and the experimental
cross variogram computed for two random functions, the variability of and relationships between the random functions are modeled
with the same basis functions. In particular, structural correlations can be defined from entries of sill matrices (coregionalization
matrices) under second-order stationarity. In this article, modified t-tests are proposed for assessing the statistical significance of estimated structural correlations. Their specific aspects
and fundamental differences, compared with an existing modified t-test for global correlation analysis with spatial data, are discussed via estimated effective sample sizes, in relation to
the superimposition of random structural components, the range of autocorrelation, the presence of correlation at another
structure, and the sampling scheme. Accordingly, simulation results are presented for one structure versus two structures
(one without and the other with autocorrelation). The performance of tests is shown to be related to the uncertainty associated
with the estimation of variogram model parameters (range, sill matrix entries), because these are involved in the test statistic
and the degrees of freedom of the associated t-distribution through the estimated effective sample size. Under the second-order stationarity and LMC assumptions, the proposed
tests are generally valid. 相似文献
4.
Normal cross-variograms cannot be estimated from data in the usual way when there are only a few points where both variables have been measured. But the experimental pseudo cross-variogram can be computed even where there are no matching sampling points, and this appears as its principal advantage. The pseudo cross-variogram may be unbounded, though for its existence the intrinsic hypothesis alone is not a sufficient stationarity condition. In addition the differences between the two random processes must be second order stationary. Modeling the function by linear coregionalization reflects the more restrictive stationarity condition: the pseudo cross-variogram can be unbounded only if the unbounded correlation structures are the same in all variograms. As an alternative to using the pseudo cross-variogram a new method is presented that allows estimating the normal cross variogram from data where only one variable has been measured at a point. 相似文献
5.
Normal cross-variograms cannot be estimated from data in the usual way when there are only a few points where both variables have been measured. But the experimental pseudo cross-variogram can be computed even where there are no matching sampling points, and this appears as its principal advantage. The pseudo cross-variogram may be unbounded, though for its existence the intrinsic hypothesis alone is not a sufficient stationarity condition. In addition the differences between the two random processes must be second order stationary. Modeling the function by linear coregionalization reflects the more restrictive stationarity condition: the pseudo cross-variogram can be unbounded only if the unbounded correlation structures are the same in all variograms. As an alternative to using the pseudo cross-variogram a new method is presented that allows estimating the normal cross variogram from data where only one variable has been measured at a point. 相似文献
6.
Variogram Fitting by Generalized Least Squares Using an Explicit Formula for the Covariance Structure 总被引：4，自引：0，他引：4
Marc G. Genton 《Mathematical Geology》1998,30(4):323-345
In the context of spatial statistics, the classical variogram estimator proposed by Matheron can be written as a quadratic form of the observations. If data are Gaussian with constant mean, then the correlation between the classical variogram estimator at two different lags is a function of the spatial design matrix and the variance matrix. When data are independent with unidimensional and regular support, an explicit formula for this correlation is available. The same is true for a multidimensional and regular support as can be shown by using Kronecker products of matrices. As variogram fitting is a crucial stage for correct spatial prediction, it is proposed to use a generalized least squares method with an explicit formula for the covariance structure (GLSE). A good approximation of the covariance structure is achieved by taking account of the explicit formula for the correlation in the independent situation. Simulations are carried out with several types of underlying variograms, as well as with outliers in the data. Results show that this technique (GLSE), combined with a robust estimator of the variogram, improves the fit significantly. 相似文献
7.
Spectral simulation of multivariable stationary random functions using covariance fourier transforms 总被引：1，自引：0，他引：1
A coregionalization simulation consists of the generation of realizations of a group of spatially related random variables. The Fourier integral method is presented, modified to carry out such a multivariable simulation. This method allows the simulation of realizations with any specified symmetrical covariance matrix and it is not limited to the classic linear model of coregionalization. The results of gaussian nonconditinal simulations from a case study modeling the spatial characteristics of a layer of coal are given. 相似文献
8.
Automatic Variogram Modeling by Iterative Least Squares: Univariate and Multivariate Cases 总被引：2，自引：1，他引：1
In this paper, we propose a new methodology to automatically find a model that fits on an experimental variogram. Starting with a linear combination of some basic authorized structures (for instance, spherical and exponential), a numerical algorithm is used to compute the parameters, which minimize a distance between the model and the experimental variogram. The initial values are automatically chosen and the algorithm is iterative. After this first step, parameters with a negligible influence are discarded from the model and the more parsimonious model is estimated by using the numerical algorithm again. This process is iterated until no more parameters can be discarded. A procedure based on a profiled cost function is also developed in order to use the numerical algorithm for multivariate data sets (possibly with a lot of variables) modeled in the scope of a linear model of coregionalization. The efficiency of the method is illustrated on several examples (including variogram maps) and on two multivariate cases. 相似文献
9.
《Mathematical Geology》1997,29(6):779-799
Generalized cross-covariances describe the linear relationships between spatial variables observed at different locations.
They are invariant under translation of the locations for any intrinsic processes, they determine the cokriging predictors
without additional assumptions and they are unique up to linear functions. If the model is stationary, that is if the variograms
are bounded, they correspond to the stationary cross-covariances. Under some symmetry condition they are equal to minus the
usual cross-variogram. We present a method to estimate these generalized cross-covariances from data observed at arbitrary
sampling locations. In particular we do not require that all variables are observed at the same points. For fitting a linear
coregionalization model we combine this new method with a standard algorithm which ensures positive definite coregionalization
matrices. We study the behavior of the method both by computing variances exactly and by simulating from various models. 相似文献
10.
Highly Robust Variogram Estimation 总被引：5，自引：0，他引：5
Marc G. Genton 《Mathematical Geology》1998,30(2):213-221
The classical variogram estimator proposed by Matheron is not robust against outliers in the data, nor is it enough to make simple modifications such as the ones proposed by Cressie and Hawkins in order to achieve robustness. This paper proposes and studies a variogram estimator based on a highly robust estimator of scale. The robustness properties of these three estimators are analyzed and compared. Simulations with various amounts of outliers in the data are carried out. The results show that the highly robust variogram estimator improves the estimation significantly. 相似文献
11.
Guillaume Larocque Pierre Dutilleul Bernard Pelletier James W. Fyles 《Mathematical Geology》2007,39(3):263-288
Coregionalization analysis has been presented as a method of multi-scale analysis for multivariate spatial data. Despite an
increasing use of this method in environmental and earth sciences, the uncertainty associated with the estimation of parameters
in coregionalization analysis (e.g., sills and functions of sills) is potentially high and has not yet been characterized.
This article aims to discuss the theory underlying coregionalization analysis and assess the robustness and limits of the
method. A theoretical framework is developed to calculate the ergodic and fluctuation variance-covariance matrices of least-squares
estimators of sills in the linear model of coregionalization. To adjust for the positive semidefiniteness constraint on estimated
coregionalization matrices, a confidence interval estimation procedure for sills and functions of sills is presented. Thereafter,
the relative importance of uncertainty measures (bias and variance) for sills and structural coefficients of correlation and
determination is assessed under different scenarios to identify factors controlling their uncertainty. Our results show that
the sampling grid density, the choice of the least-squares estimator of sills, the positive semidefiniteness constraint, the
presence of scale dependence in the correlations, and the number and range of variogram models, all affect the level of uncertainty,
sometimes through multiple interactions. The asymptotic properties of variogram model parameter estimators in a bounded sampling
domain impose a theoretical limit to their accuracy and precision. Because of this limit, the uncertainty was found to be
high for several scenarios, especially with three variogram models, and was often more dependent on the ratio of variogram
range to domain extent than on the sampling grid density. In practice, in the coregionalization analysis of a real dataset,
the circular requirement for sill estimates in the calculation of uncertainty measures makes the quantification of uncertainty
very problematic, if not impossible. The use of coregionalization analysis must be made with due knowledge of the uncertainty
levels and limits of the method. 相似文献
12.
Nonstationarity of the Mean and Unbiased Variogram Estimation: Extension of the Weighted Least-Squares Method 总被引：1，自引：0，他引：1
When concerned with spatial data, it is not unusual to observe a nonstationarity of the mean. This nonstationarity may be modeled through linear models and the fitting of variograms or covariance functions performed on residuals. Although it usually is accepted by authors that a bias is present if residuals are used, its importance is rarely assessed. In this paper, an expression of the variogram and the covariance function is developed to determine the expected bias. It is shown that the magnitude of the bias depends on the sampling configuration, the importance of the dependence between observations, the number of parameters used to model the mean, and the number of data. The applications of the expression are twofold. The first one is to evaluate a priori the importance of the bias which is expected when a residuals-based variogram model is used for a given configuration and a hypothetical data dependence. The second one is to extend the weighted least-squares method to fit the variogram and to obtain an unbiased estimate of the variogram. Two case studies show that the bias can be negligible or larger than 20%. The residual-based sample variogram underestimates the total variance of the process but the nugget variance may be overestimated. 相似文献
13.
Guillaume Larocque Pierre Dutilleul Bernard Pelletier James W. Fyles 《Mathematical Geosciences》2007,39(3):263-288
Coregionalization analysis has been presented as a method of multi-scale analysis for multivariate spatial data. Despite an increasing use of this method in environmental and earth sciences, the uncertainty associated with the estimation of parameters in coregionalization analysis (e.g., sills and functions of sills) is potentially high and has not yet been characterized. This article aims to discuss the theory underlying coregionalization analysis and assess the robustness and limits of the method. A theoretical framework is developed to calculate the ergodic and fluctuation variance-covariance matrices of least-squares estimators of sills in the linear model of coregionalization. To adjust for the positive semidefiniteness constraint on estimated coregionalization matrices, a confidence interval estimation procedure for sills and functions of sills is presented. Thereafter, the relative importance of uncertainty measures (bias and variance) for sills and structural coefficients of correlation and determination is assessed under different scenarios to identify factors controlling their uncertainty. Our results show that the sampling grid density, the choice of the least-squares estimator of sills, the positive semidefiniteness constraint, the presence of scale dependence in the correlations, and the number and range of variogram models, all affect the level of uncertainty, sometimes through multiple interactions. The asymptotic properties of variogram model parameter estimators in a bounded sampling domain impose a theoretical limit to their accuracy and precision. Because of this limit, the uncertainty was found to be high for several scenarios, especially with three variogram models, and was often more dependent on the ratio of variogram range to domain extent than on the sampling grid density. In practice, in the coregionalization analysis of a real dataset, the circular requirement for sill estimates in the calculation of uncertainty measures makes the quantification of uncertainty very problematic, if not impossible. The use of coregionalization analysis must be made with due knowledge of the uncertainty levels and limits of the method. 相似文献
14.
The experimental variogram computed in the usual way by the method of moments and the Haar wavelet transform are similar in that they filter data and yield informative summaries that may be interpreted. The variogram filters out constant values; wavelets can filter variation at several spatial scales and thereby provide a richer repertoire for analysis and demand no assumptions other than that of finite variance. This paper compares the two functions, identifying that part of the Haar wavelet transform that gives it its advantages. It goes on to show that the generalized variogram of order k=1, 2, and 3 filters linear, quadratic, and cubic polynomials from the data, respectively, which correspond with more complex wavelets in Daubechies's family. The additional filter coefficients of the latter can reveal features of the data that are not evident in its usual form. Three examples in which data recorded at regular intervals on transects are analyzed illustrate the extended form of the variogram. The apparent periodicity of gilgais in Australia seems to be accentuated as filter coefficients are added, but otherwise the analysis provides no new insight. Analysis of hyerpsectral data with a strong linear trend showed that the wavelet-based variograms filtered it out. Adding filter coefficients in the analysis of the topsoil across the Jurassic scarplands of England changed the upper bound of the variogram; it then resembled the within-class variogram computed by the method of moments. To elucidate these results, we simulated several series of data to represent a random process with values fluctuating about a mean, data with long-range linear trend, data with local trend, and data with stepped transitions. The results suggest that the wavelet variogram can filter out the effects of long-range trend, but not local trend, and of transitions from one class to another, as across boundaries. 相似文献
15.
Estimation of non-ergodic variograms and their sampling variance by design-based sampling strategies
Design-based sampling strategies based on classical sampling theory offer unprecedented potentials for estimation of non-ergodic variograms. Unbiased and uncorrelated estimates of the semivariance at the selected lags and of its sampling variance can be simply obtained. These estimates are robust against deviations from an assumed spatial autocorrelation model. The same holds for the variogram model parameters and their sampling (co)variances. Moreover, an objective measure for lack of fit of the fitted model can simply be derived. The estimators for two basic sampling designs, simple random sampling and stratified simple random sampling of pairs of points, are presented. The first has been tested in real world for estimating the non-ergodic variograms of three soil properties. The parameters of variogram models and their sampling (co)variances were estimated with 72 pairs of points distributed over six lags. 相似文献
16.
Covariance and variogram functions have been extensively studied in Euclidean space. In this article, we investigate the validity
of commonly used covariance and variogram functions on the sphere. In particular, we show that the spherical and exponential
models, as well as power variograms with 0<α≤1, are valid on the sphere. However, two Radon transforms of the exponential model, Cauchy model, the hole-effect model and
power variograms with 1<α≤2 are not valid on the sphere. A table that summarizes the validity of commonly used covariance and variogram functions on
the sphere is provided. 相似文献
17.
Application of geostatistical methods in gold geochemical anomalies identification (Montemor-O-Novo, Portugal) 总被引：1，自引：0，他引：1
The study described herein concerns the application of geostatistical methods to data soil from Montemor-O-Novo area (Southern Portugal). In the area, the gold mineralised zones (Banhos, Caeiras, Falés, Gamela, Malaca and Monfurado) are characterised by different geological settings and mineralogical assemblages. A total of 1211 soil samples were collected in Montemor-O-Novo area and analysed for Cu, Pb, Zn, As, Ba and Au by atomic absorption spectrometry.To account for spatial structure, simple and cross variograms were computed for the main directions of the grid sampling. From the experimental variograms a linear model of coregionalization composed of three structures, a nugget effect and two anisotropic spherical structures, was fitted to each of the six variables. The coregionalization matrices deduced from the theoretical model show the relationships between the variables at different scales. These matrices were compared with those obtained by principal component analysis (PCA).This methodology was the basis for estimating the corresponding spatial components (Y0, Y1 and Y2) using factorial kriging analysis (FKA). Maps of raw data, Y0, Y1 and Y2 were made for each variable.The use of multivariate analysis permit the study of the spatial structure intrinsic to geochemical data and the identification and refinement of significant anomalies related to Au-bearing mineral deposits. 相似文献
18.
The application of kriging-based geostatistical algorithms to integrate large-scale seismic data calls for direct and cross variograms of the seismic variable and primary variable (e.g., porosity) at the modeling scale, which is typically much smaller than the seismic data resolution. In order to ensure positive definiteness of the cokriging matrix, a licit small-scale coregionalization model has to be built. Since there are no small-scale secondary data, an analytical method is presented to infer small-scale seismic variograms. The method is applied to estimate the 3-D porosity distribution of a West Texas oil field given seismic data and porosity data at 62 wells. 相似文献
19.
This paper is concerned with vector random fields on spheres with second-order increments, which are intrinsically stationary and mean square continuous and have isotropic variogram matrix functions. A characterization of the continuous and isotropic variogram matrix function on a sphere is derived, in terms of an infinite sum of the products of positive definite matrices and ultraspherical polynomials. It is valid for Gaussian or elliptically contoured vector random fields, but may not be valid for other non-Gaussian vector random fields on spheres such as a χ 2, log-Gaussian, or skew-Gaussian vector random field. Some parametric variogram matrix models are derived on spheres via different constructional approaches. A simulation study is conducted to illustrate the implementation of the proposed model in estimation and cokriging, whose performance is compared with that using the linear model of coregionalization. 相似文献
20.
Kriging with imprecise (fuzzy) variograms. I: Theory 总被引：2，自引：0，他引：2
Imprecise variogram parameters are modeled with fuzzy set theory. The fit of a variogram model to experimental variograms is often subjective. The accuracy of the fit is modeled with imprecise variogram parameters. Measurement data often are insufficient to create good experimental variograms. In this case, prior knowledge and experience can contribute to determination of the variogram model parameters. A methodology for kriging with imprecise variogram parameters is developed. Both kriged values and estimation variances are calculated as fuzzy numbers and characterized by their membership functions. Besides estimation variance, the membership functions are used to create another uncertainty measure. This measure depends on both homogeneity and configuration of the data. 相似文献