翟高鵬,李文彬(1.河北省測繪產(chǎn)品質(zhì)量監(jiān)督檢驗站,河北石家莊 050031;.石家莊經(jīng)濟學(xué)院土地資源與城鄉(xiāng)規(guī)劃學(xué)院,河北石家莊 050031)

基于AutoLISP的地形圖數(shù)學(xué)精度檢測算法的設(shè)計與實現(xiàn)

翟高鵬1?,李文彬2
(1.河北省測繪產(chǎn)品質(zhì)量監(jiān)督檢驗站,河北石家莊 050031;
2.石家莊經(jīng)濟學(xué)院土地資源與城鄉(xiāng)規(guī)劃學(xué)院,河北石家莊 050031)

摘 要:介紹了測繪產(chǎn)品質(zhì)量監(jiān)督檢驗工作中對數(shù)字地形圖產(chǎn)品數(shù)學(xué)精度檢測的一種算法。該算法通過遍歷展繪的檢測點坐標(biāo),依次定位顯示檢測點視圖窗口,然后采用人機交互操作方式實現(xiàn)地形圖數(shù)學(xué)精度的快速檢測統(tǒng)計與評分。文章結(jié)合地形圖產(chǎn)品質(zhì)量檢驗項目,以LISP語言為例,對算法的基本設(shè)計思路進行了說明與實現(xiàn)。最后通過工程實例對程序的實用性和有效性進行了測試,達到了預(yù)期結(jié)果。

關(guān)鍵詞:AutoLISP;數(shù)字地形圖;數(shù)學(xué)精度;檢測

1　引 言

測繪產(chǎn)品質(zhì)量監(jiān)督檢驗中,地形圖數(shù)學(xué)精度檢測是一項十分重要的工作,檢測結(jié)果是產(chǎn)品質(zhì)量評定的主要技術(shù)指標(biāo)和重要參考依據(jù)[1]。測繪產(chǎn)品質(zhì)檢單位對測繪成果數(shù)學(xué)精度檢測的方法有很多種:目前流行使用的基于ArcGIS開發(fā)的4D Checker軟件,對于檢測入庫4D產(chǎn)品起到了很好的效果,但4D Checker軟件基于ArcGIS開發(fā)平臺,在檢測AutoAD平臺數(shù)字測繪成果前需要進行格式轉(zhuǎn)換,使得軟件的應(yīng)用受到一定限制;河南省測繪產(chǎn)品質(zhì)量監(jiān)督檢驗站開發(fā)的數(shù)字地形圖精度自動檢測系統(tǒng),在DWG格式地形圖精度檢測的應(yīng)用中發(fā)揮了很大作用。其他還有基于VB和CAD的地形圖精度檢測統(tǒng)計程序等等,這些研究工作都在一定程度上提高了地形圖數(shù)學(xué)精度檢測的效率,減輕了檢驗工作對人力因素的需求。本文所研究的是基于AutoLISP的CAD二次開發(fā)技術(shù),在大比例尺地形圖測繪中由于AutoCAD平臺應(yīng)用的廣泛性,使得本程序具有良好的通用性。該算法的具體實現(xiàn)也對提高地形圖產(chǎn)品質(zhì)檢效率具有一定的積極意義。

2　基本思路

本文通過對地形圖數(shù)學(xué)精度檢測的一般作業(yè)方式以及在精度檢測與統(tǒng)計過程中影響作業(yè)效率的關(guān)鍵步驟進行分析,有針對性的設(shè)計了展點范圍選定、檢測點定位顯示、粗差點提示和中誤差自動計算與質(zhì)量評分的算法。

2.1地形圖精度檢測作業(yè)方式

根據(jù)規(guī)范要求,進行地形圖的數(shù)學(xué)精度檢測,每幅圖需要分別采集20個～50個平面和高程檢測點[6],檢測點分布應(yīng)均勻覆蓋圖幅范圍。檢測點的采集由作業(yè)小組完成,一般使用RTK或全站儀對照地形圖在實地進行采集,采集過程中將采集的檢測點點位(高程注記檢測點為概略位置)標(biāo)繪到地形圖上。外業(yè)檢測點采集完成后,由內(nèi)業(yè)進行統(tǒng)計與評分。內(nèi)業(yè)是將采集的檢測點展繪到地形圖上,然后對照外業(yè)地形圖標(biāo)繪的檢測點位置進行一一計算與統(tǒng)計,最后基于統(tǒng)計得到的中誤差進行地形圖數(shù)學(xué)精度的質(zhì)量評分。

2.2展點范圍設(shè)定

作業(yè)小組采集的點位分布范圍比較大,一般是跨圖幅作業(yè)。由于作業(yè)習(xí)慣的不同,有些外業(yè)作業(yè)員提交的檢測點是幾個圖幅混合到一起的數(shù)據(jù),這對于內(nèi)業(yè)統(tǒng)計來說增加了必要的數(shù)據(jù)整理工作。

