俞宵,李藝,范俊甫*,韓留生,周玉科

基于穩(wěn)定夜間燈光數(shù)據(jù)的山東省城市腹地提取與變化分析

俞宵1,李藝1,范俊甫1*,韓留生1,周玉科2*

1. 山東理工大學建筑工程學院, 山東 淄博 255049 2. 中國科學院地理科學與資源研究所生態(tài)系統(tǒng)網(wǎng)絡(luò)觀測與模擬院重點實驗室, 北京 100101

城市腹地提取是探索城市服務范圍，衡量城市間聯(lián)系程度的重要方法。針對城市腹地提取，本文提出了一種新的城市腹地提取算法，該算法采用迭代閾值分割法提取城市建成區(qū)范圍，使用建成區(qū)面積作為城市質(zhì)量計算指標，并將加權(quán)Voronoi圖與引力模型方法結(jié)合；將該算法用于2010年至2013年山東省區(qū)域的4景穩(wěn)定DMSP/OLS影像數(shù)據(jù)，提取了山東省17個地級市像元級尺度的腹地空間分布；并將本文所得結(jié)果與基于交通可達性和城市競爭力腹地范圍提取結(jié)果作對比，驗證提取結(jié)果精度；進一步分析了4年間各市的腹地范圍變化，并提出腹地均衡系數(shù)，探索城市發(fā)展方向。研究表明：（1）本文方法與傳統(tǒng)經(jīng)典方法在結(jié)果上重疊度達到73%以上，且解決了存在細小斑塊空洞腹地以及腹地延道路不合理延伸問題。（2）萊蕪與泰安被濟南的城市腹地所包圍，腹地范圍打破行政邊界壁壘。（3）2010至2013年山東各市腹地均衡指數(shù)不斷增加，具有城市資源由集約化向均衡化發(fā)展的明顯地理規(guī)律。（5）山東半島藍色經(jīng)濟區(qū)各城市均衡指數(shù)明顯高于山東其他城市，2010至2013年山東半島藍色經(jīng)濟區(qū)腹地均衡指數(shù)由0.672上升為0.882，其他城市腹地均衡指數(shù)由0.145上升至0.528。

燈光數(shù)據(jù);城市建設(shè); 空間分布

城市發(fā)展離不開周邊區(qū)域的支撐，準確界定城市范圍及其腹地,有助于明確城市發(fā)展的生態(tài)環(huán)境及資源基底,對于城市可持續(xù)發(fā)展決策意義重大。因此,城市腹地的識別一直以來都是城市地理學、區(qū)域經(jīng)濟學的研究重點?！案沟亍逼鹪从趯Ω劭诔擎?zhèn)的研究，即港口集散物資的地域，現(xiàn)在被引申為位于港口或者內(nèi)陸中心城鎮(zhèn)周邊的并與其存在密切社會經(jīng)濟聯(lián)系的毗鄰地區(qū)[1]。作為空間集聚主體的城鎮(zhèn)，其發(fā)展需要周邊區(qū)域提供資源，精確識別城鎮(zhèn)與其周邊區(qū)域即腹地，有利于提升城鎮(zhèn)的功能與核心競爭力，對城鎮(zhèn)可持續(xù)發(fā)展以及行政區(qū)的劃分有著指導性意義。

