章明 李玉華 曹蕾 鐘菲

摘 要 在GIS空間分析功能的支持下，利用數(shù)字高程模型（Digital Elevation Model，DEM）提取陽(yáng)江農(nóng)場(chǎng)高程、坡度和坡向數(shù)據(jù)，計(jì)算橡膠林、農(nóng)用地及未利用地等3種土地利用類型的分級(jí)分布指數(shù)，準(zhǔn)確描述農(nóng)場(chǎng)主要用地類型的空間分布頻率，總結(jié)分布趨勢(shì)，為進(jìn)一步制定熱作區(qū)農(nóng)場(chǎng)發(fā)展規(guī)劃提供參考。

關(guān)鍵詞 DEM ；熱作區(qū) ；土地利用現(xiàn)狀

分類號(hào) F301.2

土地利用類型的空間分布格局主要受到地形條件的制約[1]，地形條件又直接或間接的作用于光、熱、水、土等自然因素的再分配[2-3]，且呈現(xiàn)出一定的規(guī)律性。有學(xué)者以DEM為基礎(chǔ)探討了地形因子與土地利用的關(guān)系[4]，利用DEM揭示土地利用類型在高程、坡度和坡向上的分布[5-6]。利用指數(shù)模型分析區(qū)域地形條件對(duì)土地利用分布的影響和評(píng)價(jià)的研究方式也受到越來(lái)越廣泛的關(guān)注[7-9]。本文使用的高分辨率DEM（5 m），通過(guò)運(yùn)用地理信息空間分析方法得出陽(yáng)江農(nóng)場(chǎng)土地利用類型在高程、坡度以及坡向的分級(jí)分布指數(shù)，準(zhǔn)確描述了陽(yáng)江農(nóng)場(chǎng)膠園空間分布的變化特征，從而為熱帶作物地區(qū)土地資源利用優(yōu)化布置和熱作經(jīng)濟(jì)的可持續(xù)發(fā)展戰(zhàn)略設(shè)計(jì)提供重要的依據(jù)。

1 研究區(qū)概況及研究方法

1.1 研究區(qū)概況

海南省國(guó)營(yíng)陽(yáng)江農(nóng)場(chǎng)位于中國(guó)唯一具有熱帶海洋氣候的島嶼——海南島，東經(jīng)109°44′01″～109°55′20″，北緯19°17′02″～19°26′24″，年平均氣溫在25℃左右，該農(nóng)場(chǎng)是海南農(nóng)墾橡膠產(chǎn)膠大場(chǎng)之一，土地面積1.19萬(wàn)hm2，0.513萬(wàn)hm2橡膠園，年產(chǎn)干膠5 000 t。區(qū)域內(nèi)地形起伏不大，最高海拔465 m，最低海拔85 m，相對(duì)高差345 m，為淺丘陵地帶，植被覆蓋達(dá)到85%以上。

1.2 研究方法

1.2.1 數(shù)據(jù)源與數(shù)據(jù)處理

本研究采用海南省獨(dú)立坐標(biāo)系1∶1萬(wàn)數(shù)字地形圖，陽(yáng)江農(nóng)場(chǎng)1∶1萬(wàn)土地利用現(xiàn)狀圖，并通過(guò)設(shè)置海南省獨(dú)立坐標(biāo)系轉(zhuǎn)換參數(shù)，將以上數(shù)據(jù)從獨(dú)立坐標(biāo)系坐標(biāo)轉(zhuǎn)換到北京1954坐標(biāo)系和1985高程系，以實(shí)現(xiàn)空間數(shù)據(jù)的疊加分析[10-11]。

（1）數(shù)字高程模型：將陽(yáng)江農(nóng)場(chǎng)矢量等高線和離散高程點(diǎn)生成不規(guī)則三角網(wǎng)（Triangulated Irregular Network，TIN），后經(jīng)TIN柵格化獲得5m分辨率DEM[12]，由所得DEM提取高程、坡度與坡向等作為實(shí)驗(yàn)數(shù)據(jù)。

