连续夜光遥感数据异常判别及插补方法技术

本发明专利技术提供连续夜光遥感数据异常判别及插补方法，包括获取连续的夜光遥感数据，还包括以下步骤：将日期连续的每幅夜光遥感数据转换为表格文件；计算每个表格文件与相邻日期的2N个表格文件平均值的相关系数；利用多步循环，每次都删除相关系数最小的表格文件，然后重新计算剩余表格文件的相关系数，直到所有相关系数均大于设定阈值；使用筛选出的合格的文件，插值得到删除日期对应的表格文件；将插值得到的表格文件按照行列对应关系转换成夜光遥感影像。本发明专利技术利用每日夜光遥感数据的全部像元值与相邻日期夜光遥感数据的亮度值的相关性剔除异常数据；利用最近日期的正常数据，逐像元进行插值，从而获取逐日连续的夜光遥感数据。数据。数据。

【技术实现步骤摘要】
连续夜光遥感数据异常判别及插补方法


[0001]本专利技术涉及降水监测的
，特别是连续夜光遥感数据异常判别及插补方法。

技术介绍

[0002]随着社会经济的发展，夜间照明设施逐渐普及。夜晚灯光的强弱，直接彰显了人类活动区域的发达程度，夜光卫星就是利用遥感技术从太空观测夜间地球的灯光，通过传感器获取夜间无云条件下地表发射的可见光
‑
近红外电磁波信息而形成的灰度影像。夜光卫星最早起源于上世纪70年代美国军事气象卫星计划，近几年，夜间灯光卫星遥感受到越来越多的关注，广泛应用于城市扩张、人口分布、社会经济指标提取、生态环境评估、光污染监测、突发事件监测等。目前在全球应用最为广泛、序列最长的夜间卫星数据为VIIRS/DNB（Visible Infrared Imaging Radiometer SuiteDay/Night Band），是NASA和NOAA根据SNPP(Suomi National PolarPartnership)卫星上的可见红外辐射套件 (VisibleInfrared Imaging Radiometer Suite (VIIRS)) 的昼夜波段 (Day/Night Band (DNB)) 收集的全球夜间灯光数据，该卫星于2011年发射成功，VIIRS/DNB的空间分辨率在500米左右。
[0003]VIIRS/DNB数据集，免费提供了2012年4月1日以来的逐日数据。逐日夜光遥感数据由于受云层、月光、雾霾、不稳定光源等因素的影响，会导致当日数据缺失或出现异常。在实际应用中进行夜光变化分析时，需要连续的逐日夜光遥感数据，因此需要从大量夜光遥感数据中自动删除异常数据，并对删除日期的夜光遥感数据进行插值。
[0004]2011年第5期《测绘通报》上公开了江振治、胡海彦、王刃和陈虹的题目为《各类卫星遥感影像间插补漏洞方法研究》的文章，该文章提出为填补某类卫星摄影覆盖区内建立控制网的局部漏洞,以该类卫星影像上已知的像点坐标为控制,先利用另一类卫星影像插补待求点在该类卫星影像中的像点坐标,再利用插补后影像数据进行构网平差。在平原地区计算过程使用9种变换公式,山区使用4种效果较好的变换公式,综合评价平面和高程精度,以选用二三次混合多项式的效果最佳。该方法的缺点是仅对影像漏洞进行插值，未考虑影像未形成漏洞但数值异常或错误的情况，二是未涉及连续多期影像的插补方法。
[0005]申请号为CN111292254A的专利技术专利申请公开了一种遥感影像中数据异常的确定方法、装置、设备及介质，所述方法包括：将待检测的遥感图像进行矩阵格式处理，得到图像信息；根据所述图像信息中的波段信息，绘制散点图，并计算所述散点图的平均值和标准差；采用预先设置的N个特定倍数对所述标准差进行计算，得到N个比较阈值；根据所述平均值和所述N个比较阈值，确定所述遥感图像是否存在数据异常。该方法的缺点是仅利用单幅影像的像元特征对单幅影像的噪声进行判别和去除，无法对连续的多期影像进行异常判别，而且采用确定阈值的判别方法可能会导致噪声的误判。

技术实现思路

