利用AutoCAD二次開發實現地形圖高程點與等高線錯誤自動查找

劉 洪，黃軍勝

(1. 廣西有色勘察設計研究院，廣西 南寧 530031； 2. 廣西壯族自治區水利電力設計研究院，廣西 南寧 530023)

利用AutoCAD二次開發實現地形圖高程點與等高線錯誤自動查找

劉 洪1，黃軍勝2

(1. 廣西有色勘察設計研究院，廣西 南寧 530031； 2. 廣西壯族自治區水利電力設計研究院，廣西 南寧 530023)

論述了數字地形圖中產生高程點與等高線錯誤的可能原因，分析了高程點和等高線在地形圖中的空間位置關系及作為判斷條件的數學關系，論述了解決問題的辦法和判斷規則，展示了通過編制計算機程序實現找出其錯誤的核心代碼，此研究對于減輕質檢人員在檢查數字地形圖的高程點與等高線錯誤時的勞動強度和提高工作效率很有意義。

數字地形圖；等高線；高程點；點線不符；C#

因為科學技術的進步，在現階段地形圖的生產能力有了很大的提高，但由于人工智能仍不夠發達且很多環節仍需要人工干預，導致生產出的地形圖仍存在著高程點與等高線的錯誤。對于質檢人員而言，如果僅憑借著眼睛，要找出這些錯誤以提高地形圖的質量水平是非常困難的，勞動強度大不說，而且效率低下、極易出現遺漏。因此，找到一種高效率的辦法解決這個問題就顯得很有必要，本文正是據此分析產生錯誤的原因，找到解決的辦法，再利用C#+ObjectARX對AutoCAD二次開發編程實現查找數字地形圖中高程點與等高線的錯誤。

1 錯誤原因分析

1.1 高程點錯誤

目前大比例尺數字地形圖主要有3種生產方式：一是老圖矢量數字化，在識別老圖中的高程注記點時，有些不能正確識別，需要人工修改，容易輸入錯誤；二是利用航片經畸變糾正等工序生成正射影像得到數字線劃圖并在立體模型中取得高程屬性值，取得高程值時讀取錯誤，與實際繪制的等高線形成矛盾；三是全野外數據采集內業編輯形成數字地形圖，其注記點的高程直接由數據文件讀取生成，這種直接注記高程點的方式不易出現錯誤，但作業員在調整自動繪制的等高線時，可能為了圓滑美觀而人為調整等高線上的節點位置，從而導致高程注記點的高程數值與等高線形成矛盾。

由以上方式生產得到的數字地形圖，都有可能因圖面高程注記點的數量達不到規范要求而人工加注高程點，而人工加點由于作業人員的疏忽，判斷等高線錯誤導致加注高程點錯誤，這是最易產生錯誤的原因。

有時，作業員在發現高程注記點的高程數值與等高線相互矛盾時，為圖方便，僅將其注記數值修改成正確數值而沒有對高程點特性的坐標Z值進行相應修改，從而導致錯誤。

1.2 等高線錯誤

在調整不合理等高線時，作業員需刪除計算機繪制的部分不合理等高線，另行繪制補充。由于粗心，忘記修改等高線的高程屬性值，使其值為零或其他錯誤值，導致新繪制的等高線與已繪制好的等高線及高程點相互矛盾；或者又由于圖面元素太密集而在瀏覽時誤刪了某一條正確的等高線。

2 解決辦法探究

2.1 錯誤歸類

根據以上原因分析，可以看出數字地形圖中高程點與等高線產生錯誤的原因可以歸為如下4類。即點的高程特性坐標Z值與點的高程注記數值不符；等高線正確，但高程點不正確；高程點正確，但等高線不正確；等高線與高程點皆不正確。

2.2 解決辦法

2.2.1 點的高程特性坐標Z值與點的高程注記值不符

對數字地形圖的*.dwg格式文件進行遍歷，在高程點圖層中找到全部高程點，讀取其高程特性坐標Z值和高程注記數值(即高程點特性中的height屬性值)，并將兩者比較，若兩不相等，即為作業員修改過高程點的高程注記數值，用圓圈將其標示出來。

如圖1所示，高程點的高程特性坐標Z值原為“61.355 0”，點的高程注記數值原為“61.36”，被作業員人為修改為“59.36”，導致兩者不符，這在入庫時或許會引起新的錯誤。

圖1 高程點自身屬性數據矛盾

2.2.2 等高線正確，但高程點不正確

如果等高線與高程點都是正確的，則高程點的

高程注記值與其最近的一根等高線的高程值之差應該在一個等高距之內。

在AutoCAD中，每一個圖元都對應著一個唯一的ID，可以從高程點開始，人為構造一條射線，與等高線相交，提取這些交點和與之對應的等高線ID，然后按高程點與這些交點的遠近距離進行排序，由此可按順序逐條找到高程點周邊的不同等高線。

圖2所示為構造線與等高線相交的可能順序。

