CN105740824B - 一种具有抗特征冲突的人口密集区识别方法 - Google Patents
一种具有抗特征冲突的人口密集区识别方法 Download PDFInfo
- Publication number
- CN105740824B CN105740824B CN201610069935.4A CN201610069935A CN105740824B CN 105740824 B CN105740824 B CN 105740824B CN 201610069935 A CN201610069935 A CN 201610069935A CN 105740824 B CN105740824 B CN 105740824B
- Authority
- CN
- China
- Prior art keywords
- dense population
- population areas
- pixel
- feature
- recognition methods
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/176—Urban or other man-made structures
Abstract
本发明涉及一种具有抗特征冲突的人口密集区识别方法。包括如下步骤:步骤1,对遥感影像进行预处理;步骤2,提取局部梯度方向密度、梯度均值、Harris特征点数和光谱同质性共四种特征;步骤3,计算各像素的特征取值并归一化;步骤4,计算各像素的特征概率函数取值;步骤5,特征融合;步骤6,保留人口密集区像素及可视化。作为人口密集区中人工地物提取的预处理步骤,可以缩小后续处理的区域面积,提高人工地物提取的准确率和算法效率。
Description
技术领域
本发明涉及一种遥感图像处理领域,具体说是一种具有抗特征冲突的人口密集区识别方法。
背景技术
由于受拍摄条件的影响,不同时期、不同地区的高分辨率遥感影像中的人口密集区表现出不同的显著特征,如果能获取显著特征的先验知识,则可以大大提高建筑物密集区的识别准确率,但这一条件往往难以满足。目前常见的方法有:(1)单一特征识别法,缺点是只对特定影像类别中的建筑物区域识别有效,(2)综合多种特征识别法,缺点是需要确定各个特征的权重,不合理的权重可能导致建筑物密集区识别的失败,而权重的确定一般采用训练或学习的方式,需要大量的样本;其次,当各个特征之间具有冲突时,识别的正确率将会大大降低。
发明内容
本发明提供了一种具有抗特征冲突的人口密集区识别方法,以遥感影像作为数据源,最大程度地挖掘遥感影像中人口密集区的有效信息,可以有效防止特征冲突导致识别正确率下降,且不涉及任何参数,达到完全自动化。
为实现本发明的目标所采用的技术方案是:方法包括以下步骤:
步骤1:对遥感影像image进行预处理;
步骤2:选取遥感影像中与人口密集区相关的局部梯度方向密度、梯度均值、Harris特征点数和光谱同质性共四种特征;
步骤3:对步骤2中四种特征在遥感影像image中各像素的取值进行归一化处理,分别得到矩阵E1、E2、E3和E4;
步骤4:利用矩阵E1、E2、E3和E4分别计算每个像素属于人口密集区、非人口密集区和(人口密集区、非人口密集区)的概率函数取值mi(b)、mi(nb)、mi(b,nb),其中i=1,2,3,4;
步骤5:利用特征融合方法对步骤4中的mi(b)、mi(nb)、mi(b,nb)进行组合,分别计算每个像素属于人口密集区、非人口密集区和(人口密集区、非人口密集区)的概率m(b)、m(nb)、m(b,nb);
步骤6:利用步骤5中的概率m(b)、m(nb)、m(b,nb)结合符合人口密集区条件的计算公式S保留相应的像素,并进行人口密集区的边界表达。
所述的局部梯度方向密度由以下公式计算:
其中,Ns为窗口内各个像素点的梯度模值之和,h为局部梯度方向密度函数的带宽,为像素点r的梯度幅值,k为核密度函数。
所述的梯度均值为所取窗口中各点的梯度模值的累加和与窗口面积的比值。
所述的Harris特征点数为落在所取窗口中的Harris特征点数量。
所述的光谱同质性所取窗口中像元灰度值的标准差。
所述的概率函数取值mi(b)、mi(nb)、mi(b,nb)由以下公式计算:
所述的特征融合方法由以下公式计算:
m(A)=p(A)+k×ε×q(A),A={b,nb,(b,nb)}
其中,
所述的符合人口密集区条件的计算公式S为:
本发明的有益效果是:作为人口密集区中人工地物提取的预处理步骤,可以缩小后续处理的区域面积,提高人工地物提取的准确率和算法效率。
附图说明
图1是本发明的总体处理流程图。
具体实施方式
下面结合附图详细描述本发明的具体实施方式。
在步骤101,输入待处理的遥感影像image为高空间分辨率遥感影像,可以是Quickbird、World view等卫星影像或各类航拍影像,空间分辨率在1米以下。
在步骤102,对image进行辐射校正、几何校正和大气校正等预处理。
在步骤103,提取image中的局部梯度方向密度、梯度均值、Harris特征点数和光谱同质性共四种特征。
在步骤104,分别以image中的各个像素为中心点,选取尺寸为11×11的窗口,分别计算步骤103中的特征值并进行归一化,将其取值赋予中心点,分别得到矩阵E1、E2、E3和E4。
在步骤105,利用公式
计算各像素的特征概率函数取值mi(b)、mi(nb)和mi(b,nb)。
在步骤106,利用特征融合公式
m(A)=p(A)+k×ε×q(A),A={b,nb,(b,nb)}
其中,
对步骤103中的四种特征进行融合。
在步骤107,利用人口密集区条件的计算公式
保留人口密集区的像素,并进行人口密集区的边界表达。
Claims (7)
1.一种具有抗特征冲突的人口密集区识别方法,其特征在于包括以下步骤:
步骤1:对遥感影像image进行预处理;
步骤2:选取遥感影像中与人口密集区相关的局部梯度方向密度、梯度均值、Harris特征点数和光谱同质性共四种特征;
步骤3:对步骤2中四种特征在遥感影像image中各像素的取值进行归一化处理,分别得到矩阵E1、E2、E3和E4;
步骤4:利用矩阵E1、E2、E3和E4分别计算每个像素属于人口密集区、非人口密集区和(人口密集区、非人口密集区)的概率函数取值mi(b)、mi(nb)、mi(b,nb),其中i=1,2,3,4;
步骤5:利用特征融合方法对步骤4中的mi(b)、mi(nb)、mi(b,nb)进行组合,分别计算每个像素属于人口密集区、非人口密集区和(人口密集区、非人口密集区)的概率m(b)、m(nb)、m(b,nb);
步骤6:利用步骤5中的概率m(b)、m(nb)、m(b,nb)结合符合人口密集区条件的计算公式S保留相应的像素,并进行人口密集区的边界表达;
所述的符合人口密集区条件,其计算公式S为:
2.根据权利要求1所述的一种具有抗特征冲突的人口密集区识别方法,其特征在于所述的局部梯度方向密度由以下公式计算:
其中,Ns为窗口内各个像素点的梯度模值之和,h为局部梯度方向密度函数的带宽,为像素点r的梯度幅值,k为核密度函数。
3.根据权利要求1所述的一种具有抗特征冲突的人口密集区识别方法,其特征在于所述的梯度均值为所取窗口中各点的梯度模值的累加和与窗口面积的比值。
4.根据权利要求1所述的一种具有抗特征冲突的人口密集区识别方法,其特征在于所述的Harris特征点数为落在所取窗口中的Harris特征点数量。
5.根据权利要求1所述的一种具有抗特征冲突的人口密集区识别方法,其特征在于所述的光谱同质性所取窗口中像元灰度值的标准差。
6.根据权利要求1所述的一种具有抗特征冲突的人口密集区识别方法,其特征在于所述的概率函数取值mi(b)、mi(nb)、mi(b,nb)由以下公式计算:
7.根据权利要求1所述的一种具有抗特征冲突的人口密集区识别方法,其特征在于所述的特征融合方法,由以下公式计算:
m(A)=p(A)+k×ε×q(A),A={b,nb,(b,nb)}
其中,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610069935.4A CN105740824B (zh) | 2016-02-01 | 2016-02-01 | 一种具有抗特征冲突的人口密集区识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610069935.4A CN105740824B (zh) | 2016-02-01 | 2016-02-01 | 一种具有抗特征冲突的人口密集区识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105740824A CN105740824A (zh) | 2016-07-06 |
CN105740824B true CN105740824B (zh) | 2019-10-22 |
Family
ID=56245680
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610069935.4A Expired - Fee Related CN105740824B (zh) | 2016-02-01 | 2016-02-01 | 一种具有抗特征冲突的人口密集区识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105740824B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112818747A (zh) * | 2020-12-31 | 2021-05-18 | 上海应用技术大学 | 一种基于空间大数据的城市特色街区人口密度估算方法和系统方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102622607B (zh) * | 2012-02-24 | 2013-09-25 | 河海大学 | 一种基于多特征融合的遥感图像分类方法 |
CN103324945B (zh) * | 2013-07-08 | 2016-12-28 | 南京大学 | 一种基于模式识别的森林点云分类方法 |
CN104463442B (zh) * | 2014-11-25 | 2018-01-26 | 中国科学院地理科学与资源研究所 | 一种城乡建设集聚性的探测方法 |
CN104680173B (zh) * | 2015-01-26 | 2018-04-27 | 河海大学 | 一种遥感图像场景分类方法 |
-
2016
- 2016-02-01 CN CN201610069935.4A patent/CN105740824B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN105740824A (zh) | 2016-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107067415B (zh) | 一种基于图像匹配的目标定位方法 | |
Zhong et al. | An adaptive subpixel mapping method based on MAP model and class determination strategy for hyperspectral remote sensing imagery | |
CN110929607B (zh) | 一种城市建筑物施工进度的遥感识别方法和系统 | |
Rao et al. | Textural analysis of IRS-1D panchromatic data for land cover classification | |
CN107358260B (zh) | 一种基于表面波cnn的多光谱图像分类方法 | |
CN108428220B (zh) | 静止轨道卫星序列遥感影像海岛礁区域自动几何校正方法 | |
CN112149594B (zh) | 基于深度学习和高分辨率卫星影像的城市建设评估方法 | |
CN106910202B (zh) | 一种遥感图像地物的图像分割方法及系统 | |
CN108960404B (zh) | 一种基于图像的人群计数方法及设备 | |
CN105354841B (zh) | 一种快速遥感影像匹配方法及系统 | |
Yuan et al. | Learning to count buildings in diverse aerial scenes | |
CN115713694B (zh) | 一种土地测绘信息管理方法 | |
CN109377511B (zh) | 基于样本组合和深度检测网络的运动目标跟踪方法 | |
CN108537169B (zh) | 一种基于中心线和路宽检测算法的高分辨率遥感影像道路提取方法 | |
CN112766184B (zh) | 基于多层级特征选择卷积神经网络的遥感目标检测方法 | |
CN107909053B (zh) | 一种基于等级学习级联卷积神经网络的人脸检测方法 | |
CN102446356A (zh) | 一种获取均匀分布匹配点的遥感影像并行自适应匹配方法 | |
CN108710909A (zh) | 一种可变形旋转不变装箱物体清点方法 | |
CN104537381B (zh) | 一种基于模糊不变特征的模糊图像识别方法 | |
CN110222609A (zh) | 一种基于图像处理的墙体裂缝智能识别方法 | |
Su et al. | Hyperspectral image classification based on volumetric texture and dimensionality reduction | |
CN106709515A (zh) | 一种适用于下视景象匹配选区准则干预方法 | |
Du et al. | Open-pit mine change detection from high resolution remote sensing images using DA-UNet++ and object-based approach | |
Norzaki et al. | A comparative study of template matching, ISO cluster segmentation, and tree canopy segmentation for homogeneous tree counting | |
CN109558803B (zh) | 基于卷积神经网络与np准则的sar目标鉴别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191022 Termination date: 20220201 |
|
CF01 | Termination of patent right due to non-payment of annual fee |