張正明，程保銀，杜陽(yáng)

(1.江蘇省測(cè)繪產(chǎn)品質(zhì)量監(jiān)督檢驗(yàn)站，江蘇 南京　210009；　2.蘇州市測(cè)繪院有限責(zé)任公司，江蘇 蘇州　215000)

?

GeoCA模型在城市土地利用預(yù)測(cè)中的應(yīng)用研究

張正明1*，程保銀2，杜陽(yáng)2

(1.江蘇省測(cè)繪產(chǎn)品質(zhì)量監(jiān)督檢驗(yàn)站，江蘇 南京210009；2.蘇州市測(cè)繪院有限責(zé)任公司，江蘇 蘇州215000)

摘要：為了彌補(bǔ)傳統(tǒng)基于BP神經(jīng)網(wǎng)絡(luò)的GeoCA模型無(wú)法兼顧土地利用變化空間分布的缺點(diǎn)，在傳統(tǒng)BP神經(jīng)網(wǎng)絡(luò)的基礎(chǔ)上，添加空間分布集中度約束建立了GeoCA模型，對(duì)城市土地利用變化進(jìn)行預(yù)測(cè)，在保證土地利用類(lèi)型轉(zhuǎn)化概率精度的前提下，實(shí)現(xiàn)了真正意義上的土地利用變化預(yù)測(cè)。

關(guān)鍵詞：GIS；GeoCA；土地利用；空間分布

1引言

土地利用及其變化是當(dāng)前全球環(huán)境變化與可持續(xù)發(fā)展的重要內(nèi)容。將影響土地利用變化格局和過(guò)程的自然、經(jīng)濟(jì)、環(huán)境等作為驅(qū)動(dòng)因子，研究土地利用變化的驅(qū)動(dòng)機(jī)制，對(duì)土地利用變化進(jìn)行動(dòng)態(tài)模擬可以揭示土地利用變化的狀況、原因、效應(yīng)等科學(xué)問(wèn)題。一直以來(lái)，該領(lǐng)域的研究吸引了許多學(xué)者注意，建立了諸多LUCC(Land use and land cover change)模型，對(duì)土地利用進(jìn)行動(dòng)態(tài)模擬。其中比較有代表性的有荷蘭瓦赫寧根農(nóng)業(yè)大學(xué)研制的CLUE-S模型及后續(xù)研究者在其基礎(chǔ)上的一些改進(jìn)模型，如CLUE-SII模型、CLUE-CH模型等。該類(lèi)模型在土地利用時(shí)空變換模擬方面得到過(guò)廣泛應(yīng)用，但都是由相互松散的模塊構(gòu)成，缺乏統(tǒng)一的思路實(shí)現(xiàn)模型的集成。由美國(guó)麻省理工學(xué)院教授創(chuàng)立的系統(tǒng)動(dòng)力學(xué)(SD：System Dynamics)模型及其后續(xù)的改進(jìn)模型充分考慮了社會(huì)經(jīng)濟(jì)因素對(duì)土地利用變化的影響，能夠從宏觀上反映土地系統(tǒng)的復(fù)雜行為，但是缺乏對(duì)土地利用空間格局特征的模擬。40年代由ULam首先提出的元胞自動(dòng)機(jī)模型(CA，Cellular automata)，被用來(lái)研究自組織系統(tǒng)的演變過(guò)程，在土地利用動(dòng)態(tài)模擬方面研究很廣泛，并成為當(dāng)前CA應(yīng)用的熱點(diǎn)[1]。它是一種空間離散的動(dòng)力學(xué)模型，具有強(qiáng)大的空間復(fù)雜系統(tǒng)運(yùn)算能力，但是其缺乏對(duì)空間信息的處理和分析能力。為此，諸多學(xué)者如黎夏、張顯峰、劉繼生等將CA原理引入到GIS領(lǐng)域，結(jié)合GIS強(qiáng)大的空間處理與分析等優(yōu)點(diǎn)建立的地理元胞自動(dòng)機(jī)(GeoCA)模型[1]，有效地模擬了復(fù)雜土地利用的變化過(guò)程，為土地利用動(dòng)態(tài)預(yù)測(cè)的研究做出了突出貢獻(xiàn)。研究表明，基于BP神經(jīng)網(wǎng)絡(luò)(Back-propagation Neural Network)的GeoCA模型特別適用于模擬非線性系統(tǒng)[2]。本文在前人研究的基礎(chǔ)上，利用BP神經(jīng)網(wǎng)絡(luò)，提出空間集中度約束性GeoCA模型，不僅滿足了土地利用類(lèi)型轉(zhuǎn)換概率上的精度要求，還彌補(bǔ)了以往做法中難以兼顧土地類(lèi)型空間分布規(guī)律方面的不足。

