李繼紅 甘依童 李文慧 貢雪紅 趙宇庭 孟康

摘要：

以哈爾濱市市區(qū)為研究對(duì)象，通過(guò)構(gòu)建基于多特征的面向?qū)ο蟮姆诸惙椒?，?duì)區(qū)域2000、2010、2015三期Landsat遙感數(shù)據(jù)進(jìn)行解譯，將土地利用類型結(jié)合分辨率特性劃分為居民地、水體濕地、林草用地、農(nóng)田和未利用地五個(gè)類別。建模時(shí)，利用構(gòu)建的坡度、坡向、交通、服務(wù)設(shè)施、河流水域和人口等多個(gè)因子，通過(guò)Markov模型計(jì)算旋轉(zhuǎn)面積矩陣和轉(zhuǎn)化概率矩陣，采用ACO（蟻群優(yōu)化算法）結(jié)合類型轉(zhuǎn)化概率矩陣構(gòu)建CA（元胞自動(dòng)機(jī)）模型土地利用類型轉(zhuǎn)化規(guī)則。然后利用土地規(guī)劃綱要，采用matlab平臺(tái)對(duì)CA-MAS（多智能體）模型進(jìn)行集成，進(jìn)而對(duì)CA模型進(jìn)行修正，最終得到模擬結(jié)果，實(shí)現(xiàn)模型耦合。最后以2015年基準(zhǔn)數(shù)據(jù)，對(duì)2020年城市空間格局進(jìn)行模擬分析。耦合模型總體精度90.88%，kappa系數(shù)0.887，一致性良好，各個(gè)地類的也達(dá)到較好模擬精度。模型總體精度、各個(gè)用戶類別精度、模型的一致性和模型具體參量的模擬上均表現(xiàn)突出，表明CA-MAS耦合模型在城市空間格局演化模擬中效果良好，具有一定的應(yīng)用潛力。

關(guān)鍵詞：

面向?qū)ο蠓诸?；CA模型；MAS模型；城市土地利用；演化模擬

中圖分類號(hào)：TU 984；P 209文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1001-005X（2018）01-0030-06

Abstract：

This paper took Harbin city as the research object，by constructing the objectoriented classification method based on multi features，combining the characteristics of resolution to interpret the Landsat remote sensing data of 2000，2010 and 2015，which were divided into residential land，wetland，water and land，farmland and unused five categories land use types individually.When modeling，utilizing multiple factors which were constructed such as slope，aspect，traffic，service facilities，water，population and so on，using Markov model to calculate rotation matrix and area conversion probability matrix，and using ACO（Ant Colony Optimization）combined with the rule of transforming type conversion probability matrix to construct the CA model of land use type conversion rules.Next，using land planning outline，the CA-MAS model was integrated to modify CA model by the MATLAB platform.Finally，based on the datum data of 2015，the urban spatial pattern in 2020 was simulated and analyzed.The overall accuracy of the coupling model was 90.88%，and the kappa coefficient was 0.887.The consistency was effective，and the various land types also achieved good simulation accuracy.Model overall accuracy，each user category accuracy，consistency and model specific parameters of the model on the simulation were prominent，which showed that the CA-MAS coupled model evolution simulation had good effect in city spatial pattern and a certain potential of application.

Keywords：

Objectoriented classification；CA model；MAS model；urban land use；evolution simulation

0引言

近年來(lái)，隨著我國(guó)城市化進(jìn)程的不斷加快，各城市用地規(guī)模不斷擴(kuò)張，城市用地類型持續(xù)變化，以城市地理空間格局演化為核心的動(dòng)態(tài)模擬分析已成為地理學(xué)研究的重要內(nèi)容[1]。城市空間格局變化受到自然、社會(huì)、經(jīng)濟(jì)、文化、政治和法律等多種因素的制約，其演化過(guò)程具有高度的復(fù)雜性，使得傳統(tǒng)的城市格局演化方法受到前所未有的挑戰(zhàn)。

