汪東川，張利輝

(天津城市建設(shè)學(xué)院 管理工程系，天津 300384)

軌跡分析與元胞自動機(jī)在土地利用動態(tài)模擬中的應(yīng)用

汪東川，張利輝

(天津城市建設(shè)學(xué)院 管理工程系，天津 300384)

元胞自動機(jī)模型是一種時間、空間、狀態(tài)都離散，空間上相互作用及時間上因果關(guān)系皆局部的網(wǎng)絡(luò)動力學(xué)模型．在該類模型研究中，經(jīng)常將其與其他類型模型相結(jié)合使用，以提高模擬有效性．軌跡分析是一種可用來分析時序動態(tài)變化的新方法．本文在總結(jié)、分析元胞自動機(jī)與軌跡分析模型各自特點的基礎(chǔ)上，探討了兩者在土地利用動態(tài)變化研究中聯(lián)合應(yīng)用的可能性，為以后建立基于軌跡分析與元胞自動機(jī)的土地利用動態(tài)模型奠定理論基礎(chǔ)．通過對天津濱海新區(qū)土地利用變化進(jìn)行初步模擬，證實了該理論的可行性.

軌跡分析；元胞自動機(jī)；土地利用動態(tài)變化；模擬

目前，土地利用模型研究越來越受到各國學(xué)者的重視．在前期土地利用/土地覆蓋變化(land use/cover change，簡稱 LUCC)機(jī)理研究的基礎(chǔ)上，研究者利用土地利用模擬模型研究土地利用/土地覆蓋的驅(qū)動機(jī)制，找出影響LUCC的驅(qū)動因素，建立LUCC的預(yù)測模型，以預(yù)測其未來變化，從而為土地利用合理評價及未來規(guī)劃提供參考與指導(dǎo)．LUCC預(yù)測模型的研究已經(jīng)成為目前LUCC研究的重點問題之一．由于土地利用/土地覆蓋本身具有較強(qiáng)的空間性，同時要研究其變化需要不同時期的土地利用數(shù)據(jù)，具有明顯的時序特征．所以LUCC預(yù)測模型研究也正朝著時空變化模型的方向發(fā)展.

元胞自動機(jī)(cellular automata，簡稱 CA)模型是20世紀(jì)40年代由Ulam提出的，具有模擬復(fù)雜系統(tǒng)時空演化過程的能力．CA的一個重要特點是復(fù)雜系統(tǒng)可以由一些簡單的局部規(guī)則來產(chǎn)生．它這種“自下而上”的研究思路充分體現(xiàn)了復(fù)雜系統(tǒng)局部個體行為產(chǎn)生全局、有秩序模式的理念，非常適用于復(fù)雜地理過程的模擬和預(yù)測[1]．根據(jù)定義的轉(zhuǎn)換規(guī)則進(jìn)行學(xué)習(xí)和反饋，以實現(xiàn)土地利用的優(yōu)化布局和空間模擬．該類模型一般要以計算機(jī)編程為基礎(chǔ)，依托GIS平臺實現(xiàn)．

在利用CA模型進(jìn)行LUCC模擬預(yù)測時，經(jīng)常將該模型與其他模型相結(jié)合來構(gòu)造元胞轉(zhuǎn)換規(guī)則，以提高模擬的有效性．CA-Markov模型是將元胞自動機(jī)模型與Markov模型相結(jié)合的、典型的LUCC模擬預(yù)測模型，許多學(xué)者運用該模型對不同區(qū)域進(jìn)行土地利用預(yù)測模擬[1-9]，提出了一種基于蟻群智能來自動獲取地理元胞自動機(jī)轉(zhuǎn)換規(guī)則的新方法；何春陽等[10]將元胞自動機(jī)與系統(tǒng)動力學(xué)相結(jié)合建立模型進(jìn)行模擬；此外，還有系統(tǒng)動力學(xué)與 CLUE-S結(jié)合的模擬[11]以及神經(jīng)網(wǎng)絡(luò)與元胞自動機(jī)結(jié)合的模擬[12]．

本文在總結(jié)、分析元胞自動機(jī)與軌跡分析模型各自特點的基礎(chǔ)上，探討了兩者在土地利用動態(tài)變化研究中聯(lián)合應(yīng)用的可能性，為日后建立基于軌跡分析與元胞自動機(jī)的LUCC動態(tài)預(yù)測模型奠定理論基礎(chǔ).