2GeoCA模型在土地利用中的研究

2.1GeoCA模型原理

從本質(zhì)上講， GeoCA模型的原理是依據(jù)一定的轉(zhuǎn)換規(guī)則，判定代表土地利用信息的柵格數(shù)據(jù)在不同時(shí)刻的狀態(tài)，從而達(dá)到預(yù)測(cè)的效果。模型可以用如下公式簡(jiǎn)單描述：

A=(C，S，N，R，T)

(1)

式中，A為地理元胞自動(dòng)機(jī)；C為描述空間實(shí)體的單個(gè)地理元胞；S為元胞的狀態(tài)；N為元胞鄰域關(guān)系；R為狀態(tài)轉(zhuǎn)換規(guī)則，轉(zhuǎn)換規(guī)則的確定是模型應(yīng)用的最核心部分；T表示模型的時(shí)間維，規(guī)定在某一時(shí)刻元胞狀態(tài)唯一。模型的工作原理可以用圖1表示：

這里以選取9個(gè)元胞為一個(gè)單元為例，對(duì)中心元胞的狀態(tài)轉(zhuǎn)換進(jìn)行研究。C對(duì)應(yīng)了圖中的一個(gè)個(gè)柵格單元，代表某一具體的地理實(shí)體；S對(duì)應(yīng)了柵格單元的顏色，不同的顏色描述了同一屬性的不同取值；N對(duì)應(yīng)了除紅色以外的其他8個(gè)元胞，其狀態(tài)是影響中心元胞狀態(tài)轉(zhuǎn)化的重要因素之一。在轉(zhuǎn)換規(guī)則R的作用下，完成中心柵格從T到T+1時(shí)刻的狀態(tài)轉(zhuǎn)換，即實(shí)現(xiàn)了對(duì)地理實(shí)體某一特征的動(dòng)態(tài)變化模擬。很顯然，在整個(gè)過(guò)程當(dāng)中，轉(zhuǎn)換規(guī)則的確定至關(guān)重要，是模型有效性的決定性因素。

具體應(yīng)用到土地利用動(dòng)態(tài)預(yù)測(cè)中，則C代表了各種土地利用類(lèi)型的地塊；N代表了中心地塊周?chē)牡貕K；S代表了某一時(shí)刻各地塊對(duì)應(yīng)的不同的土地利用類(lèi)型；而轉(zhuǎn)換規(guī)則的確定，除了受N的影響之外，還受其他空間因素的影響。

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

由GeoCA模型原理可以知道，模型的建立是基于柵格數(shù)據(jù)的。除了研究區(qū)域的土地利用信息需要用柵格數(shù)據(jù)進(jìn)行表示之外，影響轉(zhuǎn)換規(guī)則的各空間因素，也要用柵格數(shù)據(jù)來(lái)描述。研究所需的各種數(shù)據(jù)均可通過(guò)遙感影像來(lái)獲取，遙感影像的分類(lèi)圖可以作為研究區(qū)域土地利用類(lèi)型的初始狀態(tài)。