隨著城市微觀研究的深入，“自下而上”的建模思想逐步受到推崇，催生了智能化的城市空間模擬和預(yù)測(cè)方法。其中常用的模擬模型包括基于元胞自動(dòng)機(jī)模型（Cellular Automata，CA）、多智能體模型（Multi-agent Simulation，MAS）等[2]。國(guó)內(nèi)對(duì)城市擴(kuò)張的研究還主要集中在華南、華北地區(qū)，對(duì)于東北地區(qū)大型城市以及中型城市的研究相對(duì)較少[3-7]。為擴(kuò)展對(duì)東北大中型城市演化的研究，同時(shí)研究和驗(yàn)證CA-MAS耦合模型的有效性。本研究以哈爾濱市為研究對(duì)象，對(duì)其時(shí)空變化進(jìn)行模擬。

1元胞自動(dòng)機(jī)（CA）及多智能體（MAS）的分析

1.1元胞自動(dòng)機(jī)

CA模型在地理學(xué)領(lǐng)域的發(fā)展源于Tobler，他認(rèn)為類似元胞自動(dòng)機(jī)的地理模型的采用是分析模擬地理動(dòng)態(tài)現(xiàn)象的一次方法革命[8]。Couclelis也較早對(duì)元胞自動(dòng)機(jī)在城市擴(kuò)展方面的應(yīng)用潛力從理論上作了充分的論述，認(rèn)為在城市發(fā)展政策和城市發(fā)展模擬中，不確定性的特點(diǎn)決定了要用這類模型來(lái)進(jìn)行建模，尤其是對(duì)元胞自動(dòng)機(jī)模擬城市擴(kuò)散的闡述，對(duì)后來(lái)這方面的研究有深遠(yuǎn)的影響[9-10]。國(guó)內(nèi)學(xué)者黎夏、周成虎、何春陽(yáng)、龍瀛、曾輝等[11-18]應(yīng)用于不同的地區(qū)，做出了大量模擬研究。

元胞自動(dòng)機(jī)的模擬反映的主要是局部地區(qū)的繁衍、再生狀況[19]。首先與傳統(tǒng)的物理模型和統(tǒng)計(jì)模型相比CA模型不具有嚴(yán)格的模型函數(shù)表達(dá)式，其主要構(gòu)成包括元胞單元（cell）、胞體狀態(tài)（state）、領(lǐng)域（neighborhood）、轉(zhuǎn)化規(guī)則（transfer function）和時(shí)間（temporal）五個(gè)組成部分。元胞為CA模型的最基本運(yùn)算單元，胞體狀態(tài)為元胞體所具有的屬性信息的集合，胞體空間為所有胞體所在的空間，鄰域?yàn)榘w周圍的胞體空間，轉(zhuǎn)化規(guī)則是CA模型建模的核心。CA模型的建模規(guī)則的基本思想是只考慮某一時(shí)刻t元胞體狀態(tài)和領(lǐng)域內(nèi)元胞體狀態(tài)對(duì)下一時(shí)刻胞體的影響。其函數(shù)的表達(dá)為：

1.2CA-MAS模型總體架構(gòu)