1 CA的基本特征

元胞自動機(jī)模型是一種時間、空間、狀態(tài)都離散，空間上相互作用及時間上因果關(guān)系皆局部的網(wǎng)絡(luò)動力學(xué)模型．它是現(xiàn)代系統(tǒng)科學(xué)、非線性科學(xué)與人工生命、遺傳算法相互交叉、滲透產(chǎn)生的，其固有的、強(qiáng)大的并行計算能力以及時空動態(tài)特征，使得它在模擬空間復(fù)雜系統(tǒng)的時間動態(tài)演變方面具有自然性、合理性和可行性．散布在規(guī)則格網(wǎng)中的每一元胞(cell)取有限的離散狀態(tài)，遵循同樣的作用規(guī)則，依據(jù)確定的局部規(guī)則作同步更新．大量元胞通過簡單的相互作用而構(gòu)成動態(tài)系統(tǒng)的演化．元胞自動機(jī)模型(CA模型)是一種主要的智能模擬預(yù)測模型，該模型通過對其轉(zhuǎn)換規(guī)則的定義和擴(kuò)展，可以將最適宜的用地分配給優(yōu)先配置的用地類型，有效地實現(xiàn)土地利用的時空動態(tài)模擬．不同于一般的動力學(xué)模型，元胞自動機(jī)不是由嚴(yán)格定義的物理方程或函數(shù)確定，而是用一系列模型構(gòu)造的規(guī)則構(gòu)成．凡是滿足這些規(guī)則的模型都可以算作是元胞自動機(jī)模型．因此，元胞自動機(jī)是一類模型的總稱，或者說是一個方法框架.

CA模型最大的優(yōu)點是通過定義局部的元胞鄰近關(guān)系，以及使用比較簡單的作用于元胞鄰域上的局部的轉(zhuǎn)換規(guī)則，可以模擬和表示整個系統(tǒng)中復(fù)雜現(xiàn)象的時空動態(tài)變化．元胞自動機(jī)具有離散性、同步性、齊質(zhì)性、局部性、時空動態(tài)性和簡單性等特征[13].

2 軌跡分析(trajectory analysis)方法

軌跡分析是一種可用來分析時序動態(tài)變化的新方法．該方法的首要前提是認(rèn)為土地覆蓋變化中的許多現(xiàn)象在變化前后均存在特殊的時間連續(xù)性，致使光譜空間產(chǎn)生獨特的時序標(biāo)志[14]．軌跡分析方法可以同時提供非連續(xù)現(xiàn)象(如干擾時間與強(qiáng)度)與連續(xù)現(xiàn)象(如干擾后的恢復(fù))的估算．時間序列過程中的變化軌跡可以用變化代碼來表示，如AAAAA，AABEE，11112，22233等，表現(xiàn)形式多樣．其中每個代碼的字母或數(shù)字等代表每一個時期土地利用分類圖中相應(yīng)柵格上面的土地利用類型，這樣就可以表示該柵格對應(yīng)的地點在研究時序內(nèi)各個時間節(jié)點上土地利用類型的變化軌跡．軌跡分析的目標(biāo)與特點主要有以下幾點．

(1)軌跡分析是一種多變量數(shù)據(jù)分析技術(shù)，目標(biāo)是揭示自然系統(tǒng)運行的規(guī)律性，識別其驅(qū)動機(jī)制，最終找出其控制原則.

軌跡分析模型是一個簡單的概念模型．該原型由三部分組成：射彈、彈道和驅(qū)動力．當(dāng)它被用在自然的植被/環(huán)境系統(tǒng)上時，射彈就變成了植被的狀態(tài)，彈道軌跡就變成了植被演化過程，驅(qū)動力就變成了驅(qū)動環(huán)境變量[15]．模型結(jié)構(gòu)比較簡單，在自然環(huán)境下可以得到很好的限定，并且很容易被應(yīng)用于時間序列的植被數(shù)據(jù)，與數(shù)據(jù)的來源無關(guān)．其目的就是基于每一點的變化過程進(jìn)行研究，揭示自然系統(tǒng)運行的規(guī)律性，識別其驅(qū)動機(jī)制，最終找出其控制原則.

