王文倩　高明秀　朗坤

摘要：黃河三角洲作為世界上陸地資源增長(zhǎng)最快的地區(qū)之一，一直以來(lái)都是土地利用/覆被變化研究的熱點(diǎn)區(qū)域。本研究以東營(yíng)市墾利區(qū)為例，在解譯和分析研究區(qū)2005、2009和2015年3期遙感影像的基礎(chǔ)上，通過(guò)構(gòu)建土地利用動(dòng)態(tài)分析模型組，分析研究區(qū)土地利用動(dòng)態(tài)變化特征，并預(yù)測(cè)未來(lái)的土地利用狀況。結(jié)果表明：2005-2015年間墾利區(qū)各地類均處于變動(dòng)之中，2005-2009年間研究區(qū)鹽荒地所占面積比例最大，2009-2015年間耕地面積超過(guò)鹽荒地成為墾利區(qū)面積最大的地類;在空間變化劇烈程度方面，2005-2009年土地利用變化的劇烈程度稍大于后一期。耕地、水域、城鄉(xiāng)建設(shè)用地迅速增加的面積，主要由鹽荒地轉(zhuǎn)變而來(lái)。預(yù)測(cè)到2025年，建設(shè)用地進(jìn)一步擴(kuò)張，林地和草地進(jìn)一步減少，草地將成為面積最小的地類，而耕地成為面積最大的地類。本研究為預(yù)測(cè)墾利區(qū)土地利用的發(fā)展趨勢(shì)及實(shí)現(xiàn)區(qū)域土地的可持續(xù)利用提供了科學(xué)依據(jù)。

關(guān)鍵詞：黃河三角洲;土地利用/覆被變化;模型組;動(dòng)態(tài)變化特征;趨勢(shì)預(yù)測(cè)

Change Characteristics and Trend Prediction ofLand Use in the Yellow River Delta Based on Model GroupWang Wengqian， Gao Mingxiu， Lang Kun

Abstract

The Yellow River Delta is one of the fastest growing areas of land resources in the world. Andit is the focus area of researches on land use and cover changes. In this article， taking Kenli district of Dongy-ing City as the research area， we analyzed the dynamic change characteristics and predicted the land use statusin the future by constructing the dynamic analysis model group of land use based on the interpretation of remotesensing images of 2005， 2009 and 2015. The results showed that all the land use types in Kenli Districtchanged in 2005 - 2015. The salt wasteland occupied the largest area during 2005 - 2009， and the cultivatedland area was the largest during 2009 - 2015. During the research period， the spatial dynamic change ampli-tude during 2005 -2009 was greater than that during 2009 - 2015. There was rapid growth in the area of con-struction land， cultivated land and water area， which were mainly transformed from salt wasteland. It was pre-dicted that in 2025 ， the construction land would expand further， while the woodland and grassland would de-crease further; the cultivated land would be the type with the largest area， and the grassland would be onewith the least area. These results would provide scientific bases for predicting the developing trend of land useand achieving the sustainable use of land in Kenli District.

Keywords

The Yellow River Delta; Land use and cover change ; Model group; Dynamic chnange char-acteristics ; Trend prediction

20世紀(jì)90年代，在“國(guó)際地圖與生物圈計(jì)劃（IGBP）”與“全球環(huán)境變化人文計(jì)劃（IHDP）”兩大國(guó)際組織的推動(dòng)下，土地利用/土地覆被變化（ LUCC）研究成為全球環(huán)境變化研究領(lǐng)域的熱點(diǎn)[1]。利用RS技術(shù)提取土地利用動(dòng)態(tài)變化的數(shù)據(jù)，揭示其時(shí)空變化的趨勢(shì)、驅(qū)動(dòng)機(jī)制及演變特征，是該領(lǐng)域今后研究中的重點(diǎn)方向，國(guó)內(nèi)外學(xué)者已開(kāi)展大量研究并取得一系列的成果[2]。1995年“土地利用/覆被變化科學(xué)計(jì)劃（LUCC）”提出土地利用/覆被變化的3個(gè)研究重點(diǎn)，即土地利用變化機(jī)制的案例研究、土地覆被變化機(jī)制的直接觀測(cè)和診斷模型研究以及區(qū)域和全球模型研究，旨在加深對(duì)地區(qū)土地利用動(dòng)態(tài)變化的理解，分析驅(qū)動(dòng)因子，通過(guò)改進(jìn)現(xiàn)有模型和建立新模型來(lái)預(yù)測(cè)土地利用的變化[3]。2005年全球土地計(jì)劃（Global Land Project，GLP）啟動(dòng)以后，以人類一環(huán)境耦合系統(tǒng)為核心的土地利用/覆被動(dòng)態(tài)的監(jiān)測(cè)與模擬逐漸成為研究關(guān)注的焦點(diǎn)[4]。日本等國(guó)的科研部門(mén)和學(xué)者也開(kāi)展了土地利用的驅(qū)動(dòng)因子、空間特征及對(duì)策等研究[5，6]。

