CN106548141A - 一种基于三角网的面向对象耕地信息自动提取方法 - Google Patents
一种基于三角网的面向对象耕地信息自动提取方法 Download PDFInfo
- Publication number
- CN106548141A CN106548141A CN201610934973.1A CN201610934973A CN106548141A CN 106548141 A CN106548141 A CN 106548141A CN 201610934973 A CN201610934973 A CN 201610934973A CN 106548141 A CN106548141 A CN 106548141A
- Authority
- CN
- China
- Prior art keywords
- triangulation network
- cluster
- central point
- oriented
- farmland information
- 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
- 238000000605 extraction Methods 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 claims abstract description 49
- 230000011218 segmentation Effects 0.000 claims abstract description 42
- 238000005520 cutting process Methods 0.000 claims abstract description 41
- 238000010586 diagram Methods 0.000 claims abstract description 15
- 238000009826 distribution Methods 0.000 claims description 23
- 230000003595 spectral effect Effects 0.000 claims description 11
- 238000003709 image segmentation Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000012217 deletion Methods 0.000 claims description 3
- 230000037430 deletion Effects 0.000 claims description 3
- 239000004744 fabric Substances 0.000 claims description 2
- 238000002715 modification method Methods 0.000 claims description 2
- 239000000284 extract Substances 0.000 abstract description 7
- 230000008569 process Effects 0.000 abstract description 6
- 238000003973 irrigation Methods 0.000 abstract description 5
- 230000002262 irrigation Effects 0.000 abstract description 5
- 238000011156 evaluation Methods 0.000 abstract description 2
- 230000002452 interceptive effect Effects 0.000 abstract 1
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 description 10
- 241000196324 Embryophyta Species 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 238000012360 testing method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 150000002500 ions Chemical class 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 206010008190 Cerebrovascular accident Diseases 0.000 description 2
- 241001269238 Data Species 0.000 description 2
- 241001420622 Meris Species 0.000 description 2
- 208000006011 Stroke Diseases 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000003066 decision tree Methods 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 238000010191 image analysis Methods 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 241000132092 Aster Species 0.000 description 1
- NLZUEZXRPGMBCV-UHFFFAOYSA-N Butylhydroxytoluene Chemical compound CC1=CC(C(C)(C)C)=C(O)C(C(C)(C)C)=C1 NLZUEZXRPGMBCV-UHFFFAOYSA-N 0.000 description 1
- -1 Transact Ions Chemical class 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- JZCCFEFSEZPSOG-UHFFFAOYSA-L copper(II) sulfate pentahydrate Chemical compound O.O.O.O.O.[Cu+2].[O-]S([O-])(=O)=O JZCCFEFSEZPSOG-UHFFFAOYSA-L 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000009313 farming Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000000873 masking effect Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003012 network analysis Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000009331 sowing Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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/13—Satellite images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- 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
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20016—Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20021—Dividing image into blocks, subimages or windows
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30188—Vegetation; Agriculture
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (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)
- Astronomy & Astrophysics (AREA)
- Remote Sensing (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明涉及一种基于三角网的面向对象耕地信息自动提取方法，其步骤为：利用多尺度分割法对高空间分辨率影像进行分割；剔除长条状分割对象（道路，沟渠等）；提取剩余分割对象的中心点；利用中心点构建三角网；对三角网进行剥皮操作；使用AUTOCLUST聚类算法构建三角网，并聚类；利用V图约束，优化聚类结果，避免过度聚类和欠聚类；利用最大方差约束，剔除剩余的零星林地；得到耕地提取结果并进行精度评估。本发明克服了高空间分辨率遥感影像数据量大，处理难的问题，充分利用分割后的对象提供的语义信息，通过剔除居民地和道路等对耕地信息提取造成干扰的对象，从而高效地进行耕地信息的自动提取，并保证耕地提取的整体性。
Description
技术领域
本发明涉及一种耕地信息自动提取方法，特别是涉及一种面向对象的高效非监督耕地信息提取方法。
背景技术
耕地信息的自动提取对于我国农村土地资源管理和土地规划有着至关重要的作用。当前,耕地信息制图的相关技术方法主要以遥感技术为主。在不同的空间尺度上，已经有很多关于耕地信息提取的研究，主要分为局部、地区和全球区域。局部区域的研究主要涉及单个地区或灌区盆地或特定区域；地区范围主要包括大河流域和大陆性的区域；而全球尺度则主要涉及全球范围的耕地信息制图。在不同空间尺度下，为了达到最佳的耕地信息提取效果，不同的传感器数据，需要采用不同的制图方式，根据不少学者的研究成果总结，局部区域的耕地信息提取采用Landsat TM/ETM+，SPOT，LISS，ASTER，CBERS数据能够达到较好的性价比，信息提取方法主要涉及相片解译、影像数字运算、影像分类、分割、影像融合。地区范围主要采用Landsat TM/ETM+，MODIS，MERIS，AVHRR，SPOT VGT数据，通过时间序列分析方法、监督或非监督分类方法、掩膜方法识别耕地信息。全球尺度主要采用MODIS，MERIS，AVHRR，SPOT VGT数据，信息提取方法主要涉及非监督聚类、针对时间序列数据的机器学习算法，以及采用其它辅助数据(包括统计数据，地面实时数据等)。
以上研究使用的遥感影像的空间分辨率较低，当识别小范围和重要地区的耕地往往比较困难，这就需要辅以更高分辨率的影像数据。特别地，中国始于“十一五”的高分辨率对地观测***专项的实施，使得高空间分辨率影像数据(<10m)获取技术进一步得到发展。然而，由于高空间分辨率影像光谱信息更加丰富， 异物同谱现象严重，很难用基于像素的分类方法提取耕地。另外，耕地不同物候的差异，例如，在无人机高分影像中，旱地与水田的色彩相差很大，旱地为黄色，水田为绿色，很难将两者同时提取。因此，高分影像中的耕地专题信息的提取难度更大。
针对高分影像的耕地专题信息提取，已有不少学者做了相关研究工作，Liu等2008年在“GIS与建筑环境联合会议”撰文“High qual ity prime farmland extract ionpattern based on object-oriented image analys is”，利用航空影像数据(3m),通过边缘检测分割影像，随后计算对象特征值，最后利用C4.5构建耕地信息提取决策树，实现高质量标准农田的探测。Lu等2007年在《Internat ional Journal of Remote Sens ing》撰文“Comparison between several feature extract ion/classification methods formapping compl icated agricultural land use patches us ing airbornehyperspectral data”，利用航空高光谱数据(2m)提取复杂农业用地覆盖区域。Duro等2012年在《Remote Sens ing of Environment》撰文“A comparison of pixel-based andobject-based image analys is with selected machine learning algorithms forthe class ificat ion of agricultural landscapes using SPOT-5HRG imagery”，使用SPOT影像，采用多尺度分割技术，分别利用随机森、支持向量机、决策树等监督分类方法，实现耕地信息的提取。诸如此类的基于分割后的对象，利用监督方法提取耕地信息的方法还很多，但是高分影像提供更加丰富的光谱和纹理特征信息，以上方法多为黑箱操作，并没有充分利用耕地地块信息自身独有的特点，导致提取效果不好。
目前，专门针对高分影像开展的非监督耕地信息提取方法很少，自动化程度并不高，对分割后的地理对象的语义信息利用并不充分。例如，Sun和Xu 2009 年在“Transactions of the CSAE”撰文“Comer extraction algorithm for high-resolut ion imageryof agricultural land”，利用针对Quickbird全色波段影像，先对影像数据进行分割，结合角特征点所提供的形状信息(耕地外形规则，一般具有4个或4个以上的角特征点)提取分割后的耕地和坑塘等形状规则的农用地信息。该研究为高分辨率全色影像耕地信息提取开辟了一个新的思路，一定程度上利用高分辨率影像数据中耕地地块的几何和纹理特征，以及耕地地块的独特性质进行算法设计，从而快速高效的提取耕地信息。但是其仅利用了单一的形状特征，对复杂环境的提取效果并不好。考虑到耕地信息具有明显的集中连片特征，且由于高分影像分割后的耕地对象较大，而对于集中居民地区域，由于内部地物丰富，包括房屋、林地、道路、水泥地等，形成形态各异的分割对象，相对耕地区域，该分割区域看上去比较破碎。基于这一特征，我们考虑聚类的方法将破碎居民地、林地等地物剔除，从而自动提取耕地信息。
发明内容
本发明要解决的技术问题是：克服现有技术的上述不足，提供一种基于三角网的面向对象耕地信息自动提取方法，其结合了面向对象影像分类、三角网分析和AUTOCLUST聚类算法，基于高分遥感影像，自动、高效、准确地进行耕地信息的提取。
为了解决以上技术问题，本发明提供的一种基于三角网的面向对象耕地信息自动提取方法，包括以下步骤：
步骤1、利用多尺度分割法对高空间分辨率遥感影像进行分割；
步骤2、剔除长条状分割对象；
步骤3、提取剩余分割对象的中心点；
步骤4、利用中心点构建Delaunay三角网；
步骤5、对步骤4获得的Delaunay三角网进行剥皮操作，剥离Delaunay三角网的外三角；
步骤6、对步骤5中剥皮后的Delaunay三角网进行聚类，得到三角网节点分布密集和分布稀疏的两种类型的若干个聚类结果；
步骤7、利用Vornoio图约束对聚类结果进行修正，修正方法如下：利用Delaunay三角网节点构建Voronoi图，计算Voronoi图中Voronoi多边形面积A，Voronoi多边形面积小于下限阈值时，初步确定该Voronoi多边形对应的中心点属于破碎区域，Voronoi多边形面积大于上限阈值时，初步确定该多边形对应的中心点为耕地区域；将判定为破碎区域的中心点删去或归入节点分布密集的聚类中，将判定为耕地区域的中心点归入节点分布稀疏的聚类中；
步骤8、删除节点分布密集的聚类所对应的中心点；
步骤9、若所述步骤1分割后得到分割对象的最大光谱方差Max.diff大于预设阈值，则判定该分割对象为林地，剔除该分割对象所对应的中心点，最终保留下来的所有中心点所对应的分割对象为耕地，得到耕地提取结果。
本发明还具有以下进一步的特征：
1、步骤1中，使用影像分析软件对高空间分辨率影像进行反复的分割尝试，并观察分割对象的分布特征，选择对耕地地块分割较好的分割尺度进行影像分割。
2、矩形拟合度特征数值小于0.5的分割对象判定为长条状分割对象，并剔除。
3、步骤3中，计算剩余分割对象的多边形面积，并根据面积计算中心点的坐标，所述中心点为分割对象多边形的质心。
4、步骤5中，选择距离参数d，删除Delaunay三角网中三角形边长大于d的边，d的取值范围为[100m,150m]。
5、在步骤7中，下限阈值为上限阈值为其中，Voronio_mean为生成的泰森多边形的面积的平均值，和为系数，且
6、
7、在步骤7或步骤8执行之前，进行如下判断，如果第j个聚类满足则将该聚类归入节点分布稀疏的聚类，其中表示第j个聚类的三角形面积的均值，表示所有聚类的三角形集合的面积均值的平均值，表示第j个聚类的三角形面积集合的方差。
8、预设阈值为2倍所有分割对象的最大光谱方差均值与所有分割对象的最大光谱方差标准差之和。
9、步骤6中采用AUTOCLUST方法进行聚类。
本发明基于Delaunay三角网结合AUTOCLUST聚类算法和多种优化约束方法实现对耕地信息的自动提取。本发明能够在面向对象影像分类的基础上将三角网和聚类算法结合，自动准确地从中选取可靠的耕地信息，同时，简化了耕地信息提取过程中人机交互强度，保证了耕地提取的整体性。
本发明的步骤1中使用影像分割技术，利用高空分辨率影像，通过反复的分割尝试并观察分割对象的分布特征，发现耕地信息具有明显的集中连片特征，分割后的耕地块都比较大，且由于低空无人机影像的高分辨率特点，对于集中居民地区域，其内部地物丰富，包括房屋、林地、道路、水泥地等，分割后形成形态各异的分割对象，相对耕地区域，该分割区域看上去比较破碎。另一方面，对于中国广大农村，居民地区域呈棋盘状分布在耕地中，一般地，集中居 住居民地和道路区域是耕地信息提取的主要干扰因素，而道路和沟渠等对象具有明显的长条状特征，与耕地信息不符，因此在步骤2中直接剔除此类对象。
进一步地，对于非长条形的地物，本发明基于三角网分析对非耕地的干扰因素进行逐步剔除，在这个过程中将引入AUTOCLUST聚类、Voronio图约束和方差约束，多种聚类约束算法一步步提高耕地提取的精度。
可见，本发明以三角网为基础，使用AUTOCLUST聚类为主要提取算法，利用长短边聚类实现聚类信息的传递，具体地，采用局部和全局方差表征长短边信息，删除三角网中的短边集合和长边集合并进一步扩大范围。
相比现有技术，本发明充分结合了耕地遥感影像分割后的特征信息以及三角网分析，形成了一种自动、高效、准确的耕地信息提取方法。具体创新点及有益效果如下：
第一，本发明能够充分利用Delaunay三角网的能自动适应不同区域点的密度变化的特征，结合基于图论的AUTOCLUST算法，利用长短边聚类实现聚类信息的传递，具体地，采用局部和全局方差表征长短边信息，从而能够使用该方法实现三角网的聚类。
第二，本发明在AUTOCLUST算法的基础上，使用泰森多边形约束进一步处理被忽略三角网中三角形面积的特征，根据Delaunay三角网得到其对偶图Voronoi图，进一步通过计算泰森多边形面积的方法，表征点的分布密度，解决使用AUTOCLUST聚类算法时出现的欠聚类的现象；使用聚类三角网面积均值约束，解决AUTOCLUST聚类中较大三角网的误聚类问题；使用最大方差约束，进一步优化算法，剔除平原地区零星林地。
综上，本发明充分挖掘高分遥感影像的语义信息，提出了一种以三角网为基础，结合AUTOCLUST聚类的耕地提取方法，自动、高效、准确地进行耕地信 息的提取。整个方法，自动化程度较高、精度较高，稳定性较高。不仅提高了高分遥感影像的耕地信息提取效率，同时，保证了耕地提取的整体性。
附图说明
下面结合附图对本发明作进一步的说明。
图1为本发明方法流程图。
图2a为无人机试验原始影像。
图2b为无人机试验原始影像分割尺度为50的分割结果。
图2c为无人机试验原始影像分割尺度为100的分割结果。
图2d为无人机试验原始影像分割尺度为180的分割结果。
图3a为Delaunay三角网示意图。
图3b为剥离Delaunay三角网的外三角示例图。
图4a为图2d分割结果构建的三角网。
图4b为AUTOCLUST聚类结果。
图4c为Voronio图下限阈值约束的聚类图。
图4d为Voronio图上限阈值约束的聚类图。
图5a为使用提出的方法的耕地提取结果。
图5b为使用图像分析软件eCognition进行耕地提取得到的结果。
具体实施方式
下面根据附图详细说明本发明，本发明的目的和效果将变得更加明显。
如图1所示为本发明实施例基于三角网的面向对象耕地信息自动提取方法的流程图，本实施方法包括以下步骤：
步骤1、利用多尺度分割法对高空间分辨率影像进行分割。本实施例中，使用eCognit ion对高空间分辨率影像进行反复的分割尝试，并观察分割对象的分 布特征，选择对耕地地块分割较好的分割尺度进行影像分割。本发明建议分割尺度的取值范围为[150,200]，以便突出高分影像中的居民地与耕地对应的分割对象的差异。
图2为三景无人机试验影像。本步骤中，利用影像分析软件eCognition进行多尺度分割，分割参数Shape为0.4，Smoothness为0.5。分别实施三个尺度的分割(包括尺度50、100、180)，通过观察不同尺度的分割结果，选择对耕地地块分割较好的尺度进行后续分析。容易发现，随着分割尺度的增加，分割对象的面积逐渐增大。相对分割尺度50和100，当分割尺度在180时，能够更好将分辨率0.2米的无人机影像数据中具有相同像素特征的耕地地块分割出来。
图2a为原始影像，图2b是分割尺度为50的分割结果，图2c是分割尺度为100的分割结果，图2d是分割尺度为180的分割结果。
步骤2、剔除长条形状的分割对象(道路，沟渠等)。本实施例中，将矩形拟合度特征(Rectangular fit)数值小于0.5的分割对象剔除，矩形拟合度特征(Rectangular fit)可用于剔除道路、沟渠等与耕地地块特征明显差异的长条状分割对象。
成功对无人机影像数据进行分割后，将使用三角网聚类剔除破碎居民地区域，以提取耕地信息。首先利用分割对象提取中心点(步骤3)，提取的中心点个数和该区域分割对象数目相同。图2d耕地虽然分布较集中，但是难免夹杂有道路、沟渠、田埂等地物，在高分辨率的影像中，这些地物都能被分割出来，且它们都属于长条形状地物，与耕地的矩形有较大差别。在步骤3提取的中心点中，多了这些地物的中心点，它们会影响方法精度，为提高方法可靠性，在构建三角网之前，需要将这些条状地物剔除。所以，首先使用空间特征矩形拟合度(Rectangular fit)阈值剔除条状地物(Rectangular fit<0.5)对应的中心点。
步骤3、提取剩余分割对象的中心点。本例中，首先计算剩余分割对象的多边形面积，并根据面积计算其中心点的坐标。
为计算分割对象多边形的中心点(重心或质心)，首先计算多边形的面积，如果某个分割对象的分割线有N个节点(xi,yi)组成，其中0≤i≤N-1。假设最后一个节点与第一个节点相同，那么多边形封闭。那么，多边形的面积计算公式表示为：
其中xN＝x0，则表示多边形封闭。在计算得到面积以后，多边形中心点(这里实际上表示重心或质心)的坐标(Cx，Cy)能够通过下式计算：
步骤4、利用步骤3提取的中心点构建Delaunay三角网,Delaunay三角网如图3a所示。
Delaunay三角网是相互不重叠且相互邻接的三角形集合，一般需满足两个准则：Delaunay三角网中的各三角形的外接圆中不存在其它点。在所有可能的三角网中，Delaunay三角网中的各三角形的最小角一般最大。根据以上特征，许多研究提出不同的三角网构建方法，流行的方法包括：分治法、扫描线算法、增量算法、快速增量构造算法、凸包算法。本发明重点不在于三角网构建算法，本实施例利用成熟的分治算法(Divide-and-Conquer Algorithm)构建Delaunay三角网，首先，尽量划分点集到足够小的子集，对不同子集构建子三角网，随后整合各子三角网，最后，利用局部优化法保证构建的三角网为Delaunay三角网，这样构建三角网的优点是时间效率较高，缺点是需要的递归运算量大，对内存空间的要求较高。
步骤5、对步骤4获得的Delaunay三角网进行剥皮操作，剥离Delaunay三角网的外三角。本例中，利用步骤4中获得的Delaunay三角网，选择距离参数d，距离参数d应当不小于三角网内部最长边，本发明建议d的取值范围为[100m,150m]，删除三角形边长大于d的边，并将仅与一个三角形关联的边作为点群的边界。如图3b所示为剥离Delaunay三角网的外三角的结果。
点群分布范围属于不确定性的问题。在GIS的研究中，常用作替代点群分布范围。但是，当点群的分布凸凹有致，那么可能在凸壳中包含很多不存在点的凹部区域，这样形成的凸壳并不能真实的代表点群的分布范围。一般地，基于Delaunay三角网的聚类步骤为：①利用点群构建三角网，如图3a；②选择一个距离参数d，通过参数d实现对三角网的“剥皮”，删除三角形边长大于d的边；③将与一个三角形关联的边提取作为点群的边界。图3b显示，通过“剥皮”操作，删除了边长很大的三角形，使得“剥皮”后的三角网更接近真实的点群分布。
步骤6、对步骤5中剥皮后的Delaunay三角网进行AUTOCLUST聚类，得到三角网节点分布密集和分布稀疏的两种类型的若干个聚类结果。AUTOCLUST聚类算法为现有成熟算法，可参阅Estivill-Castro和Lee于2002年发表在《Computers,Environment and urbansystems》上的文献“Argument free clustering for large spat ial point data setsvia boundary extract ion from Delaunay Diagram”。如图4a中A、B、C、D、E、F四块区域为高密度聚类所在位置。
除了使用AUTOCLUST聚类方法以外，也可以使用K-means、基于密度的空间聚类等方法对中心点进行聚类，获得分布密集和分布稀疏的两种类型的若干个 聚类。
步骤7、利用Vornoio图约束对AUTOCLUST聚类结果进行修正，避免过度聚类和欠聚类，如图4b所示，图中黑色四边形和三角形内代表未被聚类的居民点区域，即欠聚类。本实施例中，利用Delaunay三角网节点(也就是第4步中获得的中心点)构建Voronoi图，计算Voronoi图中Voronoi多边形面积A，Voronoi多边形面积大于上限阈值时，初步确定该多边形对应的中心点为耕地，Voronoi多边形面积小于下限阈值时，初步确定该Voronoi多边形对应的中心点属于破碎区域(非耕地)，将破碎区域的中心点删去或归入节点分布密集的聚类中。完成对聚类结果的修正，避免过度聚类和欠聚类现象。
如图4c所示，为使用下限阈值进行约束的结果，图中线条较粗的Voronoi多边形为面积小于下限值的多边形，初步确定为非耕地。如图4d所示，为实用上限阈值进行约束的结果，图中线条较粗的Voronoi多边形为面积大于下限值的多边形，初步确定为耕地。
基于Delaunay三角网构建Voronoi图，计算多边形面积，表示为A。所以，Voronoi多边形的点密度可以通过1/A表征。一般地，该值越小，点群密度越小，该值越大，表示点群密度越大。泰森多边形面积约束策略：一般地，AUTOCLUST进行聚类时，有不少区域出现欠聚类的现象，那么在采用泰森多边形进行约束时，使用上限阈值与下限阈值综合约束，一般地，当泰森多边形面积大于上限阈值时，该多边形对应的中心点一定是耕地，当泰森多边形面积小于下限阈值时，该多边形对应的中心点一定是属于破碎区域。这里的上下限阈值尽量设置松散。用以下公式过滤高密度点区域和低密度点区域：
这里Voronio_mean为生成的泰森多边形的面积的平均值，和为系数，且明显地，这里的系数决定了点的聚类程度，通过实验，为了避免过度聚类和欠聚类，设置作为一般值，对AUTOCLUST聚类结果进行修正。
在AUTOCLUST聚类的过程中，一些聚类主要还包含较大的三角网，且其面积较周围三角网面积还大，这可能是因为极少数的两三个耕地分割对象对应的中心点距离很近，而距离相邻的其它点较远，这种点往往不是居民区域分割对象的中心点。这种聚类中的三角网个数往往少，且三角形面积的平均值比需要识别的破碎区域的三角网面积的平均值大很多。
所以进一步引入聚类三角网面积均值的概念来去除AUTOCLUST聚类过程中这类错误聚类，方程如下：
其中表示第j个聚类的三角形面积的均值，Ti为第i个三角形的面积，nj为第j个聚类的三角形个数，表示所有聚类的三角形面积均值的平均值，表示第j个聚类的三角形面积集合的方差，m表示聚类的个数，i表示第j个聚类中的三角形序号。
如果第j个聚类满足则将该聚类归入分布稀疏的聚类中去。
步骤8、删除节点分布密集的聚类所对应的中心点，至此完成了非耕地区域的中心点的剔除。
步骤9、利用最大方差约束，剔除剩余的零星林地，若所述步骤1分割后得 到分割对象的最大光谱方差Max.diff大于预设阈值，则判定该分割对象为林地，剔除该分割对象所对应的中心点，最终保留下来的所有中心点所对应的分割对象为耕地，得到耕地提取结果并进行精度评估。本实施例中，所述预设阈值为2倍最大光谱方差均值与最大光谱方差标准差之和。
本实施例提出的方法虽然能够聚类居民地聚集区域，并将其剔除，但是对于平原地区，还零星存在不少林地，由于林地的纹理特征相差不大，对应的分割对象往往较大，很难用三角网聚类方法将其剔除。而林地的纹理粗糙，方差较大，耕地较为光滑，方差较小，所以这里提出方差约束的方法用以对耕地信息提取结果做修剪。本实施例采用2倍最大光谱方差均值与最大光谱方差标准差之和作为阈值进行识别。如下式：
Max.diff>2×mean_Max.diff+Standard Deviation_Max.diff (9)
其中mean_Max.diff表示分割图层中所有对象的最大光谱方差特征的均值；特征Max.diff的标准差表示为Standard Deviation_Max.diff。
图5a为采用基于三角网的面向对象耕地信息自动提取方法提取的耕地信息，图5a中，白色表示提取的耕地，黑色表示识别的非耕地，图5b表示利用eCognit ion软件进行监督分类提取的耕地信息，图中白色为耕地，黑色为非耕地。比较发现：本发明提取方法在保证精度要求的前提下，不仅能够在整体上保持耕地的集中连片特性，还能在一定程度上避免由于播种类型差异(物候差异)，导致的零星错分现象。同时，避免监督分类选择训练样本的过程，提高了自动化处理程度。
除上述实施例外，本发明还可以有其他实施方式。凡采用分割对象中心点聚类等形式的耕地信息提取技术方案，均落在本发明要求的保护范围。
Claims (10)
1.一种基于三角网的面向对象耕地信息自动提取方法，包括以下步骤：
步骤1、利用多尺度分割法对高空间分辨率影像进行分割；
步骤2、剔除长条状分割对象；
步骤3、提取剩余分割对象的中心点；
步骤4、利用中心点构建Delaunay三角网；
步骤5、对步骤4获得的Delaunay三角网进行剥皮操作，剥离Delaunay三角网的外三角；
步骤6、对步骤5中剥皮后的Delaunay三角网进行聚类，得到三角网节点分布密集和分布稀疏的两种类型的若干个聚类结果；
步骤7、利用Vornoio图约束对聚类结果进行修正，修正方法如下：利用Delaunay三角网节点构建Voronoi图，计算Voronoi图中Voronoi多边形面积，Voronoi多边形面积小于下限阈值时，初步确定该Voronoi多边形对应的中心点属于破碎区域，Voronoi多边形面积大于上限阈值时，初步确定该多边形对应的中心点为耕地区域；将判定为破碎区域的中心点删去或归入节点分布密集的聚类中，将判定为耕地区域的中心点归入节点分布稀疏的聚类中；
步骤8、删除节点分布密集的聚类所对应的中心点；
步骤9、若所述步骤1分割后得到分割对象的最大光谱方差Max.diff大于预设阈值，则判定该分割对象为林地，剔除该分割对象所对应的中心点，最终保留下来的所有中心点所对应的分割对象为耕地，得到耕地提取结果。
2.根据权利要求1所述的基于三角网的面向对象耕地信息自动提取方法，其特征在于：步骤1中，使用影像分析软件对高空间分辨率影像进行反复的分割尝试，并观察分割对象的分布特征，选择对耕地地块分割较好的分割尺度进行影像分割。
3.根据权利要求1所述的基于三角网的面向对象耕地信息自动提取方法，其特征在于：矩形拟合度特征数值小于0.5的分割对象判定为长条状分割对象，并剔除。
4.根据权利要求1所述的基于三角网的面向对象耕地信息自动提取方法，其特征在于：步骤3中，计算剩余分割对象的多边形面积，并根据面积计算中心点的坐标，所述中心点为分割对象多边形的质心。
5.根据权利要求1所述的基于三角网的面向对象耕地信息自动提取方法，其特征在于：步骤5中，选择距离参数d，删除Delaunay三角网中三角形边长大于d的边，d的取值范围为[100m,150m]。
6.根据权利要求1所述的基于三角网的面向对象耕地信息自动提取方法，其特征在于：在步骤7中，下限阈值为上限阈值为其中，Voronio_mean为生成的泰森多边形的面积的平均值，和为系数，且
7.根据权利要求6所述的基于三角网的面向对象耕地信息自动提取方法，其特征在于：
8.根据权利要求1所述的基于三角网的面向对象耕地信息自动提取方法，其特征在于：在步骤7或步骤8执行之前，进行如下判断，如果第j个聚类满足则将该聚类归入节点分布稀疏的聚类，其中表示第j个聚类的三角形面积的均值，表示所有聚类的三角形集合的面积均值的平均值，表示第j个聚类的三角形面积集合的方差。
9.根据权利要求1所述的基于三角网的面向对象耕地信息自动提取方法，其特征在于：预设阈值为2倍所有分割对象的最大光谱方差均值与所有分割对象的最大光谱方差标准差之和。
10.根据权利要求1所述的基于三角网的面向对象耕地信息自动提取方法，其特征在于：步骤6中采用AUTOCLUST方法进行自动聚类。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610934973.1A CN106548141B (zh) | 2016-11-01 | 2016-11-01 | 一种基于三角网的面向对象耕地信息自动提取方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610934973.1A CN106548141B (zh) | 2016-11-01 | 2016-11-01 | 一种基于三角网的面向对象耕地信息自动提取方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106548141A true CN106548141A (zh) | 2017-03-29 |
| CN106548141B CN106548141B (zh) | 2019-07-26 |
Family
ID=58392223
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610934973.1A Active CN106548141B (zh) | 2016-11-01 | 2016-11-01 | 一种基于三角网的面向对象耕地信息自动提取方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106548141B (zh) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107221001A (zh) * | 2017-05-18 | 2017-09-29 | 苏州恩巨网络有限公司 | 一种区域信息获取方法及装置 |
| CN107564077A (zh) * | 2017-09-04 | 2018-01-09 | 济南浪潮高新科技投资发展有限公司 | 一种绘制道路网的方法及装置 |
| CN107705002A (zh) * | 2017-09-21 | 2018-02-16 | 中国矿业大学(北京) | 矿区土壤重金属含量采样点异常高值影响范围的确定方法 |
| CN107766496A (zh) * | 2017-10-19 | 2018-03-06 | 大连大学 | 基于MapReduce与倒排泰森多边形的大规模近邻查询*** |
| CN107766495A (zh) * | 2017-10-19 | 2018-03-06 | 大连大学 | 基于Voronoi图的分布式时空索引方法 |
| CN107958087A (zh) * | 2017-10-19 | 2018-04-24 | 大连大学 | 移动医疗的倒排泰森多边形的分布式时空索引的代表点获取方法 |
| CN108021690A (zh) * | 2017-10-19 | 2018-05-11 | 大连大学 | 倒排Voronoi索引 |
| CN108961403A (zh) * | 2018-06-22 | 2018-12-07 | 武汉大学 | 一种运用开源街道地图提取主干道的方法 |
| CN111768505A (zh) * | 2020-07-01 | 2020-10-13 | 江西师范大学 | 一种平面点集形状重建方法、装置及电子设备 |
| CN112330700A (zh) * | 2020-11-16 | 2021-02-05 | 四川航天神坤科技有限公司 | 一种卫星影像的耕地地块提取方法 |
| CN112598724A (zh) * | 2021-03-01 | 2021-04-02 | 武大吉奥信息技术有限公司 | 一种改进的基于tin的面状要素中心线提取方法 |
| CN112926029A (zh) * | 2021-01-29 | 2021-06-08 | 四川省环境政策研究与规划院 | 一种用于农村生活污水治理的聚居区识别划分方法 |
| CN114354464A (zh) * | 2021-12-22 | 2022-04-15 | 南京大学 | 一种对不同金属纳米颗粒高光谱库之间的相似性进行定量分析的方法 |
| CN114648491A (zh) * | 2022-02-22 | 2022-06-21 | 山东土地集团数字科技有限公司 | 一种用于耕地道路的遥感监测***和方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090315882A1 (en) * | 2008-06-19 | 2009-12-24 | Xerox Corporation | Reducing the number of vertices in a clustered delaunay mesh |
| CN102073981A (zh) * | 2011-02-23 | 2011-05-25 | 中国测绘科学研究院 | 一种关联要素限制下的点群地理实体选取方法 |
| CN104036024A (zh) * | 2014-06-27 | 2014-09-10 | 浙江大学 | 一种基于GACUC和Delaunay三角网的空间聚类方法 |
| JP2016062170A (ja) * | 2014-09-16 | 2016-04-25 | 三菱製紙株式会社 | 光透過性導電材料 |
| CN105608474A (zh) * | 2015-12-30 | 2016-05-25 | 二十一世纪空间技术应用股份有限公司 | 基于高分辨率影像的区域自适应耕地提取方法 |
| CN105719014A (zh) * | 2016-01-18 | 2016-06-29 | 中国农业大学 | 省域耕地后备资源开发区划方法及*** |
-
2016
- 2016-11-01 CN CN201610934973.1A patent/CN106548141B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090315882A1 (en) * | 2008-06-19 | 2009-12-24 | Xerox Corporation | Reducing the number of vertices in a clustered delaunay mesh |
| CN102073981A (zh) * | 2011-02-23 | 2011-05-25 | 中国测绘科学研究院 | 一种关联要素限制下的点群地理实体选取方法 |
| CN104036024A (zh) * | 2014-06-27 | 2014-09-10 | 浙江大学 | 一种基于GACUC和Delaunay三角网的空间聚类方法 |
| JP2016062170A (ja) * | 2014-09-16 | 2016-04-25 | 三菱製紙株式会社 | 光透過性導電材料 |
| CN105608474A (zh) * | 2015-12-30 | 2016-05-25 | 二十一世纪空间技术应用股份有限公司 | 基于高分辨率影像的区域自适应耕地提取方法 |
| CN105719014A (zh) * | 2016-01-18 | 2016-06-29 | 中国农业大学 | 省域耕地后备资源开发区划方法及*** |
Non-Patent Citations (3)
| Title |
|---|
| JINFEI LIU等: "Effective Privacy Preserved Clustering Based on Voronoi Diagram", 《2011 EIGHTH INTERNATIONAL SYMPOSIUM ON VORONOI DIAGRAMS IN SCIENCE AND ENGINEERING》 * |
| 孙小丹等: "农用地遥感影像信息的角提取方法", 《农业工程学报》 * |
| 彭仪普等: "Delaunay三角网与Voronoi图在GIS中的应用研究", 《测绘工程》 * |
Cited By (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107221001A (zh) * | 2017-05-18 | 2017-09-29 | 苏州恩巨网络有限公司 | 一种区域信息获取方法及装置 |
| CN107564077B (zh) * | 2017-09-04 | 2019-12-10 | 浪潮集团有限公司 | 一种绘制道路网的方法及装置 |
| CN107564077A (zh) * | 2017-09-04 | 2018-01-09 | 济南浪潮高新科技投资发展有限公司 | 一种绘制道路网的方法及装置 |
| CN107705002A (zh) * | 2017-09-21 | 2018-02-16 | 中国矿业大学(北京) | 矿区土壤重金属含量采样点异常高值影响范围的确定方法 |
| CN107705002B (zh) * | 2017-09-21 | 2020-07-03 | 中国矿业大学(北京) | 矿区土壤重金属含量采样点异常高值影响范围的确定方法 |
| CN107766495A (zh) * | 2017-10-19 | 2018-03-06 | 大连大学 | 基于Voronoi图的分布式时空索引方法 |
| CN107958087A (zh) * | 2017-10-19 | 2018-04-24 | 大连大学 | 移动医疗的倒排泰森多边形的分布式时空索引的代表点获取方法 |
| CN107766496A (zh) * | 2017-10-19 | 2018-03-06 | 大连大学 | 基于MapReduce与倒排泰森多边形的大规模近邻查询*** |
| CN108021690A (zh) * | 2017-10-19 | 2018-05-11 | 大连大学 | 倒排Voronoi索引 |
| CN108961403A (zh) * | 2018-06-22 | 2018-12-07 | 武汉大学 | 一种运用开源街道地图提取主干道的方法 |
| CN111768505A (zh) * | 2020-07-01 | 2020-10-13 | 江西师范大学 | 一种平面点集形状重建方法、装置及电子设备 |
| CN111768505B (zh) * | 2020-07-01 | 2023-10-17 | 江西师范大学 | 一种平面点集形状重建方法、装置及电子设备 |
| CN112330700A (zh) * | 2020-11-16 | 2021-02-05 | 四川航天神坤科技有限公司 | 一种卫星影像的耕地地块提取方法 |
| CN112926029B (zh) * | 2021-01-29 | 2022-08-26 | 四川省环境政策研究与规划院 | 一种用于农村生活污水治理的聚居区识别划分方法 |
| CN112926029A (zh) * | 2021-01-29 | 2021-06-08 | 四川省环境政策研究与规划院 | 一种用于农村生活污水治理的聚居区识别划分方法 |
| CN112598724A (zh) * | 2021-03-01 | 2021-04-02 | 武大吉奥信息技术有限公司 | 一种改进的基于tin的面状要素中心线提取方法 |
| CN114354464A (zh) * | 2021-12-22 | 2022-04-15 | 南京大学 | 一种对不同金属纳米颗粒高光谱库之间的相似性进行定量分析的方法 |
| CN114354464B (zh) * | 2021-12-22 | 2024-04-09 | 南京大学 | 一种对不同金属纳米颗粒高光谱库之间的相似性进行定量分析的方法 |
| CN114648491A (zh) * | 2022-02-22 | 2022-06-21 | 山东土地集团数字科技有限公司 | 一种用于耕地道路的遥感监测***和方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106548141B (zh) | 2019-07-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106548141B (zh) | 一种基于三角网的面向对象耕地信息自动提取方法 | |
| CN101840581B (zh) | 一种从卫星遥感影像中提取建筑物轮廓的方法 | |
| Blaschke et al. | Image segmentation methods for object-based analysis and classification | |
| CN103839267B (zh) | 一种基于形态学建筑物指数的建筑物提取方法 | |
| CN104156964B (zh) | 一种综合mrf和贝叶斯网络的遥感影像区域分割方法 | |
| CN109919944A (zh) | 一种复杂场景建筑物变化检测的联合超像素图割优化方法 | |
| CN108985247A (zh) | 多光谱图像城市道路识别方法 | |
| CN106709989B (zh) | 基于面向对象影像特征的地理本体建模与语义推理方法 | |
| CN102096816A (zh) | 基于最小生成树的多尺度多层次影像分割方法 | |
| WO2024020744A1 (zh) | 基于高分辨率影像的生态图斑提取方法 | |
| Hu et al. | Integrating CART algorithm and multi-source remote sensing data to estimate sub-pixel impervious surface coverage: a case study from Beijing Municipality, China | |
| CN109766824A (zh) | 基于模糊证据理论的主被动遥感数据融合分类方法 | |
| CN107194857B (zh) | 基于高分卫星遥感数据的可利用土地资源分析*** | |
| CN105701856B (zh) | 一种植被提取方法及*** | |
| Guo et al. | Dual-concentrated network with morphological features for tree species classification using hyperspectral image | |
| Nelson et al. | Spatial statistical techniques for aggregating point objects extracted from high spatial resolution remotely sensed imagery | |
| CN116415318B (zh) | 一种基于数学形态学的内流区湖泊水文连通性建模方法 | |
| CN112380994A (zh) | 一种基于时序分类与空间分析的人工林空间模式识别方法 | |
| CN103778413B (zh) | 一种遥感影像欠分割对象自动识别方法 | |
| Huck et al. | Delineating parameters for object-based urban structure mapping in Santiago de Chile using QuickBird data | |
| Zhang et al. | Urban road extraction from high-resolution remote sensing images based on semantic model | |
| CN112907567B (zh) | 基于空间推理方法的sar图像有序人造构筑物提取方法 | |
| Wang et al. | Combining SAR images with land cover products for rapid urban flood mapping | |
| Forghani et al. | Extracting terrain categories from multi-source satellite imagery | |
| CN112487977A (zh) | 一种基于国产gf2遥感影像的机场跑道提取方法 |
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 |