張瓊藝， 李昆， 雍志瑋， 熊俊楠,， 程維明， 肖坤洪， 劉東麗

(1.西南石油大學(xué)土木工程與測(cè)繪學(xué)院，成都 610500； 2.四川電力設(shè)計(jì)咨詢有限責(zé)任公司，成都 610041； 3.西南石油大學(xué)地球科學(xué)與技術(shù)學(xué)院，成都 610500； 4.中國(guó)科學(xué)院地理科學(xué)與資源研究所資源與環(huán)境信息系統(tǒng)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100101； 5.四川省煤田測(cè)繪工程院，成都 610072； 6.自然資源部第六測(cè)量隊(duì)，成都 610500)

0 引言

貧困是一個(gè)長(zhǎng)期的世界性難題，是造成社會(huì)不穩(wěn)定的主要因素，也是造成人類潛力喪失的最主要原因之一[1]。中國(guó)作為世界上最大的發(fā)展中國(guó)家，具有貧困人口多、分布廣的基本格局[2]。2020年11月23日，貴州省宣布其所有貧困縣脫貧摘帽，也標(biāo)志著我國(guó)消滅了絕對(duì)貧困和區(qū)域性整體貧困。2020年后，中國(guó)扶貧工作轉(zhuǎn)變?yōu)榫徑獍l(fā)展不平衡、不充分的相對(duì)貧困。相對(duì)貧困的主要特征為區(qū)域與城鄉(xiāng)收入差異、社會(huì)公共服務(wù)獲取不平等[3]。而要解決相對(duì)貧困問(wèn)題，需要精準(zhǔn)幫扶并對(duì)癥下藥，而精準(zhǔn)識(shí)別是精準(zhǔn)幫扶的前提條件[4]。因此，目前各個(gè)貧困地區(qū)的貧困測(cè)度和發(fā)展?fàn)顩r仍然是黨和政府關(guān)注的重點(diǎn)[5]，對(duì)相對(duì)貧困地區(qū)進(jìn)行精確測(cè)度和評(píng)估仍是未來(lái)很長(zhǎng)一段時(shí)間的需要[6]。

目前，貧困概念逐漸由經(jīng)濟(jì)維度擴(kuò)展到社會(huì)、文化、教育、自然等多維度層面。多維貧困測(cè)度可以綜合測(cè)量地區(qū)貧困程度及類型。1997年，聯(lián)合國(guó)開(kāi)發(fā)計(jì)劃署(united nations development programme, UNDP)結(jié)合壽命、讀寫(xiě)能力、生活水平、可支配收入水平和失業(yè)率等維度開(kāi)發(fā)了分別適合發(fā)展中國(guó)家和發(fā)達(dá)國(guó)家的人類貧困指數(shù)，來(lái)對(duì)國(guó)家的平均發(fā)展水平進(jìn)行度量[7]。2010年，UNDP與牛津大學(xué)貧困與人類發(fā)展中心共同構(gòu)建開(kāi)發(fā)了多維貧困指數(shù)(multidimensional poverty index, MPI)，該指數(shù)從健康、教育和生活水平3個(gè)維度的10個(gè)指標(biāo)反映個(gè)體的貧困度，也可以反映地區(qū)人文發(fā)展情況，更加符合實(shí)際需求[8]。國(guó)內(nèi)外學(xué)者也提出了多種多維貧困指數(shù)構(gòu)建方法，并運(yùn)用到多維貧困測(cè)度中。主要的有模糊集方法[9]、公理化方法[10]和信息論方法[11]等。Alkire等[12]提出了“A-F法”，因其更好的具體性與細(xì)致性，也成為了當(dāng)前應(yīng)用最廣泛的多維貧困測(cè)度模型。但是多維度也意味著需要更多的數(shù)據(jù)。但是目前MPI的構(gòu)建主要依靠統(tǒng)計(jì)或調(diào)查數(shù)據(jù)，此類型數(shù)據(jù)存在獲取不易、成本較高、時(shí)效不佳等不足。因此，有必要提供新型數(shù)據(jù)源來(lái)識(shí)別和評(píng)估中國(guó)的貧困狀況，以滿足貧困動(dòng)態(tài)監(jiān)測(cè)的需求。

與傳統(tǒng)的統(tǒng)計(jì)調(diào)查數(shù)據(jù)相比，遙感數(shù)據(jù)可以提供高效、準(zhǔn)確的空間數(shù)據(jù)，以便從多尺度角度觀察社會(huì)經(jīng)濟(jì)和自然現(xiàn)象。而夜光遙感數(shù)據(jù)在區(qū)域貧困研究分析中也顯示出了良好的潛力[13]。夜光實(shí)質(zhì)上是對(duì)電氣化水平的測(cè)量，而電氣化不僅使日間時(shí)間得到了擴(kuò)展，還改善了獲得醫(yī)療保健、食物、住所、教育等的方式[14]，因此夜光也被證明能夠反映區(qū)域的社會(huì)經(jīng)濟(jì)發(fā)展?fàn)顩r。已有研究通過(guò)建立夜光數(shù)據(jù)與替代表征貧困的貧困指數(shù)或貧困概率之間的關(guān)系來(lái)對(duì)貧困進(jìn)行識(shí)別測(cè)度[15-17]。目前基于夜光數(shù)據(jù)的貧困測(cè)度研究偏少，而已有研究使用傳統(tǒng)的建模方法取得了較為理想的成果。但研究中仍存在一些其他的問(wèn)題： ①多數(shù)研究使用單一的夜光指數(shù)，可能忽略了夜光中的部分有用信息； ②夜光與貧困的非線性關(guān)系可能無(wú)法由傳統(tǒng)的經(jīng)驗(yàn)?zāi)Ｐ瓦M(jìn)行較好的表達(dá)； ③貧困是一個(gè)長(zhǎng)期、動(dòng)態(tài)的過(guò)程，而多數(shù)研究?jī)H使用DMSP-OLS燈光數(shù)據(jù)或者僅使用NPP-VIIRS燈光數(shù)據(jù)開(kāi)展研究，或者在同時(shí)使用2種數(shù)據(jù)時(shí)未進(jìn)行一致性校準(zhǔn)，這會(huì)導(dǎo)致無(wú)法實(shí)現(xiàn)長(zhǎng)期連續(xù)動(dòng)態(tài)的貧困監(jiān)測(cè)及分析。因此，可以對(duì)DMSP-OLS與NPP-VIIRS燈光進(jìn)行校準(zhǔn)整合，擴(kuò)展貧困測(cè)度的時(shí)間序列，并提取夜光的數(shù)據(jù)特征，使用機(jī)器學(xué)習(xí)模型構(gòu)建基于夜間燈光特征的貧困估算模型，實(shí)現(xiàn)長(zhǎng)時(shí)間序列的貧困時(shí)空動(dòng)態(tài)監(jiān)測(cè)與分析。