本研究中采用30 m×30 m空間分辨率的柵格單元為一個(gè)元胞體。胞體的狀態(tài)屬性為地類的類別：居民地、水體、林草用地、農(nóng)田和未利用用地。采用二維四方形排列的歐幾里得空間，使用摩爾型領(lǐng)域（8領(lǐng)域），其信息的交互方向?yàn)?個(gè)方位，相比馮諾依曼型從4個(gè)方位轉(zhuǎn)化，從而提高的信息傳輸速度。根據(jù)自下而上的建模思想，采用ACO算法，利用數(shù)字高程、交通等屬性因子圖層建立局部的轉(zhuǎn)化規(guī)則；借助Markov模型進(jìn)行土地利用類型旋轉(zhuǎn)矩陣的計(jì)算，確定全局范圍內(nèi)的轉(zhuǎn)化規(guī)則。最后將局部與全局轉(zhuǎn)化規(guī)則相結(jié)合，并轉(zhuǎn)化為MAS系統(tǒng)。哈爾濱的規(guī)劃格局主要分為一主三副，以哈爾濱主城區(qū)為中心，阿城、雙城和肇東為副，進(jìn)行集群發(fā)展。工業(yè)發(fā)展集中為著力發(fā)展三大動(dòng)力、平方工業(yè)區(qū)，阿城裝備機(jī)械經(jīng)濟(jì)開(kāi)發(fā)區(qū)。生態(tài)規(guī)劃以山脈、河流、交通道路為骨架，以松花江為軸江景帶。將該規(guī)劃主題思想嵌入MAS模型的構(gòu)建中，遵循CA模型的轉(zhuǎn)化規(guī)則，哈爾濱主城區(qū)、阿城、雙城為中心居民地優(yōu)先轉(zhuǎn)化，河流及其山脈區(qū)域作用林草生態(tài)用地優(yōu)先轉(zhuǎn)化。在各地類轉(zhuǎn)化總量上遵循土地旋轉(zhuǎn)矩陣的預(yù)測(cè)面積限制。按照MAS系統(tǒng)行為規(guī)則根據(jù)起始年份解譯地類分布信息進(jìn)行模擬，利用土地規(guī)劃綱要同時(shí)采用CA-MAS模型集成的方式對(duì)CA模型進(jìn)行修正，實(shí)現(xiàn)模型的耦合。進(jìn)而預(yù)測(cè)哈爾濱市未來(lái)土地利用類型分布情況，如圖1所示。

2CA-MAS模擬城市土地利用類型變化

2.1研究區(qū)域

哈爾濱市地處三江平原與松嫩平原腹地，125°42′E～130°10′E，44°04′～46°40′N。由于地理位置的因素，使得位于北半球中溫帶的哈爾濱市大陸性季風(fēng)氣候特征顯著，表現(xiàn)為雨熱同期，夏季均溫23 ℃，且時(shí)間較短；冬季均溫-15 ℃，時(shí)間相對(duì)漫長(zhǎng)。年平均降水量為570 mm。平均海拔為160 m，地勢(shì)為中、西部較平坦，東部多丘陵和山地。哈爾濱市市轄9區(qū)7縣。本研究中是以哈爾濱市的9個(gè)區(qū)作為研究對(duì)象，研究面積超過(guò)7 000 km2。

2.2數(shù)據(jù)

本研究中除利用自然屬性數(shù)據(jù)以外，還需要用到諸多社會(huì)屬性數(shù)據(jù)。這些基礎(chǔ)數(shù)據(jù)包括：三期（2000年、2010年和2015年）Lansat系列遙感數(shù)據(jù)（30 m空間分辨率），見(jiàn)表1。在數(shù)據(jù)的篩選中，盡可能選擇無(wú)云或少云的影像，數(shù)據(jù)時(shí)相選擇地表景觀差異較大的9月份，有利于地表覆蓋類型的解譯；研究區(qū)1∶50 000地形圖用于各圖層數(shù)據(jù)的配準(zhǔn)；研究區(qū)高程數(shù)據(jù)GDEM-V2（30 m空間分辨率）以及使用DEM所生成的坡度、坡向數(shù)據(jù)；區(qū)域內(nèi)交通圖用于提取主要鐵路、高速、國(guó)道等信息；區(qū)域規(guī)劃圖用于居民點(diǎn)、車站和大型服務(wù)設(shè)施等重要信息提取；區(qū)域內(nèi)人口分布數(shù)據(jù)，通過(guò)國(guó)家人口普查數(shù)據(jù)獲取。

2.3模型