圖2 排序前等高線的順序

相交后，即可獲得各交點的坐標，據此可計算出各相交點與此高程點的距離，根據其距離值的大小，再次按近遠順序排列。

圖3所示，為構造線與各等高線相交點距離的近遠順序。

圖3 排序后等高線的順序

經過這樣排序之后，便可按順序找到相應的等高線，讀取其高程值。

圖4 高程點錯誤示例圖1

圖5 高程點錯誤示例圖2

如果H2

如果H2>H3，即H1>H3，則HA必小于H1，否則此高程點錯誤。

圖6 高程點錯誤示例圖3

2.2.3 高程點正確，但等高線不正確

(1) 等高線高程值為零這種情況存在于沿海一帶及海拔位于零左右的地方，在其他地方出現便是錯誤無疑。作業員在調整修改等高線時，對部分不合理的等高線進行了刪除操作，然后由手工繪制補充，而未對其賦予正確的高程值可能會出現零或輸入了其他錯誤值。這樣的等高線，在圖面上并不能看出不妥之處，但在建立立體模型生成DEM等產品的過程中，必定會出現與實地不符的形狀，從而導致錯誤，嚴重影響產品質量。

0米高程線在其線頭位置以圓圈標記，并注明此線高程為0。

(2) 假設A點沿構造線方向相交有3條以上多條等高線且等高線的高程值均不相等。

以此類推，可以判斷出第4條、第5條甚至更多的等高線。不過，此地形圖如果是按規范展繪高程點，其數量及位置涵蓋到了每個小區域，則不必每個點都判斷多條等高線，根據經驗，一般判斷3條就夠了。

圖7 等高線錯誤示例圖1

圖8 等高線錯誤示例圖2

2.2.4 等高線與高程點皆不正確

這種情況比較少，但如果采用以上規則進行判斷，則可以直接得出結論，無論是高程點還是等高線都將被圓圈標記，以示錯誤，需要人工干預，理清緣由。

2.2.5 特殊情況

(1) 若高程點A已在圖邊，則向圖邊構造線是無法判斷此高程點是否正確的，這種情況下，需要向另一方向構造線來與等高線相交，再用上述規則進行判斷。

如圖9所示，若右側為圖框，A點位于圖邊，向右的構造線已沒有等高線與其相交，這時，則需要向另一方向構造線段或射線來與等高線相交，如圖10所示。

圖9 構造線與等高線相交示例圖1

圖10 構造線與等高線相交示例圖2

(3) 若高程點A一側只有兩根等高線，如果兩根等高線高程不同，則可用上述規則判斷點是否錯誤或等高線是否錯誤，但如果兩根等高線高程相同，則不能確定高程點或等高線是否正確，這時，也需要向另一方向構造線來與等高線相交，再進行判斷。

(4) 若高程點A一側相交的第一條線為田埂線或陡坎(坡、崖)線，由于田的另一側地形與此高程點無太大聯系且同一條坎(坡、崖)的各點高度絕大部分都不相同，因此無法判斷，需要向另一方向構造線來與等高線相交，再進行判斷。因本程序的主要功能是查找高程點與相臨等高線的錯誤，如若點的各個方向的第一條交線都為田埂線或陡坎(坡、崖)線(如高程點在梯田中間)，則不再對此點判斷。

3 編程實現

3.1 編程環境配置

操作系統：MicrosoftWindows10中文旗艦版64位操作系統。

硬件開發環境：聯想筆記本昭陽E49A。

軟件開發環境：MicrosoftVisualStudio.NET2008+ObjectARX2010。

軟件運行環境：AutoCAD2010。

3.2 編程核心算法

選擇出數字地形圖中的田埂線和陡坎(坡、崖)線及等高線，構建選擇集；遍歷數字地形圖中的所有高程點，將高程點逐個向一側構造線與臨近的田埂線或陡坎(坡、崖)線或等高線相交，剔除田埂線和陡坎(坡、崖)線，比較高程點注記與Z值，高程點注記與鄰近等高線的高程值，以及臨近幾條等高線間的高程值，然后把錯誤的高程點在高程點位置用圓圈標記出來，把錯誤的等高線在等高線的交點處用另一半徑圓圈標記出來。部分核心代碼如下：

if(JdStruct.Length> 2)∥3條以上情況

{if(Math.Abs(JdStruct[2].Gc-JdStruct[1].Gc) ==DengGaoJu&&Math.Abs(JdStruct[1].Gc-JdStruct[0].Gc)!=DengGaoJu)

{ ∥第1條等高線有問題

CreateJianChaYuan(JdStruct[0].JD, 190,GScale,0.25);∥繪制標錯圓

CreateText(String.Format("線"),JdStruct[0].JD,GScale, 0.8, "檢查");

}

{ ∥第2條等高線有問題

CreateJianChaYuan(JdStruct[1].JD, 190, GScale, 0.25);∥繪制標錯圓

CreateText(String.Format("線"), JdStruct[1].JD, GScale, 0.8, "檢查");

}

if (Math.Abs(JdStruct[1].Gc - JdStruct[0].Gc) == DengGaoJu && Math.Abs(JdStruct[2].Gc - JdStruct[1].Gc)!= DengGaoJu)

{∥第3條等高線有問題

CreateJianChaYuan(JdStruct[2].JD, 190, GScale, 0.25);∥繪制標錯圓

CreateText(String.Format("線"), JdStruct[2].JD, GScale, 0.8, "檢查");

}

}