（2）土地利用現(xiàn)狀數(shù)據(jù)：將陽(yáng)江農(nóng)場(chǎng)土地利用現(xiàn)狀數(shù)據(jù)柵格化，分辨率設(shè)定為與DEM相同的5m。由于陽(yáng)江農(nóng)場(chǎng)主要以橡膠種植生產(chǎn)為主，其它如旱地、水田、菜地及其他熱作園所占面積較少，故將其歸并到為一種土地利用類型類——農(nóng)用地，將橡膠園和面積較大的未利用地作為獨(dú)立地類。

1.2.2 基于DEM的地形因子分級(jí)

高程分級(jí)：陽(yáng)江農(nóng)場(chǎng)的海拔高程值處于85～465 m，地形起伏較為均勻，因此參考等間距分級(jí)法將陽(yáng)江農(nóng)場(chǎng)DEM重分類為8個(gè)高程級(jí)，85～100 m為第1高程級(jí)，400～465 m為第8高程級(jí)，100～300 m按50 m等間隔分類為2至7高程級(jí)。

坡度分級(jí)：按照水土保持通用分級(jí)標(biāo)準(zhǔn)進(jìn)行分級(jí)，共分13級(jí)，依次為0°～3°，3°～5°，5°～10°，10°～15°，15°～20°，20°～25°，25°～30°，30°～35°，35°～40°，40°～45°，45°～50°，50°～55°，>55°。

坡向：該定義為坡面法線水平面投影與正北方向的夾角，通常將正北方向定為0°，按順時(shí)針?lè)较蛐D(zhuǎn)來(lái)劃分坡向，0°為平坡，0°～50°和320°～360°為陰坡，230°～320°為半陰坡，140°～230°為陽(yáng)坡，50°～140°為半陽(yáng)坡。

1.2.3 分級(jí)分布指數(shù)

本文采用分級(jí)分布指數(shù)[13]來(lái)描述不同土地利用類型在高程、坡度、坡向上的分布特征。分級(jí)分布指數(shù)是一個(gè)標(biāo)準(zhǔn)化、無(wú)量綱指數(shù)，該指數(shù)的引入將消除面積差異帶來(lái)的量綱影響[14]，計(jì)算公式如下。

P=（）·（）

P代表分布指數(shù)，e代表地形因子， 代表e地形因子特定等級(jí)下的i地類的面積，Si是i地類的面積，Se是整個(gè)區(qū)域地形因子特定等級(jí)下的總面積，S是整個(gè)區(qū)域面積。P值越大，說(shuō)明某土地利用類型出現(xiàn)頻率越高。

1.2.4 研究流程圖

利用DEM可以自動(dòng)提取地形因子，并與土地利用數(shù)據(jù)進(jìn)行空間分析，具有快速、科學(xué)、客觀等特點(diǎn)，大大方便了對(duì)土地利用的地形格局特征進(jìn)行的研究[15]。本研究基本流程見(jiàn)圖1。

2 結(jié)果與分析

2.1 分級(jí)分布指數(shù)處理結(jié)果

利用地理信息技術(shù)空間分析（Spatial Analyst）方法，將橡膠園、農(nóng)用地和未利用地柵格數(shù)據(jù)與所提取的地形因子進(jìn)行疊加分析[16]。

2.2 高程分布指數(shù)

從圖2中可以看出，膠園、農(nóng)用地、未利用地在不同高程級(jí)別中的分布趨勢(shì)。在第1高程等級(jí)（<100 m）上的分布指數(shù)均小于1。在第2、3級(jí)（100～200 m）上膠園和農(nóng)用地分布最多，即在這兩個(gè)高程級(jí)上，農(nóng)用地和膠園的出現(xiàn)頻率是最大的，但隨著高程的增加，分布逐漸下降。由于作為農(nóng)場(chǎng)主產(chǎn)業(yè)的橡膠和農(nóng)用產(chǎn)品主要分布在海拔低的地區(qū)，故絕大多數(shù)土地都得到了很好的開(kāi)發(fā)和利用。因此，在這兩個(gè)區(qū)域未利用地的分布就很小，而多數(shù)也是屬于那種利用投入過(guò)大，且不利于耕作的地塊。在第4級(jí)高程區(qū)，膠園和農(nóng)用地的分布已經(jīng)降到了0.5的水平，而未利用地開(kāi)始占據(jù)主導(dǎo)地位，膠園和農(nóng)用地的分布降到了最低點(diǎn)，且在6、7、8三個(gè)等級(jí)中，沒(méi)有膠園的分布，農(nóng)用地利用也從第7等級(jí)開(kāi)始，沒(méi)有分布。