通過(guò)Markov模型分布計(jì)算2000年-2010年和2010-2015年旋轉(zhuǎn)面積矩陣和轉(zhuǎn)化概率矩陣；采用ACO（蟻群優(yōu)化算法）結(jié)合類型轉(zhuǎn)化概率矩陣構(gòu)建CA模型土地利用類型轉(zhuǎn)化規(guī)則；利用土地規(guī)劃綱要同時(shí)采用CA-MAS模型集成的方式對(duì)CA模型進(jìn)行修正，實(shí)現(xiàn)模型的耦合。

2.3.1基于ACO算法各屬性圖層局部轉(zhuǎn)化規(guī)則

采用ACO（蟻群優(yōu)化算法思想）針對(duì)胞體的自然屬性（高程、坡度、坡向和水域）與社會(huì)屬性（交通、服務(wù)設(shè)施、人口）建立局部轉(zhuǎn)化規(guī)則。其中一個(gè)空間轉(zhuǎn)化類型為蟻群的一條路徑（農(nóng)田-居民地）。分別通過(guò)以上自然屬性和社會(huì)因子進(jìn)行路徑加權(quán)。通過(guò)按類別分區(qū)統(tǒng)計(jì)確定各個(gè)因子的判別閾值。采用參考均值、標(biāo)準(zhǔn)差和最終設(shè)定判別閾值，利用閾值依次對(duì)多個(gè)因子進(jìn)行判別。各個(gè)目標(biāo)地類在某個(gè)因子的閾值范圍內(nèi)則地類權(quán)重加一。通過(guò)各個(gè)因子的依次篩選，并將各個(gè)地類間的轉(zhuǎn)化權(quán)重與對(duì)應(yīng)的轉(zhuǎn)化概率相乘，值最大的路徑為最終轉(zhuǎn)化路徑。

2.3.2基于Markov模型土地旋轉(zhuǎn)矩陣

通過(guò)對(duì)2000年、2010年、2015年遙感數(shù)據(jù)進(jìn)行解譯得到區(qū)域地類分布圖。為了解區(qū)域土地利用類型流向變化情況，分別計(jì)算2000年-2010年和2010年-2015年土地利用類型旋轉(zhuǎn)變化矩陣。利用旋轉(zhuǎn)矩陣可以對(duì)土地類型的轉(zhuǎn)移方向、流向、流向比例進(jìn)行分析。

利用對(duì)兩期遙感解譯結(jié)果根據(jù)以下公式進(jìn)行波段運(yùn)算：

從定量的角度，可以利用土地類型旋轉(zhuǎn)矩陣推算出土地類型轉(zhuǎn)化的概率函數(shù)。t時(shí)間內(nèi)A地類轉(zhuǎn)化為t+k時(shí)間內(nèi)B地類的面積值等于A的面積與轉(zhuǎn)化概率的乘積。其公式為：

利用概率轉(zhuǎn)化函數(shù)可以從各地類面積總量上進(jìn)行預(yù)測(cè)，確定了在宏觀上的地類的轉(zhuǎn)化規(guī)則。

2.3.3CA與MAS集成

根據(jù)CA模型的的局部、全局轉(zhuǎn)化規(guī)則結(jié)合哈爾濱城鄉(xiāng)2005-2020規(guī)劃綱要，利用MAS模型將政府決策加入模型的運(yùn)算。每一個(gè)元胞體視為一個(gè)MAS系統(tǒng)，元胞的鄰域視為MAS系統(tǒng)的外部環(huán)境。MAS系統(tǒng)根據(jù)CA模型的轉(zhuǎn)化結(jié)果進(jìn)行調(diào)整，使轉(zhuǎn)化的結(jié)果符合區(qū)域規(guī)劃決策和轉(zhuǎn)化矩陣面積總量。

哈爾濱的規(guī)劃格局主要分為一主三副，以哈爾濱主城區(qū)為中心，阿城、雙城和肇東為副，進(jìn)行集群發(fā)展。工業(yè)發(fā)展集中為著力發(fā)展三大動(dòng)力、平方工業(yè)區(qū)，阿城裝備機(jī)械經(jīng)濟(jì)開(kāi)發(fā)區(qū)。生態(tài)規(guī)劃以山脈、河流、交通道路為骨架，以松花江為軸江景帶。

