一种林业图斑变化快速识别方法及系统技术方案

本发明专利技术属于图斑监测技术领域，本发明专利技术提供了一种林业图斑变化快速识别方法及系统，通过获取首个待测区域的遥感图像，根据图斑范围将所述遥感图像划分为多个图斑图像；继续接收预设时间段内的待测区域的遥感图像，并依据各个接收到的遥感图像计算各个图斑图像内的空间变化度；实时判断各个空间变化度是否超过变化阈值，如果是则将该图斑图像标记为异常图斑。根据空间变化度能够准确的监控到最近的图斑图像的邻接区域的细微的变化，能够准确的防止不明显的恶意占地，只是从遥感图像上监测位置变化，极大的减少了识别的误差率，平衡了图斑边界识别的灵敏性和地图上图斑边界变化的识别的准确性。别的准确性。别的准确性。

【技术实现步骤摘要】
一种林业图斑变化快速识别方法及系统


[0001]本专利技术属于图斑监测
，具体涉及一种林业图斑变化快速识别方法及系统。

技术介绍

[0002]林业图斑是指利用卫星遥感技术，获得某个时间内所监测的林业土地范围内用地情况的图像，结合GIS地理信息系统技术，给所监测的林业土地用一定的规则进行分区编号，而这个编号就称为图斑。利用卫星遥感技术进行林业土地资源监测越来越普及，而林业图斑在划定的过程中由于受到了诸多人为因素的影响，忽略了质量指标以及空间分布的合理性，存在林业土地和普通的工业、农业用地重叠的情况，
[0003]目前，结合GIS地理信息系统技术获得的林业图斑，将林业土地和普通的工业、农业用地进行空间叠加相交分析，从而获得两图层间重叠的林业图斑。现有技术，例如公开号为CN111797190A的专利技术专利提供了一种图斑提取方法及装置，依据约束条件选取目标重叠图斑，筛选因子满足约束条件的重叠图斑。将形状、面积、地类区域、以及地类区域的面积占比作为重叠图斑的筛选因子，依据约束条件，针对性地从各个重叠图斑中筛选出目标重叠图斑，但是，通过空间叠加相交分析之后，会产生碎小图斑、图斑形状狭长、图斑数量多等情况，无法识别出一些逐渐恶意占地，或者重叠边界的图斑区域，尤其是在一些区域经过时间较长的开发或者人为改变导致地形地貌发生改变、或者自然灾害、地质活动频繁的区域、人为恶意占地导致林业图斑无法准确的识别。

技术实现思路

