-
變異函數
鎖定
變異函數(variogram)是描述隨機場(random field)和隨機過程(random process)空間相關性的統計量,被定義為空間內兩空間點之差的方差。在實際應用中,由於無法遍歷空間內所有點,通過有限個採樣計算的變異函數被稱為經驗變異函數(empirical variogram)。變異函數有時也被稱為“變差函數”,在文獻中通常記為2γ(s1, s2)或γ(s1, s2),變異函數的一半被稱為“半變異函數(semivariogram)”,二者本質相同僅存在簡單的倍數關係。變異函數主要應用於隨機場和隨機過程的建模,在地統計學中,克里金法(Kriging)使用變異函數對空間場進行重構和插值。變異函數在穩定過程(stationary process)中存在許多經典模型,包括塊金模型(nugget effect model)、指數模型(exponential model)、高斯模型(Gaussian model)等。
變異函數定義
式中
為隨機場
中支撐集
內的一個點,
為
內任意兩點的向量(注:為與其它研究保持一致,這裏對原作的數學符號進行了修改)。馬瑟倫的定義等價於二階穩定過程中變異函數的定義,而更一般的變異函數定義如下
[2]
:
式中
和
表示數學期望和方差運算,
為隨機場內特定點的數學期望。
當變異函數應用於二階穩定過程(second-order stationary process)時,由於隨機場的數學期望處處相同,且協方差(covariance)僅與兩空間點所構成的向量
有關,此時變異函數可表示為:
變異函數性質
變異函數的數學本質是二階矩(second-order moment),因此恆為非負。且由定義可知,理論變異函數必定經過原點,為偶函數:
若隨機場的協方差
存在,則變異函數與協方差有如下關係:
變異函數參數
塊金(nugget)、基台(sill)和變程(range)是對變異函數進行描述和建模時的常見參數,它們有如下定義
[7]
:
塊金:由性質可知,變異函數必定經過原點,如果變異函數在經過原點之後跳躍至其它值,則該值被稱為塊金。
基台:在二階穩定過程下,隨着距離的增加,變異函數往往趨於平穩,在其呈平穩狀態時所達到的值稱為基台。如果變異函數在二階穩定過程下是可遍歷的(ergodic),則基台即是該隨機過程的方差:
變程:變異函數進入平穩狀態時對應的向量長度稱為變程,在實際應用中,當變異函數的值達到基台的95%以內時,可被認為進入平穩狀態。
變異函數建模
經驗變異函數(empirical variogram)
在實際應用中,隨機場或隨機過程的信息不是處處可用的,只能使用有限個樣本點計算經驗變異函數
,在二階穩定過程下,經驗變異函數有如下表示:
應用上式對格點數據進行計算時可有如下步驟:
1.計算所有樣本點的距離矩陣(distances matrix)
2.設定經驗變異函數的帶寬,並將
按帶寬等分:
3. 通過距離矩陣檢索點間距離(pairwise distance)在h_i+\delta h內的所有點對,並帶入公式計算結果。
4. 遍歷所有
並重復以上步驟
變異函數模型(variogram model)
二階穩定過程不必然是各向同性的,依據可用樣本點所構成向量的方向分別計算經驗變異函數可以觀察該隨機場是否具有各向同性。具有各向同性的隨機場由於所有方向的向量可以共用一個變異函數模型,因此易於建模。一些特定類型的各項異性,例如幾何各向異性(geometric anisotropy)可以通過座標變換轉化為各向同性並進行建模。
變異函數模型並不是任意給定的,由性質可知,在二階穩定過程下,理論變異函數和協方差直接相關(
),因此對變異函數建模時,模型必須滿足協方差矩陣
的正定性質(positive definiteness),即對任意非零向量
,
。以下在各向同性假設下給出常見的變異函數模型
[6]
:
塊金模型(nugget effect model)
塊金模型是對所有距離下的變異函數按定長塊金效應建模所得的結果,是最簡單的變異函數模型。塊金模型下空間所有點的協方差為一常數。
線性模型(linear model)
線性模型中變異函數的取值從塊金處隨距離而增加且是無邊界的(unbounded)。因為
,所以當空間的距離過遠時,協方差可能為負。線性模型無邊界的特性使得其只能在特定的距離範圍內使用。
冪模型(power model)
冪模型與線性模型一樣是無邊界的,具有與線性模型類似的性質,只能在特定的距離範圍內使用。
- 參考資料
-
- 1. Matheron, G. (1963). Principles of geostatistics. Economic geology, 58(8), 1246-1266.
- 2. Cressie, N. (1990). The origins of kriging. Mathematical geology, 22(3), 239-252.
- 3. Webster, R., & Oliver, M. A..Geostatistics for environmental scientists:John Wiley & Sons,2007
- 4. Bachmaier, M., & Backes, M. (2008). Variogram or semivariogram? Understanding the variances in a variogram. Precision Agriculture, 9(3), 173-175.
- 5. Bachmaier, M., & Backes, M. (2011). Variogram or semivariogram? Variance or semivariance? Allan variance or introducing a new term?. Mathematical Geosciences, 43(6), 735-740.
- 6. Le, N. D., & Zidek, J. V. .Statistical analysis of environmental space-time processes:Springer Science & Business Media,2006
- 7. Understanding a semivariogram: The range, sill, and nugget .ArcGIS Pro[引用日期2018-09-23]