劉石俊，薛亮

(陜西師范大學 地理科學與旅游學院，西安 710119)

0 引言

區(qū)域經(jīng)濟格局研究實質(zhì)上是對區(qū)域經(jīng)濟在時空角度上的差異研究。區(qū)域經(jīng)濟差異是指一個區(qū)域相對于另一區(qū)域在經(jīng)濟發(fā)展上的不平衡、不平等狀況。區(qū)域經(jīng)濟發(fā)展格局及區(qū)域間經(jīng)濟差異一直是學術界十分關注的問題[1-3]。改革開放以來，陜西省經(jīng)濟持續(xù)發(fā)展，總體水平不斷攀升，但由于自然資源、地理位置、政策傾斜等原因，區(qū)域間的經(jīng)濟差異始終存在[4]。區(qū)域間的這種差異會阻礙陜西省經(jīng)濟可持續(xù)發(fā)展，并由此引發(fā)一系列的社會問題[5]。實現(xiàn)區(qū)域經(jīng)濟協(xié)調(diào)發(fā)展是經(jīng)濟健康運行的前提條件，也是實施精準扶貧政策全面進入小康社會的重要基礎[6]。伴隨著陜西省經(jīng)濟實現(xiàn)跨越發(fā)展的同時，發(fā)展的不平衡不充分、資源環(huán)境破壞、產(chǎn)業(yè)結(jié)構不發(fā)達等問題日益突出，因此，研究陜西省經(jīng)濟發(fā)展時空格局變化對于推動經(jīng)濟高質(zhì)量發(fā)展、持續(xù)健康發(fā)展具有極其重要的現(xiàn)實意義。

國內(nèi)外學者對區(qū)域經(jīng)濟差異問題進行了大量研究。目前針對經(jīng)濟差異的研究大多采用以國內(nèi)生產(chǎn)總值為主的社會統(tǒng)計數(shù)據(jù)，如Yamamoto[7]采用GDP數(shù)據(jù)基于變異系數(shù)、基尼系數(shù)、泰爾指數(shù)方法分析了美國地區(qū)收入差異；Fleisher等[8]通過分析GDP數(shù)據(jù)研究了中國沿海地區(qū)和非沿海地區(qū)在收入、生產(chǎn)率、政策方面的差異；周玉翠等[9]以人均GDP為測度指標研究了中國省際經(jīng)濟差異動態(tài)變化特征。由于統(tǒng)計數(shù)據(jù)的收集需要大量的人力和時間成本，且存在數(shù)據(jù)缺失、空間尺度不統(tǒng)一等問題，隨著遙感技術的迅猛發(fā)展和廣泛應用，夜間燈光數(shù)據(jù)逐漸被許多學者引入到經(jīng)濟差異問題研究中來。美國國防氣象衛(wèi)星(defense meteorological satellite program，DMSP)搭載的業(yè)務型線掃描傳感器(operational linescan system，OLS)生產(chǎn)的DMSP/OLS夜間燈光數(shù)據(jù)是當前應用最廣泛的夜間燈光數(shù)據(jù)之一，已被廣泛應用于人口空間化[10]、城市擴張[11]、碳排放[12]等領域。在經(jīng)濟發(fā)展變化應用方面，Marx等[13]研究發(fā)現(xiàn)，GDP與DMSP/OLS夜間燈光數(shù)據(jù)DN值間存在較強的相關性；葛心韻[14]利用DMSP/OLS夜間燈光影像構建預測模型進行了江蘇各地級市GDP預測研究；王毓乾等[15]采用一種新的去飽和方法，對處理后的影像進行省域生產(chǎn)總值反演。

從現(xiàn)有文獻分析，前人基于DMSP/OLS數(shù)據(jù)所作的有關經(jīng)濟發(fā)展方面的研究多聚焦于對GDP、貧困度方面的模型構建、預測和響應，有關連續(xù)時間內(nèi)的陜西省經(jīng)濟發(fā)展時空格局特征研究還比較少，且陜西省經(jīng)過這些年的發(fā)展，經(jīng)濟格局也發(fā)生了較大變化，非常有必要從學術角度來分析研究。因此，首先，對覆蓋陜西省的2001—2013年DMSP/OLS夜間燈光數(shù)據(jù)進行預處理以提升數(shù)據(jù)質(zhì)量；然后，采用重心、標準差橢圓、基尼系數(shù)、變異系數(shù)、局部莫蘭指數(shù)對2001—2013陜西省經(jīng)濟發(fā)展時空格局特征進行綜合分析。