夜間燈光被證明與人類社會(huì)經(jīng)濟(jì)活動(dòng)高度相關(guān)，但是沒(méi)有燈光也不代表沒(méi)有人類活動(dòng)，特別是在人口稀疏的地區(qū)。而中國(guó)西南4省(市)的人口密度大，但是社會(huì)經(jīng)濟(jì)發(fā)展不均衡，半數(shù)以上的縣被確定為國(guó)家級(jí)貧困縣。因此，本文通過(guò)綜合西南4省(市)縣域內(nèi)的夜光信息，更為準(zhǔn)確地反映西南4省(市)縣域整體的發(fā)展情況。首先，本研究應(yīng)用粒子群優(yōu)化算法(particle swarm optimization,PSO)-反向傳播(back propagation, BP)神經(jīng)網(wǎng)絡(luò)模型構(gòu)建夜間燈光校準(zhǔn)整合模型，并以此實(shí)現(xiàn)了DMSP-OLS和NPP-VIIRS這2種燈光數(shù)據(jù)的整合，擴(kuò)展了數(shù)據(jù)的覆蓋周期； 從長(zhǎng)時(shí)間序列夜光數(shù)據(jù)的集中趨勢(shì)、分散程度、分布特征和空間特征中提取了12種燈光指數(shù)，構(gòu)建了基于夜間燈光的貧困估算模型，并基于估算模型得到西南4省(市)各縣域近20 a的MPI指數(shù)估算值，進(jìn)而對(duì)貧困時(shí)空特征進(jìn)行分析，以期為區(qū)域后續(xù)減貧工作的開(kāi)展以及扶貧政策的制定實(shí)施提供一定的理論基礎(chǔ)和參考。

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

1.1 研究區(qū)概況

研究區(qū)為中國(guó)西南地區(qū)，包含四川省、云南省、貴州省和重慶市4省(市)，共438個(gè)縣級(jí)行政區(qū)。研究區(qū)的地理位置在N21°08′～34°18′ ，E97°20′～110°11′ 之間，占地約為112.78×104km2(圖1)。西南地區(qū)地形可分為3個(gè)主要地形單元： 四川盆地和東北部的丘陵地區(qū)、東南部的云南—貴州高原及西北部的高海拔青藏高原。由于地質(zhì)、地形、氣候和水熱條件等差異較大，該地區(qū)的農(nóng)業(yè)條件具有明顯差異。而且該地區(qū)也是中國(guó)重要的糧食生產(chǎn)區(qū)[18]。

圖1 研究區(qū)位置示意圖Fig.1 Location of the study area

據(jù)統(tǒng)計(jì)數(shù)據(jù)可知，2019年西南地區(qū)總?cè)丝跒?.998億人，占全國(guó)總?cè)丝诩s14.2%，人口密度約為177人/km2，人口密度較大但分布不均(圖2(a))。人口密度較大的區(qū)域主要分布在四川盆地、云南滇中盆地等地勢(shì)平坦、海拔較低的位置。2019年全國(guó)人均國(guó)內(nèi)生產(chǎn)總值(gross domestic product，GDP)為70 892元，居民人均可支配收入為30 733元，而西南地區(qū)中僅有重慶市的人均GDP高于全國(guó)平均水平，其余各省都遠(yuǎn)低于平均線。同時(shí)，4省(市)間的經(jīng)濟(jì)發(fā)展并不平衡，例如2019年四川省GDP總量為46 615.82億元，占到了西南地區(qū)的42.29%。同時(shí)，各省(市)內(nèi)部各城市GDP情況分布同樣不均衡(圖2(b))。

(a) 人口分布 (b) GDP分布圖2 2019年研究區(qū)人口和GDP分布Fig.2 Population and GDP distribution of the study area in 2019

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

