基于DMSP/OLS燈光數(shù)據(jù)的川南經(jīng)濟(jì)區(qū)城鎮(zhèn)聯(lián)系研究

王俊翔， 陳一君， 郭嵐， 韓兵， 畢歡

(四川理工學(xué)院管理學(xué)院， 四川自貢643000)

引言

經(jīng)濟(jì)全球化使得城市間的競(jìng)爭(zhēng)主體由單體城市演化為區(qū)域性群體城市。目前，群體性的城市經(jīng)濟(jì)區(qū)已成為加快城鎮(zhèn)化進(jìn)程的重要引擎[1]。城鎮(zhèn)聯(lián)系是研究城市之間、區(qū)縣間的一種內(nèi)在經(jīng)濟(jì)空間聯(lián)系，從空間尺度和經(jīng)濟(jì)尺度上探尋城市發(fā)展的趨勢(shì)和內(nèi)部的聯(lián)系特點(diǎn)。研究經(jīng)濟(jì)區(qū)的城鎮(zhèn)聯(lián)系可以幫助人們更好的把握城市經(jīng)濟(jì)區(qū)的發(fā)展，構(gòu)建合理的內(nèi)部發(fā)展機(jī)制，提高群體城市的競(jìng)爭(zhēng)力，并為政府規(guī)劃和城市現(xiàn)代化提供理論指導(dǎo)。因此，研究城市經(jīng)濟(jì)區(qū)的城鎮(zhèn)聯(lián)系意義重大。近年來，經(jīng)濟(jì)社會(huì)的發(fā)展推動(dòng)著城鎮(zhèn)的迅速擴(kuò)展，城市區(qū)域的拓展一定程度上造成耕地的流失，并對(duì)于城市的可持續(xù)發(fā)展造成影響[2]。城鎮(zhèn)化的快速發(fā)展也存在著空間拓展的失控問題[3]，為了能夠?qū)?jīng)濟(jì)區(qū)域進(jìn)行有效的管控，能夠從時(shí)間尺度上客觀把握城市及城鎮(zhèn)發(fā)展的進(jìn)程也具有深遠(yuǎn)的意義。作為社會(huì)經(jīng)濟(jì)發(fā)展的引擎，城市在為居民提供利益的同時(shí)，也產(chǎn)生了許多負(fù)面效應(yīng)。因此，對(duì)目前區(qū)域城鎮(zhèn)化發(fā)展的進(jìn)程必須進(jìn)行長(zhǎng)時(shí)間深入的了解[4]。

在城鎮(zhèn)發(fā)展中，客觀把握城市發(fā)展的特征規(guī)律，是政府進(jìn)行科學(xué)規(guī)劃與制定政策的重要基礎(chǔ)。在城鎮(zhèn)發(fā)展的理論研究中，以往學(xué)者多用單一主城區(qū)用地的范圍變化來探尋城鎮(zhèn)空間的擴(kuò)展規(guī)律[5-8]。該方法雖然簡(jiǎn)便，但卻忽視了城市內(nèi)部的差異性，城市作為一個(gè)綜合的研究對(duì)象，它涵蓋了經(jīng)濟(jì)、用地、人口、環(huán)境等領(lǐng)域，僅用單一的主城區(qū)或城市用地指標(biāo)研究城鎮(zhèn)的發(fā)展有失偏頗。此外，目前的研究領(lǐng)域側(cè)重于宏觀尺度的土地利用、空間拓展、格局體系分析及模型等多方面，而對(duì)城鎮(zhèn)經(jīng)濟(jì)區(qū)空間拓展的微觀尺度研究相對(duì)比較缺乏。本文在研究過程中，用遙感影像數(shù)據(jù)的相關(guān)分析解決了單一用地變化無法描述城市內(nèi)部差異性的問題，同時(shí)從區(qū)域經(jīng)濟(jì)微觀尺度上，利用具備時(shí)間尺度的遙感影像對(duì)城市經(jīng)濟(jì)區(qū)的城鎮(zhèn)聯(lián)系及變化過程進(jìn)行了定性定量研究。

目前群體經(jīng)濟(jì)區(qū)的研究，宏觀尺度的城市經(jīng)濟(jì)區(qū)如京津冀、珠三角等，以往學(xué)者研究的都比較系統(tǒng)、透徹；而微觀尺度的城市經(jīng)濟(jì)區(qū)，其定量研究的還不夠全面。2016年“十三五規(guī)劃”關(guān)于優(yōu)化城鎮(zhèn)化布局和形態(tài)的系列政策表明：相關(guān)城市經(jīng)濟(jì)區(qū)的發(fā)展日益成為政府決策和科學(xué)界所重視的問題。近期，四川省出臺(tái)了《成渝經(jīng)濟(jì)區(qū)南部城市群發(fā)展規(guī)劃(2014—2020年)》，規(guī)劃中明確強(qiáng)調(diào)要構(gòu)建新興區(qū)域性經(jīng)濟(jì)區(qū)，尤其是川南經(jīng)濟(jì)區(qū)。因此，本文從微觀尺度選取了西南地區(qū)四川境內(nèi)的川南經(jīng)濟(jì)區(qū)作為研究范圍區(qū)。

在研究方法上，關(guān)于遙感影像數(shù)據(jù)，它能夠?qū)崿F(xiàn)區(qū)域大尺度或全球尺度的同步觀測(cè)，可以用于城鎮(zhèn)信息的快速提取和城鎮(zhèn)聯(lián)系的研究。DMSP/OLS夜間燈光影像數(shù)據(jù)是遙感影像的一種類型，DMSP(Defence Meteorological Satellite Program,DMSP)指的是美國(guó)“國(guó)防氣象衛(wèi)星計(jì)劃”，OLS(Operational Linescan System)是指DMSP搭載的傳感器掃描系統(tǒng)，主要用于燈光亮度的記錄[9]。全球范圍的DMSP/OLS夜間燈光影像，其亮斑燈光區(qū)域代表經(jīng)濟(jì)活動(dòng)較強(qiáng)的城鎮(zhèn)區(qū)域，而黑暗區(qū)域則代表經(jīng)濟(jì)活動(dòng)較弱的鄉(xiāng)村區(qū)域，可有效區(qū)分城鎮(zhèn)區(qū)和非城鎮(zhèn)區(qū)[10]。在時(shí)間尺度上，DMSP/OLS全球影像具有廣泛的時(shí)效性和經(jīng)濟(jì)實(shí)用性，且年際區(qū)域?qū)Ρ刃暂^強(qiáng)，因此廣泛用于宏觀大尺度的城市空間拓展[11]、經(jīng)濟(jì)發(fā)展[12-15]和人口估計(jì)[16-19]、城鎮(zhèn)用電與能源消耗統(tǒng)計(jì)[20]、碳排放[21]及光學(xué)污染等環(huán)境評(píng)估[22]。目前，在城市化[23-24]研究中，燈光數(shù)據(jù)能夠監(jiān)測(cè)城鎮(zhèn)化發(fā)展水平、人口密度和經(jīng)濟(jì)狀況等活動(dòng)，是人類活動(dòng)的綜合表示。國(guó)內(nèi)外研究者常用燈光數(shù)據(jù)與經(jīng)濟(jì)活動(dòng)的諸多因子進(jìn)行回歸分析，成果豐碩。DMSP/OLS數(shù)據(jù)通過區(qū)域空間的經(jīng)濟(jì)規(guī)模變化，對(duì)人類的經(jīng)濟(jì)活動(dòng)進(jìn)行預(yù)測(cè)，其應(yīng)用前景廣闊[25]。