研究表明，影響土地利用變換的因素包括一系列的距離變量、鄰近現(xiàn)有土地利用類(lèi)型的變量、單元的自然屬性等[3]。例如，某一單元距離市中心越近，其轉(zhuǎn)變?yōu)槌鞘杏玫氐母怕试礁?；若單元本身為水體，那它轉(zhuǎn)化為巖石的概率基本為0等。針對(duì)本研究選取的區(qū)域特點(diǎn)，選取了包括距市中心距離、距道路距離、鄰域用地類(lèi)型及柵格代表地塊的自然屬性等在內(nèi)的13個(gè)空間變量作為轉(zhuǎn)換規(guī)則的影響因子，并在Arc/Info Grid環(huán)境下對(duì)各變量進(jìn)行了提取。

在Arc/Info Grid環(huán)境中利用其內(nèi)置函數(shù)Eucdistance、Focal等提取空間距離變量和鄰域的用地類(lèi)型，并可直接用于后續(xù)建模之中?？紤]到數(shù)據(jù)表達(dá)和數(shù)據(jù)計(jì)算的難易程度，對(duì)提取的原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，用0～1之間的數(shù)值來(lái)表示原始數(shù)據(jù)的大小。

3模型算法研究

3.1算法的提出

模擬多種土地利用類(lèi)型之間的變化牽涉的轉(zhuǎn)換關(guān)系比較復(fù)雜，對(duì)于N種類(lèi)型的用地，最多可以有N×N種可能的轉(zhuǎn)換。單純的基于BP神經(jīng)網(wǎng)絡(luò)(Back-propagation Neural Network)的GeoCA模型雖然較好地預(yù)測(cè)了土地利用變化的趨勢(shì)，但其僅從轉(zhuǎn)換數(shù)量的角度對(duì)模擬結(jié)果進(jìn)行評(píng)定，對(duì)空間分布規(guī)律難以控制，具有一定的盲目性。土地利用變換是一個(gè)緩慢的過(guò)程，并且其變化特征多反映為中心擴(kuò)張規(guī)律[4]。因此土地利用變化往往呈現(xiàn)出在已有類(lèi)型周?chē)黾油N用地類(lèi)型的現(xiàn)象。本文加入空間集中度約束對(duì)BP算法進(jìn)行改進(jìn)，既保證了土地利用類(lèi)型的轉(zhuǎn)換概率，又使其盡量接近土地利用變化的空間分布狀況，較好地反應(yīng)了土地利用變化的規(guī)律。

3.2空間集中度約束性BP算法

(2)

(3)

BP算法是一種反向?qū)W習(xí)算法，定義誤差函數(shù)r為實(shí)際輸出與期望輸出之差的平方和的一半，即：

(4)

在傳統(tǒng)的算法中，該誤差反應(yīng)的是輸出概率的誤差。為了反映土地利用變化的空間分布規(guī)律，加入空間集中度誤差r′來(lái)共同約束土地利用的轉(zhuǎn)換規(guī)則。定義r′為實(shí)際位置與預(yù)測(cè)位置之差的平方和的一半，即

(5)

(6)

當(dāng)p越大時(shí)表示空間聚集度越高。

(7)

其中ε為學(xué)習(xí)步長(zhǎng)，且ε>0。

通過(guò)推導(dǎo)可得：

(8)

設(shè)：

(9)

(10)

則：

(11)

(12)

(13)

設(shè)：

(14)

4模型應(yīng)用與分析

4.1模型應(yīng)用

以某市沿海區(qū)域?yàn)檠芯繀^(qū)域，利用2005年及2010年的TM影像并結(jié)合該區(qū)域的土地利用變化歷史資料進(jìn)行分類(lèi)。研究區(qū)域的用地類(lèi)型分為建筑用地、鹽田、養(yǎng)殖區(qū)、耕地、林地等9種類(lèi)型。提取了距離市中心距離、距離公路距離、距離海岸帶距離、坡度及9種用地類(lèi)型在內(nèi)的13個(gè)空間變量作為影響因子。

