曹偉超，魯小丫，馬月偉 ，文學虎，王 蕾，陳 陽

(1.四川省第三測繪工程院，四川 成都 610500；2.西南民族大學 計算機科學與技術(shù)學院，四川 成都 610041；3.西南林業(yè)大學 生態(tài)旅游學院，云南 昆明 650224)

基于多時相夜間燈光數(shù)據(jù)的城市建成區(qū)提取研究
——以西南地區(qū)為例

曹偉超1，魯小丫2，馬月偉3，文學虎1，王 蕾1，陳 陽1

(1.四川省第三測繪工程院，四川 成都 610500；2.西南民族大學 計算機科學與技術(shù)學院，四川 成都 610041；3.西南林業(yè)大學 生態(tài)旅游學院，云南 昆明 650224)

城市化對地區(qū)經(jīng)濟的發(fā)展具有重要的戰(zhàn)略意義，而地理信息系統(tǒng)(GIS)和遙感技術(shù)的應(yīng)用，為城市化發(fā)展信息的獲取提供新的研究途徑。文中以2003年、2005年和2007年3期的DMSP/OLS夜間燈光數(shù)據(jù)為例，借助地理信息系統(tǒng)(GIS)平臺，提取出西南地區(qū)城市空間信息?？紤]區(qū)域差異性，將西南地區(qū)分為川渝、云貴和廣西3個區(qū)域，分別進行城市化發(fā)展過程和趨勢的分析。將夜間燈光數(shù)據(jù)應(yīng)用到經(jīng)濟不發(fā)達的地區(qū)，在數(shù)據(jù)提取精度上還需進一步結(jié)合中尺度遙感影像進行分析。

DMSP/OLS；地理信息系統(tǒng)；城市化；西南地區(qū)

城市化指“人類生產(chǎn)與生活方式由鄉(xiāng)村型向城市型轉(zhuǎn)化的歷史過程，表現(xiàn)為鄉(xiāng)村人口向城市人口轉(zhuǎn)化以及城市不斷發(fā)展和完善的過程”[1]。因此一個區(qū)域城市化的發(fā)展，不僅是該區(qū)域人口城市化過程，還包括經(jīng)濟、社會、生活方式等方面。改革開放以來，中國正經(jīng)歷著高速的城市化發(fā)展，隨著中國城市化進程的深入，大都市區(qū)和大城市群已成為目前中國城市化進程中最引人注目和最有活力的部分。

在快速城市化進程中，我國西南地區(qū)的城市化卻表現(xiàn)出整體水平相對低下、空間分布不均衡等特征，各省份大、中城市以及集中于此的城鎮(zhèn)人口與經(jīng)濟力量，成為其城市化和經(jīng)濟社會發(fā)展的核心區(qū)域和主要力量；而核心區(qū)域之外的廣大地區(qū)，城鎮(zhèn)稀少，自然條件惡劣，經(jīng)濟社會發(fā)展水平遠低于中心區(qū)域[2]。根據(jù)西南地區(qū)城市化的發(fā)展現(xiàn)狀，需要以新的理念、新的思路，在依靠中心區(qū)域吸引各種要素的同時，指導邊緣區(qū)域的城鎮(zhèn)建設(shè)，進而實現(xiàn)整個區(qū)域經(jīng)濟社會的全面發(fā)展。

GIS和遙感技術(shù)的發(fā)展為區(qū)域經(jīng)濟信息的獲取提供新的研究途徑。美國軍事氣象衛(wèi)星(Defense Meteorological Satellite Program，DMSP)搭載的(Operational Line scan System，OLS)傳感器為大尺度的城市研究提供一種新的數(shù)據(jù)獲取手段，DMSP/OLS傳感器可在夜間工作，能夠探測到城市燈光甚至小規(guī)模居民地、車流等發(fā)出的低強度燈光，比較適合大尺度城市化進程的動態(tài)監(jiān)測[3]。國外學者主要利用該數(shù)據(jù)分析與社會經(jīng)濟因子的線性相關(guān)性，這些因子包括人口密度[4]、城市面積擴張[5]以及電力能源消耗[6]等；國內(nèi)主要以卓莉等人為代表，利用夜間燈光數(shù)據(jù)對人口密度、城市化過程等進行模擬[7-9]，利用多時相數(shù)據(jù)之間的關(guān)系對未來城市的發(fā)展狀況進行估算預(yù)測，應(yīng)用前景非常廣闊。

