專利名稱:一種生成柵格插值的方法和裝置的制作方法
技術(shù)領(lǐng)域:
本申請涉及地理信息領(lǐng)域,特別是涉及一種生成柵格插值的方法和裝置。
背景技術(shù):
GIS (Geographic Information System,地理信息系統(tǒng))經(jīng)過了 40 年的發(fā)展,到今天已經(jīng)逐漸成為一門相當(dāng)成熟的技術(shù),并且得到了極廣泛的應(yīng)用。尤其是近些年,GIS更以其強(qiáng)大的地理信息空間分析功能,在GPS及路徑優(yōu)化中發(fā)揮著越來越重要的作用。GIS地理信息系統(tǒng)是以地理空間數(shù)據(jù)庫為基礎(chǔ),在計(jì)算機(jī)軟硬件的支持下,運(yùn)用系統(tǒng)工程和信息科學(xué)的理論,科學(xué)管理和綜合分析具有空間內(nèi)涵的地理數(shù)據(jù),以提供管理、決策等所需信息的技術(shù)系統(tǒng)。在GIS應(yīng)用中,需要對柵格數(shù)據(jù)進(jìn)行插值計(jì)算,即需要在較大地理空間范圍下快速進(jìn)行柵格插值計(jì)算?,F(xiàn)有技術(shù)中,地理空間范圍較大條件下進(jìn)行插值計(jì)算,目前有以下做法方法一、采用大樣本量計(jì)算。假設(shè)樣本點(diǎn)為P1,P2,...,Pn.對于每一待計(jì)算的柵格V,常用的算法都是根據(jù)樣本點(diǎn)與柵格V的距離來確定一組系數(shù)al,a2,..., an,這系數(shù)是與樣本點(diǎn)對應(yīng)的,表示每個(gè)樣本點(diǎn)在插值結(jié)果中的貢獻(xiàn)度,即Value(V)= al*Pl+a2*P2+. . . +an*Pn ;此種方法精度高,但當(dāng)樣本量η非常大、待插值的柵格非常多的時(shí)候,計(jì)算速度非常慢。方法二、采用小樣本量計(jì)算。計(jì)算原理同上。由于樣本量小,每個(gè)柵格的計(jì)算明顯加快,但是,對于地理空間范圍非常大的條件下,使用較少的采樣點(diǎn),由于插值的結(jié)果沒有實(shí)際采樣調(diào)查準(zhǔn)確可靠,插值的結(jié)果的精確度往往非常低。方法三、改進(jìn)的大樣本量計(jì)算法。采集的樣本點(diǎn)還是Pl,Ρ2,. . .,Pn, (η比較大, 譬如η > 1000),但是,在對每個(gè)柵格計(jì)算插值時(shí),采用大樣本量中較少的點(diǎn)參與實(shí)際計(jì)算, 譬如只選用(1^2,...,011(111值較小,譬如111=10)。這樣,可保證精度的同時(shí)提高速度,但是,由于只針對每個(gè)柵格選用小部分樣本點(diǎn)參與計(jì)算,就會出現(xiàn)相鄰的2個(gè)柵格點(diǎn)A和B,選擇的參與計(jì)算的樣本不同,這樣A和B插值的結(jié)果就有可能出現(xiàn)差異較大,在生成的柵格圖存在明顯的階梯斑快。
發(fā)明內(nèi)容
本申請所要解決的技術(shù)問題是提供一種生成柵格插值的方法和裝置,消除生成柵格圖的存在明顯的階梯斑塊問題。為了解決上述問題,本申請公開了一種生成柵格插值的方法,其特征在于,包括從整個(gè)柵格陣列一端的ηΧη子?xùn)鸥耜嚵虚_始,依次對所有ηΧη子?xùn)鸥耜嚵羞M(jìn)行處理,處理完畢后,將整個(gè)柵格陣列進(jìn)行輸出;其中,相鄰近的兩個(gè)柵格陣列相互交疊至少一排柵格;其中,對于每個(gè)ηΧη子?xùn)鸥耜嚵羞M(jìn)行如下處理所交疊的各柵格的高程值以相互交疊的已進(jìn)行平滑處理的子?xùn)鸥耜嚵械慕Y(jié)果為準(zhǔn);針對該nXn子?xùn)鸥耜嚵兄械奈传@得高程值的每個(gè)柵格,基于對應(yīng)各柵格所選擇的多個(gè)采樣點(diǎn)的各高程值進(jìn)行插值計(jì)算獲得每個(gè)柵格的高程值;根據(jù)所在子?xùn)鸥耜嚵械母鳀鸥竦母叱讨担错樞驅(qū)γ總€(gè)柵格進(jìn)行平滑處理。優(yōu)選的,所述的平滑處理包括如下步驟步驟111,從所述子?xùn)鸥耜嚵衅渲幸粋€(gè)柵格開始,對于該柵格,取所在子?xùn)鸥耜嚵械乃袞鸥竦母叱讨档木底鳛樵摉鸥癞?dāng)前的高程值;步驟112,針對處理順序中的下一個(gè)柵格,取所在子?xùn)鸥耜嚵械乃袞鸥癞?dāng)前的高程值的均值作為該柵格的高程值。優(yōu)選的,根據(jù)如下步驟選擇對應(yīng)所述柵格的多個(gè)采樣點(diǎn)針對子?xùn)鸥耜嚵兄械囊蛔訓(xùn)鸥馮i,計(jì)算該柵格的中心點(diǎn)到所有采樣點(diǎn)的距離Li, 取述各距離Li中最小的m個(gè)距離及其對應(yīng)的采樣點(diǎn)I。優(yōu)選的,所述的通過插值計(jì)算獲得每一子?xùn)鸥竦母叱讨捣椒ò?br>
m ^w ι根據(jù)R = (ZTT)/(Σ7Γ)計(jì)算所述柵格的高程值,其中Vi表示該柵格Ti的高程值。優(yōu)選的,所述的η至少為2。相應(yīng)的,本發(fā)明還公開了一種快速生成柵格插值的裝置,其特征在于,包括主柵格處理模塊和子?xùn)鸥裉幚砟K,所述的主柵格處理模塊用于從整個(gè)柵格陣列一端的nXn子?xùn)鸥耜嚵虚_始,依次對所有nXn子?xùn)鸥耜嚵羞M(jìn)行處理,處理完畢后,將整個(gè)柵格陣列進(jìn)行輸出;其中,相鄰近的兩個(gè)柵格陣列相互交疊至少一排柵格;所述的子?xùn)鸥裉幚砟K包括高程值獲取模塊,用于所交疊的各柵格的高程值以相互交疊的已進(jìn)行平滑處理的柵格陣列的結(jié)果為準(zhǔn);針對該nXn子?xùn)鸥耜嚵兄械奈传@得高程值的每個(gè)柵格,基于對應(yīng)各柵格所選擇的多個(gè)采樣點(diǎn)的各高程值進(jìn)行插值計(jì)算獲得每個(gè)柵格的高程值;平滑處理模塊,用于根據(jù)所在子?xùn)鸥耜嚵械母鳀鸥竦母叱讨?,按順序?qū)γ總€(gè)柵格進(jìn)行平滑處理。優(yōu)選的,所述的平滑處理模塊包括子模塊一,用于從所述子?xùn)鸥耜嚵衅渲幸粋€(gè)柵格開始,對于該柵格,取所在子?xùn)鸥耜嚵械乃袞鸥竦母叱讨档木底鳛樵摉鸥癞?dāng)前的高程值;子模塊二,用于針對處理順序中的下一個(gè)柵格,取所在子?xùn)鸥耜嚵械乃袞鸥癞?dāng)前的高程值的均值作為該柵格的高程值。優(yōu)選的,所述的高程值獲取模塊包括采樣點(diǎn)獲取模塊,用于針對子?xùn)鸥耜嚵兄械囊蛔訓(xùn)鸥瘢?jì)算該柵格的中心點(diǎn)到所有采樣點(diǎn)的距離Li,取述各距離Li中最小的m個(gè)距離及其對應(yīng)的采樣點(diǎn)I。與現(xiàn)有技術(shù)相比,本申請包括以下優(yōu)點(diǎn)本申請通過以整個(gè)柵格陣列中nXn的子?xùn)鸥耜嚵袨閱卧鶕?jù)由采樣點(diǎn)直接獲得的高程值進(jìn)行平滑處理,其中相鄰近的后一子?xùn)鸥耜嚵幸韵嗷ソ化B的已進(jìn)行平滑處理的柵格陣列的結(jié)果為準(zhǔn),循環(huán)進(jìn)行平滑處理,這樣使在生成柵格圖時(shí),速度快,效率高,并且去除了階梯斑塊問題,得到平滑的柵格圖。
圖1是本申請一種本申請一種生成柵格插值的方法的對子?xùn)鸥耜嚵械奶幚砹鞒淌疽鈭D;圖2是本申請一種本申請一種快速生成柵格插值的裝置的結(jié)構(gòu)示意圖;圖3是采用方法三大樣本量計(jì)算得到的柵格圖;圖4是采用方法一大樣本量計(jì)算得到的柵格圖;圖5是采用本申請大樣本量計(jì)算得到的柵格圖。
具體實(shí)施例方式為使本申請的上述目的、特征和優(yōu)點(diǎn)能夠更加明顯易懂,下面結(jié)合附圖和具體實(shí)施方式
對本申請作進(jìn)一步詳細(xì)的說明。本申請通過以整個(gè)柵格陣列中nXn的子?xùn)鸥耜嚵袨閱卧鶕?jù)由采樣點(diǎn)直接獲得的高程值進(jìn)行平滑處理,其中相鄰近的后一子?xùn)鸥耜嚵幸韵嗷ソ化B的已進(jìn)行平滑處理的柵格陣列的結(jié)果為準(zhǔn),循環(huán)進(jìn)行平滑處理,這樣使在生成柵格圖時(shí),速度快,效率高,并且去除了階梯斑塊問題,得到平滑的柵格圖。本申請的一種生成柵格插值的方法包括從整個(gè)柵格陣列一端的nXn子?xùn)鸥耜嚵虚_始,依次對所有nXn子?xùn)鸥耜嚵羞M(jìn)行處理,處理完畢后,將整個(gè)柵格陣列進(jìn)行輸出;其中,相鄰近的兩個(gè)柵格陣列相互交疊至少一排柵格。一個(gè)待處理的Α^Φ(Α,Β值比較大)柵格陣列包含大量的柵格,比如10000Χ 10000 的待處理柵格陣列。一般可以從左下角的一個(gè)nXn子?xùn)鸥耜嚵虚_始,從左向右以nXn子?xùn)鸥耜嚵袨閱卧M(jìn)行處理,其中相鄰近的兩個(gè)柵格陣列相互交疊至少一排柵格,比如對上述Α^Φ柵格陣列,從左下角的第一個(gè)2 X 2子?xùn)鸥耜嚵虚_始進(jìn)行處理,當(dāng)對該2 X 2子?xùn)鸥耜嚵刑幚硗戤吅?,向右選取第二個(gè)2X2子?xùn)鸥耜嚵?,該子?xùn)鸥耜嚵械淖筮呉慌排c前面的子?xùn)鸥耜嚵械挠疫呉慌沤化B,即第一個(gè)子?xùn)鸥耜嚵信c第二個(gè)子?xùn)鸥耜嚵泄靡慌艝鸥瘢绱藢Φ诙€(gè)子?xùn)鸥耜嚵羞M(jìn)行處理,如此循環(huán);當(dāng)最下面一排作為第一排的所有2X2子?xùn)鸥耜嚵刑幚硗瓿珊?,再往上進(jìn)行第二排的2X2子?xùn)鸥耜嚵械奶幚?,?dāng)從下向上對子?xùn)鸥耜嚵羞M(jìn)行處理時(shí),比如說從左下角往上的第二排第一個(gè)2 X 2子?xùn)鸥耜嚵校瑒t該2 X 2子?xùn)鸥耜嚵械南逻呉慌排c左下角的2 X 2子?xùn)鸥耜嚵械纳线呉慌沤化B,如此對第二排的第一個(gè)2X2子?xùn)鸥耜嚵羞M(jìn)行處理;從左至右,對于第二排的第二個(gè)2X2子?xùn)鸥耜嚵校渥筮叺囊慌排c第二排的第一個(gè)2X2子?xùn)鸥耜嚵械挠疫呉慌沤化B,其下邊的一排與第一排的第二個(gè)2X2子?xùn)鸥耜嚵械纳线吔化B,即整個(gè)第二排2 X 2子?xùn)鸥耜嚵械南逻吪艝鸥衽诺? X 2子?xùn)鸥耜嚵猩线呉慌艝鸥窠化B,如此循環(huán),直至將所有2X2子?xùn)鸥耜嚵刑幚硗戤?。其中,子?xùn)鸥竦奶幚眄樞蚩梢杂衅渌绞剑热?,從右上角開始,從左至右,從上之下對子?xùn)鸥耜嚵羞M(jìn)行處理等方式。所述的η至少為2,即子?xùn)鸥耜嚵兄辽贋?X2陣列,比如也可以為3Χ3、4Χ4或者其他。一般情況下,所述的nXn的大小選擇為整個(gè)柵格陣列大小的二千分之一左右,同時(shí)還可以根據(jù)實(shí)際精度的需求確定nXn的大小。
另外,相鄰近的兩個(gè)柵格陣列相互交疊的排數(shù)優(yōu)選為一排柵格,如果nXn相對較大,也可根據(jù)實(shí)際情況選擇交疊多排柵格陣列。在以nXn子?xùn)鸥耜嚵袑φ麄€(gè)柵格陣列進(jìn)行處理時(shí),當(dāng)處理到邊界而原柵格陣列缺少多列或者多行形成nXn子?xùn)鸥耜嚵?,那么可將邊界的一行或者一列?fù)制多行或者多列補(bǔ)齊,在此基礎(chǔ)上完成nXn子?xùn)鸥耜嚵刑幚砗?,將?fù)制的多行或多列的柵格刪除。比如當(dāng)從左至右,從下至上依次以3 X 3子?xùn)鸥耜嚵羞M(jìn)行處理時(shí),相鄰近兩個(gè)3 X 3子?xùn)鸥耜嚵邢嗷ソ化B一排柵格,如果從左至右處理到右邊界時(shí),原柵格陣列只剩余一列柵格,而要使最右邊形成一個(gè)3X3子?xùn)鸥耜嚵?,原柵格陣列需要剩余兩列柵格,此時(shí)可將原柵格陣列最右邊一列復(fù)制添補(bǔ)齊到該列柵格右邊,這樣即可繼續(xù)以3X3子?xùn)鸥耜嚵羞M(jìn)行處理,處理完畢后將所述補(bǔ)齊的一列柵格刪除即可;如果從下至上,也缺少一行,可通過以上相同原理進(jìn)行處理。其它情況可類似進(jìn)行處理。其中,對于每個(gè)nXn子?xùn)鸥耜嚵羞M(jìn)行如下處理參照圖1,示出了本申請一種本申請一種生成柵格插值的方法的對子?xùn)鸥耜嚵械奶幚砹鞒淌疽鈭D,包括步驟101,選定與已進(jìn)行平滑處理的nXn子?xùn)鸥耜嚵邢噜徑南嗷ソ化B至少一排柵格的nXn子?xùn)鸥耜嚵小Mㄟ^前述次序和規(guī)則選擇了 nXn子?xùn)鸥耜嚵修D(zhuǎn)入步驟102。步驟102,所交疊的各柵格的高程值以相互交疊的已進(jìn)行平滑處理的子?xùn)鸥耜嚵械慕Y(jié)果為準(zhǔn);針對該ηXη子?xùn)鸥耜嚵兄械奈传@得高程值的每個(gè)柵格,基于對應(yīng)各柵格所選擇的多個(gè)采樣點(diǎn)的各高程值進(jìn)行插值計(jì)算獲得每個(gè)柵格的高程值。其中,對于一個(gè)nXn子?xùn)鸥耜嚵兄械奈传@得高程值的柵格,可以通過如下方式計(jì)算所述柵格的高程值比如柵格陣列大小為Α*Β,如果樣本量為S1、S2、. . .,Sn,樣本量比較大,確定本方法計(jì)算使用的小樣本量m取值(例如m= 10)。對于子?xùn)鸥耜嚵兄忻總€(gè)要計(jì)算高程值的柵格Ti (1 < = i < = A*B),計(jì)算Ti到每個(gè)樣本Si的距離Li,取Li中最小的m個(gè), 并保存這最小的η個(gè)距離所對應(yīng)的樣本< j <= m)。那么,待插值點(diǎn)網(wǎng)格Ti的值就是由Sj (1 < j <=m)確定??筛鶕?jù)每個(gè)樣本點(diǎn)h到Ti的距離大小,按線性進(jìn)行插值計(jì)算,計(jì)算得到的結(jié)果為Vi,這樣就得到了該子?xùn)鸥耜嚵械母鳀鸥竦母叱讨?。其中,還可以通過其他方法獲取采樣點(diǎn)。其中,所述的通過插值計(jì)算獲得每一子?xùn)鸥竦母叱讨捣椒òǜ鶕?jù)
m ^w ι
h=(ΣΤ^)/(Στγ)計(jì)算所述柵格的高程值,其中Vi表示該柵格i的高程值。還可以通過其
他的插值方法獲得高程值,比如趨勢面法、樣條函數(shù)法、如克立格(Kriging)插值法、謝別德法(Skpard' s Method)等。步驟103,根據(jù)所在子?xùn)鸥耜嚵械母鳀鸥竦母叱讨担错樞驅(qū)γ總€(gè)柵格進(jìn)行平滑處理。對于已獲得高程值的nXn子?xùn)鸥耜嚵?,對其進(jìn)行平滑處理,優(yōu)選的,平滑處理步驟包括步驟111,從所述子?xùn)鸥耜嚵衅渲幸粋€(gè)柵格開始,對于該柵格,取所在子?xùn)鸥耜嚵械乃袞鸥竦母叱讨档木底鳛樵摉鸥癞?dāng)前的高程值;
步驟112,針對處理順序中的下一個(gè)柵格,取所在子?xùn)鸥耜嚵械乃袞鸥癞?dāng)前的高程值的均值作為該柵格的高程值。‘比如對于4個(gè)柵格Pa,Pb, Pc, Pd通過前述步驟獲得的插值結(jié)果Na,Vb, Vc, Vd,然后對這4個(gè)柵格分別進(jìn)行平滑處理,處理過程可采用均值法先對柵格1 進(jìn)行處理,計(jì)算 Va, Vb, Vc, Vd的均值Vh,Pa的高程值取Vh ;然后對柵格1 進(jìn)行處理,計(jì)算Vh,Vb, Vc, Vd 的均值Vt,Pb的高程值取Vt ;再對柵格Pc進(jìn)行處理,計(jì)算Vh,Vt,Vc,Vd的均值Vs,Pc的高程值取Vs ;再對柵格Pd進(jìn)行處理,計(jì)算Vh,Vt, Vs, Vd的均值Vg,Pd的高程值取Vg。這樣,得到比較平滑的4柵格插值結(jié)果。然后,轉(zhuǎn)入步驟101,進(jìn)入下一個(gè)子?xùn)鸥耜嚵械奶幚?。下面以前述對于依次?X2子?xùn)鸥耜嚵袑φ麄€(gè)Α^Φ柵格進(jìn)行處理加以說明對子?xùn)鸥耜嚵械难h(huán)處理過程從左下端開始,從左至右,從下至上的順序依次選擇2 X 2子?xùn)鸥耜嚵?;步驟101先選擇左下角第一個(gè)(也即最下面一排第一個(gè))2Χ2子?xùn)抨嚵校渲?個(gè)柵格均未獲得高程值,則通過上述步驟102計(jì)算4個(gè)柵格的高程值,然后轉(zhuǎn)入步驟103,對該4個(gè)柵格進(jìn)行平滑處理,得到各柵格處理后的高程值;然后進(jìn)入步驟101,選擇最下面一排第二個(gè)2X2子?xùn)鸥耜嚵?,該子?xùn)鸥耜嚵械淖筮呉慌排c第一個(gè)子?xùn)鸥耜嚵杏疫呉慌沤化B,所以該子?xùn)鸥耜嚵械淖筮呉慌诺臇鸥褚呀?jīng)獲得高程值,其余的柵格未獲得高程值,則轉(zhuǎn)入步驟102,使所述未獲得高程值的柵格獲得高程值,然后轉(zhuǎn)入步驟103,對該2X2子?xùn)鸥耜嚵羞M(jìn)行同樣的平滑處理;然后轉(zhuǎn)入下一個(gè)2 X 2子?xùn)鸥耜嚵械奶幚?,其過程如上所述;如此循環(huán)所有2 X 2子?xùn)鸥耜嚵刑幚硗戤吅?,整個(gè)柵格陣列的處理即完畢,則將整個(gè)柵格陣列進(jìn)行輸出。為了便于說明本申請的效果,以實(shí)際大小為4000*4000象素的柵格進(jìn)行地面高度插值計(jì)算,采樣點(diǎn)為130個(gè)點(diǎn)的高程值
8
權(quán)利要求
1.一種生成柵格插值的方法,其特征在于,包括從整個(gè)柵格陣列一端的nXn子?xùn)鸥耜嚵虚_始,依次對所有ηΧη子?xùn)鸥耜嚵羞M(jìn)行處理, 處理完畢后,將整個(gè)柵格陣列進(jìn)行輸出;其中,相鄰近的兩個(gè)柵格陣列相互交疊至少一排柵格;其中,對于每個(gè)ηΧη子?xùn)鸥耜嚵羞M(jìn)行如下處理所交疊的各柵格的高程值以相互交疊的已進(jìn)行平滑處理的子?xùn)鸥耜嚵械慕Y(jié)果為準(zhǔn);針對該ηΧη子?xùn)鸥耜嚵兄械奈传@得高程值的每個(gè)柵格,基于對應(yīng)各柵格所選擇的多個(gè)采樣點(diǎn)的各高程值進(jìn)行插值計(jì)算獲得每個(gè)柵格的高程值;根據(jù)所在子?xùn)鸥耜嚵械母鳀鸥竦母叱讨?,按順序?qū)γ總€(gè)柵格進(jìn)行平滑處理。
2.根據(jù)權(quán)利要求1所述的方法,其特征在于,所述的平滑處理包括如下步驟步驟111,從所述子?xùn)鸥耜嚵衅渲幸粋€(gè)柵格開始,對于該柵格,取所在子?xùn)鸥耜嚵械乃袞鸥竦母叱讨档木底鳛樵摉鸥癞?dāng)前的高程值;步驟112,針對處理順序中的下一個(gè)柵格,取所在子?xùn)鸥耜嚵械乃袞鸥癞?dāng)前的高程值的均值作為該柵格的高程值。
3.根據(jù)權(quán)利要求1所述的方法,其特征在于,根據(jù)如下步驟選擇對應(yīng)所述柵格的多個(gè)采樣點(diǎn)針對子?xùn)鸥耜嚵兄械囊蛔訓(xùn)鸥馮i,計(jì)算該柵格的中心點(diǎn)到所有采樣點(diǎn)的距離Li,取述各距離Li中最小的m個(gè)距離及其對應(yīng)的采樣點(diǎn)I。
4.根據(jù)權(quán)利要求3所述的方法,其特征在于所述的通過插值計(jì)算獲得每一子?xùn)鸥竦母叱讨捣椒ò╩ ^w ι根據(jù)h = (ΣΤ^)/(ΣΤΓ)計(jì)算所述柵格的高程值,其中Vi表示該柵格Ti的高程值。J=I L1 J=I lI
5.根據(jù)權(quán)利要求1所述的方法,其特征在于 所述的η至少為2。
6.一種快速生成柵格插值的裝置,其特征在于,包括 主柵格處理模塊和子?xùn)鸥裉幚砟K,所述的主柵格處理模塊用于從整個(gè)柵格陣列一端的ηΧη子?xùn)鸥耜嚵虚_始,依次對所有ηΧη子?xùn)鸥耜嚵羞M(jìn)行處理,處理完畢后,將整個(gè)柵格陣列進(jìn)行輸出;其中,相鄰近的兩個(gè)柵格陣列相互交疊至少一排柵格; 所述的子?xùn)鸥裉幚砟K包括高程值獲取模塊,用于所交疊的各柵格的高程值以相互交疊的已進(jìn)行平滑處理的柵格陣列的結(jié)果為準(zhǔn);針對該ηΧη子?xùn)鸥耜嚵兄械奈传@得高程值的每個(gè)柵格,基于對應(yīng)各柵格所選擇的多個(gè)采樣點(diǎn)的各高程值進(jìn)行插值計(jì)算獲得每個(gè)柵格的高程值;平滑處理模塊,用于根據(jù)所在子?xùn)鸥耜嚵械母鳀鸥竦母叱讨担错樞驅(qū)γ總€(gè)柵格進(jìn)行平滑處理。
7.根據(jù)權(quán)利要求6所述的裝置,其特征在于 所述的平滑處理模塊包括子模塊一,用于從所述子?xùn)鸥耜嚵衅渲幸粋€(gè)柵格開始,對于該柵格,取所在子?xùn)鸥耜嚵械乃袞鸥竦母叱讨档木底鳛樵摉鸥癞?dāng)前的高程值;子模塊二,用于針對處理順序中的下一個(gè)柵格,取所在子?xùn)鸥耜嚵械乃袞鸥癞?dāng)前的高程值的均值作為該柵格的高程值。
8.根據(jù)權(quán)利要求6所述的裝置,其特征在于 所述的高程值獲取模塊包括采樣點(diǎn)獲取模塊,用于針對子?xùn)鸥耜嚵兄械囊蛔訓(xùn)鸥?,?jì)算該柵格的中心點(diǎn)到所有采樣點(diǎn)的距離Li,取述各距離Li中最小的m個(gè)距離及其對應(yīng)的采樣點(diǎn)I。
全文摘要
本申請?zhí)峁┝艘环N生成柵格插值的方法和裝置,涉及地理信息領(lǐng)域。所述的方法包括從整個(gè)柵格陣列一端的n×n子?xùn)鸥耜嚵虚_始,依次對所有n×n子?xùn)鸥耜嚵羞M(jìn)行處理,處理完畢后,將整個(gè)柵格陣列進(jìn)行輸出;其中,相鄰近的兩個(gè)柵格陣列相互交疊至少一排柵格,并以已進(jìn)行平滑處理子?xùn)鸥耜嚵袨榛鶞?zhǔn)對下一個(gè)子?xùn)鸥耜嚵羞M(jìn)行平滑處理。申請通過以整個(gè)柵格陣列中n×n的子?xùn)鸥耜嚵袨閱卧鶕?jù)由采樣點(diǎn)直接獲得的高程值進(jìn)行平滑處理,其中相鄰近的后一子?xùn)鸥耜嚵幸韵嗷ソ化B的已進(jìn)行平滑處理的柵格陣列的結(jié)果為準(zhǔn),循環(huán)進(jìn)行平滑處理,這樣使在生成柵格圖時(shí),速度快,效率高,并且去除了階梯斑塊問題,得到平滑的柵格圖。
文檔編號G06F17/30GK102339310SQ201110274280
公開日2012年2月1日 申請日期2011年9月15日 優(yōu)先權(quán)日2011年9月15日
發(fā)明者史明昌, 孫成寶, 曹剛, 李團(tuán)宏, 李嵩, 黃兆偉 申請人:北京地拓科技發(fā)展有限公司