在本研究中，2種夜光產(chǎn)品(DMSP-OLS和NPP-VIIRS NTL)被用作主要的夜光數(shù)據(jù)源。1992—2013年的第四版DMSP-OLS年際穩(wěn)定夜光合成數(shù)據(jù)由科羅拉多礦業(yè)大學(xué)地球觀測(cè)組提供(https: //eogdata.mines.edu/products/dmsp/)。2012年以來(lái)的NPP-VIIRS夜光數(shù)據(jù)也是從科羅拉多礦業(yè)大學(xué)地球觀測(cè)組獲得的(https: //eogdata.mines.edu/products/vnl/)，此外，本研究中使用的數(shù)據(jù)源還包括行政邊界數(shù)據(jù)、社會(huì)經(jīng)濟(jì)統(tǒng)計(jì)數(shù)據(jù)和數(shù)字高程模型(digital elevation model，DEM)數(shù)據(jù)。各個(gè)數(shù)據(jù)源的詳細(xì)信息見(jiàn)表1。

表1 本文主要數(shù)據(jù)匯總Tab.1 Description of each data sources used in this study

2 數(shù)據(jù)處理及分析方法

2.1 長(zhǎng)時(shí)間序列夜光數(shù)據(jù)集的構(gòu)建

由于DMSP-OLS 數(shù)據(jù)和 NPP-VIIRS 數(shù)據(jù)之間存在傳感器設(shè)計(jì)和空間分辨率的差異，所以建立2000—2019年的長(zhǎng)時(shí)間序列夜間燈光數(shù)據(jù)需要對(duì)這2種數(shù)據(jù)進(jìn)行一致性校準(zhǔn)預(yù)處理； 首先，通過(guò)投影、重采樣進(jìn)行預(yù)處理。然后，采用不變目標(biāo)區(qū)域法、年內(nèi)校正和年際校正，以減少亮度飽和效應(yīng)，提高DMSP-OLS時(shí)間序列的一致性和可比性[19-21]。此外，本研究采取了4個(gè)步驟來(lái)減少NPP-VIIRS每月數(shù)據(jù)中的負(fù)面因素[22]： ①應(yīng)用從校正的DMSP-OLS數(shù)據(jù)中提取的掩碼來(lái)去除NPP-VIIRS數(shù)據(jù)的背景噪聲，同時(shí)將負(fù)值替換為零； ②將北京、上海、廣州和深圳的最大像素值作為閾值，以過(guò)濾每個(gè)月夜光圖像中的突變亮度，并將成都、重慶、昆明和貴陽(yáng)的像素值作為補(bǔ)充，以處理一些未解決的像素； ③為了消除雜散光污染的影響，采用除4—8月外的NPP-VIIRS 每月夜光數(shù)據(jù)來(lái)合成年度數(shù)據(jù)[23-24]； ④采用與DMSP-OLS數(shù)據(jù)相同的方法對(duì)NPP-VIIRS夜光年度數(shù)據(jù)進(jìn)行剪裁、重采樣和重新投影。

通過(guò)數(shù)據(jù)預(yù)處理，可以獲得2000—2013年的DMSP-OLS校正數(shù)據(jù)和2012—2019年的NPP-VIIRS校正數(shù)據(jù)。隨后本研究通過(guò)實(shí)施從NPP-VIIRS數(shù)據(jù)到DMSP-OLS數(shù)據(jù)的像素級(jí)轉(zhuǎn)換來(lái)拓展夜光數(shù)據(jù)。2個(gè)數(shù)據(jù)集的重疊期(2012年和2013年)為2種夜光數(shù)據(jù)的整合提供了基礎(chǔ)。目前，偽不變特征(Pseudo invariant features, PIF)已被證明在多時(shí)空衛(wèi)星圖像的相互校準(zhǔn)中非常有用[25-26]。因此，考慮到世界上沒(méi)有專門(mén)為夜光數(shù)據(jù)設(shè)計(jì)的校準(zhǔn)點(diǎn)，在本研究中可以以PIF的形式提供用于相互校準(zhǔn)的樣品[27]。根據(jù)以往的研究[27-28]，以20%為閾值的變異系數(shù)(coefficient of variation, CV)來(lái)選擇合適的PIF，其過(guò)程包括3個(gè)方面： ①用3×3窗口計(jì)算DMSP-OLS柵格和NPP-VIIRS柵格中所有像素的CV值； ②提取DMSP-OLS和NPP-VIIRS中CV值小于20%的交叉區(qū)域，這些區(qū)域被認(rèn)為在時(shí)間和空間上是穩(wěn)定的； ③將交叉區(qū)域的正像素組成PIFs。鑒于2012年夜光數(shù)據(jù)的不完整性，本研究采用了2013年的PIF。CV的定義為：

(1)

式中：σ為3×3窗口中像素值的標(biāo)準(zhǔn)偏差；μ為3×3窗口中像素值的平均值。

在本研究中，通過(guò)將PSO與BP人工神經(jīng)網(wǎng)絡(luò)相結(jié)合，采用PSO-BP混合算法來(lái)構(gòu)建人工神經(jīng)網(wǎng)絡(luò)，以探索DMSP-OLS數(shù)據(jù)與NPP-VIIRS數(shù)據(jù)之間的潛在關(guān)系，構(gòu)建長(zhǎng)時(shí)間序列的夜間燈光數(shù)據(jù)集。

2.2 MPI的構(gòu)建

貧困是一種多維現(xiàn)象，涉及人類發(fā)展和自然環(huán)境的諸多方面?；诮y(tǒng)計(jì)的多維貧困評(píng)估已被學(xué)者廣泛采用[12, 29-31]?；谘芯繀^(qū)的數(shù)據(jù)可用性，本研究采用MPI作為縣域?qū)用娑嗑S貧困的指標(biāo)。參考以往研究[16, 32-33]提出的可持續(xù)生計(jì)框架，以科學(xué)性、全面性、可度量性和可操作性為原則，從經(jīng)濟(jì)維度、社會(huì)維度和自然維度提取了11個(gè)社會(huì)經(jīng)濟(jì)和自然變量，建立MPI作為探索使用夜光數(shù)據(jù)進(jìn)行貧困評(píng)估的參考。此外，MPI的計(jì)算是一個(gè)多指標(biāo)綜合評(píng)價(jià)的過(guò)程，而各指標(biāo)的權(quán)重對(duì)于最終的MPI起著至關(guān)重要的決定性作用。熵值法是一種通過(guò)判斷信息熵的大小來(lái)進(jìn)行權(quán)重賦值的方法[34]。在信息論中認(rèn)為信息量的大小與不確定性和熵值成反比，可以通過(guò)熵值來(lái)判斷指標(biāo)的離散程度，而指標(biāo)的離散程度通常被認(rèn)為與其影響力成正比。因此熵值法僅通過(guò)指標(biāo)自身的特征來(lái)確定權(quán)重，在一定程度上避免了主觀因素的干擾，所得到的結(jié)果也更加客觀。11個(gè)指標(biāo)及其權(quán)重如表2所示。