在區(qū)域城鎮(zhèn)聯(lián)系的對(duì)象研究方面，國(guó)外相關(guān)學(xué)者[26-28]重點(diǎn)分析國(guó)家和省域級(jí)別范圍的城市等級(jí)體系及空間結(jié)構(gòu)方面的研究，而對(duì)區(qū)域上地級(jí)市及內(nèi)部區(qū)縣的定量研究較少。模型數(shù)據(jù)方面，學(xué)者多用引力模型和統(tǒng)計(jì)數(shù)據(jù)分析，但統(tǒng)計(jì)數(shù)據(jù)存在搜索難度大、不齊全、指標(biāo)多的缺點(diǎn)，難以大范圍的分析城鎮(zhèn)化區(qū)位及空間聯(lián)系。而本文選取了區(qū)域地級(jí)市及區(qū)縣，并用引力模型進(jìn)行定量研究，解決了統(tǒng)計(jì)數(shù)據(jù)難以大范圍分析城鎮(zhèn)化區(qū)位及空間聯(lián)系的難題。國(guó)內(nèi)學(xué)者對(duì)城市聯(lián)系的研究多采用數(shù)學(xué)模型和計(jì)量方法[29]，綜觀文獻(xiàn)資料，可以發(fā)現(xiàn)在以往的研究中，多采用引力模型分析聯(lián)系緊密且經(jīng)濟(jì)水平高的東部地區(qū)，而對(duì)于聯(lián)系松散且經(jīng)濟(jì)水平低的西南地區(qū)研究比較缺乏。因此，本文重點(diǎn)選取川南經(jīng)濟(jì)區(qū)作為研究對(duì)象，用2002年、200年7和2013年的DMSP/OLS數(shù)據(jù)表征城市發(fā)展規(guī)模[30-32]，根據(jù)DMSP/OLS影像數(shù)據(jù)分析川南地區(qū)的城鎮(zhèn)發(fā)展總體情況、經(jīng)濟(jì)區(qū)位強(qiáng)度及空間聯(lián)系，在城鎮(zhèn)聯(lián)系上分析內(nèi)部之間的拓展差異性，并輔助以統(tǒng)計(jì)數(shù)據(jù)進(jìn)行區(qū)位強(qiáng)度的驗(yàn)證對(duì)比，在區(qū)縣層次上做進(jìn)一步探討，以期探尋城鎮(zhèn)化發(fā)展的總體趨勢(shì)及城鎮(zhèn)聯(lián)系特點(diǎn)，并為川南經(jīng)濟(jì)區(qū)的建設(shè)發(fā)展提供理論指導(dǎo)。

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

1.1 數(shù)據(jù)來源

DMSP/OLS燈光影像數(shù)據(jù)，由美國(guó)NOAA：(National Oceanic and Atmospheric Administration)下屬的NGDC(National Geophysical Data Center)發(fā)布。已有數(shù)據(jù)主要包括六代DMSP衛(wèi)星F10(1992～1994)、F12(1994～1999)、F14(1997～2003)、F15(2000～2007)、F16(2004～2009)、F18(2010～2013)，涉及1992～2013年共22年的34期影像，該系列影像可在NGDC官方網(wǎng)站中下載(http://ngdc.noaa.gov/eog/dmsp/downloadV4composites.html)。全年的平均燈光影像主要包含無云觀測(cè)的頻數(shù)狀態(tài)、穩(wěn)定狀態(tài)和平均狀態(tài)3種類型，其中的像元DN值(灰度值)表示平均燈光強(qiáng)度，范圍為0～63，該值與燈光強(qiáng)度成正比。穩(wěn)定狀態(tài)的影像排除了背景噪聲的影響，它囊括了城市、縣鎮(zhèn)和其余持久發(fā)光的光源區(qū)域。而且大部分學(xué)者研究用到的是穩(wěn)定燈光產(chǎn)品數(shù)據(jù)，因?yàn)檫@套產(chǎn)品數(shù)據(jù)能有效區(qū)分城市區(qū)和非城鎮(zhèn)區(qū)，且處理了噪聲的影響；其他兩種產(chǎn)品有部分欠缺而相對(duì)使用較少，所以在本文的研究中采用此數(shù)據(jù)。

此外，影像的選取時(shí)段應(yīng)該便于分析研究區(qū)的城鎮(zhèn)差異性，而連續(xù)年份的影像變化不明顯，以往學(xué)者們多采用間隔五年左右的影像，因此，根據(jù)研究時(shí)段選取F142002、F162007、F182013共3期影像。

1.2 研究區(qū)介紹

本次研究范圍以川南經(jīng)濟(jì)區(qū)為研究區(qū)(圖1)，川南經(jīng)濟(jì)區(qū)主要包括內(nèi)江、自貢、瀘州、宜賓四個(gè)地級(jí)市及樂山五縣(井研縣、犍為縣、沐川縣、馬邊彝族自治縣、峨邊彝族自治縣)。對(duì)于樂山五縣，在研究市級(jí)層次時(shí)，為了方便，單獨(dú)處理成一個(gè)獨(dú)立的行政區(qū)域，共計(jì)5個(gè)單元；在研究區(qū)縣級(jí)層次時(shí)，將市轄區(qū)和下轄各縣級(jí)行政區(qū)作為獨(dú)立的研究單元，共計(jì)26個(gè)單元。行政區(qū)邊界方面，中國(guó)省級(jí)、市級(jí)、縣級(jí)行政區(qū)劃圖來自國(guó)家基礎(chǔ)地理信息中心(http://ngcc.sbsm.gov.cn/)，比例尺為1∶400萬，坐標(biāo)系為WGS-84。社會(huì)經(jīng)濟(jì)統(tǒng)計(jì)數(shù)據(jù)方面，本文用到的2002年、2007年、2013年川南經(jīng)濟(jì)區(qū)統(tǒng)計(jì)數(shù)據(jù)來源于四川省歷年統(tǒng)計(jì)年鑒、各地市統(tǒng)計(jì)年鑒，并根據(jù)地方政府統(tǒng)計(jì)公報(bào)進(jìn)行插補(bǔ)完善。

圖1 川南經(jīng)濟(jì)區(qū)地理區(qū)位圖

2 研究方法

2.1 城鎮(zhèn)聯(lián)系理論基礎(chǔ)

城鎮(zhèn)聯(lián)系及發(fā)展的理論研究，多從單一主城區(qū)用地的變化分析城鎮(zhèn)的擴(kuò)展特點(diǎn)，忽視了城市的內(nèi)部差異化；而且研究偏向于國(guó)家及省域?qū)哟蔚某鞘畜w系、空間結(jié)構(gòu)、空間用地拓展等宏觀尺度上的城市經(jīng)濟(jì)區(qū)領(lǐng)域，缺乏區(qū)域?qū)哟紊系丶?jí)市區(qū)縣的城鎮(zhèn)體系、區(qū)位差異和空間聯(lián)系等微觀尺度上的城鎮(zhèn)經(jīng)濟(jì)區(qū)定量研究[33]。本文在選取川南經(jīng)濟(jì)區(qū)研究時(shí)，利用長(zhǎng)時(shí)間尺度的燈光影像數(shù)據(jù)進(jìn)行系列分析，在充分考慮城市內(nèi)部差異性的前提下，從區(qū)域?qū)哟蔚奈⒂^經(jīng)濟(jì)尺度上，對(duì)城鎮(zhèn)經(jīng)濟(jì)區(qū)的城鎮(zhèn)體系、區(qū)位強(qiáng)度差異、城鎮(zhèn)空間聯(lián)系及總體變化趨勢(shì)進(jìn)行定性、定量研究。

