1 / 39

Geo479/579: Geostatistics Ch17. Cokriging

Geo479/579: Geostatistics Ch17. Cokriging. Introduction. Data sets often contain more than one variable of interest These variables are usually spatially cross-correlated. The Cokriging System.

jafari
Download Presentation

Geo479/579: Geostatistics Ch17. Cokriging

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Geo479/579: GeostatisticsCh17. Cokriging

  2. Introduction Data sets often contain more than one variable of interest These variables are usually spatially cross-correlated

  3. The Cokriging System A method for estimation that minimizes the variance of the estimation error by exploiting the cross-correlation between several variables Cross-correlated information contained in the secondary variable should help reduce the variance of the estimation errors

  4. The Cokriging System When is the secondary variable useful in estimates? Primary variable of interest is under sampled then the only information we have is the cross correlated information

  5. The Cokriging System The cokriging estimate is a linear combination of both primary and secondary data values This is the Equation used in Ordinary Kriging pg 279

  6. The Cokriging System The development of the cokriging system is identical to the development of ordinary kriging system Estimation Error R can be defined as This is a modification of the error estimation in Ordinary Kriging( pg 279)

  7. The Cokriging System Using matrix notation we can write w = { a1, a2, a3,…an, b1, b2, b3,…bm} Z = { U1, U2…..Ui, V1,….Vj}

  8. The Cokriging System Using Equation 9.14 (p216), 12.6 (p283), 16.3 (p372) we can write

  9. This is similar to Chapter 16

  10. The Cokriging System 1) Unbiasedness condition • Note: Other nonbias conditions are also possible

  11. It is similar to unbiasedness in Ordinary Kriging (p281) • We set error at as 0:

  12. The Cokriging System 2) Minimizing error variance Lagrangean Relaxation:

  13. The Cokriging System Lagrangean Relaxation: Original Lagrange parameter: (12.9)

  14. The Cokriging System Equating n+m+2 partial derivatives of Var{R} to zero, we get the following system of equations

  15. This is similar to minimizing the varianves of error in Ordinary Kriging • The set of weights that minimize the error variance under the unbiasedness condition satisfies the following n+1 equations - ordinary kriging system: (12.11) (12.12)

  16. Minimization of the Error Variance • The ordinary kriging system expressed in matrix (12.13) (12.14)

  17. The Cokriging System Positive definiteness must hold for the set of auto- and cross-variograms (Eq16.44, p391).

  18. The Cokriging System If the primary and secondary variables both exist at all data locations and the auto- and cross-variograms are proportional to the same basic model then the cokriging estimates will be identical to the ordinary kriging estimates

  19. A Cokriging Example

  20. A Case Study • Compares cokriging and ordinary kriging • Undersampled variable U is estimated using 275 U & 470 V sample data for cokriging and only the 275 U data for ordinary kriging

  21. Ordinary kriging 275 U values Using eq 17.11 for the variogram model

  22. Cokriging 275 U and 470 V values Using eq 17.11 for the variogram model Two non-bias conditions 1) uses the initial conditions 2) uses only one nonbias condition

  23. In the alternate unbiased condition, the unknown U value is now estimated as a weighted linear combination of nearby V values adjusted by a constant so that their mean is equal to the mean of the U values

  24. Negative estimates occur because of the nonbias condition

  25. Cokriging with two nonbias conditions is less than satisfactory A physical process with both negative and positive weighting scheme is difficult to imagine Cokriging with one nonbias condition considerably improved the spread of errors and bias Though we had to calculate global means of U and V

  26. If the spatial continuity is modeled using semivariograms then they can be converted to covariance values for cokriging matrices by following equation:

  27. The Cokriging System If we want an estimate over a local area A, there are two options: 1) Average of point estimations within A

  28. The Cokriging System 2) Replace all the covariance terms and in point cokriging system, with average covariance values and

  29. The Cokriging System With the unbiasedness conditions, we can calculate error variance as follows

More Related