黃河三角洲作為中國(guó)最年輕的土地，是世界上造陸速度最快的河口三角洲之一，其地域遼闊，自然資源十分豐富，因此它的開(kāi)發(fā)和利用備受矚目。東營(yíng)市墾利區(qū)是黃河三角洲的核心區(qū)域，其土地利用類型的演化過(guò)程對(duì)黃河三角洲而言既有典型的代表性又有其獨(dú)有的特殊性和復(fù)雜性，已有不少學(xué)者對(duì)黃河三角洲的土地利用/覆被變化特征進(jìn)行研究[7-9]。在此背景下，及時(shí)、準(zhǔn)確地掌握這一地區(qū)的土地利用信息，摸清土地利用動(dòng)態(tài)變化的特征和規(guī)律，對(duì)于采取措施合理配置和優(yōu)化利用區(qū)域土地和各種資源、推進(jìn)科學(xué)發(fā)展意義重大[10]。

本研究以墾利區(qū)作為研究區(qū)，基于RS技術(shù)獲取和處理數(shù)據(jù)，將傳統(tǒng)數(shù)量分析模型、綜合動(dòng)態(tài)分析模型與空間動(dòng)態(tài)分析模型相結(jié)合，構(gòu)建能夠表達(dá)土地利用系統(tǒng)復(fù)雜微觀格局的動(dòng)態(tài)分析模型組，全面、深入地分析區(qū)域土地利用的時(shí)空特征，并預(yù)測(cè)土地利用的發(fā)展趨勢(shì)，為區(qū)域土地的可持續(xù)利用提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

墾利區(qū)是山東省東營(yíng)市的轄區(qū)，位于山東省域的北部、黃河下游的人?？谔帲靥巹倮吞锔沟?。地理坐標(biāo)37024'～ 38010'N、118 015～llg。19'E。區(qū)域地形呈西南一東北走向，南北縱距55.5 km，東西橫距96.2km。該區(qū)處于暖溫帶季風(fēng)性氣候區(qū)，雖瀕臨渤海，但受大陸性季風(fēng)影響大，四季分明，冬季干冷，夏季濕熱，降水主要集中在夏季，春秋干旱比較嚴(yán)重[2]。地下水的礦化程度比較高，埋藏淺，極易隨毛管上升引起表層土壤的鹽漬化[11.12]，形成了以鹽化潮土為主的土壤，各地類插花分布。

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

選取研究區(qū)2005年5月24日、2009年5月19日的Landsat 4-5的TM影像和2015年5月4日的Landsat 8的OLI影像以及墾利區(qū)行政區(qū)劃矢量圖作為基礎(chǔ)數(shù)據(jù)源（表1），三期遙感圖像時(shí)間間隔約為5年，成像時(shí)間均處于5月，云量均低于0.6%，圖像清晰，影像波段所含信息豐富?；诟鞑ǘ螖?shù)據(jù)進(jìn)行統(tǒng)計(jì)特征分析、主成分分析和相關(guān)性分析[9]，最佳波段選擇波段5、波段4、波段3。

參考土地利用現(xiàn)狀分類標(biāo)準(zhǔn)（GB/T 21010-2007）和全國(guó)遙感監(jiān)測(cè)土地利用/覆被分類體系，結(jié)合研究區(qū)土地利用的實(shí)際情況，將墾利區(qū)的土地利用類型劃分為耕地、林地、草地、城鄉(xiāng)建設(shè)用地、水域、灘涂和鹽荒地七種類型。其中，灘涂即沿海灘涂，鹽荒地則為除灘涂外其他用地里的鹽堿地，包括未經(jīng)開(kāi)發(fā)的鹽堿荒地及雖經(jīng)開(kāi)發(fā)整理后又撂荒或退化的耕地、草地等。

1.3 研究方法

