一種緊支撐權函數性質的推導及空間插值實驗

李　佳，段　平，盛業華，呂海洋，張思陽

(1.云南師范大學　旅游與地理科學學院，云南 昆明 650050；2.南京師范大學　虛擬地理環境教育部重點實驗室，江蘇 南京 210023；3.江蘇省地理環境演化國家重點實驗室培育建設點，江蘇 南京 210023；4.江蘇省地理信息資源開發與利用協同創新中心，江蘇 南京 210023)

?

一種緊支撐權函數性質的推導及空間插值實驗

李佳1,2,3,4，段平1,2,3,4，盛業華2,3,4，呂海洋2,3,4，張思陽2,3,4

(1.云南師范大學旅游與地理科學學院，云南 昆明 650050；2.南京師范大學虛擬地理環境教育部重點實驗室，江蘇 南京 210023；3.江蘇省地理環境演化國家重點實驗室培育建設點，江蘇 南京 210023；4.江蘇省地理信息資源開發與利用協同創新中心，江蘇 南京 210023)

在空間插值方法中，權函數是影響空間插值精度的重要因素之一，通常需要選擇具有性質較好的權函數進行空間插值。推導一種緊支撐權函數的性質，以鐵礦品位為數據源，采用緊支撐權函數進行空間插值，最終驗證該緊支撐權函數的優良性質：具有較高的插值精度。推導的緊支撐權函數性質可為后期構建新的緊支撐權函數作為一種借鑒。

緊支撐權函數；空間插值；反距離權重；鐵礦品位

地理空間現象是連續變化的，通常為了表達連續的地理空間現象，需要建立連續的空間場。受觀測成本、觀測不可重復性等原因，只能在研究區域中采集到有限的散亂樣本數據，空間插值方法就是利用這些有限的散亂樣本數據建立一個連續函數模型，并計算地理空間任意位置上的特征屬性，最終重建整個空間場的空間分析方法[1-2]。空間插值方法主要確定參考點的選擇和權函數。當參考點確定后，權函數決定插值的精度。

國內外學者對空間插值方法中權函數進行大量的研究，如在經典的Shepard插值方法中權函數是以冪函數為倒數的函數[3]，Bekele，Ping，Lloyd分別選取冪函數的不同次數(1～5)進行插值實驗[4-6]，Rühaak等人認為冪函數的次數應該根據觀測點的可信度來動態變化，Lu 等人則認為冪函數次數應該根據參考點的密度自適應變化[7]，史利民從數學曲面的角度定義幾類新的權函數[8]，徐遵義從權函數的衰減性上對Shepard插值方法中的權函數進行改進[9]，李正泉認為權函數的定義應該顧及插值參考點的方向遮蔽性，提出權函數應該考慮到反映方向的調和權重系數[10]。上述各種權函數的研究集中在構建新的權函數和對原始權函數改進，而對于權函數數學性質的研究較少，實際上權函數的數學性質指導著權函數的構建和改進。基于此，本文證明了一種常見的緊支撐權函數的性質，并采用該緊支撐權函數對模擬數據和實際的采樣數據進行空間插值實驗，實驗表明本文采用的緊支撐權函數是一種插值精度較高的權函數，本文證明的緊支撐函數性質可為后期構建新的緊支撐函數作為一種指導。

1　空間插值方法中常見的權函數

(1)

式中：wk表達式常見的有3種，均是與距離有關的函數：

(2)

式中：μ為冪函數的次數，可選候選值為1、2、3、4、5等，但一般取2。

(3)

(4)

其中

(5)

式中：Rw為設置的緊支撐半徑。

從數學建模上分析，只要權函數的衰減性和光滑性越好，則插值效果越理想。而采用空間插值方法重建整個連續的空間場過程可以理解為數學建模過程。因此，在空間插值方法中，選擇的權函數應該具備較好的衰減性和光滑性。上述各式中權函數是一個衰減性函數，隨著距離的增加權函數對其插值點的貢獻越小，一般的權函數都具備衰減性，而光滑性表現在權函數的偏導數相等，式(4)是一種具有較好的光滑性的緊支撐權函數，給出權函數光滑性的證明。

2　一種緊支撐權函數性質的證明

以二維情況為例，給出式(4)緊支撐權函數的基本性質，在此基礎上給出偏導數相等的證明。

