CN115049834B - 基于夜间灯光数据与高分辨率影像的城市建成区提取方法 - Google Patents

基于夜间灯光数据与高分辨率影像的城市建成区提取方法 Download PDF

Info

Publication number
CN115049834B
CN115049834B CN202210971201.0A CN202210971201A CN115049834B CN 115049834 B CN115049834 B CN 115049834B CN 202210971201 A CN202210971201 A CN 202210971201A CN 115049834 B CN115049834 B CN 115049834B
Authority
CN
China
Prior art keywords
area
remote sensing
built
sensing image
night light
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210971201.0A
Other languages
English (en)
Other versions
CN115049834A (zh
Inventor
杨柳林
柳敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangsu Dianboshi Energy Equipment Co ltd
Original Assignee
Nantong Electric Doctor Automation Equipment Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nantong Electric Doctor Automation Equipment Co ltd filed Critical Nantong Electric Doctor Automation Equipment Co ltd
Priority to CN202210971201.0A priority Critical patent/CN115049834B/zh
Publication of CN115049834A publication Critical patent/CN115049834A/zh
Application granted granted Critical
Publication of CN115049834B publication Critical patent/CN115049834B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/20—Image preprocessing
    • G06V10/26—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
    • G06V10/267—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion by performing operations on regions, e.g. growing, shrinking or watersheds
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/20—Image preprocessing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/20—Image preprocessing
    • G06V10/30—Noise filtering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
    • G06V10/774—Generating sets of training patterns; Bootstrap methods, e.g. bagging or boosting
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00—Scenes; Scene-specific elements
    • G06V20/70—Labelling scene content, e.g. deriving syntactic or semantic representations

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Multimedia (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computational Linguistics (AREA)
  • Computing Systems (AREA)
  • Databases & Information Systems (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Image Analysis (AREA)

Abstract

本发明涉及城市遥感信息技术应用领域，具体涉及一种基于夜间灯光数据与高分辨率影像的城市建成区提取方法，包括：采集同一区域的夜间灯光遥感图像和高分辨率全色遥感图像，并分别对这两种图像进行预处理；对高分辨率全色遥感图像样本内的建筑区域与其它区域进行分割标注；将同一区域的高分辨率全色遥感图像样本和夜间灯光遥感图像的分辨率进行调整，使两种图像的分辨率相同；利用同一区域的两种图像中每行每列的像素及其权重构建交叉熵损失函数；利用交叉熵损失函数得到训练好的语义分割网络；将预处理后的两种图像输入到训练好的语义分割网络中，完成城市建成区的提取。上述方法用于城市建成区的提取，可提高建成区的提取准确性。

Description

基于夜间灯光数据与高分辨率影像的城市建成区提取方法
技术领域
本发明涉及城市遥感信息技术应用领域，具体涉及一种基于夜间灯光数据与高分辨率影像的城市建成区提取方法。
背景技术
NPP/VIIRS夜间灯光数据由于其较高的分辨率以及覆盖面广、成本低廉、效率高等特点，被广泛用来衡量人类活动的强度与广度，主要应用于城市建成区提取、区域经济、地缘政治等研究领域。在目前利用夜间灯光数据提取城市建成区的研究中，常用的最佳阈值分割法主要有四种：经验阈值法、突变检测法、统计数据验证法和较高分辨率影像比较法。关于高分辨率影像提取城市建成区的研究大致可以分为两类：一类是基于区域增长方法的半自动提取方法，另一类方法则主要是基于分类的思想，根据光谱和纹理特征进行分类并作进一步处理以达到提取城市建成区的目的。
然而，利用NPP/VIIRS夜间灯光数据提取城市建成区仍然存在挑战：主要是NPP/VIIRS夜间灯光数据受到溢出效应的影响。由于非相干光从光源向各个方向辐射，比如，灯光亮度在周围区域的弥散，造成灯光亮度的溢出效应，导致对城市土地用地面积的高估，限制了NPP/VIIRS夜间灯光数据在城市土地用地精确提取的广泛应用。
发明内容
为了克服上述现有技术的缺点，本发明提供了一种基于夜间灯光数据与高分辨率影像的城市建成区提取方法，以实现对城市建成区的精确提取。
为达到上述目的，本发明采用以下技术方案，一种基于夜间灯光数据与高分辨率影像的城市建成区提取方法，包括：
S1：采集同一区域的NPP/VIIRS夜间灯光遥感图像和高分辨率全色遥感图像，并分别对这两种图像进行预处理；
S2：构建语义分割网络模型：
S201：对预处理后的高分辨率全色遥感图像样本内的建筑区域与其它区域进行分割标注；
S202：将同一区域的高分辨率全色遥感图像样本和夜间灯光遥感图像的分辨率进行调整，使同一区域的高分辨率全色遥感图像样本和夜间灯光遥感图像的分辨率相同；
S203：利用同一区域的高分辨率全色遥感图像样本和夜间灯光遥感图像中每行每列的像素及其权重构建交叉熵损失函数；
S204：利用构建的交叉熵损失函数对语义分割网络模型进行监督训练，直至语义分割网络收敛为止，得到训练好的语义分割网络；
S3：将预处理后的同一区域的NPP/VIIRS夜间灯光遥感图像和高分辨率全色遥感图像输入到训练好的语义分割网络模型中，完成城市建成区的提取。
所述一种基于夜间灯光数据与高分辨率影像的城市建成区提取方法，所述交叉熵损失函数的表达式为：
Figure DEST_PATH_IMAGE001
式中，
Figure 324475DEST_PATH_IMAGE002
表示交叉熵损失函数，M、N分别为图像的宽、高，
Figure 321250DEST_PATH_IMAGE003
表示图像中第i行第j列 像素的权重，
Figure 111351DEST_PATH_IMAGE004
表示语义分割网络预测的第i行第j列像素的概率，
Figure 927997DEST_PATH_IMAGE005
表示语义分割网络的 第i行第j列像素的标签值。
所述一种基于夜间灯光数据与高分辨率影像的城市建成区提取方法，所述图像中第i行第j列像素的权重的表达式为：
Figure 228791DEST_PATH_IMAGE006
=
Figure 814493DEST_PATH_IMAGE007
式中，
Figure 572234DEST_PATH_IMAGE003
表示图像中第i行第j列像素的权重，
Figure 876176DEST_PATH_IMAGE008
分别表示三种不同情 况得到的第i行第j列像素的初始权重。
所述一种基于夜间灯光数据与高分辨率影像的城市建成区提取方法，所述
Figure 712152DEST_PATH_IMAGE009
是 按照以下方式获取：
将夜间灯光遥感图像分为三个区域，该三个区域包括：建成区、过渡区和非建成区；
对建成区和非建成区进行中值滤波去噪，然后在两个区域中各选取一个中心邻域，并对中心邻域内各个像素的灰度值按大小进行排序，选取排序序列中的中值作为该中心邻域中心像素的灰度值；
对建成区和非建成区进行形态重构，获取形态重构后的夜间灯光遥感图像；
获取形态重构后的夜间灯光遥感图像的梯度图像；
利用分水岭分割算法提取梯度图像中的建成区；
将梯度图像中的建成区和非建成区进行标记，得到城市建成区的二值图，将该城 市建成区的二值图的像素值作为
Figure 417939DEST_PATH_IMAGE009
。
所述一种基于夜间灯光数据与高分辨率影像的城市建成区提取方法，所述
Figure 815423DEST_PATH_IMAGE010
的 计算公式为：
Figure 606661DEST_PATH_IMAGE011
式中，
Figure 249257DEST_PATH_IMAGE012
表示取最优阈值时所有城市斑块的周长之和，
Figure 809552DEST_PATH_IMAGE013
表示每个城市斑块的平 均灯光强度，
Figure 643515DEST_PATH_IMAGE014
表示高斯核函数，
Figure 656471DEST_PATH_IMAGE015
表示每个城市斑块中像素的坐标；
所述
Figure 670187DEST_PATH_IMAGE012
的表达式为：
Figure 84987DEST_PATH_IMAGE016
式中，
Figure 355432DEST_PATH_IMAGE017
表示取最优阈值时第n个城市斑块的周长，Q表示取最优阈值时城市斑块 的数量；
所述
Figure 855683DEST_PATH_IMAGE013
的表达式为：
Figure 840082DEST_PATH_IMAGE018
=
Figure 374968DEST_PATH_IMAGE019
式中：
Figure 550735DEST_PATH_IMAGE020
表示取最优阈值时第
Figure 803862DEST_PATH_IMAGE021
个城市斑块的年平均夜间灯光总量；
所述
Figure 589022DEST_PATH_IMAGE014
的表达式为：
Figure 447256DEST_PATH_IMAGE022
式中，
Figure 59503DEST_PATH_IMAGE023
是带宽。
所述一种基于夜间灯光数据与高分辨率影像的城市建成区提取方法，所述最优阈值是按照以下方式获取：
设置初始阈值，并将该初始阈值为基础，以间隔权重为间隔进行迭代增加；
利用迭代后的结果获得城市斑块的周长随阈值迭代的分布情况；
找出城市斑块的周长突变时所对应的辐射值，该辐射值减掉间隔权重即为最优阈值。
所述一种基于夜间灯光数据与高分辨率影像的城市建成区提取方法，所述
Figure 534347DEST_PATH_IMAGE024
是 按照以下方式获取：
从高分辨率全色遥感图像样本中提取NDVI和NDBI，其中NDVI为归一化植被指数， NDBI为归一化建筑指数，
Figure DEST_PATH_IMAGE025
和
Figure 188444DEST_PATH_IMAGE026
的表达式分别为：
Figure 432344DEST_PATH_IMAGE027
Figure 684334DEST_PATH_IMAGE028
式中，
Figure 145009DEST_PATH_IMAGE029
为近红外波段，
Figure 39015DEST_PATH_IMAGE030
为红外波段，
Figure DEST_PATH_IMAGE031
为中红外波段；
采用改进的城市夜间灯光指数法VBANUI提取城市建成区，
Figure 465317DEST_PATH_IMAGE032
的计算公式如 下：
Figure DEST_PATH_IMAGE033
式中，
Figure 451990DEST_PATH_IMAGE034
为夜间灯光遥感图像；
在完成
Figure 901426DEST_PATH_IMAGE032
指标的计算后，令
Figure 474490DEST_PATH_IMAGE035
式中，
Figure 427402DEST_PATH_IMAGE036
为每个城市斑块的
Figure 785309DEST_PATH_IMAGE032
最大值。
本发明的有益效果是：通过结合NPP/VIIRS夜间灯光数据和高分辨率影像提供了一种城市建成区的提取方法，该提取方法可以有效减少NPP/VIIRS夜间灯光数据溢出效应的影响，尤其在水体、建筑区内的植被区域等的夜间灯光亮度溢出现象，弥补当前已有方法的不足。同时本发明构建一种语义分割网络模型：利用夜间灯光数据提取不同形态城市建成区的最佳阈值，获取加权交叉熵。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明的城市建成区提取方法的流程示意图；
图2是本发明的语义分割网络模型的构建过程示意图；
图3是本发明的城市建成区的最优阈值获取过程示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例
本实施例提供了一种基于夜间灯光数据与高分辨率影像的城市建成区提取方法，如图1所示，包括：
S1：采集同一区域的NPP/VIIRS夜间灯光遥感图像和高分辨率全色遥感图像，并分别对这两种图像进行预处理。
1.对NPP/VIIRS夜间灯光遥感图像进行预处理。
首先用ArcMap完成NPP/VIIRS夜间灯光遥感图像的裁剪、重采样和重投影，裁剪、重采样和重投影属于常见的遥感图像预处理，过程细节不再赘述。NPP/VIIRS夜间灯光遥感图像的原始数据采用WGS1984坐标系，将其转换为Albers等积投影坐标系。然后将其重采样到0.5km×0.5km的网格中，本步骤主要是为了避免由坐标系引起的图像失真并简化图像区域的计算。
对NPP/VIIRS夜间灯光遥感图像进行年度均值图像的合成。在下载时月合成数据出现5、6月份缺失的情况，这可能是因为极光等污染光源的影响，且这种现象在每年的月合成数据中基本都会出现，因此采取了剔除这两个月的数据然后把剩余10个月数据求平均的方法来合成年度均值图像。合成公式如下所示：
Figure 456462DEST_PATH_IMAGE037
式中：
Figure 957850DEST_PATH_IMAGE038
表示第
Figure DEST_PATH_IMAGE039
个月的夜间灯光总量，
Figure 594630DEST_PATH_IMAGE040
表示年平均夜间灯光总量。
对NPP/VIIRS夜间灯光遥感图像进行不稳定光源和背景噪声消除、极值消除，把小于零的值设为0.001（近似为0，对统计分析没影响），把背景值设为0；把大于235的数值设为235（根据推理经验值），以删除一些异常值。
2.对高分辨率全色遥感图像进行预处理。
首先对高分辨率全色遥感图像进行辐射校正、大气校正，然后利用小波变换算法对校正后的图像进行去云处理，对图像进行直方图均衡化拉伸，其作用是通过拉伸像素强度分布范围来增强图像对比度。辐射校正、大气校正和直方图均衡化拉伸属于常见的遥感图像预处理，过程细节不再赘述。
S2：构建语义分割网络模型。
构建的网络模型采用基于深度学习的语义分割方法，具体有Unet、SegNet和FCN等模型，采用编码器-解码器结构。
如图2所示，本实施例所述语义分割网络模型的构建过程具体为：
S201：对预处理后的高分辨率全色遥感图像样本内的建筑区域与其它区域进行分割标注。
对高分辨率全色遥感图像样本进行标注，像素位置为建筑区域的标注为0，其它标注为1，最终得到一张像素值包含0,1的标注图像。
S202：将同一区域的高分辨率全色遥感图像样本和夜间灯光遥感图像的分辨率进行调整，使同一区域的高分辨率全色遥感图像样本和夜间灯光遥感图像的分辨率相同。
对同一区域的夜间灯光遥感图像和高分辨率全色遥感图像样本进行分块处理并调整分辨率。由于遥感图像数据分辨率较大，输入网络时需要进行分块处理，例如一个10000*10000的图像，可以分块成400块500*500大小的。同时由于两种图像分辨率不同，需要采用图像插值的方法把两种图像调到相同分辨率，常用的图像插值方法有：最邻近元法、双线性内插法和三次内插法。本实施例图像插值方法采用三次内插法，该方法利用三次多项式求逼近理论上最佳插值函数。
S203：利用同一区域的高分辨率全色遥感图像样本和夜间灯光遥感图像中每行每列的像素及其权重构建交叉熵损失函数。
语义分割网络模型采用交叉熵，它是对给定随机变量或事件集的两个概率分布之间的差异的度量，分类任务采用softmax激活函数+交叉熵损失函数，softmax激活函数将一个向量进行“归一化”成概率分布的形式，再采用交叉熵损失函数计算Loss。
交叉熵损失函数
Figure 624903DEST_PATH_IMAGE041
的具体表达为：
Figure 783352DEST_PATH_IMAGE001
其中M、N分别为图像的宽、高，
Figure 822852DEST_PATH_IMAGE003
表示第i行第j列像素的权重，加权交叉熵基于图 像位置进行加权，可以使得网络更注意VBANUI大的区域，且可以学习到DN与建成区的关系， 提高提取城市建成区的准确性。
Figure 983313DEST_PATH_IMAGE004
表示语义分割网络预测的第i行第j列像素的概率。
Figure 450066DEST_PATH_IMAGE005
表 示语义分割网络的第i行第j列像素的标签值。
上述加权交叉熵是基于图像位置进行加权，可以使得网络更注意VBANUI大的区域，且可以学习到DN与建成区的关系，提高城市建成区的提取准确性，每个像素的权重的获取方法具体如下：
（1）将夜间灯光遥感图像分为三个区域：建成区、非建成区和过渡区。将亮度辐射值大于c的像元作为建成区，像元亮度辐射值介于a和c的作为过渡区，亮度辐射值低于a的像元作为非建成区。将整幅夜间灯光遥感图像输入做如下处理：
①对建成区和非建成区进行中值滤波去噪，各选取一个中心邻域，然后对中心邻域内各个像素的灰度值按大小进行排序，选取排序序列中的中值作为该中心邻域中心像素点的灰度值；
②对建成区和非建成区进行形态重构，选择3维矩阵
Figure 95811DEST_PATH_IMAGE042
作为结构元， 对得到的图像进行闭运算，此操作可以较好地提高图像的轮廓清晰度，进而得到三个分区 的轮廓信息；
③利用Sobel算子计算梯度数据，该算子是一种1阶微分算子，它利用像素邻近区域的梯度值来计算1个像素的梯度，然后根据一定的绝对值来取舍。
④利用分水岭分割算法提取建成区，通常把梯度图像作为输入图像。
通过形态重构和标记得到二值图，其像素值记为
Figure 407844DEST_PATH_IMAGE009
，有效地抑制了提取建成区的 过程中存在欠分割的问题。
（2）对预处理后的夜间遥感图像采用周长突变检测法，由于城市建成区是聚合人类活动的整体，且NPP/VIIRS夜间灯光遥感图像分辨率较低，所以城市建成区内部应该很少有破碎斑块，且其夜间灯光数值应该比非城市建成区和城郊都更高。因此，随着设置的阈值不断增大，划归城市建成区的像元不断减少，城市建成区的周长相应也变小。但是，当阈值增大到一定程度时，城市建成区内部开始破碎，进而导致城市建成区周长没有变小，反而骤然增大。正确的阈值应该能够保证提取出来的城市结构的整体性，使之鲜少存在破碎区，所以将突变点之前的灯光值识别为最优阈值。
如图3所示，获取最优阈值的具体实现步骤如下：
1.设置初始阈值为
Figure 956899DEST_PATH_IMAGE043
，将辐射值不小于
Figure 594554DEST_PATH_IMAGE043
的像元划分为城市建成区像元；
2.将连续的城市像元聚合成城市斑块，计算城市斑块的总数量，记为Q；
3.从1到Q标识城市斑块，计算每一个城市斑块的周长，并求和；
4.城市斑块的周长计算公式如下：
Figure 727595DEST_PATH_IMAGE044
式中，
Figure 843319DEST_PATH_IMAGE017
表示取最优阈值时第n个城市斑块的周长，Q表示取最优阈值时城市斑块 的数量。求周长是斑块连通域的一个属性，通过OpenCV直接调用，这里不再赘述。
5.回到第一步，将阈值
Figure 243951DEST_PATH_IMAGE043
以0.01（根据推理经验值）为间隔，迭代增加；观察并统计 随着阈值的增加，城市斑块的周长的分布情况。
6.找出最优阈值。找到城市斑块的周长突然剧烈增加时对应的辐射值，这个辐射值减掉0.01就是最优阈值；
7.基于城市斑块的高斯赋权，采用高斯核函数，得到
Figure 786928DEST_PATH_IMAGE010
，高斯核函数可以将有限 维数据映射到高维空间，高斯核函数圆心为每个斑块的质心，其定义为：
Figure DEST_PATH_IMAGE045
其中，
Figure 344948DEST_PATH_IMAGE046
是带宽（根据经验值，在二维情况下
Figure 264362DEST_PATH_IMAGE046
=1），控制径向作用范围，换句话说，
Figure 522431DEST_PATH_IMAGE046
控制高斯核函数的局部作用范围。在图形上，正态分布是一种钟形曲线，越接近中心，取 值越大，越远离中心，取值越小。
Figure 236309DEST_PATH_IMAGE015
表示每个城市斑块中像素的坐标，当
Figure 609521DEST_PATH_IMAGE047
和
Figure 67047DEST_PATH_IMAGE048
的欧式距 离处于某一个区间范围内的时候，假设固定
Figure 911113DEST_PATH_IMAGE048
，
Figure 327051DEST_PATH_IMAGE049
随
Figure 656401DEST_PATH_IMAGE047
的变化而变化，且变化程度相当 显著。
8.提取每一个城市斑块内的年平均夜间灯光总量，最终每个城市斑块都得到一个权重，则每个城市斑块的平均灯光强度计算公式如下：
Figure 917618DEST_PATH_IMAGE018
=
Figure 884700DEST_PATH_IMAGE019
Figure 471539DEST_PATH_IMAGE020
表示取最优阈值时第
Figure 288185DEST_PATH_IMAGE021
个城市斑块的年平均夜间灯光总量。
Figure 87514DEST_PATH_IMAGE011
Figure 171751DEST_PATH_IMAGE012
表示取最优阈值时所有城市斑块的周长之和，
Figure 663912DEST_PATH_IMAGE013
表示每个城市斑块的平均灯光 强度，
Figure 702275DEST_PATH_IMAGE014
表示高斯核函数。
（3）从高分辨率全色遥感图像样本中提取NDVI和NDBI，其中NDVI为归一化植被指数，NDBI为归一化建筑指数，NDVI和NDBI的表达式如下：
Figure 305295DEST_PATH_IMAGE050
Figure 512548DEST_PATH_IMAGE028
其中，
Figure 910031DEST_PATH_IMAGE029
为近红外波段，
Figure 966849DEST_PATH_IMAGE030
为红外波段，
Figure 107980DEST_PATH_IMAGE031
为中红外波段。
采用改进的城市夜间灯光指数法VBANUI提取城市建成区，
Figure 190247DEST_PATH_IMAGE032
的计算公式如 下：
Figure 758632DEST_PATH_IMAGE033
其中，
Figure 37166DEST_PATH_IMAGE034
为夜间灯光遥感图像，
Figure 981988DEST_PATH_IMAGE025
为归一化植被指数，
Figure 537735DEST_PATH_IMAGE026
为归一化建筑 指数。
NDVI取值为[-1,1]，NDVI>0为植被覆盖区，NDVI<0为非植被覆盖区。表达式(1-NDVI)则表征了较大的非植被覆盖区权重在城市核心地区，(1-NDVI)与NTL的结合可减少城市核心区夜间灯光的饱和现象，并增加城市内部核心区变化特征的快速识别。
其次，NDBI取值范围也在[-1，1]之间。研究表明NDBI正值为城市用地区域，NDBI负值为非城市用地区域。融入NDBI后，在水体对城市建成区混淆或者NDVI和NTL有溢出效应时，都可以有效控制并减弱其影响。
在完成
Figure 44065DEST_PATH_IMAGE032
指标的计算后，令
Figure 544316DEST_PATH_IMAGE051
其中，
Figure 27250DEST_PATH_IMAGE036
为每个城市斑块的
Figure 827716DEST_PATH_IMAGE032
最大值。
至此，对于图像中的每个像素都有三个初始权重，分别为
Figure 236438DEST_PATH_IMAGE052
。
（4）由上述公式计算可知：
Figure 223986DEST_PATH_IMAGE053
Figure 510611DEST_PATH_IMAGE054
Figure 900004DEST_PATH_IMAGE055
最终得到每个像素的权重的表达式为：
Figure 748136DEST_PATH_IMAGE006
=
Figure 222980DEST_PATH_IMAGE007
将
Figure 47716DEST_PATH_IMAGE003
代入上述交叉熵损失函数中进行计算。
S204：利用构建的交叉熵损失函数对语义分割网络模型进行监督训练，直至语义分割网络收敛为止，得到训练好的语义分割网络。
S3：将预处理后的同一区域的NPP/VIIRS夜间灯光遥感图像和高分辨率全色遥感图像输入到训练好的语义分割网络模型中，完成城市建成区的提取。
将上述两种图像输入编码器，经过编码器特征提取，输出为特征图，解码器输入为特征图，经过解码器的上采样和拟合，得到目标分割图。
对语义分割网络进行优化，常用优化器有：SGD、Adam和Lookahead优化器，本实施例优选使用Lookahead优化器。
综上所述，本实施例提出的方法能够结合夜间灯光数据和高分辨率数据，提出加权交叉熵的方法，再引入NDVI和NDBI指标，有效减少了灯光亮度与NDVI溢出效应的影响，达到提高城市建成区提取准确性的目的。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (2)

1.一种基于夜间灯光数据与高分辨率影像的城市建成区提取方法，其特征在于，包括：
S1：采集同一区域的NPP/VIIRS夜间灯光遥感图像和高分辨率全色遥感图像，并分别对这两种图像进行预处理；
S2：构建语义分割网络模型：
S201：对预处理后的高分辨率全色遥感图像样本内的建筑区域与其它区域进行分割标注；
S202：将同一区域的高分辨率全色遥感图像样本和夜间灯光遥感图像的分辨率进行调整，使同一区域的高分辨率全色遥感图像样本和夜间灯光遥感图像的分辨率相同；
S203：利用同一区域的高分辨率全色遥感图像样本和夜间灯光遥感图像中每行每列的像素及其权重构建交叉熵损失函数；
所述交叉熵损失函数的表达式为：
Figure DEST_PATH_IMAGE002
式中，
Figure DEST_PATH_IMAGE004
表示交叉熵损失函数，M、N分别为图像的宽、高，
Figure DEST_PATH_IMAGE006
表示图像中第i行第j列像素的权重，
Figure DEST_PATH_IMAGE008
表示语义分割网络预测的第i行第j列像素的概率，
Figure DEST_PATH_IMAGE010
表示语义分割网络的第i行第j列像素的标签值；
所述图像中第i行第j列像素的权重的表达式为：
Figure DEST_PATH_IMAGE012
=
Figure DEST_PATH_IMAGE014
式中，
Figure 313137DEST_PATH_IMAGE006
表示图像中第i行第j列像素的权重，
Figure DEST_PATH_IMAGE016
分别表示三种不同情况得到的第i行第j列像素的初始权重；
所述
Figure DEST_PATH_IMAGE018
是按照以下方式获取：
将夜间灯光遥感图像分为三个区域，该三个区域包括：建成区、过渡区和非建成区；
对建成区和非建成区进行中值滤波去噪，然后在两个区域中各选取一个中心邻域，并对中心邻域内各个像素的灰度值按大小进行排序，选取排序序列中的中值作为该中心邻域中心像素的灰度值；
对建成区和非建成区进行形态重构，获取形态重构后的夜间灯光遥感图像；
获取形态重构后的夜间灯光遥感图像的梯度图像；
利用分水岭分割算法提取梯度图像中的建成区；
将梯度图像中的建成区和非建成区进行标记，得到城市建成区的二值图，将该城市建成区的二值图的像素值作为；
所述
Figure DEST_PATH_IMAGE020
的计算公式为：
Figure DEST_PATH_IMAGE022
式中，
Figure DEST_PATH_IMAGE024
表示取最优阈值时所有城市斑块的周长之和，
Figure DEST_PATH_IMAGE026
表示每个城市斑块的平均灯光强度，
Figure DEST_PATH_IMAGE028
表示高斯核函数，
Figure DEST_PATH_IMAGE030
表示每个城市斑块中像素的坐标；
所述
Figure 715431DEST_PATH_IMAGE024
的表达式为：
Figure DEST_PATH_IMAGE032
式中，
Figure DEST_PATH_IMAGE034
表示取最优阈值时第n个城市斑块的周长，Q表示取最优阈值时城市斑块的数量；
所述
Figure 69402DEST_PATH_IMAGE026
的表达式为：
Figure DEST_PATH_IMAGE036
=
Figure DEST_PATH_IMAGE038
式中，
Figure DEST_PATH_IMAGE040
表示取最优阈值时第
Figure DEST_PATH_IMAGE042
个城市斑块的年平均夜间灯光总量；
所述
Figure 841442DEST_PATH_IMAGE028
的表达式为：
Figure DEST_PATH_IMAGE044
式中，
Figure DEST_PATH_IMAGE046
是带宽；
所述
Figure DEST_PATH_IMAGE048
是按照以下方式获取：
从高分辨率全色遥感图像样本中提取NDVI和NDBI，其中NDVI为归一化植被指数，NDBI为归一化建筑指数，
Figure DEST_PATH_IMAGE050
和
Figure DEST_PATH_IMAGE052
的表达式分别为：
Figure DEST_PATH_IMAGE054
Figure DEST_PATH_IMAGE056
式中，
Figure DEST_PATH_IMAGE058
为近红外波段，
Figure DEST_PATH_IMAGE060
为红外波段，
Figure DEST_PATH_IMAGE062
为中红外波段；
采用改进的城市夜间灯光指数法VBANUI提取城市建成区，
Figure DEST_PATH_IMAGE064
的计算公式如下：
Figure DEST_PATH_IMAGE066
式中，
Figure DEST_PATH_IMAGE068
为夜间灯光遥感图像；
在完成
Figure 433311DEST_PATH_IMAGE064
指标的计算后，令
Figure DEST_PATH_IMAGE070
式中，
Figure DEST_PATH_IMAGE072
为每个城市斑块的最大值；
S204：利用构建的交叉熵损失函数对语义分割网络模型进行监督训练，直至语义分割网络收敛为止，得到训练好的语义分割网络；
S3：将预处理后的同一区域的NPP/VIIRS夜间灯光遥感图像和高分辨率全色遥感图像输入到训练好的语义分割网络模型中，完成城市建成区的提取。
2.根据权利要求1所述的一种基于夜间灯光数据与高分辨率影像的城市建成区提取方法，其特征在于，所述最优阈值是按照以下方式获取：
设置初始阈值，并将该初始阈值为基础，以间隔权重为间隔进行迭代增加；
利用迭代后的结果获得城市斑块的周长随阈值迭代的分布情况；
找出城市斑块的周长突变时所对应的辐射值，该辐射值减掉间隔权重即为最优阈值。
CN202210971201.0A 2022-08-15 2022-08-15 基于夜间灯光数据与高分辨率影像的城市建成区提取方法 Active CN115049834B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210971201.0A CN115049834B (zh) 2022-08-15 2022-08-15 基于夜间灯光数据与高分辨率影像的城市建成区提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210971201.0A CN115049834B (zh) 2022-08-15 2022-08-15 基于夜间灯光数据与高分辨率影像的城市建成区提取方法

Publications (2)

Publication Number Publication Date
CN115049834A CN115049834A (zh) 2022-09-13
CN115049834B true CN115049834B (zh) 2022-11-11

Family

ID=83167533

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210971201.0A Active CN115049834B (zh) 2022-08-15 2022-08-15 基于夜间灯光数据与高分辨率影像的城市建成区提取方法

Country Status (1)

Country Link
CN (1) CN115049834B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117612003A (zh) * 2023-11-27 2024-02-27 通友微电(四川)有限公司 一种基于多源遥感影像的城市建成区绿地变化识别方法
CN117495425B (zh) * 2023-12-29 2024-04-12 武汉大学 一种基于多维夜光特征的资产财富估算方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108596139A (zh) * 2018-05-03 2018-09-28 武汉大学 一种基于Gabor特征显著性的遥感影像城市区域提取方法
CN111597949A (zh) * 2020-05-12 2020-08-28 中国科学院城市环境研究所 一种基于npp-viirs夜间灯光数据的城市建成区提取方法
CN112989985A (zh) * 2021-03-08 2021-06-18 武汉大学 一种融合夜间灯光数据与Landsat8 OLI影像的城市建成区提取方法
US11189034B1 (en) * 2020-07-22 2021-11-30 Zhejiang University Semantic segmentation method and system for high-resolution remote sensing image based on random blocks

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108596139A (zh) * 2018-05-03 2018-09-28 武汉大学 一种基于Gabor特征显著性的遥感影像城市区域提取方法
CN111597949A (zh) * 2020-05-12 2020-08-28 中国科学院城市环境研究所 一种基于npp-viirs夜间灯光数据的城市建成区提取方法
US11189034B1 (en) * 2020-07-22 2021-11-30 Zhejiang University Semantic segmentation method and system for high-resolution remote sensing image based on random blocks
CN112989985A (zh) * 2021-03-08 2021-06-18 武汉大学 一种融合夜间灯光数据与Landsat8 OLI影像的城市建成区提取方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于NPP-VIIRS夜间灯光数据和Landsat-8数据的城镇建筑用地提取方法改进――以广州市为例;唐梁博等;《测绘与空间地理信息》;20170925(第09期);第79-83页 *

Also Published As

Publication number Publication date
CN115049834A (zh) 2022-09-13

Similar Documents

Publication Publication Date Title
CN115049834B (zh) 基于夜间灯光数据与高分辨率影像的城市建成区提取方法
CN107045629B (zh) 一种多车道线检测方法
CN109977801B (zh) 一种联合光学和雷达的区域水体快速动态提取方法及系统
Zhang et al. A fusion algorithm for infrared and visible images based on saliency analysis and non-subsampled Shearlet transform
CN109472200B (zh) 一种智能的海面垃圾检测方法、系统和存储介质
CN112101271A (zh) 一种高光谱遥感影像分类方法及装置
CN111027446B (zh) 一种高分辨率影像的海岸线自动提取方法
CN111881816B (zh) 一种长时序的河湖围埂养殖区域监测方法
CN110070545B (zh) 一种城镇纹理特征密度自动提取城镇建成区的方法
CN107992856B (zh) 城市场景下的高分遥感建筑物阴影检测方法
CN110176005B (zh) 基于归一化指数和多尺度模型的遥感图像分割方法
CN115861359B (zh) 一种水面漂浮垃圾图像自适应分割提取方法
CN112200083B (zh) 一种基于多元高斯混合模型的机载多光谱LiDAR数据分割方法
CN111310771B (zh) 遥感影像的道路图像提取方法、装置、设备及存储介质
CN115631372B (zh) 基于土壤遥感数据的土地信息分类管理方法
CN114842262A (zh) 融合线路通道正射影像的激光点云地物自动识别方法
CN116091937A (zh) 基于深度学习的高分辨率遥感影像地物识别模型计算方法
CN117079117B (zh) 水下图像处理及目标识别的方法及装置、存储介质及电子设备
CN113705433A (zh) 一种基于可见光航拍图像的电力线检测方法
Ju et al. A novel fully convolutional network based on marker-controlled watershed segmentation algorithm for industrial soot robot target segmentation
CN109165590A (zh) 利用稀疏形态成分的高分辨率遥感图像道路提取方法
Subhashini et al. An innovative hybrid technique for road extraction from noisy satellite images
Wang et al. Simultaneous extracting area and quantity of agricultural greenhouses in large scale with deep learning method and high-resolution remote sensing images
CN115294322A (zh) 水下船底可疑目标检测方法、装置、电子设备及可读介质
Borra et al. Satellite image enhancement and analysis

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP03 Change of name, title or address

Address after: 226200 group 11, Chengbei Village, Huilong Town, Qidong City, Nantong City, Jiangsu Province

Patentee after: Jiangsu Dianboshi Energy Equipment Co.,Ltd.

Address before: 226000 group 11, Chengbei Village, Huilong Town, Qidong City, Nantong City, Jiangsu Province

Patentee before: Nantong Electric doctor automation equipment Co.,Ltd.

CP03 Change of name, title or address