(2)軌跡分析主要用來識別景觀隨時間變化的軌跡過程.

通過將變化歷史與具體變化位置相結(jié)合，不僅可以描述格局的時序變化規(guī)律，而且可以同時考慮他們的連續(xù)性和方向性．所以，強(qiáng)調(diào)變化軌跡很有用．例如，可以重建過去的變化過程，分析其變化的主要驅(qū)動力，或者基于相關(guān)的政治、經(jīng)濟(jì)和技術(shù)條件發(fā)生變化的情況下景觀變化趨勢預(yù)測等[16-17]．過去的變化規(guī)律影響當(dāng)前不同景觀元素的生態(tài)特征，同時影響其在規(guī)劃策略制定方面所起的作用，因此這種動態(tài)過程重建變得更加重要[18].

(3)軌跡分析將時間維與景觀格局相結(jié)合，同時，軌跡分析方法的多樣性展現(xiàn)了這種結(jié)合所面臨的方法的挑戰(zhàn).

在軌跡分析中，每個軌跡是一條單維度的線，在定義了任何時間和空間上的點的軌跡狀態(tài)的多維度空間中，表示狀態(tài)隨時間的變化趨勢．軌跡的概念意味著要使用過去的歷史序列數(shù)據(jù)，利用調(diào)查數(shù)據(jù)找出植被變化過程定量的模擬分析，進(jìn)而揭示相關(guān)過程規(guī)律及管理的信息.

(4)軌跡分析是一種綜合的方法框架，在該框架下，可以將不同的具體方法集成應(yīng)用.

在軌跡分析框架下，各種不同的分析方法可以被集成應(yīng)用．如空間分析與多變量統(tǒng)計方法的合成等，表明了軌跡分析方法在景觀動態(tài)變化分析方面具有較好的發(fā)展?jié)摿Γ梢詰?yīng)用各種光譜技術(shù)來探測植被變化的各方面信息．只要特有的時間或空間維度的最小要求得到滿足，數(shù)據(jù)就可以用來分析處理，當(dāng)然數(shù)據(jù)的詳盡程度經(jīng)常對結(jié)果具有限制作用.

(5)軌跡分析的前提假設(shè)是變化由隨機(jī)變化和規(guī)律性變化兩部分組成.

一般的生態(tài)過程表明隨機(jī)變化可能控制了短期變化過程，特別是在具體的小尺度高度同質(zhì)的環(huán)境下，例如林隙動態(tài)(gap dynamics)；規(guī)律性變化是指環(huán)境變化在一定的時間和空間尺度內(nèi)的總體趨勢，表示其變化的真實類型和強(qiáng)度．時空尺度并不會隔離變化的隨機(jī)與規(guī)律性部分．隨機(jī)變化與規(guī)律性變化在自然界內(nèi)以波的層疊方式回旋前進(jìn)，在整個系統(tǒng)內(nèi)是密不可分的[15]．與其他模型相比，軌跡分析模型的主要優(yōu)點是該模型能夠從最高的過程完整性上來探求其變化規(guī)律性，而不是像其他模型那樣割裂其基本變化過程，在初級水平上探討植被動態(tài)來獲取信息規(guī)律．

3 軌跡分析與CA結(jié)合進(jìn)行土地利用動態(tài)模擬的可能性

景觀格局演化過程模擬是一項十分復(fù)雜的系統(tǒng)模擬，而從微觀的角度入手建構(gòu)新的方法，從而揭示系統(tǒng)的規(guī)律性，一直是復(fù)雜系統(tǒng)研究較為活躍的領(lǐng)域．元胞自動機(jī)與軌跡分析的各自特點決定了它們都適合于空間信息的時空動態(tài)分析，尤其是時空動態(tài)過程的模擬.

(1)在結(jié)構(gòu)上，元胞自動機(jī)采用“自下而上”的構(gòu)模方式，符合景觀作為復(fù)雜系統(tǒng)的形成規(guī)律及其研究方法.“復(fù)雜來自于簡單”，只有從景觀單元的狀態(tài)和行為入手，模擬它們的相互作用，才能從根本上解決復(fù)雜性問題．而且，它沒有一個既定的數(shù)學(xué)方程，只是一個建模原則，因此具有很好的開放性和靈活性，更符合人們認(rèn)識復(fù)雜事物的思維方式[19]；軌跡分析模型結(jié)構(gòu)也比較簡單，在自然環(huán)境下可以得到很好的限定，它是一種綜合的方法框架，在該框架下，可以將不同的具體方法集成應(yīng)用，與將完整的變化過程分開進(jìn)行處理的抽象數(shù)學(xué)模型形成了鮮明的對比.