表2 多維貧困測(cè)度指標(biāo)體系及權(quán)重Tab.2 Evaluation indices and weight distribution of multidimensional poverty

對(duì)于數(shù)據(jù)的標(biāo)準(zhǔn)化處理，多維貧困指標(biāo)體系中各指標(biāo)原始數(shù)據(jù)的量綱和數(shù)量級(jí)不盡相同，因此在數(shù)據(jù)分析前需要進(jìn)行數(shù)據(jù)標(biāo)準(zhǔn)化處理以消除量綱不同而造成的影響。由于涉及到正負(fù)2種屬性的指標(biāo)，其中正屬性指標(biāo)的數(shù)據(jù)值大小與其對(duì)綜合指數(shù)的“促進(jìn)”作用成正比，而負(fù)屬性數(shù)據(jù)值與“抑制”作用成正比，所以需要采取不同方式對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理。假設(shè)共計(jì)k個(gè)年份，m個(gè)縣域和n項(xiàng)指標(biāo)，則對(duì)于正負(fù)屬性指標(biāo)的標(biāo)準(zhǔn)化方式分別為：

(2)

(3)

然后計(jì)算樣本縣的區(qū)域MPI，其計(jì)算公式為：

(4)

式中：MPIi為第i個(gè)縣在某一年的數(shù)值；wj為第j個(gè)指數(shù)的權(quán)重；Xij為第i個(gè)縣在同一年的第j個(gè)指數(shù)的標(biāo)準(zhǔn)值。本研究采用熵值法來(lái)確定指數(shù)的權(quán)重(表2)，相應(yīng)的過(guò)程與前人的工作一致[35]。因此，獲得了38個(gè)樣本縣2000—2019年的MPI。理論上，MPI指數(shù)值越大代表對(duì)應(yīng)縣域經(jīng)濟(jì)發(fā)展水平高，居民生活水平和生產(chǎn)條件高，貧困度低； 而MPI指數(shù)值越小代表社會(huì)經(jīng)濟(jì)發(fā)展速度緩慢，生活和生產(chǎn)條件落后，貧困度高。因此，MPI指數(shù)可以用于度量縣域綜合貧困度。

2.3 基于夜光的貧困測(cè)度方法

本研究的基本假設(shè)是夜光信息可以作為分析區(qū)域貧困的準(zhǔn)確和有利的識(shí)別特征，這已經(jīng)被現(xiàn)有研究的發(fā)現(xiàn)所證明[33]。并通過(guò)挖掘MPI和夜光特征變量之間的相關(guān)性建立了相應(yīng)的評(píng)價(jià)模型，以評(píng)價(jià)研究區(qū)域的多維度貧困。

首先，參考相關(guān)研究的方法[17, 32, 36]，本研究選取了12個(gè)夜光特征，其中不僅包括平均亮度指數(shù)等簡(jiǎn)單統(tǒng)計(jì)特征，還包括空間特征。并且可以通過(guò)這些選定的指數(shù)來(lái)描述各縣夜光分布的集中趨勢(shì)、分散程度、分布特征和空間特征，從而更全面地反映各縣之間在不同方面的差異和相似之處。表3描述了12個(gè)夜光特征變量用于評(píng)價(jià)貧困。

表3 區(qū)域夜光特征變量及描述Tab.3 Feature variables of regional NTL and its description

考慮到貧困的復(fù)雜性以及缺乏連接貧困和夜光特征的先驗(yàn)知識(shí)，還使用了PSO-BP算法來(lái)擬合和訓(xùn)練模型。具體來(lái)說(shuō)，根據(jù)樣本地區(qū)的數(shù)據(jù)，選擇了12個(gè)夜光特征變量和年份作為輸入?yún)?shù)，而縣級(jí)MPI作為輸出參數(shù)。最后將基于PSO-BP神經(jīng)網(wǎng)絡(luò)的貧困估算模型推廣到西南4省(市)，可以得到2000—2019年研究區(qū)中各縣區(qū)的多維貧困估算值。因此，將整理得到的西南地區(qū)2000—2019年435個(gè)縣域的12種燈光指數(shù)，結(jié)合時(shí)間因子作為已訓(xùn)練好的PSO-BP神經(jīng)網(wǎng)絡(luò)模型的輸入樣本，共計(jì)8 700組。經(jīng)由模型輸出得到逐年各縣域的夜光估算MPI指數(shù)(NLMPI)。

2.4 時(shí)空演變分析方法

為了調(diào)查和量化貧困縣的時(shí)空動(dòng)態(tài)，計(jì)算了Global Moran’s I和局部空間關(guān)聯(lián)指標(biāo)(local indicators of spatial association, LISA)。Global Moran’s I評(píng)估了整體的空間相關(guān)現(xiàn)象，范圍從[-1，1][17, 32]。此外，正值表示正空間相關(guān)性，其絕對(duì)值越接近1，空間相關(guān)性越強(qiáng)； 0值表示隨機(jī)分布。其計(jì)算公式為：

(5)

(6)

式中LISAi為地理現(xiàn)象內(nèi)部單元i的局部自相關(guān)指數(shù)。

3 結(jié)果與分析

3.1 夜間燈光指數(shù)與MPI的關(guān)聯(lián)性分析