城鎮(zhèn)模型數(shù)據(jù)通過引力模型對(duì)影像數(shù)據(jù)處理計(jì)算，克服了統(tǒng)計(jì)數(shù)據(jù)搜索難度大、指標(biāo)多，難以大范圍分析城鎮(zhèn)區(qū)位及空間聯(lián)系的弱點(diǎn)。

理論及實(shí)踐創(chuàng)新上運(yùn)用燈光分布圖與DN值指標(biāo)、燈光強(qiáng)度的空間聯(lián)系相結(jié)合的方法，分析城鎮(zhèn)聯(lián)系是本文的創(chuàng)新之處，區(qū)縣領(lǐng)域的分析，也豐富了區(qū)域性微觀尺度的定量研究。

城鎮(zhèn)聯(lián)系的研究設(shè)計(jì)步驟。

(1)根據(jù)川南行政區(qū)劃裁剪的燈光影像分布圖，對(duì)影像進(jìn)行系統(tǒng)研究，分析川南城鎮(zhèn)的總體分布情況，以此從時(shí)間尺度上對(duì)城鎮(zhèn)體系發(fā)展過程進(jìn)行整體客觀把握。

(2)以引力模型為基礎(chǔ)，根據(jù)傳統(tǒng)統(tǒng)計(jì)數(shù)據(jù)和DMSP/OLS影像數(shù)據(jù)計(jì)算兩兩城市間引力值，由引力矩陣推算川南地區(qū)的區(qū)位強(qiáng)度和差異，同時(shí)以統(tǒng)計(jì)數(shù)據(jù)進(jìn)行區(qū)位強(qiáng)度的驗(yàn)證對(duì)比，以探尋其城鎮(zhèn)化發(fā)展的內(nèi)部差異性。

(3)根據(jù)引力模型計(jì)算川南城際、區(qū)縣單元的引力值，然后在行政區(qū)劃上畫出對(duì)應(yīng)的最大引力分布圖和城鎮(zhèn)聯(lián)系分布圖，以此分析川南城鎮(zhèn)化發(fā)展的總體趨勢(shì)及城鎮(zhèn)聯(lián)系特點(diǎn)，并服務(wù)于川南經(jīng)濟(jì)區(qū)的建設(shè)發(fā)展。

2.2 夜間燈光數(shù)據(jù)處理

DMSP/OLS夜間燈光數(shù)據(jù)，包括很多柵格影像，影像幅寬為3000 km，空間分辨率為1 km，光普分辨率為6 bit。由于不同衛(wèi)星、不同年際的影像數(shù)據(jù)具有差異性和中心飽和性，所以要進(jìn)行數(shù)據(jù)校正[34]。本文應(yīng)用不變目標(biāo)法，選用鶴崗市作為不變目標(biāo)區(qū)，并以無飽和的F162006數(shù)據(jù)為參考對(duì)象，構(gòu)建校正模型進(jìn)行校正?？紤]到川南經(jīng)濟(jì)區(qū)的地理行政位置，并減小地理變形誤差，故采用WGS84的地理坐標(biāo)系，投影方式是蘭勃特等角圓錐投影(Lambert Conformal Conic Projection)。同時(shí)，為了避免中心城市燈光影像的溢出效應(yīng)，這里采取經(jīng)驗(yàn)閥值法進(jìn)行DN總值的提取。

具體操作：首先對(duì)2002年、2007年和2013年的DMSP/OLS燈光數(shù)據(jù)及中國(guó)行政區(qū)劃進(jìn)行蘭勃特等角圓錐投影，轉(zhuǎn)換為直角坐標(biāo)系，然后用川南地區(qū)的行政區(qū)邊界對(duì)3年的DMSP/OLS燈光影像進(jìn)行裁剪，經(jīng)過投影變換和重采樣后，即得到2002年、2007年和2013年的川南經(jīng)濟(jì)區(qū)DMSP/OLS燈光數(shù)據(jù)。

全球的夜間燈光數(shù)據(jù)相互校正是由Elvidge等(2009)提出，本文使用不變目標(biāo)法，具體校正方法如下：

第一步，確定參考區(qū)，根據(jù)年際統(tǒng)計(jì)數(shù)據(jù)變化小且DN值分布范圍廣的雙項(xiàng)原則，即城鎮(zhèn)發(fā)展水平較低，選定鶴崗市[35]。

第二步，確定參考基準(zhǔn)數(shù)據(jù)，根據(jù)DN值累積最高的原則，選取F162006年數(shù)據(jù)。

第三步，構(gòu)建校正模型，以參考影像數(shù)據(jù)值為準(zhǔn)，與待校正的影像數(shù)據(jù)回歸求得相關(guān)參數(shù)(表1)。

DNcorrect=a×DNb

(1)

式中，DNcorrect為校正后的影像灰度值，a和b是常數(shù)。

表1 夜間燈光數(shù)據(jù)冪數(shù)模型的參數(shù)

第四步，根據(jù)校正模型對(duì)影像數(shù)據(jù)進(jìn)行校正。

調(diào)整后的燈光影像如圖2～圖4所示。

圖2 川南經(jīng)濟(jì)區(qū)2002年DMSP/OLS穩(wěn)定燈光強(qiáng)度圖

圖3 川南經(jīng)濟(jì)區(qū)2007年DMSP/OLS穩(wěn)定燈光強(qiáng)度圖

圖4 川南經(jīng)濟(jì)區(qū)2013年DMSP/OLS穩(wěn)定燈光強(qiáng)度圖

2.3 引力模型

現(xiàn)有研究中，多數(shù)學(xué)者選用城際間的功能流向用于分析城市的相互聯(lián)系，一般用人口、經(jīng)濟(jì)、信息、物流等指標(biāo)[36]。關(guān)于城鎮(zhèn)聯(lián)系模型的研究，傳統(tǒng)方法是以統(tǒng)計(jì)數(shù)據(jù)建立模型分析城市間相互聯(lián)系，這種方式整體上比較準(zhǔn)確，并且實(shí)際操作性比實(shí)地調(diào)查相對(duì)簡(jiǎn)便。借助于引力模型，基于夜間燈光影像的定量研究則相對(duì)較少。

引力模型是現(xiàn)有模型中最常用的一種，通常用于研究城市間的作用聯(lián)系。引力模型的含義是城市作用聯(lián)系與經(jīng)濟(jì)規(guī)模成正比，而與城市距離成反比。

(2)

式中，F(xiàn)ij表示聯(lián)系強(qiáng)度，Pi,Pj代表經(jīng)濟(jì)規(guī)模，Dij是城市距離。模型的形式一般不變，不同的經(jīng)濟(jì)問題可用對(duì)應(yīng)的參數(shù)定義。對(duì)于表示區(qū)域城市規(guī)模的變量，考慮到部分?jǐn)?shù)據(jù)難以獲取，在選擇變量和參數(shù)時(shí)，一般選取與經(jīng)濟(jì)活動(dòng)直接相關(guān)的變量，如地區(qū)GDP，人均GDP，人均收入，人口數(shù)量等。