1 研究區(qū)概況與數(shù)據(jù)來源

1.1 研究區(qū)概況

陜西東鄰山西、河南，西連寧夏、甘肅，南抵四川、重慶、湖北，北接內(nèi)蒙古，總面積20.56×104km2。由于南北方向橫跨范圍大，陜西省內(nèi)部地形氣候差異極大，陜北為黃土高原區(qū)，屬北亞熱帶氣候；關中為平原區(qū)，屬暖溫帶氣候；陜南為秦巴山區(qū)，屬中溫帶氣候。自改革開放以來，陜西省經(jīng)濟發(fā)展水平不斷提高，1979—2005年間年平均經(jīng)濟增速一直高于10%，且從1996年后經(jīng)濟增速一直高于全國平均水平。2001年，陜西省GDP總量僅1 768億元，到2013年，陜西省GDP總量已達到16 132億元。在陜西省經(jīng)濟快速發(fā)展的同時，陜西省107個縣級行政單元，各縣之間經(jīng)濟發(fā)展極不平衡。以2013年數(shù)據(jù)為例，西安市雁塔區(qū)GDP為988.2億元，漢中市佛坪縣GDP僅為5.79億元。另外，由于自然條件、經(jīng)濟發(fā)展、政策傾向等原因，陜西省是全國貧困面大、貧困人口多、貧困程度深的省份之一。在全國14個集中連片特困地區(qū)中，陜西省涉及秦巴山片區(qū)、六盤山片區(qū)、呂梁山片區(qū)。

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

本文數(shù)據(jù)主要包括統(tǒng)計數(shù)據(jù)、矢量數(shù)據(jù)、柵格數(shù)據(jù)。統(tǒng)計數(shù)據(jù)為2001—2013年陜西省縣域GDP和人口數(shù)據(jù)，來源于《陜西統(tǒng)計年鑒》，部分缺失數(shù)據(jù)由陜西各市統(tǒng)計年鑒補齊。矢量數(shù)據(jù)為2015年陜西省各縣行政區(qū)劃矢量，來源于全國地理信息資源目錄服務系統(tǒng)(http：//www.webmap.cn)中的1∶100萬全國基礎地理數(shù)據(jù)庫。柵格數(shù)據(jù)分為DMSP/OLS夜間燈光數(shù)據(jù)和MODIS數(shù)據(jù)。DMSP/OLS夜間燈光數(shù)據(jù)時間尺度為2001—2013年(https：//ngdc.noaa.gov/eog/dmsp)。MODIS數(shù)據(jù)為EVI數(shù)據(jù)(http：//www.gscloud.cn/)和MCD12Q1產(chǎn)品數(shù)據(jù)(https：//ladsweb.modaps.eosdis.nasa.gov)。

1.3 數(shù)據(jù)預處理

DMSP/OLS夜間燈光數(shù)據(jù)由于傳感器原因，導致長時間序列的各期數(shù)據(jù)之間缺乏連續(xù)性和可比性，且影像DN值存在飽和現(xiàn)象，因此，要對數(shù)據(jù)進行預處理，主要包括針對性的不連續(xù)性校正和飽和校正。在這之前，先根據(jù)研究區(qū)裁剪出陜西省夜間燈光影像，為盡可能減少面積變形，選取橢球體D_Krassovsky_1940，選取Albers投影，重采樣數(shù)據(jù)分辨率為1 000 m。

1)陜西省DMSP/OLS夜間燈光數(shù)據(jù)不連續(xù)性校正。由于存在同一年份有多個傳感器數(shù)據(jù)和不同年份傳感器不同的情況，且傳感器在獲取數(shù)據(jù)的過程中會受到多種因素的影響，因此，要對陜西省DMSP/OLS夜間燈光數(shù)據(jù)進行不連續(xù)性校正，以達到獲得隨著時間連續(xù)變化的夜間燈光數(shù)據(jù)的目的。參考文獻[16-17]的研究，不連續(xù)性校正具體可分為年內(nèi)融合和年際間校正兩步。

