杜佳威,武 芳,李靖涵,行瑞星,鞏現(xiàn)勇
信息工程大學(xué)地理空間信息學(xué)院,河南 鄭州 450001
線要素化簡(jiǎn)是自動(dòng)制圖綜合中的重要研究?jī)?nèi)容和經(jīng)典研究問(wèn)題[1],當(dāng)前線化簡(jiǎn)方法主要包括頂點(diǎn)壓縮式[2-8]、彎曲化簡(jiǎn)式[9-14]、智能算法[15-17]等。海岸線要素是一種形態(tài)復(fù)雜且具有特定地理意義的線要素,通用線化簡(jiǎn)算法直接用于海岸線化簡(jiǎn)并不合適[18]。
海岸線是海圖綜合中重點(diǎn)關(guān)注的要素之一(本文海圖主要指航海圖),顧及相關(guān)綜合原則[19],雙向緩沖區(qū)算法[20-21]常用于海岸線化簡(jiǎn)應(yīng)用。一般形態(tài)海岸線化簡(jiǎn)相對(duì)容易,如文獻(xiàn)[22]改進(jìn)單調(diào)鏈法識(shí)別彎曲,以骨架線長(zhǎng)度作為彎曲取舍指標(biāo)化簡(jiǎn)海岸線;而復(fù)雜海岸線化簡(jiǎn)往往需要多算子協(xié)同操作[23],研究難度較大。針對(duì)形態(tài)復(fù)雜的溺谷灣海岸線,文獻(xiàn)[18]基于約束Delaunay三角網(wǎng)構(gòu)建河口層次樹(shù)模型,通過(guò)層次剪枝實(shí)現(xiàn)海岸線化簡(jiǎn),但將狹長(zhǎng)的沉溺河谷視為整體直接取舍,忽略了尺度效應(yīng)影響下要素漸進(jìn)變化的自然規(guī)律[5],不符合其隨尺度變化的漸變過(guò)程;改進(jìn)的Visvalingam算法[24]能夠漸進(jìn)化簡(jiǎn)狹長(zhǎng)河谷,但該方法與比例尺間關(guān)系并不明確,且化簡(jiǎn)后容易出現(xiàn)尖角、鋸齒等形態(tài)問(wèn)題及自相交的拓?fù)鋯?wèn)題。此外,以上方法都缺少視覺(jué)約束,化簡(jiǎn)結(jié)果中仍會(huì)存在人眼無(wú)法分辨[3,27]的局部細(xì)節(jié)。
河口灣海岸線是最具代表性的復(fù)雜海岸線之一。結(jié)合現(xiàn)有研究不足,設(shè)計(jì)河口灣海岸線化簡(jiǎn)算法應(yīng)顧及:①建立算法參數(shù)與比例尺間關(guān)系,增強(qiáng)算法對(duì)化簡(jiǎn)結(jié)果的可控性;②總結(jié)人工綜合知識(shí)、經(jīng)驗(yàn),模擬人工化簡(jiǎn)過(guò)程,漸進(jìn)化簡(jiǎn)河口灣海岸線;③引入視覺(jué)約束,確?;?jiǎn)結(jié)果符合目標(biāo)尺度下視覺(jué)要求;④協(xié)調(diào)形態(tài)化簡(jiǎn)與地理特征保持[19]間關(guān)系,既要充分化簡(jiǎn)細(xì)小彎曲及局部細(xì)節(jié),又要顧及地理要素及其對(duì)應(yīng)地理特征保持,避免化簡(jiǎn)過(guò)度。
基于以上考慮,本文提出一種尺度驅(qū)動(dòng)[14]的河口灣海岸線漸進(jìn)化簡(jiǎn)方法:首先分析河口灣海岸線形態(tài)特征、化簡(jiǎn)約束及化簡(jiǎn)過(guò)程,然后以約束Delaunay三角網(wǎng)為支撐構(gòu)建河口灣骨架線二叉樹(shù)模型,在此模型基礎(chǔ)上通過(guò)細(xì)小彎曲刪除或彎曲細(xì)小部分刪除、漸進(jìn)化簡(jiǎn)、局部夸大等步驟實(shí)現(xiàn)河口灣海岸線漸進(jìn)化簡(jiǎn)。
存在河流注入的海灣均屬河口灣范疇[25],河口灣內(nèi)徑流與潮流交匯,環(huán)流結(jié)構(gòu)、水團(tuán)組成、泥沙沉積等復(fù)雜多變[26]。受復(fù)雜地理環(huán)境影響,河口灣海岸線形態(tài)表達(dá)具有以下特點(diǎn):①水流分汊是河口灣常見(jiàn)現(xiàn)象,河口灣海岸線多呈裂葉狀、樹(shù)枝狀、喇叭狀[26]。②從彎曲構(gòu)成上看,河口灣海岸線主要由用于表達(dá)徑流河道、沉溺河谷、峽江河槽等水流流路的明顯條帶狀彎曲構(gòu)成,也存在部分用于表達(dá)次生河口、淺灣等非條帶狀彎曲。蜿蜒、狹長(zhǎng)的條帶狀彎曲使河口灣海岸線明顯區(qū)別于其他線要素,具有更加復(fù)雜的幾何形態(tài)。③河口灣海岸受侵蝕作用與沉積作用影響顯著,河口灣海岸線中存在局部狹窄部分[23],如圖1所示。
圖1 河口灣海岸線表達(dá)特點(diǎn)Fig.1 Cartographic representation of estuary coastline
河口灣海岸線化簡(jiǎn)應(yīng)遵循以下約束:
(1) 應(yīng)用約束:顧及應(yīng)用安全性,海圖綜合中海岸線化簡(jiǎn)應(yīng)遵循“擴(kuò)陸縮?!盵19]原則。
(2) 地理約束:海岸線充分化簡(jiǎn)后應(yīng)保持準(zhǔn)確、豐富的地理特征[10,18-19];河口灣內(nèi)海水流通,海岸線化簡(jiǎn)后不能導(dǎo)致河口灣被“淤堵”,確?;?jiǎn)結(jié)果符合地理規(guī)律。
(3) 視覺(jué)約束:視覺(jué)約束是受尺度效應(yīng)影響的關(guān)鍵約束,化簡(jiǎn)結(jié)果應(yīng)符合目標(biāo)尺度下視覺(jué)要求,即不存在圖上不可視的細(xì)小彎曲或局部細(xì)節(jié)。結(jié)合制圖經(jīng)驗(yàn)[27],可利用以下閾值(均為圖上距離)量化衡量化簡(jiǎn)后的河口灣海岸線是否符合視覺(jué)要求:①不存在小于深度閾值(thd)的條帶狀彎曲;②不存小于面積閾值(ths)的非條帶狀彎曲;③不存在小于寬度閾值(thw)的狹窄部分。
河口灣海岸線化簡(jiǎn)以彎曲作為化簡(jiǎn)單元[19],其化簡(jiǎn)漸進(jìn)性主要表現(xiàn)為兩方面:
(1) 彎曲化簡(jiǎn)的漸進(jìn)性。主要表現(xiàn)為條帶狀狹長(zhǎng)河谷隨尺度變化的漸變過(guò)程,如圖2所示。當(dāng)條帶狀彎曲中既包含目標(biāo)尺度下可視部分又包含不可視部分時(shí),既不能直接保留(不符合視覺(jué)要求),也不能直接刪除(導(dǎo)致大量地理特征丟失),往往對(duì)其中不可視部分有選擇地刪除或夸大,漸進(jìn)化簡(jiǎn)條帶狀彎曲。
(2) 化簡(jiǎn)過(guò)程的漸進(jìn)性。條帶狀彎曲漸進(jìn)化簡(jiǎn)后可能轉(zhuǎn)化為非條帶狀彎曲,非條帶狀彎曲刪除后也可能影響其他彎曲形態(tài)。因此,需要重復(fù)彎曲化簡(jiǎn)過(guò)程,直至充分化簡(jiǎn)為止。
河口灣海岸線漸進(jìn)化簡(jiǎn)后,還需對(duì)其狹窄部分進(jìn)行適當(dāng)夸大以及次生沖突消除等,才得到最終化簡(jiǎn)結(jié)果。
圖2 狹長(zhǎng)河谷漸進(jìn)化簡(jiǎn)Fig.2 The progressive simplification results of a narrow river valley
利用約束Delaunay三角網(wǎng)[9,18]能夠較好地識(shí)別河口灣海岸線中的條帶狀彎曲。依文獻(xiàn)[18]方法對(duì)河口灣海岸線加密、構(gòu)建約束Delaunay三角網(wǎng),提取海岸線海洋一側(cè)三角形并構(gòu)建骨架線網(wǎng)絡(luò)。
采用文獻(xiàn)[18]中三角形分類方法,將三角形分為入??谌切?、分叉三角形、源頭三角形和連接三角形4類?;谝陨戏诸?,本文將源頭三角形非約束邊對(duì)應(yīng)的頂點(diǎn)稱為源頭,分叉三角形重心稱為分叉點(diǎn),入??谌切沃胁慌c其他三角形鄰接的非約束邊稱為入海口。對(duì)圖1中河口灣海岸線加密、構(gòu)建約束Delaunay三角網(wǎng)、構(gòu)建骨架線網(wǎng)絡(luò)后如圖3所示,其中連接三角形、源頭三角形、分叉三角形、入海口三角形分別被染為白色、綠色、藍(lán)色、淺藍(lán)色。
骨架線網(wǎng)絡(luò)是對(duì)河口灣的降維表達(dá),能夠清晰反映河口灣海岸線結(jié)構(gòu)特點(diǎn)。文獻(xiàn)[18]利用骨架線深度多叉樹(shù)對(duì)具有明顯主干特征的樹(shù)枝狀溺谷灣進(jìn)行結(jié)構(gòu)化表達(dá);而形態(tài)多樣的河口灣未必存在明顯主干特征,本文提出一種骨架線二叉樹(shù)模型,用于河口灣海岸線的結(jié)構(gòu)化表達(dá)。按如下步驟初步構(gòu)建骨架線二叉樹(shù):
步驟1:對(duì)骨架線網(wǎng)絡(luò)進(jìn)行拓?fù)浠幚?,得到所有骨架線段,且標(biāo)記為0。其中,與入??谙噙B的骨架線段為根節(jié)點(diǎn),根節(jié)點(diǎn)作為第i(i=1)層節(jié)點(diǎn)執(zhí)行步驟2。
步驟2:將第i層節(jié)點(diǎn)骨架線段均標(biāo)記為-1,若所有骨架線段都標(biāo)記為-1,則二叉樹(shù)構(gòu)建完成;否則,將第i層所有節(jié)點(diǎn)依次執(zhí)行步驟3,得到的所有孩子節(jié)點(diǎn)作為第i+1層節(jié)點(diǎn),令i=i+1,重復(fù)步驟2。
步驟3:若存在與該節(jié)點(diǎn)骨架線段共用相同分叉點(diǎn)且標(biāo)記為0的骨架線段(有且只有兩條),則這兩條骨架線段即為該節(jié)點(diǎn)的孩子節(jié)點(diǎn);否則,該節(jié)點(diǎn)不存在孩子節(jié)點(diǎn)。
本文將葉子骨架線記為L(zhǎng)0(j),L0(j)貫穿的三角形集合稱為葉子流路,記為{T(j)};將非葉子骨架線記為L(zhǎng)1(k),L1(k)貫穿的三角形集合稱為過(guò)渡流路,記為{T(k)}。葉子骨架線與葉子流路、非葉子骨架線與過(guò)渡流路一一對(duì)應(yīng),所有葉子流路與過(guò)渡流路構(gòu)成河口灣。
河口灣海岸線通過(guò)刪除細(xì)小彎曲、部分刪除條帶狀彎曲、漸進(jìn)化簡(jiǎn)、局部夸大等步驟實(shí)現(xiàn)化簡(jiǎn)?;诠羌芫€二叉樹(shù)結(jié)構(gòu)模型,具體實(shí)現(xiàn)步驟如下:①刪除目標(biāo)尺度下圖上不可視的葉子流路以實(shí)現(xiàn)細(xì)小彎曲刪除;②通過(guò)葉子流路“退化”實(shí)現(xiàn)條帶狀彎曲刪除或部分刪除,即從源頭開(kāi)始沿葉子流路向入??谧匪荩匪葸^(guò)程中不可視的狹窄部分隨之刪除,直至足夠可視為止;③重復(fù)骨架線二叉樹(shù)構(gòu)建、葉子流路刪除、退化等過(guò)程實(shí)現(xiàn)漸進(jìn)化簡(jiǎn);④對(duì)漸進(jìn)化簡(jiǎn)后的所有葉子流路、過(guò)渡流路中狹窄部分適當(dāng)拓寬,實(shí)現(xiàn)局部夸大。
為敘述方便,將化簡(jiǎn)前后比例尺分別記為Sa、Sf;線要素L的圖上長(zhǎng)度記為L(zhǎng)en(L);面要素X的圖上面積記為Area(X)。再定義以下兩個(gè)函數(shù):
(1)X=Union({Xi}):返回多個(gè)面要素{Xi}求并后的結(jié)果X。
(2) {Xi}=Diff(X1,X2):返回屬于面要素X1卻不屬于X2的面要素集合{Xi}。
首先,確定葉子流路是條帶狀還是非條帶狀。利用骨架線長(zhǎng)度衡量葉子流路深度,再用葉子流路面積和深度反算其平均寬度,由經(jīng)驗(yàn)可知深度是寬度2倍以上的葉子流路可以視為條帶狀,否則為非條帶狀;然后,基于各自化簡(jiǎn)閾值確定目標(biāo)尺度下葉子流路是否可視。由于彎曲化簡(jiǎn)具有漸進(jìn)性,本節(jié)將目標(biāo)尺度下不可視的葉子流路標(biāo)記為“待刪除”,取舍操作將在3.4節(jié)中介紹。按以下步驟遍歷所有{L0(j)}:
葉子骨架線L0(j)是葉子流路{T(j)}的中心線,Buf(L0(j))的圖上覆蓋范圍即為目標(biāo)尺度下不可視部分。條帶狀葉子流路退化過(guò)程如下:基于Buf(L0(j))與T(j)間關(guān)系判斷三角形是否可視,進(jìn)而確定條帶狀葉子流路中的可視部分與不可視部分;顧及最小可視距離[3](記為svo),將視覺(jué)上近似不可視的可視部分調(diào)整為不可視部分;將包含源頭的不可視部分刪除,實(shí)現(xiàn)葉子流路退化。具體步驟如下:
步驟1:遍歷{T(j)},若T(j)包含于Buf(L0(j)),則T(j)標(biāo)記為“待化簡(jiǎn)”;否則,T(j)標(biāo)記為“保留”。鄰接的待化簡(jiǎn)三角形構(gòu)成待化簡(jiǎn)區(qū)域,記為{D(d)};鄰接的保留三角形構(gòu)成保留區(qū)域,記為{P(p)}。執(zhí)行步驟2。
步驟2:保留區(qū)域P(p)中不被Buf(L0(j))覆蓋的區(qū)域?yàn)槟繕?biāo)尺度下可視部分,記為Vis(p),Vis(p)=Union(Diff(P(p),Buf(L0(j))))。包含于P(p)內(nèi)的部分L0(j)記為L(zhǎng)0(jp),若Len(L0(jp)) 步驟3:遍歷{D(d)},若{D(d)}中存在包含源頭的D(d),則刪除該D(d),即實(shí)現(xiàn)葉子流路退化;否則,不作處理。 按以上方法,圖5(a)中紅邊三角形即為待化簡(jiǎn)三角形,綠邊為保留三角形,連續(xù)的綠邊三角形構(gòu)成保留區(qū)域{P(1)},連續(xù)的紅邊三角形構(gòu)成待化簡(jiǎn)區(qū)域{D(1),D(2)},格網(wǎng)區(qū)域?yàn)镻(1)中的可視部分Vis(1)。若P(1)足夠可視,則刪除D(1)即實(shí)現(xiàn)該葉子流路退化;若P(1)近似不可視,則P(1)與D(1)、D(2)合并成新的D(1),刪除D(1)實(shí)現(xiàn)該葉子流路退化,如圖5(b)所示。 骨架線二叉樹(shù)中互為兄弟節(jié)點(diǎn)的骨架線對(duì)應(yīng)的葉子流路(或過(guò)渡流路)互為兄弟流路,兄弟節(jié)點(diǎn)的父親節(jié)點(diǎn)的骨架線對(duì)應(yīng)的過(guò)渡流路為父親流路。葉子流路刪除后其兄弟流路將與其父親流路構(gòu)成新葉子流路,互為兄弟的葉子流路不應(yīng)同時(shí)刪除。葉子流路取舍、退化后重新構(gòu)建骨架線二叉樹(shù)并重復(fù)此過(guò)程,直至化簡(jiǎn)結(jié)果中不存在需要?jiǎng)h除或退化的葉子流路為止,完成漸進(jìn)化簡(jiǎn),具體步驟如下: 步驟1:按3.2節(jié)中方法遍歷骨架線二叉樹(shù)中所有葉子流路,若存在被標(biāo)記為“待刪除”的{T(j)},則執(zhí)行步驟2;否則,執(zhí)行步驟3。 步驟2:待刪除葉子流路{T(j)}對(duì)應(yīng)骨架線L0(j)的兄弟骨架線記為L(zhǎng)0(j′),{T(j′)}為L(zhǎng)0(j′)對(duì)應(yīng)的葉子流路。當(dāng){T(j′)未被標(biāo)記時(shí),刪除{T(j)}及其標(biāo)記;當(dāng){T(j′)}標(biāo)記為待刪除時(shí),依深度最大準(zhǔn)則進(jìn)行取舍:若Len(L0(j)) 步驟3:按3.3節(jié)方法化簡(jiǎn)葉子流路,若骨架線二叉樹(shù)中所有葉子流路均未發(fā)生退化,則漸進(jìn)化簡(jiǎn)結(jié)束;否則,重新構(gòu)建骨架線二叉樹(shù),執(zhí)行步驟1。 骨架線也是對(duì)應(yīng)葉子流路(或過(guò)渡流路)的中心線,對(duì)葉子流路與過(guò)渡流路的中心線構(gòu)建視覺(jué)緩沖區(qū),以緩沖區(qū)邊界代替狹窄部分岸線實(shí)現(xiàn)局部夸大。具體實(shí)現(xiàn)步驟如下: 步驟1:遍歷3.4節(jié)中的漸進(jìn)化簡(jiǎn)結(jié)果,對(duì)過(guò)渡流路{T(k)},執(zhí)行步驟2;對(duì)葉子流路{T(j)},執(zhí)行步驟3。 圖3 約束Delaunay三角網(wǎng)、骨架線網(wǎng)絡(luò)Fig.3 Constrained Delaunay triangulation and skeletons network 圖4 河口灣海岸線骨架線二叉樹(shù)模型Fig.4 The binary tree structure of estuary skeletons 圖5 葉子流路“退化”Fig.5 Partial deletion of the leaf stream 步驟2:令{Xi}=Diff( Buf(L1(k)),Union({T(k)})),若{Xi}不為空,則Union({Xi},{T(k)}),實(shí)現(xiàn)過(guò)渡流路局部夸大;否則,無(wú)須夸大。 步驟3:令{Xi}=Diff( Buf(L0(j)),Union({T(j)})),易知{Xi}一定不為空,且至少存在包含該葉子流路源頭的Xi′。將Xi′從{Xi}中剔除后,若{Xi}仍不為空,則Union({Xi},{T(k)}),實(shí)現(xiàn)葉子流路局部夸大;否則,無(wú)須夸大。 圖6(a)中{X1,X2}=Diff(Buf(L1(1)),Union({T(1)})),該過(guò)渡流路局部夸大為Union({X1,X2},{T(1)}),如圖6(b)所示。圖7(a)中{X1,X2,X3}=Diff(Buf(L0(1)),Union({T(1)})),除去包含源頭的X1外,該葉子流路局部夸大為Union({X1,X2},{T(1)}),如圖7(b)所示。 圖6 過(guò)渡流路局部夸大Fig.6 Exaggeration of the non-leaf stream 圖7 葉子流路局部夸大Fig.7 Exaggeration of the leaf stream 對(duì)夸大后的所有葉子流路及過(guò)渡流路求并后提取邊界線,剔除入??谕獾倪吔缇€,即為夸大后的河口灣海岸線。 河口灣海岸線漸進(jìn)化簡(jiǎn)過(guò)程中不會(huì)產(chǎn)生自相交的拓?fù)鋯?wèn)題;局部夸大能夠有效消除漸進(jìn)化簡(jiǎn)后河口灣內(nèi)視覺(jué)沖突。通過(guò)判斷夸大后的河口灣海岸線是否貫穿葉子流路、過(guò)渡流路的夸大部分{Xi},即可確定對(duì)應(yīng)夸大處理是否產(chǎn)生拓?fù)鋯?wèn)題。由于夸大程度較小,且夸大處理前河口灣海岸線已經(jīng)過(guò)充分化簡(jiǎn),一般不會(huì)產(chǎn)生拓?fù)鋯?wèn)題。少數(shù)情況,若產(chǎn)生拓?fù)鋯?wèn)題及其他次生沖突,可參照文獻(xiàn)[6,28]中的處理方法進(jìn)行處理。 通常認(rèn)為圖上最小可視距離svo為0.2 mm,海岸線線劃寬度b介于0.08~0.1 mm之間[27]。本文試驗(yàn)令b=0.27磅(約0.094 5 mm),圖上視覺(jué)指標(biāo)設(shè)置如下:thw=2·b+svo;thd=2·thw;ths=thw·thd?;贏rcEngine二次開(kāi)發(fā)實(shí)現(xiàn)本文試驗(yàn)。 圖8為利用本文方法將1∶25萬(wàn)河口灣海岸線化簡(jiǎn)至1∶50萬(wàn)的主要過(guò)程:圖8(a)為裂葉狀河口灣海岸線,圖8(b)為構(gòu)建約束Delaunay三角網(wǎng)、三角形分類、構(gòu)建骨架線網(wǎng)絡(luò)的結(jié)果,圖8(c)為漸進(jìn)化簡(jiǎn)后的結(jié)果,圖8(d)為局部夸大后的最終化簡(jiǎn)結(jié)果。利用本文方法化簡(jiǎn)喇叭狀、樹(shù)枝狀河口灣海岸線,化簡(jiǎn)結(jié)果(藍(lán)色)與原線(黑色)疊加顯示分別如圖9、圖10所示。 參數(shù)設(shè)置上看:比例尺是本文方法參數(shù)中的唯一變量,調(diào)整目標(biāo)比例尺即可得到滿足該尺度視覺(jué)要求的化簡(jiǎn)結(jié)果,符合尺度驅(qū)動(dòng)的制圖規(guī)律,方法可控性強(qiáng)。 化簡(jiǎn)過(guò)程上看:重復(fù)葉子流路退化、取舍的漸進(jìn)化簡(jiǎn)過(guò)程能夠充分化簡(jiǎn)細(xì)小彎曲或彎曲細(xì)小部分,如圖8(e)中1—3處所示;通過(guò)局部夸大,消除河口灣內(nèi)視覺(jué)沖突,如圖8(e)中4—7處所示。多算子協(xié)調(diào)下由局部向整體的漸進(jìn)變化過(guò)程與人工化簡(jiǎn)過(guò)程相符合。 化簡(jiǎn)結(jié)果上看:各化簡(jiǎn)結(jié)果符合目標(biāo)尺度下視覺(jué)要求,且具有良好的幾何相似性和地理一致性;均未產(chǎn)生次生沖突,各化簡(jiǎn)結(jié)果不存在自相交或視覺(jué)上的自相交,具有良好的拓?fù)湟恢滦裕缓涌跒郴?jiǎn)后未被淤堵(或視覺(jué)上被淤堵),化簡(jiǎn)結(jié)果符合地理規(guī)律;狹長(zhǎng)河谷漸進(jìn)化簡(jiǎn),利于地理特征保持,如圖10中區(qū)域1、2被充分化簡(jiǎn)的同時(shí),主要地理要素及其對(duì)應(yīng)地理特征得到保留;海岸線化簡(jiǎn)后(除局部夸大部分外)位于原線海洋一側(cè),符合應(yīng)用約束。 方法適用性上看,本文方法化簡(jiǎn)裂葉狀、喇叭狀、樹(shù)枝狀河口灣海岸線均能得到良好的化簡(jiǎn)結(jié)果,分別如圖8、圖9、圖10所示。本文方法適用于不同形態(tài)河口灣海岸線多尺度化簡(jiǎn)應(yīng)用。 圖8 河口灣海岸線化簡(jiǎn)過(guò)程Fig.8 Simplification process of the estuary coastline by the proposed method 圖9 利用本文方法化簡(jiǎn)喇叭狀河口灣海岸線結(jié)果Fig.9 Trumpet-like estuary coastline simplified results by the proposed method 圖10 利用本文方法化簡(jiǎn)樹(shù)枝狀河口灣海岸線結(jié)果Fig.10 Dendritic estuary coastline simplified results by the proposed method 以圖8中1∶25萬(wàn)河口灣海岸線為化簡(jiǎn)對(duì)象,對(duì)比試驗(yàn)結(jié)果如圖11—13所示。分析可知: D-P算法化簡(jiǎn)閾值設(shè)置繁瑣,化簡(jiǎn)結(jié)果形態(tài)尖銳且存在自相交的拓?fù)鋯?wèn)題,也存在目標(biāo)尺度下不可視的局部細(xì)節(jié),不符合視覺(jué)約束;此外,化簡(jiǎn)結(jié)果也不滿足應(yīng)用約束。D-P算法化簡(jiǎn)效果較差,不適于河口灣海岸線化簡(jiǎn)應(yīng)用。 圖11 利用D-P算法化簡(jiǎn)河口灣海岸線Fig.11 Results simplified by D-P algorithm 本文方法與雙向緩沖區(qū)算法化簡(jiǎn)效果相對(duì)較好,都不存在自相交的拓?fù)鋯?wèn)題,也都滿足應(yīng)用約束。幾何層次上看,本文方法化簡(jiǎn)結(jié)果形態(tài)良好,較雙向緩沖區(qū)算法更好地保持了原海岸線裂葉狀形態(tài)特征,圖12、13中區(qū)域1對(duì)比明顯;雙向緩沖區(qū)算法化簡(jiǎn)結(jié)果中包含諸多連續(xù)細(xì)小彎曲構(gòu)成的“鋸齒”(如圖12中區(qū)域3所示),即不美觀也不符合人工綜合規(guī)律。地理層次上看,這些細(xì)小鋸齒會(huì)引發(fā)錯(cuò)誤地理認(rèn)知,破壞化簡(jiǎn)前后地理一致性;且雙向緩沖區(qū)算法較本文方法同尺度化簡(jiǎn)結(jié)果丟失了顯著河谷、海灣等重要地理要素,圖12、13中區(qū)域2對(duì)比明顯。 計(jì)算雙向緩沖區(qū)法與本文方法化簡(jiǎn)結(jié)果的位置誤差[30],用于表征線化簡(jiǎn)后的平均位移和整體變形,分別記為Eb、Em。設(shè)計(jì)參數(shù)Ebm=(Eb-Em)/Eb,若Ebm>0,則本文方法化簡(jiǎn)后產(chǎn)生的平均位移小于雙向緩沖區(qū)法,且Ebm越大,本文方法平均位移相對(duì)越?。蝗鬍bm<0,則本文方法產(chǎn)生的平均位移大于雙向緩沖區(qū)法;若Ebm=0,則兩方法產(chǎn)生的平均位移相等。 圖8中1∶25萬(wàn)河口灣海岸線多尺度化簡(jiǎn)結(jié)果的位置誤差如表1所示,分析可知:比例尺跨度在4倍以內(nèi)時(shí),各結(jié)果均符合目標(biāo)尺度下視覺(jué)要求而本文方法產(chǎn)生的平均位移及整體變形更小,這是由雙向緩沖區(qū)算法刪除了明顯海灣、溺谷等地理要素引起的,如圖12、13中1∶50萬(wàn)~1∶100萬(wàn)化簡(jiǎn)結(jié)果中區(qū)域2所示;隨比例尺跨度增大,兩方法化簡(jiǎn)結(jié)果均趨于概略,產(chǎn)生的平均位移與整體變形也趨于相同;當(dāng)比例尺跨度較大時(shí),可能產(chǎn)生本文方法的位置誤差略大于雙向緩沖區(qū)算法的情況,這主要是由雙向緩沖區(qū)算法化簡(jiǎn)不徹底、化簡(jiǎn)后仍存在細(xì)小鋸齒等引起的,如圖12、13中1∶125萬(wàn)化簡(jiǎn)結(jié)果中區(qū)域3等所示,結(jié)合幾何次、地理層次分析可知,此種情況下本文方法仍具有一定優(yōu)勢(shì)。 表1雙向緩沖區(qū)算法與本文方法化簡(jiǎn)結(jié)果的位置誤差 Tab.1Locationerrorsofresultssimplifiedbydouble-buffermethodandtheproposedmethod 比例尺Eb/mEm/mEbm/(%)1∶30萬(wàn)116.6028.2075.811∶40萬(wàn)236.7273.6168.901∶50萬(wàn)247.2994.8261.661∶60萬(wàn)262.50139.6946.781∶70萬(wàn)292.95192.9533.141∶75萬(wàn)296.14207.8229.821∶80萬(wàn)303.82215.7628.981∶100萬(wàn)359.30290.7219.091∶125萬(wàn)389.58410.36-5.33 圖12 利用雙向緩沖區(qū)算法化簡(jiǎn)河口灣海岸線Fig.12 Results simplified by double-buffer algorithm 本文基于Delaunay三角網(wǎng)構(gòu)建骨架線二叉樹(shù)模型,在此模型基礎(chǔ)上提出一種河口灣海岸線漸進(jìn)化簡(jiǎn)方法,充分化簡(jiǎn)目標(biāo)尺度下不可視的細(xì)小彎曲或彎曲細(xì)小部分,保留主要彎曲及彎曲主要部分。本文方法具有以下特點(diǎn):方法參數(shù)設(shè)置與比例尺間關(guān)系明確,方法可控性強(qiáng);方法設(shè)計(jì)模擬人工化簡(jiǎn)的漸進(jìn)過(guò)程,符合認(rèn)知規(guī)律;顧及視覺(jué)約束,局部細(xì)節(jié)充分化簡(jiǎn)的同時(shí)保持了豐富、準(zhǔn)確地理特征,化簡(jiǎn)結(jié)果符合目標(biāo)尺度下視覺(jué)要求,且具有良好的幾何相似性、拓?fù)湟恢滦院偷乩硪恢滦?。本文方法可用于多種形態(tài)河口灣海岸線化簡(jiǎn)應(yīng)用。 圖13 利用本文方法化簡(jiǎn)河口灣海岸線Fig.13 Results simplified by the proposed method 此外,如何更好地處理次生沖突及要素間沖突還有待進(jìn)一步研究,海岸帶多要素協(xié)同綜合將作為接下來(lái)的研究重點(diǎn)。本文方法不局限于航海圖,也可用于海底地形圖等其他海圖河口灣海岸線化簡(jiǎn)應(yīng)用,但針對(duì)不同應(yīng)用范疇,如何量化衡量、協(xié)調(diào)化簡(jiǎn)中地理約束、應(yīng)用約束及視覺(jué)約束間關(guān)系仍有待進(jìn)一步研究。 參考文獻(xiàn): [1] 武芳, 鞏現(xiàn)勇, 杜佳威. 地圖制圖綜合回顧與前望[J]. 測(cè)繪學(xué)報(bào), 2017, 46(10): 1645-1664. DOI: 10.11947/j.AGCS.2017.20170287. WU Fang, GONG Xianyong, DU Jiawei. Overview of the Research Progress in Automated Map Generalization[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(10): 1645-1664. DOI: 10.11947/j.AGCS.2017.20170287. [2] DOUGLAS D H, PEUCKER T K. Algorithms for the Reduction of the Number of Points Required to Represent a Digitized Line or its Caricature[J]. Cartographica: The International Journal for Geographic Information and Geovisualization, 1973, 10(2): 112-122. [3] LI Zhilin, OPENSHAW S. A Natural Principle for the Objective Generalization of Digital Maps[J]. Cartography and Geographic Information Systems, 1993, 20(1): 19-29. [4] VISVALINGAM M, WHYATT J D. Line Generalization by Repeated Elimination of Points[J]. The Cartographic Journal, 1993, 30(1): 46-51. [5] GUO Qingsheng, BRANDENBERGER C, HURNI L. A Progressive Line Simplification Algorithm[J]. Geo-spatial Information Science, 2002, 5(3): 41-45. [6] RAPOSO P. Scale-specific Automated Line Simplification by Vertex Clustering on a Hexagonal Tessellation[J]. Cartography and Geographic Information Science, 2013, 40(5): 427-443. [7] TIENAAH T, STEFANAKIS E, COLEMAN D. Contextual Douglas-Peucker Simplification[J]. Geomatica, 2015, 69(3): 327-338. [8] 李成名, 郭沛沛, 殷勇, 等. 一種顧及空間關(guān)系約束的線化簡(jiǎn)算法[J]. 測(cè)繪學(xué)報(bào), 2017, 46(4): 498-506. DOI: 10.11947/j.AGCS.2017.20160546. LI Chengming, GUO Peipei, YIN Yong, et al. A Line Simplification Algorithm Considering Spatial Relations between Two Lines[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(4): 498-506. DOI: 10.11947/j.AGCS.2017.20160546. [9] 艾廷華, 郭仁忠, 劉耀林. 曲線彎曲深度層次結(jié)構(gòu)的二叉樹(shù)表達(dá)[J]. 測(cè)繪學(xué)報(bào), 2001, 30(4): 343-348. AI Tinghua, GUO Renzhong, LIU Yaolin. A Binary Tree Representation of Curve Hierarchical Structure in Depth[J]. Acta Geodaetica et Cartographica Sinica, 2001, 30(4): 343-348. [10] 翟仁健, 武芳, 朱麗, 等. 利用地理特征約束進(jìn)行曲線化簡(jiǎn)[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2009, 34(9): 1021-1024, 1089. ZHAI Renjian, WU Fang, ZHU Li, et al. Structured Representation of Curve Shape[J]. Geomatics and Information Science of Wuhan University, 2009, 34(9): 1021-1024, 1089. [11] 杜佳威, 武芳, 李靖涵, 等. 采用多元彎曲組劃分的線要素化簡(jiǎn)方法[J]. 計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào), 2017, 29(12): 2189-2196. DU Jiawei, WU Fang, LI Jinghan, et al. Line Simplification Method Based on Multi-bends Groups Division[J]. Journal of Computer-Aided Design & Computer Graphics, 2017, 29(12): 2189-2196. [12] 彭東亮, 鄧敏, 劉慧敏. 更充分利用獨(dú)立彎曲結(jié)構(gòu)的線狀要素Morphing變換方法[J]. 測(cè)繪學(xué)報(bào), 2014, 43(6): 637-644, 652. DOI: 10.13485/j.cnki.11-2089.2014.0100. PENG Dongliang, DENG Min, LIU Huimin. Morphing Transformation of Linear Features by Using Independent Bend Structures More Sufficiently[J]. Acta Geodaetica et Cartographica Sinica, 2014, 43(6): 637-644, 652. DOI:10.13485/j.cnki.11-2089.2014.0100. [13] QIAN Haizhong, ZHANG Meng, WU Fang. A New Simplification Approach Based on the Oblique-dividing-curve Method for Contour Lines[J]. International Journal of Geo-Information, 2016, 5(9): 153. [14] AI Tinghua, KE Shu, YANG Min, et al. Envelope Generation and Simplification of Polylines Using Delaunay Triangulation[J]. International Journal of Geographical Information Science, 2017, 31(2): 297-319. [15] 武芳, 鄧紅艷. 基于遺傳算法的線要素自動(dòng)化簡(jiǎn)模型[J]. 測(cè)繪學(xué)報(bào), 2003, 32(4): 349-355. WU Fang, DENG Hongyan. Using Genetic Algorithms for Solving Problems in Automated Line Simplification[J]. Acta Geodaetica et Cartographica Sinica, 2003, 32(4): 349-355. [16] JIANG Bin, NAKOS B. Line Simplification Using Self-organizing Maps[C/DK]∥A Working Paper Presented at ISPRS Workshop on Spatial Analysis and Decision Making. Hong Kong, China: ISPRS, 2003. [17] 鄭春燕, 郭慶勝, 胡華科. 基于蟻群優(yōu)化算法的線狀目標(biāo)簡(jiǎn)化模型[J]. 測(cè)繪學(xué)報(bào), 2011, 40(5): 635-638. ZHENG Chunyan, GUO Qingsheng, HU Huake. The Simplification Model of Linear Objects Based on Ant Colony Optimization Algorithm[J]. Acta Geodaetica et Cartographica Sinica, 2011, 40(5): 635-638. [18] AI Tinghua, ZHOU Qi, ZHANG Xiang, et al. A Simplification of Ria Coastline with Geomorphologic Characteristics Preserved[J]. Marine Geodesy, 2014, 37(2): 167-186. [19] 王厚祥, 李進(jìn)杰. 海圖制圖綜合[M]. 北京: 測(cè)繪出版社, 1999. WANG Houxiang, LI Jinjie. Chart Cartographic Generalization[M]. Beijing: Surveying and Mapping Press, 1999. [20] CHRISTENSEN A H J. Cartographic Line Generalization with Waterlines and Medial-axes[J]. Cartography and Geographic Information Science, 1999, 26(1): 19-32. [21] 王亮, 吳艷蘭, 湯軍, 等. 改進(jìn)的雙向緩沖海岸線綜合方法研究[J]. 測(cè)繪科學(xué), 2015, 40(1): 18-21, 13. WANG Liang, WU Yanlan, TANG Jun, et al. A Coastline Generalization Method for Marine Delimitation[J]. Science of Surveying and Mapping, 2015, 40(1): 18-21, 13. [22] 陳惠榮, 彭認(rèn)燦, 鄭義東, 等. 以彎曲骨架線為化簡(jiǎn)指標(biāo)的海岸線綜合方法[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2011, 36(12): 1418-1422. CHEN Huirong, PENG Rencan, ZHENG Yidong, et al. Coastline Generalization Based on Skeleton Line of Curve Bends[J]. Geomatics and Information Science of Wuhan University, 2011, 36(12): 1418-1422. [23] WANG Zeshen, MULLER J C. Complex Coastline Generalization[J]. Cartography and Geographic Information Systems, 1993, 20(2): 96-106. [24] VISVALINGAM M, WHELAN J C. Implications of Weighting Metrics for Line Generalization with Visvalingam’s Algorithm[J]. The Cartographic Journal, 2016, 53(3): 253-267. [25] 金長(zhǎng)茂, 眭良仁. 試論河口灣[J]. 海洋學(xué)報(bào), 1989, 11(3): 378-384. JIN Changmao, SUI Liangren. Discussion on the Estuary[J]. Acta Oceanologica Sinica, 1989, 11(3): 378-384. [26] 夏東興, 邊淑華, 豐愛(ài)平, 等. 海岸帶地貌學(xué)[M]. 北京: 海洋出版社, 2014. XIA Dongxing, BIAN Shuhua, FENG Aiping, et al. Coastal Zone Geomorphology[M]. Beijing: China Ocean Press, 2014. [27] 王家耀, 范亦愛(ài), 韓同春, 等. 普通地圖制圖綜合原理[M]. 北京: 測(cè)繪出版社, 1993. WANG Jiayao, FAN Yiai, HAN Tongchun, et al. Cartographic Generalization Theory of General Map[M]. Beijing: Surveying and Mapping Press, 1993. [28] 李振豪, 楊春成, 魏斌, 等. 道路與建筑物空間沖突的幾何相似性移位算法[J]. 測(cè)繪學(xué)報(bào), 2016, 45(6): 747-755. DOI: 10.11947/j.AGCS.2016.20150559. LI Zhenhao, YANG Chuncheng, WEI Bin, et al. A Displacement Algorithm Based on Geometry Similarity for Spatial Conflicts Between Roads and Buildings[J]. Acta Geodaetica et Cartographica Sinica, 2016, 45(6): 747-755. DOI: 10.11947/j.AGCS.2016.20150559. [29] 董箭, 彭認(rèn)燦, 張立華, 等. 顧及"保真性"原則的雙向滾動(dòng)球變換DDM多尺度表達(dá)算法[J]. 測(cè)繪學(xué)報(bào), 2017, 46(6): 789-801. DOI: 10.11947/j.AGCS.2017.20160558. DONG Jian, PENG Rencan, ZHANG Lihua, et al. Multi-scale Representation of Digital Depth Model Based on Double Direction Rolling Ball Transform According to the Reality Principle[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(6): 789-801. DOI: 10.11947/j.AGCS.2017.20160558. [30] 朱鯤鵬, 武芳, 王輝連, 等. Li-Openshaw算法的改進(jìn)與評(píng)價(jià)[J]. 測(cè)繪學(xué)報(bào), 2007, 36(4): 450-456. DOI: 10.3321/j.issn:1001-1595.2007.04.015. ZHU Kunpeng, WU Fang, WANG Huilian, et al. Improvement and Assessment of Li-Openshaw Algorithm[J]. Acta Geodaetica et Cartographica Sinica, 2007, 36(4): 450-456. DOI: 10.3321/j.issn:1001-1595.2007.04.015.3.4 基于骨架線二叉樹(shù)的漸進(jìn)化簡(jiǎn)
3.5 局部夸大
4 試驗(yàn)與分析
4.1 方法有效性驗(yàn)證
4.2 方法優(yōu)越性驗(yàn)證
5 結(jié) 論