將以上規(guī)劃主題思想嵌入MAS模型的構(gòu)建中。遵循CA模型的轉(zhuǎn)化規(guī)則，哈爾濱主城區(qū)、阿城、雙城為中心在調(diào)整轉(zhuǎn)化階段居民地優(yōu)先轉(zhuǎn)化，河流及其山脈區(qū)域采用林草生態(tài)用地優(yōu)先轉(zhuǎn)化。在各地類轉(zhuǎn)化總量上遵循土地旋轉(zhuǎn)矩陣的預(yù)測(cè)面積限制。

得到轉(zhuǎn)化結(jié)果后與初始地類計(jì)算得到轉(zhuǎn)化后的旋轉(zhuǎn)矩陣XZ_conversion，并將之與預(yù)測(cè)旋轉(zhuǎn)矩陣XZ_forecast進(jìn)行比對(duì)。根據(jù)比對(duì)結(jié)果利用決策信息建立的緩沖區(qū)，在緩沖區(qū)內(nèi)根據(jù)政策進(jìn)行引導(dǎo)轉(zhuǎn)化。緩沖區(qū)外區(qū)域，將轉(zhuǎn)化結(jié)果帶入轉(zhuǎn)化函數(shù)進(jìn)行迭代，直到接近預(yù)測(cè)旋轉(zhuǎn)矩陣結(jié)束迭代。

2.4技術(shù)路線

CA-MAS橫型模擬預(yù)測(cè)城市用地?cái)U(kuò)張流程如圖2所示。

3結(jié)果與分析

3.1解譯結(jié)果

研究中以2015年7月的實(shí)地調(diào)查樣本點(diǎn)67個(gè)野外實(shí)地調(diào)查樣本補(bǔ)充點(diǎn)以及從Google Earth上隨機(jī)生成的100個(gè)點(diǎn)作為精度評(píng)價(jià)的基礎(chǔ)數(shù)據(jù)（基本符合均勻分布），其中野外調(diào)查樣本已經(jīng)明確其各自地類，隨機(jī)生成的樣本類型通過(guò)Google Earth和現(xiàn)有的試驗(yàn)區(qū)土地利用現(xiàn)狀圖中目視確定。精度見(jiàn)表2，2000年、2010年、2015年研究區(qū)域遙感圖像分類結(jié)果如圖3所示。

3.2模型評(píng)價(jià)

利用2015年遙感數(shù)據(jù)解譯結(jié)果對(duì)2015年模型模擬結(jié)果進(jìn)行檢驗(yàn)。精度檢驗(yàn)采用混淆矩陣、Kappa系數(shù)進(jìn)行評(píng)定?；煜仃囀且环N精度評(píng)定常用方式。其對(duì)角線數(shù)值為各地類正確模擬的像元數(shù)，其余位置為某類別模擬結(jié)果錯(cuò)分到其他類別的像元數(shù)。利用2015年模擬結(jié)果和2015年解譯結(jié)果構(gòu)建混淆矩陣（表3）與獲得模型用戶精度（表4）。通過(guò)計(jì)算模型的總體精度為90.88%。模型模擬總體精度較高，對(duì)于居民地、林草地和農(nóng)田具有很好的擬合性。但在水體濕地和未利用用地模擬上精度稍低。因?yàn)樗w濕地受年降雨量、氣候條件等諸多自然因素的影響年變化較大，而研究中又沒(méi)有將氣象數(shù)據(jù)作為因子考慮在內(nèi)，因此不具有顯著的規(guī)律性。對(duì)于未利用用地而言，從分類解譯的角度來(lái)看，誤差較大，易于與其他地類混淆。其次，未利用用地總體數(shù)量較小，一定數(shù)量上的變更對(duì)于總體比例的變化影響較大，加上未利用用地受政策和土地規(guī)劃主觀性影響較大，隨著分類誤差和模型誤差的傳遞和累加造成未利用用地模型模擬精度較差。經(jīng)過(guò)計(jì)算2015年模擬結(jié)果和分類結(jié)果的卡帕系數(shù)為0.887。通過(guò)一致性檢驗(yàn)，可知模型模擬結(jié)果與實(shí)際數(shù)據(jù)具有較好的一致性。