中國利用夜間燈光數(shù)據(jù)在區(qū)域尺度對城市化發(fā)展過程的研究多以經(jīng)濟較發(fā)達區(qū)為主，對于經(jīng)濟欠發(fā)達地區(qū)涉及較少。本文以西南地區(qū)為例，借助GIS平臺，基于多時相的夜間燈光數(shù)據(jù)，獲取2003年、2005年以及2007年的城市空間分布信息。

1 研究區(qū)域和數(shù)據(jù)

本文以西南地區(qū)為研究區(qū)，包括四川、重慶、云南、貴州和廣西5個行政單元。作為一個具有很強的地理、經(jīng)濟、人文特點的區(qū)域，該研究區(qū)在經(jīng)濟和地理上構(gòu)成一個整體，同時在城市化模式上也相對一致，是中國城市化中一個獨特區(qū)域[4]。研究區(qū)面積137萬km，約占全國總面積的14.27%， 2007年末總?cè)丝谶_23 987萬人，其中城鎮(zhèn)人口有8470萬人，城鎮(zhèn)化率為35.31%，低于中國44.94%的平均水平，區(qū)域城市化水平亟待提高。

夜間燈光數(shù)據(jù)是美國國家地球物理數(shù)據(jù)中心提供的DMSP/OLS夜間燈光遙感影像，城鎮(zhèn)建設(shè)用地、人口以及土地面積等統(tǒng)計數(shù)據(jù)來源于2004-2008年《中國統(tǒng)計年鑒》；輔助數(shù)據(jù)主要包括來自國家基礎(chǔ)地理信息中心的中國1:400萬的全國行政區(qū)劃圖、交通線路圖，用于驗證夜間燈光數(shù)據(jù)反映城市化過程及其特征的合理性。

2 研究方法

2.1 數(shù)據(jù)預(yù)處理

1)投影轉(zhuǎn)換。獲取的原始夜間燈光數(shù)據(jù)是WGS84坐標系統(tǒng)的數(shù)據(jù)，考慮到數(shù)據(jù)分析和處理的需要，坐標系統(tǒng)轉(zhuǎn)換為Albers投影，D_Krasovsky_1940橢球體。

2)誤差消除。結(jié)合2003-2007年中國城市化發(fā)展的特點，假設(shè)后期城鎮(zhèn)建設(shè)用地均以前期為基礎(chǔ)，沒有前期存在而后期消失的城鎮(zhèn)建設(shè)用地；同時一些多云的地區(qū)也常表現(xiàn)為高灰度值區(qū)域，這些誤差都應(yīng)該消除。結(jié)合城鎮(zhèn)居民點數(shù)據(jù)，通過GIS柵格運算，可以較好地消除該兩項誤差，提高數(shù)據(jù)精度。

2.2 信息提取方法

提取燈光數(shù)據(jù)前，再次假設(shè)現(xiàn)有土地利用的統(tǒng)計數(shù)據(jù)可以基本反映中國城鎮(zhèn)用地的總量特征，利用DMSP/OLS數(shù)據(jù)提取的各省區(qū)城鎮(zhèn)用地的總量應(yīng)充分接近統(tǒng)計數(shù)據(jù)[8]。基于第二種假設(shè)，對預(yù)處理過的數(shù)據(jù)利用二分法確定出各行政單元燈光閾值DNT，從而提取各省區(qū)城鎮(zhèn)燈光數(shù)據(jù)。

具體方法就是設(shè)定夜間燈光亮度閾值DNT，統(tǒng)計各燈光斑塊面積與各行政單元城鎮(zhèn)建設(shè)用地面積進行比較，直到某一閾值條件下利用夜間數(shù)據(jù)提取的城鎮(zhèn)建設(shè)用地面積與統(tǒng)計數(shù)據(jù)充分接近為止，具體流程[10]如圖1所示。

圖1 城鎮(zhèn)夜間燈光閾值確定流程

圖1中：DNmax和DNmin分別是DMSP/OLS影像的最大和最小亮度值；f(DNT)為閾值DNT下提取的城鎮(zhèn)建設(shè)面積；US為實際統(tǒng)計的城鎮(zhèn)建設(shè)面積；△f(DNT)表示在閾值DNT條件下，計算得到的城鎮(zhèn)建設(shè)面積與US的差值；△f(DNT)min表示△f(DNT)與△f(DNT-1)和△f(DNT+1)的絕對值相比是否最小，如果最小，DNT便是閾值；如果不是，再比較f(DNT)與US的關(guān)系，如果大于US，便用DNT值代替DNmin，相反則用DNT代替DNmax，重新比較，直到得出最佳閾值。