構(gòu)建MPI時(shí)用到的統(tǒng)計(jì)指標(biāo)通常較多，由于各地區(qū)的統(tǒng)計(jì)普查方式和指標(biāo)項(xiàng)目不盡相同，且各項(xiàng)統(tǒng)計(jì)數(shù)據(jù)的獲取時(shí)間節(jié)點(diǎn)也無(wú)法完全保持一致，導(dǎo)致全區(qū)域性的數(shù)據(jù)獲取不易。一些學(xué)者則通過(guò)選取典型研究區(qū)來(lái)進(jìn)行建模，然后進(jìn)行推廣。所以，本研究同樣通過(guò)選取典型研究區(qū)來(lái)分析多維貧困指標(biāo)體系的合理性，并挖掘夜間燈光指數(shù)與MPI間的關(guān)系。圖3展示了2019年重慶市各縣域MPI指數(shù)分布狀況。縣域多維貧困的空間分布模式基本滿足“一圈兩翼”的模式： MPI值較高的區(qū)域主要分布在重慶市的“一小時(shí)經(jīng)濟(jì)圈”區(qū)域； MPI值較低的區(qū)域主要位于重慶市的東北翼和東南翼?；谧匀粩帱c(diǎn)法，將重慶市38個(gè)縣域的MPI指數(shù)值分為了5類： ①極低[1.489 4,1.925 7]; ②較低(1.925 7,2.554 1]; ③中等(2.554 1,3.116 8]; ④較高(3.116 8,4.550 1]; ⑤極高M(jìn)PI(4.550 1,7.025 5]。

圖3 2019年重慶市MPI指數(shù)分布Fig.3 Spatial patterns of MPI in Chongqing in 2019

從圖3中可知，渝中區(qū)、江北區(qū)、九龍坡區(qū)、南岸區(qū)和沙坪壩區(qū)等縣區(qū)的MPI指數(shù)最高，同時(shí)作為重慶市的文化和經(jīng)濟(jì)中心，具有良好的社會(huì)經(jīng)濟(jì)基礎(chǔ)和發(fā)展環(huán)境，也是重慶市經(jīng)濟(jì)最繁華的縣區(qū)，使得縣區(qū)貧困度普遍較低。而城口縣、巫溪縣、巫山縣和酉陽(yáng)土家族苗族自治縣等MPI值最低。這些縣區(qū)由于自然環(huán)境惡劣、交通條件不便和少數(shù)民族聚集等原因，導(dǎo)致經(jīng)濟(jì)基礎(chǔ)薄弱，而且其社會(huì)經(jīng)濟(jì)發(fā)展主要依賴于第一產(chǎn)業(yè)，同時(shí)第二和第三產(chǎn)業(yè)的發(fā)展速度緩慢，使得區(qū)域發(fā)展受到限制，貧困程度加深。以上結(jié)果與實(shí)際情況基本相符，也說(shuō)明了MPI指數(shù)可以較好地反映縣區(qū)的貧困狀況。

2000—2019年重慶市各縣區(qū)的MPI指數(shù)如圖4所示，可以看出2019年重慶市各縣區(qū)的MPI指數(shù)對(duì)于2000年的MPI指數(shù)有明顯提高，說(shuō)明近20 a來(lái)重慶市各縣區(qū)得到了較大的發(fā)展。不難發(fā)現(xiàn)，渝中區(qū)、江北區(qū)和九龍波區(qū)等綜合經(jīng)濟(jì)發(fā)展水平較高，而城口縣、巫溪縣、酉陽(yáng)土家族苗族自治縣和巫山縣等經(jīng)濟(jì)發(fā)展水平較低，且都是重慶市的貧困縣。

(a) 2000年 (b) 2005年

(c) 2010年 (d) 2015年

(e) 2019年圖4 重慶市各縣區(qū)MPI指數(shù)Fig.4 MPI of various county in Chongqing

基于本文提出的夜間燈光指數(shù)的構(gòu)建方法，從2000—2019年夜間燈光影像中提取了典型研究區(qū)各縣區(qū)的12種燈光指數(shù)，并用于分析夜間燈光指數(shù)與MPI的關(guān)聯(lián)性。圖5展示了2019年重慶市F1-F12指數(shù)的空間分布示意圖。其中F1即為ANLI，F(xiàn)6即為T(mén)NLI。已有研究表明，ANLI和TNLI與MPI指數(shù)存在較高的關(guān)聯(lián)性，通過(guò)分別將ANLI，TNLI，MPI劃分為極低、較低、中等、較高和極高等5個(gè)等級(jí)，而對(duì)比發(fā)現(xiàn)，它們?cè)谕粋€(gè)縣區(qū)中的等級(jí)差距并不明顯。結(jié)合圖5與圖3進(jìn)行分析，F(xiàn)1，F(xiàn)2與MPI的分布存在較大的相似度： 高值區(qū)域集中在“一小時(shí)經(jīng)濟(jì)圈”，而低值區(qū)域分布主要分布于渝東北翼和渝東南翼。并且通過(guò)對(duì)比分級(jí)差異，可以發(fā)現(xiàn)F1，F(xiàn)6和MPI指數(shù)的主要在類別相差一個(gè)等級(jí)之內(nèi)。同時(shí)，F(xiàn)1與MPI處于同一等級(jí)間的縣域數(shù)為13個(gè)，占總數(shù)約34%，而F6與MPI處于同一等級(jí)間的縣域數(shù)為16個(gè)，占總數(shù)約42%。ANLI，TNLI和MPI的關(guān)系符合已有研究的結(jié)論，也表明本研究中提取的夜光指數(shù)是有效的。此外，F(xiàn)2，F(xiàn)3，F(xiàn)4，F(xiàn)5和F8的分布也大致符合與MPI類似的“一圈兩翼”的分布模式，表明這些指數(shù)同樣與MPI指數(shù)存在關(guān)聯(lián)性。