2.3 膠園分布隨坡度的變化

從圖3可以看出，膠園主要分布于1～4坡度級(jí)，農(nóng)用地主要分布在1～3級(jí)，坡度分布指數(shù)均接近和超過(guò)1。隨著坡度的增加，膠園分布逐漸增加，農(nóng)用地分布逐漸減小，未利用地呈現(xiàn)出明顯增多的趨勢(shì)。由于不同坡度坡面，物質(zhì)交流與能量轉(zhuǎn)換的方式與效率有著明顯的差異，并且直接影響到土壤的母質(zhì)組成、土層厚度、土壤水分、無(wú)機(jī)鹽淋失和土壤酸堿度等土壤的基本屬性，坡度不同，也產(chǎn)生了坡面接受太陽(yáng)輻射量的差異，這些都將進(jìn)一步影響到該坡面的植被類型與分布頻度。

2.4 膠園分布隨坡向的變化

坡面的朝向，決定了該坡面接受太陽(yáng)輻射以及地面水分的量值與變化程度，因而在很大程度上影響植被的類型與分布，在陰陽(yáng)坡之間出現(xiàn)明顯的差異。但是由于陽(yáng)江農(nóng)場(chǎng)地處熱帶，主要以低丘為主，而日照和降雨充足，這點(diǎn)從圖4也反映出來(lái)。坡向?qū)δz園及農(nóng)用地分布影響不大，即使在陰坡中分布指數(shù)達(dá)到1，且在半陽(yáng)坡出現(xiàn)頻度更大。說(shuō)明橡膠樹(shù)這種高度在20～40 m的高大的喬木樹(shù)種，對(duì)坡向的因子并不敏感，可以利用這一特點(diǎn)，排除該因子在新建橡膠園規(guī)劃中的權(quán)重。

3 結(jié)論與討論

（1）本文以DEM和土地利用數(shù)據(jù)為主要信息源，在地理信息空間分析技術(shù)的支持下，較準(zhǔn)確的評(píng)價(jià)了熱作區(qū)農(nóng)場(chǎng)空間格局。熱作區(qū)有效土地利用類型分布主要受高程和坡度的影響，基本可不考慮坡向影響；陽(yáng)江農(nóng)場(chǎng)膠園和農(nóng)用地主要集中于300 m高程以下，坡度小于25°，分布較為均勻，且利用程度較好；隨著高程的增加和坡度的提高，未利用地分布頻度逐漸增加。

（2）結(jié)合對(duì)陽(yáng)江農(nóng)場(chǎng)的實(shí)際調(diào)查表明，陽(yáng)江農(nóng)場(chǎng)土地利用類型空間分布特征結(jié)果的可靠性和準(zhǔn)確性較高，驗(yàn)證了橡膠樹(shù)這種高度在20～40 m的高大的喬木樹(shù)種對(duì)坡向的因子影響并不敏感，可以利用這一特點(diǎn)，排除該因子在新建橡膠園規(guī)劃中的布局權(quán)重。

（3）通過(guò)對(duì)陽(yáng)江農(nóng)場(chǎng)土地利用空間分布特征研究，可以看出，人工開(kāi)墾農(nóng)場(chǎng)受地形因子中高程和坡度的限制較為明顯，應(yīng)在現(xiàn)有土地利用空間格局的基礎(chǔ)上，進(jìn)一步優(yōu)化土地利用結(jié)構(gòu)，實(shí)現(xiàn)橡膠種植區(qū)的合理布局，提高植膠面積，實(shí)現(xiàn)熱作區(qū)農(nóng)場(chǎng)的可持續(xù)發(fā)展。