(2)在時空動態(tài)模擬方面，元胞自動機(jī)是一個基于微觀個體相互作用的時空動態(tài)模擬模型．將地理實體的空間和時間特性統(tǒng)一在模型中，通過劃分研究對象的元胞空間和研究初始狀態(tài)及狀態(tài)轉(zhuǎn)換規(guī)則，元胞自動機(jī)就可以自行迭代運算，模擬景觀的演化過程；同樣，軌跡分析主要用來識別景觀隨時間變化的軌跡過程，具有時空動態(tài)特征.

(3)從對數(shù)據(jù)對象要求的角度來看，元胞自動機(jī)具有不依賴數(shù)據(jù)比例尺的概念．元胞只是提供了一個行為空間，本身不受元胞空間測度和時間測度的影響，時空測度的影響通過轉(zhuǎn)換規(guī)則體現(xiàn)．因此，元胞自動機(jī)可以用來模擬局部的、區(qū)域的景觀演化過程；軌跡分析很容易被應(yīng)用于時間序列數(shù)據(jù)，與數(shù)據(jù)的來源無關(guān).

(4)從數(shù)據(jù)模型的角度看，元胞自動機(jī)中的元胞和基于柵格GIS中的柵格一樣，所以元胞自動機(jī)易于和 GIS、遙感數(shù)據(jù)處理等系統(tǒng)集成；軌跡分析通過基于每一點(可以是GIS中的柵格，或者CA中的元胞)的變化過程進(jìn)行研究，揭示自然系統(tǒng)運行的規(guī)律性，識別其驅(qū)動機(jī)制，最終找出其控制原則.

(5)在對連續(xù)時空的模擬方法上，元胞自動機(jī)將空間和時間離散化，適合于建立計算機(jī)模型．但是在元胞自動機(jī)中，元胞在t+1的時刻只直接決定于t時刻該元胞及鄰居元胞的狀態(tài)，與 t以前的時刻無關(guān)；軌跡分析將時間維與景觀格局相結(jié)合，其對客觀世界的表示也是離散的，但該模型能夠從最高的過程完整性上來探求其變化的規(guī)律性，這對于元胞自動機(jī)來說是一個有益的補(bǔ)充.

4 軌跡分析與CA結(jié)合進(jìn)行土地利用動態(tài)模擬的關(guān)鍵問題

在軌跡分析與 CA模型結(jié)合進(jìn)行土地利用動態(tài)模擬應(yīng)用中，需要解決的關(guān)鍵問題有兩點：①模擬步長的確定．元胞自動機(jī)的模擬時間步長是對時間序列連續(xù)模擬成功的一個前提條件．利用軌跡分析找出適合于模擬區(qū)域的演變步長，可以為模擬結(jié)果的可信度奠定基礎(chǔ)；②元胞自動機(jī)的元胞轉(zhuǎn)換規(guī)則的建立．元胞轉(zhuǎn)換規(guī)則是土地利用動態(tài)模擬的關(guān)鍵．軌跡分析的優(yōu)點之一就是可以通過時間序列連續(xù)變化過程的時空分析，探討土地利用的時空演變規(guī)律，從而確定元胞轉(zhuǎn)換規(guī)則．所以兩者的結(jié)合能夠很好地解決這一關(guān)鍵問題.

5 基于軌跡分析與CA模型的濱海新區(qū)土地利用變化初步模擬試驗

在研究CA模型的基礎(chǔ)上，結(jié)合LUCC軌跡特點以及軌跡分析思想，對濱海新區(qū)土地利用發(fā)展變化情況進(jìn)行了初步模擬試驗.

通過對三期遙感影像(包括1993年Landsat TM影像、2001年Landsat ETM+影像和2009年HJ1A影像)進(jìn)行計算機(jī)分類，在子類提取的基礎(chǔ)上進(jìn)行類別綜合，得到幾個綜合的土地利用類型，分別用代碼1—4代替(見表1)，得到各個時間節(jié)點的土地利用分類圖.

