基于夜間燈光數(shù)據(jù)的山東省城市擴(kuò)展時(shí)空格局

張智璇單寶艷王俊凝劉洋洋林琪凱

(山東建筑大學(xué)測(cè)繪地理信息學(xué)院,山東濟(jì)南250101)

0 引言

城市是人口、資源、環(huán)境和各種社會(huì)經(jīng)濟(jì)因素高度集合的綜合體,在區(qū)域經(jīng)濟(jì)發(fā)展進(jìn)程中發(fā)揮著核心作用,城市空間擴(kuò)展直接帶動(dòng)區(qū)域空間結(jié)構(gòu)的改變。國(guó)內(nèi)外學(xué)者對(duì)城市擴(kuò)展的研究一般基于傳統(tǒng)統(tǒng)計(jì)數(shù)據(jù)以及Landsat、SPOT等遙感數(shù)據(jù)[1-2],主要運(yùn)用神經(jīng)網(wǎng)絡(luò)分類、監(jiān)督分類、歸一化建筑指數(shù)、極限學(xué)習(xí)機(jī)、構(gòu)造CLUE-S模型、城市年輪模型等方法對(duì)城市擴(kuò)張進(jìn)行模擬分析[3-14],研究?jī)?nèi)容主要集中于城市擴(kuò)展的特征、驅(qū)動(dòng)力、空間格局和演變格局等方面。基于傳統(tǒng)統(tǒng)計(jì)數(shù)據(jù)的城市擴(kuò)展研究存在著時(shí)效性較差等問(wèn)題。1978年,Croft首次提出運(yùn)用DMSP/OLS夜間燈光數(shù)據(jù)進(jìn)行城市建成區(qū)提取[15]。自此以后,夜間燈光數(shù)據(jù)在城市用地提取方面的研究受到了廣泛關(guān)注。夜間燈光數(shù)據(jù)為提取城市建設(shè)用地信息[16]、測(cè)算城市發(fā)展空間格局演變[17-22]以及推測(cè)城市城鎮(zhèn)化水平[23]等提供了可靠的數(shù)據(jù)來(lái)源。

夜間燈光數(shù)據(jù)既能提取城市建成區(qū),也可反映區(qū)域經(jīng)濟(jì)發(fā)展水平。有關(guān)學(xué)者采用相關(guān)分析和回歸分析等方法,對(duì)夜間燈光總輻射率與人口總數(shù)[24-31]、GDP數(shù)據(jù)[32-33]之間的關(guān)系進(jìn)行分析,以城市為統(tǒng)計(jì)單元,發(fā)現(xiàn)夜間燈光亮度總值和統(tǒng)計(jì)人口總數(shù)、GDP等經(jīng)濟(jì)因子之間呈現(xiàn)顯著的相關(guān)性,夜間燈光影像已成為反演社會(huì)經(jīng)濟(jì)發(fā)展水平的良好數(shù)據(jù)源。因此,與傳統(tǒng)方法相比,采用不同年份夜間燈光數(shù)據(jù)提取城市建成區(qū)面積不僅可以反映城市建成區(qū)面積的變化,還可以反映城市建設(shè)質(zhì)量、城鎮(zhèn)人口、經(jīng)濟(jì)發(fā)展水平等內(nèi)涵特征。目前,山東省正處于城市化高速發(fā)展階段,城市范圍不斷擴(kuò)展,由于各市縣區(qū)位條件和經(jīng)濟(jì)發(fā)展水平的差異,各地城市化水平存在非均質(zhì)性,城市擴(kuò)展規(guī)模、強(qiáng)度和質(zhì)量存在一定差異。文章基于傳統(tǒng)統(tǒng)計(jì)數(shù)據(jù)和夜間燈光數(shù)據(jù),對(duì)山東省各市建成區(qū)面積擴(kuò)展與內(nèi)涵增長(zhǎng)進(jìn)行綜合分析,分析其擴(kuò)展的空間格局和發(fā)展演變特點(diǎn),對(duì)全面分析山東省各地城鎮(zhèn)化水平和區(qū)域空間結(jié)構(gòu)具有重要意義。