(1)年內(nèi)融合。需要進行年內(nèi)融合的年份為2001—2007年，年內(nèi)融合的計算如式(1)所示。

(1)

(2)年際間校正。在進行年內(nèi)融合處理后，年份間仍存在非連續(xù)性的情況，因此要進行年際間校正。其核心思想是同一像元位置后一年的DN值一定不小于前一年的DN值。需要進行年際間校正的年份為2001—2013年。年際間校正的計算如式(2)所示。

(2)

式中：n為夜間燈光數(shù)據(jù)的年份；i為像元；DN(n-1，i)、DN(n，i)、DN(n+1，i)分別為年際間校正后的i像元在第n-1、n、n+1年的DN值。

陜西省DMSP/OLS夜間燈光數(shù)據(jù)在不連續(xù)性校正前后累計DN值對比如圖1所示?？梢园l(fā)現(xiàn)，在進行不連續(xù)性校正處理前，同一年份內(nèi)數(shù)據(jù)上下差距較大，2001、2002、2007三年最為明顯，不同年份之間數(shù)據(jù)起伏明顯，而且在2009—2011年間存在明顯的陡增和陡降的趨勢，這些現(xiàn)象的出現(xiàn)都會極大地影響后續(xù)處理結(jié)果的準確性。在進行處理后，隨著時間的推移，數(shù)據(jù)呈連續(xù)性逐漸增大的趨勢，基本解決了不連續(xù)性的問題。

圖1 陜西省DMSP/OLS夜間燈光數(shù)據(jù)不連續(xù)性校正前后累計DN值

2)陜西省DMSP/OLS夜間燈光數(shù)據(jù)飽和校正。由于DMSP/OLS夜間燈光數(shù)據(jù)存在“天花板問題”，其取值范圍為0～63，這就使得城市中心、大型建筑場所等燈光聚集地燈光亮度達到飽和狀態(tài)無法繼續(xù)增加，燈光實際值將遠小于真實值，無法有效與經(jīng)濟數(shù)據(jù)相互結(jié)合研究，因此要進行飽和校正。考慮到植被與人類活動具有負相關的趨勢[18]，許多學者利用植被指數(shù)數(shù)據(jù)構建相關指數(shù)對夜間燈光數(shù)據(jù)進行飽和校正。植被指數(shù)是對地表植被狀況簡單、有效和經(jīng)驗的度量。歸一化植被指數(shù)(normalized vegetation index，NDVI)是最常用的植被指數(shù)，增強型植被指數(shù)(enhanced vegetation index，EVI)是在其基礎上改善而來的。與NDVI相比，EVI可以在一定程度上克服NDVI較容易飽和的缺陷，而且可以減弱土壤、大氣的影響，因此，本文選用EVI數(shù)據(jù)進行飽和校正。目前采用的方法主要有差值夜間燈光降飽和指數(shù)(difference nighttime-light vegetation index，DNVI)[19]、歸一化夜間燈光降飽和指數(shù)(normalized difference nighttime-light vegetation index，NDNVI)[20]、改進型夜間燈光降飽和指數(shù)(modified difference nighttime-light vegetation index，MDNVI)[21]等。分別用三種方法對數(shù)據(jù)進行處理，選擇效果最好的方法作為最終處理方法。三種夜間燈光降飽和指數(shù)的計算如式(3)至式(5)所示。

DNVI=NTL-EVI

(3)

(4)

(5)

式中：NTL為歸一化處理后的夜間燈光數(shù)據(jù)；EVI為增強型植被指數(shù)。

