vlisp地塊自動分割算法在土地確權(quán)中的應(yīng)用研究

葛俊濤，高彥濤，賈復(fù)生，耿艷輝

（1.河南省地質(zhì)礦產(chǎn)勘查開發(fā)局測繪地理信息院，河南 鄭州 450006；2.河南省天空地遙感智能監(jiān)測工程技術(shù)研究中心，河南 鄭州 450006；3.河南省自然資源天空地遙感智能監(jiān)測研究科技創(chuàng)新中心，河南 鄭州 450006）

土地確權(quán)工作中，很多針對確權(quán)工作的建庫軟件應(yīng)運而生，但多數(shù)測量工作人員還是習(xí)慣用CAD或者南方CASS軟件，根據(jù)數(shù)字測圖的方法繪制地塊和填寫地塊屬性。在CAD中將外業(yè)采集獲得的數(shù)據(jù)進行預(yù)處理后，導(dǎo)入國土資源等其他農(nóng)經(jīng)權(quán)建庫軟件。然而數(shù)字測圖正常流程在清晰度高的基礎(chǔ)上矢量化繪制出地塊，鍵入承包方姓名和合同面積工作時，存在大量重復(fù)繁瑣的步驟，耗費大量的人力、物力、財力，在地塊勾繪、承包方姓名以及合同面積填寫方面的工作效率也有待提高[1-8]。

1 土地確權(quán)工作數(shù)學(xué)基礎(chǔ)與調(diào)查任務(wù)

1.1 土地確權(quán)內(nèi)容

根據(jù)GB/T1.1-2009要求，農(nóng)村土地承包經(jīng)營權(quán)確權(quán)采用2000國家大地坐標系，1985國家高程基準，調(diào)查比例尺為1:2 000選擇高斯-克呂格投影，采用標準3°分帶平面直角坐標系統(tǒng)。采用航測圖解和R TKGPS實測相結(jié)合的方式，輔以農(nóng)戶調(diào)研，對農(nóng)村土地各類權(quán)屬界線全部現(xiàn)場調(diào)查指界。由內(nèi)業(yè)人員在GIS平臺上繪制權(quán)屬線，標注地塊編號和權(quán)屬，記錄土地所有權(quán)權(quán)利人、界線的類型、位置、范圍、歷史沿革及權(quán)屬爭議等情況，調(diào)查清楚承包方土地承包經(jīng)營權(quán)權(quán)屬信息等問題，查清承包地塊名稱、面積、四至、空間位置以及土地用途等信息[5]，并在重要界址點埋設(shè)界樁，將權(quán)屬矢量線劃數(shù)據(jù)疊加到工作底圖上，形成地籍正射影像圖；再將權(quán)屬信息導(dǎo)入地理數(shù)據(jù)庫，最終完成地塊地籍權(quán)屬的確定，有效提升了地籍登記管理的效率，為鞏固我國農(nóng)村土地基本經(jīng)營制度提供了強有力的技術(shù)支持[9]。

1.2 數(shù)字測圖技術(shù)在土地確權(quán)中的應(yīng)用

數(shù)字測圖技術(shù)分為兩個步驟，第一步，外業(yè)數(shù)據(jù)采集：主要利用全站儀采集數(shù)據(jù)，電子手簿記錄，同時人工繪制標注測點點號的草圖，到室內(nèi)將測量數(shù)據(jù)直接由記錄器傳輸?shù)接嬎銠C中，再由人工按草圖編輯圖形文件[10]。由于外業(yè)基本都是空曠地域的農(nóng)田開展，并無高大樹木遮擋屏蔽，可用GPS-RTK進行外業(yè)采集代替全站儀，工作效率髙、定位精度高、全天候作業(yè)。只要滿足基本工作條件，在一定的作業(yè)半徑范圍內(nèi)的平面精度和高程精度都能達到厘米級[1]。在保證外業(yè)采集的點位數(shù)據(jù)精度達到農(nóng)業(yè)部規(guī)定的誤差要求（如表1所示）情況下，才能進行下一步工作。

表1 農(nóng)經(jīng)權(quán)實測法點位精度

普通地區(qū)的精度達到二級以上；特殊地區(qū)精度達到三級以上。第二步，內(nèi)業(yè)數(shù)據(jù)處理。將外業(yè)采集的每一個地塊的界址點坐標，通過“展點號”在南方CASS繪圖系統(tǒng)中顯示。將單個地塊的界址點勾繪成閉合的圖形，在地塊質(zhì)心位置標注承包方姓和合同面積（如圖1所示）。

圖1 地塊繪制示意圖

地塊繪制完成后，為了能夠?qū)肫渌r(nóng)經(jīng)權(quán)建庫軟件，還需將數(shù)據(jù)進一步處理，即界址點放到JZD圖層；圖形放置到JZX圖層；承包方姓名放置到CBFXM圖層；合同面積放置到HTMJ圖層；圖元分層放置以便轉(zhuǎn)入其他建庫軟件時查找（如表2）。