腹地的研究最早可追溯至1826年，杜能首次闡述了農(nóng)業(yè)區(qū)位論的思想，以圈層方式劃分空間農(nóng)業(yè)結(jié)構(gòu)[2]。而后韋伯提出工業(yè)區(qū)位論，表明區(qū)位因子決定生產(chǎn)場所，將企業(yè)吸引到生產(chǎn)費用最小、節(jié)約費用最大的地點，為區(qū)域規(guī)劃與城市規(guī)劃提供了理論基礎(chǔ)[3]。1933年克里斯塔勒，他首次指定腹地邊界，城鎮(zhèn)的服務上限即消費者愿意去一個中心地得到貨物或者服務的最遠距離，服務下限即門檻距離是維持一項中心地職能經(jīng)營所必需的腹地的最短距離[4]。20世紀50年代之后腹地研究由定性描述過渡到定量的數(shù)學模擬法。1955年Green通過交通通勤駛向、電話通信流向、報紙銷售區(qū)域、公司及銀行決策者工作點5項指標探究了紐約與波士頓兩個大城鎮(zhèn)的城鎮(zhèn)相互作用，并計算了兩者之間影響范圍并劃定影響邊界[5]。20世紀中期以后，隨著經(jīng)濟社會急遽發(fā)展，腹地有關(guān)的理論爆發(fā)式呈現(xiàn)：增長中心理論、核心-邊緣理論、都市-腹地理論、經(jīng)濟輻射理論等從城鎮(zhèn)等級結(jié)構(gòu)、集聚擴散、不平等發(fā)展等多角度探究城鎮(zhèn)與腹地的聯(lián)系，并逐步形成腹地界定方法[6-8]。

我國開展腹地研究始于20世紀80年代，腹地界定方法主要分為：（1）基于場強模型與引力模型測度方法。（2）基于交通網(wǎng)絡(luò)最短時間距離測度方法。（3）基于人流、物流、技術(shù)流、信息流、資本流等“流”的測度方法。其中場強模型與引力模型方法采用行政區(qū)經(jīng)濟指標作為城市質(zhì)量，行政區(qū)間的通行時間作為距離計算不同區(qū)域兩兩之間的引力值，依據(jù)最大值原則判斷腹地歸屬[4,9-11]。交通網(wǎng)絡(luò)時間距離是衡量兩個地區(qū)聯(lián)系強度的重要指標，基于交通網(wǎng)絡(luò)的腹地劃分方法通過計算不同城市間的通行距離，按照依據(jù)就近原則劃分腹地范圍[12]。信息是“網(wǎng)絡(luò)城市”背景下城市聯(lián)系的重要度量指標，“信息流”不僅在表征“可達性”、對場強模型進行距離修正上具有重要的作用，而且“信息腹地”也日益成為腹地識別的重要方面[1]。信息流方法通過信息的通信強度作為衡量兩個城市的聯(lián)系強度指標，依據(jù)聯(lián)系強度判斷腹地歸屬?？傮w來看，多源數(shù)據(jù)對于信息腹地的細致考量是對傳統(tǒng)城市腹地定量識別的有效修正與補充。以上方法以點代面，忽視了城市的面域性，城市發(fā)展的方向性；還忽略了城市發(fā)達區(qū)域可能存在多個發(fā)達中心等特征。

遙感技術(shù)憑借著數(shù)據(jù)的客觀性，多源性，現(xiàn)勢性，動態(tài)性等特點，在自然資源普查，社會經(jīng)濟發(fā)展，城市擴張等領(lǐng)域應用廣泛。夜間燈光數(shù)據(jù)作為一種間接衡量社會經(jīng)濟發(fā)展的特殊遙感數(shù)據(jù)，被廣泛使用在社會經(jīng)濟發(fā)展以及戰(zhàn)爭破壞程度評估等領(lǐng)域[13-15]。

夜間燈光數(shù)據(jù)作為一種融合了經(jīng)濟與位置且更加準確以及客觀的遙感數(shù)據(jù)，很適合用于解決腹地劃分問題，然而選取何種算法將同時包含屬性與空間信息的夜間燈光數(shù)據(jù)應用在腹地提取中，目前沒有相應的研究作為參考。

基于以上問題，本文提出一種新的基于遙感的腹地提取算法，該算法基于穩(wěn)定夜間燈光數(shù)據(jù)，并將引力模型與面加權(quán)Voronoi圖結(jié)合，劃分像元尺度的腹地范圍，并提出腹地均衡指數(shù)衡量腹地面積與行政區(qū)面積的偏離程度。以期為山東省空間結(jié)構(gòu)優(yōu)化整合提供相應的理論依據(jù)和實踐參考。

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

1.1 研究區(qū)域