選取三層BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)，即一層隱藏層的神經(jīng)網(wǎng)絡(luò)來(lái)建立模型。輸入層13個(gè)神經(jīng)元對(duì)應(yīng)選取的13個(gè)空間變量；輸出層的9個(gè)神經(jīng)元對(duì)應(yīng)了9種用地類(lèi)型。在遙感影像上隨機(jī)選取兩組包含 1 000個(gè)點(diǎn)的樣本，分別進(jìn)行神經(jīng)網(wǎng)絡(luò)的訓(xùn)練與驗(yàn)證，獲取模型參數(shù)。為了防止出現(xiàn)誤差太小而無(wú)法終止的情況，增加閾值來(lái)控制迭代次數(shù)，這里設(shè)置轉(zhuǎn)換概率和空間集中度閾值分別為0.75和0.5，迭代次數(shù)設(shè)置為 1 000，誤差期望值設(shè)置為0.3作為訓(xùn)練的終止條件。另外將學(xué)習(xí)速率設(shè)置為0.2，若學(xué)習(xí)速率設(shè)置太大，會(huì)導(dǎo)致準(zhǔn)確率難以收斂。

利用訓(xùn)練好的神經(jīng)網(wǎng)絡(luò)，在Arc/InfoGrid環(huán)境下，建立了研究區(qū)域的GeoCA模型。基于研究區(qū)域用地類(lèi)型的初始狀態(tài)，對(duì)未來(lái)幾年的土地利用變化進(jìn)行動(dòng)態(tài)模擬，與未加約束的模擬結(jié)果及實(shí)際用地情況進(jìn)行對(duì)比，驗(yàn)證研究方案的有效性。

4.2結(jié)果分析

將研究區(qū)域2005年的遙感影像分類(lèi)圖作為研究區(qū)域土地利用的初始狀態(tài)，對(duì)2010年的用地情況進(jìn)行模擬，并與傳統(tǒng)不加空間集中度約束的預(yù)測(cè)進(jìn)行結(jié)果進(jìn)行對(duì)比，分析本研究的可行性。圖2、圖3分別表示了基于不帶約束的BP神經(jīng)網(wǎng)絡(luò)GeoCA模型和論文研究的GeoCA模型模擬的統(tǒng)計(jì)結(jié)果，模擬精度分別為66.1%和64.1%。圖標(biāo)的X軸表示迭代次數(shù)，Y軸代表誤差大小。結(jié)果表明兩者在進(jìn)行模擬的過(guò)程中，類(lèi)型的轉(zhuǎn)化精度相差不大。

圖4為2010年某市未加約束情況下的土地覆蓋預(yù)測(cè)圖，圖5為2010年某市加約束情況下的土地覆蓋預(yù)測(cè)圖，圖6為2010年某市實(shí)際土地利用分布圖。

通過(guò)用地類(lèi)型分類(lèi)模擬圖與實(shí)際分類(lèi)圖的對(duì)比可以看出，基于空間集中約束的BP神經(jīng)網(wǎng)絡(luò)建立的GeoCA模型模擬的結(jié)果在空間分布上與實(shí)際狀況比較一致，而未加空間分布約束的模型模擬的結(jié)果在空間上的分布相比于實(shí)際情況較為分散。例如建筑用地和鹽田的分布就很明顯，在實(shí)際分類(lèi)圖中呈集中分布，加約束的模型模擬結(jié)果與實(shí)際情況明顯優(yōu)于未加約束的模擬結(jié)果。雖然加空間約束的模型模擬的結(jié)果不是每一種用地類(lèi)型都能優(yōu)于未加約束的模擬結(jié)果，例如鹽堿地的分布情況，但是從整體上看，其模擬結(jié)果更接近用地類(lèi)型的分布情況。

綜上可知，加空間分布約束的模型能夠在保證和傳統(tǒng)模型相同水平轉(zhuǎn)換率的情況下，使模擬結(jié)果在空間分布上更接近實(shí)際的分布情況。從而驗(yàn)證了論文研究模型在預(yù)測(cè)土地利用變化應(yīng)用中的有效性。

5結(jié)論