區(qū)位強(qiáng)度的統(tǒng)計(jì)數(shù)據(jù)引力計(jì)算以經(jīng)濟(jì)區(qū)中整體的行政區(qū)域城市為研究對(duì)象，需要注意的是，對(duì)于樂山五縣，在研究市級(jí)層次時(shí)，為了方便計(jì)算，這里單獨(dú)處理成一個(gè)獨(dú)立的行政區(qū)域。引力模型公式為：

(3)

式中，F(xiàn)ij表示城市引力，Gi,Gj代表城市的GDP，Pi,Pj為該地區(qū)的非農(nóng)業(yè)人口，Dij是城市間距離。

區(qū)位強(qiáng)度的燈光數(shù)據(jù)引力計(jì)算參考已有的研究成果，建立城市間空間聯(lián)系強(qiáng)度。引力模型公式為：

(4)

式中，Iij代表該地區(qū)的引力值，Dij代表兩城市中的距離，Pi,Pj代表城市的燈光強(qiáng)度[37]。通過統(tǒng)計(jì)年鑒獲取統(tǒng)計(jì)數(shù)據(jù)，使用模型計(jì)算出川南經(jīng)濟(jì)區(qū)各城市之間的引力值。

區(qū)位強(qiáng)度的綜合計(jì)算通過計(jì)算各城市引力矩陣，即可求得區(qū)位強(qiáng)度。區(qū)位強(qiáng)度的含義是單體城市與領(lǐng)域內(nèi)其它所有城市引力總和占領(lǐng)域內(nèi)全部城市的引力總和百分比。

(5)

F總=∑Fij(i∈[1,n-1]；j∈[i+1,n])

(6)

式中，Ri為城市i的區(qū)位強(qiáng)度，n為區(qū)域內(nèi)的城市數(shù)量。

城鎮(zhèn)空間聯(lián)系的分析主要包括城市最大引力線和引力線分布圖的繪制和經(jīng)濟(jì)擴(kuò)展的趨勢(shì)分析兩部分。本文從城市和區(qū)縣兩個(gè)層次上進(jìn)行分析。

首先畫最大引力連接線，根據(jù)兩兩城市的引力畫出最大引力的連接線然后根據(jù)畫出的引力分布圖，就可以觀察其演變規(guī)律，并歸納總結(jié)空間聯(lián)系特征及聯(lián)系分布圖，分析經(jīng)濟(jì)擴(kuò)展的趨勢(shì)。

整體而言，引力模型在一定的范圍內(nèi)，可使城鎮(zhèn)結(jié)構(gòu)精確化，在城鎮(zhèn)聯(lián)系及發(fā)展的研究中，應(yīng)用廣泛。它能夠精確刻畫城市聯(lián)系的強(qiáng)弱，在管理學(xué)意義上具有重要的區(qū)域經(jīng)濟(jì)價(jià)值。本文根據(jù)DMSP/OLS燈光影像數(shù)據(jù)，運(yùn)用引力模型分析川南經(jīng)濟(jì)區(qū)的城鎮(zhèn)區(qū)位強(qiáng)度、空間聯(lián)系。在空間聯(lián)系分析上，主要通過畫出引力分布圖分析川南經(jīng)濟(jì)區(qū)城鎮(zhèn)聯(lián)系擴(kuò)展規(guī)律和變化特征。

3 結(jié)果與分析

3.1 川南經(jīng)濟(jì)區(qū)燈光影像整體分析

城鎮(zhèn)總體分布：以川南四市為“核心”的多中心網(wǎng)絡(luò)格局基本形成(圖5、圖6)。

圖5 川南經(jīng)濟(jì)區(qū)2002年燈光分布圖

圖6 川南經(jīng)濟(jì)區(qū)2013年燈光分布圖

由圖5和圖6可知，川南經(jīng)濟(jì)區(qū)總體形成以內(nèi)江、自貢、宜賓和瀘州為核心的多中心星狀式結(jié)構(gòu)。川南經(jīng)濟(jì)區(qū)總體分布是一種以自貢、內(nèi)江、瀘州、宜賓為“核心”，以其周邊縣鎮(zhèn)區(qū)域?yàn)椤熬W(wǎng)絡(luò)連接點(diǎn)”的M-C-N結(jié)構(gòu)。區(qū)域內(nèi)的燈光密度呈現(xiàn)出顯著的差異性，按照四川省的行政地區(qū)邊界范圍，其中2013年自貢、內(nèi)江、宜賓和瀘州的燈光密度分別為4.75/km2、6.49/km2、2.64/km2和2.34/km2，是典型的密集區(qū)，而樂山五縣城鎮(zhèn)，其燈光密度僅為0.9/km2，是明顯的低密度區(qū)，川南整體的燈光密度為2.91/km2，從近年來的燈光占比知(表2)，川南四市區(qū)域燈光占比達(dá)94.1%，樂山五縣區(qū)域占比僅為5.9%。盡管瀘州、內(nèi)江兩市占比有下降，但仍然不影響川南四市的主體地位。

表2 川南經(jīng)濟(jì)區(qū)域燈光數(shù)據(jù)

由圖5、圖6可知，目前該區(qū)域呈現(xiàn)顯著的城鎮(zhèn)經(jīng)濟(jì)區(qū)集群特征，以榮縣-自貢市區(qū)-富順縣、內(nèi)江市區(qū)、宜賓縣-宜賓市轄區(qū)-南溪縣和瀘州市區(qū)為核心，以綿遂內(nèi)宜鐵路(銀昆高速)、渝昆鐵路(成渝環(huán)線高速)的長(zhǎng)江段與蓉遵高速公路為主軸發(fā)展線的“點(diǎn)軸格局”特征顯著。點(diǎn)軸外的燈光密集區(qū)僅有隆昌、資中、威遠(yuǎn)和犍為等地零星分布。瀘州市轄區(qū)、宜賓縣-宜賓市轄區(qū)-南溪縣區(qū)域即成渝環(huán)線高速公路及渝昆鐵路川內(nèi)區(qū)域?yàn)殚L(zhǎng)江上游經(jīng)濟(jì)區(qū)的“核心”，該地為連接滇黔渝的主要交通物流中轉(zhuǎn)樞紐，并承擔(dān)川南經(jīng)濟(jì)區(qū)的港口運(yùn)輸和白酒基地的物流經(jīng)濟(jì)功能。

3.2 基于DMSP/OLS和統(tǒng)計(jì)數(shù)據(jù)的川南區(qū)位強(qiáng)度分析

3.2.1基于統(tǒng)計(jì)數(shù)據(jù)的引力模型和城市引力矩陣

川南經(jīng)濟(jì)區(qū)各城市地區(qū)GDP和非農(nóng)業(yè)人口見表3～表5。

表3 川南經(jīng)濟(jì)區(qū)2002年各城市GDP與非農(nóng)業(yè)人口數(shù)據(jù)

注：統(tǒng)計(jì)數(shù)據(jù)來源：四川省統(tǒng)計(jì)年鑒2003，川南經(jīng)濟(jì)區(qū)各城市統(tǒng)計(jì)年鑒2003。

表4 川南經(jīng)濟(jì)區(qū)2007年各城市GDP與非農(nóng)業(yè)人口數(shù)據(jù)