選取陜西省關中、陜北、陜南三個地區(qū)各一個市作為處理實驗區(qū)域，分別為西安市、延安市、漢中市，時間選擇2010年。處理前可以發(fā)現(xiàn)，原始DMSP/OLS夜間燈光數(shù)據(jù)中心區(qū)域亮度范圍比城市建成區(qū)要大，且飽和現(xiàn)象突出。用三種方法處理后與天地圖影像數(shù)據(jù)和原始DMSP/OLS夜間燈光數(shù)據(jù)作比較，處理結(jié)果如圖2所示。取值范圍在-1到1之間的DNVI和NDNVI方法突出了真實夜間燈光區(qū)域與溢出夜間燈光區(qū)域之間的差異，但相比取值范圍大于0的MDNVI方法在突出真實夜間燈光區(qū)域內(nèi)部差異有所不足。同時，MDNVI更強調(diào)地物對真實夜間燈光區(qū)域內(nèi)部差異的影響[22]。通過目視對比并與天地圖影像比較，發(fā)現(xiàn)MDNVI的校正效果最佳，因此，最終選用該方法。

圖2 陜西省DMSP/OLS夜間燈光數(shù)據(jù)飽和校正效果對比

3)夜間燈光數(shù)據(jù)與陜西省各縣全域和建設用地經(jīng)濟相關性。國外學者通過研究發(fā)現(xiàn)，DMSP/OLS夜間燈光數(shù)據(jù)的累計DN值與GDP數(shù)據(jù)存在高度相關關系，借助夜間燈光數(shù)據(jù)可以估算GDP數(shù)據(jù)。國內(nèi)學者也利用其進行了有關GDP數(shù)據(jù)反演、貧困度計算、經(jīng)濟區(qū)域差異等方面的研究。因此，將DMSP/OLS夜間燈光數(shù)據(jù)作為經(jīng)濟發(fā)展特征研究的替代量是合理和可行的。由于DMSP/OLS能夠捕捉到來自城市、鄉(xiāng)村、漁船等發(fā)出的微弱的燈光，但同時也包含了地面、湖面反射光芒，會對夜間燈光數(shù)據(jù)與經(jīng)濟的相關性研究產(chǎn)生干擾。為了盡可能減小這種干擾，考慮到經(jīng)濟的發(fā)展與人類活動息息相關，而基于人口主要分布在建設用地上這一認識，考慮提取研究陜西省各縣建設用地與夜間燈光相關性。以2013年為例，進行MCD12Q1數(shù)據(jù)預處理后，得到陜西省建設用地數(shù)據(jù)，利用ArcGIS軟件分別計算各縣全域和建設用地上夜間燈光數(shù)據(jù)DN 累計值，并與GDP相互擬合，結(jié)果如圖3所示?？梢园l(fā)現(xiàn)，建設用地上二者相關性高于各縣全域，說明該種方法對于減少無關燈光的干擾有一定的效果。

圖3 夜間燈光數(shù)據(jù)累計DN值與陜西省各縣全域和建設用地GDP相關性

2 研究方法

2.1 重心計算

經(jīng)濟的發(fā)展會促使經(jīng)濟中心的形成以及經(jīng)濟重心的移動，通過研究重心的位置和移動軌跡可以度量經(jīng)濟發(fā)展方向。重心模型[23]的計算如式(6)所示。

(6)

2.2 標準差橢圓

標準差橢圓是對點的方向和分布進行分析的一種算法。標準差橢圓的基本要素由長半軸、短半軸、扁率、中心點構成，可以用來表征數(shù)據(jù)分布的方向性、分布范圍、離散程度、數(shù)據(jù)變異程度等。其計算如式(7)所示[24]。

(7)

2.3 基尼系數(shù)

基尼系數(shù)是衡量區(qū)域收入差距的指標，其值域為0到1，數(shù)值越大表明區(qū)域經(jīng)濟發(fā)展越不平衡?，F(xiàn)構建陜西省夜間燈光數(shù)據(jù)的基尼模型，計算如式(8)所示。

(8)

式中：G表示陜西省夜間燈光數(shù)據(jù)的基尼系數(shù)；B表示洛倫茲曲線右下方部分面積。

2.4 變異系數(shù)

變異系數(shù)反映的是數(shù)據(jù)偏移平均值的程度，數(shù)值越大偏移程度越高。在基于陜西省夜間燈光數(shù)據(jù)計算過程中，變異系數(shù)的計算如式(9)所示。

(9)

2.5 局部莫蘭指數(shù)