從數(shù)據(jù)可以得出：伴隨著城鎮(zhèn)化的進(jìn)程農(nóng)田不斷減少，居民地城鎮(zhèn)用地隨之增加。隨著環(huán)境意識(shí)的增強(qiáng)，水體濕地，林地草地等生態(tài)用地，呈現(xiàn)增加態(tài)勢(shì)。隨著規(guī)劃用地日趨合理化，未利用用地呈現(xiàn)快速減少態(tài)勢(shì)。根據(jù)統(tǒng)計(jì)結(jié)果無(wú)論從各個(gè)類型的用地的變化趨勢(shì)上還是各地類實(shí)際變化率上模型預(yù)測(cè)結(jié)果與實(shí)際較為接近，居民地、水體、林地和農(nóng)田地類的土地變化率誤差均控制在1%以內(nèi)，模型在未利用地類的模擬精度稍差。

通過(guò)2010年的基點(diǎn)數(shù)據(jù)運(yùn)行得到2015年的模擬結(jié)果，經(jīng)檢驗(yàn)?zāi)Ｐ涂傮w精度較高，各個(gè)地類的增長(zhǎng)率也達(dá)到較好模擬精度，模擬效果理想。之后對(duì)2020年哈爾濱市的發(fā)展情況進(jìn)行模擬，得到模擬結(jié)果如圖4所示，對(duì)模擬結(jié)果進(jìn)行統(tǒng)計(jì)分析，得到統(tǒng)計(jì)數(shù)據(jù)見(jiàn)表7。

4結(jié)論

采用CA-MAS集成方式利用土地規(guī)劃信息主動(dòng)干預(yù)土地轉(zhuǎn)化的方向。有效彌補(bǔ)了CA模型未考慮到人為政策因素主動(dòng)干預(yù)所產(chǎn)生的影響，從而提高模型模擬精度。模型總體精度各個(gè)用戶類別精度、模型的一致性和模型具體參量的模擬上均表現(xiàn)突出，表明CA-MAS耦合模型在城市擴(kuò)展模擬上具有一定的應(yīng)用潛力。

從模型對(duì)2020年預(yù)測(cè)結(jié)果可知，2015-2020年哈爾濱居民用地繼續(xù)增長(zhǎng)，預(yù)計(jì)增長(zhǎng)40.256 km2，水體濕地和林地也呈現(xiàn)出增長(zhǎng)態(tài)勢(shì)，農(nóng)田和未利用用地繼續(xù)減少，但總體趨勢(shì)呈現(xiàn)出變化減緩趨勢(shì)。在一定程度上表明研究區(qū)城鎮(zhèn)化進(jìn)程減緩。在空間中呈現(xiàn)出哈爾濱主城區(qū)周邊至雙城、阿城、肇東沿道路發(fā)展由線到面擴(kuò)展。小城鎮(zhèn)擴(kuò)展較明顯，全區(qū)小城鎮(zhèn)呈現(xiàn)出環(huán)狀擴(kuò)展的現(xiàn)象。在平房、阿城經(jīng)濟(jì)開(kāi)發(fā)區(qū)、雙城、松北食品加工工業(yè)園區(qū)呈現(xiàn)出居民用地聚集明顯。林草生態(tài)用地沿松花江、阿什河和呼蘭河流域呈現(xiàn)增加態(tài)勢(shì)。在哈雙、哈阿、哈五、哈同、哈伊哈大等高速公路兩側(cè)出現(xiàn)小斑塊林草用地增長(zhǎng)。模擬結(jié)果與哈爾濱2005-2020規(guī)劃綱要發(fā)展方向相符，綜合注重經(jīng)濟(jì)區(qū)的快速發(fā)展，小城鎮(zhèn)的城鎮(zhèn)化建設(shè)和綠水藍(lán)天的生態(tài)發(fā)展方向。