(a) F1指數(shù) (b) F2指數(shù) (c) F3指數(shù)

(d) F4指數(shù) (e) F5指數(shù) (f) F6指數(shù)

(g) F7指數(shù) (h) F8指數(shù) (i) F9指數(shù)

(j) F10指數(shù) (k) F11指數(shù) (l) F12指數(shù)圖5 2019年重慶市12種燈光指數(shù)空間分布Fig.5 Spatial distribution of 12 lighting indices in Chongqing in 2019

同時(shí)，為進(jìn)一步分析各個(gè)燈光指數(shù)與MPI指數(shù)的關(guān)聯(lián)性，以及評(píng)價(jià)通過(guò)夜光指數(shù)對(duì)MPI指數(shù)進(jìn)行估算的可行性，本研究使用Pearson相關(guān)系數(shù)法對(duì)F1—F12與MPI指數(shù)進(jìn)行了相關(guān)性分析。2000—2019年重慶市夜間燈光F1—F12指數(shù)與MPI指數(shù)的相關(guān)系數(shù)如表4所示，所有結(jié)果均通過(guò)了顯著性檢驗(yàn)。從表4中可以看出，夜間燈光的集中特征、分散程度和空間特征等6個(gè)指數(shù)，以及分布特征的5個(gè)指數(shù)均與MPI指數(shù)呈顯著正相關(guān)； 分布特征的F7指數(shù)與MPI指數(shù)呈顯著負(fù)相關(guān)。F1和F6指數(shù)與MPI指數(shù)的相關(guān)系數(shù)較其他大部分燈光指數(shù)更高，表明ANLI和TNLI指數(shù)能夠較好地反映貧困程度，符合已有研究結(jié)論，同時(shí)其他指數(shù)也存在與MPI指數(shù)的明顯相關(guān)性，說(shuō)明可以作為夜間燈光信息的補(bǔ)充。

表4 MPI和夜光特征變量的相關(guān)系數(shù)及顯著性Tab.4 Correlation analysis between MPI and NTL feature variables

綜上所述，本研究所選取的12種燈光指數(shù)與貧困程度存在明顯關(guān)聯(lián)性，并且使用這些燈光指數(shù)對(duì)MPI進(jìn)行估算是合理且可行的。

3.2 西南地區(qū)縣域多維貧困測(cè)度

根據(jù)夜光數(shù)據(jù)提取12種夜光指數(shù)結(jié)合PSO-BP神經(jīng)網(wǎng)絡(luò)構(gòu)建立MPI的估算模型，然后借助估算結(jié)果開(kāi)展基于夜光數(shù)據(jù)的貧困測(cè)度。同時(shí)，經(jīng)過(guò)驗(yàn)證表明該方法具有較強(qiáng)的可靠性和普適性。所以該模型可以擴(kuò)展到西南地區(qū)來(lái)對(duì)各縣域進(jìn)行長(zhǎng)期的貧困測(cè)度。圖6為西南地區(qū)NLMPI指數(shù)及國(guó)家級(jí)貧困縣的分布，由圖6可知，研究區(qū)中各縣域夜光估算NLMPI的空間分布具有明顯的隨時(shí)間變化的動(dòng)態(tài)特征。在整體分布上，研究區(qū)中的多維貧困分布呈現(xiàn)兩極分化的特征： 低值NLMPI指數(shù)主要分布在川西北、滇西北和云貴高原地區(qū)，伴隨自然環(huán)境和交通條件的不利，且具有連片分布的特征； 高值NLMPI則主要集中位于成都平原、滇中、渝中和黔中等地勢(shì)相對(duì)平坦的城市群，屬于政府優(yōu)先發(fā)展其經(jīng)濟(jì)以對(duì)周邊城市發(fā)揮帶動(dòng)作用的區(qū)域。此外，可以發(fā)現(xiàn)大部分貧困縣的NLMPI指數(shù)分級(jí)為最低，也表明本研究中NLMPI指數(shù)是合理的。

(a) 2000年 (b) 2005年 (c) 2010年

(d) 2015年 (e) 2019年 (f) 國(guó)家級(jí)貧困縣圖6 西南地區(qū)NLMPI指數(shù)及國(guó)家級(jí)貧困縣的分布Fig.6 Spatial patterns of NLMPI and national poverty counties in southwest China

從時(shí)間變化特征上看，在國(guó)家開(kāi)展扶貧工作和經(jīng)濟(jì)快速發(fā)展的背景下，整體NLMPI指數(shù)有顯著提高，2000年、2005年、2010年、2015年、2019年研究區(qū)的NLMPI指數(shù)的中位數(shù)分別為： 0.33，0.44，0.87，1.64和2.21(圖7)。自“新綱要”頒布實(shí)施以來(lái)，各地區(qū)的脫貧攻堅(jiān)任務(wù)取得了優(yōu)秀的成果。四川省、重慶市、云南省和貴州省在2013年的貧困發(fā)生率分別為： 9.6%，6.0%，17.8%和21.3%； 到2019年底，貧困率分別下降到了0.30%，0.12%，1.32%和0.85%，累計(jì)實(shí)現(xiàn)209個(gè)國(guó)家貧困縣的脫貧摘帽，而剩余的25個(gè)貧困縣在2020年也實(shí)現(xiàn)全部摘帽，脫貧攻堅(jiān)任務(wù)取得重大勝利，并標(biāo)志著區(qū)域性貧困在我國(guó)的消除。

圖7 不同時(shí)期西南地區(qū)NLMPI指數(shù)Fig.7 NLMPI of southwest China in different years

