CN115409742A - 一种基于园林绿化的植被覆盖密度评估方法 - Google Patents

一种基于园林绿化的植被覆盖密度评估方法 Download PDF

Info

Publication number
CN115409742A
CN115409742A CN202211358958.9A CN202211358958A CN115409742A CN 115409742 A CN115409742 A CN 115409742A CN 202211358958 A CN202211358958 A CN 202211358958A CN 115409742 A CN115409742 A CN 115409742A
Authority
CN
China
Prior art keywords
noise
image
area
window size
suspicious
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
Application number
CN202211358958.9A
Other languages
English (en)
Other versions
CN115409742B (zh
Inventor
袁传镇
张海燕
于丽霞
张宏智
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jinxiang County Forestry Protection And Development Service Center Jinxiang County Wetland Protection Center Jinxiang County Wildlife Protection Center Jinxiang County State Owned Baiwa Forest Farm
Original Assignee
Jinxiang County Forestry Protection And Development Service Center Jinxiang County Wetland Protection Center Jinxiang County Wildlife Protection Center Jinxiang County State Owned Baiwa Forest Farm
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jinxiang County Forestry Protection And Development Service Center Jinxiang County Wetland Protection Center Jinxiang County Wildlife Protection Center Jinxiang County State Owned Baiwa Forest Farm filed Critical Jinxiang County Forestry Protection And Development Service Center Jinxiang County Wetland Protection Center Jinxiang County Wildlife Protection Center Jinxiang County State Owned Baiwa Forest Farm
Priority to CN202211358958.9A priority Critical patent/CN115409742B/zh
Publication of CN115409742A publication Critical patent/CN115409742A/zh
Application granted granted Critical
Publication of CN115409742B publication Critical patent/CN115409742B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/62Analysis of geometric attributes of area, perimeter, diameter or volume
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/764Arrangements for image or video recognition or understanding using pattern recognition or machine learning using classification, e.g. of video objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10032Satellite or aerial image; Remote sensing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30181Earth observation
    • G06T2207/30188Vegetation; Agriculture

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Computing Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Geometry (AREA)
  • Databases & Information Systems (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Multimedia (AREA)
  • Image Processing (AREA)

Abstract

本发明涉及图像处理技术领域,具体涉及一种基于园林绿化的植被覆盖密度评估方法,该方法包括:获取园林高空遥感图像,将初步降噪后的图像中的背景部分进行掩膜遮盖,进而获取图像中所有闭合的连通域的面积,根据面积对所有闭合的连通域进行分类;根据连通域的数量和面积得到类别特征值,进而计算评价特征值,确定可疑噪声区域和非可疑噪声区域;根据可疑噪声区域内像素点到非可疑噪声区域边缘的最短距离得到像素点的特征长度,确定最大窗口尺寸;根据最大窗口尺寸确定最佳窗口尺寸;利用最佳窗口尺寸对初始遥感图像进行二次降噪,根据降噪后的图像获得植被的分布密度。本发明能够获得更加准确的植被分布密度的评估结果。

Description

一种基于园林绿化的植被覆盖密度评估方法
技术领域
本发明涉及图像处理技术领域,具体涉及一种基于园林绿化的植被覆盖密度评估方法。
背景技术
城市园林的绿植及其绿植覆盖率估算是对于城市生态以及居住环境的重要参考。早期利用常规技术进行树木株数及其树冠覆盖率的测算比较费工、费时。现有利用无人机遥感技术对城市园林绿化植被进行密度评估,这对于园林植被的养护、布局、生长监测等相关领域都具有重要意义。
无人机在高空拍摄园林遥感图像时,总会受到风力、光线等因素干扰,进而造成图像质量参差不齐,所以需要先进行一系列的预处理操作提高图像的可信度和清晰度,为后期的特征提取、分类等提供了更精确可靠的图像信息,但常规的滤波、均衡化等预处理操作存在丢失图像细节的问题,虽然对整体绿化覆盖率的评估结果影响不大,但对于不同种类绿植的覆盖密度进行评估时,图像细节丢失也会带来不同种类绿植之间层次感的丢失,对评估结果影响较大,使得评估结果并不准确。
发明内容
为了解决上述技术问题,本发明的目的在于提供一种基于园林绿化的植被覆盖密度评估方法,所采用的技术方案具体如下:
获取园林高空遥感图像,对该图像进行初步降噪,将初步降噪后的图像中的背景部分进行掩膜遮盖得到园林高空的初始遥感图像;获取初始遥感图像中所有闭合的连通域的面积,根据面积对所有闭合的连通域进行分类;
根据分类结果中各类别内闭合的连通域的数量和面积得到类别特征值,根据最大类别特征值对应的类别内连通域的面积得到各闭合的连通域的评价特征值,根据所述评价特征值确定可疑噪声区域和非可疑噪声区域;
根据可疑噪声区域内像素点到非可疑噪声区域边缘的最短距离得到可疑噪声区域内像素点对应的特征长度,根据特征长度确定像素点对应的最大窗口尺寸;在小于或等于最大窗口尺寸的范围内选择不同的尺寸,根据不同的尺寸对应的窗口内属于可疑噪声区域的像素点的数量确定最佳窗口尺寸;
利用最佳窗口尺寸对初始遥感图像进行二次降噪,根据降噪后的图像获得植被的分布密度。
优选地,所述将初步降噪后的图像中的背景部分进行掩膜遮盖得到园林高空的初始遥感图像具体为:
对初步降噪后的图像进行边缘检测,获取各闭合边缘的边缘像素点,并获取各边缘像素点的梯度方向,根据边缘像素点的梯度方向上相邻的边缘像素点之间的像素值差值得到边缘像素点的像素差异值,根据闭合边缘的边缘像素点对应的像素差异值确定背景区域;对背景区域进行掩膜遮盖得到园林高空的初始遥感图像。
优选地,所述类别特征值的获取方法具体为:
Figure DEST_PATH_IMAGE001
其中,
Figure 823741DEST_PATH_IMAGE002
表示类别
Figure DEST_PATH_IMAGE003
对应的类别特征值,
Figure 374064DEST_PATH_IMAGE004
表示类别
Figure 447063DEST_PATH_IMAGE003
内闭合的连通域的总数量,
Figure DEST_PATH_IMAGE005
表示类别
Figure 759226DEST_PATH_IMAGE003
内所有闭合的连通域的面积均值,
Figure 652096DEST_PATH_IMAGE006
表示所有类别中所有闭合的连通域的面积均值,
Figure DEST_PATH_IMAGE007
表示类别
Figure 737602DEST_PATH_IMAGE003
内所有闭合的连通域的面积的方差,exp()表示以自然常数e为底的指数函数。
优选地,所述评价特征值的获取方法具体为:
Figure 415839DEST_PATH_IMAGE008
其中,
Figure DEST_PATH_IMAGE009
表示第v个闭合的连通域的评价特征值,
Figure 931527DEST_PATH_IMAGE010
表示最大类别特征值对应的类别内所有闭合的连通域的面积均值,
Figure DEST_PATH_IMAGE011
表示第v个闭合的连通域的面积。
优选地,所述根据所述评价特征值确定可疑噪声区域和非可疑噪声区域具体为:
设置特征阈值,若闭合的连通域对应的评价特征值大于或等于特征阈值,则该闭合的连通域为可疑噪声区域,否则为非可疑噪声区域。
优选地,述根据特征长度确定该像素点对应的最大窗口尺寸具体为:
若可疑噪声区域内像素点的特征长度的取值为奇数,则该像素点对应的最大窗口尺寸与特征长度相同;若可疑噪声区域内像素点的特征长度的取值为偶数,则该像素点对应的最大窗口尺寸为特征长度的取值加1。
优选地,所述最佳窗口尺寸的获取方法具体为:
若可疑噪声区域内像素点对应的最大窗口尺寸的取值为3,则该像素点对应的最佳窗口尺寸的大小为3*3;若像素点可疑噪声区域内像素点对应的最大窗口尺寸的取值大于3,则在小于最大窗口尺寸的范围内选择不同的尺寸,计算不同的尺寸对应的窗口内属于可疑噪声区域内的像素点数量与窗口内像素点的总数量之间的比值,获取比值的最小值对应的尺寸记为像素点的最佳窗口尺寸。
本发明实施例至少具有如下有益效果:
本发明在对园林高空遥感图像进行初步降噪后,将图像中的背景部分进行掩膜遮盖得到初始遥感图像,考虑到了背景部分对图像进行二次滤波降噪的影响,能够在一定程度避免出现二次滤波降噪时出现图像失真的问题;根据分类结果中各类别内闭合的连通域的数量和面积得到类别特征值,通过分析各类别内闭合的连通域的数量特征以及面积特征,能够筛选出符合噪声部分特征的类别,进而可以获取噪声部分对应的连通域面积特征,能够进一步确认图像中可能属于噪声的区域以及不可能属于噪声的区域;进而根据可疑噪声区域内像素点到非可疑噪声区域边缘的最短距离得到可疑噪声区域内像素点对应的特征长度,充分考虑对像素点进行滤波降噪时,自适应选择的尺寸对应的窗口内应当尽可能多的包含非噪声部分,使得最终二次降噪时能够保留图像中细节纹理的同时,降噪效果更好,进而利用二次降噪后的图像获得植被的分布密度时,能够获得更加准确的评估结果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案和优点,下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它附图。
图1是本发明的一种基于园林绿化的植被覆盖密度评估方法的方法流程图。
具体实施方式
为了更进一步阐述本发明为达成预定发明目的所采取的技术手段及功效,以下结合附图及较佳实施例,对依据本发明提出的一种基于园林绿化的植被覆盖密度评估方法,其具体实施方式、结构、特征及其功效,详细说明如下。在下述说明中,不同的“一个实施例”或“另一个实施例”指的不一定是同一实施例。此外,一或多个实施例中的特定特征、结构或特点可由任何合适形式组合。
除非另有定义,本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。
下面结合附图具体的说明本发明所提供的一种基于园林绿化的植被覆盖密度评估方法的具体方案。
实施例:
本发明所针对的情景具体为:无人机在高空拍摄园林遥感图像时,总会受到风力、光线等因素干扰造成图像质量参差不齐,所以需要先进行一系列的预处理操作提高图像的可信度和清晰度,为后期的特征提取、分类等提供了更精确可靠的图像信息,但常规的滤波、均衡化等预处理操作存在丢失图像细节的问题,虽然对整体绿化覆盖率的评估结果影响不大,但对于不同种类绿植的覆盖密度进行评估时,图像细节丢失也会带来不同种类绿植之间层次感的丢失,对评估结果影响较大。
请参阅图1,其示出了本发明一个实施例提供的一种基于园林绿化的植被覆盖密度评估方法的方法流程图,该方法包括以下步骤:
步骤一,获取园林高空遥感图像,对该图像进行初步降噪,将初步降噪后的图像中的背景部分进行掩膜遮盖得到园林高空的初始遥感图像;获取初始遥感图像中所有闭合的连通域的面积,根据面积对所有闭合的连通域进行分类。
首先,利用无人机在高空拍摄园林的图像,并对该图像进行灰度化处理得到园林高空遥感图像,对图像进行灰度化处理是为了减少过多色彩信息的干扰,并且在进行预处理操作时相较于彩色图像具有更快的速度。
需要说明的是,在遥感数据采集过程中,由于大气、光线、反射光等干扰因素,经常会导致园林高空遥感图像内含有大量椒盐噪声,因此在进行绿植覆盖密度分析之前需要对园林高空遥感图像进行噪声去除的操作。一般情况下,噪声是比较小的图斑,类似于在一张面饼上随机分布的椒盐,所以可以将其看作椒盐点去除。
一般中值滤波对椒盐噪声点的去除效果最佳,但是不同尺度下的滤波核去噪效果明显不同,滤波尺度较小会导致噪声点去除不完全,尺度过大则又会丢失图像原本的细节信息。并且由于种植密度、光线角度的不同,园林遥感图像上不同区域的光照程度存在差异,在绿植较密集的区域,地面阴影连成一片,使其呈现较低的灰度值区域,这部分区域极容易与绿植部分发生混淆,在对图像进行降噪处理时,若滤波核尺寸选择不当,则极易导致原始图像上较弱的纹理消失。因此,需要选择合适尺寸的滤波核对图像进行滤波处理,使得图像中的噪声能够被去除的同时不会丢失图像原始的细节纹理。
然后,为了避免滤波核尺寸选择不当使得图像失真,首先对园林高空遥感图像进行初步降噪处理,在本实施例中采用中值滤波对园林高空遥感图像进行初步降噪处理,由于滤波核尺寸过小可能会导致图像中的噪声去除不完全,滤波核尺寸过大可能会丢失图像中的细节纹理,因此在对园林高空遥感图像进行初步降噪时,选择最小的滤波核尺寸,即尺寸大小为3*3。利用该尺寸大小对园林高空遥感图像进行降噪处理时,可以将边缘平滑的影响程度降到最低,但是噪声不会被完全去除,得到一幅降噪效果并不理想的遥感图像,此时图像上大量噪声可以去除,但仍有噪声残留。进而需要选择合适的滤波核尺寸对初步降噪后的遥感图像进行二次滤波处理。
进一步的,获取初步降噪后的图像中的背景部分,将初步降噪后的图像中的背景部分进行掩膜遮盖得到园林高空的初始遥感图像。对图像中的背景部分进行掩膜遮盖,使得图像进行二次滤波时仅对图像中的前景部分进行处理,可以仅消除噪声,并且不对背景部分进行滤波进而避免出现前景部分的边缘过度平滑的问题,进而在一定程度上避免二次滤波时图像失真的问题。其中,对图像中的背景区域进行掩膜遮盖为公知技术,在此不再过多介绍。
具体地,对于园林绿植的遥感图像而言,所有植物的图像为前景信息,地面及其他背景部分属于背景信息。提取图像中的背景部分为公知技术,实施者可根据实际情况选择合适的方法进行获取。
在本实施例中,获取图像中的背景部分的方法具体为:对初步降噪后的图像进行边缘检测,获取各闭合边缘的边缘像素点,利用sobel算子对初步降噪后的图像进行梯度检测,获得各边缘像素点的梯度方向,其中,sobel算子是利用两个卷积核分别对图像水平和垂直方向进行卷积,并结合两个卷积结果求出每个边缘像素点的梯度方向。同时,利用sobel算子对图像进行边缘检测为公知技术,在此不再过多介绍。
由于图像中的前景部分相较于背景部分之间的差异较大,则前景部分的区域的边缘像素点的像素值,相对于背景部分的区域内的像素点的像素值,都是统一高于或统一低于背景部分的区域内的像素点的像素值,即闭合边缘上各边缘像素点与其梯度方向上相邻的像素点的灰度值差异性保持一致。故根据边缘像素点的梯度方向上相邻的边缘像素点之间的像素值差值得到边缘像素点的像素差异值,根据闭合边缘的边缘像素点对应的像素差异值的正负性是否统一确定背景区域,若统一,则对应的闭合边缘内的区域为前景区域,其他闭合边缘内的区域为背景区域。
最后,由于本实施例对园林中的植被图像进行预处理操作,则图像中的前景部分可能包含了多种植被,同时存在噪声点,则每个植被在图像中均对应一个局部最大的闭合的连通域,在植被对应的局部最大的闭合的连通域内,可能存在较多的噪声点,且噪声点表现在图像中也为闭合的连通域,在本实施例中,以一个植被对应的局部最大的闭合的连通域为例进行说明,图像中所有植被对应的局部最大的闭合的连通域均作相同的处理。即首先获取初始遥感图像中的任意一个局部最大连通域。
具体地,所述局部最大连通域的获取方法具体为:判断闭合连通域的包含关系,当一个连通域无法被其余连通域所包含,则该连通域为初始遥感图像内部所包含的局部最大连通域。获取局部最大连通域中所有闭合的连通域的面积,根据面积对所有闭合的连通域进行分类,得到至少两个类别。其中,局部最大连通域内各闭合的连通域可能为噪声部分,也可能为植被的纹理部分。
步骤二,根据分类结果中各类别内闭合的连通域的数量和面积得到类别特征值,根据最大类别特征值对应的类别内连通域的面积得到各闭合的连通域的评价特征值,根据所述评价特征值确定可疑噪声区域和非可疑噪声区域。
首先,由于椒盐噪声具有随机分布、面积较小且均匀的特征,对各类别内闭合的连通域进行分析,若类别内闭合的连通域的面积越小、越均匀,且类别内闭合的连通域的数量相对较多,则该类别内各闭合的连通域越可能属于椒盐噪声。即可根据分类结果中各类别内闭合的连通域的数量和面积得到类别特征值,用公式表示为:
Figure 655901DEST_PATH_IMAGE012
其中,
Figure 217201DEST_PATH_IMAGE002
表示类别
Figure 264791DEST_PATH_IMAGE003
对应的类别特征值,
Figure 387599DEST_PATH_IMAGE004
表示类别
Figure 520640DEST_PATH_IMAGE003
内闭合的连通域的总数量,
Figure 888561DEST_PATH_IMAGE005
表示类别
Figure 790658DEST_PATH_IMAGE003
内所有闭合的连通域的面积均值,
Figure 84367DEST_PATH_IMAGE006
表示所有类别中所有闭合的连通域的面积均值,
Figure 439125DEST_PATH_IMAGE007
表示类别
Figure 607807DEST_PATH_IMAGE003
内所有闭合的连通域的面积的方差,exp()表示以自然常数e为底的指数函数。
Figure DEST_PATH_IMAGE013
表示类别
Figure 177460DEST_PATH_IMAGE014
内所有闭合的连通域的面积均值与所有类别中所有闭合的连通域的面积均值之间的比值,能够反应类别
Figure 156917DEST_PATH_IMAGE014
内所有闭合的连通域的面积均值与所有闭合的连通域的面积均值之间的差异关系,当类别内所有闭合的连通域的面积均偏小时,该比值的取值越小,对应的类别特征值的取值越大,说明该类别内各闭合的连通域越可能属于椒盐噪声,进而利用指数函数进行反比例归一化并矫正逻辑关系。
Figure DEST_PATH_IMAGE015
能够反应类别
Figure 832925DEST_PATH_IMAGE014
内所有闭合的连通域的面积的波动情况,该取值越小,说明类别
Figure 306763DEST_PATH_IMAGE014
内所有闭合的连通域的面积之间的波动越小,即差异越小,进而说明类别
Figure 917873DEST_PATH_IMAGE014
内所有闭合的连通域的面积越均匀,则对应的类别特征值的取值越大,说明该类别内闭合的连通域越符合椒盐噪声点分布具有普遍性、大小均匀的特征,该类别内闭合的连通域越可能属于椒盐噪声,而不同的前景图像的内部纹理是不具备普遍性,且面积可能各不相同,同时利用指数函数进行反比例归一化并矫正逻辑关系。
Figure 51920DEST_PATH_IMAGE016
表示对类别
Figure 646849DEST_PATH_IMAGE014
内所有闭合的连通域的面积的均值参数和方差参数进行L2范数的计算,该均值参数和方差参数均属于类别
Figure 658799DEST_PATH_IMAGE014
内闭合的连通域的面积的数值特征,因此利用L2范数进行整合。考虑到噪声点的分布数量比图像中纹理信息的分布数量较多的情况,当
Figure 858836DEST_PATH_IMAGE017
的取值越大时,说明类别内闭合的连通域的数量越多,则对应的类别特征值的取值越大,说明该类别越可能属于椒盐噪声。
由于椒盐噪声具有随机分布、面积较小且均匀,同时数量相对较多的特征,故当类别特征值的取值越大时,说明对应类别内闭合的连通域的面积越小且越均匀,闭合的连通域的数量越少,则该类别越可能为椒盐噪声所在类别,故选择所有类别中类别特征值最大的类别,该类别即为椒盐噪声所在类别,将该类别内所有闭合的连通域的面积均值作为判断噪声分布的可参考标准。
进一步的,将局部最大连通域内各闭合的连通域的面积分别与椒盐噪声所在类别内所有闭合的连通域的面积均值进行比较,当两者之间的差异越小时,则说明该闭合的连通域越可能属于椒盐噪声。基于此,根据最大类别特征值对应的类别内连通域的面积得到各闭合的连通域的评价特征值,用公式表示为:
Figure 432293DEST_PATH_IMAGE018
其中,
Figure 248939DEST_PATH_IMAGE009
表示第v个闭合的连通域的评价特征值,
Figure 799001DEST_PATH_IMAGE010
表示最大类别特征值对应的类别内所有闭合的连通域的面积均值,即噪声所在类别内所有闭合的连通域的面积均值,
Figure 650282DEST_PATH_IMAGE011
表示第v个闭合的连通域的面积。
Figure 391711DEST_PATH_IMAGE019
表示椒盐噪声所在类别内所有闭合的连通域的面积均值与第v个闭合的连通域的面积之间的比值,该比值能够反应两者之间的相似程度,该比值的取值越接近于1,则说明第v个闭合的连通域的面积越符合噪声的面积特征,第v个闭合的连通域属于噪声的可能性越大。
Figure 446385DEST_PATH_IMAGE020
表示数值1与两者之间相似程度的差异,该差异越小,说明第v个闭合的连通域属于噪声的可能性越大。
Figure 49405DEST_PATH_IMAGE021
表示第v个闭合的连通域的概率值,即所述差异越小,对应的闭合的连通域为噪声的概率值越大,即评价特征值的取值越大。
最后,按照上述方法获得所有闭合的连通域对应的评价特征值,设置特征阈值,若闭合的连通域对应的评价特征值大于或等于特征阈值,则该闭合的连通域为可疑噪声区域,否则为非可疑噪声区域。其中,非可疑噪声区域可能为图像中植被自身的纹理特征所在区域。
步骤三,根据可疑噪声区域内像素点到非可疑噪声区域边缘的最短距离得到可疑噪声区域内像素点对应的特征长度,根据特征长度确定该像素点对应的最大窗口尺寸;在小于或等于最大窗口尺寸的范围内选择不同的尺寸,根据不同的尺寸对应的窗口内属于可疑噪声区域的像素点的数量确定最佳窗口尺寸。
具体地,按照上述方法已经获得可疑噪声区域以及非可疑噪声区域,将非可疑噪声区域进行掩膜遮盖处理,即认为这部分区域内不存在噪声点,则不需要对其进行二次滤波处理,仅对可疑噪声区域中的各像素点进行处理即可。其中,对图像中的区域进行掩膜遮盖为公知技术,在此不再过多介绍。
在对可疑噪声区域中的可疑噪声部分进行去除时,需要根据局部最大连通域内包含的非可疑噪声区域中的纹理信息,判断将不同的可疑噪声区域去除时需要的滤波核尺寸。由于在对图像进行滤波时,滤波核尺寸对应的窗口内属于非噪声的像素点越多,去除噪声的效果越好,滤波核尺寸对应的窗口内各种元素越复杂,即滤波核尺寸对应的窗口内既含有较多非噪声的像素点,也含有较多噪声的像素点,则通过滤波进行去除噪声的效果越差。
基于此,在对可疑噪声区域内像素点进行滤波处理时,所采用的滤波核尺寸对应的窗口内应该尽可能的包含较多的非噪声部分的像素点,同时,该窗口内应当尽可能的不包含植被自身纹理对应的像素点。由于非可疑噪声区域可能为植被自身的纹理特征所在区域,因此,通过可疑噪声区域内像素点到其周围非可疑噪声区域的边缘线的最短距离获得滤波核尺寸,在该滤波核尺寸对应的窗口内,尽可能的包含较多的非噪声部分的像素点,又不包含非可疑噪声区域内的像素点。
在本实施例中,获取非可疑噪声区域的边缘像素点,根据可疑噪声区域内像素点与非可疑噪声区域的边缘像素点之间的距离最小值得到可疑噪声区域内像素点对应的特征长度,用公式表示为:
Figure 272969DEST_PATH_IMAGE022
其中,
Figure 404873DEST_PATH_IMAGE023
表示可疑噪声区域内第c个像素点对应的特征长度,以图像左上角为原点建立以像素点为单位的直角坐标系。像素点的横坐标与纵坐标分别是在其图像中所在的列数与所在行数。
Figure 946844DEST_PATH_IMAGE024
Figure 353555DEST_PATH_IMAGE025
分别表示可疑噪声区域内第c个像素点在图像上的横坐标和纵坐标,
Figure 163117DEST_PATH_IMAGE026
Figure 997081DEST_PATH_IMAGE027
分别表示非可疑噪声区域的边缘像素点的横坐标和纵坐标,min()表示求最小值函数。
同时,实施者也可根据实际情况选择其他合适的方法获取可疑噪声区域内像素点到其周围非可疑噪声区域的边缘线的最短距离。
进一步的,由于中值滤波算法中滤波核的尺寸大小均为奇数,因此需要根据特征长度确定可疑噪声区域内像素点对应的最大窗口尺寸。具体地,若可疑噪声区域内像素点的特征长度的取值为奇数,则该像素点对应的最大窗口尺寸与特征长度相同;若可疑噪声区域内像素点的特征长度的取值为偶数,则该像素点对应的最大窗口尺寸为特征长度的取值加1。
在获取可疑噪声区域内像素点对应的最大窗口尺寸后,若像素点可疑噪声区域内像素点对应的最大窗口尺寸的取值为3,则该像素点对应的最佳窗口尺寸的大小为3*3。即当最大窗口尺寸的取值为3时,滤波核的尺寸大小此时最小,且仅能使用尺寸大小为3*3的滤波核进行滤波降噪处理。
若像素点可疑噪声区域内像素点对应的最大窗口尺寸的取值大于3,则在小于或等于最大窗口尺寸的范围内选择不同的尺寸,在各尺寸对应的窗口内,根据窗口内属于可疑噪声区域内的像素点数量与窗口内像素点的总数量的比值,确定最佳窗口尺寸。
具体地,计算不同的尺寸对应的窗口内属于可疑噪声区域内的像素点数量与窗口内像素点的总数量之间的比值,获取比值的最小值对应的尺寸记为像素点的最佳窗口尺寸。即该比值表示不同尺寸对应的窗口内属于噪声部分所占窗口整体的比例,该比例越小,则说明在对应窗口内属于非噪声部分的像素点较多,则利用该尺寸进行滤波时的效果更好。
本实施例在满足可疑噪声区域内每个像素点的最大窗口尺寸内调整最终的滤波核尺寸,当某个滤波核尺寸对应的窗口内的可疑噪声部分的占比越小,则说明在该窗口内属于非噪声部分的像素点越多,则利用对应的滤波核尺寸进行降噪处理时的噪声去除效果越好。
步骤四,利用最佳窗口尺寸对初始遥感图像进行二次降噪,根据降噪后的图像获得植被的分布密度。
具体地,根据上述步骤获取的可疑噪声区域内像素点的最佳窗口尺寸对初始遥感图像进行二次降噪,在本实施例中利用中值滤波算法进行降噪处理。
需要说明的是,本发明首先通过二次去噪来避免单次去噪时滤波核尺寸选取不当的问题,其次通过分割前景、背景区域来避免二次分割造成图像细节大幅丢失的问题,然后对残留在前景区域内的噪声进行识别,并获取可疑噪声的分布、位置信息,并针对不同可疑噪声与非可疑噪声纹理的间隔距离获取自适应滤波核尺寸,利用中值滤波核进行去噪滤波,得到前景部分较为清晰的遥感图像,背景部分即使含有少量噪声,也不会对后续植被密度的评价结果产生影响。
进一步的,对二次降噪后的图像进行特征匹配,将其按照同类型的植物进行标记,以每个植物在图像中的质心点作为目标点,计算不同类型植物的分布密度,目标点的分布密度的获取方法为公知技术,在此不再过多介绍。
以上所述实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请各实施例技术方案的范围,均应包含在本申请的保护范围之内。

Claims (7)

1.一种基于园林绿化的植被覆盖密度评估方法,其特征在于,该方法包括以下步骤:
获取园林高空遥感图像,对该图像进行初步降噪,将初步降噪后的图像中的背景部分进行掩膜遮盖得到园林高空的初始遥感图像;获取初始遥感图像中所有闭合的连通域的面积,根据面积对所有闭合的连通域进行分类;
根据分类结果中各类别内闭合的连通域的数量和面积得到类别特征值,根据最大类别特征值对应的类别内连通域的面积得到各闭合的连通域的评价特征值,根据所述评价特征值确定可疑噪声区域和非可疑噪声区域;
根据可疑噪声区域内像素点到非可疑噪声区域边缘的最短距离得到可疑噪声区域内像素点对应的特征长度,根据特征长度确定像素点对应的最大窗口尺寸;在小于或等于最大窗口尺寸的范围内选择不同的尺寸,根据不同的尺寸对应的窗口内属于可疑噪声区域的像素点的数量确定最佳窗口尺寸;
利用最佳窗口尺寸对初始遥感图像进行二次降噪,根据降噪后的图像获得植被的分布密度。
2.根据权利要求1所述的一种基于园林绿化的植被覆盖密度评估方法,其特征在于,所述将初步降噪后的图像中的背景部分进行掩膜遮盖得到园林高空的初始遥感图像具体为:
对初步降噪后的图像进行边缘检测,获取各闭合边缘的边缘像素点,并获取各边缘像素点的梯度方向,根据边缘像素点的梯度方向上相邻的边缘像素点之间的像素值差值得到边缘像素点的像素差异值,根据闭合边缘的边缘像素点对应的像素差异值确定背景区域;对背景区域进行掩膜遮盖得到园林高空的初始遥感图像。
3.根据权利要求1所述的一种基于园林绿化的植被覆盖密度评估方法,其特征在于,所述类别特征值的获取方法具体为:
Figure 405686DEST_PATH_IMAGE001
其中,
Figure 866754DEST_PATH_IMAGE002
表示类别
Figure 220900DEST_PATH_IMAGE003
对应的类别特征值,
Figure 665788DEST_PATH_IMAGE004
表示类别
Figure 862283DEST_PATH_IMAGE003
内闭合的连通域的总数量,
Figure 443437DEST_PATH_IMAGE005
表示类别
Figure 293450DEST_PATH_IMAGE003
内所有闭合的连通域的面积均值,
Figure 946673DEST_PATH_IMAGE006
表示所有类别中所有闭合的连通域的面积均值,
Figure 963170DEST_PATH_IMAGE007
表示类别
Figure 382519DEST_PATH_IMAGE003
内所有闭合的连通域的面积的方差,exp()表示以自然常数e为底的指数函数。
4.根据权利要求1所述的一种基于园林绿化的植被覆盖密度评估方法,其特征在于,所述评价特征值的获取方法具体为:
Figure 341117DEST_PATH_IMAGE008
其中,
Figure 495018DEST_PATH_IMAGE009
表示第v个闭合的连通域的评价特征值,
Figure 315206DEST_PATH_IMAGE010
表示最大类别特征值对应的类别内所有闭合的连通域的面积均值,
Figure 589061DEST_PATH_IMAGE011
表示第v个闭合的连通域的面积。
5.根据权利要求1所述的一种基于园林绿化的植被覆盖密度评估方法,其特征在于,所述根据所述评价特征值确定可疑噪声区域和非可疑噪声区域具体为:
设置特征阈值,若闭合的连通域对应的评价特征值大于或等于特征阈值,则该闭合的连通域为可疑噪声区域,否则为非可疑噪声区域。
6.根据权利要求1所述的一种基于园林绿化的植被覆盖密度评估方法,其特征在于,所述根据特征长度确定该像素点对应的最大窗口尺寸具体为:
若可疑噪声区域内像素点的特征长度的取值为奇数,则该像素点对应的最大窗口尺寸与特征长度相同;若可疑噪声区域内像素点的特征长度的取值为偶数,则该像素点对应的最大窗口尺寸为特征长度的取值加1。
7.根据权利要求1所述的一种基于园林绿化的植被覆盖密度评估方法,其特征在于,所述最佳窗口尺寸的获取方法具体为:
若可疑噪声区域内像素点对应的最大窗口尺寸的取值为3,则该像素点对应的最佳窗口尺寸的大小为3*3;若像素点可疑噪声区域内像素点对应的最大窗口尺寸的取值大于3,则在小于最大窗口尺寸的范围内选择不同的尺寸,计算不同的尺寸对应的窗口内属于可疑噪声区域内的像素点数量与窗口内像素点的总数量之间的比值,获取比值的最小值对应的尺寸记为像素点的最佳窗口尺寸。
CN202211358958.9A 2022-11-02 2022-11-02 一种基于园林绿化的植被覆盖密度评估方法 Active CN115409742B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211358958.9A CN115409742B (zh) 2022-11-02 2022-11-02 一种基于园林绿化的植被覆盖密度评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211358958.9A CN115409742B (zh) 2022-11-02 2022-11-02 一种基于园林绿化的植被覆盖密度评估方法

Publications (2)

Publication Number Publication Date
CN115409742A true CN115409742A (zh) 2022-11-29
CN115409742B CN115409742B (zh) 2023-01-31

Family

ID=84169161

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211358958.9A Active CN115409742B (zh) 2022-11-02 2022-11-02 一种基于园林绿化的植被覆盖密度评估方法

Country Status (1)

Country Link
CN (1) CN115409742B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116740580A (zh) * 2023-08-16 2023-09-12 山东绿博园市政工程有限公司 基于遥感技术的园林工程数据处理方法
CN117132423A (zh) * 2023-08-22 2023-11-28 珠海市经典电子有限公司 园区管理系统
CN117496448A (zh) * 2024-01-02 2024-02-02 深圳市旭仓科技有限公司 一种建筑施工安全智能监控系统及方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105243649A (zh) * 2015-11-09 2016-01-13 天津大学 基于二次噪点检测的图像降噪方法
CN106373096A (zh) * 2016-08-30 2017-02-01 电子科技大学 一种多特征权值自适应的阴影消除方法
CN107038688A (zh) * 2017-03-01 2017-08-11 武汉科技大学 基于Hessian矩阵的图像噪点检测及去噪方法
US20210090219A1 (en) * 2020-01-13 2021-03-25 Jiangxi University Of Science And Technology Salt and pepper noise filtering method and device based on morphological component analysis
CN113313641A (zh) * 2021-04-28 2021-08-27 北京理工大学 一种自适应中值滤波的ct图像去噪方法
CN113724167A (zh) * 2021-09-03 2021-11-30 中国矿业大学 一种高清视频信息的自适应采集及图像处理方法
CN115239736A (zh) * 2022-09-23 2022-10-25 海安玻克超硬材料有限公司 用于金刚石磨轮生产的磨料层混料质量监测方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105243649A (zh) * 2015-11-09 2016-01-13 天津大学 基于二次噪点检测的图像降噪方法
CN106373096A (zh) * 2016-08-30 2017-02-01 电子科技大学 一种多特征权值自适应的阴影消除方法
CN107038688A (zh) * 2017-03-01 2017-08-11 武汉科技大学 基于Hessian矩阵的图像噪点检测及去噪方法
US20210090219A1 (en) * 2020-01-13 2021-03-25 Jiangxi University Of Science And Technology Salt and pepper noise filtering method and device based on morphological component analysis
CN113313641A (zh) * 2021-04-28 2021-08-27 北京理工大学 一种自适应中值滤波的ct图像去噪方法
CN113724167A (zh) * 2021-09-03 2021-11-30 中国矿业大学 一种高清视频信息的自适应采集及图像处理方法
CN115239736A (zh) * 2022-09-23 2022-10-25 海安玻克超硬材料有限公司 用于金刚石磨轮生产的磨料层混料质量监测方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
ALESSIO BARDUCCI ET AL: "Object Recognition in Floor Plans by Graphs of White Connected Components", 《21ST INTERNATIONAL CONFERENCE ON PATTERN RECOGNITION (ICPR 2012)》 *
于泳波等: "基于图像连通域的桥梁裂缝提取方法", 《公路交通科技》 *
周佳慧等: "基于方向信息自适应2次噪点检测的降噪方法", 《南开大学学报(自然科学版)》 *
姚薇等: "矿山遥感图像自适应加权改进中值滤波算法", 《金属矿山》 *
王昕等: "视网膜血管分割算法研究", 《吉林大学学报》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116740580A (zh) * 2023-08-16 2023-09-12 山东绿博园市政工程有限公司 基于遥感技术的园林工程数据处理方法
CN117132423A (zh) * 2023-08-22 2023-11-28 珠海市经典电子有限公司 园区管理系统
CN117132423B (zh) * 2023-08-22 2024-04-12 深圳云创友翼科技有限公司 园区管理系统
CN117496448A (zh) * 2024-01-02 2024-02-02 深圳市旭仓科技有限公司 一种建筑施工安全智能监控系统及方法
CN117496448B (zh) * 2024-01-02 2024-04-26 深圳市旭仓科技有限公司 一种建筑施工安全智能监控系统及方法

Also Published As

Publication number Publication date
CN115409742B (zh) 2023-01-31

Similar Documents

Publication Publication Date Title
CN115409742B (zh) 一种基于园林绿化的植被覆盖密度评估方法
CN110414411B (zh) 基于视觉显著性的海面船只候选区域检测方法
CN111260616A (zh) 一种基于Canny算子二维阈值分割优化的绝缘子裂纹检测方法
CN109191432B (zh) 基于域变换滤波多尺度分解的遥感图像云检测方法
CN108389175B (zh) 融合变差函数和颜色衰减先验的图像去雾方法
CN109255757B (zh) 一种机器视觉自然放置葡萄串果梗区域分割方法
CN114118144A (zh) 抗干扰的航空遥感图像阴影精准检测方法
CN108229342B (zh) 一种海面舰船目标自动检测方法
CN112149543B (zh) 一种基于计算机视觉的建筑扬尘识别系统与方法
CN114372983A (zh) 一种基于图像处理的屏蔽盒涂装质量检测方法及系统
CN109389609B (zh) 基于fart神经网络的交互自反馈红外目标检测方法
CN113221881B (zh) 一种多层级的智能手机屏幕缺陷检测方法
CN109871900A (zh) 一种基于图像处理的复杂背景下苹果的识别定位方法
CN111353371A (zh) 基于星载sar影像的海岸线提取方法
CN110660065A (zh) 一种红外故障检测识别算法
CN114842337A (zh) 基于深度学习与多维信息融合聚类的果实采摘点识别方法
CN117575953B (zh) 一种高分辨率林业遥感图像细节增强方法
CN109446896B (zh) 基于图像处理的无人机信号识别检测算法
CN111667509B (zh) 目标与背景颜色相似下的运动目标自动跟踪方法及系统
CN110765887A (zh) 隧道衬砌裂纹自动识别技术及检测方法
CN113705441A (zh) 协同多光谱和sar影像的高时空分辨率地表水体提取方法
CN113096103A (zh) 一种面向放空火炬的烟气智能影像感知方法
CN110047064B (zh) 一种马铃薯疮痂检测方法
CN111882529B (zh) 显示屏Mura缺陷检测方法及装置
CN114119634A (zh) 结合植被剔除和图像特征一致性约束的建筑物自动提取方法及系统

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