2.1緊支撐權函數的基本性質

當插值點(x, y)趨近于某個已知點(xj, yj)時，即(x,y)→(xj,yj)分為兩種情況：

1)當j=k時，可以認為第k個點在wk處的wk的值為無窮大。

wk(xj,yj)→∞，而權重值為

(6)

分母：

w1+w2+…+wj+…+wk+…+

wn～wj(xj,yj)→∞.

(7)

(8)

(9)

(10)

2.2緊支撐權函數偏導數相等的證明

1)若(x,y)→(xj,yj)=(xk,yk)，則：

wk(xk,yk)→∞，wj(xk,yk)→有界量

(11)

(12)

其中

(13)

(14)

其中

(15)

(16)

(17)

(18)

其中

(19)

(20)

(21)

(22)

2)若(x,y)→(xj,yj)≠(xk,yk)，

(23)

M1，M2與上述相同，均為有界量。則式(23)中各分量為

(24)

(25)

3　鐵礦品位的空間插值實驗與空間場重建

3.1鐵礦品位的空間插值實驗分析

為了驗證本文推導的緊支撐權函數的性質，采用與緊支撐權函數同一類型的其它權函數進行反距離權重(Inverse Distance Weighted，IDW)空間插值實驗對比，插值模型為

(26)

圖1　礦體水平巷道鉆孔數據空間分布

為了消除各種插值方法中參考點不同而影響插值誤差分析，因此需要各種插值方法所采用的參考點數量相同，本文采用最鄰近的30個點作為參考點，插值點與最鄰近的第30個點的距離作為緊支撐權重的半徑Rw。采用交叉驗證統計各種插值方法的最小誤差(Min)、最大誤差(Max)、平均誤差(ME)、均方根誤差(RMSE)，如表1所示。

表1中緊支撐權函數的Min較大，但實際數值也趨近于0，認為緊支撐權函數在Min方面與另外幾種插值方法無差異，而在Max方面，緊支撐權函數插值與式(2)權函數插值誤差近似，但是這二種權函數插值都不如式(3)權函數插值的誤差大，說明緊支撐權函數可能出現某個局部點的誤差較大。就ME和RMSE而言，緊支撐權函數插值方法的這二項誤差指標最小，而這二種誤差指標是評價插值方法總體性能的最重要指標[13]，說明緊支撐權函數插值的總體誤差最小，精度最高。從插值類型分析，與緊支撐權函數相似的幾種權函數插值誤差精度均優于克里金插值和薄板樣條插值，4種誤差最大的是克里金插值，其次是薄板樣條插值，其原因是這二種插值方法所涉及的插值模型參數需要一定的人工經驗，如果設置的不合理可能引起較大誤差。

表1　各種插值方法誤差統計指標

3.2鐵礦品位的空間場重建

分別采用5種插值方法對鐵礦品位進行空間場重建，如圖2所示中(a)、(b)、(c)、(d)、(e)分別為采用式(2)權函數的IDW插值、式(3)權函數的IDW插值、緊支撐權函數插值、薄板樣條插值、克里金插值的鐵礦品位空間場重建結果，各插值結果的邊界出現鋸齒狀，其原因是裁剪的多邊形有鋸齒狀，但不影響空間場重建的分析。

從圖2中可知，采用同一類型權函數插值方法(式(2)和式(3)權函數及緊支撐權函數)建立的空間場整體形態近似，但是采用式(2)和式(3)權函數空間插值方法的空間場呈現出整體“模糊”，而緊支撐權函數構建的空間場呈現出的光滑性和連續性是優于其他二種權函數插值方法，這也與3.1節中分析的插值誤差結論相對應，插值方法誤差越小表現在構建的空間場越光滑和越連續。與其他不同類型的樣條函數插值和克里金插值比較，從空間場形態上分析，樣條函數插值方法與其他方法的形態相差較大，主要原因與在3.1節中的各項誤差最大相應，但是其空間場表明的光滑性和連續性較好，但是從光滑性和連續性不能判定插值方法的好壞，要與精度一起綜合評價插值方法，而采用克里金插值方法重建的空間場形態與緊支撐權函數及其他權函數的空間場形態大體一樣，其空間場整體上光滑性和連續性也較好，但是在3.1節中的誤差不及緊支撐權函數空間插值。

