林春峰 閔世平 黃華平

(中鐵二院工程集團(tuán)有限責(zé)任公司，四川成都 610031)

1 概述

為了提高生產(chǎn)效率，在數(shù)字化制圖中需要利用已知的等高線進(jìn)行內(nèi)插，目前常用方法有網(wǎng)格內(nèi)插、三角網(wǎng)內(nèi)插、輔助線內(nèi)插等［1-7］。網(wǎng)格內(nèi)插和三角網(wǎng)內(nèi)插的計(jì)算量較大，在實(shí)際應(yīng)用中，會(huì)顯得較遲鈍;輔助線內(nèi)插方法計(jì)算速度快，但由于采用的數(shù)據(jù)點(diǎn)較少，數(shù)據(jù)點(diǎn)的分布位置不佳，在一些等高線急劇變化的地方，內(nèi)插出的等高線存在較大失真［8］。結(jié)合輔助線方法的特點(diǎn)，提出一種新的內(nèi)插方法，首先利用特征參數(shù)對(duì)已知的等高線節(jié)點(diǎn)進(jìn)行描述;再利用特征參數(shù)進(jìn)行等高線節(jié)點(diǎn)匹配，將等高線分割成不同的區(qū)域，分區(qū)域進(jìn)行輔助線構(gòu)建，利用構(gòu)建的輔助線內(nèi)插等高線;最后對(duì)內(nèi)插出的等高線進(jìn)行節(jié)點(diǎn)簡化，得到最終結(jié)果。實(shí)驗(yàn)結(jié)果表明，利用該方法內(nèi)插的等高線能夠很好的顧及局部地形，有效地解決等高線內(nèi)插的問題。

2 算法介紹

等高線內(nèi)插要解決的問題可以轉(zhuǎn)化成:在相鄰兩條已知等高線l1和l2中，內(nèi)插N條等高線。其中，假定l1和l2節(jié)點(diǎn)排列順序一致，N為已知常數(shù)。

算法首先按照給定的規(guī)則計(jì)算等高線l1和l2上每個(gè)節(jié)點(diǎn)的特征參數(shù);然后利用匹配規(guī)則進(jìn)行節(jié)點(diǎn)匹配，將l1和l2分割成一一對(duì)應(yīng)的區(qū)域;然后對(duì)各個(gè)區(qū)域進(jìn)行輔助線構(gòu)建，利用構(gòu)建后的輔助線內(nèi)插等高線;最后進(jìn)行節(jié)點(diǎn)簡化，得到最終結(jié)果。

2.1 節(jié)點(diǎn)特征參數(shù)計(jì)算

節(jié)點(diǎn)特征參數(shù)的定義為

Pi為某條等高線上第i個(gè)節(jié)點(diǎn)的特征參數(shù)。其中:i∈［0，n］，n 為等高線節(jié)點(diǎn)的個(gè)數(shù);Xi、Yi和 Zi為節(jié)點(diǎn)的坐標(biāo)值;αi為前一節(jié)點(diǎn)和當(dāng)前節(jié)點(diǎn)連線與水平方向夾角角度值，αi∈［0，360);βi為當(dāng)前節(jié)點(diǎn)和后一節(jié)點(diǎn)連線與水平方向夾角角度值，βi∈［0，360)，如圖1所示。

圖1 角度參數(shù)計(jì)算示意

αi和βi的計(jì)算公式如下

其中，函數(shù)Fun(Math．Atan(k))表示計(jì)算斜率為k的向量與水平軸正方向夾角，夾角數(shù)值大于0、小于360。

2.2 節(jié)點(diǎn)匹配

對(duì)已知兩條等高線l1和l2所有節(jié)點(diǎn)，按照2.1小節(jié)對(duì)節(jié)點(diǎn)特征參數(shù)的定義以及公式(1)、(2)進(jìn)行計(jì)算，得到等高線l1和l2所有節(jié)點(diǎn)的特征參數(shù)集合

其中，m和n為等高線l1和l2節(jié)點(diǎn)個(gè)數(shù)。

以兩條等高線中節(jié)點(diǎn)個(gè)數(shù)多的那條為基準(zhǔn)，遍歷該條等高線上的每一個(gè)節(jié)點(diǎn)，假設(shè)當(dāng)前節(jié)點(diǎn)為Pc，然后在另外一條等高線中尋找距離該點(diǎn)最近的M個(gè)點(diǎn)集合R，計(jì)算點(diǎn)Pc與點(diǎn)集合R中每個(gè)點(diǎn)Pi(Pi∈R，i∈［0，M］)的特征距離 Di，計(jì)算公式為

其中，Pi∈R，i∈［0，M］。

設(shè)定距離閾值d，當(dāng)Di大于d時(shí)，進(jìn)行下一點(diǎn)的判斷;否則認(rèn)為Pi是Pc的變化趨勢一致點(diǎn)，記錄Pc和Pi兩點(diǎn)的坐標(biāo)，這兩點(diǎn)的連線即能將等高線分割成兩部分。

對(duì)基準(zhǔn)等高線上所有的節(jié)點(diǎn)匹配完成之后，就能得到一組匹配結(jié)果，將其圖形化表示，如圖2所示。

圖2 節(jié)點(diǎn)特征匹配

節(jié)點(diǎn)特征匹配結(jié)果可以將兩條等高線分割成多個(gè)匹配區(qū)域，但是通過圖形化結(jié)果發(fā)現(xiàn):匹配區(qū)域有交叉重疊部分和共點(diǎn)部分，需要按照一定的規(guī)則對(duì)匹配結(jié)果進(jìn)行篩除，即剔除圖2中相交和共點(diǎn)的線。篩除規(guī)則如下:

①依次遍歷匹配結(jié)果集合，取出當(dāng)前匹配結(jié)果，獲得匹配結(jié)果中兩點(diǎn)分別在兩條等高線中節(jié)點(diǎn)的序號(hào)index11、index21;

②取出前一匹配結(jié)果，獲得前一匹配結(jié)果中兩點(diǎn)分別在兩條等高線中節(jié)點(diǎn)的序號(hào)index10、index20;

③當(dāng)index21＜=index20，則剔除該匹配結(jié)果;否則繼續(xù)執(zhí)行①、②、③，直到所有的匹配結(jié)果遍歷完成，得到如圖3所示的結(jié)果。

圖3 特征匹配錯(cuò)誤剔除

2.3 輔助線構(gòu)造

相鄰兩組匹配結(jié)果可以將兩條等高線分割成變化趨勢一致的區(qū)域，依次遍歷篩選后的匹配結(jié)果，取當(dāng)前和下一組匹配結(jié)果，得到其分割兩條已知等高線的區(qū)域，進(jìn)行定數(shù)等分，然后依次匹配兩條等高線上的定數(shù)等分點(diǎn)，存入輔助線數(shù)組。圖形化的輔助線如圖4所示。

2.4 等高線生成

假設(shè)每根輔助線的起止點(diǎn)坐標(biāo)為(X1，Y1，Z1)和(X2，Y2，Z2)，可以得到待內(nèi)插的N條等高線的節(jié)點(diǎn)坐標(biāo)

以此遍歷每根輔助線，計(jì)算內(nèi)插點(diǎn)坐標(biāo)，將內(nèi)差點(diǎn)連接起來，即可得到等高線，如圖5所示。

圖4 構(gòu)造輔助線

圖5 通過構(gòu)造線內(nèi)插等高線

上述方法會(huì)造成等高線線上節(jié)點(diǎn)過多，道格拉斯－普克算法能夠準(zhǔn)確的刪除線段上小彎曲的節(jié)點(diǎn)，在整體上有效保持了線要素的形態(tài)特征［9］。因此，在得到等高線后，利用道格拉斯－普克算法進(jìn)行簡化，得到最終結(jié)果，如圖6所示。

圖6 道格拉斯－普克簡化后的等高線

3 實(shí)驗(yàn)及分析

利用等高線上當(dāng)前節(jié)點(diǎn)坐標(biāo)值，以及當(dāng)前節(jié)點(diǎn)和前后兩個(gè)相鄰節(jié)點(diǎn)向量與水平軸正方向的兩個(gè)夾角數(shù)值當(dāng)作描述當(dāng)前節(jié)點(diǎn)的特征參數(shù)，一方面能夠準(zhǔn)確地確定節(jié)點(diǎn)的空間位置;同時(shí)，兩個(gè)角度數(shù)值能夠很好的確定當(dāng)前節(jié)點(diǎn)在空間中的方位，便于在另一條等高線上尋找變化趨勢一致點(diǎn)。尋找到了變化趨勢一致點(diǎn)后，有利于將已知的等高線進(jìn)行區(qū)域分割。同一分割區(qū)域內(nèi)，兩條已知等高線變化趨勢是一致的。最后，在變化趨勢已知的分割區(qū)域內(nèi)，進(jìn)行輔助線構(gòu)造、等高線生成、節(jié)點(diǎn)簡化等操作，得到最后的結(jié)果。

利用本文提出的新方法進(jìn)行等高線的生成，并利用文獻(xiàn)［4］提出的算法進(jìn)行對(duì)比實(shí)驗(yàn)，對(duì)比結(jié)果如圖7所示，其中中間部分實(shí)現(xiàn)為本算法內(nèi)插結(jié)果，虛線為文獻(xiàn)［4］提出算法的結(jié)果。

圖7 內(nèi)插等高線結(jié)果對(duì)比

經(jīng)過對(duì)比分析可知，本文提出的新方法能夠很好的顧及局部地形，可以更好地解決等高線內(nèi)插的問題。

［1］ 龔有亮，何玉華，付子傲，等．一種實(shí)用的等高線內(nèi)插算法［J］．測繪學(xué)院學(xué)報(bào)，2002，19(1):36-37

［2］ 胡衛(wèi)明，吳兵，凌海濱．地圖等高線內(nèi)插算法［J］．計(jì)算機(jī)學(xué)報(bào)，2000，23(8):847-851

［3］ 楊長強(qiáng)，鄭永果，鄭作亞．一種等高線自動(dòng)加密算法［J］．信息技術(shù)與信息化，2006(2):56-58

［4］ 李洪玉，唐詩華，李景文，等．在CAD中實(shí)現(xiàn)等高線自動(dòng)內(nèi)插的一種方法［J］．海洋測繪，2004，24(3):52-54

［5］ 張渭軍，王文科，孔金玲，等．等雨量線的自動(dòng)繪制方法［J］．水土保持通報(bào)，2006，25(5):54-57

［6］ 曹穎，趙牡丹．DLG內(nèi)插等高線算法實(shí)現(xiàn)［J］．測繪科學(xué)，2007，32(2):67-68

［7］ 黃培之．等高線形狀變化規(guī)律的研究［J］．中國圖象圖形學(xué)報(bào)，2006，11(1):103-106

［8］ 姜波．一種等高線內(nèi)插的新思路［J］．科技通報(bào)，2010，26(5):780-781

［9］ 彭認(rèn)燦，董箭，鄭義東，等．垂距法與道格拉斯－普克法刪除冗余頂點(diǎn)效率的比較［J］．測繪通報(bào)，2010(3):66-67