表2 數(shù)據(jù)圖元分層

在實際的繪圖過程中，繪制界址線和鍵入地塊屬性信息是一個重復(fù)繁瑣的過程，本文將主要研究運用VLISP進行程序設(shè)計，開發(fā)出CAD插件，從而實現(xiàn)高效繪制承包地塊和鍵入承包方姓名以及合同面積，以解決重復(fù)繁瑣的工作步驟。

2 vlisp地塊自動分割算法的程序設(shè)計

2.1 LISP語言

AutoCAD是使用非常廣泛的計算機輔助設(shè)計繪圖軟件，但大量的重復(fù)性工作又使得工程設(shè)計制圖繁瑣。用Visual LISP則可以很好地解決這一問題，在編程過程中，CAD系統(tǒng)提供軟件的界面、環(huán)境、核心算法、數(shù)據(jù)庫及設(shè)計需要相關(guān)的、較專業(yè)的支持軟件，而在工程測量中CAD計算導(dǎo)線，不能達到導(dǎo)線、展點、圖形一次性完成，不利于需用其他軟件來輔助完成，這樣不僅費時費力，而且也增加了出錯的可能性。Visual LISP可以直接調(diào)用執(zhí)行AutoCAD所有命令，是強化AutoCAD的程序語言[12]。

2.2 程序流程

本文的程序算法是根據(jù)多邊形固定一邊的方法設(shè)計，執(zhí)行本文程序時首先將要處理的大地塊輪廓勾繪完成，然后在程序中導(dǎo)入承包方姓名和合同面積，程序會根據(jù)讀取的大地塊面積進行差額分配，將面積差額按比例分配到各承包方面積中，將集體大地塊分割成承包戶二輪承包的承包地塊，最后導(dǎo)入承包方的姓名和合同面積標注在分割完成的地塊的質(zhì)心處，并自動將數(shù)據(jù)分圖層放置，如圖2所示。

圖2 程序流程圖

如圖3所示，將A1D1作為固定的地塊的一邊，讀取要分割地塊長邊和短邊夾角，根據(jù)要分出的合同面積Sh，運用梯形的面積計算公式和三角函數(shù)進行推導(dǎo)可得分割寬度H，計算公式如下：

圖3 地塊分割示意圖

捕捉固定邊的中點為起點，根據(jù)求得的分割寬度H，繪制出分割線，運用“boundary”命令進行構(gòu)面，分割出二輪承包地塊，將多余的構(gòu)面線刪減，得到封閉的承包方地塊，將承包方姓名和合同面積標注在封閉面域的質(zhì)心處，并分圖層放置。

3 vlisp地塊自動分割算法的程序?qū)崿F(xiàn)

3.1 實例驗證

將目標村的影像底圖導(dǎo)入到CAD中，將界線清晰的集體大地塊勾繪出來，如圖4所示，將項目前期搜集的承包方信息表整理后在Excel中編輯，保存成csv格式，或保存成txt格式，數(shù)據(jù)用逗號或者空格分開，導(dǎo)入小程序，實驗數(shù)據(jù)如表3所示。

表3 測試數(shù)據(jù)

圖4 待分割地塊

繪制地塊草圖的目的是為了將集體大地塊分割成承包戶二輪承包的地塊。

標注承包方姓名和合同面積時，該程序可直接允許Excel表格中編輯完成的承包方姓名和合同面積粘貼在對話框中，執(zhí)行后承包姓名和合同面積標注在承包戶小地塊的中心位置。

點擊執(zhí)行，如圖5所示。根據(jù)南方CASS命令行的提示進行操作，即可一鍵繪出地塊草圖，如圖6所示。

圖5 程序界面

圖6 分割結(jié)果

3.2 野外數(shù)據(jù)采集和界址點矯正

運用小程序?qū)⒉輬D勾繪完成后，進行數(shù)字測圖的第一階段和第二階段，將用GPS-RTK采集獲得上述承包方地塊的界址點勾繪到草圖中，依次檢查草圖地塊界址點和實測點的偏差，如發(fā)現(xiàn)偏差，僅需將對象捕捉開啟，拖動小地塊公共節(jié)點與實測點重合即可。

4 結(jié)語

繪制地塊及其屬性工作是農(nóng)經(jīng)權(quán)內(nèi)業(yè)數(shù)據(jù)處理的繁瑣任務(wù)之一，承包地塊及其屬性信息繪制具有大批量、同步驟等特性，針對這種情況，利用vlisp開發(fā)CAD插件，用于農(nóng)經(jīng)權(quán)批量繪制地塊及權(quán)屬信息賦值。這種方法已經(jīng)在河南省周口市淮陽縣大連鄉(xiāng)的農(nóng)村土地承包經(jīng)營權(quán)登記頒證工作中使用，結(jié)果表明：自動分割插件在農(nóng)經(jīng)權(quán)工作中提高效率、節(jié)省了大量的人力、物力，在地塊勾繪、承包方姓名以及合同面積填寫方面的工作效率有了顯著的提高。