[0006]为了解决上述的技术问题，本专利技术提出的连续夜光遥感数据异常判别及插补方法，将逐日夜光遥感数据转换为表格文件，利用每日夜光遥感数据的全部像元值与相邻日期夜光遥感数据的亮度值的相关性剔除异常数据；利用最近日期的正常数据，逐像元进行插值，获取每个像元的正常值，从而获取逐日连续的夜光遥感数据。
[0007]本专利技术的目的是提供连续夜光遥感数据异常判别及插补方法，包括获取连续的夜光遥感数据，还包括以下步骤：步骤1：将日期连续的每幅夜光遥感数据转换为表格文件；步骤2：计算每个表格文件与相邻日期的2N个表格文件平均值的相关系数；步骤3：利用多步循环，每次都删除相关系数最小的表格文件，然后重新计算剩余表格文件的相关系数，直到所有相关系数均大于设定阈值；步骤4：使用筛选出的合格的文件，插值得到删除日期对应的表格文件；步骤5：将步骤4得到的表格文件按照行列对应关系转换成夜光遥感影像文件。
[0008]优选的是，第p行第q列的影像像元值来自表格中第p*q行的数据。
[0009]在上述任一方案中优选的是，所述步骤2包括以下子步骤：步骤21：选择每个表格文件对应日期的前后各N日，共2N日的表格文件；步骤22：计算日期为d的表格文件的相邻2N日的像元平均值；步骤23：计算每个表格文件与相邻2N日均值的相关系数。
[0010]在上述任一方案中优选的是，日期为d的表格文件的相邻2N日的像元平均值计算结果中第i行的值的计算公式为：
[0011]其中，A
d,i
表示日期为d的表格文件的相邻2N日的像元平均值计算结果中第i行的值，v
d
‑
n,i
表示第d
‑
n日的表格文件中第i行的值，v
d+n,i
表示第d+n日的表格文件中第i行的值，1≤n≤N。
[0012]在上述任一方案中优选的是，日期为d的相关系数r
d
的计算公式为，，，其中，m表示表格文件的行数，v
d,i
表示日期为d的表格文件第i行的值，表示日期为d的所有行的平均值，表示相邻2N日表格文件平均值中所有行的平均值。
[0013]在上述任一方案中优选的是，所述步骤3还包括每次删除相关系数最小的表格文件之后，选择相邻表格时，跳过删除文件对应的日期，直到找到日期最近的2N个表格文件。
[0014]在上述任一方案中优选的是，所述插值公式为，其中，v
D,i
表示日期为D的待插值表格第i行的值，D1表示日期D之前离D最近的合格文件的日期，D2表示日期D之后离D最近的合格文件的日期，v
D1,i
表示日期为D1的表格文件中第i行的值；v
D2,i
表示日期为D2的表格文件中第i行的值，D－D1和D2－D1表示的是相差的天数。
[0015]本专利技术提出了连续夜光遥感数据异常判别及插补方法，可以从多期连续的逐日夜光遥感数据中自动剔除掉异常数据，并对剔除掉日期的夜光遥感数据进行插值，从而形成日期连续的逐日夜光遥感数据。
附图说明
[0016]图1为按照本专利技术的连续夜光遥感数据异常判别及插补方法的一优选实施例的流程图。
[0017]图2为按照本专利技术的连续夜光遥感数据异常判别及插补方法的另一优选实施例的流程图。
[0018]图3为按照本专利技术的连续夜光遥感数据异常判别及插补方法的部分影像原始效果的一实施例的示意图。
[0019]图4为按照本专利技术的连续夜光遥感数据异常判别及插补方法的夜光遥感数据插值前后对比的一实施例的示意图。
具体实施方式
[0020]下面结合附图和具体的实施例对本专利技术做进一步的阐述。
实施例
[0021]如图1所示，连续夜光遥感数据异常判别及插补方法，执行步骤100，获取连续的夜光遥感数据。
[0022]执行步骤110，将日期连续的每幅夜光遥本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.连续夜光遥感数据异常判别及插补方法，包括获取连续的夜光遥感数据，还包括以下步骤：步骤1：将日期连续的每幅夜光遥感数据转换为表格文件；步骤2：计算每个表格文件与相邻日期的2N个表格文件平均值的相关系数；步骤3：利用多步循环，每次都删除相关系数最小的表格文件，然后重新计算剩余表格文件的相关系数，直到所有相关系数均大于设定阈值；步骤4：使用筛选出的合格的文件，插值得到删除日期对应的表格文件；步骤5：将步骤4得到的表格文件按照行列对应关系转换成夜光遥感影像文件。2.如权利要求1所述的连续夜光遥感数据异常判别及插补方法，其特征在于，第p行第q列的影像像元值来自表格中第p*q行的数据。3.如权利要求2所述的连续夜光遥感数据异常判别及插补方法，其特征在于，所述步骤2包括以下子步骤：步骤21：选择每个表格文件对应日期的前后各N日，共2N日的表格文件；步骤22：计算日期为d的表格文件的相邻2N日的像元平均值；步骤23：计算每个表格文件与相邻2N日均值的相关系数。4.如权利要求3所述的连续夜光遥感数据异常判别及插补方法，其特征在于，日期为d的表格文件的相邻2N日的像元平均值计算结果中第i行的值的计算公式为：，其中，A
d,i
表示日期为d的表格文件的相邻2N日的像元平均值计算结果中第i行的值，v
d
‑

【专利技术属性】
技术研发人员：刘业森，刘媛媛，刘舒，李敏，臧文斌，郝晓丽，黄耀欢，
申请(专利权)人：中国水利水电科学研究院，
类型：发明
国别省市：