卫星遥感夜光辐亮度的观测角度归一化方法技术

本发明专利技术涉及一种卫星遥感夜光辐亮度的观测角度归一化方法。首先建立夜光辐亮度与观测天顶角之间的天顶

【技术实现步骤摘要】
卫星遥感夜光辐亮度的观测角度归一化方法


[0001]本专利技术属于统计回归模型在夜光遥感影像各向异性消除领域的应用，特别是涉及一种卫星遥感夜光辐亮度的观测角度归一化方法。

技术介绍

[0002]夜光遥感是获取夜间无云时地表发射的可见光信息，通过夜间灯光可以发现人类聚集区并记录地表灯光强度信息，进而可以通过夜光遥感影像间接反映居民在城镇中的行为活动与空间分布。深入挖掘夜光辐亮度时间序列的动态变化规律，对于准确监测城市化进程、电力供应以及评估自然灾害等重大突发事件有着重要的作用。
[0003]随着遥感技术的不断发展，夜光遥感产品逐渐丰富。特别是新一代夜光辐亮度数据NPP/VIIRS(National Polar
‑
Orbiting Partnership's Visible Infrared Imaging Radiometer Suite)不仅具有检测微弱夜间灯光的能力，而且弥补了DMSP/OLS数据低空间分辨率与低辐射分辨率的缺点，能更多地应用在建立夜光辐亮度与人类活动模式经验关系的研究中。由于夜间灯光信号会受到杂散光、大气成分、季节性等噪声信号干扰，从而造成夜光辐亮度时间序列的波动。为此，NASA团队生产出了一套全球可用的、包含丰富时间与空间信息的黑色大理石夜光影像产品(VNP46)。该产品去除了受云污染的像素，并对受大气、地形、植被、雪、月球以及杂散光影响的夜光辐亮度进行了校正，使得产品的夜光辐亮度更为可靠。但是Li等学者发现，夜光辐亮度与卫星传感器观测天顶角存在非线性关系，特别是在高建筑物区域，受到城市三维形态的影响，夜光辐亮度会存在明显的冷点和热点效应。因此，在利用长时间序列的夜光影像进行研究时，为了更准确地反映城市的动态变化，消除夜光辐亮度的角度效应十分必要。
[0004]目前许多学者在光学各向异性的研究主要包括基于双向反射率分布函数(BRDF)表征地球表面反射率各向异性的研究与模拟城市上空光分布的光污染研究。地球表面反射光会受到地物的物理结构、传感器视角以及太阳光照角度变化的影响，产生地球表面反射率各向异性特征。Landsat卫星通过最低点一定范围的视角获取地球表面影像，会导致地表反射率产生较小的方向性影响。Roy等人基于全球高质量的MODIS BRDF产品，得出一组固定的BRDF光谱模型参数，从而完成Landsat数据地表反射率角度校正，但作为描述叶子空间分布格局的关键参数，聚集指数(CI)会受到BRDF模型与太阳天顶角的影响。Wei等人使用归一化冷热点差异(NDHD)评估BRDF模型的不同参数配置与太阳天顶角对CI估计的影响，结果表明太阳天顶角为0
°
时会低估CI值，而60
°
时会高估CI值。此外，Jiao等人在线性核驱动的BRDF模型基础上，提出了一种表征纯雪反射率各向异性的雪核，从而更好地模拟雪的散特性。
[0005]在主动光源各向异性的研究领域，Li等人利用卫星传感器提供的多角度观测信息，探究卫星传感器观测角度与夜光辐亮度之间的关系，结果表明，在城市区域二次模型可以较好地描述夜光辐亮度与观测天顶角之间的关系，并且城市内部空间结构会影响角度与辐亮度之间的关系。由于受到观测角度效应的影响，夜光辐亮度会造成错误检测，而无法反
映灯光的实际变化。Elvidge等人针对夜光辐亮度随观测天顶角变化的四种趋势(平稳、凸起、凹陷、最低观测点峰值)，采用最接近观测天顶角0
°
簇内夜光辐亮度中位数与当前校正角度簇内夜光辐亮度中位数的比值作为系数，分别乘以校正角度簇内的夜光辐亮度观测值从而完成夜光辐亮度角度效应的归一化，进而有益于电力可靠性的时空模式挖掘。在光污染研究中，Tong等人基于无暮光与月光的VIIRS夜光辐亮度数据绘制了北欧和北非向上人造光的角度分布，借此更好地了解光污染在夜间环境中的传播机制。Kocifaj等人利用观测角度模拟夜间环境的人造光，并引入纠正功能改进城市发射函数，从而提高了夜晚城市上空光分布的估算结果。
[0006]综合前人研究，目前关于角度各向异性的研究主要针对地表反射率的遥感产品，很少考虑由主动光源生成的夜光影像角度各向异性，并且目前关于夜光遥感的研究多数聚集在城市化进程与社会经济参数估算等领域的应用，忽略了数据本身的质量存在不确定性与无效性，从而导致分析结果的偏差。

技术实现思路

[0007]本专利技术针对现有技术的不足，提供一种卫星遥感夜光辐亮度的观测角度归一化方法。首先通过建立夜光辐亮度与观测天顶角之间的天顶
‑
辐亮度二次(ZRQ)模型，采用最小二乘多项式曲线拟合方式求解二次函数的系数，然后基于ZRQ模型，将观测到的夜光辐亮度分解为天顶角为0
°
时的夜光辐亮度与二次函数值的乘积，最后采用最小二乘法求解分解后的二次函数系数，并计算当前传感器获取的夜光辐亮度与二次函数值的比值，求解天顶角为0
°
时的夜光辐亮度，进而消除夜光辐亮度的角度效应。
[0008]为了达到上述目的，本专利技术提供的技术方案是一种卫星遥感夜光辐亮度的观测角度归一化方法，包括以下步骤：
[0009]步骤1，获取研究区域黑色大理石(Black Marble)夜光影像产品中VNP46A1与VNP46A2数据。
[0010]步骤2，对步骤1获取到的所有VNP46A1与VNP46A2数据进行数据格式转换，并得到需要的质量文件与影像文件数据。
[0011]步骤3，对步骤2所得到的质量标志数据，设置一定的筛选条件，得到符合需求的高质量夜光辐亮度数据。
[0012]步骤4，以像素为研究单位，从步骤3得到的高质量夜光辐亮度数据中提取像素的夜光辐亮度时间序列和传感器观测天顶角时间序列。
[0013]步骤5，对步骤4提取的夜光辐亮度与卫星传感器观测天顶角进行统计建模，确立两者之间的量化关系，并计算多项式系数，包括以下子步骤：
[0014]步骤5.1，以观测天顶角作为自变量，夜光辐亮度作为因变量，建立两者之间的天顶
‑
辐亮度二次(ZRQ)模型。
[0015]步骤5.2，计算ZRQ模型的二次函数系数。
[0016]步骤6，基于步骤5建立的天顶
‑
辐亮度二次(ZRQ)统计模型，消除夜光辐亮度角度效应，包括以下子步骤：
[0017]步骤6.1，基于步骤5建立ZRQ统计模型，将模型公式分解为天顶角为0
°
时的夜光辐亮度(常数项系数)与二次函数值的乘积形式。
[0018]步骤6.2，将步骤4提取的夜光辐亮度时间序列{R1,R2,...,R
n
}除以步骤5得到的二次模型的常数项系数c，获得新的夜光辐亮度时间序列并用新的夜光辐亮度时间序列与对应时刻的卫星传感器观测天顶角时间序列{z1,z2,...,z
n
}进行二次建模，使用如步骤5.2中的最小二乘法多项式曲线拟合求解二次函数系数a
′
、b
′
。
[0本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种卫星遥感夜光辐亮度的观测角度归一化方法，其特征在于，包括如下步骤：步骤1，获取研究区域Black Marble夜光影像产品中VNP46A1与VNP46A2数据；步骤2，对步骤1获取到的所有VNP46A1与VNP46A2数据进行数据格式转换，并得到需要的质量文件与影像文件数据；步骤3，对步骤2所得到的质量标志数据，设置一定的筛选条件，得到符合需求的高质量夜间灯光数据；步骤4，以像素为研究单位，从步骤3得到的高质量夜间灯光数据中提取像素的夜光辐亮度时间序列和传感器观测天顶角时间序列；步骤5，对步骤4提取的夜光辐亮度与卫星传感器观测天顶角进行统计建模，确立两者之间的量化关系，并计算多项式系数；步骤6，基于步骤5建立的天顶
‑
辐亮度二次模型，消除夜光辐亮度角度效应；步骤7，计算夜光辐亮度时间序列角度归一化后的变异系数与归一化冷热点夜光辐亮度，量化角度效应在夜间灯光时间序列中的影响程度。2.如权利要求1所述的一种卫星遥感夜光辐亮度的观测角度归一化方法，其特征在于：所述步骤2中将下载获取的HDF5格式科学数据集转换成GeoTiff格式，选取VNP46A1产品中传感器观测天顶角、太阳天顶角、月光照度分数文件，以及VNP46A2产品中经过BRDF校正的夜间灯光影像、云层掩码、强制性质量指标文件，其中太阳天顶角、月光照度分数、云层掩码以及强制性质量指标文件作为后续像素筛选的质量标志数据，经过BRDF校正的夜间灯光影像与传感器观测天顶角将用于天顶角与夜光辐亮度的建模。3.如权利要求1所述的一种卫星遥感夜光辐亮度的观测角度归一化方法，其特征在于：所述步骤3中设置云层掩码指标、强制质量指标、太阳天顶角和月光照度分数作为逐像素夜间灯光质量的筛选条件，但是在云层筛选的基础上，仍然存在部分像素值异常的情况，因此，将强制性质量指标设置为0，以保证所选像素均为高质量，从而提高数据质量，剔除太阳天顶角大于108
°
或月光照度分数小于1％的像素，以保证获取像素的有效性和降低云层掩码标志数据质量错误的可能性。4.如权利要求1所述的一种卫星遥感夜光辐亮度的观测角度归一化方法，其特征在于：所述步骤5包括以下子步骤：步骤5.1，以观测天顶角作为自变量，夜光辐亮度作为因变量，建立两者之间的天顶
‑
辐亮度二次模型；步骤5.2，计算天顶
‑
辐亮度二次模型的二次函数系数。5.如权利要求4所述的一种卫星遥感夜光辐亮度的观测角度归一化方法，其特征在于：所述步骤5.1是以观测天顶角作为自变量，夜光辐亮度作为因变量，建立两者之间的天顶
‑
辐亮度二次模型拟合关系如下：R＝az2+bz+c
ꢀꢀꢀꢀꢀꢀꢀꢀ
(1)式中，z为卫星观测天顶角，R为夜光辐亮度，a、b、c分别为二次模型系数。6.如权利要求4所述的一种卫星遥感夜光辐亮度的观测角度归一化方法，其特征在于：所述步骤5.2中中假设夜光辐亮度时间序列为{R1,R2,...,R
n
}，传感器观测天顶角时间序列为{z1,z2,...,z
n
}，在夜光辐亮度与观测天顶角拟合二次曲线过程中定义拟合的目标函数使得拟合偏差最小，目标函数公式如下：
式中，z
i
与R
i
分别为时间序列中第i时刻的观测天顶角与夜光辐亮度，为拟合的夜光辐亮度，n为时间序列长度；建立误差方程，并求解二次曲线...

【专利技术属性】
技术研发人员：李熙，贾明晖，
申请(专利权)人：武汉大学，
类型：发明
国别省市：