1 數(shù)據(jù)來(lái)源及研究方法介紹

山東省是經(jīng)濟(jì)大省,自2007年以來(lái),經(jīng)濟(jì)總量穩(wěn)居全國(guó)第3位。山東省內(nèi)城市發(fā)展較快,2018年,山東省城市化率達(dá)到61.18%,山東半島部分地區(qū)達(dá)到了70%以上。據(jù)統(tǒng)計(jì),2001—2015年間,山東省城市建成區(qū)面積由119 200 km2增長(zhǎng)到339 800 km2,年均增長(zhǎng)率超過(guò)7%,明顯超過(guò)全國(guó)的平均增長(zhǎng)率6.21%。

1.1 數(shù)據(jù)來(lái)源

2007、2012和2017年的山東省各城市建成區(qū)面積及各經(jīng)濟(jì)指標(biāo)來(lái)源于對(duì)應(yīng)年份的《山東統(tǒng)計(jì)年鑒》和《山東省建設(shè)統(tǒng)計(jì)年報(bào)》,地圖矢量文件來(lái)源于山東天地圖(www.sdmap.gov.cn)。研究區(qū)以2017年底的山東省行政區(qū)劃為基準(zhǔn),包括99個(gè)市縣。各年份城市建成區(qū)統(tǒng)計(jì)面積數(shù)據(jù)處理時(shí),涉及行政區(qū)域變更的全部統(tǒng)一到2017年底的行政區(qū)劃。膠南市調(diào)整為黃島區(qū)歸并入青島市,即墨市調(diào)整為即墨區(qū)歸并入青島市,章丘市調(diào)整為章丘區(qū)歸并入濟(jì)南市,墾利縣調(diào)整為墾利區(qū)歸并入東營(yíng)市,兗州市調(diào)整為兗州區(qū)歸并入濟(jì)寧市,文登市調(diào)整為文登區(qū)歸并入威海市,陵縣調(diào)整為陵城區(qū)歸并入德州市,沾化市調(diào)整為沾化區(qū)歸并入濱州市,定陶縣調(diào)整為定陶區(qū)歸并入菏澤市,蒼山縣更名為蘭陵縣,文中地圖審圖號(hào)為:魯SG(2019)066號(hào)。

研究使用的兩種夜間燈光數(shù)據(jù)均來(lái)源于NOAA網(wǎng)站(http://www.ngdc.noaa.gov/),分別為美國(guó)軍事氣象衛(wèi)星計(jì)劃發(fā)射的DMSP衛(wèi)星通過(guò)搭載OLS傳感器獲得的全球夜間燈光影像數(shù)據(jù),以及NPOESS Preparatory Project衛(wèi)星系統(tǒng)搭載傳感器得到的全球夜間燈光數(shù)據(jù),時(shí)間序列為2007、2012和2017年。2007和 2012年為 DMSP/OLS影像,2017年為NPP/VIIRS影像,設(shè)定投影坐標(biāo)系為 WGS 84/UTM zone 51N,以2017年底山東省行政區(qū)劃為邊界對(duì)全球夜間燈光數(shù)據(jù)進(jìn)行裁剪和地理配準(zhǔn)。

1.2 研究方法

1.2.1 相關(guān)分析

相關(guān)分析是研究變量之間相關(guān)性的分析方法,是使用適當(dāng)?shù)慕y(tǒng)計(jì)指標(biāo)來(lái)描述客觀事物之間關(guān)系緊密程度的過(guò)程。對(duì)于變量X=(x1,x2,x3,…,xi)、Y=(y1,y2,y3,…,yi),其相關(guān)系數(shù)rxy由式(1)[34]表示為

式中:、分別為變量X、Y的平均值。相關(guān)系數(shù)rxy介于-1和1之間,為正值時(shí),兩個(gè)變量為正相關(guān)；為負(fù)值時(shí),兩個(gè)變量為負(fù)相關(guān),r的絕對(duì)值越接近1,兩個(gè)變量之間的相關(guān)性越強(qiáng),越接近0,兩個(gè)變量之間的關(guān)聯(lián)程度越弱；而為0時(shí),兩變量無(wú)線性相關(guān)。