參考文獻(xiàn)

[1] Ning Wu， Cheng qi， Shi deGuo， et al. A study of farmland landscape patterm with TM&DEM-Case study in Yongsheng county， Yunnan province， P.R.China[J]. Internation Geoscience and Remote Sensing Symposium， 2004， 5： 3 420-3 423.

[2] 陳利頂，張淑榮，傅伯杰，等. 流域尺度土地利用與土壤類型空間分布的相關(guān)性研究[J]. 生態(tài)學(xué)報(bào)，2003，23（12）：2 497-2 505.

[3] 閆小培，毛蔣興，普 軍. 巨型城市區(qū)域土地因變化的人文因素分析[J]. 地理學(xué)報(bào)，2006，61（6）：613-623.

[4] James P. McNamara， Alan D. Ziegler， S H. Wood， et al.Channel head locations with respect to geomorphologic thresholds derived from a digital elevation model： A cas study in northern Thainland[J]. Forest Ecology and Management， 2006， 224： 147-156.

[5] 代冬梅，吳國(guó)平，司惠超. 基于DEM的瀏陽(yáng)市土地利用空間格局分析[J]. 安徽農(nóng)業(yè)科學(xué)，2009，37（3）：1 231-1 234.

[6] 崔衛(wèi)國(guó)，文 倩，劉艷艷，等. 基于DEM的醴陵市土地利用空間格局分析[J]. 資源科學(xué)，2008，30（2）：228-234.

[7] 陳利頂，楊 爽，馮曉明，等. 土地利用變化的地形梯度特征與空間擴(kuò)展——以北京市海淀區(qū)和延慶縣為例[J]. 地理研究，2008，27（6）：1 225-1 234.

[8] Shen Huaifei. Land use spatial pattern characteristics along the terrain gradient in Yellow River Basin in west Henan province， China[C]. 2011 19th International Conference on Geoinformatics，Geoinformatics 2011.

[9] 陳 楠，楊武年，李 娟. 巴中市丘陵山區(qū)土地利用在地形梯度上的分布特征[J]. 水土保持通報(bào)，2012，32（1）：185-188.

[10] 程三友，李英杰，劉少峰. 基礎(chǔ)DEM的大別山地區(qū)地貌特征研究[J]. 東華理工大學(xué)學(xué)報(bào)（自然科學(xué)版），2010，33（3）：270-271.

[11] 童紹玉. 云南山區(qū)與壩區(qū)農(nóng)業(yè)利用劃分方法探析[J]. 貴州農(nóng)業(yè)科學(xué)，2011，39（11）：89-91.

[12] 唐新明，李 莉，季小燕，等. 全國(guó)七大江河流域重點(diǎn)防范區(qū)1∶1萬(wàn)數(shù)字高程模型（DEM）數(shù)據(jù)庫(kù)的建立[J]. 測(cè)繪通報(bào)，2002（6）：19-22.

[13] 孫 麗，陳煥偉，潘家文. 運(yùn)用DEM剖析土地利用類型的分布及時(shí)空變化——以北京延慶縣為例[J]. 山地學(xué)報(bào)，2004（6）：762-766.

[14] 喻 紅，曾 輝，江子瀛，等. 快速城市化地區(qū)景觀組分在地形梯度上的分布特征研究[J]. 地理科學(xué)，2001，21（2）：64-69.

[15] 湯國(guó)安，劉學(xué)軍，閭過(guò)年. 數(shù)字高程模型及地學(xué)分析的原理與方法[M]. 北京：科學(xué)出版社，2005.

[16] 馮險(xiǎn)峰，汪 閩，孟雪蓮，等. ARCGIS空間分析使用指南[Z]. ArcInfo中國(guó)技術(shù)咨詢與培訓(xùn)中心，2002.