CN110472559A - 一种遥感影像建筑区域土地利用属性空间迁移方法 - Google Patents
一种遥感影像建筑区域土地利用属性空间迁移方法 Download PDFInfo
- Publication number
- CN110472559A CN110472559A CN201910743346.3A CN201910743346A CN110472559A CN 110472559 A CN110472559 A CN 110472559A CN 201910743346 A CN201910743346 A CN 201910743346A CN 110472559 A CN110472559 A CN 110472559A
- Authority
- CN
- China
- Prior art keywords
- building
- attribute
- spot
- poi
- land use
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/24323—Tree-organised classifiers
-
- 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/25—Determination of region of interest [ROI] or a volume of interest [VOI]
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Multimedia (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明揭示了一种遥感影像土地利用属性空间迁移方法,包括获取遥感影像建筑区域的POI数据,并对POI数据进行建筑属性清洗,进一步将获得的建筑属性与土地利用属性建立映射关系;将遥感影像建筑区域的建筑图斑与属性清洗后的POI数据进行落点分析,获取具有POI落点的建筑图斑的建筑土地利用属性和特征;将具有建筑土地利用属性和特征的建筑图斑输入随机森林分类器中进行训练,并通过随机森林分类器对不具有POI落点的建筑图斑进行预测,使不具有POI落点的建筑图斑获得建筑土地利用属性。本发明实现土地利用属性的空间迁移,提高了遥感影像土地利用属性的赋值效率。
Description
技术领域
本发明涉及一种遥感影像图像处理技术领域,尤其是涉及一种遥感影像建筑区域土地利用属性空间迁移方法。
背景技术
随着人口的急剧增长,可利用的土地资源越来也少,土地利用问题也逐渐被人们所重视。在遥感技术领域,土地利用一直是遥感影像解译的重要工作,也是资源合理分配调度及市政规划的重要参考依据。与传统的通过每十年一次的土地调查来探知全国土地利用情况相比,通过遥感影像解译土地利用情况可为全国土地调查工作提供可靠高效的解决方案。
然而,现有技术中,在利用遥感影像解译土地利用时还通过人工对遥感影像中的土地利用属性进行判读,并进一步对土地利用属性进行赋值,然而人工进行土地利用属性的赋值存在如下问题:(1)人工进行土地利用属性的赋值需耗费大量的人力物力,且效率低下;(2)人工进行土地利用属性的赋值的过程中容易出错。
发明内容
本发明的目的在于克服现有技术的缺陷,提供一种遥感影像建筑区域土地利用属性空间迁移方法,提高土地利用属性赋值效率。
为实现上述目的,本发明提出如下技术方案:一种遥感影像建筑区域土地利用属性空间迁移方法,包括如下步骤:
S100,获取遥感影像建筑区域的至少一种POI数据,并根据所述POI数据进行建筑属性清洗,使每个建筑图斑对应唯一的建筑属性;
S200,将建筑属性与土地调查获得的土地利用属性建立映射关系,使每个建筑属性对应建筑土地利用属性;
S300,获取遥感影像建筑区域已有的土地覆盖图层中的建筑图斑,并根据所述建筑图斑与属性清洗后的POI数据进行落点分析,获取具有POI落点的建筑图斑的建筑土地利用属性,进一步获取所述具有POI落点的建筑图斑的特征;
S400,将具有建筑土地利用属性和特征的建筑图斑输入随机森林分类器中进行训练,并通过随机森林分类器对不具有POI落点的建筑图斑进行预测,使不具有POI落点的建筑图斑获得建筑土地利用属性。
优选地,在步骤S100中,根据所述POI数据进行建筑属性清洗包括:
建立属性关键词样本库;
计算待匹配的建筑属性与属性关键词样本库中每个属性关键词的相似度;
选取相似度最大的待匹配的建筑属性作为建筑图斑的建筑属性。
优选地,通过如下步骤计算相似度:
将待匹配的建筑属性及属性关键词转换为向量,并根据如下公式计算相似度:
其中,A为待匹配的建筑属性转换后的向量,B为属性关键词转换后的向量。
优选地,通过词袋模型One-hot文本表示方法将待匹配的建筑属性和属性关键词转换为向量。
优选地,在步骤S300中,具有POI落点的建筑图斑的建筑土地利用属性通过如下步骤获得:
将POI数据转换为POI点图层;
将POI点图层与土地覆盖图层进行匹配处理,并在匹配过程中获取落在建筑图斑内的POI点,进一步通过最邻近算法对落在建筑图斑内的POI点进行筛选,将筛选出的POI点对应的建筑土地利用属性作为建筑图斑的建筑土地利用属性。
优选地,具有POI落点的建筑图斑的特征通过如下步骤获取:
针对根据遥感影像的建筑特点进行特征设计;
根据特征设计提取建筑图斑的特征。
优选地,还包括S500,将随机森林分类器的训练结果与验证集进行匹配,验证建筑土地利用属性迁移精度。
本发明的有益效果是:
本发明所述的遥感影像建筑区域土地利用属性空间迁移方法,通过POI数据获取建筑区域内具有POI落点的建筑图斑的建筑属性和特征,进一步通过具有建筑属性和特征的建筑图斑训练随机森林分类器,使随机森林分类器预测不具有POI落点的建筑图斑的建筑属性,最终实现土地利用属性的空间迁移,提高了遥感影像土地利用属性的赋值效率。
附图说明
图1是本发明的流程图示意图;
图2是本发明的相似度计算流程图示意图;
图3是本发明的遥感影像建筑区域部分影像图;
图4是图3中土地利用属性空间迁移前的影像图;
图5是图4中土地利用属性空间迁移后的影像图。
具体实施方式
下面将结合本发明的附图,对本发明实施例的技术方案进行清楚、完整的描述。
本发明所揭示的一种遥感影像建筑区域土地利用属性空间迁移方法,通过POI(Point Of Interest,兴趣点)数据获取建筑区域内具有POI落点的建筑图斑的建筑属性和特征,进一步通过具有建筑属性和特征的建筑图斑训练随机森林分类器(RF,RandomForest),使随机森林分类器预测不具有POI落点的建筑图斑的建筑属性,最终实现土地利用属性的空间迁移,提高了遥感影像土地利用属性的赋值效率。
如图1所示,为本发明所揭示的一种遥感影像建筑区域土地利用属性空间迁移方法,包括如下步骤:
S100,获取遥感影像建筑区域的至少一种POI数据,并对所述POI数据进行建筑属性清洗,使每个建筑图斑对应唯一的建筑属性;
具体地,在地里信息系统中,每个POI(Point Of Interest,兴趣点)通常包括名称、类别、经度、纬度信息,其可以是一个商店、酒吧、加油站等。实施时,POI数据可通过主流地图服务商获取,如通过网络爬虫从百度地图、高德地图、腾讯地图中的一个或多个获取POI数据。
本实施例中,以通过网络爬虫分别从百度地图、高德地图和腾讯地图获取POI数据为例,对如何进行建筑属性清洗进行详细的说明。
通过网络爬虫百度地图、高德地图和腾讯地图获取三种POI数据,POI数据包括建筑属性信息,如商店、酒吧、加油站等,则遥感影像建筑区域已有的土地覆盖图层中的每个建筑图斑具有三个临时建筑属性,如小学、中学和大学,可通过关键词匹配的方法对多个临时建筑属性进行清洗,使每个建筑图斑对应唯一的建筑属性,也即建筑区域的POI数据一般都是建筑相关的属性,如团体、机关、社等等,将这些关键字作为语义因子,并通过字符串匹配的分词方法将三种临时建筑属性整合为唯一的一种建筑属性,如将小学、中学和大学整合为学校。
进一步地,通过关键字匹配方法对POI数据进行建筑属性清洗包括如下步骤:首先,建立属性关键词样本库;其次,计算待匹配的建筑属性与属性关键词样本库中每个建筑属性的相似度;最后,选取相似度最大的待匹配的建筑属性作为建筑图斑的建筑属性。
实施时,首先将待匹配的建筑属性转换为向量,同时将属性关键词样本库中的建筑属性也转换为向量,进一步根据如下公式计算两个词的相似度:
其中,A为待匹配的建筑属性转换后的向量,B为属性关键词转换后的向量。
余弦值可表征两个向量的相似性,夹角越小,余弦值越接近于1,两个向量的方向更加吻合,也即两个词相似度高。
具体实施时,可通过词袋模型One-hot(独热编码)文本表示方法将待匹配的关键词转换为向量,及通过词袋模型One-hot文本表示方法将属性关键词样本库中的属性关键词转换为向量。
步骤S200,将建筑属性与土地调查获得的土地利用属性建立映射关系,使每个建筑属性对应建筑土地利用属性;
具体地,将步骤S100获得的建筑属性与全国第三次土地调查二级类属性进行映射,如学校对应科教文卫用地,编号为08H2等。具体实施时,并不限于与全国第三次土地调查二级类属性进行映射,可根据实际需求进行映射,如与全国第二次土地调查二级类属性进行映射等。
步骤S300,获取遥感影像建筑区域已有的土地覆盖图层中的建筑图斑,并根据所述建筑图斑与属性清洗后的POI数据进行落点分析,获取具有POI落点的建筑图斑的建筑土地利用属性,进一步获取所述具有POI落点的建筑图斑的特征。
具体地,建筑图斑可从遥感影像建筑区域已有的土地覆盖图层中提取。进一步地,具有POI落点的建筑图斑的建筑土地利用属性通过如下步骤获得:
首先,将属性清洗后的POI数据转换为POI点图层;
其次,将POI点图层与土地覆盖图层进行匹配处理,并在匹配过程中获取落在建筑图斑内的POI点,通过最邻近算法对落在建筑图斑内的POI点进行筛选,将筛选出的POI点对应的建筑土地利用属性作为建筑图斑的建筑土地利用属性。对于未落入POI点的建筑图斑,则等待土地利用属性空间迁移。
遥感影像建筑区域的建筑具有特定的建筑特点,为获取具有POI落点的建筑图斑的特征,首先根据遥感影像建筑区域的建筑特点进行特征设计,以从遥感影像中识别相应的建筑,最后,根据特征设计结果获取建筑图斑的特征。本实施例中,根据建筑特点设计了包括但不限于图像特征、矢量特征和空间关系特征,图像特征包括建筑图斑的标准差、对比度、熵、相异性和相关性,矢量特征包括形状指数、主方向、向量长、向量宽和向量长宽比;空间关系特征包括相邻边界、相对边界比例和城郊率。
其中,标准差通过如下公式计算:
Pi为像素值,n为像素个数,l为波段,u为均值;
对比度通过如下公式计算:
相异性通过如下公式计算:
熵通过如下公式计算:
相关性通过如下公式计算:
在对比度、相异性、熵、相关性公式中,Pi,j为归一化灰度共生矩阵GLCM的值,i,j为GLCM的行、列值,N是GLCM的行或列的数目;
形状指数通过如下公式计算:
其中,Pe为矢量图形的周长,Area为矢量图形的面积;
主方向通过如下公式计算:
其中,λ为特征值,VarX和VarY分别为X,Y的方差;
长度和宽度分别通过如下公式计算:
其中,Area为矢量图形的面积,λ为矢量图形的长宽比;
长度和宽度分别通过如下公式计算:
λ=λ1/λ2
其中,λ1为方差矩阵特征值的较大值,λ2为方差矩阵特征值的较小值;
相邻边界通过如下公式计算:
其中,b(u,v)为对象u和v共同边界的长度;
相邻边界比例通过如下公式计算:
其中,b(u,v)为对象u和v共同边界的长度,Pe(v)为对象v的边界总长度;
城郊率通过如下公式计算:
ratio=dcity/rcity
其中,dcity为城市距离,rcity为城市半径。
具体实施时可根据上述图像特征、矢量特征或空间特征计算获取建筑图斑的特征。
S400,将具有建筑土地利用属性和特征的建筑图斑输入随机森林分类器中进行训练,并通过随机森林分类器对不具有POI落点的建筑图斑进行预测,使不具有POI落点的建筑图斑获得建筑土地利用属性。
具体地,将多个具有建筑土地利用属性和特征的建筑图斑建立样本库。由于建筑图斑的不规则形状,则将样本库输入随机分类器中进行训练,并通过训练后的随机森林分类器对不具有POI落点的建筑图斑进行预测,使不具有POI落点的建筑图斑获得建筑土地利用属性,也即完成土地利用属性空间属性的迁移的同时提高预测的准确性。
当完成土地利用属性空间迁移后,还需进行精度评定,以验证随机森林分类器预测结果的准确性。实施时,可随机森林分类器的训练结果与验证集进行匹配验证土地利用属性空间迁移的精度。当不满足精度要求时,继续使用样本库训练随机森林分类器,直至土地利用属性空间迁移的精度满足要求。
本发明所述的遥感影像建筑区域土地利用属性空间迁移方法,有效提高了土地利用属性的赋值效率。如图3~5所示,图3为遥感影像建筑区域部分影像图,图4为图3中土地利用属性空间迁移前的影像图,该图中包含多个不具有POI落点的建筑图斑;图5为利用本发明的方法进行土地利用属性空间迁移后的示意图,由图可知,不具有POI落点的建筑图斑均具有了建筑土地利用属性。
同时,本发明还通过选取一总体面积约390平方千米的区域对本发明土地利用属性空间迁移方法进行验证。通过本发明所述的方法,仅需约6个小时便完成建筑区域土地利用属性赋值,而通过人工进行建筑区域土地利用属性赋值,则需要约60个小时完成建筑区域土地利用属性赋值,由此可见本发明大大缩短了建筑区域土地利用属性赋值效率。
本发明的技术内容及技术特征已揭示如上,然而熟悉本领域的技术人员仍可能基于本发明的教示及揭示而作种种不背离本发明精神的替换及修饰,因此,本发明保护范围应不限于实施例所揭示的内容,而应包括各种不背离本发明的替换及修饰,并为本专利申请权利要求所涵盖。
Claims (7)
1.一种遥感影像建筑区域土地利用属性空间迁移方法,其特征在于,包括如下步骤:
S100,获取遥感影像建筑区域的至少一种POI数据,并根据所述POI数据进行建筑属性清洗,使每个建筑图斑对应唯一的建筑属性;
S200,将建筑属性与土地调查获得的土地利用属性建立映射关系,使每个建筑属性对应建筑土地利用属性;
S300,获取遥感影像建筑区域已有的土地覆盖图层中的建筑图斑,并根据所述建筑图斑与属性清洗后的POI数据进行落点分析,获取具有POI落点的建筑图斑的建筑土地利用属性,进一步获取所述具有POI落点的建筑图斑的特征;
S400,将具有建筑土地利用属性和特征的建筑图斑输入随机森林分类器中进行训练,并通过随机森林分类器对不具有POI落点的建筑图斑进行预测,使不具有POI落点的建筑图斑获得建筑土地利用属性。
2.根据权利要求1所述的方法,其特征在于,在步骤S100中,根据所述POI数据进行建筑属性清洗包括:
建立属性关键词样本库;
计算待匹配的建筑属性与属性关键词样本库中每个属性关键词的相似度;
选取相似度最大的待匹配的建筑属性作为建筑图斑的建筑属性。
3.根据权利要求2所述的方法,其特征在于,通过如下步骤计算相似度:
将待匹配的建筑属性及属性关键词转换为向量,并根据如下公式计算相似度:
其中,A为待匹配的建筑属性转换后的向量,B为属性关键词转换后的向量。
4.根据权利要求3所述的方法,其特征在于,通过词袋模型One-hot文本表示方法将待匹配的建筑属性和属性关键词转换为向量。
5.根据权利要求1所述的方法,其特征在于,在步骤S300中,具有POI落点的建筑图斑的建筑土地利用属性通过如下步骤获得:
将属性清洗后的POI数据转换为POI点图层;
将POI点图层与土地覆盖图层进行匹配处理,并在匹配过程中获取落在建筑图斑内的POI点,进一步通过最邻近算法对落在建筑图斑内的POI点进行筛选,将筛选出的POI点对应的建筑土地利用属性作为建筑图斑的建筑土地利用属性。
6.根据权利要求1所述的方法,其特征在于,具有POI落点的建筑图斑的特征通过如下步骤获取:
针对根据遥感影像的建筑特点进行特征设计;
根据特征设计提取建筑图斑的特征。
7.根据权利要求1所述的方法,其特征在于,还包括S500,将随机森林分类器的训练结果与验证集进行匹配,验证建筑土地利用属性迁移精度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910743346.3A CN110472559B (zh) | 2019-08-13 | 2019-08-13 | 一种遥感影像建筑区域土地利用属性空间迁移方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910743346.3A CN110472559B (zh) | 2019-08-13 | 2019-08-13 | 一种遥感影像建筑区域土地利用属性空间迁移方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110472559A true CN110472559A (zh) | 2019-11-19 |
CN110472559B CN110472559B (zh) | 2023-06-23 |
Family
ID=68511750
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910743346.3A Active CN110472559B (zh) | 2019-08-13 | 2019-08-13 | 一种遥感影像建筑区域土地利用属性空间迁移方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110472559B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112925773A (zh) * | 2019-12-10 | 2021-06-08 | 中国再保险(集团)股份有限公司 | 一种构建行业风险暴露数据库的poi数据清洗与融合方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001075790A2 (en) * | 2000-04-03 | 2001-10-11 | 3-Dimensional Pharmaceuticals, Inc. | Method, system, and computer program product for representing object relationships in a multidimensional space |
CN107247938A (zh) * | 2017-06-08 | 2017-10-13 | 中国科学院遥感与数字地球研究所 | 一种高分辨率遥感影像城市建筑物功能分类的方法 |
CN108038448A (zh) * | 2017-12-13 | 2018-05-15 | 河南理工大学 | 基于加权熵的半监督随机森林高光谱遥感影像分类方法 |
-
2019
- 2019-08-13 CN CN201910743346.3A patent/CN110472559B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001075790A2 (en) * | 2000-04-03 | 2001-10-11 | 3-Dimensional Pharmaceuticals, Inc. | Method, system, and computer program product for representing object relationships in a multidimensional space |
CN107247938A (zh) * | 2017-06-08 | 2017-10-13 | 中国科学院遥感与数字地球研究所 | 一种高分辨率遥感影像城市建筑物功能分类的方法 |
CN108038448A (zh) * | 2017-12-13 | 2018-05-15 | 河南理工大学 | 基于加权熵的半监督随机森林高光谱遥感影像分类方法 |
Non-Patent Citations (1)
Title |
---|
赵海莹等: "遥感影像更新土地利用数据库的方法与流程", 《城市勘测》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112925773A (zh) * | 2019-12-10 | 2021-06-08 | 中国再保险(集团)股份有限公司 | 一种构建行业风险暴露数据库的poi数据清洗与融合方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN110472559B (zh) | 2023-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Song et al. | Faster region convolutional neural network for automated pavement distress detection | |
CN108764193B (zh) | 融合poi和遥感影像的城市功能区划分方法 | |
CN110348506B (zh) | 基于遥感图像的土地利用分类方法、存储介质和计算设备 | |
Lloyd | Spatial data analysis: an introduction for GIS users | |
Qi et al. | Knowledge discovery from soil maps using inductive learning | |
Rashed et al. | Measuring temporal compositions of urban morphology through spectral mixture analysis: toward a soft approach to change analysis in crowded cities | |
Wu et al. | Using geometrical, textural, and contextual information of land parcels for classification of detailed urban land use | |
Zhang et al. | Heuristic sample learning for complex urban scenes: Application to urban functional-zone mapping with VHR images and POI data | |
Liu et al. | Road selection based on Voronoi diagrams and “strokes” in map generalization | |
CN105912643A (zh) | 一种基于内容改进型均值哈希的图像检索方法 | |
de Souza Pereira et al. | Heterogeneity and homogeneity of Southern Brazilian beaches: A morphodynamic and statistical approach | |
CN110399819A (zh) | 一种基于深度学习的遥感影像居民区自动提取方法 | |
Gibril et al. | New semi-automated mapping of asbestos cement roofs using rule-based object-based image analysis and Taguchi optimization technique from WorldView-2 images | |
CN108280396A (zh) | 基于深度多特征主动迁移网络的高光谱图像分类方法 | |
Zhang et al. | How do land-use patterns influence residential environment quality? A multiscale geographic survey in Beijing | |
CN111814596B (zh) | 一种遥感影像与出租车轨迹融合的自动城市功能分区方法 | |
CN106846322A (zh) | 基于曲线波滤波器和卷积结构学习的sar图像分割方法 | |
Rosier et al. | Fusing Earth observation and socioeconomic data to increase the transferability of large-scale urban land use classification | |
Piyoosh et al. | Semi-automatic mapping of anthropogenic impervious surfaces in an urban/suburban area using Landsat 8 satellite data | |
Staab et al. | Predicting traffic noise using land-use regression—a scalable approach | |
Nduwayezu et al. | Urban growth and land use/land cover changes in the post-genocide period, Kigali, Rwanda | |
Zhang et al. | Using street view images to identify road noise barriers with ensemble classification model and geospatial analysis | |
Ai et al. | Construction of small sample seismic landslide susceptibility evaluation model based on Transfer Learning: a case study of Jiuzhaigou earthquake | |
Kundu et al. | Urban change detection analysis during 1978–2017 at Kolkata, India, using multi-temporal satellite data | |
Jianwen et al. | Land-use classification using ASTER data and self-organized neutral networks |
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 |