利用互联网地图数据进行地表覆盖变化检测的方法及系统技术方案

本发明专利技术提出了利用互联网地图数据进行地表覆盖变化检测的方法及系统，首先，依照需要进行变化检测的地表覆盖产品的分辨率，划分地表覆盖区域的栅格单元；继而，提取每个栅格单元的中互联网地图数据的文本信息，利用层次狄利克雷过程算法，对每个栅格单元进行地表覆盖分类；最后，将分类后每个栅格单元的地表覆盖类型与需要进行变化检测的地表覆盖产品的类型进行对比，以达到地表覆盖变化检测的目的。

【技术实现步骤摘要】
利用互联网地图数据进行地表覆盖变化检测的方法及系统
本专利技术涉及利用互联网地图数据进行地表覆盖变化检测的方法及系统。
技术介绍
地表覆盖及变化是环境变化研究、地理世情监测、可持续发展规划等不可或缺的重要基础信息和关键参量，对于研究地球系统的能量平衡、碳循环及其他生物地球化学循环、气候变化等有着十分重要的意义。传统的地表覆盖变化检测方法主要基于遥感影像，通过对比不同时期遥感影像中的光谱信息，提取地表覆盖类型及其变化信息。该方法耗时耗力、生产周期长，无法满足地表覆盖产品快速更新的需求。因此，需要发展一种新的地表覆盖变化快速检测方法。目前，已经出现了诸多应用众源地理数据进行地表覆盖变化检测的方法。例如，“Johnson,BrianA.,andKotaroIizuka."IntegratingOpenStreetMapcrowdsourceddataandLandsattime-seriesimageryforrapidlanduse/landcover(LULC)mapping:CasestudyoftheLagunadeBayareaofthePhilippines."AppliedGeography67(2016):140-149.”中利用OpenStreetMap中的建筑物、水体、绿地数据，作为地表覆盖变化的参考数据，并用于遥感数据中地表覆盖类型的变化检测。“Fonte,CidáliaCosta,etal."GeneratingUp-to-DateandDetailedLandUseandLandCoverMapsUsingOpenStreetMapandGlobeLand30."ISPRSInternationalJournalofGeo-Information6.4(2017):125.”中将所有OSM中的数据类型分为不同的地表覆盖类型，对GlobeLand30中2010年的地表覆盖数据进行变化检测。上述研究通过利用众源地理数据，实现了地表覆盖的快速变化检测。然而，由于众源地理数据大多来自未经过训练的志愿者，其数据质量无法保证，并且其覆盖范围为志愿者主要活动区域，数据更新也无法保证。因此，将其做为数据源进行地表覆盖变化检测，在一定程度上影响了精度。近年来，随着互联网地图公司的兴起，以及商用互联网地图的广泛应用，出现了诸如百度地图、高德地图等互联网地图数据。这些数据来源于专业的商业公司，其精度高、更新实时、覆盖面广，也包含了大量能够反映地表覆盖类型的信息。因此，将商用互联网地图数据应用于地表覆盖变化检测是一条可行途径。
技术实现思路
本专利技术的目的就是为了解决上述问题，提供利用互联网地图数据进行地表覆盖变化检测的方法及系统，由于互联网数据包含了点、线、面等多种不同的数据类型，同时其文本信息类型多样，无法直接应用于地表覆盖变化检测。为了实现上述目的，本专利技术采用如下技术方案：利用互联网地图数据进行地表覆盖变化检测的方法，包括如下步骤：步骤A：将地表覆盖区域划分为若干个栅格单元，利用每个栅格单元中内互联网地图数据的位置信息，将互联网地图数据与划分后的栅格单元进行空间位置叠加；步骤B：对互联网地图数据中的文本信息进行分类，使互联网地图数据的文本信息与地表覆盖类型相对应；步骤C：对每个栅格单元进行地表覆盖类型的变化检测。优选的，所述步骤A中地表覆盖区域的栅格单元划分：依据需要进行变化检测的地表覆盖产品的分辨率，将地表覆盖区域划分n个与分辨率相同的栅格单元Ri，j，其中，R为划分的栅格单元，i为栅格单元R的行数，j为栅格单元R的列数。优选的，所述步骤B中，对互联网地图数据中的文本信息进行分类：提取每个栅格单元中互联网地图数据的文本信息，利用层次狄利克雷过程算法，对每个栅格单元内的文本信息进行自动分类。优选的，所述步骤C：首先，依据需要进行变化检测的地表覆盖产品，提取每个栅格单元中包含的原有地表覆盖类型Cl；继而，依据步骤B中的文本信息分类结果，提取每个栅格单元内文本信息类型Ct；如果Cl＝Ct，则保留需要进行变化检测的地表覆盖产品的类型；如果Cl≠Ct，则保留互联网地图数据的文本信息类型，作为新的地表覆盖类型，去除原有地表覆盖产品的类型，从而完成地表覆盖的变化检测。步骤B例如，互联网地图数据中的“购物商场”归类为地表覆盖中的“人造地表”类型，“高尔夫球场”归类为“草地”。所述步骤B包括：步骤B1：提取每个栅格单元中互联网地图数据的文本信息，将每个栅格单元内的文本信息作为一个文档di；所有文档定义为D∈{d1，d2，d3....dn}，其中，n为所有地表覆盖栅格单元的数量；步骤B2：利用层次狄利克雷过程算法，对每个栅格单元内的文本信息进行自动分类，使其与地表覆盖类型相对应：利用层次狄利克雷过程算法计算得到每种地表覆盖类型的分类概率Gj，取概率最大的类型为最终的分类类型：G0～DP(γ，H)；(1)Gj～DP(α，G0)；(2)其中，G0、Gj为两个狄利克雷过程，H为假设的度量空间的一个分布，γ、α分别为满足公式(1)(2)中狄利克雷分布的两个参数，由经验值取得；其中，狄利克雷过程G的计算公式如下：其中，Dir(·)为狄利克雷分布，计算公式为：其中，B(α)为对于参数α的B(·)函数，计算公式为：其中，Γ(·)为gamma函数，表示为αi的阶乘。利用互联网地图数据进行地表覆盖变化检测的系统，包括：叠加模块：将地表覆盖区域划分为若干个栅格单元，利用每个栅格单元中内互联网地图数据的位置信息，将互联网地图数据与划分后的栅格单元进行空间位置叠加；分类模块：对互联网地图数据中的文本信息进行分类，使互联网地图数据的文本信息与地表覆盖类型相对应；检测模块：对每个栅格单元进行地表覆盖类型的变化检测。优选的，所述叠加模块中地表覆盖区域的栅格单元划分：依据需要进行变化检测的地表覆盖产品的分辨率，将地表覆盖区域划分n个与分辨率相同的栅格单元Ri，j，其中，R为划分的栅格单元，i为栅格单元R的行数，j为栅格单元R的列数。优选的，所述分类模块中，对互联网地图数据中的文本信息进行分类：提取每个栅格单元中互联网地图数据的文本信息，利用层次狄利克雷过程算法，对每个栅格单元内的文本信息进行自动分类。优选的，所述检测模块：首先，依据需要进行变化检测的地表覆盖产品，提取每个栅格单元中包含的原有地表覆盖类型Cl；继而，依据步骤B中的文本信息分类结果，提取每个栅格单元内文本信息类型Ct；如果Cl＝Ct，则保留需要进行变化检测的地表覆盖产品的类型；如果Cl≠Ct，则保留互联网地图数据的文本信息类型，作为新的地表覆盖类型，去除原有地表覆盖产品的类型，从而完成地表覆盖的变化检测。本专利技术的有益效果是：首先，依照需要进行变化检测的地表覆盖产品的分辨率，划分地表覆盖区域的栅格单元；继而，提取每个栅格单元的中互联网地图数据的文本信息，利用层次狄利克雷过程算法，对每个栅格单元进行地表覆盖分类；最后，将分类后每个栅格单元的地表覆盖类型与需要进行变化检测的地表覆盖产品的类型进行对比，以达到地表覆盖变化检测的目的。附图说明图1为地表覆盖变化检测方法示意图；图2为北京市部分区域的高德地图数据与地表覆盖栅格单元叠加示意图；图3(a)-图3(b)为北京市部分地区地表覆盖单元及地本文档来自技高网...

