利用MapBasic7.0制作基于MapInfo的廣西森林資源連續(xù)清查控制點(diǎn)生成程序

黃甫權(quán)

摘要：在廣西連續(xù)清查控制點(diǎn)坐標(biāo)的基礎(chǔ)上，重點(diǎn)闡述了利用MapBasic7.0制作的廣西連續(xù)清查控制點(diǎn)生成程序。

關(guān)鍵詞：MapBasic;MapInfo;廣西連續(xù)清查控制

廣西森林資源連續(xù)清查體系始建于1997年，是以掌握宏觀森林資源現(xiàn)狀與動(dòng)態(tài)為目的，利用固定樣地為主進(jìn)行定期復(fù)查的森林資源調(diào)查方法（一類(lèi)調(diào)查），廣西已于2005年完成了第7次復(fù)查的任務(wù)[2]。

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

1 研究方法

1.1 研究區(qū)域概況

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

廣西是中國(guó)南方的重要林區(qū)之一，據(jù)2005年廣西森林資源連續(xù)清查第7次復(fù)查統(tǒng)計(jì)結(jié)果顯示，廣西土地面積2376.0萬(wàn)hm²，其中：林業(yè)用地面積1509.44萬(wàn)hm²，占63.53%。

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

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

1.2 數(shù)據(jù)來(lái)源

從2005年度廣西第七次森林資源連續(xù)清查外業(yè)調(diào)查數(shù)據(jù)（主要分為樣地表和樣本表），提取縱橫坐標(biāo)及樣點(diǎn)號(hào)等信息。

1.3 技術(shù)路線

通過(guò)對(duì)廣西第七次森林資源連續(xù)清查數(shù)據(jù)的提取，生成新的廣西連續(xù)清查控制點(diǎn)坐標(biāo)表;其次，將控制點(diǎn)坐標(biāo)表導(dǎo)入到MapInfo軟件中，生成Map-Info的*.Tab文件;最后，通過(guò)MapBasic編程，在MapInfo軟件中生成連續(xù)清查控制點(diǎn)，并與屬性數(shù)據(jù)相對(duì)接。

2 方法實(shí)現(xiàn)

2.1 確定連續(xù)清查控制點(diǎn)的度帶

廣西森林資源連續(xù)清查體系是以1：50000地形圖為基礎(chǔ)建立的，采用北京54坐標(biāo)系，使得4948個(gè)連清控制點(diǎn)中存在有分度帶的問(wèn)題。廣西經(jīng)度范圍是東經(jīng)104°28′～122°04′，其中：104°28′～108°部分為18°帶;108°～122°04′部分為19°帶。

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

2.2 確定連續(xù)清查控制點(diǎn)的投影

大地基準(zhǔn)：采用北京54坐標(biāo)系，橢球體采用Krasovsky（1940）。投影方式：采用高斯一克呂格分帶投影。分帶方式：6°分帶。高程基準(zhǔn)：采用1985國(guó)家高程基準(zhǔn)。中央經(jīng)線：105°（18°帶）、111°（19°帶）。

2.3 確定數(shù)據(jù)庫(kù)表結(jié)構(gòu)

數(shù)據(jù)庫(kù)表主要是用于存儲(chǔ)連清控制點(diǎn)數(shù)據(jù)，而且對(duì)數(shù)據(jù)進(jìn)行相應(yīng)的規(guī)范。數(shù)據(jù)庫(kù)表結(jié)構(gòu)主要決定于外業(yè)調(diào)查因子及統(tǒng)計(jì)表的需要。

2.4 編程代碼實(shí)現(xiàn)

2.4.1 生成18°度帶連續(xù)清查控制點(diǎn)。通過(guò)技術(shù)路線、投影、數(shù)據(jù)庫(kù)表的結(jié)構(gòu)分析，在MapBasic7.0的編譯器中，編寫(xiě)模塊代碼。在編寫(xiě)模塊代碼之后，可以進(jìn)行編譯，生成*.MBX可執(zhí)行文件。

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

'生成控制點(diǎn)過(guò)程：

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

'生成樣地點(diǎn)

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°帶連續(xù)清查控制點(diǎn)。生成190帶連續(xù)清查控制點(diǎn)只需要改動(dòng)生成18°帶程序，將Set MapXY Units“m”CoordSys Earth Projection 8，1001，“m”，105，0，1，500000，0中的“105”（中央經(jīng)線）改為“111”，即能生成19°帶連清控制點(diǎn)。

2.5 合并2個(gè)度帶連續(xù)清查控制點(diǎn)

在MapInfo下，使用“向表中添加行”的辦法，將2個(gè)度帶連續(xù)清查控制點(diǎn)進(jìn)行合并，即為最終的全區(qū)連續(xù)清查控制點(diǎn)分布圖。

3 結(jié)論

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

參考文獻(xiàn)：

[1]劉光.地理信息系統(tǒng)二次開(kāi)發(fā)教程.

[2]廣西森林資源連續(xù)清查第七次復(fù)查操作細(xì)則.