2.3 城市空間信息提取

在夜間燈光數(shù)據(jù)中，燈光強度(圖像灰度值)越高的地方，城鎮(zhèn)建設(shè)用地的可能性越大。采用二分法獲取各省區(qū)的城鎮(zhèn)建設(shè)用地，將獲取的城鎮(zhèn)建設(shè)用地面積與統(tǒng)計數(shù)據(jù)進行比較，直到提取的各省城鎮(zhèn)用地總量與統(tǒng)計數(shù)據(jù)充分接近為止。最終得到2003年、2005年和2007年的城市空間信息(見圖2)，并提取最終的效果評價(見表1)。從表1中可以看出，提取的城鎮(zhèn)建設(shè)用地面積與統(tǒng)計值的相對誤差絕對值，該值除2005年在廣西超過8%，其余均在8%以內(nèi)，每年總體誤差絕對值均在1%之內(nèi)。該結(jié)果表明利用DMSP/OLS數(shù)據(jù)提取的研究區(qū)城市空間分布信息具有一定可信度，在此基礎(chǔ)上進行城市化空間過程的研究是合理的。

表1 基于DMSP/OLS的西南地區(qū)城市空間信息提取結(jié)果分析

圖2 2003-2007年西南地區(qū)城市空間格局

3 結(jié)果分析

結(jié)合提取的三期西南地區(qū)的城鎮(zhèn)建設(shè)用地空間分布柵格數(shù)據(jù)，借助GIS平臺，將3個時期用地擴展結(jié)果進行疊加，恢復(fù)2003-2007年西南地區(qū)城鎮(zhèn)建設(shè)用地的擴展過程(見圖3)。在擴展過程圖中，可以比較清晰的看到2003-2005年城鎮(zhèn)擴展區(qū)域，而代表2005-2007年城鎮(zhèn)擴展的區(qū)域則不是非常明顯。從統(tǒng)計數(shù)據(jù)(見表1)中可知：2003-2005年的城鎮(zhèn)面積增長率為3.14%，而2005-2007年的增長率只有2.57%，可見在城鎮(zhèn)面積擴展的速率上有所減慢，也證實利用夜間燈光數(shù)據(jù)分析城市化過程的可行性。

借助交通線路和行政區(qū)劃數(shù)據(jù)，以2007年夜間燈光數(shù)據(jù)為基礎(chǔ)，得到西南地區(qū)城市的空間分布情況(見圖4)。圖中表明西南地區(qū)城鎮(zhèn)擴展除以大城市為中心，如圍繞成都、重慶、昆明、貴陽、南寧省會城市逐漸發(fā)展成為西南地區(qū)的城市群外，還存在沿交通干線發(fā)展起來的多條小有規(guī)模的城市帶(如廣西的桂林-柳州-南寧-北海)，這些城市帶已成為區(qū)域城市發(fā)展的主體，城鎮(zhèn)建設(shè)用地在此呈現(xiàn)較明顯的集聚。

圖3 2003-2007年西南地區(qū)城市化過程

圖4 2007年西南地區(qū)城市空間分布

4 結(jié)束語

利用DMSP/OLS夜間燈光數(shù)據(jù)，結(jié)合統(tǒng)計數(shù)據(jù)，提取西南地區(qū)5省(自治區(qū)、市)2003-2007年的城市化空間信息。從空間分析的角度出發(fā)，根據(jù)城市化發(fā)展現(xiàn)狀，把西南地區(qū)分為川渝地區(qū)、云貴地區(qū)和廣西3大區(qū)域，系統(tǒng)地對西南地區(qū)2003-2007年的空間演化過程進行分析。基本結(jié)論如下：

1)西南地區(qū)城市化過程中，圍繞省會等特大型、大型城市以及主要交通干線，呈現(xiàn)出較明顯的集聚現(xiàn)象。集聚對于經(jīng)濟相對落后的西南地區(qū)而言，將增強城鎮(zhèn)集中區(qū)域的吸引力，促使鄉(xiāng)村人口和社會生產(chǎn)要素向區(qū)位優(yōu)越、空間可達性強、適宜居住、產(chǎn)業(yè)和城鎮(zhèn)集聚程度較高的大中城市集中，最終形成城鄉(xiāng)一體化的空間結(jié)構(gòu)體系，擴展城市化范圍。