注：統(tǒng)計(jì)數(shù)據(jù)來源：四川省統(tǒng)計(jì)年鑒2008，川南經(jīng)濟(jì)區(qū)各城市統(tǒng)計(jì)年鑒2008。

根據(jù)川南經(jīng)濟(jì)區(qū)區(qū)位圖，用ArcGIS軟件測(cè)算各城市的距離，計(jì)算結(jié)果見表6。

表5 川南經(jīng)濟(jì)區(qū)2013年各城市GDP與非農(nóng)業(yè)人口數(shù)據(jù)

注：統(tǒng)計(jì)數(shù)據(jù)來源：四川省統(tǒng)計(jì)年鑒2014，川南經(jīng)濟(jì)區(qū)各城市統(tǒng)計(jì)年鑒2014。

表6 川南經(jīng)濟(jì)區(qū)各城市間直線距離(單位：公里)

將表3～表5中的GDP與非農(nóng)業(yè)人口統(tǒng)計(jì)數(shù)據(jù)與表6中的各城市距離代入式(3)計(jì)算得到各城市的引力值，最后得到的引力值矩陣(表7～表9)。

表7 川南經(jīng)濟(jì)區(qū)2002年各城市間引力值矩陣(統(tǒng)計(jì)數(shù)據(jù))

表8 川南經(jīng)濟(jì)區(qū)2007年各城市間引力值矩陣(統(tǒng)計(jì)數(shù)據(jù))

表9 川南經(jīng)濟(jì)區(qū)2013年各城市間引力值矩陣(統(tǒng)計(jì)數(shù)據(jù))

3.2.2基于DMSP/OLS燈光數(shù)據(jù)的引力模型和城市間引力矩陣

已有研究證明：DN總值與不同類型產(chǎn)業(yè)數(shù)量之間的相關(guān)性都較高(房地產(chǎn)業(yè)最高)，DN總值與城市發(fā)展水平、區(qū)域地理?xiàng)l件、區(qū)域產(chǎn)業(yè)政策相關(guān)性較大，隨著產(chǎn)業(yè)結(jié)構(gòu)的升級(jí)，DN總值也在增大。DN總值在某種水平上可反映出區(qū)域產(chǎn)業(yè)的結(jié)構(gòu)發(fā)展、演變以及城鎮(zhèn)化水平，它可用于區(qū)域產(chǎn)業(yè)結(jié)構(gòu)及其演變分析、城市區(qū)域發(fā)展和城鎮(zhèn)化規(guī)劃等。

通過對(duì)行政邊界范圍內(nèi)燈光影像的像元DN值進(jìn)行統(tǒng)計(jì)，加權(quán)求和即得到整個(gè)區(qū)域的燈光強(qiáng)度(表10)。

表10 川南經(jīng)濟(jì)區(qū)各城市區(qū)域燈光強(qiáng)度統(tǒng)計(jì)表

將表10中的燈光強(qiáng)度數(shù)據(jù)與表6中的各城市間直線距離代入式(4)計(jì)算各城市之間的引力值，得到的引力值矩陣(表11～表13)。

表11 川南經(jīng)濟(jì)區(qū)2002年各城市間引力值矩陣(燈光數(shù)據(jù))

表12 川南經(jīng)濟(jì)區(qū)2007年各城市間引力值矩陣(燈光數(shù)據(jù))

表13 川南經(jīng)濟(jì)區(qū)2013年各城市間引力值矩陣(燈光數(shù)據(jù))

3.2.3川南區(qū)位強(qiáng)度分析

根據(jù)區(qū)位強(qiáng)度公式(5)式和(6)式，分別對(duì)表7～表9中的統(tǒng)計(jì)引力值和表11～表13的燈光引力值進(jìn)行計(jì)算，計(jì)算結(jié)果見表14、表15。

表14 川南經(jīng)濟(jì)區(qū)各城市間區(qū)位強(qiáng)度統(tǒng)計(jì)表(統(tǒng)計(jì)數(shù)據(jù))

表15 川南經(jīng)濟(jì)區(qū)各城市間區(qū)位強(qiáng)度統(tǒng)計(jì)表(燈光數(shù)據(jù))

通過對(duì)比傳統(tǒng)統(tǒng)計(jì)數(shù)據(jù)和燈光數(shù)據(jù)的區(qū)位強(qiáng)度，可知二者整體上都維持著一致的總體趨勢(shì)，誤差比較小。這表明：夜間燈光數(shù)據(jù)在一定程度上與經(jīng)濟(jì)活動(dòng)方面的參數(shù)有著較高的相關(guān)性，可定量研究城市發(fā)展、城鎮(zhèn)化進(jìn)程、城鎮(zhèn)聯(lián)系、空間結(jié)構(gòu)演變等領(lǐng)域。

自貢和內(nèi)江在三個(gè)年份中利用燈光數(shù)據(jù)和統(tǒng)計(jì)數(shù)據(jù)計(jì)算的值相較其它城市而言，有較小的偏差，變化范圍分別在0.6至0.7，0.5至0.6之間，偏差都在0.1以內(nèi)。其中自貢、宜賓的燈光計(jì)算值始終比統(tǒng)計(jì)計(jì)算值偏小，內(nèi)江和瀘州的燈光計(jì)算值與統(tǒng)計(jì)計(jì)算值時(shí)高時(shí)低，而樂山五縣的燈光計(jì)算值始終比統(tǒng)計(jì)計(jì)算值偏大。近十年各個(gè)城市之間的經(jīng)濟(jì)區(qū)位強(qiáng)度變化不大，相對(duì)都比較穩(wěn)定。按數(shù)值劃分可分為三級(jí)：自貢和內(nèi)江的區(qū)位強(qiáng)度始終處在第一級(jí)，都在0.5以上，瀘州和宜賓處于第二級(jí)，在0.2到0.4之間，第三級(jí)即0至0.2之間的城市為樂山五縣。

綜上所述：川南地區(qū)整體形成自貢、內(nèi)江為核心、宜賓、瀘州為次級(jí)、樂山五縣為三等級(jí)的區(qū)位特征。

川南經(jīng)濟(jì)區(qū)的三個(gè)等級(jí)劃分，主要是受地理上城市距離分布的影響。雖然從GDP、人口和燈光總亮度上自貢、內(nèi)江、瀘州、宜賓相差不大，且瀘州、宜賓的總量大于自貢、內(nèi)江，但是從引力模型上分析，自貢和內(nèi)江的較近距離對(duì)兩個(gè)城市的經(jīng)濟(jì)區(qū)位強(qiáng)度起著決定性作用。而瀘州、宜賓的城市地理分布上較為偏遠(yuǎn)，因此經(jīng)濟(jì)區(qū)位強(qiáng)度相對(duì)較低。