傳統(tǒng)的單一數(shù)量模型雖然計(jì)算方便，但忽視了土地利用的空間區(qū)位性和獨(dú)特性[13]，沒(méi)有考慮土地利用過(guò)程中的空間變化和相關(guān)屬性。為了能夠更全面地分析墾利區(qū)土地利用變化現(xiàn)狀，預(yù)測(cè)將來(lái)的土地利用變化情況，本研究參考多個(gè)空間分析模型，應(yīng)用傳統(tǒng)土地利用變化幅度模型、空間動(dòng)態(tài)模型和轉(zhuǎn)移矩陣模型[7]構(gòu)建成一個(gè)鏈狀模型組（表2），用以分析土地利用變化的特征。

1.4 數(shù)據(jù)處理

對(duì)遙感圖像首先進(jìn)行幾何精校正[14，15]，然后利用基于面向?qū)ο蟮闹С窒蛄繖C(jī)方法[16]并結(jié)合目視解譯，提取土地利用及其覆被信息特征，經(jīng)分類處理后作出精度評(píng)價(jià)。經(jīng)計(jì)算，本研究所選墾利區(qū)三幀遙感影像均符合總體精度達(dá)到80%以上的要求，其土地利用/覆被結(jié)果如圖1所示。

2 結(jié)果與分析

2.1 墾利區(qū)土地利用結(jié)構(gòu)變化

利用ENVI軟件，基于3期遙感圖像統(tǒng)計(jì)墾利區(qū)的土地利用類型，并分別計(jì)算10年間各土地利用類型面積的變化幅度，結(jié)果見(jiàn)表3。

2005-2015年間，墾利區(qū)的主要土地利用類型所占比重都發(fā)生了不同程度的變化，耕地、水域、城鄉(xiāng)建設(shè)用地增加，而鹽荒地、草地、灘涂、林地面積縮小。未利用地（灘涂、鹽荒地）是墾利區(qū)面積最大的土地類型，各期的面積分別占全區(qū)總面積的52.20%、46.57%和44.63%。鹽荒地是時(shí)段內(nèi)減少幅度最大的土地類型，由2005年的68 789.72 h㎡縮減至2015年的53194.05 h㎡，總體降幅達(dá)到6.69%;灘涂也相應(yīng)減少了2065. 86h㎡，總體降幅為0.88%。耕地是僅次于未利用地的主要用地類型，占全區(qū)總面積的比例由2005年的18.73%增長(zhǎng)至2015年的23.17%，總體增幅4.44%.2009-2015年間耕地的增幅相對(duì)于2005-2009年間更為突出，是變化幅度第二的土地利用類型。水域是增加幅度僅次于耕地的土地利用類型，共增長(zhǎng)了3.51%。表明，近幾年墾利區(qū)十分重視鹽荒地的改造和耕地保護(hù)，加大了對(duì)鹽荒地、灘涂等多種難利用地的開(kāi)墾，為墾利區(qū)的發(fā)展提供了豐富的后備土地資源，但其難利用性也是限制墾利區(qū)經(jīng)濟(jì)發(fā)展的重要因素之一。

2.2 墾利區(qū)土地利用空間變化

傳統(tǒng)的土地利用變化分析模型（變化幅度模型）僅限于對(duì)土地利用面積進(jìn)行單純的數(shù)量分析，動(dòng)態(tài)度模型只考慮土地變化的某一方面，無(wú)法測(cè)算土地利用變化的綜合活躍度?？臻g分析模型是由劉盛和等[13]提出的一種基于轉(zhuǎn)移矩陣的測(cè)算模型，是在土地動(dòng)態(tài)度模型的基礎(chǔ)上進(jìn)行修正，加入了對(duì)土地新增部分的測(cè)算。

利用ENVI軟件，從3期土地利用現(xiàn)狀圖中提取出轉(zhuǎn)移矩陣，基于模型2計(jì)算2005-2015年間各種土地利用類型的動(dòng)態(tài)變化情況，如表4所示：

（1）草地是10年間總變化速率排名第一的土地利用類型，但因其初期基數(shù)小，總體變化面積并不是最高的。其轉(zhuǎn)出速率較高，明顯高于其他土地利用類型，但轉(zhuǎn)入速率較低，僅1.34%，屬于高速衰減型。