為此,本文設(shè)計了按圖幅范圍展繪檢測點的算法:在被檢測圖幅上選取圖框左下角(x1,y1)和右上角(x2,y2)角點坐標(biāo),以(x1＜x＜x2) and(y1＜y＜y2)作為篩選條件,依次讀取檢測點數(shù)據(jù)文件中的點,并將滿足條件的點展繪到圖幅上。檢測點展點調(diào)用的命令代碼為(command“circle”zb r) (command "text" tzb 0.8 0h),即調(diào)用CAD繪圓與添加文本命令。展繪的檢測點樣式是以點位為圓心,點位精度允許誤差為半徑的圓,檢測點高程值作為點注記。如圖1所示。

圖1　展點樣式

圖1中展繪檢測點示例所對應(yīng)的DXF組碼為: ((-1.＜圖元名:7ef76948＞) (0."CIRCLE") (330.＜圖元名:7ef71c20＞) (5." FC581") (100." AcD-bEntity") (67.0) (410."Model") (8."0") (62.230) (6."Continuous") (370.0) (100."AcDbCircle") (10 3.94676e+007 4.49896e+006684.026) (40 .0.5) (210 0.0 0.0 1.0))。

2.3檢測點定位顯示

影響地形圖數(shù)學(xué)精度統(tǒng)計效率的關(guān)鍵步驟是建立外業(yè)檢測點和與之相對應(yīng)圖上點相關(guān)信息的對應(yīng)關(guān)系。受計算機數(shù)字地形圖顯示范圍和人眼分辨能力所限,在地形圖中快速找出檢測點的識別范圍耗時費力。本文據(jù)此設(shè)計了檢測點定位顯示的算法,其基本思想如下:

(1)以展繪檢測點的DXF組碼“(0."CIRCLE") (40.0.5)”為篩選條件,建立以檢測點為成員的圖元選擇集。建立圖元選擇集的LISP代碼為(setq ss(ssget "X" lst))。

(2)遍歷選擇集圖元,獲取圖元坐標(biāo)。根據(jù)圖元坐標(biāo)并結(jié)合地形圖比例尺大小計算視窗角坐標(biāo)。設(shè)計出合理的視窗顯示范圍。定位顯示調(diào)用的CAD命令為(command ".zoom" " w" wpt1 wpt2),其中wpt1, wpt2為計算機視窗角坐標(biāo)。

(3)定位到檢測點后,通過人機交互實現(xiàn)檢測點與圖上點相關(guān)信息的匹配。主要包括平面檢測點與地物點、高程檢測點與高程注記點或高程值內(nèi)插點的匹配。

2.4 粗差記錄提示與評分

根據(jù)國家標(biāo)準(zhǔn)[6]的規(guī)定,合格地形圖產(chǎn)品數(shù)學(xué)精度檢測的粗差率不高于5%,因此統(tǒng)計過程中粗差的出現(xiàn)需要經(jīng)過謹(jǐn)慎判斷。地形圖數(shù)學(xué)精度評分依照國家標(biāo)準(zhǔn)規(guī)定的數(shù)學(xué)精度評分方法進行評分。

本文在數(shù)學(xué)精度檢測點點位信息計算完成后,設(shè)置了粗差提示。即根據(jù)匹配后計算出的差值與允許差值進行比較,當(dāng)差值超過允許值時判定為粗差并自動給出提示,根據(jù)實際情況判定是否對該檢測點位進行記錄。點位匹配完成后,通過程序?qū)崿F(xiàn)中誤差計算公式與數(shù)學(xué)精度評分方法所定義的分段函數(shù)公式[6],完成單幅地形圖成果中誤差計算和數(shù)學(xué)精度質(zhì)量評分。

3　算法實現(xiàn)與效果

該算法可以用多種編程語言實現(xiàn),如AutoLISP、ObjectARX、VBA等,本文采用AutoLISP編寫并在AutoCAD 2010下運行通過。目前作業(yè)單位普遍使用的數(shù)字成圖軟件,如南方CASS、北京威遠(yuǎn)圖易SV等軟件,均是在AutoCAD平臺下開發(fā),所以該程序具有良好的通用性。下面說明算法的主要步驟,其中算法思路相類似步驟受篇幅所限本文不再細(xì)述。最后對該算法的實用效果進行說明。

3.1算法主要步驟

(1)人機交互方式確定展點的范圍。交互操作只需提取范圍角點坐標(biāo),相應(yīng)代碼:(setq pt1(getpoint " n請點取展點范圍左下角點"))(setq pt2(getpoint " 請點取展點范圍右上角點"))

(2)篩選檢測點文件中圖幅范圍內(nèi)的點并以設(shè)定的圓形半徑進行展繪。主要程序代碼:

(if(and (＜x1 (atof x) x2) (＜y1 (atof y) y2)) (progn

(command "circle" zb r) ;繪點

(command "text" tzb 0.8 0 h) ;注記

) ;end progn

);end if

(3)高程精度檢測

①建立以展繪點圖元為單位的選擇集ss,(setq ss (ssget "X" lst));

②獲取選擇集圖元個數(shù)num,(setq num(sslength ss));

③依圖元個數(shù)num遍歷圖元,并據(jù)圖元坐標(biāo)定位顯示交互操作視圖窗口。