山東省研究區(qū)域概況：山東省，省會濟南，位于中國東部沿海，界于北緯34°22.9′-38°24.01′，東經(jīng)114°47.5′-122°42.3′之間，山東省總面積15.71萬km2，山東省是中國經(jīng)濟最發(fā)達的省份之一，也是中國的農(nóng)業(yè)與工業(yè)大省。2011年1月4日，國務院批復《山東半島藍色經(jīng)濟區(qū)發(fā)展規(guī)劃》，標志著全國海洋經(jīng)濟發(fā)展試點工作進入實施階段，也標志著山東半島藍色經(jīng)濟區(qū)建設(shè)正式上升為國家戰(zhàn)略，成為國家層面海洋發(fā)展戰(zhàn)略和區(qū)域協(xié)調(diào)發(fā)展戰(zhàn)略的重要組成部分。山東半島藍色經(jīng)濟區(qū)，是中國第一個以海洋經(jīng)濟為主題的區(qū)域發(fā)展戰(zhàn)略，是中國區(qū)域發(fā)展從陸域經(jīng)濟延伸到海洋經(jīng)濟、積極推進陸海統(tǒng)籌的重大戰(zhàn)略舉措。規(guī)劃主體區(qū)范圍包括山東全部海域和青島、煙臺、威海，濰坊，淄博、東營，日照等八市及濱州的無棣、沾化2個沿海縣所屬陸域，海域面積15.95萬km2，陸域面積6.4萬km2 [16]。本研究以山東省17個地級市為研究對象，研究范圍如圖1：

圖1 山東省市級行政邊界

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

依據(jù)業(yè)務型線掃描傳感器，美國國防氣象衛(wèi)星計劃(DMSP)提供了獲取城市燈光信息的新數(shù)據(jù)源。該衛(wèi)星傳感器能靈敏地探測到包含城市燈光，汽車，船只等發(fā)出的燈光，并準確有效地區(qū)分城市燈光和黑暗農(nóng)村，本文使用2010-2013年共四年的DMSP/OLS的非輻射定標穩(wěn)定夜間燈光數(shù)據(jù)，數(shù)據(jù)像元值范圍為1-63，分辨率為1 km。2010-2013年四年山東省各市建成區(qū)面積以及城市競爭力指標來自《山東省統(tǒng)計年鑒》以及各市國民經(jīng)濟與社會發(fā)展統(tǒng)計公報。用于腹地劃分的交通網(wǎng)絡(luò)數(shù)據(jù)來源于全國地理信息資源目錄服務系統(tǒng)提供的全國1:1 000 000基礎(chǔ)地理數(shù)據(jù)中的道路矢量線數(shù)據(jù)，包括；標準軌鐵路、窄軌鐵路等鐵路，以及國道、省道、縣道、鄉(xiāng)道、其它公路、街道、鄉(xiāng)村道路等公路。由于高速公路和鐵路具有一定的封閉性, 只有互通口或站點位置才能與周邊地區(qū)發(fā)生聯(lián)系，將高速以及鐵路除互通口以及站點進行封閉處理。

2 研究方法

2.1 基于夜間燈光的腹地劃分算法流程

腹地劃分基本流程如圖所示：

圖2 腹地劃分基本流程

依據(jù)夜間燈光數(shù)據(jù)以及統(tǒng)計年鑒中四年間的建成區(qū)面積采用迭代閾值分割法獲取建成區(qū)，結(jié)合建成區(qū)數(shù)據(jù)采用公式（1）計算每個像元與17個市建成區(qū)的引力值，依據(jù)面加權(quán)Voronoi圖原則，得到17個市的腹地范圍。將基于夜間燈光的腹地劃分結(jié)果與基于交通網(wǎng)絡(luò)數(shù)據(jù)的腹地劃分結(jié)果進行相對合理性評價。使用公式（5）將行政區(qū)劃數(shù)據(jù)與腹地數(shù)據(jù)進行均衡指數(shù)計算。