1.2.2 熱點(diǎn)分析

熱點(diǎn)分析又稱Getis-OrdGi?統(tǒng)計(jì),其計(jì)算由式(2)[35]表示為式中:xj為要素j的屬性值；wi,j為要素i和j之間的空間權(quán)重；n為要素總數(shù)；和S分別由式(3)和(4)表示為

Getis-OrdGi?統(tǒng)計(jì)可使用ArcGIS10.3計(jì)算,計(jì)算結(jié)果輸出Z得分和P值,統(tǒng)計(jì)即為Z得分[36]。通過(guò)2個(gè)數(shù)值可以得到高值和低值要素在空間上發(fā)生聚類的位置。Z為正值表現(xiàn)為熱點(diǎn),為負(fù)值表現(xiàn)為冷點(diǎn),并且其絕對(duì)值越大,在空間上聚類效果越顯著。Z＜-1.65或>1.65時(shí),其要素反映置信度為90% 的統(tǒng)計(jì)顯著性；Z＜-1.96或>1.96時(shí),其要素反映置信度為95% 的統(tǒng)計(jì)顯著性；Z＜-2.58或>2.58時(shí),其要素反映置信度為99% 的統(tǒng)計(jì)顯著性；Z得分接近0時(shí),表示沒(méi)有顯著的空間聚類。

2 基于夜間燈光數(shù)據(jù)的城市建成區(qū)面積提取

2.1 城市建成區(qū)面積提取

采用夜間燈光數(shù)據(jù)提取城市建成區(qū)面積,經(jīng)常采用閾值法[37]。基于夜間燈光數(shù)據(jù)提取建成區(qū)面積的標(biāo)準(zhǔn)為:使其滿足某一閾值的夜間燈光數(shù)據(jù)提取的總面積與當(dāng)年的《山東統(tǒng)計(jì)年鑒》和《山東省建設(shè)統(tǒng)計(jì)年報(bào)》統(tǒng)計(jì)的全省城市建成區(qū)面積總和的差值最小,并且滿足兩者之間的誤差平方和達(dá)到最小。具體提取方法為:

(1)假設(shè)某年山東省夜間燈光度的灰度最大值為DNmax、最小值為DNmin,山東省建成區(qū)面積的潛在變量閾值DNx,定義變量DNx為灰度值范圍的平均值；

(2)提取閾值為DNx時(shí)的面積S1,將其與當(dāng)年統(tǒng)計(jì)年鑒統(tǒng)計(jì)的建成區(qū)面積總和S2作比較,若S1>S2,則將上述變量DNx作為第2次計(jì)算的最大值,替換原DNmax,若S1＜S2,則將上述變量DNx作為第2次計(jì)算的最小值,替換原DNmin；

(3)重復(fù)以上步驟進(jìn)行計(jì)算,直到S1≈S2,且滿足誤差平方和最小為止,則該DNx作為該年份建成區(qū)面積提取的最佳閾值,在此閾值下提取的面積即作為建成區(qū)面積。

基于以上方法,選取2017年的提取閾值為8.9,即以DNx>8.9提取2017年的山東省各城市建成區(qū)面積。為分析各城市建成區(qū)提取面積(M1)與統(tǒng)計(jì)面積(M2)的誤差及其空間差異,計(jì)算提取誤差W=(M1-M2)/M2,并運(yùn)用ArcGIS 10.3做空間分布分析,采用分位數(shù)法進(jìn)行分類,結(jié)果如圖1所示。采用同樣方法提取2007和2012年的建成區(qū)面積,并分析其提取誤差和空間差異,結(jié)果與2017年類似。

圖1 2017年山東省城市建成區(qū)面積提取誤差圖

由圖1可見(jiàn),山東省各城市建成區(qū)提取面積和統(tǒng)計(jì)面積均存在一定差異,其中城市建成區(qū)提取面積明顯大于統(tǒng)計(jì)面積的主要為地級(jí)以上城市,建成區(qū)提取面積小于統(tǒng)計(jì)面積的全部為縣級(jí)市或縣城。