4 驗證效果

圖11所示為查找程序運行的開始界面，此時可在對話框中輸入當前數字地形圖的等高距，其默認值為1 m，使用時可根據實際情況填寫。

圖11 程序運行開始界面

圖12所示為查找程序運行的結束界面，點擊確定后即顯示檢查結果。

圖12 程序運行結束確認界面

圖13所示為程序運行后結果顯示。

圖13 錯誤查找結果

5 結 語

經過以上論述過程可以看出，本文的研究結果是正確的，可以利用前述判斷規則編制的程序來完成數字地形圖中高程點與等高線的錯誤自動查找，據此改正即可提高數字地形圖的質量，可以明顯減輕質檢人員在檢查數字地形圖的高程點與等高線錯誤時的勞動強度，提高工作效率，解決了質檢人員人工查找點線矛盾的難題。

[1] 洪濤，譚仁春，李寧. AutoCAD數據自動創建ArcGIS“地圖瓦片”服務的方法研究[J]. 測繪通報，2015(12): 85-87.

[2] 余代俊，蒲朝旭，朱逍賢. 基于AutoCAD二次開發的高斯投影轉等積割圓錐投影[J]. 測繪通報，2013(7): 71-73.

[3] 胡紹輝， 陳海英， 魏葉敏. AutoCAD 2004平臺下農村宅基地宗地圖批量生產軟件的研發[J]. 測繪通報，2016(9): 83-86.

[4] 詹必偉，高穎，翟翔，等. 基于AutoCAD的DLG自動化質量檢查探究[J]. 測繪通報，2014(10): 111-113.

[5] 郭沛沛，李成名，殷勇.建筑物合并的Delaunay三角網分類過濾法[J]. 測繪學報，2016，45(8): 1001-1007.

[6] 竇世卿，趙學勝，劉成軍，等. 河網線要素與DEM綜合的三維Douglas-Peucker算法[J]. 測繪學報，2016，45(4): 450-457.

[7] 李國輝，許文帥，龍毅，等. 面向等高線與河流沖突處理的多約束移位方法[J]. 測繪學報， 2014，43(11): 1204-1210.

[8] 馮向前. 使用ObjectARX實現自動檢查點線矛盾[J]. 城市勘測，2008(2): 115-118.

[9] ROBINSON S， NAGEL C. C#高級編程[M].3版.北京：清華大學出版社，2005.

[10] 李冠億. 深入淺出AutoCAD.NET二次開發[M].北京：中國建筑工業出版社，2012.

[11] 楊坤. 可視域分析算法的并行化與容錯研究[D].南京：南京師范大學，2014.

[12] 蔡偉明.一種基于AutoCAD平臺高程值檢查功能的開發[J].礦山測量，2010(3)：26-28.

[13] 董玉德，趙韓. CAD二次開發理論與技術[D].合肥：合肥工業大學出版社，2009.

[14] 陸貽豪.ObjectARX在城市地下排水管道上的開發[J].科技風，2016(14)：154-155.

[15] 秦洪現，崔惠嵐，孫劍.Autodesk 系列產品開發培訓教程[M]. 北京：化學工業出版社, 2008.

Using Secondary Development of AutoCAD to Automatically Find Errors of Elevation Points and Contour Lines in Digital Topographic Map

LIU Hong1,HUANG Junsheng2

(1. Guangxi Nonferrous Survey and Design Institute, Nanning 530031, China; 2. Guangxi Water and Power Design Institute, Nanning 530023, China)

This paper discusses the possible causes of the elevation point and contour errors in digital topographic map, analyzes the spatial location relationship and the mathematical relationship of the elevation point and contour in the topographic map, with the mathematical relationship working as a judgment condition. It discusses the solutions and judgment rules, and shows the realization of finding out its mistakes’core codes through computer program. This research is very meaningful to reduce the labor intensity of inspection personnel while checking elevation point and contour errors in the digital topographic map so as to improve the working efficiency.

digital topographic map； contour； elevation point； no conformation between the point and line； C#

劉洪，黃軍勝.利用AutoCAD二次開發實現地形圖高程點與等高線錯誤自動查找[J].測繪通報，2017(5)：100-104.

10.13474/j.cnki.11-2246.2017.0164.

2016-11-22；

2017-02-06

劉 洪(1974—)，男，碩士，高級工程師，從事大地測量與測繪工程工作。E-mail: 36682165@qq.com

P208

A

0494-0911(2017)05-0100-05