表1 土地利用分類代碼表

利用軌跡分析對各期土地利用圖進(jìn)行軌跡計算，得到土地利用變化軌跡空間分布圖(見圖 1)．圖 1代碼中的每個數(shù)字代表三個時間節(jié)點土地利用類型，代碼整體代表土地利用變化過程.

圖1 土地利用變化軌跡模擬

通過與影響土地利用變化的自然、社會因子相疊加進(jìn)行軌跡分析，獲取土地利用時空演變規(guī)律，制定元胞自動機(jī)的元胞轉(zhuǎn)換規(guī)則．元胞轉(zhuǎn)換規(guī)則中主要考慮四種轉(zhuǎn)換概率：①先驗概率(PA)；②研究區(qū)內(nèi)軌跡轉(zhuǎn)移概率(PB)；③元胞空間內(nèi)轉(zhuǎn)移概率(PC)；④元胞空間內(nèi)地類概率統(tǒng)計(PD)．通過上述各個影響因素來確定模擬期土地利用各個元胞上每一種地類在前一步長中的轉(zhuǎn)移概率．

通過 VC++ 編程語言實現(xiàn)預(yù)測轉(zhuǎn)換規(guī)則(部分代碼見圖2)，對研究區(qū)土地利用時空變化進(jìn)行模擬．

將野外調(diào)查數(shù)據(jù)與歷史數(shù)據(jù)相結(jié)合獲取50個點的土地利用變化真實數(shù)據(jù)，對模擬結(jié)果進(jìn)行驗證，模擬精度達(dá)到80%以上，效果較好.

圖2 基于軌跡分析和元胞自動機(jī)的土地利用模擬程序代碼(部分)

6 結(jié) 論

根據(jù)以上分析，軌跡分析方法能夠很好地與元胞自動機(jī)模型相結(jié)合，突出時空動態(tài)特征．通過對濱海新區(qū)模擬過程與實際演變過程進(jìn)行比較，進(jìn)一步證明元胞自動機(jī)與軌跡分析具有較好的相互結(jié)合潛力．軌跡分析通過與土地利用影響因子相結(jié)合進(jìn)行時空動態(tài)分析，從而為元胞自動機(jī)提供元胞演變規(guī)則，使其能夠進(jìn)行土地利用的動態(tài)模擬，并為元胞自動機(jī)對土地利用未來演變進(jìn)行預(yù)測提供幫助.

在以后的研究中，將繼續(xù)為元胞自動機(jī)與軌跡分析在不同區(qū)域聯(lián)合應(yīng)用進(jìn)行土地利用時空動態(tài)模擬提供實證，對基于軌跡分析的元胞演變準(zhǔn)則的設(shè)置進(jìn)行深入研究，使其更具科學(xué)性，能夠準(zhǔn)確地描述土地利用變化規(guī)律，并根據(jù)該準(zhǔn)則進(jìn)行未來土地利用變化情況的預(yù)測，以得到更具科學(xué)性的預(yù)測結(jié)果.

［1］劉小平，黎 夏，葉嘉安，等. 利用蟻群智能挖掘地理元胞自動機(jī)的轉(zhuǎn)換規(guī)則[J]. 中國科學(xué) D 輯：地球科學(xué)，2007，37(6)：824-834.

［2］張顯峰，崔偉宏. 集成 GIS和細(xì)胞自動機(jī)模型進(jìn)行地理時空過程模擬與預(yù)測的新方法[J]. 測繪學(xué)報，2001，30(2)：148-155.

［3］黎 夏，葉嘉安. 基于神經(jīng)網(wǎng)絡(luò)的單元自動機(jī)CA及真實和優(yōu)化的城市模擬[J]. 地理學(xué)報，2002，57(2)：159-166.

［4］陳龍泉，鄭海金. 基于 Markov-CA 的土地利用/土地覆被變化動態(tài)模型研究[J]. 測繪信息與工程，2004，29(1)：108-113.

［5］邱炳文，陳崇成. 基于多目標(biāo)決策和CA模型的土地利用變化預(yù)測模型及其應(yīng)用[J]. 地理學(xué)報，2008，63(2)：165-174.

