鑒于傳統(tǒng)各向異性二維半變異函數(shù)各向同性化方法未充分考慮或無力精確描述其內(nèi)部結(jié)構(gòu)信息的缺陷,本研究通過引入線性廣義尺度不變（GSI）模型,以DEM數(shù)據(jù)作為驗(yàn)證對象,對二維半變異函數(shù)各向異性結(jié)構(gòu)信息進(jìn)行多尺度建模,并采用旋轉(zhuǎn)橢圓法、兩步搜索作圖法等方法對系統(tǒng)參數(shù)進(jìn)行估計(jì),最后以球狀模型為例對理論半變異函數(shù)的估計(jì)精度,及其在空間數(shù)據(jù)插值中的應(yīng)用效果進(jìn)行對比研究。結(jié)果表明:各向異性普遍存在于地形數(shù)據(jù)的空間變異中,有證據(jù)表明,這種各向異性結(jié)構(gòu)中處處顯現(xiàn)出不同的變形特征,但是也存在著某種規(guī)則性的成分,如各向同性圓形或近圓形等值線,因此,在對坐標(biāo)進(jìn)行各向同性化處理時(shí)不適合采用"一刀切"的方式去處理;GSI系統(tǒng)參數(shù)皆能得到較高精度的估計(jì),如決定系數(shù)R²普遍達(dá)到了0.99以上,間接證明了GSI模型對地形數(shù)據(jù)各向異性結(jié)構(gòu)處理的有效性和適用性;通過理論模型估計(jì)和插值結(jié)果對比,線性GSI坐標(biāo)轉(zhuǎn)換法比傳統(tǒng)坐標(biāo)轉(zhuǎn)換法有了明顯的精度提升,并且展現(xiàn)出了較高的邊緣信息恢復(fù)能力,但也表現(xiàn)出了一定的局限性和不穩(wěn)定性。 

【文章來源】：地理研究. 2020,39(11)北大核心CSSCICSCD

【文章頁數(shù)】：19 頁

【部分圖文】：

傳統(tǒng)坐標(biāo)轉(zhuǎn)換法和GSI坐標(biāo)轉(zhuǎn)換法對各向異性半變異函數(shù)等值線處理對比

GSI模型（Generalized Scale Invariance），也叫廣義尺度不變模型，是一種冪律模型，它是Schertzer D等在使用分形理論研究湍流現(xiàn)象中的尺度不變特征提出來的，他們認(rèn)為在自然界中除了一些簡單的自相似現(xiàn)象以外（如兩個(gè)物體經(jīng)過放大和縮小變得相同），還存在一些復(fù)雜的自相似現(xiàn)象（經(jīng)過旋轉(zhuǎn)、拉伸和縮放后變得相同）[20,21]。如前所述，坐標(biāo)轉(zhuǎn)換的核心是旋轉(zhuǎn)和尺度縮放操作，同樣這兩種操作也是GSI的核心思想，它能夠通過這兩種操作兼顧全部等值線的變形信息，從而使半變異函數(shù)結(jié)構(gòu)從整體上得到控制。GSI模型中比較難理解的是尺度概念，同其他冪律模型相比，GSI中尺度定義的方式有所不同，一般冪律模型中尺度是定義在某個(gè)支撐區(qū)間上，如污染物濃度值采樣的某個(gè)時(shí)間段、礦石品位估計(jì)的礦塊體積等物理量上，而GSI將位于同一條等值線（如圖1中的圓和橢圓，也可以是很復(fù)雜的曲線）上的點(diǎn)視為同一尺度，即一條等值線代表一個(gè)尺度，相同尺度上點(diǎn)位的函數(shù)值都是相等的。那么這種尺度應(yīng)如何度量呢？對于分布在圓上的點(diǎn)，可以利用該點(diǎn)到原點(diǎn)的距離作為其尺度的度量，即半徑（圓上每個(gè)點(diǎn)都有一個(gè)尺度值，且都相等）。按照這種定義，對于一些復(fù)雜幾何曲線如橢圓上的點(diǎn)來說，它們到原點(diǎn)的距離是不等的，顯然這種定義在GSI模型中是行不通的，因此需要提出一種新的尺度度量函數(shù)或范式。通常在操作歐氏距離中的點(diǎn)時(shí)，關(guān)注的重點(diǎn)是坐標(biāo)值和其坐標(biāo)構(gòu)成的向徑距離（尺度），不同的是，在這里操作這些點(diǎn)時(shí)，既要習(xí)慣于其在歐氏坐標(biāo)系中的坐標(biāo)，同時(shí)還要關(guān)注其位于某條等值線上，因?yàn)楹笳邲Q定了其尺度的大�。ǚ菤W氏距離）。顯然，可以很輕易地將這種尺度定義與整個(gè)區(qū)域中的單元格一一對應(yīng)起來，對于一個(gè)研究區(qū)域來說，用無數(shù)條等值線去覆蓋整個(gè)區(qū)域，每條曲線代表了一個(gè)尺度，也就意味著對于研究區(qū)域中的任意一點(diǎn)，都能找到劃過該點(diǎn)的等值線，同時(shí)也就確定了它的尺度。如果該系統(tǒng)是自相似的，那么這些不同尺度上的函數(shù)值就可以通過冪律準(zhǔn)則建立起相互之間的關(guān)聯(lián)，意味著一旦已知其中某個(gè)尺度上的函數(shù)值，就能計(jì)算出所有尺度上的值。下面給出GSI模型的符號(hào)描述，如上所述，要建立一個(gè)GSI模型，則需要確定三個(gè)要素：(1)尺度改變算子Tλ；(2)一個(gè)基準(zhǔn)圓或者球E1和一組由其生成的等值線Eλ；(3)等值線尺度的度量函數(shù)φ。冪律關(guān)系是湍流等地球物理領(lǐng)域中的一個(gè)基本準(zhǔn)則，例如在湍流域中，結(jié)構(gòu)函數(shù)滿足S(λ-1x)=λ-2HS(x),S為結(jié)構(gòu)函數(shù)，λ是尺度比（注意與Λ區(qū)分），有λ=φ(E1)/φ(Eλ)，H代表冪階數(shù)，x為位置或者距離向量。上述模型可視為各向同性的，Schertzer D等在此基礎(chǔ)上，將其調(diào)整為S(Tλx)=λ-2HS(x)，其中Tλ=λ-G是尺度算子，G是一矩陣，λ-G為一矩陣函數(shù)，需要將其化為指數(shù)形式然后通過冪級(jí)數(shù)展開進(jìn)行計(jì)算[20,21]�？梢钥闯觯话憬Y(jié)構(gòu)函數(shù)中的各向同性尺度被一個(gè)更一般的尺度算子λ-G所取代，意味著一組同心圓等值線被一組可以任意復(fù)雜的曲線所代替。GSI模型本質(zhì)上就是將一個(gè)尺度轉(zhuǎn)換算子λ-G不斷地作用在一個(gè)圓或者其他曲線上，隨著λ的變化從而生成一組同心等值線的過程，這組等值線僅僅用G中少數(shù)的幾個(gè)參數(shù)就關(guān)聯(lián)了起來。下面看一下φ的定義，對于一個(gè)圓或者球尺度φ(E1)（基準(zhǔn)尺度）來說，很容易找到一個(gè)正數(shù)來定義，如圓或者球的半徑，而對于一個(gè)不規(guī)則形狀的曲線，應(yīng)如何定義呢？從圓或者球半徑的定義受到啟發(fā)，已知圓或者球半徑與圓面積或者球體積之間滿足冪函數(shù)關(guān)系，如其可以表示為φD(E1),D為該半徑φ(E1)的冪。因此，這里可以采用對不規(guī)則曲線的面積或者體的體積進(jìn)行冪運(yùn)算來獲取其尺度，且結(jié)果肯定是一正數(shù)，如果該不規(guī)則曲線為一橢圓，其面積為A，則其尺度可以定義為A1/2。

