DMSP/OLS夜間燈光影像由于傳感器未經過星上定標，導致長時間序列的各期影像之間缺乏連續性和可比性，且由于光譜分辨率的限制，使得大量的城市中心像元出現DN值為63的集聚飽和現象等問題，因此學者們對以上存在的問題提出相應的解決辦法。2009年，文獻[12]提出了不變目標區域的方法對全球的夜間燈光影像的不連續和不可比問題進行了校正，該方法以Sicily為不變動的參考區域，以F121999為參考數據集，將待校正影像與F121999參考數據集進行比較建立回歸方程。該方法較適用于全球夜間燈光數據的校正，對于某一國家尺度的適用性不強。文獻[3]在Elvidge的基礎上，提出了適合中國區域的夜間燈光影像的校正方法。首先統計分析了中國主要城市1992—2008年的GDP統計數據和建成區面積數據，確定黑龍江省的雞西市為校正中國區域的不變動參考區域，選擇F16的2007年數據作為參考數據集，實現長時間序列DMSP/OLS夜間燈光數據的相互校正、年內融合和年際間校正，卻忽略了影像的飽和問題。文獻[13]基于地理探測器模型驗證分析，提出利用單元路網長度去飽和算法。結合2006年的輻射定標產品數據對2013年DMSP/OLS穩定燈光數據進行飽和校正。文獻[14]結合DMSP/OLS數據和Landsat數據，建立去除水體并且緩解燈光強度飽和的指數RwNTLI，結果表明該方法在一定程度上降低了影像的飽和問題。以上幾位學者提出的方法能夠在一定程度上緩解影像的飽和問題，但是僅局限在單期影像的校正，無法在長時間序列上進行研究。因此，本文旨在對1992—2013年長時間序列的DMSP/OLS夜間燈光影像進行飽和校正、相互校正、年內融合及年際間校正等一系列校正內容，從而解決夜間燈光影像的飽和、影像間不連續及不可比等的問題。

1 校正方法

具體校正思路為：①確定黑龍江省的雞西市為參考區域；②選取已經經過飽和校正的輻射定標產品數據中的2006年F16傳感器的數據作為對連續穩定燈光影像進行飽和校正、相互校正的參考數據集；③對1992年至2013年34期夜間燈光影像進行相同年份的年內融合及不同年份的年際間校正，從而解決燈光影像的飽和、影像間不連續及不可比等問題。

1.1 影像預處理

原始穩定燈光影像的坐標為WGS-84，空間分辨率為30″。研究區域為中國區域，因此根據原始影像裁剪出中國區域影像。為了使投影面積變形最小，對坐標系進行變換，橢球體為D_Krassovsky_1940，投影坐標為Krassovsky_1940_Albers。統一影像分辨率為1000 m。

1.2 中國區域的DMSP/OLS夜間燈光數據校正1.2.1 影像數據分析

圖 1所示為未經校正的中國區域的所有亮值像元(像元DN值大于0)的數量，圖 2所示為未經校正的中國區域的所有亮值像元的DN值總和(TDV)。可以看出：①同一年份的數據可由不同傳感器進行獲取，且不同傳感器獲取的亮值像元總數量不一致、總DN值不一致；②同一傳感器獲取的不同年份的像元總數量、總DN值存在異常波動；③由于以上兩個問題的存在，導致多期夜間燈光影像間數據不連續、不可比；④由于傳感器光譜分辨率的限制，使得城市中心區域的像元值集聚飽和，DN值最大為63，使得城市中心地帶的差異不明顯。TDV的計算公式為

圖 1 校正前的亮值像元總數

圖 2 校正前的亮值像元總DN值

式中，DNi表示i像元的DN值；Ni表示i像元的數量。

1.2.2 相互校正與飽和校正

將1992—2013年的34期待校正的穩定燈光影像與輻射定標產品數據中的2006年F16傳感器的數據進行相互校正、飽和校正。不變動的參考區域選擇：黑龍江省的雞西市。分別對待校正影像和輻射定標產品數據裁剪出雞西市的數據，分別統計DN值，進行如式(2)的二次回歸模型計算，將模型結果運用到待校正的34期影像，完成影像的飽和校正、相互校正。

式中，DN表示校正前的亮度值；a、b、c為回歸系數；DNcorrect表示校正后的DN值。所求結果見表 1 。

表 1 穩定燈光影像二次回歸模型的模型參數

衛星

年份

a

b

c

R2

F10

1992

0.036

0.337

2.581

0.850

1993

0.045

0.334

2.150

0.850

1994

0.039

0.518

1.252

0.850

F12

1994

0.039

0.599

0.277

0.854

1995

0.034

0.513

0.485

0.851

1996

0.037

0.479

0.511

0.862

1997

0.030

0.560

0.413

0.877

1998

0.027

0.511

0.481

0.876

1999

0.029

0.502

0.454

0.866

F14

1997

0.035

0.893

0.165

0.863

1998

0.036

0.731

0.229

0.865

1999

0.028

0.901

0.124

0.863

2000

0.024

0.606

0.346

0.873

2001

0.021

0.780

0.188

0.862

2002

0.017

1.018

-0.113

0.876

2003

0.017

1.096

-0.101

0.890

F15

2000

0.028

0.578

