利用MapBasic7.0制作基于MapInfo的廣西森林資源連續清查控制點生成程序

黃甫權

摘要：在廣西連續清查控制點坐標的基礎上，重點闡述了利用MapBasic7.0制作的廣西連續清查控制點生成程序。

關鍵詞：MapBasic;MapInfo;廣西連續清查控制

廣西森林資源連續清查體系始建于1997年，是以掌握宏觀森林資源現狀與動態為目的，利用固定樣地為主進行定期復查的森林資源調查方法（一類調查），廣西已于2005年完成了第7次復查的任務[2]。

廣西森林資源連續清查體系，主要由按8km×6km布設的4948個地面固定樣地組成，因此對這4948個地面固定樣地進行計算機生成及存檔有著極大的重要意義，而利用MapInfo對這些固定樣地的GIS化管理提供了極大的方便。

1 研究方法

1.1 研究區域概況

廣西地處中國南部，位于北緯20°54′～26°23′，東經104°28′～122°04′。南臨北部灣，與海南省隔海相望，東連廣東省，東北接湖南省，西北靠貴州省，西與云南省接壤，西南與越南社會主義共和國毗鄰，具有沿海、沿江、沿邊的區位優勢。

廣西是中國南方的重要林區之一，據2005年廣西森林資源連續清查第7次復查統計結果顯示，廣西土地面積2376.0萬hm²，其中：林業用地面積1509.44萬hm²，占63.53%。

林業用地中，有林地面積924.35萬hm²，占61.24%;疏林地面積15.37萬hm²，占1.02%;灌木林地面積356.01萬hm²，占23.58%;未成林造林地49.48萬hm²，占3.28%;無立木林地面積48.01萬hm²，占3.18%;宜林地面積116.22萬hm²，占7.7%。森林（含國家特別規定灌木林）覆蓋率為52.71%.

全區活立木蓄積 51056.78萬m³。其中：有林地蓄積46875.18萬m³，占91.18%;疏林蓄積176.44萬m³，占0.35%;散生木蓄積3069.74萬m³，占6.01%;四旁樹蓄積935.42萬m³，占1.83%。

1.2 數據來源

從2005年度廣西第七次森林資源連續清查外業調查數據（主要分為樣地表和樣本表），提取縱橫坐標及樣點號等信息。

1.3 技術路線

通過對廣西第七次森林資源連續清查數據的提取，生成新的廣西連續清查控制點坐標表;其次，將控制點坐標表導入到MapInfo軟件中，生成Map-Info的*.Tab文件;最后，通過MapBasic編程，在MapInfo軟件中生成連續清查控制點，并與屬性數據相對接。

2 方法實現

2.1 確定連續清查控制點的度帶

廣西森林資源連續清查體系是以1：50000地形圖為基礎建立的，采用北京54坐標系，使得4948個連清控制點中存在有分度帶的問題。廣西經度范圍是東經104°28′～122°04′，其中：104°28′～108°部分為18°帶;108°～122°04′部分為19°帶。

18°帶與19°帶存在著中央經線的不同，它們的中央經線分別為105°和111°。因此，要分別生成18°帶控制點與19°帶控制點，最后合并成為全廣西連續清查控制點圖。

2.2 確定連續清查控制點的投影

大地基準：采用北京54坐標系，橢球體采用Krasovsky（1940）。投影方式：采用高斯一克呂格分帶投影。分帶方式：6°分帶。高程基準：采用1985國家高程基準。中央經線：105°（18°帶）、111°（19°帶）。

2.3 確定數據庫表結構

數據庫表主要是用于存儲連清控制點數據，而且對數據進行相應的規范。數據庫表結構主要決定于外業調查因子及統計表的需要。

2.4 編程代碼實現

2.4.1 生成18°度帶連續清查控制點。通過技術路線、投影、數據庫表的結構分析，在MapBasic7.0的編譯器中，編寫模塊代碼。在編寫模塊代碼之后，可以進行編譯，生成*.MBX可執行文件。

下面列出部分主要模塊代碼（以開頭的部分為代碼注譯部分）：

'生成控制點過程：

Set Map XY Units“m”CoordSys Earth Projection8，1001，“m”，105，0，1，500000，0Fetch First From t1

Do While Not EOT（t1）

x1=t1.x1

y1=t1.y1

no1_1=t1.no1

'生成樣地點

Create Point

Into Window map win_id

（y1，x1）

Symbol（49，red，20）

selectfrom t2 where x01=0

update Selection set x01

=x1，y01 =y1，no01=not 1

2.4.2 生成19°帶連續清查控制點。生成190帶連續清查控制點只需要改動生成18°帶程序，將Set MapXY Units“m”CoordSys Earth Projection 8，1001，“m”，105，0，1，500000，0中的“105”（中央經線）改為“111”，即能生成19°帶連清控制點。

2.5 合并2個度帶連續清查控制點

在MapInfo下，使用“向表中添加行”的辦法，將2個度帶連續清查控制點進行合并，即為最終的全區連續清查控制點分布圖。

3 結論

通過MapBasic編程生成的辦法，能夠很好地處理批量的控制點數據。同時，為今后廣西連續清查控制點的使用提供了方便，即廣西連續清查控制點可以通過MapInfo輸出到任何比例的地形圖上，而不是只能借助于原來的紙質圖紙。其次，通過修改程序中的中央經線，可以生成不同度帶的控制點，如需要生成3°帶（36帶），則將其改為“108”。因此，本程序具有很強的擴展性。（收稿：2019-07-23）

參考文獻：

[1]劉光.地理信息系統二次開發教程.

[2]廣西森林資源連續清查第七次復查操作細則.