綜上所述，從插值方法的各種誤差指標的總體性能及重建的空間場形態的連續性和光滑性分析，可知本文推導的緊支撐權函數在空間插值中表現較優。

圖2　采用5種插值方法構建的鐵礦品位空間場

4　結　論

緊支撐權函數的數學性質影響插值精度，本文推導一種緊支撐權函數的數學性質，通過實驗驗證采用該緊支撐權函數插值具有較高的插值精度，同時推導的緊支撐權函數性質可為后期構建新的權函數作為一種指導。

[1]段平,盛業華,李佳,等.自適應的IDW插值方法及其在氣溫場中的應用[J].地理研究,2014,33(8):1417-1426.

[2]張思陽,盛業華,段平,等.顧及各向異性的CSRBF空間插值及其在氣溫場重建中的應用[J].地理與地理信息科學,2014,30(7):117-121.

[3]RANKE R.Scattered data interpolation: tests of some methods[J].Mathematics of computation,1982,38(157): 181-200.

[4]BEKELE A,DOWNER R G,WOLCOTT M C,el al.Comparative evaluation of spatial prediction methods in a fieldexperiment for mapping soil potassium [J].Soil Science,2003,168(1):15-28.

[5]PING J L,GREEN C J,ZARTMAN R.Exploring spatial dependence of cotton yield using global and local autocorrelation statistics [J].Field Crops Research,2004,89(2): 219-236

[6]LLOYD C D.Assessing the effect of integrating elevation data intothe estimation of monthly precipitation in Great Britain[J].Journal ofHydrology,2005,308(1): 128-150.

[7]LU G Y,WONG D W.An adaptive inverse-distance weighting spatial interpolation technique[J].Computers & Geosciences,2008,34(9): 1044-1055.

[8]史利民,王仁宏.幾種基于散亂數據擬合的局部插值方法[J].數學研究與評論,2006,26(2): 283-291.

[9]徐遵義,姜玉祥,趙亮,等.改進的Shepard算法及其在重力異常插值中的應用[J].武漢大學學報(信息科學版),2010,35(4):477-480.

[10] 李正泉,吳堯詳.顧及方向遮蔽性的反距離權重插值法[J].測繪學報,2015,44(1): 91-98.

[11] 黃靜,王愛倩,翟世龍.基于Kriging方法的負偏態分布數據的空間插值研究[J].測繪工程,2015,24(8):16-19

[12] 張子昕，周強波.空間插值算法在GIS中的應用[J].測繪與空間地理信息，2015,38(2)：103-107.

[13] 焦永清,李斌,張坤.基于空間可視化的氣溫插值方法比較[J].測繪工程,2013,22(5):62-64.

[責任編輯：李銘娜]

Theoretical derivation of compactly supported weight function andits application to spatial interpolation

LI Jia1,2,3,4,DUAN Ping1,2,3,4, SHENG Yehua2,3,4,LYU Haiyang2,3,4, ZHANG Siyang2,3,4

(1. College of Tourism and Geographical Sciences, Yunnan Normal University, Kunming 650500, China;2. Key Laboratory of Virtual Geographic Environment (Nanjing Normal University), Ministry of Education, Nanjing, 210023, China;3. State Key Laboratory Cultivation Base of Geographical Environment Evolution (Jiangsu Province), Nanjing, 210023, China;4. Jiangsu Center for Collaborative Innovation in Geographical Information Resource Development and Application, Nanjing 210023, China)

The weight function is one of the important factors that influence the spatial interpolation precision in many spatial interpolation methods. Weight function with better mathematical property is selected with spatial interpolation methods. Based on it, the mathematical property of compactly supported weight function is deduced in this paper. Iron ore grade sampling data as data source is used in compactly supported weight function interpolation experiment. By experimental validation, compactly supported weight function with fine properties has better precision and versatility. The mathematical property derivation of compactly supported weight function can be a guide for building a new compactly supported weight function.

compactly supported weight function; spatial interpolation; inverse distance weighted; iron ore grade

10.19349/j.cnki.issn1006-7949.2016.11.002

2016-01-23

國家自然科學基金資助項目(41271383)；云南師范大學博士基金資助項目(01300205020503128)

李佳(1984-)，女，講師，博士.

P223

A

1006-7949(2016)11-0007-05