［6］唐寬金．基于景觀生態(tài)學(xué)的土地利用時空格局變化研究——以長清區(qū)為例[D]. 濟(jì)南：山東師范大學(xué)，2009.

［7］趙建軍，張洪巖，喬志和，等. 基于CA-Markov模型的向海濕地土地覆被變化動態(tài)模擬研究[J]. 自然資源學(xué)報，2009，24(12)：2 178-2 186．

［8］鄭燕鳳. 基于 GIS的CA-MARKOV模型的土地利用變化研究——以招遠(yuǎn)市為例[D]. 泰安：山東農(nóng)業(yè)大學(xué)，2009.

［9］黃海洋. DEM修正下的CA-MARKOV模型在土地利用變化中的應(yīng)用研究[D]. 重慶：西南大學(xué)，2009．

［10］何春陽，史培軍，陳 晉，等. 基于系統(tǒng)動力學(xué)模型和元胞自動機(jī)模型的土地利用情景模型研究[J]. 中國科學(xué)D輯：地球科學(xué)，2005，35(5)：464-473．

［11］鄧祥征，劉紀(jì)遠(yuǎn). 太仆寺旗土地利用變化時空格局的動態(tài)模擬[J]. 地理研究，2004，23(2)：147-156．

［12］黎 夏，葉嘉安. 基于神經(jīng)網(wǎng)絡(luò)的元胞自動機(jī)及模擬復(fù)雜土地利用系統(tǒng)[J]. 地理研究，2005，24(1)：19-27．

［13］李宗花，葉正偉. 基于元胞自動機(jī)的洪澤湖洪水蔓延模型研究[J]. 計算機(jī)應(yīng)用，2007(3)：718-720.

［14］KENNEDY R E，COHEN W B， SCHROEDER T A.Trajectory-based change detection for automated characterization of forest disturbance dynamics[J]. Remote Sensing of Environment，2007，110：370-386.

［15］ORLOCI László. Multi-scale trajectory analysis ：powerful conceptual tool for understanding ecological change[J]. Front Biol China，2009，4(2)：158-179．

［16］DOMON G，BOUCHARD A. The landscape history of Godmanchester(Quebec，Canada)，two centuries of shifting relationships between anthropic and biophysical factors[J]. Landscape ecology，2007，22：1,201-1,214.

［17］RUIZ Julie，DOMON G . Analysis of landscape pattern change trajectories within areas of intensive agricultural use：case study in a watershed of southern Quebec，Canada[J]. ,Landscape ecology，2009，24：419-432.

［18］HAINES-YOUNG R. Landscape pattern：context and process. In：Wiens JA，Moss MR(eds)Issues and perspectives in landscape ecology[M]. Cambridge：Cambridge University Press，2005：103-111.

［19］周成虎，孫戰(zhàn)利，謝一春. 地理元胞自動機(jī)研究[M].北京：科學(xué)出版社，1999.

Primary Exploration of Trajectory Analysis and Cellular Automaton in the Dynamic Simulation of Land Use Changes

WANG Dong-chuan，ZHANG Li-hui
(Department of Management Engineering，TIUC，Tianjin 300384，China)

Cellular automaton model is a dynamic network model, which is discrete in time, space and state, and local in spatial interaction and temporal causality. In the research of this model, it is often combined with other types of models to improve the simulation efficiency. The trajectory analysis method is a new method for analyzing dynamic changes in the time series. After summarizing and analyzing the characteristics of cellular automaton model and trajectory analysis method,this paper discusses the possibility of joint application of the two methods in the simulation of land use changes, laying a theoretical foundation for the establishment of land use dynamic model. The feasibility of the combination of the two methods is further confirmed by a preliminary simulation of land use changes in Binhai New Area in Tianjin.

trajectory analysis；cellular automaton；dynamic changes of land use；simulation

F301

A

1006-6853(2011)02-0135-05

2011-03-25；

2011-04-10

國家科技支撐計劃項目(2008BAJ08B16)；國家自然科學(xué)基金(70803032)；天津市哲學(xué)社會科學(xué)研究規(guī)劃項目(TJGL10-922)；天津市高等學(xué)?？萍及l(fā)展基金計劃項目(20090819)

汪東川（1972—），男，河北唐山人，天津城市建設(shè)學(xué)院講師，博士.