（2） 2005-2015年間鹽荒地轉(zhuǎn)移速度較快，轉(zhuǎn)出面積及其所占比例最大，且轉(zhuǎn)出速率遠(yuǎn)大于轉(zhuǎn)入速率，也屬于高速衰減型，表明鹽荒地是其他用地類型的主要新增來(lái)源;但由于其初期總體基數(shù)較大，總體的變化速率較低。林地、灘涂的變化趨勢(shì)同鹽荒地一致，也是轉(zhuǎn)出部分大于轉(zhuǎn)入部分，均屬于衰減型，但轉(zhuǎn)出面積占轉(zhuǎn)移總面積的比例均不足10%。

（3）城鄉(xiāng)建設(shè)用地、耕地、水域三種用地類型新增部分均遠(yuǎn)遠(yuǎn)大于轉(zhuǎn)出部分，且轉(zhuǎn)出速率遠(yuǎn)小于新增速率，均屬于高速擴(kuò)展型。耕地是新增面積及其占比最高的土地利用類型，其新增面積是轉(zhuǎn)出面積的3倍以上，是研究區(qū)其他用地轉(zhuǎn)出的主要接收者;水域作為新增面積排名第二的土地利用類型，其新增面積是同期轉(zhuǎn)出面積的6. 25倍，新增速率也遠(yuǎn)大于轉(zhuǎn)出速率;城鄉(xiāng)建設(shè)用地的新增面積占區(qū)域新增總面積的12.24%，遠(yuǎn)遠(yuǎn)大于同期的轉(zhuǎn)出面積，轉(zhuǎn)入速率是其轉(zhuǎn)出速率的34.58倍。

墾利區(qū)2005-2015年土地利用變化的空間分析模型測(cè)算結(jié)果如表5所示，可以看出，在2005-2009年間，墾利區(qū)土地利用變動(dòng)的劇烈程度由大到小依次為水域>鹽荒地>林地>草地>城鄉(xiāng)建設(shè)用地>耕地>灘涂;而在2009-2015年，變化劇烈程度的排序略有不同，依次為草地>林地>鹽荒地>耕地>城鄉(xiāng)建設(shè)用地>水域>灘涂。

兩階段中，林地和草地的空間動(dòng)態(tài)度不降反升，尤其是草地，由前一期的5.39%增長(zhǎng)至7.65%;其余類型用地的空間動(dòng)態(tài)度則均降低，表明這些土地利用類型在2005 -2009年階段的利用變化更為強(qiáng)烈，其中最明顯的是水域，由6.58%降至1.79%，反映出水域經(jīng)過(guò)2005-2009年間的劇烈變動(dòng)后，在2009-2015年間逐漸趨于穩(wěn)定。此外，鹽荒地的動(dòng)態(tài)度也由6.10%降至4.22%，其開(kāi)發(fā)利用的困難性是造成其動(dòng)態(tài)度降低的主要原因。

2.3 墾利區(qū)土地利用流轉(zhuǎn)分析

為了更清楚地了解各類土地利用類型之間的流轉(zhuǎn)關(guān)系，本研究引入了轉(zhuǎn)移矩陣模型[17]，轉(zhuǎn)移矩陣能夠清晰地反映出某種土地利用類型的增加來(lái)源和消失去向，可為土地規(guī)劃和決策提供數(shù)據(jù)來(lái)源和依據(jù)。

由表6可見(jiàn)，研究區(qū)城鄉(xiāng)建設(shè)用地面積的增加主要來(lái)源于鹽荒地，2005-2015年間城鄉(xiāng)建設(shè)用地的轉(zhuǎn)入面積為4606.97h㎡，其中由鹽荒地轉(zhuǎn)入的面積為3 088.45h㎡，占總體的67.04%;其次來(lái)源于耕地，占總體的25.98%。耕地面積增加的主要來(lái)源同樣為鹽荒地，10年間由鹽荒地轉(zhuǎn)入的面積達(dá)9450.63 h㎡，占總體的68.02%，另有少量水域（內(nèi)陸?zhàn)B殖用地）、少數(shù)農(nóng)村居民點(diǎn)、部分自然生長(zhǎng)的林草地被開(kāi)墾為耕地;在轉(zhuǎn)出方面.耕地主要與城鄉(xiāng)建設(shè)用地、鹽荒地之間存在流轉(zhuǎn)關(guān)系，城市發(fā)展占用耕地，不當(dāng)?shù)耐恋乩梅绞狡茐母貙?dǎo)致土地退化。水域（水體及水利設(shè)施用地）面積的增加主要來(lái)自鹽荒地和沿海灘涂，10年間由這兩種地類轉(zhuǎn)入的面積占總體的84.35%，還有少部分耕地、林草地被開(kāi)發(fā)用于內(nèi)陸?zhàn)B殖業(yè)。鹽荒地作為轉(zhuǎn)移面積最大的用地類型，主要轉(zhuǎn)移至耕地、水域、城鄉(xiāng)建設(shè)用地，占轉(zhuǎn)移總體的92.33%，還有少部分被整理為林地和草地;在轉(zhuǎn)入方面，草地、林地、耕地與鹽荒地的流轉(zhuǎn)比較密切，可見(jiàn)，鹽荒地的整改在近些年初見(jiàn)成效，但難度增加，同期的土地退化現(xiàn)象同樣嚴(yán)重。