空間相關性表明的是事物在空間上是否集聚的非隨機性，莫蘭指數(shù)是用來衡量空間相關性的一個重要指標，可分為全局莫蘭指數(shù)和局部莫蘭指數(shù)。由于本研究中夜間燈光數(shù)據(jù)的分布與局部莫蘭指數(shù)相比，注重屬性變量的區(qū)域差異的特性相符合，因此選用局部莫蘭指數(shù)進行計算。計算結(jié)果可分為高高聚集、高低聚集、低高聚集、低低聚集、不顯著五種，計算如式(10)、式(11)所示[25]。

(10)

(11)

3 結(jié)果與分析

3.1 陜西省經(jīng)濟發(fā)展方向特征分析

根據(jù)重心和標準差橢圓的計算方法，分別基于夜間燈光數(shù)據(jù)和行政區(qū)劃計算陜西省經(jīng)濟發(fā)展的重心和標準差橢圓，結(jié)果如圖4所示。對比基于夜間燈光和行政區(qū)劃數(shù)據(jù)計算的經(jīng)濟重心遷移和標準差橢圓變化情況，基于行政區(qū)劃數(shù)據(jù)計算的經(jīng)濟重心遷移與基于夜間燈光的相比規(guī)律性較弱，且基于行政區(qū)劃數(shù)據(jù)計算的標準差橢圓變化規(guī)律不明顯。經(jīng)分析主要原因是：首先，數(shù)據(jù)樣本較少，陜西省共107個縣級行政區(qū)所有樣本僅107例，而夜間燈光是1 km柵格數(shù)據(jù)樣本，高達205 555例；其次，以行政區(qū)劃為單元基于GDP計算時，將GDP視為在行政區(qū)劃內(nèi)均勻分布，但實際上GDP與土地利用類型關系較大，而且主要分布在建設用地上。夜間燈光數(shù)據(jù)與傳統(tǒng)數(shù)據(jù)相比，數(shù)據(jù)樣本大，考慮了空間異質(zhì)性，克服了上述不足，因此，選用基于夜間燈光數(shù)據(jù)計算的重心遷移和標準差橢圓進行分析。

注：該圖基于陜西測繪地理信息局標準地圖服務下載的審圖號為陜S(2018)006號的標準地圖制作，底圖無修改。圖4 基于夜間燈光和行政區(qū)劃計算的陜西省經(jīng)濟重心和標準差橢圓

從重心遷移的角度來看，2001—2013年間陜西省經(jīng)濟重心大致沿自南往北方向移動，水平位移較小，垂直位移較大，直線位移距離44.15 km，說明隨著時間的推移，陜西北部經(jīng)濟發(fā)展狀況優(yōu)于南部。陜西南北跨度極大，可分為陜北、陜南、關中三大塊。一個地區(qū)的經(jīng)濟發(fā)展要受到自然條件、地理位置、資源稟賦、政策傾斜等諸多因素的影響。由于陜北地區(qū)農(nóng)村經(jīng)濟繁榮發(fā)展、礦產(chǎn)資源的發(fā)現(xiàn)和大規(guī)模開發(fā)、交通網(wǎng)絡體系基本建成等原因，逐漸拉開了與交通不便、自然災害頻繁的陜南地區(qū)的差距。

同時，經(jīng)濟重心年間移動距離大致呈現(xiàn)先增加后減少的趨勢，具體來看，在2001—2007年間移動距離較大，在2007—2013年間移動距離較小。經(jīng)分析，這主要與陜北陜南地區(qū)自然資源、政策傾斜和基礎設施建設三方面因素有關。2007年之前，陜北一直在這三方面占有較大優(yōu)勢，與陜南經(jīng)濟發(fā)展差距較大；2007年之后，陜南由于發(fā)展動植物資源產(chǎn)業(yè)，“精準扶貧”政策實施，交通體系開始逐步完善，經(jīng)濟取得了一定的發(fā)展，使得與陜北的相對差距縮小。