2.2 城市建成區(qū)提取面積與經(jīng)濟(jì)指標(biāo)的關(guān)系

為分析山東省各城市夜間燈光數(shù)據(jù)和城市經(jīng)濟(jì)發(fā)展水平的相關(guān)性,對(duì)2007、2012、2017年基于夜間燈光數(shù)據(jù)提取的建成區(qū)面積與當(dāng)年的一般公共預(yù)算收入、一般公共預(yù)算支出、規(guī)模以上工業(yè)企業(yè)總產(chǎn)值、社會(huì)消費(fèi)品零售額等經(jīng)濟(jì)指標(biāo)做相關(guān)分析。采用式(1),運(yùn)用Matlab2014計(jì)算相關(guān)系數(shù)R及其對(duì)應(yīng)的P值(兩者不相關(guān)的概率),結(jié)果見(jiàn)表1。

表1 2007、2012、2017年山東省各市縣夜間燈光提取面積與經(jīng)濟(jì)指標(biāo)的相關(guān)系數(shù)表

由表1可知,2007、2012、2017年的基于夜間燈光提取的建成區(qū)面積與各經(jīng)濟(jì)指標(biāo)均為正相關(guān),且高度顯著。因此,山東省各城市夜間燈光數(shù)據(jù)不僅可以反映建成區(qū)面積,而且還能夠反映城市經(jīng)濟(jì)和城市化質(zhì)量。

3 山東省城市擴(kuò)展的時(shí)空格局分析

城市擴(kuò)展既包括量的擴(kuò)展,表現(xiàn)為城市建成區(qū)面積的擴(kuò)張；也包括質(zhì)的提升,表現(xiàn)為城市的內(nèi)涵增長(zhǎng)、功能提升、結(jié)構(gòu)優(yōu)化和建設(shè)質(zhì)量與水平的提高等。城市發(fā)展過(guò)程中,在擴(kuò)張面積相同的情況下,城市建設(shè)投入和開(kāi)發(fā)強(qiáng)度不同,城市擴(kuò)展質(zhì)量必然存在差異。對(duì)于城市擴(kuò)展水平的定量分析,僅基于單一的城市空間擴(kuò)展面積進(jìn)行分析就比較片面,而構(gòu)建綜合指標(biāo),從空間擴(kuò)展面積和擴(kuò)展質(zhì)量等方面進(jìn)行綜合分析,可以較為全面地反映城市擴(kuò)展綜合水平。借鑒已有研究成果[37-38],引入城市擴(kuò)展綜合指標(biāo),對(duì)城市擴(kuò)展質(zhì)和量進(jìn)行定量表達(dá),使其不僅能夠反映建成區(qū)的空間擴(kuò)展,也能反映城鎮(zhèn)人口、經(jīng)濟(jì)發(fā)展水平和建設(shè)質(zhì)量等方面。

3.1 城市擴(kuò)展綜合指標(biāo)構(gòu)建

夜間燈光數(shù)據(jù)能夠反映城市經(jīng)濟(jì)和城市質(zhì)量等方面,也能夠表達(dá)建成區(qū)面積,但各城市均存在一定誤差。建成區(qū)統(tǒng)計(jì)數(shù)據(jù)能夠比較精確地表達(dá)建成區(qū)面積,但難以反應(yīng)城市經(jīng)濟(jì)、城市建設(shè)質(zhì)量。因此,文章取兩者之長(zhǎng),以夜間燈光數(shù)據(jù)提取的建成區(qū)面積M1和建成區(qū)面積統(tǒng)計(jì)數(shù)據(jù)M2和構(gòu)建城市擴(kuò)展綜合指標(biāo)M,由式(5)表示為

式中:λ1和λ2分別為M1和M2的權(quán)重,采用專家調(diào)查法和靈敏度分析法綜合確定。