2.1.1 建成區(qū)范圍提取為提取腹地中心的范圍，考慮到各個城市年度城區(qū)范圍的不確定性，以及夜間燈光的飽和特征，參照文獻[13]，以山東省統(tǒng)計年鑒的各市建成區(qū)的面積為基礎(chǔ)，采用燈光DN遞增方法提取山東省17個市級行政區(qū)的建成區(qū)作為城市腹地中心，方法流程如圖3所示：

圖3 腹地中心范圍提取算法流程圖

其中：()為燈光值，初始值為1，()為像素值大于()的柵格單元總面積，為各市的建成區(qū)面積。

2.1.2 引力模型 引力模型來源于牛頓1687年提出的萬有引力定律，用于測度城鎮(zhèn)空間聯(lián)系強弱。本研究通過改進引力模型，測度不同城市間的引力強度時引入夜間燈光強度作為權(quán)重系數(shù)，公式如下：

其中：Q為城市質(zhì)量，計算公式為：

其中，F為柵格單元與城市之間的引力；其中Q表示需要被劃分的城市柵格單元質(zhì)量，Q為中心城市的城市質(zhì)量，即中心城市夜間燈光閾值分割后的面積；d表示地到地的距離，本研究采用歐式距離；為系數(shù)，本文參考文獻取1；為距離摩擦系數(shù)，本文取2，D為城市的夜間燈光值總和，為整個研究區(qū)域的燈光值總和；S為城市腹地中心面積。

2.1.3 面加權(quán)Voronoi圖加權(quán)Voronoi圖能劃分最合理服務范圍，可以為腹地劃分提供科學的理論依據(jù)?？紤]到各城市建成區(qū)為面數(shù)據(jù)，因此本文采用面加權(quán)Voronoi圖作為腹地劃分算法。面加權(quán)Voronoi圖常用于劃分服務區(qū)范圍，定義如下：

其中Vor表示當前中心城區(qū)面域A的Voronoi區(qū)域，為歐幾里得平面上的任意一點，d為任意點與中心城區(qū)面域之間的最短歐式距離，W，W分別為A與A的權(quán)重，其中W表示點到A的引力大小，W為點到A的引力大小。

2.2 腹地準確度評價指標

腹地雖然沒有準確的范圍，但是依據(jù)不同的方法所得的范圍相互驗證有助于判斷不同方法的合理性。因此本文采用像素精度這一評價指標定量衡量本文方法的適宜性，像素精度定義如下：

式中為像素精度。為總的像素個數(shù)。為分類正確的像素個數(shù)，本文表示與其他方法腹地范圍一致的像素個數(shù)。

2.3 腹地均衡指數(shù)

收縮與擴張是城市空間結(jié)構(gòu)的演變規(guī)律[17]。城市收縮會伴隨著城市資源集約化分布，即資源向發(fā)達城市集中，再由發(fā)達城市向周邊城市提供資源。城市擴張同時伴隨著資源均衡化分布，即資源較均勻分配到每個城市，每個城市主要依靠自身資源均衡發(fā)展。

為探測城市發(fā)展規(guī)律，本文依據(jù)各市行政區(qū)范圍與各市腹地范圍差異程度，并結(jié)合相關(guān)系數(shù)量算方法，提出腹地均衡指數(shù)，用于探測城市集約與均衡度，公式如下：

3 結(jié)果分析

3.1 山東省各市腹地范圍

利用公式（3）逐個計算柵格像元到17個市級行政區(qū)的引力值，依據(jù)面Voronoi圖原理，取對其吸引力最大值像元所在的建成區(qū)作為該像元的腹地中心，單個像元的計算過程示例如圖4所示。

圖4 單個像元劃分示例

山東省各市2010、2011、2012、2013年17個市級行政區(qū)腹地范圍，如圖5所示。其中濟南、青島、煙臺、濰坊、臨沂相對于其他城市具有更廣闊的腹地。而萊蕪、泰安、德州、聊城、日照、威海等城市腹地范圍相對較小。隨著時間的推移，腹地范圍較大的城市腹地面積開始慢慢減小。相反，腹地較小的城市腹地范圍開始增大。按行政區(qū)劃分原則，各個行政區(qū)為自身范圍提供資源供應。然而2010至2013年濟南市的城市腹地范圍跨越了萊蕪，泰安，聊城三個市的反?，F(xiàn)象，打破各市的行政壁壘，主要是由于濟南周圍城市綜合實力較弱，導致其周邊城市的腹地范圍遠小于其自身行政區(qū)范圍。這與2018年萊蕪市被劃入濟南市的政策相符。

