珞珈一号夜间灯光数据的福建省人为热通量估算

珞珈一号夜间灯光数据的福建省人为热通量估算 Estimation of anthropogenic heat flux of Fujian Province （China） based on Luojia 1-01 nighttime light data 1 first-author 林中立，研究方向为城市环境遥感研究。E-mail：linzl@fjut.edu.cn 林中立，研究方向为城市环境遥感研究。E-mail：linzl@fjut.edu.cn linzl@fjut.edu.cn 2 3 corresp 徐涵秋，研究方向为环境遥感应用研究。E-mail：hxu@fzu.edu.cn 徐涵秋，研究方向为环境遥感应用研究。E-mail：hxu@fzu.edu.cn hxu@fzu.edu.cn 1 福建工程学院 建筑与城乡规划学院, 福州 350118 福建工程学院 建筑与城乡规划学院, 福州 350118 College of Architecture and Urban Planning, Fujian University of Technology, Fuzhou 350118, China College of Architecture and Urban Planning, Fujian University of Technology, Fuzhou 350118, China 福州大学 环境与安全工程学院 空间数据挖掘与信息共享教育部重点实验室, 福州 350116 福州大学 环境与安全工程学院 空间数据挖掘与信息共享教育部重点实验室, 福州 350116 Ministry of Education Key Laboratory of Spatial Data Mining & Information Sharing, College of Environmental and Safety Engineering, Fuzhou University, Fuzhou 350116, China Ministry of Education Key Laboratory of Spatial Data Mining & Information Sharing, College of Environmental and Safety Engineering, Fuzhou University, Fuzhou 350116, China 福州大学 遥感信息工程研究所 福建省水土流失遥感监测评估与灾害防治重点实验室, 福州 350116 福州大学 遥感信息工程研究所 福建省水土流失遥感监测评估与灾害防治重点实验室, 福州 350116 Fujian Provincial Key Laboratory of Remote Sensing Soil Erosion and Disaster Prevention, Institute of Remote Sensing Information Engineering, Fuzhou University, Fuzhou 350116, China Fujian Provincial Key Laboratory of Remote Sensing Soil Erosion and Disaster Prevention, Institute of Remote Sensing Information Engineering, Fuzhou University, Fuzhou 350116, China 夜间灯光数据是估算人为热通量（AHF）的重要数据，但当前应用最广的DMSP/OLS和Suomi-NPP/VIIRS夜间灯光数据由于受限于粗糙的空间分辨率，而无法刻画城市内部的AHF分布细节。中国2018年6月发射的Luojia 1-01卫星所获取的130 m高空间分辨率夜间灯光数据，则有望解决这一问题。因此本文利用Luojia 1-01夜间灯光数据，通过将统计年鉴中的能源统计数据细化至福建省84个县（市、区），然后与3个夜间灯光指数（NTL nor 、HSI、VANUI）进行回归分析，分别构建了基于这3个指数的福建省AHF空间估算模型，并采用交叉验证法对其进行筛选。结果显示：（1）在3个指数中，基于VANUI的乘幂估算模型的 R 2 最高，且RMSE最小，因此精度最高；（2）利用VANUI乘幂估算模型反演得到的2018年福建省年均AHF为0.88 W/m 2 ，其中厦门市的年均AHF最高，达10.98 W/m 2 ，泉州、莆田、福州、漳州等沿海城市次之，年均值在0.98—1.95 W/m 2 ，而宁德、龙岩、三明、南平等城市的AHF则较低，均值在0.38—0.46 W/m 2 ；（3）Luojia 1-01夜间灯光数据可以揭示城市内部的AHF分异细节。根据用地属性和功能的不同，AHF数值表现为：城市集中商业区 > 大型市政公共设施区 > 城市主干道 > 城市住宅区 > 近郊住宅区。研究表明，基于Luojia 1-01夜间灯光数据建立的AHF估算模型可以较好地揭示城市尺度AHF的空间分异情况。 Nighttime light (NTL) data are important for estimating Anthropogenic Heat Flux (AHF). However, the commonly used DMSP/OLS and Suomi-NPP/VIIRS NTL data are restricted by their coarse spatial resolution and therefore, cannot exhibit the spatial details of AHF at city scale. The 130 m high-resolution NTL data obtained by the Luojia 1-01 satellite launched in June 2018 show potential to solve this problem. Therefore, this study aims to construct an AHF estimation model using the NTL data of Luojia 1-01 for Fujian Province based on three indexes, namely, normalized nighttime light data (NTL nor ), Human Settlement Index (HSI), and Vegetation Adjusted NTL Urban Index (VANUI). To determine the best estimation model of AHF, the AHF of 84 county-level cities of Fujian Province has also been calculated using energy-consumption statistics data and then correlated with the corresponding data of three indexes. Results show that (1) based on a five-fold cross validation approach, VANUI power estimation model achieves the highest R 2 along with the smallest RMSE; therefore, it has the highest accuracy among the three indexes; (2) according to the VANUI power estimation model, the average annual AHF of Fujian Province in 2018 is 0.88 W/m 2 , of which Xiamen has the highest average annual AHF of 10.98 W/m 2 , followed by Quanzhou, Putian, Fuzhou, and Zhangzhou, with the annual average of 0.98—1.95 W/m 2 , whereas the figures of Ningde, Longyan, Sanming, and Nanping are relatively low, ranging from 0.38—0.46 W/m 2 ; (3) Luojia 1-01 NTL data can reveal the AHF differentiation details at a city scale. The AHF values of different land properties and functions show the following order: urban commercial area > large municipal public facility area > urban main road > urban residential area > suburban residential area. Studies have shown that the AHF estimation model developed by Luojia 1-01 NTL data can achieve high accuracy of the city-scale estimation of AHF. 遥感 人为热 AHF 珞珈一号01星 夜间灯光影像 福建省 remote sensing anthropogenic heat AHF Luojia 1-01 nighttime light imagery Fujian Province 1　引言 1 当前，城市化进程正在全球范围内以前所未有的速度进行。城市化的快速发展伴随着城市空间的不断扩张，自然地表被大面积的人工建筑所取代，城市热岛效应逐步加剧（ Rizwan等，2008 Rizwan等，2008 ； Imhoff等，2010 Imhoff等，2010 ）。同时，城市人口数量激增，根据联合国人口司统计数据显示，世界城市人口比例从1950年的30%增加至2018年的55.3%，预计到2050年，世界城市人口比例将会达到68%（United Nations，2019）。密集的人类生产和生活消耗了大量不同形式的能源，最终除了向大气排放温室气体外，还释放大量热量，这种热量即为人为热（Anthropogenic Heat）（ Oke，1976 Oke，1976 ； 陈兵 等，2016 陈兵 等，2016 ）。城市地区大量的人为热排放可直接加热地表和近地表大气，影响局地能量平衡，加剧城区及其周边区域温度的升高，对局地气候，尤其是城市热岛效应产生不可忽视的影响（ Block，2004 Block，2004 ； 何晓凤 等，2007 何晓凤 等，2007 ； Pal等，2012 Pal等，2012 ）。 Torrance和Shun（1976） Torrance和Shun（1976） 是最早关注城市人为热排放及其对气候影响的学者。随后，许多学者也投入到人为热的研究当中，对城市和区域的人为热排放进行估算（ Ichinose等，1999 Ichinose等，1999 ； 佟华 等，2004 佟华 等，2004 ； Bohnenstengel等，2014 Bohnenstengel等，2014 ）。研究者通常以人为热通量AHF（Anthropogenic Heat Flux），即在单位时间和单位面积内所产生的人为热排放总通量（W/m 2 ）来定量衡量人为热排放的数量（ Iamarino等，2012 Iamarino等，2012 ）。 Wang等（2019） Wang等（2019） 对2016年中国地区的AHF进行估算，结果表明，在省级水平上，上海的年均AHF最高，达12.53 W/m 2 ，天津、北京和广东次之，分别为6.91 W/m 2 、5.84 W/m 2 和4.53 W/m 2 。 Flanner（2009） Flanner（2009） 根据研究的模型预测，到2040年美国大陆、西欧以及中国地区的年均AHF将分别达到0.59 W/m 2 、0.89 W/m 2 和0.76 W/m 2 ，同时他还指出全球年均AHF的不断增加，将会加剧全球平均温度的升高，并且致使行星边界层内的气溶胶等大气污染物释放进入AHF较高的地区。可见，人为热排放是当前全球城市，尤其是快速发展大背景下的中国城市的温度升高、环境质量下降等问题的不可忽视的重要原因。 夜间灯光数据能在夜间探测到地球表面微弱的近红外辐射，包括城市甚至小规模居民地的灯光、车流和其他持久性光源，以及云、油气燃烧等短暂性光源（ Elvidge等，1997 Elvidge等，1997 ）。基于这一特性，夜间灯光数据能够较好地反映社会经济活动与能源消耗的空间分布情况，已有许多学者发现了夜间灯光亮度值与AHF之间存在显著的相关关系（ Chen等，2012 Chen等，2012 ； 马盼盼 等，2016 马盼盼 等，2016 ； Wang等，2019 Wang等，2019 ），利用这一特性可进行AHF的空间估算。 Wang等（2019） Wang等（2019） 利用NPP/VIIRS夜间灯光正规化亮度值NTL nor （Normalized nighttime light data）获得了与AHF较高的拟合度（ R 2 =0.95）。一些学者发现，将夜间灯光数据结合植被指数能够在一定程度上消除夜间灯光数据过饱和的问题，有效提高城市夜间灯光亮度的分异性与区分度，在城市开发较早、发展速度平稳的地区可以有效凸显城市核心区的灯光亮度与周边临近区域的差异（ 陈颖彪 等，2019 陈颖彪 等，2019 ），从而提高夜间灯光与AHF间的拟合度。 Lu等（2008） Lu等（2008） 结合DMSP/OLS夜间灯光数据与MODIS NDVI植被产品，构建了人居指数HSI（Human Settlement Index）。 Chen和Hu（2017） Chen和Hu（2017） 利用HSI对北京市县区级尺度的AHF进行拟合，得到了较为显著的拟合度（ R 2 =0.989）。 Zhang等（2013） Zhang等（2013） 则指出在灯光饱和区，可能出现NDVI接近0，从而导致HSI呈指数级增长，造成矫枉过正的问题，因此他们在HSI的成果基础上，构建了植被调节城市夜间灯光指数VANUI（Vegetation Adjusted NTL Urban Index）。 夜间灯光数据是AHF定量估算的重要数据之一，但一直以来，应用最为广泛的夜间灯光数据为美国前后两代的夜间灯光数据（DMSP/OLS和Suomi-NPP/VIIRS），但其分辨率仅为30 arc-second（≈1000 m）和15 arc-second（≈500 m），基本无法揭示城市内部AHF的空间分异细节。2018年6月武汉大学成功发射了Luojia 1-01卫星（ http：//www.lmars.whu.edu.cn/index.php/xwzx/2168.html http://www.lmars.whu.edu.cn/index.php/xwzx/2168.html ［2020-06-22］），提供了空间分辨率为130 m，辐射分辨率为14 bit的夜间灯光影像，显著改善了当前夜间灯光影像空间分辨率低的问题（ Jiang等，2018 Jiang等，2018 ； 李德仁 等，2019 李德仁 等，2019 ； 钟亮和刘小生，2019 钟亮和刘小生，2019 ； Wang等，2020 Wang等，2020 ），运用Luojia 1-01夜间灯光数据有望显著提高AHF估算的空间分辨率，然而针对这一方面的研究和应用迄今还鲜有报道。福建省是中国东部沿海的重要省份，当前针对福建省人为热排放的相关研究仍鲜有报道。鉴于此，本文基于Luojia 1-01夜间灯光数据，构建2018年福建省AHF的遥感估算模型，实现省域范围的AHF高空间分辨率估算，并对福建省AHF的空间分布规律进行分析和探讨。因此，本研究对Luojia 1-01夜间灯光数据的科学应用及其省域尺度人为热排放的研究都具有重要意义。 2　研究区和数据源 1 2.1　研究区概况 2 福建省位于中国东南沿海，地理位置为23°33′ N—28°19′ N，115°50′ E—120°43′ E，全省陆地总面积为124000 km 2 。福建省地形以山地丘陵为主，占全省总面积的80%以上，森林覆盖率达67%，长期居全国首位。福建省属亚热带海洋性季风气候，温暖湿润。2010年—2018年，总人口从3693万人增至3941万人，地区生产总值从14737亿元增至35804亿元，能源消费方面也出现较大的增幅，从9189万吨标准煤增至13131万吨标准煤（ 福建省统计局和国家统计局福建调查总队，2019 福建省统计局和国家统计局福建调查总队，2019 ）。 福建省共有9个地级市（福州、厦门、漳州、泉州、三明、莆田、南平、龙岩和宁德），下辖85个县（市、区），其中泉州市的金门县，以及一些面积较小的沿海岛礁不在此次研究区范围内，因此研究区总面积为121518 km 2 （ 图1（a） 图1（a） ）。 1,2; 2.2　遥感影像数据源 2 本文所用Luojia 1-01夜间灯光影像来自高分辨率对地观测系统湖北数据与应用网（ http：//www.hbeos.org.cn http://www.hbeos.org.cn ［2020-06-22］）发布的GEC系统几何纠正产品。Luojia 1-01于2018年6月2日发射升空后，2018年全年先后过境福建省10次，经逐影像检查，剔除其中4次质量不佳（大面积云对信息的遮挡、薄云造成的影像模糊与灯光溢出）的过境影像，本次采用的影像日期分别为：7月31日、8月21日、9月15日、9月26日、10月5日、10月30日、11月23日。 由于本文采用的一些AHF反演算法需要借助植被指数，所以进一步选用了2018年的Landsat 8影像。数据来自美国地质调查局（USGS）官方网站（ https：//earthexplorer. https://earthexplorer. ‍usgs.‍gov/［2020-06-22］‍）发布的Level 1T产品（path/row：118—121/40—44）。影像获取时间尽量选取7—8月植被长势最好的夏季，但因福建多云雨，部分影像质量欠佳，因此使用9—10月影像或相近年份的数据进行补充，并利用时间相近的数据对少量的空洞进行填补（ 表1 表1 ）。由于福建地区以常绿植被为主，所以7—10月的植被长势仍较为相近。 行列号 影像获取日期 云量/% 行列号 影像获取日期 云量/% 118/041 2018-10-30 2 119/043 2019-09-22 1 118/042 2018-10-30 5 120/040 2018-10-28 0 2 119/041 2018-10-05 12 120/041 2018-10-28 0 2019-11-09 0 120/042 2018-10-28 0 4 119/042 2018-10-05 27 120/043 2018-10-28 0 2019-09-22 20 120/044 2018-10-28 10 2018-03-11 0 121/041 2019-09-20 0 2016-06-25 13 121/042 2019-09-20 0 119/043 2018-10-05 0 121/043 2019-09-20 0 3　研究方法 1 鉴于Luojia 1-01夜间灯光数据尚无计算AHF的先例，本文提出一套利用该数据反演AHF影像的技术流程（ 图2 图2 ）。 3.1　几何校正 2 所获取的Luojia 1-01夜间灯光影像虽然已进行系统几何纠正，但根据产品信息显示的定位精度，仍存在150—400 m不等的误差，而且不同景的影像其几何偏移量也不同。因此，本文将几何定位较为精确的Landsat 8作为基准影像，结合夜间灯光影像的成像特性，选取具有空间特征的地面控制点，如道路的交叉路口、桥梁的端点、建筑与水体等自然地物的边界折点等，将每一景Luojia 1-01逐一与Landsat 8影像进行配准，误差控制在0.5个像元内。 3.2　辐射校正 2 （1）Luojia 1-01夜间灯光影像的拼接与辐射亮度转换。由于覆盖福建省全境范围的夜间灯光数据需由多景影像进行拼接，并且存在同一区域覆盖多景影像的情况，为了确保夜间灯光数据的空间连续性，并消除夜间灯光数据存在的短暂灯光噪声，将相互重叠的区域进行灯光亮度均值化的处理，同时规定若有一景影像的像元值为0，则该像元所代表的这一位置则为0（ 曹子阳 等，2015 曹子阳 等，2015 ）。将处理后的影像进行拼接，并裁剪得到研究区的夜间灯光影像。进一步将影像的DN值转换为辐射亮度值，转换公式源于Luojia 1-01官方网站（ http：//www.hbeos.org.cn http://www.hbeos.org.cn ［2020-06-22］）： <math display="block" id="M1"><mi>L</mi><mo>=</mo><mi mathvariant="normal">D</mi><msup><mrow><mi mathvariant="normal">N</mi></mrow><mrow><mfrac><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac></mrow></msup><mo>×</mo><msup><mrow><mn mathvariant="normal">10</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">10</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770364&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770362&type= 24.12999916 5.24933338 式中， L 为夜间灯光绝对辐射校正后辐射亮度值（W/（m 2 ·sr·μm））。 （2）Luojia 1-01夜间灯光影像的异常值校正。夜间灯光异常值是指那些具有异常高亮度值的像元，可由火光或油气燃烧产生（ Shi等，2014 Shi等，2014 ）。本文根据 Shi等（2014） Shi等（2014） 提出的异常值校正方法，将福建地区城市发展程度最高的省会福州和经济特区厦门的夜间灯光进行统计，得到其最大辐射亮度值为9814.52 W/（m 2 ·sr·μm），将该值作为异常值校正阈值，若研究区内有高于这一阈值的像元，则查看该像元相邻的8个像元，并选择8个相邻像元中最大的DN值重新赋值给这一异常像元。 （3）Landsat 8的辐射校正则采用Landsat用户手册中的模型和Chavez的COST大气校正模型，将原始影像的DN值转换为传感器处反射率（at-sensor reflectance）（ Chavez，1996 Chavez，1996 ； Chander和Markham，2003 Chander和Markham，2003 ； USGS，2019 USGS，2019 ）。 最后将经过影像配准和辐射校正的Luojia 1-01夜间灯光影像和Landsat 8的投影转换为Alberts等面积投影，并将Luojia 1-01空间分辨率重采样为与Landsat 8相同的30 m（ 图1（b） 图1（b） ）。 3.3　AHF空间估算指数的计算 2 本文选取了当前常用的3种夜间灯光数据AHF估算指数进行实验，分别为NTL nor 、HSI和VANUI。这3种指数都是基于美国前后两代夜间灯光数据DMSP/OLS和NPP/VIIRS提出的，其中HSI和VANUI还需综合MODIS植被指数产品进行计算。目前针对Luojia 1-01夜间灯光数据的AHF估算指数还未见报道。因此，本文尝试将这3个指数移植到Luojia 1-01数据中，并结合Landsat影像的归一化植被指数NDVI（Normalized Difference Vegetation Index）进行AHF的反演。NTL nor 、HSI和VANUI指数的计算公式分别为（ Lu等，2008 Lu等，2008 ； Zhang等，2013 Zhang等，2013 ）： <math display="block" id="M2"><mi mathvariant="normal">N</mi><mi mathvariant="normal">T</mi><msub><mrow><mi mathvariant="normal">L</mi></mrow><mrow><mi mathvariant="normal">n</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">r</mi></mrow></msub><mo>=</mo><mfrac><mrow><mi mathvariant="normal">N</mi><mi mathvariant="normal">T</mi><mi mathvariant="normal">L</mi><mo>-</mo><mi mathvariant="normal">N</mi><mi mathvariant="normal">T</mi><msub><mrow><mi mathvariant="normal">L</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub></mrow><mrow><mi mathvariant="normal">N</mi><mi mathvariant="normal">T</mi><msub><mrow><mi mathvariant="normal">L</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi></mrow></msub><mo>-</mo><mi mathvariant="normal">N</mi><mi mathvariant="normal">T</mi><msub><mrow><mi mathvariant="normal">L</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770368&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770367&type= 38.60800171 8.72066593 <math display="block" id="M3"><mi mathvariant="normal">H</mi><mi mathvariant="normal">S</mi><mi mathvariant="normal">I</mi><mo>=</mo><mfrac><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>-</mo><mi mathvariant="normal">N</mi><mi mathvariant="normal">D</mi><mi mathvariant="normal">V</mi><msub><mrow><mi mathvariant="normal">I</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><mi mathvariant="normal">N</mi><mi mathvariant="normal">T</mi><msub><mrow><mi mathvariant="normal">L</mi></mrow><mrow><mi mathvariant="normal">n</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">r</mi></mrow></msub></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>-</mo><mi mathvariant="normal">N</mi><mi mathvariant="normal">T</mi><msub><mrow><mi mathvariant="normal">L</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><mi mathvariant="normal">N</mi><mi mathvariant="normal">D</mi><mi mathvariant="normal">V</mi><msub><mrow><mi mathvariant="normal">I</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi></mrow></msub><mo>+</mo><mi mathvariant="normal">N</mi><mi mathvariant="normal">T</mi><msub><mrow><mi mathvariant="normal">L</mi></mrow><mrow><mi mathvariant="normal">n</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">r</mi></mrow></msub><mo>×</mo><mi mathvariant="normal">N</mi><mi mathvariant="normal">D</mi><mi mathvariant="normal">V</mi><msub><mrow><mi mathvariant="normal">I</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770372&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770370&type= 73.23667145 9.14400005 <math display="block" id="M4"><mi mathvariant="normal">V</mi><mi mathvariant="normal">A</mi><mi mathvariant="normal">N</mi><mi mathvariant="normal">U</mi><mi mathvariant="normal">I</mi><mo>=</mo><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>-</mo><mi mathvariant="normal">N</mi><mi mathvariant="normal">D</mi><mi mathvariant="normal">V</mi><msub><mrow><mi mathvariant="normal">I</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi></mrow></msub><mo stretchy="false">)</mo><mo>×</mo><mi mathvariant="normal">N</mi><mi mathvariant="normal">T</mi><msub><mrow><mi mathvariant="normal">L</mi></mrow><mrow><mi mathvariant="normal">n</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">r</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770375&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770374&type= 49.86866760 4.57200003 式中，NTL、NTL max 和NTL min 分别为夜间灯光辐射亮度值及其最大值和最小值；NDVI max 为NDVI在研究年份中的最大值。本文利用所获取的各景Landsat 8影像计算NDVI值，然后通过逐景NDVI数据的叠加分析，从中求出NDVI的最大值（ Lu等，2008 Lu等，2008 ），其计算公式为 <math display="block" id="M5"><mi mathvariant="normal">N</mi><mi mathvariant="normal">D</mi><mi mathvariant="normal">V</mi><msub><mrow><mi mathvariant="normal">I</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi></mrow></msub><mo>=</mo><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi><mfenced separators="|"><mrow><mi mathvariant="normal">N</mi><mi mathvariant="normal">D</mi><mi mathvariant="normal">V</mi><msub><mrow><mi mathvariant="normal">I</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>,</mo><mtext> </mtext><mi mathvariant="normal">N</mi><mi mathvariant="normal">D</mi><mi mathvariant="normal">V</mi><msub><mrow><mi mathvariant="normal">I</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mtext> </mtext><mo>,</mo><mo>…</mo><mo>,</mo><mtext> </mtext><mi mathvariant="normal">N</mi><mi mathvariant="normal">D</mi><mi mathvariant="normal">V</mi><msub><mrow><mi mathvariant="normal">I</mi></mrow><mrow><mi>n</mi></mrow></msub></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770379&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770377&type= 64.34667206 5.24933338 式中，NDVI 1 ，NDVI 2 ，…，NDVI n 为各景影像的NDVI数据。 对于NTL min 的获取，本文采用常用的指数阈值法（ Imhoff等，1997 Imhoff等，1997 ； Ma等，2012 Ma等，2012 ； Yue等，2014 Yue等，2014 ）。利用Google Earth高分影像对Landsat 8 NDVI max 数据纯植被像元提取效果的目视比对，确定将NDVI max 大于0.52设定为纯植被像元，其上叠加的灯光辐射亮度值如大于0，即认为这一区域为灯光溢出。据此统计得到灯光溢出区域的夜间灯光辐射亮度均值为2.02 W/（m 2 ·sr·μm），因此将NTL min 设为2.02。 3.4　基于统计年鉴的AHF数据统计计算 2 基于统计年鉴计算的AHF数据是夜间灯光数据反演结果的对比验证依据。因此，本文采用至上而下源清单法，根据福建省统计年鉴数据（ 福建省统计局和国家统计局福建调查总队，2019 福建省统计局和国家统计局福建调查总队，2019 ），对2018年福建省84个县（市、区）的工业、交通、建筑（商业和生活建筑）以及人体新陈代谢（ 马盼盼 等，2016 马盼盼 等，2016 ； Chen和Hu，2017 Chen和Hu，2017 ； Wang等，2019 Wang等，2019 ），这4个主要人为热排放源的年均AHF进行计算，公式为 <math display="block" id="M6"><msub><mrow><mi>Q</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo>=</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi mathvariant="normal">I</mi></mrow></msub><mo>+</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi mathvariant="normal">V</mi></mrow></msub><mo>+</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi mathvariant="normal">B</mi></mrow></msub><mo>+</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi mathvariant="normal">M</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770383&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770381&type= 36.99933624 4.57200003 式中， <math id="M7"><msub><mrow><mi>Q</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770386&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770384&type= 3.64066648 4.48733330 为地区年均AHF总量（W/m 2 ）； Q I 为工业热通量（W/m 2 ）； Q V 为交通热通量（W/m 2 ）； Q B 为建筑热通量（W/m 2 ）； Q M 为人体新陈代谢热通量（W/m 2 ）。 Q I 、 Q V 、 Q B 和 Q M 的计算公式分别为 <math display="block" id="M8"><msub><mrow><mi>Q</mi></mrow><mrow><mi mathvariant="normal">I</mi></mrow></msub><mo>=</mo><mfrac><mrow><msub><mrow><mi>E</mi></mrow><mrow><mi mathvariant="normal">I</mi></mrow></msub><mo>×</mo><mi>C</mi></mrow><mrow><mi>A</mi><mo>×</mo><mi>T</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770390&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770388&type= 18.45733452 7.61999989 式中， E I 为工业标准煤热消耗，用每吨标准煤当量表示（tce）； C 为标准煤热值，中国规定的标准煤发热值为29306 kJ/kg； A 为研究区土地面积（m 2 ）； T 为时间，本次研究为1年。 <math display="block" id="M9"><msub><mrow><mi>Q</mi></mrow><mrow><mi mathvariant="normal">V</mi></mrow></msub><mo>=</mo><mfrac><mrow><mi>C</mi><mi>M</mi><mi>V</mi><mo>×</mo><mi>D</mi><mo>×</mo><mi>E</mi><mo>×</mo><mi>ρ</mi><mo>×</mo><mi>N</mi><mi>H</mi><mi>C</mi></mrow><mrow><mi>A</mi><mo>×</mo><mi>T</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770393&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770391&type= 48.09066391 7.36600018 式中， CMV 为民用机动车保有量； D 为机动车年均行驶距离（km）； E 为汽油消耗率（L/km）； ρ 为汽油消耗密度（kg/L）； NHC 为净热量消耗（kJ/g）（ Chen和Hu，2017 Chen和Hu，2017 ）。 <math display="block" id="M10"><msub><mrow><mi>Q</mi></mrow><mrow><mi mathvariant="normal">B</mi></mrow></msub><mo>=</mo><mfrac><mrow><mo stretchy="false">(</mo><msub><mrow><mi>E</mi></mrow><mrow><mi mathvariant="normal">B</mi><mi mathvariant="normal">C</mi></mrow></msub><mo>+</mo><msub><mrow><mi>E</mi></mrow><mrow><mi mathvariant="normal">B</mi><mi mathvariant="normal">R</mi></mrow></msub><mo stretchy="false">)</mo><mo>×</mo><mi>C</mi></mrow><mrow><mi>A</mi><mo>×</mo><mi>T</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770398&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770396&type= 33.02000046 7.87400007 式中， E BC 和 E BR 分别为商业建筑和生活建筑标准煤热消耗（tce）。 <math display="block" id="M11"><msub><mrow><mi>Q</mi></mrow><mrow><mi mathvariant="bold-italic">M</mi></mrow></msub><mo>=</mo><mfrac><mrow><mo stretchy="false">(</mo><msub><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>×</mo><msub><mrow><mi>t</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>+</mo><msub><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>×</mo><msub><mrow><mi>t</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">)</mo><mo>×</mo><mi>P</mi></mrow><mrow><mo stretchy="false">(</mo><msub><mrow><mi>t</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>+</mo><msub><mrow><mi>t</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">)</mo><mo>×</mo><mi>A</mi><mo>×</mo><mi>T</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770402&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32770399&type= 41.82533264 9.14400005 式中， s 1 和 s 2 分别为人体活动阶段（7：00—23：00）和人体睡眠阶段（23：00—7：00）的人均代谢热消耗（W）（ Grimmond，1992 Grimmond，1992 ； Quah和Roth，2012 Quah和Roth，2012 ）； t 1 和 t 2 分别为人体活动时长和睡眠时长（h）； P 为人口数量。 3.5　AHF遥感估算模型的构建与验证 2 前文述及，夜间灯光亮度值与AHF有着显著的相关性，利用这一特性可将夜间灯光估算指数结果与统计年鉴数据的AHF计算结果进行关联，构建基于夜间灯光数据的AHF遥感空间估算模型。根据统计年鉴，通过式（ 6 6 ）—（ 10 10 ）可计算得到福建省84个县域（不含金门县）的年均AHF统计值（记为：AHF sta ），并将其作为因变量（ y ）；将NTL nor 、HSI和VANUI估算指数计算得到的各县域均值作为自变量（ x ），然后进行多种函数的回归分析，建立AHF的遥感估算模型。 通过对比不同估算模型的回归方程拟合度（ R 2 ），初步筛选出 R 2 较高的估算模型。为了进一步评价AHF估算模型的结果，采用统计学上常用的五折交叉验证法（5-fold cross-validation），该交叉验证是将样本数据随机分成5份，将其中4份作为训练数据，1份作为测试数据，最终得到5次验证结果，并将其作为对估算模型精度的结果，使用均方根误差（RMSE）进行定量表示（ 杨立娟 等，2018 杨立娟 等，2018 ； Yang等，2019 Yang等，2019 ）。 4　结果与讨论 1 4.1　AHF最佳遥感估算模型的确定与验证 2 为了从NTL nor 、HSI和VANUI这3个指数中确定出最佳估算模型，在回归时分别采用线性、多项式、指数、对数、乘幂等5种函数对3个指数（ x ）和基于统计年鉴计算的AHF（记为：AHF sta ）进行回归分析。各估算模型的方程拟合都通过了1%的显著性检验（ p < 0.01），从 表2 表2 可以看出，乘幂函数估算模型在3个指数中的 R 2 都是最高的，而其中NTL nor 、VANUI的 R 2 都大于0.9，高于HSI的0.8。 估算指数 函数 估算模型 R 2 5 NTL nor 线性 AHF sta =2102.8332 x +0.7801 0.6234 二次多项式 AHF sta =-51032.11 x 2 +4286.88 x -1.7702 0.7020 指数 AHF sta =0.6776 e 141.8742 x 0.4950 对数 AHF sta =6.9995 ln（ x ）+62.3435 0.3524 乘幂 AHF sta =822.8343 x 0.8482 0.9027 5 HSI 线性 AHF sta =126.6906 x -47.3019 0.4417 二次多项式 AHF sta =396.48 x 2 -294.87 x +54.9677 0.5194 指数 AHF sta =0.0041 e 12.8702 x 0.7952 对数 AHF sta =59.9527 ln（ x ）+59.6828 0.3916 乘幂 AHF sta =318.9966 x 6.5477 0.8149 5 VANUI 线性 AHF sta =2520.4561 x +0.8851 0.6233 二次多项式 AHF sta =-75973.01 x 2 +5215.28 x -1.6148 0.7007 指数 AHF sta =0.6835 e 169.4279 x 0.4914 对数 AHF sta =6.7345 ln（ x ）+62.3494 0.3493 乘幂 AHF sta =848.2905 x 0.8198 0.9028 虽然VANUI乘幂估算模型的 R 2 略高于NTL nor 乘幂估算模型，但优势不明显。因此，将二者乘幂函数估算模型进行五折交叉验证。从 图3 图3 的验证结果可知，2个估算模型的估算值与统计的AHF sta 验证值之间的拟合度很高， R 2 均达到0.98。就RMSE误差而言，VANUI乘幂估算模型的2.1786 W/m 2 小于NTL nor 的2.2411 W/m 2 ，因此本文选用VANUI乘幂估算模型作为福建省AHF的最佳遥感估算模型。 1,2; 从上述NTL nor 、HSI和VANUI的估算模型构建结果可以看出，这3个指数与AHF sta 的回归方程都有着较好的拟合度，其中NTL nor 和VANUI乘幂估算模型验证有着较高的精度。但VANUI从估算模型拟合度和模型估算精度方面都要略高于NTL nor ，说明在单一夜间灯光数据得到的NTL nor 中引入植被指数，对模型估算精度的提高确实有一定作用。然而，对于HSI，虽然也同样引入植被指数，但估算模型拟合效果仍不够理想，通过考察HSI的计算公式（ 式（3） 式（3） ）可以发现，由于NTL nor 与NDVI max 有着较显著的乘幂关系，因此当NTL nor 在高亮灯光区域出现最大值，或接近最大值1时，NDVI max 则将接近最小值0，这将造成HSI值呈指数级快速增加，出现极大值的极端情况，导致HSI计算结果的失真（ Zhang等，2013 Zhang等，2013 ）。 Wang等（2019） Wang等（2019） 也在研究中指出，运用NTL nor 和VANUI所构建的估算模型其AHF估算结果要明显优于HSI。 4.2　AHF空间估算结果 2 根据基于VANUI的乘幂估算模型（AHF sta =848.2905 x 0.8198 ）反演出2018年福建省AHF分布图（ 图4（a） 图4（a） ）和各县（市、区）的AHF均值分布图（ 图4（b） 图4（b） ）。可以看出，福建省AHF整体上呈现明显的空间分异性，经济水平较高的东南沿海地区，其AHF也较高，在空间上形成了从福州到厦门的沿海人为热排放集聚带；而中、西部地区多为山地，森林覆盖度高，城市较为分散，因此AHF整体较低，分布也较为稀疏。 对 图4 图4 进行统计可知，2018年全省AHF均值为0.88 W/m 2 。9个地级市中，厦门的年均AHF最高，达10.98 W/m 2 ，泉州、莆田、福州、漳州等沿海城市次之，年均AHF在0.98—1.95 W/m 2 ，而北部的宁德以及中西部内陆山区的龙岩、三明、南平等城市则较低，均值在0.38—0.46 W/m 2 （ 图5 图5 ），在空间上明显体现出东西分异的状况。 1,2; 为探究Luojia 1-01高分辨率数据在确定城市内部AHF空间细节方面是否有优势，对福建省东部沿海AHF高值集聚最明显的福州市区、泉州市区和厦门岛进行城市AHF空间细节呈现（ 图6 图6 ）。可以看出，Luojia 1-01影像具有较强的城市内部空间刻画能力，能够清晰地揭示城区人为热排放的分异细节。进一步对Luojia 1-01所呈现的AHF空间细节进行分析，在城市空间尺度上，3个城市的AHF均呈现在老城区或大型商业区的块状集聚，如福州东街口和江滨商业区、泉州城市中心区、厦门的老城区和集中商业区等区域，这些区域的AHF多在150 W/m 2 以上，局部最高值区域可超过500 W/m 2 。此外，一些大型的市政公共设施区也有较大面积的高值集聚，如厦门高崎国际机场、货运码头、厦门国际会展中心的AHF都超过100 W/m 2 。相比于商业区和大型公共设施区，一些连片居住片区，如福州金山住宅区、泉州东海湾片住宅区、厦门岛内住宅区等，其AHF则有一定程度的降低，可保持在10—50 W/m 2 。而一些城市近郊的村庄和社区的AHF则更低，一般为2—10 W/m 2 。从 图6 图6 还可以发现，人为热排放除了呈块状集聚外，同时沿城市主干道呈网状蔓延，根据主干道的车流情况与繁华程度的不同，AHF呈现出不同的数值，一般为30—100 W/m 2 ，而在一些道路的交汇点，AHF可达200 W/m 2 （ 表3 表3 ）。 1. 福州东街口商业区；2. 江滨商业区；3. 金山住宅区；4. 近郊村庄；5. 泉州商业区；6. 东海湾住宅区；7. 近郊社区；8. 厦门老城区和商业区；9. 住宅区；10. 厦门国际机场；11. 货运码头；12. 厦门国际会展中心 1,2,3;4; 用地类型 AHF均值/（W/m 2 ） 城市集中商业区 169.081 大型市政公共设施区 136.836 城市主干道 35.937 城市住宅区 23.644 近郊住宅区 8.264 以上分析表明Luojia 1-01高空间分辨率数据有助于城市内部AHF空间细节的辨析，并可获取各土地覆盖类型的AHF定量数值。根据用地属性和用地功能的不同，AHF数值可表现为：城市集中商业区>大型市政公共设施区>城市主干道>城市住宅区>近郊住宅区，这一发现将有助于人为热的减排和治理。但是，研究也发现作为城市白天关键热源之一的工业区未能在AHF结果中很好地体现，许多工业区在反演的AHF影像中的值并不高。这主要是由于工业区的生产活动及其热排放主要是在白天进行，因此夜间灯光数据未能捕捉到夜间工业区的信息，造成了工业区AHF数值偏低的结果。在今后的AHF研究中可加入白天过空的热红外遥感数据，以弥补这一缺陷。这是由于热红外数据可以较全面地揭示日间人为活动热源的空间分异规律，特别是对工业区的生产及其热排放有明显的热辐射信息反馈，因此，在后续研究中可尝试耦合夜间灯光与热红外数据，以期更加全面、精准地定量估算AHF的空间分布。 5　结论 1 本文运用Luojia 1-01夜间灯光数据，辅以Landsat 8植被指数数据，构建了AHF空间估算的遥感模型，并反演得到了福建省AHF空间分布影像，所获得主要结论如下： （1）Luojia 1-01夜间灯光可以有效地反演AHF数据，揭示其空间分布细节。在所采用的3个指数中，VANUI指数与福建省2018年各县域的AHF统计数据的拟合度最高，基于其乘幂函数构建的AHF遥感估算模型的RMSE误差也最小。 （2）从VANUI估算模型反演得到的AHF结果可知，2018年福建省年均AHF为0.88 W/m 2 ，且在空间上表现为经济水平较高的东南沿海地区的AHF明显高于中、西部内陆地区，并形成了东部沿海热排放集聚带。由于AHF与工业产值和城市人口密度的相关度最高，因此，在今后的城市发展中，应避免人口和工业过度集中，逐步对高热集聚区进行有序疏散。 （3）Luojia 1-01相比于前两代夜间灯光DMSP/OLS和Suomi-NPP/VIIRS在影像的空间分辨率上有了较大的提升，这显著提高了夜间灯光数据对AHF空间细节的辨析能力。本文也因此发现了研究区城市内部因用地属性和功能的不同而表现出的AHF不同，具体表现为：城市集中商业区>大型市政公共设施区>城市主干道>城市住宅区>郊区住宅区，这一发现将有助于城市人为热的减排和治理。 参考文献（References） Block A , Keuler K and Schaller E . 2004 . Impacts of anthropogenic heat on regional climate patterns . Geophysical Research Letters , 31 ( 12 ): L 12211 [ DOI: 10.1029/2004GL019852 http://dx.doi.org/10.1029/2004GL019852 ] Impacts of anthropogenic heat on regional climate patterns Bohnenstengel S I , Hamilton I , Davies M and Belcher S E . 2014 . Impact of anthropogenic heat emissions on London's temperatures . Quarterly Journal of the Royal Meteorological Society , 140 ( 679 ): 687 - 698 [ DOI: 10.1002/qj.2144 http://dx.doi.org/10.1002/qj.2144 ] Impact of anthropogenic heat emissions on London's temperatures Cao Z Y , Wu Z F , Kuang Y Q and Huang N S . 2015 . Correction of DMSP/OLS night-time light images and its application in China . Journal of Geo-information Science , 17 ( 9 ): 1092 - 1102 Correction of DMSP/OLS night-time light images and its application in China 曹子阳 , 吴志峰 , 匡耀求 , 黄宁生 . 2015 . DMSP/OLS夜间灯光影像中国区域的校正及应用 . 地球信息科学学报 , 17 ( 9 ): 1092 - 1102 [ DOI: 10.3724/SP.J.1047.2015.01092 http://dx.doi.org/10.3724/SP.J.1047.2015.01092 ] DMSP/OLS夜间灯光影像中国区域的校正及应用 Chander G and Markham B . 2003 . Revised Landsat-5 TM radiometric calibration procedures and postcalibration dynamic ranges . IEEE Transactions on Geoscience and Remote Sensing , 41 ( 11 ): 2674 - 2677 [ DOI: 10.1109/TGRS.2003.818464 http://dx.doi.org/10.1109/TGRS.2003.818464 ] Revised Landsat-5 TM radiometric calibration procedures and postcalibration dynamic ranges Chavez P S . 1996 . Image-based atmospheric corrections revisited and improved . Photogrammetric Engineering and Remote Sensing , 62 ( 9 ): 1025 - 1035 Image-based atmospheric corrections revisited and improved Chen B , Chen L F , Dong L and Shi G Y . 2016 . Estimating the global distribution of anthropogenic heat release and exploring its possible climatic effect . Chinese Journal of Atmospheric Sciences , 40 ( 2 ): 289 - 295 Estimating the global distribution of anthropogenic heat release and exploring its possible climatic effect 陈兵 , 陈良富 , 董理 , 石广玉 . 2016 . 人为热释放: 全球分布的估算及其气候效应的探索 . 大气科学 , 40 ( 2 ): 289 - 295 [ DOI: 10.3878/j.issn.1006-9895.1504.14268 http://dx.doi.org/10.3878/j.issn.1006-9895.1504.14268 ] 人为热释放: 全球分布的估算及其气候效应的探索 Chen B , Shi G Y , Wang B , Zhao J Q and Tan S C . 2012 . Estimation of the anthropogenic heat release distribution in China from 1992 to 2009 . Acta Meteorologica Sinica , 26 ( 4 ): 507 - 515 [ DOI: 10.1007/s13351-012-0409-y http://dx.doi.org/10.1007/s13351-012-0409-y ] Estimation of the anthropogenic heat release distribution in China from 1992 to 2009 Chen S S and Hu D Y . 2017 . Parameterizing anthropogenic heat flux with an energy-consumption inventory and multi-source remote sensing data . Remote Sensing , 9 ( 11 ): 1165 [ DOI: 10.3390/rs9111165 http://dx.doi.org/10.3390/rs9111165 ] Parameterizing anthropogenic heat flux with an energy-consumption inventory and multi-source remote sensing data Chen Y B , Zheng Z H , Wu Z F and Qian Q L . 2019 . Review and prospect of application of nighttime light remote sensing data . Progress in Geography , 38 ( 2 ): 205 - 223 Review and prospect of application of nighttime light remote sensing data 陈颖彪 , 郑子豪 , 吴志峰 , 千庆兰 . 2019 . 夜间灯光遥感数据应用综述和展望 . 地理科学进展 , 38 ( 2 ): 205 - 223 [ DOI: 10.18306/dlkxjz.2019.02.005 http://dx.doi.org/10.18306/dlkxjz.2019.02.005 ] 夜间灯光遥感数据应用综述和展望 Elvidge C D , Baugh K E , Kihn E A , Kroehl H W and Davis E R . 1997 . Mapping city lights with nighttime data from the DMSP Operational Linescan System . Photogrammetric Engineering and Remote Sensing , 63 ( 6 ): 727 - 734 Mapping city lights with nighttime data from the DMSP Operational Linescan System Flanner M G . 2009 . Integrating anthropogenic heat flux with global climate models . Geophysical Research Letters , 36 ( 2 ): L 02801 [ DOI: 10.1029/2008GL036465 http://dx.doi.org/10.1029/2008GL036465 ] Integrating anthropogenic heat flux with global climate models Fujian Province Bureau of Statistics and Fujian Survey Team of National Bureau of Statistics . 2019 . Fujian Statistical Yearbook-2019 . Beijing : China Statistics Press Fujian Statistical Yearbook-2019 福建省统计局, 国家统计局福建调查总队 . 2019 . 福建统计年鉴—2019 . 北京 : 中国统计出版社 福建统计年鉴—2019 Grimmond C S B . 1992 . The suburban energy balance: methodological considerations and results for a mid-latitude west coast city under winter and spring conditions . International Journal of Climatology , 12 ( 5 ): 481 - 497 [ DOI: 10.1002/joc.3370120506 http://dx.doi.org/10.1002/joc.3370120506 ] The suburban energy balance: methodological considerations and results for a mid-latitude west coast city under winter and spring conditions He X F , Jiang W M , Chen Y and Liu G . 2007 . Numerical simulation of the impacts of anthropogenic heat on the structure of the urban boundary layer . Chinese Journal of Geophysics , 50 ( 1 ): 74 - 82 Numerical simulation of the impacts of anthropogenic heat on the structure of the urban boundary layer 何晓凤 , 蒋维楣 , 陈燕 , 刘罡 . 2007 . 人为热源对城市边界层结构影响的数值模拟研究 . 地球物理学报 , 50 ( 1 ): 74 - 82 [ DOI: 10.3321/j.issn:0001-5733.2007.01.010 http://dx.doi.org/10.3321/j.issn:0001-5733.2007.01.010 ] 人为热源对城市边界层结构影响的数值模拟研究 Iamarino M , Beevers S and Grimmond C S B . 2012 . High-resolution (space, time) anthropogenic heat emissions: London 1970 - 2025 . International Journal of Climatology , 32 ( 11 ): 1754 - 1767 [ DOI: 10.1002/joc.2390 http://dx.doi.org/10.1002/joc.2390 ] High-resolution (space, time) anthropogenic heat emissions: London Ichinose T , Shimodozono K and Hanaki K . 1999 . Impact of anthropogenic heat on urban climate in Tokyo . Atmospheric Environment , 33 ( 24/25 ): 3897 - 3909 [ DOI: 10.1016/S1352-2310(99)00132-6 http://dx.doi.org/10.1016/S1352-2310(99)00132-6 ] Impact of anthropogenic heat on urban climate in Tokyo Imhoff M L , Lawrence W T , Stutzer D C and Elvidge C D . 1997 . A technique for using composite DMSP/OLS “city lights” satellite data to map urban area . Remote Sensing of Environment , 61 ( 3 ): 361 - 370 [ DOI: 10.1016/S0034-4257(97)00046-1 http://dx.doi.org/10.1016/S0034-4257(97)00046-1 ] A technique for using composite DMSP/OLS “city lights” satellite data to map urban area Imhoff M L , Zhang P , Wolfe R E and Bounoua L . 2010 . Remote sensing of the urban heat island effect across biomes in the continental USA . Remote Sensing of Environment , 114 ( 3 ): 504 - 513 [ DOI: 10.1016/j.rse.2009.10.008 http://dx.doi.org/10.1016/j.rse.2009.10.008 ] Remote sensing of the urban heat island effect across biomes in the continental USA Jiang W , He G J , Long T F , Guo H X , Yin R Y , Leng W C , Liu H C and Wang G Z . 2018 . Potentiality of using Luojia 1 - 01 nighttime light imagery to investigate artificial light pollution. Sensors , 18 ( 9 ): 2900 [ DOI: 10.3390/s18092900 http://dx.doi.org/10.3390/s18092900 ] Potentiality of using Luojia Li D R , Zhang G , Shen X , Zhong X , Jiang Y H , Wang T Y , Tu J G and Li Z J . 2019 . Design and processing night light remote sensing of LJ-1 01 satellite . Journal of Remote Sensing , 23 ( 6 ): 1011 - 1022 Design and processing night light remote sensing of LJ-1 01 satellite 李德仁 , 张过 , 沈欣 , 钟兴 , 蒋永华 , 汪韬阳 , 涂建光 , 李治江 . 2019 . 珞珈一号01星夜光遥感设计与处理 . 遥感学报 , 23 ( 6 ): 1011 - 1022 [ DOI: 10.11834/jrs.20199327 http://dx.doi.org/10.11834/jrs.20199327 ] 珞珈一号01星夜光遥感设计与处理 Lu D S , Tian H Q , Zhou G M and Ge H L . 2008 . Regional mapping of human settlements in southeastern China with multisensor remotely sensed data . Remote Sensing of Environment , 112 ( 9 ): 3668 - 3679 [ DOI: 10.1016/j.rse.2008.05.009 http://dx.doi.org/10.1016/j.rse.2008.05.009 ] Regional mapping of human settlements in southeastern China with multisensor remotely sensed data Ma P P , Wu J J , Yang X C and Qi J G . 2016 . Spatialization of anthropogenic heat using multi-sensor remote sensing data: a case study of Zhejiang Province, East China . China Environmental Scienc , 36 ( 1 ): 314 - 320 Spatialization of anthropogenic heat using multi-sensor remote sensing data: a case study of Zhejiang Province, East China 马盼盼 , 吾娟佳 , 杨续超 , 齐家国 . 2016 . 基于多源遥感信息的人为热排放量空间化——以浙江省为例 . 中国环境科学 , 36 ( 1 ): 314 - 320 [ DOI: 10.3969/j.issn.1000-6923.2016.01.052 http://dx.doi.org/10.3969/j.issn.1000-6923.2016.01.052 ] 基于多源遥感信息的人为热排放量空间化——以浙江省为例 Ma T , Zhou C H , Pei T , Haynie S and Fan J F . 2012 . Quantitative estimation of urbanization dynamics using time series of DMSP/OLS nighttime light data: a comparative case study from China’s cities . Remote Sensing of Environment , 124 : 99 - 107 [ DOI: 10.1016/j.rse.2012.04.018 http://dx.doi.org/10.1016/j.rse.2012.04.018 ] Quantitative estimation of urbanization dynamics using time series of DMSP/OLS nighttime light data: a comparative case study from China’s cities Oke T R . 1976 . The distinction between canopy and boundary-layer urban heat islands . Atmosphere , 14 ( 4 ): 268 - 277 [ DOI: 10.1080/00046973.1976.9648422 http://dx.doi.org/10.1080/00046973.1976.9648422 ] The distinction between canopy and boundary-layer urban heat islands Pal S , Xueref-Remy I , Ammoura L , Chazette P , Gibert F , Royer P , Dieudonné E , Dupont J C , Haeffelin M , Lac C , Lopez M , Morille Y and Ravetta F . 2012 . Spatio-temporal variability of the atmospheric boundary layer depth over the Paris agglomeration: an assessment of the impact of the urban heat island intensity . Atmospheric Environment , 63 : 261 - 275 [ DOI: 10.1016/j.atmosenv.2012.09.046 http://dx.doi.org/10.1016/j.atmosenv.2012.09.046 ] Spatio-temporal variability of the atmospheric boundary layer depth over the Paris agglomeration: an assessment of the impact of the urban heat island intensity Quah A K L and Roth M . 2012 . Diurnal and weekly variation of anthropogenic heat emissions in a tropical city, Singapore . Atmospheric Environment , 46 : 92 - 103 [ DOI: 10.1016/j.atmosenv.2011.10.015 http://dx.doi.org/10.1016/j.atmosenv.2011.10.015 ] Diurnal and weekly variation of anthropogenic heat emissions in a tropical city, Singapore Rizwan A M , Dennis L Y C and Liu C . 2008 . A review on the generation, determination and mitigation of Urban Heat Island . Journal of Environmental Sciences , 20 ( 1 ): 120 - 128 [ DOI: 10.1016/S1001-0742(08)60019-4 http://dx.doi.org/10.1016/S1001-0742(08)60019-4 ] A review on the generation, determination and mitigation of Urban Heat Island Shi K F , Yu B L , Huang Y X , Hu Y J , Yin B , Chen Z Q , Chen L J and Wu J P . 2014 . Evaluating the ability of NPP-VIIRS nighttime light data to estimate the gross domestic product and the electric power consumption of China at multiple scales: a comparison with DMSP-OLS data . Remote Sensing , 6 ( 2 ): 1705 - 1724 [ DOI: 10.3390/rs6021705 http://dx.doi.org/10.3390/rs6021705 ] Evaluating the ability of NPP-VIIRS nighttime light data to estimate the gross domestic product and the electric power consumption of China at multiple scales: a comparison with DMSP-OLS data Tong H , Liu H Z , Sang J G and Hu F . 2004 . The impact of urban anthropogenic heat on Beijing heat environment . Climatic and Environmental Research , 9 ( 3 ): 409 - 421 The impact of urban anthropogenic heat on Beijing heat environment 佟华 , 刘辉志 , 桑建国 , 胡非 . 2004 . 城市人为热对北京热环境的影响 . 气候与环境研究 , 9 ( 3 ): 409 - 421 [ DOI: 10.3969/j.issn.1006-9585.2004.03.001 http://dx.doi.org/10.3969/j.issn.1006-9585.2004.03.001 ] 城市人为热对北京热环境的影响 Torrance K E and Shun J S W . 1976 . Time-varying energy consumption as a factor in urban climate . Atmospheric Environment ( 1967 ), 10 ( 4 ): 329 - 337 [ DOI: 10.1016/0004-6981(76)90174-8 http://dx.doi.org/10.1016/0004-6981(76)90174-8 ] Time-varying energy consumption as a factor in urban climate United Nations . 2019 . The 2018 revision of World urbanization prospects [EB/OL]. [ 2020-06-22 ]. https://population.un.org/wup/Publications/Files/ WUP2018-Report.pdf https://population.un.org/wup/Publications/Files/WUP2018-Report.pdf The 2018 revision of World urbanization prospects USGS . 2019 . Landsat 8 data users handbook (Version 2.0) [EB/OL]. [ 2020-06-22 ]. https://www.usgs.gov/media/files/landsat-8-data-users-handbook https://www.usgs.gov/media/files/landsat-8-data-users-handbook Landsat 8 data users handbook (Version 2.0) Wang C X , Chen Z Q , Yang C S , Li Q X , Wu Q S , Wu J P , Zhang G and Yu B L . 2020 . Analyzing parcel-level relationships between Luojia 1 - 01 nighttime light intensity and artificial surface features across Shanghai, China: a comparison with NPP-VIIRS data. International Journal of Applied Earth Observation and Geoinformation , 85 : 101989 [ DOI: 10.1016/j.jag.2019.101989 http://dx.doi.org/10.1016/j.jag.2019.101989 ] Analyzing parcel-level relationships between Luojia Wang S S , Hu D Y , Chen S S and Yu C . 2019 . A partition modeling for anthropogenic heat flux mapping in China . Remote Sensing , 11 ( 9 ): 1132 [ DOI: 10.3390/rs11091132 http://dx.doi.org/10.3390/rs11091132 ] A partition modeling for anthropogenic heat flux mapping in China Yang L J , Xu H Q and Jin Z F . 2018 . Estimation of ground-level PM 2.5 concentrations using MODIS satellite data in Fuzhou, China . Journal of Remote Sensing , 22 ( 1 ): 64 - 75 Estimation of ground-level PM2.5 concentrations using MODIS satellite data in Fuzhou, China 杨立娟 , 徐涵秋 , 金致凡 . 2018 . MODIS卫星遥感估计福州地区近地面PM 2.5 浓度 . 遥感学报 , 22 ( 1 ): 64 - 75 [ DOI: 10.11834/jrs.20186501 http://dx.doi.org/10.11834/jrs.20186501 ] MODIS卫星遥感估计福州地区近地面PM2.5浓度 Yang L J , Xu H Q and Jin Z F . 2019 . Estimating ground-level PM 2.5 over a coastal region of China using satellite AOD and a combined model . Journal of Cleaner Production , 227 : 472 - 482 [ DOI: 10.1016/j.jclepro.2019.04.231 http://dx.doi.org/10.1016/j.jclepro.2019.04.231 ] Estimating ground-level PM2.5 over a coastal region of China using satellite AOD and a combined model Yue W Z , Gao J B and Yang X C . 2014 . Estimation of gross domestic product using multi-sensor remote sensing data: a case study in Zhejiang Province, East China . Remote Sensing , 6 ( 8 ): 7260 - 7275 [ DOI: 10.3390/rs6087260 http://dx.doi.org/10.3390/rs6087260 ] Estimation of gross domestic product using multi-sensor remote sensing data: a case study in Zhejiang Province, East China Zhang Q L , Schaaf C and Seto K C . 2013 . The vegetation adjusted NTL Urban Index: a new approach to reduce saturation and increase variation in nighttime luminosity . Remote Sensing of Environment , 129 : 32 - 41 [ DOI: 10.1016/j.rse.2012.10.022 http://dx.doi.org/10.1016/j.rse.2012.10.022 ] The vegetation adjusted NTL Urban Index: a new approach to reduce saturation and increase variation in nighttime luminosity Zhong L and Liu X S . 2019 . Application potential analysis of L J1 - 01 new nighttime light data. Bulletin of Surveying and Mapping , ( 7 ): 132 - 137 Application potential analysis of L 钟亮 , 刘小生 . 2019 . 珞珈一号新型夜间灯光数据应用潜力分析 . 测绘通报 , ( 7 ): 132 - 137 [ DOI: 10.13474/j.cnki.11-2246.2019.0235 http://dx.doi.org/10.13474/j.cnki.11-2246.2019.0235 ] 珞珈一号新型夜间灯光数据应用潜力分析 10.11834/jrs.20210295 1007-4619（2022）06-1236-11 林中立，徐涵秋，林从华.2022.珞珈一号夜间灯光数据的福建省人为热通量估算.遥感学报，26（6）： 1236-1246 林中立，徐涵秋，林从华.2022.珞珈一号夜间灯光数据的福建省人为热通量估算.遥感学报，26（6）： 1236-1246 Lin Z L，Xu H Q and Lin C H. 2022. Estimation of anthropogenic heat flux of Fujian Province （China） based on Luojia 1-01 nighttime light data. National Remote Sensing Bulletin， 26（6）：1236-1246 Lin Z L，Xu H Q and Lin C H. 2022. Estimation of anthropogenic heat flux of Fujian Province （China） based on Luojia 1-01 nighttime light data. National Remote Sensing Bulletin， 26（6）：1236-1246 国家自然科学基金（编号：31971639）；福建省自然科学基金（编号：2020J05193）；福建省教育厅中青年教师教育科研项目（编号：JAT190403）；福建工程学院本科教学改革研究项目（编号：2022JG021） Supported by National Natural Science Foundation of China (No. 31971639); Natural Science Foundation of Fujian Province, China (No. 2020J05193); Undergraduate Teaching Reform Research Project of Fujian University of Technology (No. 2022JG021) 2020-07-22 2022-06-07 2022-08-04 2021-01-05 遥感学报 6 26