【參考文獻(xiàn)】

[1]胡曉明，李月臣，黃孝艷，等.城市空間擴(kuò)展研究及進(jìn)展[J].現(xiàn)代城市研究，2013，28（6）：60-82.

[2]閆梅，黃金川.國(guó)內(nèi)外城市空間擴(kuò)展研究評(píng)析[J].地理科學(xué)進(jìn)展，2013，32（7）：1039-1050.

[3]樸妍，馬克明.北京城市建成區(qū)擴(kuò)張的經(jīng)濟(jì)驅(qū)動(dòng)：1978-2002[J].中國(guó)國(guó)土資源經(jīng)濟(jì)，2006，20（7）：33-37.

[4]李曉文，方精云，樸世龍.上海及周邊主要城鎮(zhèn)用地?cái)U(kuò)展空間特征及其比較[J].地理研究，2003，22（6）：769-779.

[5]李開(kāi)宇，曹天艷.大都市邊緣地區(qū)城市化與城市空間擴(kuò)展研究——以廣州市番禺區(qū)為例[J].西北大學(xué)學(xué)報(bào)（自然科學(xué)版），2014，40（3）：523-526.

[6]吝濤，李新虎，張國(guó)欽，等.廈門島城市空間擴(kuò)張?zhí)卣骷捌溆绊懸蛩胤治鯷J].地理學(xué)報(bào)，2010，65（6）：715-726.

[7]張恒.基于MODIS遙感影像及元胞自動(dòng)機(jī)的京津冀地區(qū)城市模擬研究[D].北京：中國(guó)地質(zhì)大學(xué)，2009.

[8]Tobler W R.A computer movie simulating urban growth in the Detroit region[J].Economic Geography，1970，46（1）：234-240.

[9]Couclelis H M.Cellular Worlds：A Framework for Modeling MicroMacro Dynamics[J].Environment & Planning A，1985，17（5）：585-596.

[10]Couclelis H.Macrostructure and microbehavior in a metropolitan area[J].Environment & Planning B Planning & Design，1989，16（2）：141-154.

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

[12]黎夏，葉嘉安.約束性單元自動(dòng)演化CA模型及可持續(xù)城市發(fā)展形態(tài)的模擬[J].地理學(xué)報(bào)，1999，54（4）：159-166.

[13]陶海燕，黎夏，陳曉翔，等.基于多智能體的地理空間分異現(xiàn)象模擬——以城市居住空間演變?yōu)槔齕J].地理學(xué)報(bào)，2007，62（6）：579-588.

[14]劉小平，黎夏，葉嘉安.基于多智能體系統(tǒng)的空間決策行為及土地利用格局演變的模擬[J].地球科學(xué)，2006，36（11）：1027-1036.

[15]楊青生，黎夏.多智能體與元胞自動(dòng)機(jī)結(jié)合及城市用地?cái)U(kuò)張模擬[J].地理科學(xué)，2007，27（4）：542-548.

[16]何春陽(yáng)，陳晉，史培軍，等.基于CA的城市空間動(dòng)態(tài)模型研究[J].地球科學(xué)進(jìn)展，2002，17（2）：188-195.

[17]龍瀛，韓昊英，毛其智.利用約束性 CA 制定城市增長(zhǎng)邊界[J].地理學(xué)報(bào)，2009，64（8）：999-1008.

[18]曾輝.深圳市龍華地區(qū)城鎮(zhèn)用地動(dòng)態(tài)模型建設(shè)及模擬研究[J].生態(tài)學(xué)報(bào)，2000，20（4）：545-551.

[19]楊青生，黎夏，劉小平.基于Agent和CA的城市土地利用變化研究[J].地球信息科學(xué)，2005，27（2）：78-81.