圖5 2010-2013年度腹地范圍分布圖

3.2 腹地劃分算法對比

基于交通可達性就近原則以及城市競爭力引力模型的腹地劃分方法[12]被廣泛使用在腹地提取中。但是該方法，在一些封閉網(wǎng)絡(luò)，如：鐵路，高度公路，中存在不合理的通行，進行封閉處理方法復雜。且該方法的腹地劃分結(jié)果明顯的細小空洞斑塊腹地（如圖6 d的局部放大圖所示），以及存在著腹地隨道路延伸等現(xiàn)象（如圖6 a局部放大圖所示），與實際情況不符。

圖6 基于交通可達性與城市競爭力腹地提取結(jié)果

為評價本文方法的合理性，將本文方法所得結(jié)果與僅有的像元尺度腹地劃分方法－基于交通可達性與城市競爭力腹地提取方法所得結(jié)果進行定量對比。依據(jù)公式（4）進行定量精度評價。精度結(jié)果如表1所示。結(jié)果表明本文方法與經(jīng)典方法在重疊度上均達到73%以上。由此說明本文方法具有一定的合理性。基于此結(jié)論以下章節(jié)將對結(jié)果進行分析。

表1 2010至2013年腹地范圍準確度

3.3 山東省各市腹地流失與擴張

腹地流失與腹地擴張能體現(xiàn)各個城市的經(jīng)濟實力與城市間的經(jīng)濟來往強度。本文基于腹地劃分的結(jié)果，分別統(tǒng)計2010年至2013年度的各市腹地流失與腹地擴張的百分比情況，結(jié)果如圖7。像元值為橫坐標城市行政區(qū)范圍流入縱坐標城市的面積占橫坐標城市的百分比。各列和約為100(考慮到邊界像元值的歸屬情況和小數(shù)位舍棄情況，存在小于等于1%的誤差)。結(jié)果表明：發(fā)達城市行政區(qū)范圍均位于該城市的腹地范圍，如：濟南2010-2013年度腹地流出百分比僅為1%至3%。而經(jīng)濟規(guī)模較小城市腹地范圍僅占該行政區(qū)范圍的40%上下，如萊蕪，日照等城市。

其中腹地流失最為嚴重的城市萊蕪2010年至2013年間分別有44%、62%、45%、38%的行政區(qū)范圍流入濟南市。表明這兩個城市間具有較強的經(jīng)濟聯(lián)系強度，濟南對萊蕪具有很強的輻射作用，這與2018年萊蕪市劃分到濟南市，設(shè)立濟南萊蕪區(qū)的政策相符。其次為德州，因與濟南交界，且經(jīng)濟較為落后，經(jīng)濟中心位于該市邊緣地區(qū)，導致其中心城區(qū)吸引力范圍較小，因此2010年至2013年間腹地流出至濟南市分別為自身范圍的58%、56%、51%、41%。相反，濟南市作為省會城市，因其較強的經(jīng)濟實力，在2012年達到腹地面積的最大值，侵占了其周邊城市，如：泰安行政區(qū)面積的49%，萊蕪行政區(qū)面積的62%，德州行政區(qū)面積的56%，聊城行政區(qū)面積的42%之多。僅次于濟南的臨沂，因其周邊城市的經(jīng)濟實力較弱，導致其腹地范圍也十分廣闊，并于2013年達到最大，其侵占了棗莊行政區(qū)面積的34%，日照行政區(qū)面積的59%之多。

圖7 2010-2013年度腹地面積占行政區(qū)面積百分比

3.4 腹地均衡指數(shù)變化分析