式中：Z(xi)為某點(diǎn)的函數(shù)值，這里指的是區(qū)域某點(diǎn)的高程值（m);N為某確定間隔的樣點(diǎn)對數(shù)量；xi為某點(diǎn)的向量坐標(biāo)（xi,yi);h為兩點(diǎn)之間的增量，有h=(△xi,△yi)。如果半變異函數(shù)具有冪律特征，那么變換后的公式為γ(Tλh)=λ-sγ(h),s為冪階數(shù)。總之，GSI模型使用少量參數(shù)來關(guān)聯(lián)各級(jí)等值線，將它們納入到一個(gè)統(tǒng)一的框架中，使從全局上控制半變異函數(shù)成為可能。這里需要強(qiáng)調(diào)的是，GSI模型不僅僅是一種冪律模型，更重要的是，它提供了一種坐標(biāo)轉(zhuǎn)換的思路和方法，一種全局各向異性結(jié)構(gòu)各向同性化的方法。如前所述，GSI系統(tǒng)結(jié)構(gòu)是尺度算子λ-G不斷地作用在一個(gè)圓形曲線上從而生成的一組各向異性復(fù)雜曲線，那么它的逆運(yùn)算就是將一組各向異性的復(fù)雜曲線轉(zhuǎn)換成了同一個(gè)圓，然后再分別除以各級(jí)等值線的尺度比，從而就實(shí)現(xiàn)了各向同性化，后面將會(huì)把這種坐標(biāo)轉(zhuǎn)換方法應(yīng)用于理論半變異函數(shù)模型的擬合和空間數(shù)據(jù)插值中。

【參考文獻(xiàn)】：
期刊論文
[1]基于兩種插值算法的三維地質(zhì)建模對比[J]. 馮波,陳明濤,岳冬冬,李勝濤,賈小豐,宋丹.  吉林大學(xué)學(xué)報(bào)(地球科學(xué)版). 2019(04)
[2]復(fù)雜地形地區(qū)月平均氣溫(混合)地理加權(quán)回歸克里格插值[J]. 聶磊,舒紅,劉艷.  武漢大學(xué)學(xué)報(bào)(信息科學(xué)版). 2018(10)
[3]基于RASM的緊支撐徑向基函數(shù)自適應(yīng)并行地形插值方法[J]. 呂海洋,盛業(yè)華,李佳,段平,張思陽.  武漢大學(xué)學(xué)報(bào)(信息科學(xué)版). 2017(09)
[4]各向異性的分形地形生成方法研究[J]. 夏偉杰,周建江,姚楠.  中國圖象圖形學(xué)報(bào). 2009(11)



本文編號(hào)：3501261