大部分城市在近十年里的影響強(qiáng)度都有略微下降，自貢和內(nèi)江的經(jīng)濟(jì)區(qū)位強(qiáng)度始終保持較高的水平，分別占60%～70%、50%～60%，表明自貢和內(nèi)江處于經(jīng)濟(jì)區(qū)的中心位置，在地理上對(duì)經(jīng)濟(jì)區(qū)中其它城市的總體聯(lián)系影響較大。宜賓在近十年中略有波動(dòng)，其經(jīng)濟(jì)區(qū)位強(qiáng)度占30%～40%，但距離周邊城市相對(duì)較遠(yuǎn)，地理上與自貢和內(nèi)江的聯(lián)系受到一定程度的影響，且南部有較多山區(qū)，在一定程度上對(duì)其它地區(qū)的輻射影響強(qiáng)度有所阻礙。瀘州在近十年中變化比較穩(wěn)定，其經(jīng)濟(jì)區(qū)位強(qiáng)度占25%～35%，距離周邊城市都很遠(yuǎn)，地理上與其他城市的聯(lián)系受到很大程度的影響，在一定程度上對(duì)其它地區(qū)的輻射影響強(qiáng)度阻礙較大。樂山五縣是從樂山分割出的幾個(gè)縣城，與川南四市相比，區(qū)域面積、經(jīng)濟(jì)總量和人口都很小，而且地理分布呈現(xiàn)由南向北發(fā)展，其經(jīng)濟(jì)區(qū)位強(qiáng)度占7%～16%，雖然距離自貢、內(nèi)江、宜賓相對(duì)較近，但其經(jīng)濟(jì)總量和人口都很小，因此在一定程度上對(duì)其它地區(qū)的輻射影響強(qiáng)度最弱。

3.3 川南經(jīng)濟(jì)區(qū)城鎮(zhèn)空間聯(lián)系強(qiáng)度分析

3.3.1川南經(jīng)濟(jì)區(qū)城市空間聯(lián)系強(qiáng)度分析

圖7 2002年川南經(jīng)濟(jì)區(qū)最大引力線分布圖

圖8 2007年川南經(jīng)濟(jì)區(qū)最大引力線分布圖

圖9 2013年川南經(jīng)濟(jì)區(qū)最大引力線分布圖

圖10 川南經(jīng)濟(jì)區(qū)城市聯(lián)系分布圖

將城市聯(lián)系的引力分布圖按聯(lián)系強(qiáng)弱劃分等級(jí)，引力值的總和越大，聯(lián)系越強(qiáng)。城市的引力線越粗，數(shù)目越多，該城市的引力和聯(lián)系就越強(qiáng)。由圖10可知，川南地區(qū)城市一直保持“雙核雙副”的多中心城市團(tuán)體系，城際聯(lián)系已形成特殊的點(diǎn)軸狀結(jié)構(gòu)，在各時(shí)間節(jié)點(diǎn)均顯著；內(nèi)江、自貢、宜賓、瀘州始終是川南經(jīng)濟(jì)區(qū)的多核心群團(tuán)。隨著城市規(guī)模的拓展，城市的結(jié)構(gòu)聯(lián)系軸帶化趨勢(shì)十分明顯。

2002年，“雙核雙副”城市團(tuán)結(jié)構(gòu)已經(jīng)突顯，其中自貢和內(nèi)江引力值最強(qiáng)，城市聯(lián)系最為緊密，其次是宜賓和自貢，瀘州和宜賓、樂山五縣和宜賓聯(lián)系較小。2007年，除原有的自貢和內(nèi)江引力最強(qiáng)，聯(lián)系最為緊密外，宜賓和自貢的聯(lián)系逐漸加強(qiáng)并升級(jí)為一級(jí)聯(lián)系，瀘州和宜賓之間的聯(lián)系、樂山五縣和宜賓之間的聯(lián)系都不斷增強(qiáng)和升級(jí)，軸帶化趨勢(shì)開始顯現(xiàn)。2013年，自貢和內(nèi)江兩市的影響范圍繼續(xù)擴(kuò)大，各市的輻射能力均有增強(qiáng)，其中瀘州和宜賓聯(lián)系已升級(jí)為一級(jí)，輻射能力明顯增強(qiáng)，樂山五縣的總體影響范圍也略有顯現(xiàn)。

從燈光斑塊分布分析，2002年，川南四市只是點(diǎn)狀發(fā)展，城市聯(lián)系較少，川南各縣以四市為核心發(fā)展；2007年，川南四市點(diǎn)狀發(fā)展不斷增強(qiáng)，其中宜賓市轄區(qū)開始向西南宜賓縣與屏山縣交界處發(fā)展，瀘州市轄區(qū)也沿著長(zhǎng)江一帶向西南延伸，瀘州、宜賓的長(zhǎng)江沿岸軸帶化趨勢(shì)加深；2013年，四市斑塊有顯著的的擴(kuò)張發(fā)展，其中自貢市轄區(qū)向東南延伸，和富順連接成片；內(nèi)江市轄區(qū)開始向自貢方向延伸；宜賓向西南宜賓縣和屏山縣交界處延伸，宜賓市轄區(qū)東部也有所擴(kuò)張，軸帶化趨勢(shì)很強(qiáng)；瀘州市轄區(qū)分別向西南和東部延伸，這和宜賓的趨勢(shì)相同；在自然環(huán)境、交通通道和區(qū)位條件上，自貢市區(qū)向東南發(fā)展，主要是得益于成瀘高速和省道的交通優(yōu)勢(shì)；內(nèi)江向自貢方向延伸則突顯兩市聯(lián)系緊密，經(jīng)濟(jì)一體化的利益需求；宜賓和瀘州的軸帶化趨勢(shì)主要是由于兩市依附于長(zhǎng)江一帶臨港產(chǎn)業(yè)，并且有成渝環(huán)線高速和渝昆鐵路的物流區(qū)位優(yōu)勢(shì)所致。

3.3.2川南經(jīng)濟(jì)區(qū)各區(qū)縣空間聯(lián)系強(qiáng)度分析

各區(qū)縣空間聯(lián)系分析主要包括區(qū)縣最大引力線和引力線分布圖的繪制和經(jīng)濟(jì)擴(kuò)展的趨勢(shì)分析兩部分。類比城市空間聯(lián)系的方法，區(qū)縣的引力分布圖如圖11～圖13所示。根據(jù)引力圖即可得到區(qū)縣層級(jí)的城鎮(zhèn)聯(lián)系分布圖(圖14)。

圖11 2002年川南各區(qū)縣最大引力線分布圖

圖12 2007年川南各區(qū)縣最大引力線分布圖

圖13 2013年川南各區(qū)縣最大引力線分布圖

圖14 川南區(qū)縣城鎮(zhèn)聯(lián)系分布圖