0.485

0.878

2001

0.032

0.560

0.487

0.869

2002

0.023

0.631

0.356

0.888

2003

0.020

1.173

-0.215

0.855

2004

0.024

1.146

-0.153

0.901

2005

0.026

0.872

0.123

0.855

2006

0.015

1.076

-0.081

0.894

2007

0.018

0.977

0.086

0.892

F16

2004

0.022

0.809

0.231

0.866

2005

0.022

0.919

0.096

0.887

2006

0.019

0.920

0.100

0.924

2007

0.019

0.701

0.310

0.888

2008

0.012

0.625

0.381

0.868

2009

0.018

0.512

0.482

0.872

F18

2010

0.021

0.146

1.525

0.851

2011

0.020

0.510

0.901

0.851

2012

0.017

0.483

0.508

0.852

2013

0.016

0.271

2.680

0.850

1.2.3 年內融合

不同傳感器獲取的同一年份的數據不一致，為了充分利用各獨立傳感器獲取的數據，同時為了能夠解決傳感器獲取數據的不連續問題，研究按照式(3)對相互校正及飽和校正后的部分影像進行年內融合。需要進行年內融合的年份包括：1994、1997—2007年。年內融合包括：①兩幅影像對應的某像元的DN值，均為0，則該像元的DN值為0；②否則，取二者的平均值作為該像元的DN值。

式中，n=1994, 1997, 1998, …, 2007;DN(n, i)a、DN(n, i)b分別表示相互校正、飽和校正后n年的兩個不同傳感器獲取的i像元的DN值；DN(n, i)表示影像年內融合校正后n年的i像元的DN值。

1.2.4 年際間校正

經過年內融合之后，仍然存在不同年份之間的影像不可比的現象，因此對其進行年際間校正。校正依據是：后一年的像元DN值應該不小于前一年同一位置的像元DN值。如式(4)所示，年際間校正包括：①后一年的某像元DN值為0，則前一年的同一位置的像元DN值也為0；②前一年的某像元的DN值大于后一年同位置的像元DN值，則將前一年的DN值賦值給后一年該位置的像元DN值；③否則，后一年的像元DN值為本身的DN值。

式中，n=1992, 1993, 1994, …, 2013；DN(n-1, i)、DN(n, i)、DN(n+1, i)分別表示影像的像元i在n-1、n、n+1年的DN值。

1.2.5 影像校正結果

對1992—2013年的夜間燈光影像進行以上相互校正、飽和校正，年內融合及年際間校正，校正結果如圖 3所示(由于篇幅限制，文中僅列出中國區域1995、2000、2005及2010年的影像的校正結果)。

圖 3 1995、2000、2005及2010年夜間燈光影像校正結果

2 檢驗校正結果

定量檢驗對中國區域夜間燈光影像校正的合理性，分別從2個角度進行校驗：①對校正后影像的亮值像元總數和亮值像元總DN值進行統計，校正后的像元總數和像元總DN值變化趨勢如圖 4、圖 5所示；②早期文獻[15—16]的研究結果表明，GDP數據與夜間燈光數據在國家或大區域空間尺度上存在線性相關關系，因此，本文研究在省級尺度上對2013年校正前后影像的TDV值進行對比，比較省級尺度GDP與校正前后的TDV值之間的線性相關性，校正前后的相關性如圖 6、圖 7所示。

圖 4 校正后的亮值像元總數

圖 5 校正后的亮值像元總DN值

圖 6 2013年夜間燈光數據校正前的TDV值

圖 7 2013年夜間燈光數據校正后的TDV值

校正后的夜間燈光數據每一年的亮值像元總數和總DN值較未校正前的數據均呈現逐漸增長的趨勢。經過校正，充分利用各傳感器獲取的數據，每一年均只有一份有效數據，影像的一致性和連貫性均得到提升。

以2013年的夜間燈光數據為例，對校正前后的夜間燈光TDV值進行統計，并與中國各個省GDP統計數據(除香港、澳門、臺灣的統計數據)進行線性關系比較，發現校正前的TDV值與GDP之間的相關性并不高，且數據點雜亂分布，不呈現明顯的相關性，R2僅為0.333 3。校正后的TDV值與GDP的相關性較高，R2為0.791 8，影像數據得到了較好的校正。

3 結 語

本文研究對中國區域1992—2013年的DMSP/OLS夜間燈光影像進行了校正，包括：相互、飽和校正，年內融合，年際間校正等步驟，并對校正結果從一致性、連貫性及與GDP的相關性等角度進行驗證，校正結果較好，可使DMSP/OLS夜間燈光數據充分發揮其功能性和時效性，后續可用于GDP、人口等統計數據的長時間序列擬合，用于城市群演化、城鎮空間擴展等研究。

盧秀, 李佳, 段平, 等. 中國區域DMSP/OLS夜間燈光影像的校正[J]. 測繪通報，2019(7)：127-131, 159.

DOI:10.13474/j.cnki.11-2246.2019.0234.

作者簡介:盧秀,女,博士生,主要從事資源環境遙感、夜間燈光遙感研究。E-mail:lx_rsgis@163.com

總結

以上是生活随笔為你收集整理的夜间灯光数据dn值_中国区域DMSP/OLS夜间灯光影像的校正的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。