這樣構(gòu)建的城市擴(kuò)展綜合指標(biāo),既能較精確地表達(dá)建成區(qū)面積擴(kuò)展,同時(shí)還能反映城市經(jīng)濟(jì)和城市建設(shè)質(zhì)量等方面。分別采用λ1/λ2為0.5/0.5、0.4/0.6、0.3/0.7進(jìn)行靈敏度分析,計(jì)算結(jié)果差別較小,經(jīng)綜合分析,最終確定權(quán)重λ1和λ2分別為0.4和0.6。

3.2 山東省城市擴(kuò)展的空間格局

依據(jù)城市擴(kuò)展綜合指標(biāo)數(shù)據(jù)計(jì)算2007—2012年、2012—2017年的山東省各城市擴(kuò)展綜合指標(biāo)的增長(zhǎng)量和增長(zhǎng)率,并運(yùn)用ArcGIS10.3進(jìn)行空間格局分析。將增長(zhǎng)量和增長(zhǎng)率按分位法各分為3類,每一類各占1/3,分別表示低、中、高三級(jí)增長(zhǎng)率和低、中、高三級(jí)增長(zhǎng)量。其中,2007—2012年間,-0.32～0.19定義為低增長(zhǎng)率、0.19～0.48為中增長(zhǎng)率、0.48～1.42為高增長(zhǎng)率；-6.03～3.19定義為低增長(zhǎng)量、3.19～7.82為中增長(zhǎng)量、7.82～132.37為高增長(zhǎng)量。2012—2017年間,-0.34～0.14定義為低增長(zhǎng)率、0.14～0.37為中增長(zhǎng)率、0.37～1.69為高增長(zhǎng)率；-9.48～3.74定義為低增長(zhǎng)量、3.74～9.30為中增長(zhǎng)量、9.30～212.2為高增長(zhǎng)量。增長(zhǎng)量采用分級(jí)設(shè)色表示,增長(zhǎng)率采用線狀圖形表示,將圖層疊加得到增長(zhǎng)率和增長(zhǎng)量的綜合分析圖,既可分別表達(dá)山東省各城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量和增長(zhǎng)率的空間格局,也可表達(dá)增長(zhǎng)量和增長(zhǎng)率9種組合類型的空間格局,結(jié)果如圖2、3所示。

圖2 2007—2012年山東省城市擴(kuò)展綜合指標(biāo)增長(zhǎng)率和增長(zhǎng)量綜合圖

圖3 2012—2017年山東省城市擴(kuò)展綜合指標(biāo)增長(zhǎng)率和增長(zhǎng)量綜合圖

由圖2可知,2007—2012年,山東省城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量和增長(zhǎng)率的空間格局存在著明顯差異,增長(zhǎng)量較大的主要為山東省的各設(shè)區(qū)市,多為中增長(zhǎng)率,說(shuō)明各城市擴(kuò)展的增長(zhǎng)量主要取決于城市的原有規(guī)模。設(shè)區(qū)市原有城市規(guī)模基數(shù)較大,所以城市擴(kuò)展增長(zhǎng)量較大,而增長(zhǎng)率相對(duì)較小。城市擴(kuò)展增長(zhǎng)率高的市縣主要位于魯西北和魯西南地區(qū),其原有城市規(guī)模基數(shù)較小,所以城市擴(kuò)展率相對(duì)較高。其中,山東省東部的膠州市,中部的臨朐縣、昌樂(lè)縣、安丘市,北部的濱州市、博興縣、廣饒縣,西部的德州市、齊河縣、鄆城縣,南部的曹縣、濟(jì)寧市、臨沂市、曲阜市的城市擴(kuò)展綜合指標(biāo)增長(zhǎng)率和增長(zhǎng)量都較高,城市擴(kuò)展相對(duì)迅猛。南部的諸城市、莒縣、莒南縣、五蓮縣、臨沭縣、新泰市、寧陽(yáng)市、平邑縣、費(fèi)縣、滕州市、微山縣、成武縣、蘇明縣、郯城縣,東部的海陽(yáng)市、乳山市、萊陽(yáng)市、龍口市,北部的樂(lè)陵市、慶云縣、利津縣,西部的高唐縣、陽(yáng)谷縣、東阿縣的城市擴(kuò)展綜合指標(biāo)增長(zhǎng)率和增長(zhǎng)量都較低,城市擴(kuò)展相對(duì)緩慢。