川南經(jīng)濟(jì)區(qū)城鎮(zhèn)內(nèi)部空間形成了以自貢、內(nèi)江市區(qū)為核心，瀘州、宜賓市區(qū)為副核心的“雙核雙副”的網(wǎng)絡(luò)化圈層結(jié)構(gòu)體系，川南區(qū)縣內(nèi)部分化為以自貢-內(nèi)江市區(qū)、瀘州市區(qū)和宜賓市區(qū)為中心的3個(gè)城鎮(zhèn)聯(lián)系經(jīng)濟(jì)圈(圖15)。在城鎮(zhèn)帶聯(lián)系發(fā)展中形成“一軸四帶”的結(jié)構(gòu)體系，在各時(shí)間節(jié)點(diǎn)均顯著。內(nèi)江、自貢、宜賓與瀘州市區(qū)始終是川南經(jīng)濟(jì)區(qū)內(nèi)的多核心團(tuán)體。隨著城鎮(zhèn)規(guī)模的擴(kuò)張，區(qū)縣結(jié)構(gòu)網(wǎng)絡(luò)化趨勢(shì)明顯。2002年自貢市區(qū)和內(nèi)江市區(qū)的網(wǎng)絡(luò)化輻射已比較明顯，瀘州市區(qū)-瀘縣-隆昌縣開始軸帶化，宜賓市區(qū)開始星狀輻射發(fā)展，其中自貢市區(qū)-威遠(yuǎn)縣、宜賓市區(qū)-宜賓縣的引力值最強(qiáng)，區(qū)縣聯(lián)系最為緊密；其次是自貢市區(qū)-內(nèi)江市區(qū)、自貢市區(qū)-富順縣、內(nèi)江市區(qū)-資中縣、內(nèi)江市區(qū)-威遠(yuǎn)縣、內(nèi)江市區(qū)-隆昌縣、榮縣-威遠(yuǎn)縣，瀘州市區(qū)-瀘縣、瀘縣-隆昌縣；樂山五縣與其它地區(qū)聯(lián)系較少，內(nèi)部聯(lián)系則以犍為縣為中心。自貢、宜賓、瀘州市區(qū)間聯(lián)系較弱。2007年，在2002年原有的區(qū)縣聯(lián)系體系基礎(chǔ)上，宜賓市區(qū)和宜賓縣聯(lián)系逐漸加強(qiáng)并升級(jí)為一級(jí)聯(lián)系；此外，瀘州市區(qū)-合江縣、高縣-珙縣、宜賓市區(qū)-長(zhǎng)寧縣/屏山縣的聯(lián)系也逐步的加強(qiáng)；川南區(qū)縣網(wǎng)絡(luò)化趨勢(shì)開始顯現(xiàn)出來。但是，自貢、宜賓、瀘州市區(qū)間聯(lián)系依然較弱，樂山五縣與其它地區(qū)聯(lián)系依然較少，內(nèi)部聯(lián)系沒有變化。2013年與2007年相比，川南四市城區(qū)的影響范圍繼續(xù)擴(kuò)大，各市的輻射能力均有增強(qiáng)，自貢市區(qū)和威遠(yuǎn)縣升級(jí)為一級(jí)聯(lián)系，自貢市區(qū)-富順縣、自貢市區(qū)-內(nèi)江市區(qū)、內(nèi)江市區(qū)-資中縣、內(nèi)江市區(qū)-威遠(yuǎn)縣、內(nèi)江市區(qū)-隆昌縣、榮縣-威遠(yuǎn)縣、瀘州市區(qū)-瀘縣及合江縣、瀘縣-隆昌縣、長(zhǎng)寧縣-江安縣、興文縣-敘永縣、犍為縣-沐川縣聯(lián)系均得到增強(qiáng)；自貢、宜賓、瀘州市區(qū)間聯(lián)系增強(qiáng)，川南區(qū)縣網(wǎng)絡(luò)化趨勢(shì)已十分明顯。不過，樂山五縣的總體對(duì)外聯(lián)系仍沒有顯現(xiàn)出來。

圖15 川南區(qū)縣城鎮(zhèn)帶分布圖

4 結(jié)束語

本文基于DMSP/OLS夜間燈光數(shù)據(jù)，以傳統(tǒng)統(tǒng)計(jì)數(shù)據(jù)為輔助，運(yùn)用引力模型對(duì)川南經(jīng)濟(jì)區(qū)的城鎮(zhèn)區(qū)位強(qiáng)度、空間聯(lián)系進(jìn)行了研究，并探討了城鎮(zhèn)化發(fā)展趨勢(shì)及城鎮(zhèn)聯(lián)系特點(diǎn)。通過計(jì)算分析發(fā)現(xiàn)，川南經(jīng)濟(jì)區(qū)總體形成以自貢、內(nèi)江為核心，瀘州、宜賓為副核心的“雙核雙副”多中心網(wǎng)絡(luò)結(jié)構(gòu)。受地理影響，自貢、內(nèi)江的經(jīng)濟(jì)區(qū)位強(qiáng)度最強(qiáng)，川南內(nèi)部分化為三個(gè)區(qū)位等級(jí)；空間聯(lián)系上，川南經(jīng)濟(jì)區(qū)空間內(nèi)部長(zhǎng)期保持著“雙核雙副”的城市團(tuán)體系，城際聯(lián)系形成特殊的點(diǎn)軸狀結(jié)構(gòu)，城際軸帶化趨勢(shì)明顯，在各時(shí)間節(jié)點(diǎn)均顯著；區(qū)縣聯(lián)系層次上，川南經(jīng)濟(jì)區(qū)內(nèi)形成了“雙核雙副一軸四帶”的區(qū)縣團(tuán)體系，隨著城鎮(zhèn)規(guī)模的擴(kuò)張，城鎮(zhèn)聯(lián)系不斷的深化，各縣鎮(zhèn)圍繞中心城區(qū)發(fā)展，區(qū)縣間結(jié)構(gòu)網(wǎng)絡(luò)化趨勢(shì)明顯。區(qū)縣領(lǐng)域的燈光數(shù)據(jù)研究為城鎮(zhèn)空間結(jié)構(gòu)及經(jīng)濟(jì)聯(lián)系提供了一種更為細(xì)致的研究方法，豐富了城鎮(zhèn)聯(lián)系的定量研究。

此外，在計(jì)算區(qū)縣層次的城鎮(zhèn)聯(lián)系時(shí)還發(fā)現(xiàn)，根據(jù)城鎮(zhèn)結(jié)構(gòu)、經(jīng)濟(jì)發(fā)展情況，自貢、內(nèi)江、宜賓、瀘州四個(gè)城市周邊已形成三個(gè)城鎮(zhèn)聯(lián)系經(jīng)濟(jì)圈，周邊區(qū)縣與中心市區(qū)聯(lián)系十分緊密，亟需以區(qū)域的發(fā)展理念指導(dǎo)城市規(guī)劃，以期推動(dòng)區(qū)域城市化的發(fā)展。目前，川南城鎮(zhèn)化正在城市范疇內(nèi)擴(kuò)散，政府在規(guī)劃時(shí)可著手從政策上給予自貢和內(nèi)江的合作環(huán)境，并促進(jìn)宜賓、瀘州兩地的聯(lián)系，加快構(gòu)建自貢與內(nèi)江、瀘州、宜賓等區(qū)域性城市團(tuán)，促進(jìn)城鎮(zhèn)一體化的發(fā)展。

運(yùn)用引力模型分析城市空間聯(lián)系時(shí)，以夜間燈光亮度替代傳統(tǒng)統(tǒng)計(jì)數(shù)據(jù)的研究方法還有很多細(xì)節(jié)值得商榷，例如提取燈光數(shù)據(jù)的閥值選取，以歐式空間距離表述城市通達(dá)性的實(shí)效程度，不同地區(qū)利用引力模型所達(dá)到的實(shí)際效果。還有，目前學(xué)者們多數(shù)側(cè)重于省級(jí)大尺度城市聯(lián)系的定量及定性分析，而缺乏小尺度區(qū)縣層次的定量分析，這方面以后可作為研究者的一個(gè)研究方向。

[1] 方創(chuàng)琳,周成虎,顧朝林,等.特大城市群地區(qū)城鎮(zhèn)化與生態(tài)環(huán)境交互耦合效應(yīng)解析的理論框架及技術(shù)路徑[J].地理學(xué)報(bào),2016,71(4):531-550.

[2] 蔡凌雁.城市化后期城鎮(zhèn)建設(shè)用地空間演變模擬研究[D].南京:南京大學(xué),2017.

[3] 方創(chuàng)琳.改革開放30年來中國(guó)的城市化與城鎮(zhèn)發(fā)展[J].經(jīng)濟(jì)地理,2009,29(1):19-25.