從標準差橢圓變化的角度來看，相較于基于行政區(qū)劃數(shù)據(jù)計算的標準差橢圓，基于夜間燈光數(shù)據(jù)計算的標準差橢圓更加狹長，詳細數(shù)據(jù)如表1所示。長軸長度從2001年的509 km逐漸增加到2013年的579 km，為南北方向分布；短軸長度從2001年的209 km逐漸增加到2013年的223 km。相應地，隨著長軸和短軸長度的增加，橢圓面積也從2001年的83 380 km2增加到了2013年的131 044 km，標準差橢圓面積的增大反映了隨著時間的增加，陜西省經(jīng)濟不斷取得進度與發(fā)展，因此使得陜西省夜間燈光數(shù)據(jù)的廣度和亮度逐漸加大。方向角與扁率的變化正好相反，方向角逐漸變小，扁率逐漸變大，方向角變小說明了橢圓的方向逐漸由東南向西北方向偏轉(zhuǎn)，這與圖4中基于夜間燈光計算的經(jīng)濟重心偏移方向一致；長軸與短軸之比即扁率從2.44增長至2.60，表明陜西省經(jīng)濟格局在南北方向上有擴張趨勢，在東西方向有收縮趨勢，這也進一步反映了陜西省南北經(jīng)濟發(fā)展不平衡。

表1 基于夜間燈光計算的陜西省經(jīng)濟標準差橢圓變化

3.2 陜西省經(jīng)濟發(fā)展差異性特征分析

以陜西省10個市級行政區(qū)為基本單位，將這10個市級行政區(qū)的平均夜間燈光數(shù)據(jù)按從小到大順序排列，分別對其進行疊加累計，根據(jù)洛倫茲曲線求得2001—2013年陜西省夜間燈光數(shù)據(jù)的基尼系數(shù)，同理求得變異系數(shù)，具體結(jié)果如圖5所示?？梢园l(fā)現(xiàn)，2001—2013年間，這一時間段內(nèi)基尼系數(shù)都處于0.3～0.4的區(qū)間范圍內(nèi)，說明陜西省經(jīng)濟收入分配相對合理。陜西省基尼系數(shù)和變異系數(shù)的變化趨勢基本一致，從整體上來看，隨著時間的推移，都呈波動性減小的趨勢，說明陜西省經(jīng)濟呈現(xiàn)向平衡的趨勢發(fā)展。具體來看，在2008年前后，陜西省經(jīng)濟發(fā)展平衡性有了較大改變，由逐漸不平衡轉(zhuǎn)變?yōu)槠胶獍l(fā)展，這證明了在2008年之后，由夜間燈光數(shù)據(jù)模擬的陜西省城市間的區(qū)域經(jīng)濟差異開始越來越小，經(jīng)濟發(fā)展越來越平衡。經(jīng)分析，這主要得益于國家實施的區(qū)域協(xié)調(diào)發(fā)展戰(zhàn)略和關中-天水經(jīng)濟區(qū)的成立，共6市1區(qū)，形成了中國西部發(fā)達的城市群和產(chǎn)業(yè)集群帶，這一地區(qū)經(jīng)濟的迅速發(fā)展，使得陜西省經(jīng)濟相對差異開始逐漸減小。

圖5 2001—2013年陜西省經(jīng)濟發(fā)展時間平衡性變化

3.3 陜西省經(jīng)濟發(fā)展空間集聚性特征分析