整個研究區(qū)域每個市的腹地面積與行政區(qū)面積差異程度能體現(xiàn)該研究區(qū)的經(jīng)濟發(fā)展情況，以及資源分配情況。為此本文提出腹地均衡指數(shù)來衡量整個研究區(qū)域的腹地平衡指數(shù)，腹地指數(shù)計算方式如公式(5)。

分別統(tǒng)計各市2010年、2011年、2012年、2013年四個年度行政區(qū)面積與各市腹地面積，結(jié)果如圖8。濟南作為省會城市腹地面積遠遠高于其行政區(qū)面積，這與實際情況相符。其次青島、濰坊、臨沂等城市腹地面積也略高于其行政區(qū)面積，這表面這些城市相對于其周邊城市具有更強的資源輻射程度。而德州、泰安、日照等城市行政區(qū)面積卻遠大于其腹地面積。表明德州等城市資源增長量與本市需求增加量不符，需要其他城市提供相應的資源供應。其中，濟南因其發(fā)達的經(jīng)濟情況，以及周邊城市經(jīng)濟較為落后等原因，腹地范圍十分廣闊，其中2012年腹地范圍達到最大，為自身行政區(qū)范圍的367%，而后兩年腹地范圍逐漸變小，分別下降至自身行政區(qū)范圍283%和198%。與此同時青島雖經(jīng)濟發(fā)達，但其周邊沿海城市經(jīng)濟也較為發(fā)達，因此青島的腹地范圍于2010年達到最大，為自身范圍的164%。而后幾年分別在123%至127%之間。

圖8 2010-2013年行政區(qū)面積與腹地面積散點圖

腹地均衡指數(shù)的變化趨勢體現(xiàn)了該研究區(qū)的發(fā)展均衡情況。本文基于腹地范圍與行政區(qū)范圍，結(jié)合公式（5），計算2010年至2013年度整個山東省腹地均衡指數(shù)，結(jié)果如圖9。山東省的腹地指數(shù)經(jīng)歷了穩(wěn)定小幅度上升，到急劇上升的階段。其中2010年至2011年度腹地均衡指數(shù)較為穩(wěn)定，而2011-2013年腹地均衡指數(shù)呈現(xiàn)出較強的增長趨勢，分別上升至0.262和0.635。表明山東省各市隨著時間推移，腹地均衡系數(shù)不斷增加。由此可知各市的資源不斷均衡化，發(fā)達城市腹地范圍不斷縮小，落后城市腹地范圍不斷擴大，各城市之間的差距逐漸變小，城市趨向于資源均衡化發(fā)展。上述結(jié)果表明山東半島藍色經(jīng)濟區(qū)的劃定對山東各市的腹地均衡性具有顯著的影響。

圖9 2010-2013年度腹地均衡指數(shù)

山東半島藍色經(jīng)濟區(qū)各市經(jīng)濟較為發(fā)達，因此本文分別統(tǒng)計了山東半島藍色經(jīng)濟區(qū)與山東省其他城市的腹地均衡系數(shù)，結(jié)果如圖10。結(jié)果表明山東半島藍色經(jīng)濟區(qū)的腹地均衡指數(shù)遠遠高于山東省其它城市。其中山東半島藍色經(jīng)濟區(qū)腹地均衡指數(shù)經(jīng)歷了“中速上升-低速穩(wěn)步上升”的動態(tài)發(fā)展過程。在2010年至2013年間，由0.672上升至0.882，其中2010年至2011年上升速率較大，達到0.149/年，而后2011年至2013年間上升平緩，上升速率僅為0.305/年。除山東半島藍色經(jīng)濟區(qū)以外其他城市的均衡指數(shù)呈現(xiàn)出“小幅下降-小幅上升-高速上升”的動態(tài)發(fā)展過程，其中2010年至2011年呈現(xiàn)小幅度下降趨勢，由0.145下降至0.1，2011年至2012年間由0.1小幅上升至0.17，而后2012年至2013年間指數(shù)高速上升至0.528，與山東半島藍色經(jīng)濟區(qū)的指數(shù)值差距達到最小，僅為0.354，2011年腹地指數(shù)差異度達到最大值0.721。無論山東半島藍色經(jīng)濟區(qū)還是山東其它城市，2011年之后較之前都呈現(xiàn)不同的發(fā)展趨勢。因此表明，山東半島藍色經(jīng)濟區(qū)的劃定對山東各市的腹地范圍變化具有很顯著的影響。

