SNPP VIIRS DBN日尺度夜间灯光数据下载

数据介绍

VNP46（NASA’s Black Marble nighttime lights product suite）数据是由美国国家航空航天局土地科学研究人员领导的处理系统(SIPS)为常规全球处理而开发和实施，利用所有高质量、无云、大气、地形、植被、雪、月和杂散光校正辐射来估计每日夜间光线(NTL)和其他内在表面光学特性。该数据集是一个逐日的夜间辐射产品。

该产品包括VNP46A1和VNP46A2两个数据产品。 其中VNP46A1是逐日到达传感器的夜间辐射，该数据更新较快，去除了大气，地形和云层影响，但没有消除季节性和与月球相关变化的影响，受月亮阴晴圆缺影响非常大，基本没法直接用；而VNP46A2则是校正了月光影响的逐日数据，下载下来可以直接用，但该数据目前只有2012-2018年的，更新较慢，VNP46产品的介绍可以参考这篇文章（Román, M.O., Wang, Z., Sun, Q., Kalb, V., Miller, S.D., Molthan, A., Schultz, L., Bell, J., Stokes, E.C., Pandey, B. and Seto, K.C., et al. (2018). NASA’s Black Marble nighttime lights product suite. Remote Sensing of Environment 210, 113-143. doi:10.1016/j.rse.2018.03.017.）。接下来，我将介绍VNP46A1数据集及其下载过程，VNP46A2下载方式相同。

VNP46A1数据集包含26个科学数据集（SDS），包括传感器辐射，天顶和方位角（在传感器，太阳和月球），云遮罩标志，时间，短波红外辐射， 亮度温度，VIIRS质量标志，月相角和月光照度。 它还提供特定于云掩码，VIIRS中分辨率波段M10，M11，M12，M13，M15，M16和DNB的质量标志（QF）信息。产品信息如图1。该数据更新非常快，如2020年11月27日就已经可以下载26日的数据了。该数据介绍见网址[https://ladsweb.modaps.eosdis.nasa.gov/missions-and-measurements/products/VNP46A1]。

图1 数据介绍

如何下载该数据

1、在主页中点击“search for Product Files”，即可跳转到数据检索页面，且定位到了该产品(https://ladsweb.modaps.eosdis.nasa.gov/search/order/2/VNP46A1–5000/2020-11-13…2020-11-26/DB/World)

图2 数据介绍

2、选择你所需要的时间，这里我选择的是11月13日到11月26日（我写博客的日期是11月27日），然后点Add Date，则该日期被添加到了右边。点击右上角的箭头符号，进入到区域选择“Location”界面。

图3 数据时间选择

3、选择你所需要的区域，区域选择有多重方式。“World”表示下载全球的，如图4；“country”代表按国家来选择，如点击中国，则会检索出整个中国的数据；“Tiles”则表示选择某一瓦片，如图5，如下面第二张图所示；“Validation sites”表示验证点，如图6；“draw custom box”表示按自己绘制区域，如图7；“Enter coordinates ”表示手动输入经纬度，如图8.完成后，点击右上角的箭头符号，进入到files界面，查看检索出来的文件。本次采用的是自己绘制区域（大致为北京）

图4 区域选择页面，“World”

图5 Tiles选择

图6 “Validation sites”选择

图7“draw custom box” 选择

图 8 Enter coordinates 选择

4、查看检索的文件。我所选的区域是北京，下图显示了所有的文件。这里可以点击每个文件后面的下载按钮一个一个下载（文件少的时候推荐），也可以选择“select all”全部选中，则界面变成图9所示。再点击右上角的箭头符号进入下一步的Review & Order界面。

图9 查看检索的文件

5、当你点击下载的时候或者进入Review & Order界面之前会让你登陆Earthdata的账户，如图10界面所示。如果你没有账户的话则需要注册一个。

图10 登陆Earthdata的账户

6、登陆完成之后则会进入Review & Order界面，点击submit order，则可提交订单。等几分钟后便会收到邮件，此时则表示该数据可以下载了。

图11 Review & Order界面

7、点击“Past orders”，则进入所有的订单。

图12 “Past orders” 界面

8、这里我的订单是501520114。

图13 “Past orders” 界面的订单展示

9、点击该订单将其展开，即可看到详细信息，archive/orderd/501520114，即可跳转到下载页面。

图14 将订单展开图

10、在该下载页面中，可以选中所有数据，然后点击“Download All”，然后数据会逐条开始下载，不过下载每条数据的时候都会弹出“另存为”对话框，导致我们还是得守在电脑前，不太方便。下面还有另一种方法。

图15 订单中的数据下载页面

11、第二种下载方法是利用wget代码。点击“See wget Download Command”，可以看到下面有两行代码，这行代码则是用于批量下载该订单数据的。如果你没有“wget.exe”的话，首先要下载该程序。你可以点击页面中的蓝色字体“wget”（如图16红色圈圈里），即可跳转到wget工具介绍页面，里面有工具下载的网址以及安装方式。

图16 点击“See wget Download Command”后的显示

图17 wget工具介绍页面

12、嫌麻烦可以直接点击wget网址，即可下载。我下载的是版本1.20.3,64位的exe，其他人根据你的电脑是32位还是64位来选择，不过一定要下载exe的那个。下载完成后，把它放到C:\Windows\System32里（需要管理员权限），后续则可以使用该工具了。

图18 wget下载页面

图19 wget程序放到C:\Windows\System32里

13、将订单页面中“See wget Download Command”下面的两行代码复制到记事本，将APP_KEY 改为你的Earthdata账户的APP_KEY（是一串数字加字母），将TARGET_DIRECTORY_ON_YOUR_FILE_SYSTEM改为数据要存在的文件夹路径，我这里是E\NTLdata。代码为 wget -e robots=off -m -np -R .html,.tmp -nH --cut-dirs=3 “https://ladsweb.modaps.eosdis.nasa.gov/archive/orders/501520114/” --header “Authorization: Bearer APP_KEY” -P TARGET_DIRECTORY_ON_YOUR_FILE_SYSTEM

图20 代码修改

14、改完之后将代码复制到命令提示符里，运行即可自动下载订单里的数据，没有出错你就可以做其他的了，不用一直等着。

图21 运行wget代码

15、有一个小问题就是，我订单里有26条数据，但我的下载完19条后就退出下载了，这是因为你的earthdata账户退出登陆了，导致你的下载也终止了，解决方法就是刷新一下数据下载页面，让他再自动登录一次，然后再运行那段wget代码就可以了（注意，不要改动目标文件夹的任何文件，因为程序首先会判断目标文件夹里是否已经存在了该数据，存在的话就跳过该数据的下载。如果你改了文件名的话，就又会下载该数据）。

图22 只下载完19条

16、最后当文件夹里有checksums文件时，则表示全部数据下载完成。

图23 下载完成