由圖3可知,2012—2017年,山東省城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量和增長(zhǎng)率的空間格局仍然存在明顯差異,增長(zhǎng)量較大的仍為山東省的各設(shè)區(qū)市,多為低中增長(zhǎng)率,城市擴(kuò)展增長(zhǎng)率高的市縣主要位于魯西和魯西南地區(qū)。其中,山東省東北部的蓬萊市,南部的平邑縣、滕州市、沂南縣、莒南縣、郯城縣、滕州市,西部的菏澤市、鄆城縣、濟(jì)寧市、東明縣、曹縣、陽(yáng)谷縣、冠縣、聊城市、荏平縣、高唐縣等市縣的城市擴(kuò)展綜合指標(biāo)增長(zhǎng)率和增長(zhǎng)量都較高,城市擴(kuò)展相對(duì)迅猛。西北部的無(wú)棣縣、慶云縣、樂(lè)陵縣、商河縣、寧津縣、利津縣、高唐縣、惠民縣、夏津縣、武城縣、平原縣、商河縣,巨野縣,中部的壽光市、昌邑市、臨朐縣、曲阜市、肥城市,東部的高密市、萊州市、招遠(yuǎn)市、日照市、棲霞市、海陽(yáng)市、乳山市等市縣的城市擴(kuò)展綜合指標(biāo)增長(zhǎng)率和增長(zhǎng)量都較低,城市擴(kuò)展相對(duì)緩慢。

3.3 山東省城市擴(kuò)展的時(shí)空熱點(diǎn)分析

采用式(4),運(yùn)用Arcgis10.3的空間統(tǒng)計(jì)工具箱對(duì)山東省城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量和增長(zhǎng)率進(jìn)行熱點(diǎn)分析,空間權(quán)重采用“Queen”相鄰,即共享邊或點(diǎn)即為相鄰,選擇熱點(diǎn)分析工具的 CONTIGUITY_EDGES_CORNERS實(shí)現(xiàn)。 2007—2012年、2012—2017年山東省城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量和增長(zhǎng)率的熱點(diǎn)分析結(jié)果如圖4～7所示。

由圖4、5可知,2007—2012年、2012—2017年,山東省城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量的熱點(diǎn)分布較為明顯,且存在明顯的空間差異,均無(wú)顯著的冷點(diǎn)。2007—2012年,山東省城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量的熱點(diǎn)主要分布在省會(huì)濟(jì)南及其周邊市縣。其中,置信水平99%的熱點(diǎn)位于萊蕪市；置信水平95%的熱點(diǎn)位于淄博市、濟(jì)南市、泰安市、濟(jì)陽(yáng)縣和鄒平縣；置信水平90%的熱點(diǎn)位于齊河縣和平陰縣。說(shuō)明此期間濟(jì)南及周邊市縣的城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量呈現(xiàn)高值集聚。置信水平90%的熱點(diǎn)位于膠州市,其余城市冷熱點(diǎn)均不顯著。2012—2017年,山東省城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量的熱點(diǎn)主要分布在東部沿海的青島市、日照市及其相鄰市縣,其中,置信水平99%的熱點(diǎn)位于萊西市、膠州市、諸城市、五蓮縣和日照市；置信水平95%的熱點(diǎn)位于青島市、平度市、萊陽(yáng)市；其余市縣冷熱點(diǎn)均不顯著。說(shuō)明此期間青島市、日照市及其周邊市縣的城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量呈現(xiàn)高值集聚。

圖4 2007—2012年山東省城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量熱點(diǎn)分析圖

圖5 2012—2017年山東省城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量熱點(diǎn)分析圖