圖10 山東半島藍色經(jīng)濟區(qū)及其他城市2010-2013腹地指數(shù)

4 結(jié)論

本文提出一種基于遙感數(shù)據(jù)的腹地劃分新算法，該算法采用夜間燈光數(shù)據(jù)，借助引力模型，面加權(quán)Voronoi圖等方法，劃分2010至2013年四年間山東省各市腹地范圍，以及分析各市腹地范圍變化情況。并創(chuàng)新性的提出腹地均衡指數(shù)衡量腹地均衡程度，主要結(jié)論包括：

（1）本文方法和與基于交通可達性與城市競爭力腹地提取方法結(jié)果進行對比，結(jié)果顯示4年間兩總方法的腹地范圍分別達到73%至75%之間。并且本文方法所提取的腹地結(jié)果不存在細小斑塊空洞腹地以及按道路方向不合理延伸的情況；

（2）2010年至2013年間濟南、青島、煙臺、臨沂等城市腹地范圍遠遠大于自身的行政區(qū)范圍，其中2010年濟南腹地為自身范圍的3倍之多。萊蕪、泰安、德州、聊城、日照、威海等城市腹地被周邊發(fā)達省市侵占嚴重，其中萊蕪2011年有62%行政區(qū)范圍被濟南侵占。城市自身經(jīng)濟發(fā)達以及周邊城市經(jīng)濟較弱是城市腹地擴張的重要條件。設(shè)立山東半島藍色經(jīng)濟區(qū)的劃定使得山東省腹地的均衡度呈現(xiàn)顯著性的上升趨勢；

（3）山東半島藍色經(jīng)濟區(qū)腹地平衡指數(shù)明顯高于山東省其他城市，2011年至2013年間山東半島藍色經(jīng)濟區(qū)均衡指數(shù)上升趨勢放緩，增速由0.15/年降至0.06/年上下。然而，除山東半島藍色經(jīng)濟區(qū)以外其他市級行政區(qū)2010年至2011年腹地均衡指數(shù)有緩慢的下降趨勢，2011年至2013年卻以0.214/年的速率上升。山東省整個地區(qū)腹地平衡指數(shù)呈現(xiàn)上升趨勢。表明山東各市差距不斷縮小，各市發(fā)展呈現(xiàn)出越來越均衡的趨勢。其中2011年至2013年均衡發(fā)展趨勢最為顯著。

[1] 彭建,陳云謙,胡智超,等.城市腹地定量識別研究進展與展望[J].地理科學進展,2016,35(1):14-24

[2] 約翰·馮·杜能.孤立國同農(nóng)業(yè)和國民經(jīng)濟的關(guān)系[M]北京:商務印書館,1986:103

[3] 韋伯.工業(yè)區(qū)位論[M]北京:商務印書館,2010:47

[4] 周一星,許學強,寧越敏.城市地理學[M]第2版.北京:高等教育出版社,2009:39

[5] Green H. Hinterland boundaries of New York City and Boston in Southern New England [J]. Economic Geography, 1955,31(4):283-300

[6] Simon XB Zhao. Spatial restructuring of financial centers in mainland China and Hong Kong: A geography of finance perspective [J]. Urban Affairs Review, 2003,38(4):535-571

[7] Wheeler JO, Mitchelson RL. Information flows among major metropolitan areas in the United States [J]. Annals of the Association of American Geographers, 1989,79(4):523-543

[8] Dennis C, Marsland D, Cockett T. Central place practice: shopping centre attractiveness measures, hinterland boundaries and the UK retail hierarchy [J].Journal of Retailing and Consumer Services, 2002,9(4):185-199