為驗(yàn)證NLMPI對(duì)貧困的識(shí)別準(zhǔn)確性，將研究區(qū)NLMPI指數(shù)利用自然斷點(diǎn)法分為了5級(jí)，分別為： 極高、較高、中等、較低和極低。對(duì)國(guó)家劃定貧困縣進(jìn)行了統(tǒng)計(jì)比較，結(jié)果如表5。在極低和較低2個(gè)類型區(qū)劃所占比例中，2015 年為 92.9%和60.6%； 所占比例比較大。而在其他類型中，貧困縣所占數(shù)目都較少?？梢?jiàn)，國(guó)家劃定的貧困區(qū)域，其主要分布在NLMPI值極低和較低縣域。2019年貧困縣所占比例大幅下降，這是由于到2019年為止，西南地區(qū)大部分貧困縣均已脫離絕對(duì)貧困，此時(shí)NLMPI的等級(jí)劃分主要反映為相對(duì)貧困。這進(jìn)一步說(shuō)明西南地區(qū)的絕對(duì)貧困狀況已經(jīng)基本消除，但是仍存在發(fā)展不均衡的現(xiàn)象。

表5 國(guó)家級(jí)貧困縣在每一等級(jí)縣域數(shù)量中所占數(shù)目及比例Tab.5 Number and proportion of national poverty counties in the counties of each grade

3.3 西南地區(qū)縣域多維貧困時(shí)序演進(jìn)特征

為便于分析NLMPI時(shí)間變化，統(tǒng)計(jì)了各年中NLMPI各級(jí)縣域的數(shù)量和比例(圖8)。由圖8可知，2000—2019年NLMPI極高和較高等級(jí)縣域個(gè)數(shù)增加到61個(gè)，綜合占比在8.74%～16.09%之間，而NLMPI極低、較低以及中等縣域綜合占比則在83.91%～91.26%之間，差距較大，也表明研究區(qū)中多維貧困狀況分化較為嚴(yán)重。2000—2019年近10 a來(lái)，多維貧困等級(jí)極低和較低的縣域占比下降，中等縣域占比提高，表明期間扶貧工作取得了成效。西南地區(qū)2000—2019年平均NLMPI指數(shù)變化如圖9所示。由圖9中可以看出，2000—2019年研究區(qū)整體NLMPI指數(shù)呈上升趨勢(shì)，說(shuō)明在扶貧開(kāi)發(fā)工作的開(kāi)展下，區(qū)域的發(fā)展得到較好的支持，多維貧困的整體狀況有較大改善。此外，如圖10所示，整個(gè)研究區(qū)NLMPI指數(shù)增長(zhǎng)速度在2011年前呈加快趨勢(shì)，而在2011—2019年時(shí)間段增長(zhǎng)速度放緩，造成這一現(xiàn)象的原因在于《中國(guó)農(nóng)村扶貧開(kāi)發(fā)綱要(2001—2010)》實(shí)施期間，扶貧工作的方式是“大水漫灌”，同時(shí)扶貧開(kāi)發(fā)的工作強(qiáng)度相較上一階段有明顯提高，使得這一階段的貧困人口有明顯減少，人民基本解決了溫飽問(wèn)題，所以2000—2010年NLMPI指數(shù)增長(zhǎng)較快，且增長(zhǎng)速率逐年提高； 而2011—2019年是我國(guó)扶貧工作的脫貧攻堅(jiān)階段，在經(jīng)過(guò)上一階段的開(kāi)發(fā)式扶貧工作后，貧困分布特征有顯著變化，剩余貧困人口趨于分散化，并主要集中在自然環(huán)境條件惡劣、交通不便的地區(qū)，同時(shí)這一時(shí)期政府實(shí)施了“精準(zhǔn)扶貧”政策，將扶貧對(duì)象瞄準(zhǔn)到戶和人，提高了減貧成本，導(dǎo)致脫貧難度加大，所以這一階段NLMPI指數(shù)增長(zhǎng)相較前一時(shí)間階段有所放緩。分省來(lái)看，研究區(qū)4省(市)NLMPI值呈現(xiàn)相同的先加速后減緩的變化趨勢(shì)，四川省、重慶市、云南省和貴州省的NLMPI在2000—2019年中的增長(zhǎng)速率分別為： 11.67%，13.23%，10.59%和11.21%，其中重慶市的增長(zhǎng)速度最高，同時(shí)在研究區(qū)4省(市)中最先完成國(guó)家貧困縣摘帽。此外，圖9中顯示2013—2014年NLMPI指數(shù)的斜率有所下降，分析可能是2種燈光相互校準(zhǔn)造成的誤差，因?yàn)橄嗷バ?zhǔn)只能作為一種提升燈光連續(xù)性的手段，無(wú)法實(shí)現(xiàn)對(duì)真值的完全模擬。

圖8 西南地區(qū)不同時(shí)期各貧困等級(jí)縣域數(shù)量比Fig.8 Percentage of counties at different levels of poverty in southwestern China

圖9 2000—2019年西南地區(qū)平均NLMPI指數(shù)的時(shí)間變化特征Fig.9 Temporal variations of the average NLMPI in southwest China from 2000 to 2019

(a) 2000—2010年 (b) 2011—2019年圖10 西南地區(qū)不同時(shí)期NLMPI指數(shù)變化趨勢(shì)Fig.10 Trend of MPI in southwest China during the different periods

3.4 西南地區(qū)縣域貧困空間分異特征