[4] 李娜.基于夜間燈光數(shù)據(jù)對(duì)我國(guó)城鎮(zhèn)發(fā)展的探討[D].蘭州:蘭州大學(xué),2016.

[5] 陳玉光.大城市空間擴(kuò)展方式與空間結(jié)構(gòu)研究的歷史演進(jìn)及未來展望[J].江淮論壇,2017,285(5):60-65.

[6] 李保杰,頤和和,紀(jì)亞洲,等.基于地學(xué)信息圖譜的礦業(yè)城市空間擴(kuò)展研究——以徐州市為例[J].地域研究與開發(fā),2012,31(1):50-54.

[7] 陳東湘,周生路,吳紹華.基于遙感評(píng)價(jià)城市擴(kuò)張對(duì)耕地質(zhì)量等級(jí)結(jié)構(gòu)及產(chǎn)能的影響[J].農(nóng)業(yè)工程學(xué)報(bào),2017,33(13):264-269.

[8] 李全林,馬曉冬,朱傳耿,等.基于GIS的鹽城城市空間結(jié)構(gòu)演化分析[J].地理與地理信息科學(xué),2007,23(3):69-73.

[9] 董晨煒,曹宇,譚永忠.基于夜間燈光數(shù)據(jù)的環(huán)杭州灣城市擴(kuò)張及植被變化[J].應(yīng)用生態(tài)學(xué)報(bào),2017,28(1):231-238.

[10] 王曉慧.基于DMSP/OLS夜間燈光數(shù)據(jù)的中國(guó)近30年城鎮(zhèn)擴(kuò)展研究[D].南京:南京大學(xué),2013.

[11] YI N H,JIAN P,YANXU L,et al.Mapping Development Pattern in Beijing-Tianjin-Hebei Urban Agglomeration Using DMSP/OLS Nighttime Light Data[J].Remote Sensing,2017,9(7):760.

[12] RICHARD F,CHARLOTTA M,TIM G.Global Metropolis:Assessing Economic Activity in Urban Centers Based on Night-time Satellite Images[J].Professional Geographer,2010,64(2):178-187.

[13] LIU H,LIN M A,GUOPING L I,et al.Pattern evolution and its contributory factor of cold spots and hot spots of economic development in Beijing-Tianjin-Hebei region[J].Geographical Research,2017(1):97-108.

[14] QI K,HU Y,CHENG C,et al.Transferability of Economy Estimation Based on DMSP/OLS Night-Time Light[J].Remote Sensing,2017,9(8):786.

[15] HENDERSON J V,STOREYGARD A,WEIL D N.Measuring Economic Growth from Outer Space[J].American Economic Review,2012,102(2):994-1028.

[16] TRIPATHY B R,TIWARI V,PANDEY V,et al.Estimation of Urban Population Dynamics Using DMSP/OLS Night-Time Lights Time Series Sensors Data[J].IEEE Sensors Journal,2017,17(4):1013-1020.

[17] HUANG Q,YANG Y,LI Y,et al.A Simulation Study on the Urban Population of China Based on Nighttime Light Data Acquired from DMSP/OLS[J].Sustainability,2016,8(6):521.

[18] TRIPATHY B R,TIWARI V,PANDEY V,et al.Estimation of Urban Population Dynamics Using DMSP/OLS Night-Time Lights Time Series Sensors Data[J].IEEE Sensors Journal,2017,17(4):1013-1020.

[19] SUTTON P C,TAYLOR M J,ELVIDGE C D.Using DMSP/OLS Imagery to Characterize Urban Populations in Developed and Developing Countries[M].Remote Sensing of Urban and Suburban Areas,2010,329-348.

[20] XIE Y,WENG Q.World energy consumption pattern as revealed by DMSP-OLS nighttime light imagery[J].Mapping Sciences & Remote Sensing,2016,53(2):265-282.

[21] ZHANG X,WU J,PENG J,et al.The Uncertainty of Nighttime Light Data in Estimating Carbon Dioxide Emissions in China:A Comparison between DMSP/OLS and NPP-VIIRS[J].Remote Sensing,2017,9(8):797.

[22] 鄒進(jìn)貴,陳艷華,田徑,等.基于ArcGIS的DMSP/OLS夜間燈光影像校正模型的構(gòu)建[J].測(cè)繪地理信息,2014,39(4):33-37.

[23] LI X,ZHOU Y.A Stepwise Calibration of Global DMSP/OLS Stable Nighttime Light Data (1992-2013)[J].Remote Sensing,2017,9(6):63701-63714.

[24] JIA T,CHEN K,WANG J.Characterizing the Growth Patterns of 45 Major Metropolitans in Mainland China Using DMSP/OLS Data[J].Remote Sensing,2017,9(6):571.

[25] 楊眉,王世新,周藝,等.DMSP/OLS夜間燈光數(shù)據(jù)應(yīng)用研究綜述[J].遙感技術(shù)與應(yīng)用,2011,26(1):45-51.

[26] SUTTON P,ROBERTS D,ELVIDGE C,et al.Census from Heaven:An estimate of the global human population using night-time satellite imagery[J].International Journal of Remote Sensing,2001,22(16):3061-3076.

[27] LIU Y,YANG Y,JING W,et al.A new urban index for expressing inner-city patterns based on MODIS LST and EVI regulated DMSP/OLS NTL[J].Remote Sensing,2017,9(8):777.

[28] BENNETT M M,SMITH L C.Advances in using multitemporal night-time lights satellite imagery to detect,estimate,and monitor socioeconomic dynamics[J].Remote Sensing of Environment,2017,192:176-197.

[29] 張仲元.基于GIS的中原城市群空間結(jié)構(gòu)與布局研究[D].開封:河南大學(xué),2008.

[30] 顧朝林,龐海峰.基于重力模型的中國(guó)城市體系空間聯(lián)系與層域劃分[J].地理研究,2008,27(1):1-12.

[31] 張陸,高素英.多中心視角下的京津冀都市圈空間聯(lián)系分析[J].城市發(fā)展研究,2014,21(5):49-54.

[32] 劉建朝,高素英.基于城市聯(lián)系強(qiáng)度與城市流的京津冀城市群空間聯(lián)系研究[J].地域研究與開發(fā),2013,32(2):57-61.

[33] 吳健生,赫勝彬,彭建,等.基于DMSP-OLS數(shù)據(jù)的城市發(fā)展空間特征研究[J].地理與地理信息科學(xué),2014,30(2):20-25.

[34] 曹子陽,吳志峰,匡耀求,等.DMSP/OLS夜間燈光影像中國(guó)區(qū)域的校正及應(yīng)用[J].地球信息科學(xué)學(xué)報(bào),2015,17(9):1092-1102.

[35] WU J,HE S,PENG J,et al.Intercalibration of DMSP/OLS night-time light data by the invariant region method[J].International Journal of Remote Sensing,2013,34(20):7356-7368.

[36] 徐慧超,韓增林,趙林,等.中原經(jīng)濟(jì)區(qū)城市經(jīng)濟(jì)聯(lián)系時(shí)空變化分析——基于城市流強(qiáng)度的視角[J].經(jīng)濟(jì)地理,2013,33(6):53-58.

[37] 張雨晨.基于DMSP/OLS夜間燈光數(shù)據(jù)的中原城市群時(shí)空格局分析[D].成都:成都理工大學(xué),2016.