2.4 未來(lái)10年墾利區(qū)的土地利用趨勢(shì)預(yù)測(cè)

土地利用轉(zhuǎn)移矩陣除了能夠分析土地利用類型的空間變化和流轉(zhuǎn)情況，還能夠根據(jù)矩陣算出它的轉(zhuǎn)移概率矩陣，進(jìn)而預(yù)測(cè)今后土地利用可能的變化結(jié)果。利用軟件導(dǎo)出2005-2015年墾利區(qū)的土地利用轉(zhuǎn)移矩陣（表6），并通過(guò)模型4計(jì)算出10年間土地利用的轉(zhuǎn)移概率矩陣，據(jù)此來(lái)推測(cè)2015-2025年墾利區(qū)的土地利用狀況。由表7可見(jiàn)，2015-2025階段，墾利區(qū)的城鄉(xiāng)建設(shè)用地進(jìn)一步擴(kuò)展，由2015年的6.68%上漲至2025年的8.40%，增加了4013.03h㎡;耕地?cái)U(kuò)張速度明顯減慢，2015 -2025年的增加量比2005-2015年少3889. 26h㎡，但耕地將成為墾利區(qū)所占比例最大的土地利用類型，占比達(dá)25.94%;林地和草地進(jìn)一步減少，到2025年，草地僅占2.64%，林地占4.33%，成為墾利區(qū)面積最少的兩種土地利用類型;鹽荒地進(jìn)一步減少，占比由2015年的22.82%下降到2025年的17.86%，表明未來(lái)鹽荒地的開(kāi)發(fā)利用前景仍十分廣闊。上述結(jié)果僅僅是以2005-2015年間的情況為基準(zhǔn)推測(cè)出來(lái)的，而新時(shí)期政策、科學(xué)技術(shù)等的變化都會(huì)引起未來(lái)土地利用類型改變，因此本研究結(jié)果僅作為對(duì)未來(lái)用地規(guī)劃的參考使用。

3 結(jié)論

本研究利用2005、2009、2015年三期遙感影像，通過(guò)簡(jiǎn)單的變化幅度模型及復(fù)雜的空間分析和轉(zhuǎn)移矩陣模型，分別從時(shí)間和空間的角度上分析了黃河三角洲核心地墾利區(qū)近10年間土地利用/覆被的變化特征，主要結(jié)論如下：

（1）相較于傳統(tǒng)的單一模型，基于模型組對(duì)黃河三角洲墾利區(qū)的分析能更加詳細(xì)地反映當(dāng)?shù)氐耐恋乩每臻g變化特征及變化趨勢(shì)，可以有效兼顧土地利用的空間區(qū)位性和獨(dú)特性，并能直觀反映土地利用過(guò)程中的空間變化和相關(guān)屬性。

（2）該地區(qū)在2005-2009年間以鹽荒地為主要用地類型，且很大一部分鹽荒地被開(kāi)墾整理成耕地，還有部分被開(kāi)發(fā)為建設(shè)用地，到2015年，耕地超越鹽荒地成為最主要的用地類型，說(shuō)明鹽荒地得到了有效地開(kāi)發(fā)與治理。目前該區(qū)仍有大面積的鹽荒地未得到有效利用，同時(shí)又有大量土地鹽荒化，說(shuō)明墾利區(qū)土地荒蕪與鹽堿化問(wèn)題仍十分嚴(yán)峻，但也反映了該區(qū)土地利用的較大潛力。

（3）墾利區(qū)城鄉(xiāng)建設(shè)用地一直表現(xiàn)為高速擴(kuò)張的態(tài)勢(shì)，土地利用空間動(dòng)態(tài)度一直處于較高位置，體現(xiàn)了研究區(qū)快速城鎮(zhèn)化的趨勢(shì)。城市擴(kuò)張對(duì)土地的需求及粗放式的土地管理是目前保護(hù)耕地所要面臨的首要問(wèn)題。