由圖6、7可知,2007—2012年、2012—2017年,山東省城市擴(kuò)展綜合指標(biāo)增長(zhǎng)率的冷熱點(diǎn)分布比較分散,空間差異明顯。2007—2012年,山東省城市擴(kuò)展綜合指標(biāo)增長(zhǎng)率的熱點(diǎn)主要分布在魯西北和魯西南的部分市縣,冷點(diǎn)主要分布在日照市、臨沂市及其周邊部分市縣。其中,置信水平95%的熱點(diǎn)位于武城縣、濟(jì)陽(yáng)縣、惠民縣、濟(jì)寧市、嘉祥縣、魚(yú)臺(tái)縣；置信水平90%的熱點(diǎn)位于臨邑縣、高青縣、高唐縣、鄆城縣、金鄉(xiāng)縣、鄄城縣；置信水平95%的冷點(diǎn)位于臨沭縣,置信水平90%的冷點(diǎn)位于東南部的五蓮縣、日照市、臨沂市、郯城縣；其余城市冷熱點(diǎn)均不顯著。說(shuō)明此期間魯西北和魯西南部分市縣的城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量呈現(xiàn)高值集聚,日照市、臨沂市及其周邊部分市縣的城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量呈現(xiàn)低值集聚。2012—2017年,山東省城市擴(kuò)展綜合指標(biāo)增長(zhǎng)率的熱點(diǎn)主要分布在魯西的聊城市及其周邊部分市縣、魯西南的部分市縣和臨沂市及其周邊的部分市縣,冷點(diǎn)主要分布在煙臺(tái)市及其周邊部分市縣及魯西北的陽(yáng)信縣。其中,置信水平99%的熱點(diǎn)位于冠縣、聊城市、莘縣、臨沂市和臨沭縣；置信水平95%的熱點(diǎn)位于臨清縣、陽(yáng)谷縣、梁山縣；置信水平90%的熱點(diǎn)位于東阿縣、鄆城縣、嘉祥縣、東明縣、棗莊市、莒南縣和郯城縣；置信水平95%的冷點(diǎn)位于山東省東部的乳山市和海陽(yáng)市；置信水平90%的冷點(diǎn)位于煙臺(tái)市；其余市縣冷熱點(diǎn)均不顯著。兩段時(shí)期冷熱點(diǎn)變化最明顯的是臨沂市、臨沭縣和郯城縣,城市擴(kuò)展綜合指標(biāo)增長(zhǎng)率由2007—2012年的冷點(diǎn)轉(zhuǎn)變?yōu)?012—2017年的熱點(diǎn)。

圖7 2012—2017年山東省城市擴(kuò)展綜合指標(biāo)增長(zhǎng)率熱點(diǎn)分析圖

4 結(jié)論

通過(guò)研究,得到以下結(jié)論:

(1)2007—2012、2012—2017 年,山東省城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量和增長(zhǎng)率的空間格局存在明顯差異,增長(zhǎng)量高低主要取決于城市的原有規(guī)模,增長(zhǎng)量較高的主要為設(shè)區(qū)城市,其城市原有規(guī)模較大,城市擴(kuò)展增長(zhǎng)量也較大。增長(zhǎng)率較高的市縣主要位于山東省經(jīng)濟(jì)發(fā)展水平相對(duì)落后的魯西、魯西北和魯西南地區(qū)。這些地區(qū)城市規(guī)模基數(shù)較小,增長(zhǎng)量較小,增長(zhǎng)率相對(duì)較高。2007—2012年,山東省城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量高的市縣,其增長(zhǎng)率主要處于中等水平；2012—2017年,城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量高的市縣,位于魯西和魯西南的部分市縣中其增長(zhǎng)率仍處于高水平。

(2)2007—2012、2012—2017年,山東省城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量和增長(zhǎng)率的熱點(diǎn)和冷點(diǎn)的地區(qū)分布明顯不同。城市擴(kuò)展綜合指標(biāo)增長(zhǎng)量的熱點(diǎn)一直存在集聚特征,隨著時(shí)間推移,在空間上發(fā)生了明顯變化,發(fā)展路徑為“中部地區(qū)—東部沿海地區(qū)”。增長(zhǎng)率的熱點(diǎn)和冷點(diǎn)一直較為分散,熱點(diǎn)主要由山東省西北部、西南部轉(zhuǎn)移到西部、南部,冷點(diǎn)主要由西南部轉(zhuǎn)移到東北部。