④人機交互,實現(xiàn)點位信息匹配:

(setq pt1 (getpoint " 點擊檢測高程點:"))

(if點擊空格,直接跳到第③步進行下一個檢測點視圖窗口顯示)

(progn

(高程值獲取方法選擇,“內(nèi)插高程值”或“點取高程點”) (setq gch(abs gch1)) ;計算高程差值。(if(＜mg gch);判斷高程誤差是否限差,(如果超限進行提示及相應(yīng)操作,否則計算并記錄到輸出文件)

);end if );end if

⑤調(diào)用預(yù)定函數(shù),進行中誤差統(tǒng)計與評分;

⑥完成統(tǒng)計成果輸出。

(4)平面精度檢測實現(xiàn)與高程精度檢測實現(xiàn)類似,在此不再細(xì)述。

3.2應(yīng)用示例

本文以某1∶500地形圖測繪產(chǎn)品質(zhì)量檢驗項目為例說明算法及程序的有效性。根據(jù)GB/ T 24356-2009測繪成果質(zhì)量查檢與驗收,外業(yè)采集的檢測點個數(shù)介于20個～50之間均可。示例中實際采集到的檢測點個數(shù)為高程和平面檢測點各24個,其中包括公共點。如圖2所示,展繪的檢測點顯示為紅色圓形并以高程作注記。

圖2　展繪點示意圖

展繪點完成后分別進行地形圖高程與平面精度的檢測與統(tǒng)計。在計算機命令行顯示操作提示,通過人機交互實現(xiàn)點位信息的匹配。高程和平面精度檢測的交互界面分別如圖3、圖4所示。

圖3　高程精度檢測人機交互界面

圖4　平面精度檢測人機交互界面

通過人機交互操作將點位信息匹配完成后,分別輸出高程和平面精度統(tǒng)計的評分成果。輸出格式為TXT,如圖5、圖6所示。

4　結(jié) 語

本文針對測繪產(chǎn)品質(zhì)量監(jiān)督檢驗工作中的地形圖數(shù)學(xué)精度檢測與統(tǒng)計,基于AutoCAD二次開發(fā)技術(shù)的LISP計算機語言,設(shè)計并實現(xiàn)了地形圖數(shù)學(xué)精度檢測與統(tǒng)計的算法,最后將程序應(yīng)用于地形圖質(zhì)檢項目。

圖5　高程精度檢測與統(tǒng)計結(jié)果輸出

圖6　平面精度檢測與統(tǒng)計結(jié)果輸出

通過實現(xiàn)點位信息人工匹配與粗差提示等算法,有效地降低了檢測統(tǒng)計工作中出錯的概率;檢測點定位顯示與自動評分,提高了地形圖數(shù)學(xué)精度檢測與統(tǒng)計的效率。

參考文獻

[1] 鄧思勝,章亞平,呂冰.基于VB和CAD的地形圖精度檢測統(tǒng)計程序[J].地理空間信息,2008(5),124～126.

[2] 蔡建德,張福利.數(shù)字地形圖精度自動檢查統(tǒng)計系統(tǒng)的開發(fā)和應(yīng)用[J].測繪通報,2008(6),41～44.

[3] 胡友明.AutoCAD VBA在數(shù)字地形圖精度檢測統(tǒng)計中的應(yīng)用[J].城市勘測,2009(4),123～125.

[4] 建國.基于AutoCAD平臺運用VBA語言設(shè)計編寫地形圖高程精度檢測軟件[J].工程勘察,2011(11),85～89.

[5] 李學(xué)智,方戈亮,孫力紅.Visual LISP程序設(shè)計(第二版) [M].北京:清華大學(xué)出版社,2010,3～53.

[6] GB/ T 24356-2009.測繪成果質(zhì)量檢查與驗收[S].

Design and Implementation of an Algorithm for Topographic Map Mathematical Precision Detection Based on AutoLISP

Zhai Gaopeng1,Li Wenbin2
(1.The Station of Quality Supervision and Inspection about Surveying and Mapping Products,Shijiazhuang 050031,China; 2.College of land resources and rural planning,Shijiazhuang University of Economics,Shijiazhuang 050031,China)

Abstract:Introduce an algorithm for mathematical precision detection of Digital Topographic Map products in the surveying and mapping product quality supervision and inspection.Mathematical precision of topographic map shall be efficiently detected,statistic analyzed and scored,which is by traversing the detecting points’coordinate,displaying the view window of their position in turn,and then operating by human-computer interaction mode till completed.Based on topographic map quality inspection items,taking LISP as an illustration,the basic idea of the algorithm was illustrated and implemented.Finally through an example of test about the practicability and validity of application,which achieved the desired effect.

Key words:AutoLISP;digital topographic map;mathematical precision;detection

文章編號:1672-8262(2015)01-37-04中圖分類號:P208.1

文獻標(biāo)識碼:A

收稿日期:?2014—06—29

作者簡介:翟高鵬(1985—),男,碩士,助理工程師,主要從事工程測量與測繪產(chǎn)品質(zhì)檢工作。