【技术保护点】
利用互联网地图数据进行地表覆盖变化检测的方法，其特征是，包括如下步骤：步骤A：将地表覆盖区域划分为若干个栅格单元，利用每个栅格单元中内互联网地图数据的位置信息，将互联网地图数据与划分后的栅格单元进行空间位置叠加；步骤B：对互联网地图数据中的文本信息进行分类，使互联网地图数据的文本信息与地表覆盖类型相对应；步骤C：对每个栅格单元进行地表覆盖类型的变化检测。

【技术特征摘要】
1.利用互联网地图数据进行地表覆盖变化检测的方法，其特征是，包括如下步骤：步骤A：将地表覆盖区域划分为若干个栅格单元，利用每个栅格单元中内互联网地图数据的位置信息，将互联网地图数据与划分后的栅格单元进行空间位置叠加；步骤B：对互联网地图数据中的文本信息进行分类，使互联网地图数据的文本信息与地表覆盖类型相对应；步骤C：对每个栅格单元进行地表覆盖类型的变化检测。2.如权利要求1所述的方法，其特征是，所述步骤A中地表覆盖区域的栅格单元划分：依据需要进行变化检测的地表覆盖产品的分辨率，将地表覆盖区域划分n个与分辨率相同的栅格单元Ri，j，其中，R为划分的栅格单元，i为栅格单元R的行数，j为栅格单元R的列数。3.如权利要求1所述的方法，其特征是，所述步骤B中，对互联网地图数据中的文本信息进行分类：提取每个栅格单元中互联网地图数据的文本信息，利用层次狄利克雷过程算法，对每个栅格单元内的文本信息进行自动分类。4.如权利要求1所述的方法，其特征是，所述步骤C：首先，依据需要进行变化检测的地表覆盖产品，提取每个栅格单元中包含的原有地表覆盖类型Cl；继而，依据步骤B中的文本信息分类结果，提取每个栅格单元内文本信息类型Ct；如果Cl＝Ct，则保留需要进行变化检测的地表覆盖产品的类型；如果Cl≠Ct，则保留互联网地图数据的文本信息类型，作为新的地表覆盖类型，去除原有地表覆盖产品的类型，从而完成地表覆盖的变化检测。5.如权利要求1所述的方法，其特征是，所述步骤B包括：步骤B1：提取每个栅格单元中互联网地图数据的文本信息，将每个栅格单元内的文本信息作为一个文档di；所有文档定义为D∈{d1，d2，d3....dn}，其中，n为所有地表覆盖栅格单元的数量；步骤B2：利用层次狄利克雷过程算法，对每个栅格单元内的文本信息进行自动分类，使其与地表覆盖类型相对应：利用层次狄利克雷过程算法计算得到每种地表覆盖类型的分类概率Gj，取概率最大的类型为最终的分类类型：G0～DP(γ，H)；(1)Gj～DP(α，G0)；(2)其中，G0、Gj为两个狄利克雷过程，H为假设的度量空间的一个分布，γ、α分别为满足公式(1)(2)中狄利克雷分布的两个参数，由经验值取得；其中，狄利克雷过程G的计算公式如下：

【专利技术属性】
技术研发人员：邢汉发，吴泉源，徐海滨，孟媛，秦嘉徽，王梓璇，
申请(专利权)人：山东师范大学，
类型：发明
国别省市：山东,37