將CA技術(shù)引入到GIS中，建立GeoCA模型對(duì)土地利用變化進(jìn)行預(yù)測(cè)，能夠揭示土地利用變化的科學(xué)規(guī)律。傳統(tǒng)的基于BP神經(jīng)網(wǎng)絡(luò)建立的GeoCA模型雖然能夠在轉(zhuǎn)換概率上滿足預(yù)測(cè)的要求，但是其無(wú)法兼顧轉(zhuǎn)換類(lèi)型在空間上的分布，不能真正地實(shí)現(xiàn)土地利用變化的預(yù)測(cè)。本研究提出空間分布集中約束對(duì)BP神經(jīng)網(wǎng)絡(luò)進(jìn)行改進(jìn)，建立了能夠同時(shí)兼顧轉(zhuǎn)化概率和空間分布規(guī)律的GeoCA模型，并通過(guò)實(shí)際應(yīng)用驗(yàn)證了模型的有效性。

參考文獻(xiàn)

[1]周成虎，孫戰(zhàn)利，謝一春. 地理元胞自動(dòng)機(jī)研究[M]. 北京：北京出版社，1999：128～137.

[2]黎夏，劉小平，李少英. 智能式GIS與空間優(yōu)化[M]. 北京：科學(xué)出版社，2010：37～130.

[3]BattyM，XieY.Fromcellstocities.[J].EnvironmentandPlanningB，1994，21：531～548.

[4]王立波. 基于GIS和Geo-CA的營(yíng)口城市用地?cái)U(kuò)展模擬預(yù)測(cè)研究[D]. 長(zhǎng)春：東北師范大學(xué)，2009.

[5]李鳳霞，石輝，馮曉剛等. 西安市土地利用格局動(dòng)態(tài)演變及其驅(qū)動(dòng)力研究.[J]. 測(cè)繪通報(bào)，2015，12：41～56.

[6]邱小平，于丹，孫若曉等. 考慮坡道的元胞自動(dòng)機(jī)交通流模型研究遙.[J]. 計(jì)算機(jī)應(yīng)用研究，2016，33(9).

[7]陳建平，丁火平，王功文等. 基于GIS和元胞自動(dòng)機(jī)的荒漠化演化預(yù)測(cè)模型[J]. 遙感學(xué)報(bào)，2004，8(3)：254～260.[8]劉衛(wèi)芳，夏斌，聶云峰. 基于GIS與多智能體的建設(shè)用地變化仿真研究-以廣州番禺區(qū)為例[J]. 安徽農(nóng)業(yè)科學(xué)，2012，40(21)：11083～11085.

[9]王福增，尤麗敏，王凡. 大連市土地利用現(xiàn)狀及變化分析[J]. 管理科學(xué)，2008，29(6)：123～125.

[10]李德仁，張良培，夏桂松. 遙感大數(shù)據(jù)自動(dòng)分析與數(shù)據(jù)挖掘.[J]. 測(cè)繪學(xué)報(bào)，2014，43(12)：1211～1216.

Research of GeoCA Model Applied in Urban Land Use Prediction

Zhang Zhengming1，Cheng Baoyin2，Du Yang2

(1.Jiangsu Praince Surveying and Mapping Product Quality Supervion and Luspection Station，Nanjing 210009，China；2.Suzhou Surveying & Mapping Institute Co.，Ltd. Suzhou 215006，China)

Key words：GIS；GeoCA；land use；spatial distribution

Abstract：In order to make up for the shortcomings that the traditional method could not take into account the spatial distribution of the land use change，this paper based on the traditional BP neural network，GeoCA model was built by adding the Constraint conditions of spatial distribution concentration degree to predicte the change of city land use.Under the conditions of ensuring the accuracy of land use type conversion probability，the method studied in this paper truly realized the prediction of Land use change.

文章編號(hào)：1672-8262(2016)03-55-04

中圖分類(lèi)號(hào)：P208.2

文獻(xiàn)標(biāo)識(shí)碼：B

*收稿日期：2016—04—05

作者簡(jiǎn)介：張正明(1963—)，男，高級(jí)工程師，主要從事測(cè)繪產(chǎn)品質(zhì)檢方面的工作。