[0004]本专利技术的目的在于提出一种林业图斑变化快速识别方法及系统，以解决现有技术中所存在的一个或多个技术问题，至少提供一种有益的选择或创造条件。
[0005]为了实现上述目的，根据本专利技术的一方面，提供一种林业图斑变化快速识别方法，所述方法包括以下步骤：
[0006]S100，获取首个待测区域的遥感图像，根据图斑范围将所述遥感图像划分为多个图斑图像；
[0007]S200，继续接收预设时间段内的待测区域的遥感图像，并依据各个接收到的遥感图像计算各个图斑图像内的空间变化度；
[0008]S300，实时判断各个空间变化度是否超过变化阈值，如果是则将该图斑图像标记为异常图斑；
[0009]S400，将所有异常图斑影像输出到客户端进行显示。
[0010]进一步地，在S100中，所述待测区域为林地，所述林地是生长有乔木、竹类、灌木、沿海红树林的土地,所述林地不包括居民绿化用地,以及铁路、公路、河流沟渠的护路、护草林。
[0011]进一步地，在S100中，根据图斑范围将所述遥感图像划分为多个图斑图像的方法
为：将地貌、土地利用类型相同、水土流失类型一致的林业土地地块分为一类，以该一类林业土地地块作为图斑范围，以图斑范围内的林业土地地块作为基础调查单元，然后将图斑范围勾绘到遥感图像上成为图斑；(其中，遥感图像上的最小图斑面积不小于0.5cm2，实地面积0.5km2；最大不超过50cm2，实地面积50km2)，通过图像分割或者矢量图套合，将整个遥感图像根据图斑划分成若干个小的区域，以所述区域作为图斑图像。
[0012]进一步地，在S200中，接收预设时间段内的待测区域的遥感图像为：以1～3天为时间间隔在预设时间段中进行采集待测区域的遥感图像；所述预设时间段为最近120～360天的时间，以1～3天为时间间隔进行采集待测区域的遥感图像；所述遥感图像通过无人机或卫星的机载遥感设备采集待测区域的遥感图像，所述机载遥感设备为SAR遥感设备。
[0013]进一步地，在S200中，依据各个接收到的遥感图像计算各个图斑图像内的空间变化度的方法为：
[0014]以各个接收到的遥感图像构成的集合记为Image＝{Img
i1
},Img
i1
是集合Image中的第i1个遥感图像，i1是集合Image中遥感图像的序号，Sum1是集合Image中遥感图像的总数；
[0015]取各个图斑图像在遥感图像上的对应区域的集合作为监测集合Suv＝{Suvy
i2
}，Suvy
i2
是集合Suv中的第i2个图斑图像在遥感图像上的对应区域，i2是集合Suv中图斑图像在遥感图像上的对应区域的序号，Sum2是监测集合Suv中元素的总数；
[0016]将Image中所有遥感图像Img
i1
灰度化后经过边缘检测得到多个边缘线，由各个边缘线构成多个封闭区域；筛选出各个封闭区域中的封闭区域在遥感图像上的对应区域的几何重心和Suvy
i2
的几何重心之间的欧氏距离最短的封闭区域记为Suvy
i2
在遥感图像Img
i1
上的映射区域Suvy
i1,i2
，则Suvy
i2
在遥感图像Img
i1
上有一个对应的映射区域Suvy
i1,i2
；
[0017]在i1∈[2,Sum1
‑
1]、i2∈[2,Sum2
‑
1]的范围内，遍历映射区域Suvy
i1,i2
，如果首次遍历到映射区域Suvy
i1,i2
、Suvy
i1
‑
1,i2
和Suvy
i1+1,i2
之间存在补集，则以Suvy
i1,i2
的几何重心为锚定点Anchor，以补集区域的边缘线上的各个点中与Anchor之间距离最大的点为Pmax，以补集区域的边缘线上的各个点中与Anchor之间距离最小的点为Pmin，同时记对应的遥感图像Img
i1
的采集时间为T1，以T1到当前时间T0之间的时间段为空间变化时间tG；(补集区域是映射区域Suvy
i1,i2
、Suvy
i1
‑
1,i2
和Suvy
i1+1,i2
之间做求补集运算得到的区域)；
[0018]其中，Suvy
i1
‑
1,i2
和Suvy
i1+1,i2
分别是Suvy
i2
在遥感图像Img
i1
‑1上的映射区域和Suvy
i2
在遥感图像Img
i1+1
上的映射区域；
[0019]对于在一些区域经过时间较长的开发或者人为改变导致地形地貌发生改变、或者自然灾害、地质活动频繁的区域、人为恶意占地发生变化后的待测区域，由于已经监测到地图的图斑区域(图斑也可以为根据行政界线，土地权属界线或线状地物所分割的地块)由于遥感图像中图斑边界发生了变化，所以在它们之间存在补集后(即图斑图像的邻接区域Suvy
i1,i2
、Suvy
i1
‑
1,i2
和Suvy
i1+1,i2
的边界发生变化后)，点Pmax和点Pmin之间的线段的长度会随时间发生变化；
[0020]以点P本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种林业图斑变化快速识别方法，其特征在于，所述方法包括以下步骤：S100，获取首个待测区域的遥感图像，根据图斑范围将所述遥感图像划分为多个图斑图像；S200，继续接收预设时间段内的待测区域的遥感图像，并依据各个接收到的遥感图像计算各个图斑图像内的空间变化度；S300，实时判断各个空间变化度是否超过变化阈值，如果是则将该图斑图像标记为异常图斑；S400，将所有异常图斑影像输出到客户端进行显示。2.根据权利要求1所述的一种林业图斑变化快速识别方法，其特征在于，在S100中，根据图斑范围将所述遥感图像划分为多个图斑图像的方法为：将地貌、土地利用类型相同、水土流失类型一致的林业土地地块分为一类，以该一类林业土地地块作为图斑范围，以图斑范围内的林业土地地块作为基础调查单元，然后将图斑范围勾绘到遥感图像上成为图斑，通过图像分割或者矢量图套合，将整个遥感图像根据图斑划分成若干个小的区域，以所述区域作为图斑图像。3.根据权利要求1所述的一种林业图斑变化快速识别方法，其特征在于，在S200中，接收预设时间段内的待测区域的遥感图像为：以1～3天为时间间隔在预设时间段中进行采集待测区域的遥感图像；所述预设时间段为最近120～360天的时间，以1～3天为时间间隔进行采集待测区域的遥感图像；所述遥感图像通过无人机或卫星的机载遥感设备采集待测区域的遥感图像，所述机载遥感设备为SAR遥感设备。4.根据权利要求1所述的一种林业图斑变化快速识别方法，其特征在于，在S200中，依据各个接收到的遥感图像计算各个图斑图像内的空间变化度的方法为：以各个接收到的遥感图像构成的集合记为Image＝{Img
i1
},Img
i1
是集合Image中的第i1个遥感图像，i1是集合Image中遥感图像的序号，Sum1是集合Image中遥感图像的总数；取各个图斑图像在遥感图像上的对应区域的集合作为监测集合Suv＝{Suvy
i2
}，Suvy
i2
是集合Suv中的第i2个图斑图像在遥感图像上的对应区域，i2是集合Suv中图斑图像在遥感图像上的对应区域的序号，Sum2是监测集合Suv中元素的总数；将Image中所有遥感图像Img
i1
灰度化后经过边缘检测得到多个边缘线，由各个边缘线构成多个封闭区域；筛选出各个封闭区域中的封闭区域在遥感图像上的对应区域的几何重心和Suvy
i2
的几何重心之间的欧氏距离最短的封闭区域记为Suvy
i2
在遥感图像Img
i1
上的映射区域Suvy
i1,i2
，则Suvy
i2
在遥感图像Img
i1
上有一个对应的映射区域Suvy
i1,i2
；在i1∈[2,Sum1
‑
1]、i2∈[2,Sum2
‑
1]的范围内，遍历映射区域Suvy
i1,i2
，如果首次遍历到映射区域Suvy
i1,i2
...

【专利技术属性】
技术研发人员：骆水灵，彭威雄，卓国军，梁志明，袁东，赖娇娇，袁超岸，邓学启，袁桂鸿，钟昌栋，袁倩雯，
申请(专利权)人：广州市绿之城园林绿化工程有限公司，
类型：发明
国别省市：