本研究計(jì)算了西南地區(qū)2000—2019年全局Moran’s I指數(shù)，以此測(cè)量整體空間相關(guān)關(guān)系。其中所有結(jié)果均通過(guò)了0.01的顯著性檢驗(yàn)，結(jié)果如圖11所示。從整體上看，西南地區(qū)2000—2019年的Moran’s I值的平均值為0.478，表明西南地區(qū)的貧困指數(shù)存在正的空間自相關(guān)性，即表現(xiàn)為高值(低值)NLMPI空間聚集的狀態(tài)。另一方面，研究區(qū)整體Moran’s I值還呈現(xiàn)了一個(gè)先下降后上升的趨勢(shì)，這主要由于2000—2010年期間我國(guó)實(shí)施了《中國(guó)農(nóng)村扶貧開(kāi)發(fā)綱要(2001—2010年)》，明確以減少貧困人口作為工作重點(diǎn)，并推行整村推進(jìn)政策，所以這一階段貧困聚集現(xiàn)象有所減弱； 而2010—2019年階段Moran’s I值上升的主要原因?yàn)椋?經(jīng)過(guò)大規(guī)模開(kāi)發(fā)式扶貧后，遺留的貧困區(qū)域逐漸分散化，且由于致貧原因較為多樣化，增加了減貧成本，脫貧難度加大。

圖11 2000—2019年西南地區(qū)的Morans’ I指數(shù)Fig.11 Morans’I index for southwest China during the period of 2000—2019

西南地區(qū)各縣域2000—2019年NLMPI指數(shù)的局部空間自相關(guān)分布狀況如圖12所示。整體而言，西南地區(qū)多維貧困度分布整體呈現(xiàn)改善但分布不均衡的狀態(tài)。2000—2019年H-H，L-L，H-L和L-H這4種類別縣域的數(shù)量明顯失衡，其中H-H類別縣域數(shù)量呈增多趨勢(shì)，H-L和L-H類別縣域數(shù)量較少，L-L類別數(shù)量先增加后減少，總體保持穩(wěn)定。但L-L類別中貧困縣的占比在2016年之后保持下降趨勢(shì)，與貧困縣摘帽的時(shí)間段吻合，說(shuō)明西南地區(qū)的貧困縣整體貧困度正在下降。2000—2019年研究區(qū)H-H區(qū)域主要集中在以成都市、昆明市以及貴陽(yáng)市等省會(huì)城市和重慶市主城區(qū)為中心的城市群，區(qū)域范圍逐年增長(zhǎng)，符合實(shí)際情況。L-L區(qū)域分布較廣，主要分布在川西北、滇西北、重慶市兩翼和貴州東部等地區(qū)，多為貧困縣。H-L區(qū)域分布較為零散，多為位于貧困縣周邊的非貧困縣，如景洪市、大理市和西昌市等，因此這些縣域被具有較低NLMPI指數(shù)的縣域包圍，形成高低聚集的現(xiàn)象。L-H區(qū)域中的縣域則主要分布于H-H區(qū)域周邊，因?yàn)橹行膮^(qū)域NLMPI值相對(duì)較高，所以呈一種“包圍”的狀態(tài)。

(a) 2000年 (b) 2005年 (c) 2010年圖12-1 2000—2019年西南地區(qū)各縣域NLMPI指數(shù)的局部空間同相關(guān)分布狀況Fig.12-1 Distribution of NLMPI index of local spatial autocorrelation in southwest China from 2000 to 2019

(d) 2015年 (e) 2019年圖12-2 2000—2019年西南地區(qū)各縣域NLMPI指數(shù)的局部空間同相關(guān)分布狀況Fig.12-2 Distribution of NLMPI index of local spatial autocorrelation in southwest China from 2000 to 2019

4 結(jié)論

本研究建立了一種基于夜間燈光的貧困測(cè)度方法： 通過(guò)整合DMSP-OLS和NPP-VIIRS這2種夜間燈光數(shù)據(jù)構(gòu)建長(zhǎng)時(shí)間序列夜間燈光數(shù)據(jù)集，并基于整合的長(zhǎng)時(shí)間序列夜光數(shù)據(jù)開(kāi)展貧困測(cè)度。同時(shí)，經(jīng)過(guò)驗(yàn)證表明該方法可以更好地探究區(qū)域的貧困時(shí)空動(dòng)態(tài)，所以該模型可以擴(kuò)展到西南地區(qū)來(lái)對(duì)各縣域進(jìn)行長(zhǎng)期的貧困測(cè)度。研究結(jié)果表明：

1)在2000年，西南地區(qū)各縣域的多維貧困狀況分化較為嚴(yán)重，但隨國(guó)家扶貧工作的開(kāi)展，極低和較低等級(jí)縣域占比下降，中等縣域占比提高，表明期間扶貧工作取得了顯著的成效。

2)2000—2019年各縣域的Moran’s I指數(shù)呈現(xiàn)了先下降后上升的趨勢(shì)，這反映出在2000—2010年，貧困聚集現(xiàn)象有所減弱，而在之后進(jìn)入了較為分散的脫貧攻堅(jiān)階段。

3)局部空間自相關(guān)結(jié)果表明，研究區(qū)的多維貧困度分布整體呈現(xiàn)改善但分布不均衡的狀態(tài)，表現(xiàn)為以成渝片區(qū)、昆明市、貴陽(yáng)市為主的高高聚集和以川西北、滇西地區(qū)為主的低低聚集的分布格局。

基于夜間燈光數(shù)據(jù)，本文監(jiān)測(cè)了西南地區(qū)各縣域的貧困時(shí)空變化特征。但是由于夜光遙感數(shù)據(jù)的限制性，未來(lái)的貧困與社會(huì)經(jīng)濟(jì)發(fā)展研究應(yīng)該加入高分辨率的輔助數(shù)據(jù)，增強(qiáng)夜光數(shù)據(jù)在不同區(qū)域尺度研究中的應(yīng)用能力。