[9] 王法輝,金鳳君,曾光.區(qū)域人口密度函數(shù)與增長模式:兼論城市吸引范圍劃分的GIS方法[J].地理研究,2004,23(1):97-103

[10] 閆衛(wèi)陽,秦耀辰,郭慶勝,等.城市斷裂點理論的驗證、擴展及應用[J].人文地理,2004(2):12-16

[11] 鐘業(yè)喜,陸玉麒.基于空間聯(lián)系的城市腹地范圍劃分——以江蘇省為例[J].地理科學,2012,32(5):536-543

[12] 曾冰,張艷,胡亞光.基于交通可達性與城市競爭力的城市腹地范圍識別——以長江中游城市群為例[J].熱帶地 理,2020,40(1):119-127

[13] 趙曉偉,孟治國,鄭興明,等.夜間燈光數(shù)據(jù)的東三省城市群發(fā)展特征研究[J].測繪科學,2019,44(4):72-80

[14] 范俊甫,何惠馨,胡桃英,等.一種敘利亞戰(zhàn)爭城市破壞程度的評估方法[J].測繪科學,2018,43(9):129-134

[15] Ji GX, Tian L, Zhao JC,. Detecting spatiotemporal dynamics of PM 2.5 emission data in China using DMSP-OLS nighttime stable light data [J]. Journal of Cleaner Production, 2019,209(3):363-370

[16] 國家發(fā)展改革委.山東半島藍色經(jīng)濟區(qū)規(guī)劃[EB]. 2011. http://wwwndrcgovcn/fzgggz/dqjj/zhdt/201408/t20140827_623621html.

[17] 班茂盛,方創(chuàng)琳.國內(nèi)城市邊緣區(qū)研究進展與未來研究方向[J].城市規(guī)劃學刊,2007(3):49-54

Urban Hinterland Extraction and Change Analysis in Shandong Province Based on Stable Night Light Data

YU Xiao1, LI Yi1, FAN Jun-fu1*, HAN Liu-sheng1, ZHOU Yu-ke2*

1.255049,2.100101,

Urban hinterland extraction is an important method to explore the scope of urban services and measure the degree of inter-city linkages. For the urban hinterland extraction, this paper propose an new method which based on the stable version of DMSP/OLS data, using iterative threshold segmentation method to extract high-quality city range and night light area as the city quality calculation index, night light total value proportion as the weight coefficient, combined with surface weighting Voronoi Graphs, gravitational models and other methods to extract the hinterland range of cities from 17 prefecture-level cities in Shandong Province from 2010 to 2013 on a pixel scale. Analysis of changes in the hinterland range of each city over the four years. And propose the hinterland balance coefficient to explore the direction of urban development. The research shows that: ① The overlap of the results of this method and the traditional method is over 73%, beside more solve some small patch problem and hinterland extend with road problem. ②Laiwu and Taian are surrounded by the hinterland of Jinan, and the hinterland breaks the administrative boundary barrier. ③From 2010 to 2013, the balance index of the hinterland of each city in Shandong has been increasing, with obvious geographical rules of urban resources from intensification to equilibrium. ④The balance index of each city in the Shandong Peninsula Blue Economic Zone is significantly higher than that of other cities in Shandong. From 2010 to 2013, the balance index of the hinterland of Shandong Peninsula increased from 0.672 to 0.882, and the balance index of other cities in the hinterland increased from 0.145 to 0.528.

Light data; Urban construction; spatial distribution

TU984.11+3

A

1000-2324(2021)05-0870-10

2020-06-09

2020-11-23

國家自然科學基金項目（42171413）;山東省自然科學基金項目（ZR2020MD0115、ZR2020MD018）;廣東省科學院建設(shè)國內(nèi)一流研究機構(gòu)行動專項資金項目(2019GDASYL-0103003);山東理工大學青年教師發(fā)展支持計劃項目(4072-115016)

俞宵(1994-),男,碩士研究生,主要從事智能地理計算研究. E-mail:1126876778@qq.com

通訊作者:Author for correspondence. E-mail:fanjf@sdut.edu.cn; zhouyk@igsnrr.ac.cn