2)川渝、云貴以及廣西三大區(qū)域呈現(xiàn)出不同的城市化過程：川渝地區(qū)憑借經(jīng)濟基礎(chǔ)和地理區(qū)位優(yōu)勢，逐漸發(fā)展成為一個以成都和重慶為中心，附帶眾多中小城鎮(zhèn)的城市群，成為西南地區(qū)城市化水平最高的區(qū)域；云貴和廣西兩個區(qū)域缺少特大型城市，總體經(jīng)濟和社會發(fā)展水平相對落后，其城市化過程主要是依托交通干線的集聚能力，吸引更多的元素在沿線聚集，形成具有一定規(guī)模的大中城市，并隨著城市的發(fā)展，其功能不斷完善，輻射能力增強，進而推動沿交通干線城市帶的發(fā)展。

借助夜間燈光數(shù)據(jù)對城市化過程的研究還有很多，本文只從城鎮(zhèn)建設(shè)用地擴展的角度對西南地區(qū)的城市化過程進行分析。城市化過程除表現(xiàn)為非城鎮(zhèn)建設(shè)用地向城鎮(zhèn)用地轉(zhuǎn)化之外，還涉及到城市內(nèi)部自我功能的完善以及城市三維立體方向上的發(fā)展等一系列的問題。

[1]王莉，衛(wèi)海燕，宋雪娟.中國大陸城市化水平的省際差異及分類研究[J].資源開發(fā)與市場，2010,26(9):782-785.

[2]柏貴喜，羅義云.西南地區(qū)城市化發(fā)展戰(zhàn)略研究[J].中南民族大學學報：人文社會科學版，2004,24(4)：37-42.

[3]陳晉，卓莉，史培軍,等.基于DMSP/OLS數(shù)據(jù)的中國城市化過程研究[J].遙感學報，2003,7(3)：168-175.

[4]SOTTON P. Modeling population density with nighttime satellite imagery and GIS [J].Computers Environment and Urban Systems, 1997, 21(3-4):227-244.

[5]IMHOFF M L, LAWRENCE W T, STUTZER D C, et al. A technique for using composite DMSP/OLS city lights satellite data to map urban area [J].Remote Sensing of Environment, 1997(61):261-370.

[6]AMARAL S, CAMARA G, MIGUEL A,et al. Estimating population and energy consumption in Brazilian Amazonia using DMSP night-time satellite data [J].Computers, Environment and Urban Systems,2005(29):179-195.

[7]宋仁波，李文慧，王細元，等.基于雙張影像的城市建筑物高度提取方法[J].測繪工程，2014,23(9)：66-69.

[8]何春陽，李景剛.基于夜間燈光數(shù)據(jù)的環(huán)渤海地區(qū)城市化過程[J].地理學報，2005,60(3)：409-417.

[9]李景剛，何陽春.基于DMSP/OLS燈光數(shù)據(jù)的快速城市化過程的生態(tài)效應(yīng)評價研究[J].遙感學報，2007,11(1)：115-126.

[10]曹麗琴,李平湘,張良培.基于DMSP/OLS夜間燈光數(shù)據(jù)的城市人口估算——以湖北省各縣市為例[J].GIS技術(shù),2009(1):83-87.

[責任編輯：張德福]

A study of extracting urban built-up area by using multitemporal night-time light data——Taking Southwest China as An Example

CAO Wei-chao1, LU Xiao-ya2，MA Yue-wei3，WEN Xue-hu1，WANG Lei1，CHEN Yang1

(1.The Third Mapping of Sichuan Academy of Engineering, Chengdu 610500, China;2. School of Computer Science and Technology, Southwest University For Nationalities,Chengdu 610041,China;3. College of Ecotourism, Southwest Forestry University, Kunming 650224, China)

The urbanization has in important strategic sense in the economic development in Southwest China, and the application of GIS and remote sensing provides a new research approach for getting the information about urbanization. With GIS platform, the urban information is derived from night-time light (DMSP/OLS) Data in 2003, 2005 and 2007 with the support of statistical data. Considering the regional differences, Southwest China is divided into three regions for analyzing the process and trend of urbanization. In addition, the data’s precision shall be improved by combining mesoscale image data in underdeveloped regions .

DMSP/OLS; Geographic Information System(GIS)；urbanization; Southwest China

2014-05-01;

2014-10-10

四川測繪地理信息局2013年科技支撐項目(J2013ZC03)；四川省地理國情監(jiān)測工程技術(shù)研究中心資助項目(GC201414;GC201413)

曹偉超(1986-)，男，碩士.

P208

：A

：1008-5696(2014)12-0048-04