本研究為及時(shí)掌握黃河三角洲地區(qū)土地利用情況及地類變化的動(dòng)態(tài)信息提供了有效的技術(shù)手段與方法，同時(shí)系統(tǒng)探討了墾利區(qū)的土地利用空間變化特征及變化趨勢(shì)，為進(jìn)一步土地規(guī)劃與可持續(xù)利用提供了科學(xué)依據(jù)。

參考文獻(xiàn)：

[1] 張成揚(yáng)，趙智杰.近10年黃河三角洲土地利用/覆蓋時(shí)空變化特征與驅(qū)動(dòng)因素定量分析[J].北京大學(xué)學(xué)報(bào)（自然科學(xué)版），2015，51（1）：151 -158.

[2] 李百紅，趙庚星，董超，等.基于遙感和GIS的黃河三角洲鹽化土地動(dòng)態(tài)及其驅(qū)動(dòng)力分析[J].自然資源學(xué)報(bào)，2011，26（2）：310 -318.

[3] 羅湘華，倪晉仁.土地利用/土地覆蓋變化研究進(jìn)展[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào)，2000，8（3）：262 -272.

[4] 劉紀(jì)遠(yuǎn)，張?jiān)鱿椋煨铝?，?21世紀(jì)初中國(guó)土地利用變化的空間格局與驅(qū)動(dòng)力分析[J].地理學(xué)報(bào)，2009，64 （12）：1411 - 1420.

[5] 王昭生，潘愛(ài)民，肖志坤，等.基于異速增長(zhǎng)模型的土地利用變化驅(qū)動(dòng)因子分析——以湘潭市為例[J].水土保持研究，2010，17（5）：40 - 44，50.

[6] 尹慧君.天津市塘沽地區(qū)土地利用空間結(jié)構(gòu)優(yōu)化研究[D].天津：天津大學(xué)，2008.

[7] 王宏亮，郝晉珉，高陽(yáng)，等.基于多模型測(cè)度的內(nèi)蒙古土地利用動(dòng)態(tài)變化分析[J].中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào)，2017，22（4）：59 -66.

[8] 趙志杰，張成揚(yáng).黃河三角洲土地利用/覆蓋時(shí)空變化研究[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào)，2016，24（4）：731 -740.

[9]李暢，李桂娥，謝志文，等.近15年來(lái)黃河三角洲土地利用演化及成因分析[J].科技導(dǎo)報(bào)，2015，33（22）：37 -44.

[10]謝興震，史同廣.東營(yíng)市土地利用/覆被變化分析及驅(qū)動(dòng)力研究[J].綠色科技，2013，50（2）：6-8.

[11]常春艷，趙庚星，王凌，等.黃河口生態(tài)脆弱區(qū)土地利用時(shí)空變化及驅(qū)動(dòng)因素分析[J].農(nóng)業(yè)工程學(xué)報(bào)，2012，28 （24）：226 -234.

[12]李靜，趙庚星，楊配國(guó).基于知識(shí)的墾利區(qū)土地利用/覆被遙感信息提取技術(shù)研究[J].科學(xué)通報(bào)，2006，51（增刊）：183 -188.

[13]劉盛和，何書(shū)金.土地利用動(dòng)態(tài)變化的空間分析測(cè)算模型[J].自然資源學(xué)報(bào)，2002，17（5）：533 -540.

[14]趙庚星，林戈.黃河口耕地遙感動(dòng)態(tài)監(jiān)測(cè)及其生態(tài)環(huán)境安全分析[J].水土保持學(xué)報(bào)，2004，18（2）：37 -40.

[15]畢如田，白中科，李華，等.基于RS和GIS技術(shù)的露天礦區(qū)土地利用變化分析[J].農(nóng)業(yè)工程學(xué)報(bào)，2008，24（ 12）：201 - 204.

[16]常春艷，趙庚星，王凌，等.濱海光譜混淆區(qū)面向?qū)ο蟮耐恋乩眠b感分類[J].農(nóng)業(yè)工程學(xué)報(bào)，2012，28（5）：226 -231.

[17]陶蘊(yùn)之，張?zhí)穑虹?，?基于轉(zhuǎn)移矩陣的環(huán)渤海土地利用研究[J].綠色科技，2016，53（1）：4-7.