通過夜間燈光累計DN值計算2001—2013年陜西省縣域經(jīng)濟發(fā)展空間集聚性變化，可以避免傳統(tǒng)統(tǒng)計年鑒數(shù)據(jù)存在的數(shù)據(jù)缺失、GDP均勻分布等缺陷，計算結(jié)果如表2所示?？偟貋砜矗吒呔奂?、高低聚集、低高聚集、低低聚集的縣域個數(shù)都呈逐漸減少或不變的趨勢。具體來看，高高聚集的縣域主要分布在關中中部和陜北的部分地區(qū)，包括府谷縣、安塞縣、灞橋區(qū)等。府谷縣、安塞縣一直屬于陜西省經(jīng)濟10強縣，灞橋區(qū)隸屬于西安市，這些區(qū)域與經(jīng)濟發(fā)展條件較好的縣域基本吻合，構成了陜西省經(jīng)濟重要的發(fā)展軸。2001—2013年間陜北地區(qū)的高高聚集的縣域表現(xiàn)出明顯的輻射帶動作用，極大地帶動了周圍地區(qū)的經(jīng)濟發(fā)展，從2001年的2個增長到了2013年的7個，這與上文研究中陜西省經(jīng)濟重心逐漸向北轉(zhuǎn)移相符合。高低聚集分布著的主要為經(jīng)濟發(fā)展?jié)摿^大、速度較快的縣域，包括綏德縣、彬縣、洋縣等地，這些縣域由于自然資源人才等因素優(yōu)于周圍縣域，更易于吸引資本發(fā)展相關產(chǎn)業(yè)，經(jīng)濟發(fā)展也就更為迅速。以2013年綏德縣為例，GDP總量為53.94億元，GDP增速為7.4%，都要比周圍清澗縣、吳堡縣要高出不少。低高聚集主要分布在陜北地區(qū)，與周圍縣域相比，經(jīng)濟發(fā)展速度較慢，處于相對衰退趨勢。2001—2013年間陜北地區(qū)低高集聚縣域數(shù)量從6個下降為2個，但位于延安市的集中分布態(tài)勢并未發(fā)生改變。低低聚集縣域一直只有1個，分布在陜南地區(qū)安康市的漢濱區(qū)，直到2020年，經(jīng)陜西省政府同意，漢濱區(qū)才正式退出貧困縣序列。

表2 2001—2013年陜西省經(jīng)濟發(fā)展空間集聚性變化

4 結(jié)束語

本文以2001—2013年陜西省縣域GDP數(shù)據(jù)、夜間燈光數(shù)據(jù)、MODIS數(shù)據(jù)等多源數(shù)據(jù)作為數(shù)據(jù)源，在對夜間燈光進行預處理后，采用重心、標準差橢圓、基尼系數(shù)、變異系數(shù)、局部莫蘭指數(shù)的計算方法，從經(jīng)濟發(fā)展方向、經(jīng)濟發(fā)展差異性以及空間集聚性三個方面對陜西省2001—2013經(jīng)濟發(fā)展特征進行了分析比較，主要結(jié)論如下。

1)2001—2013年，陜西省北部經(jīng)濟總量比南部大、南部增長速度比北部快，陜西省經(jīng)濟總量增大的同時內(nèi)部差異縮小。具體來看，2001—2013年間，陜西省經(jīng)濟重心大致沿自南往北方向移動，且移動距離為先增加后減少的趨勢，這表明陜西省經(jīng)濟發(fā)展總量北部優(yōu)于南部，增長速度南部優(yōu)于北部；陜西省經(jīng)濟標準差橢圓面積和扁率都逐漸變大，這表明陜西省經(jīng)濟總量在不斷增加的同時內(nèi)部差異逐漸減小。

2)2001—2013年，陜西省經(jīng)濟發(fā)展時間平衡性為波動性減小的趨勢，收入分配相對合理。基于夜間燈光數(shù)據(jù)計算的基尼系數(shù)和變異系數(shù)變化趨勢基本一致，2008年是陜西省內(nèi)部區(qū)域經(jīng)濟發(fā)展由逐漸不平衡轉(zhuǎn)向平衡的拐點。

3)2001—2013年，高高聚集和高低聚集為最常見的兩種經(jīng)濟聚集類型。高高聚集主要分布在關中和陜北的發(fā)達縣域，陜北地區(qū)高高聚集縣域輻射帶動作用明顯；高低聚集為經(jīng)濟發(fā)展?jié)摿^大、速度較快縣域；低高聚集主要分布在陜北的不發(fā)達縣域；低低聚集僅分布在陜南的較為貧困的漢濱區(qū)。

本研究還有兩個方向可以做進一步研究。第一，由于DMSP/OLS夜間燈光數(shù)據(jù)2013年已停止更新，因此，如果結(jié)合后續(xù)發(fā)布的NPP/VIIRS數(shù)據(jù)，將兩種數(shù)據(jù)進行整合后，可以將研究時間延長到2019年，具有更大的現(xiàn)實意義；第二，不同的空間鄰接矩陣可能會對空間集聚性的結(jié)果產(chǎn)生較大的影響，在與結(jié)果進行精度驗證和比較后，可以探究不同的空間鄰接矩陣哪種效果最優(yōu)。