留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

1种基于全局对比度的害虫收集方法

柳懿祥 汪杭军 徐铁平

方庆, 谭菊荣, 许惠春, 等. 珍稀濒危植物细果秤锤树群落物种组成与生态位分析[J]. 浙江农林大学学报, 2022, 39(5): 931-939. DOI: 10.11833/j.issn.2095-0756.20220197
引用本文: 柳懿祥, 汪杭军, 徐铁平. 1种基于全局对比度的害虫收集方法[J]. 浙江农林大学学报, 2021, 38(2): 369-376. DOI: 10.11833/j.issn.2095-0756.20200318
FANG Qing, TAN Jurong, XU Huichun, et al. Species composition and niche of Sinojackia microcarpa, a rare and endangered plant[J]. Journal of Zhejiang A&F University, 2022, 39(5): 931-939. DOI: 10.11833/j.issn.2095-0756.20220197
Citation: LIU Yixiang, WANG Hangjun, XU Tieping. A pest collection method based on global contrast[J]. Journal of Zhejiang A&F University, 2021, 38(2): 369-376. DOI: 10.11833/j.issn.2095-0756.20200318

1种基于全局对比度的害虫收集方法

DOI: 10.11833/j.issn.2095-0756.20200318
基金项目: 浙江省基础公益研究计划项目(LGN19C140006);浙江省自然科学基金资助项目(LY16C160007);绍兴市科技计划项目(2018C20013)
详细信息
    作者简介: 柳懿祥(ORCID: 0000-0002-7561-7334),从事图像处理与模式识别研究。E-mail: 616731542@qq.com
    通信作者: 汪杭军(ORCID: 0000-0002-7698-1226),教授,博士,从事林业图像处理与模式识别研究。E-mail: whj@zafu.edu.cn
  • 中图分类号: S763

A pest collection method based on global contrast

  • 摘要:   目的  目前利用测报灯,通过灯光诱捕昆虫,并由计算机完成昆虫图像的采集、计数和识别已逐步成为害虫测报的重要方法。为了减少昆虫在采样盘上重叠造成的计数和识别误差,基于害虫图像,根据昆虫密度研究采样盘中昆虫的收集方法,从而提高采集效率和精度。  方法  根据昆虫在采样盘上姿态特点,提出基于全局对比度的图像分割方法,结合阈值迭代分割获得昆虫区域,计算昆虫比例,并控制采样盘翻转完成对昆虫的收集。  结果  通过对5种害虫的实际图像进行的试验表明:与水平集、大津法(OTSU)、阈值迭代法和基于直方图对比度的显著性检测(HC)4种算法相比,本研究方法在准确率和召回率上均提高10%以上,取得了较好的结果;同时,在分割速度上比水平集快3倍,与阈值和HC算法基本持平。  结论  基于全局对比度的分割方法简单、高效,在害虫自动测报中具有较高的实际应用价值。图7表1参17
  • 开展珍稀濒危植物的群落生态学研究有助于野生植物资源的保护、恢复和可持续更新。群落生态学研究一般通过探究物种的分布范围、群落结构及种内与种间联结关系等,揭示群落生活史、适应性、生长趋势等[1-3]。物种组成与群落结构在一定程度上展现植物对资源的利用能力和群落的稳定程度[4]。汪国海等[5]通过研究濒危植物单性木兰Kmeria septentrionalis的群落结构与空间分布格局,探究其聚集方式和传播途径。濒危物种的生态位宽度与群落总体关联度能够反映物种间的相互关系(竞争或促进作用)及对生境条件的适应状况和资源利用情况等[6-8]。刘万德等[9]对藤枣Eleutharrhena macrocarpa的生境特征和种间联结研究发现:藤枣与下层木呈极显著负相关,减少群落内下层木可以促进藤枣群落可持续生长[3, 9-11]。杨国平等[12]通过建立预测景东翅子树Pterospermum kingtungense群落动态的Lefkovitch矩阵模型,探究濒危物种在特定的小生境片段中的分布区间。因此,基于群落生态学的研究方法,有助于全面评估珍稀濒危物种的内外致濒因子,缓解其濒危态势,实现有效的拯救保护[10-11]

    细果秤锤树Sinojackia microcarpa为中国特有的极小群落野生植物,多分布在浙江临安、建德等地,处于极度濒危和受胁迫状态[13-17]。目前,对秤锤属Sinojackia的研究相对较多。杨国栋等[18]采用生态学理论结合自组织特征映射网络(SOM)方法,划分了野生秤锤树群落的群丛类型。徐惠明等[19]分析了狭果秤锤树S. rehderiana的群落年龄结构,发现该群落具有良好的更新潜力。周赛霞等[20]研究发现:受密度制约或种子扩散限制等,狭果秤锤树的空间聚集分布趋势逐渐减弱。秤锤属物种多表现出竞争能力相对较弱,对外界干扰的响应较为显著[18-19]。本研究通过对细果秤锤树群落的长期动态监测,分析细果秤锤树群落的物种组成、生态位宽度及其与主要树种的种间关联,揭示细果秤锤树的生境适应性与竞争强度,有助于在就地、迁地保护回归实践中建立适宜的生存环境。

    浙江省建德市属亚热带北缘季风气候,雨量充沛,四季分明,年平均气温为17.4 ℃。土壤类型以凝灰岩发育的红壤、黄棕色壤土为主,土层浅薄且质地较为疏松,钱塘江水系中上游,境内以低山丘陵地貌为主。细果秤锤树集中分布于浙江省建德市建德林场乌石滩林区(29°32′56″~29°35′43″N,119°33′08″~119°34′05″E),主要分布在林区乌石滩、富家坞和灵山顶,海拔为23~429 m。多生长在岩石裸露率较大的山谷溪沟边的灌丛林中,呈条带状分布,群落生境数年前遭受人为砍伐干扰较严重。

    细果秤锤树为典型极小群落野生植物,残存数量较少,因适存的小流域生境使得群落呈带状分布,样地设置受限。2020年8—9月,在全面踏查细果秤锤树野生群落的基础上,参照热带森林科学研究中心(CTFS)的样地建设技术规程,建立0.18 hm2的固定监测样地。使用全站仪在乌石滩、富家坞和灵山顶分别设置3个典型样方开展群落调查,共计9个10 m×20 m样方;在每个样方内设置3个5 m×5 m的下层木样地以及3个1 m×1 m的草本层样地。开展树种定位、地形测定(海拔、经纬度、坡向坡位等)、生境因子测定(土壤理化性质等)。

    本研究计算上层木与下层木的物种重要值。上层木重要值=(相对多度+相对频度+相对显著度)/3;下层木重要值=(相对多度+相对频度)/2;相对多度=(某种植物的数量/样地植物的总数量)×100%;相对频度=(某种植物的频度/样地所有植物物种的频度总和)×100%;相对优势度=(某种植物的胸高断面积之和/样地所有物种的胸高断面积之和)×100%。

    物种生态位特征主要采用Levins指数、Shannon-Wiener指数[21-23]反映生态位宽度,Schoener生态位相似性[24-25]与Pianka生态位重叠指数[26]反映生态相似与重叠程度。种间联结分析主要采用总体联结指数[6, 8]、卡方检验(χ2)、联结系数(AC)[24]和Pearson相关系数[8, 22]探究物种间关联性。采用R 4.1.0中spaa包计算生态位宽度、生态位相似性和生态位重叠程度、χ2检验、Pearson相关系数检验结果。

    细果秤锤树总计509株,其中富家坞分布个体数量最多(243株),灵山顶最少(71株)。群落里单丛萌蘖枝干中的最大胸径为8.10 cm,平均树高为5.40 m(表1)。乌石滩、富家坞、灵山顶细果秤锤树群落的胸径变异系数分别为34%、33%和33%,均表现为较低变异性。

    表 1  细果秤锤树群落资源组成
    Table 1  Composition of population resources of S. microcarpa
    分布区数量/
    胸径/
    cm
    树高/
    m
    胸径变异
    系数/%
    树高变异
    系数/%
    乌石滩1953.07±1.055.00±1.873438
    富家坞2433.05±1.025.40±1.983341
    灵山顶712.95±0.984.90±2.413354
      说明:胸径和树高数值为平均值±标准差
    下载: 导出CSV 
    | 显示表格

    细果秤锤树分布在海拔23~429 m的区域(表2表3),乌石滩和富家坞受人工干预程度较高,存在人为滥砍及割灌除草等抚育过程。土壤呈较疏松多孔的黏质土,土壤容重为1.06~1.19 g·cm−3,pH为4.72~5.79,偏酸性土壤,有效磷和速效钾偏低。细果秤锤树群落土壤有机质、氮、磷、钾及其速效成分中等,土壤养分条件一般。

    表 2  细果秤锤树群落生境调查
    Table 2  Environmental survey of S. microcarpa population
    分布区样地海拔/m纬度(N)经度(E)坡向群落特征
    乌石滩 P1 58 29°34′16″ 119°33′10″ 西 樟树Cinnamomum camphora-板栗Castanea mollissima混交林
    P2 45 29°34′18″ 119°3360″ 西 板栗林
    P3 64 29°34′17″ 119°33′00″ 东北 板栗林
    富家坞 P4 58 29°34′57″ 119°33′42″ 东南 柏木Cupressus funebris-南酸枣Choerospondias axiliaris混交林
    P5 95 29°34′57″ 119°33′36″ 东南 柏木林
    P6 128 29°35′20″ 119°33′24″ 柏木-拟赤杨Alniphyllum fortunei混交林
    灵山顶 P7 190 29°35′35″ 119°33′52″ 东北 樟树林
    P8 384 29°35′11″ 119°33′11″ 东北 毛竹Phyllostachys edulis
    P9 396 29°35′40″ 119°33′10″ 东北 毛竹林
    下载: 导出CSV 
    | 显示表格
    表 3  细果秤锤树群落的生境因素
    Table 3  Habitat factors of S. microcarpa
    分布区海拔/m土壤容重/
    (g·cm−3)
    土壤pH土壤有机
    质/(g·kg−1)
    土壤总孔
    隙度/%
    土壤碱解氮/
    (mg·kg−1)
    土壤有效磷/
    (mg·kg−1)
    土壤速效钾/
    (mg·kg−1)
    乌石滩 70±26 a 1.01±0.10 a 5.46±0.20 a 38.84±3.66 a 61.74±3.67 a 103.41±3.08 a 6.23±0.82 a 82.46±3.22 a
    富家坞 109±39 a 1.12±0.06 a 5.47±0.43 a 40.76±1.22 a 57.72±2.25 a 97.61±6.90 a 5.79±1.26 a 82.93±6.82 a
    灵山顶 370±110 a 1.07±0.09 a 5.23±0.15 a 45.74±3.42 a 59.72±3.44 a 107.71±8.72 a 5.54±1.45 a 95.48±14.02 a
    变化范围 23~429 1.00~1.19 4.72~5.79 36.81~48.38 55.20~62.42 91.04~113.67 5.30~7.84 75.69~102.80
      说明:数值为平均值±标准差。同列不同小写字母表示同一指标不同分布区之间差异显著(P<0.05)
    下载: 导出CSV 
    | 显示表格

    细果秤锤树样地内共记录到胸径≥1 cm的木本植物401株,隶属于35科50属51种。其中优势科有樟科Lauraceae (5属6种)、山茶科Theaceae (3属4种)、壳斗科Fagaceae (3属3种)、马鞭草科Verbenaceae (3属3种)、安息香科Styracaceae (2属3种)、大戟科Euphorbiaceae (2属2种)、金缕梅科Hamamelidaceae (2属2种)、漆树科Anacardiaceae (2属2种)、茜草科Rubiaceae (2属2种)、榆科Ulmaceae (2属2种)。樟树的平均胸径最大,达30.8 cm,有22株;平均胸径较大的树种有臭椿Ailanthus altissima、枫香Liquidambar formosana、柏木、南酸枣和毛竹。

    样地中重要值≥1%的上层木物种共16种,重要值排前4位的物种是毛竹、柏木、板栗和细果秤锤树,这4个物种重要值之和为49.85%,是群落优势树种(表4)。下层中阔叶箬竹Indocalamus latifolius的重要值最高,为15.48%;重要值排前3位的物种有水团花Adina pilulifera、毛花连蕊茶Camellia fraterna和细果秤锤树(表5)。细果秤锤树在上、下木层中重要值分别为9.50%和4.60%,是主要建群种之一。

    表 4  细果秤锤树群落上层木主要物种的重要值和生态位宽度
    Table 4  Important values and niche breadth of the dominant species in upper wood layer of S. microcarpa community
    编号物种重要值/
    %
    生态位宽度编号物种重要值/
    %
    生态位宽度
    Levins
    指数
    Shannon-Wiener
    指数
    Levins
    指数
    Shannon-Wiener
    指数
    1 毛竹 19.63 1.96 0.68 11 杉木 2.00 1.78 0.63
    2 柏木 10.84 2.48 1.00 12 黄檀 1.95 2.29 0.90
    3 板栗 9.88 2.80 1.13 13 白花泡桐 1.70 1.00 0.00
    4 细果秤锤树 9.50 5.87 1.92 14 盐肤木 1.51 1.00 0.00
    5 樟树 8.44 1.82 0.64 15 木油桐 1.27 1.96 0.68
    6 南酸枣 2.75 1.83 0.80 16 大叶白纸扇 1.21 2.00 0.69
    7 拟赤杨 2.34 1.95 0.68 17 厚壳树 0.99 1.00 0.00
    8 枫香 2.32 1.00 0.00 18 臭椿 0.96 1.00 0.00
    9 木蜡树 2.18 2.70 1.05 19 檵木 0.88 1.63 0.00
    10 棕榈 2.09 2.78 1.06
      说明:木蜡树Toxicodendron sylvestr;棕榈Trachycarpus fortunei;杉木Cunninghamia lanceolata;黄檀Dalbergia hupeana;白花泡桐Paulownia fortunei;盐肤木Rhus chinensis;木油桐Vernicia montana;大叶白纸扇Mussaenda shikokiana;厚壳树Ehretia thysiflora;檵木Loropetalum chinensis
    下载: 导出CSV 
    | 显示表格
    表 5  细果秤锤树群落下层木主要物种的重要值和生态位宽度值
    Table 5  Important value and niche breadth of the dominant species in lower wood layer of S. microcarpa community
    编号物种重要值/%生态位宽度编号物种重要值/%生态位宽度
    Levins
    指数
    Shannon-Wiener
    指数
    Levins
    指数
    Shannon-Wiener
    指数
    1 阔叶箬竹 15.48 1.98 0.84 9 短柄枹栎 2.41 1.84 0.65
    2 水团花 8.45 3.43 1.30 10 紫麻 2.30 1.08 0.16
    3 细果秤锤树 4.60 6.82 2.00 11 木荷 1.86 1.00 0.00
    4 毛花连蕊茶 4.58 4.95 1.77 12 华箬竹 1.63 1.00 0.00
    5 4.44 4.09 1.73 13 杉木 1.59 1.28 0.38
    6 檵木 3.05 4.39 1.60 14 海金子 1.58 1.92 0.74
    7 窄基红褐柃 2.98 1.00 0.00 15 黄檀 1.54 2.81 1.06
    8 杭州榆 2.69 1.00 0.00 16 朱砂根 1.45 3.90 1.57
      说明:窄基红褐柃Eurya rubiginosa var. attenuata;杭州榆Ulmus changii;短柄枹栎Quercus glandulifera;木荷Schima superba;华箬竹Sasa sinica;朱砂根Ardisia crenata
    下载: 导出CSV 
    | 显示表格

    细果秤锤树具有最大的生态位宽度,Levins的生态位宽度指数及Shannon-Wiener的生态位宽度指数在上层木中分别为5.87%和1.92%(表5),板栗、棕榈、木蜡树与柏木的生态位宽度依次降低。细果秤锤树在上层木林层与下层木林层中生态位宽度差异不明显,说明细果秤锤树的种对竞争具有一定优势,在所调查的小流域生境中具有较强的适应能力,分布幅度较广。

    细果秤锤树群落上层木物种生态位相似性和生态位重叠值最大均为盐肤木-臭椿(表6)。细果秤锤树与上层优势树种樟树生态相似性值最高(0.62),白花泡桐次之(0.59)。生态位宽度较大的柏木和黄檀的生态位相似性达0.65,而生态位宽度较窄的枫香和臭椿的生态位相似性为0,说明生态位相似性与生态位宽度有一定关联。生态位重叠值在0.8~1.0的种对有杉木-盐肤木和南酸枣-枫香,大于0.5的种对有39对(占20.53%),其中生态位重叠值小于0.1的种对共有90对(占47.37%)。上层木树种间生态位重叠值总体偏低,对资源利用的利用策略存在差异。细果秤锤树与樟树(0.62)和黄檀(0.59)具有较大的生态位重叠,存在较大的生态和资源利用相似性。

    表 6  细果秤锤树群落上层木主要优势种间的生态位相似性比例和生态位重叠指数
    Table 6  Niche similarity and niche overlap of dominant plant species in S. microcarpa community in the upper wood layer
    编号12345678910111213141516171819
    1000.06000.14000.270.5800.070.420000.420
    200.040.3400.230.410.450.260.4300.52000.190.740.1900.55
    300.020.350.150.410.040.0400.0400.040.5200.040.040.0400
    40.090.540.470.420.430.180.150.130.250.200.210.380.040.020.210.160.040.19
    5000.110.620.4000000.4200.340000.5200
    600.150.490.480.680.230.230.090.140.400.140.2700.140.140.5400
    70.120.460.010.3300.220.570.710.290.140.1500.140.150.150.150.140
    800.490.050.2200.320.720.420.4900.58000.580.420.4800
    900.4100.3200.170.960.590000000000
    100.270.580.050.2500.220.240.6500.270.7300.270.490.660.480.270.24
    110.64000.320.620.540.12000.2800.070.38000.420.380
    1200.650.060.2400.250.180.7300.890000.680.750.4800.32
    130.0900.630.590.300.3900000.040000000
    140.59000.09000.19000.450.63000001.000
    1500.300.060.0400.270.200.8100.8000.90000.420.4800
    1600.870.040.4000.160.120.4800.7900.88000.590.4200.58
    1700.200.040.280.660.750.130.5500.540.520.61000.670.4000
    180.59000.09000.19000.450.63001.000000
    1900.8600.46000000.3900.430000.8000
      说明:编号所代表物种见表4。对角线下方为生态位相似性,对角线上方为生态位重叠值
    下载: 导出CSV 
    | 显示表格

    下层木物种生态位相似性为0~0.96,生态位重叠为0~0.10,最大值种对均为海金子Pittosporum illiciodes-紫麻Oreocnide frutescens。细果秤锤树与下层优势树种檵木生态相似性值最高(0.86);与水团花(0.51)和茶Camellia sinensis (0.48)具有较大生态重叠(表7)。下层木主要物种生态位重叠平均值为0.23,且多数种对的生态位重叠在其平均值附近,表明下层木主要物种的竞争关系相对稳定。

    表 7  细果秤锤树群落下层木主要优势种间的生态位相似性比例和生态位重叠指数
    Table 7  Niche similarity and niche overlap of dominant plant species in S.microcarpa community in the lower wood layer
    编号12345678910111213141516
    10.090.320.210.160.200.340.020.25000.090.09000
    20.040.510.430.530.520.290.46000.130.560.19000
    30.300.650.360.480.630.330.380.230.010.100.0700.0100
    40.340.580.380.590.530.460.5400.080.390.240.110.110.170.08
    50.180.580.530.640.550.270.4700.340.190.210.110.340.060
    60.190.650.860.520.620.440.340.040.040.240.130.130.070.090.12
    70.310.270.320.470.210.360.140.3400.390.290.240.040.170.16
    80.020.560.470.750.480.410.1500.050.180.1700.060.030.01
    90.3500.41000.090.7000000000
    10000.020.200.750.0700.0800.01000.9600
    1100.120.030.410.100.290.440.1000.010.1100.050.220.64
    120.060.770.120.440.320.120.330.24000.140.35000
    130.130.3500.290.230.260.500.010000.48000
    14000.020.210.750.080.010.0801.000.050000.04
    150000.450.130.190.350.05000.320000
    160000.2100.240.330.02000.93000.040
      说明:编号所代表物种见表5。对角线下方为生态位相似性,对角线上方为生态位重叠值
    下载: 导出CSV 
    | 显示表格

    细果秤锤树群落上层木12个优势种间总体联结性方差比率为1.23,大于1,即种间存在一定程度正联结;其显著检验统计量为11.05,高于χ2分布临界值,表明上层木群落间总体上呈显著的正联结关系。下层木12个优势种间总体联结性方差比率为0.58,小于1,即种间存在一定程度负联结;其显著检验统计量为5.19,介于χ2分布临界值之间,即下层木12个优势种间呈不显著负联结关系。

    χ2检验主要反映不同种对之间联结的显著度。联结系数检验结果显示:上层和下层各12个优势木中,正、负联结种对数相接近。细果秤锤树群落上层木中正、负联结的种对分别为27和28个(各占种对数的40.91%和42.42%),正负关联比为0.96∶1.00。种对间总体显著率为12.12%,种间联结较松散,无联结的种对占16.67%,细果秤锤树与其他种之间都不存在联结性。下层木种对联结显著度的分布大致与上层木相似,正负关联比0.83∶1.00。细果秤锤树与水团花呈显著正联结关系。细果秤锤树-阔叶箬竹、细果秤锤树-茶、细果秤锤树-檵木、细果秤锤树-窄基红褐柃表现出极显著负关联(表8)。

    表 8  细果秤锤树群落12个优势种χ2检验、联结系数(AC)及Pearson相关检验结果
    Table 8  Result of χ2 test, association coefficient (AC) and Pearson correlation coefficient of the 12 dominant species in S. microcarpa community
    检验方法检验结果数值范围上层木下层木检验方法检验结果数值范围上层木下层木
    种对数占比/%种对数占比/%种对数占比/%种对数占比/%
    χ2正相关P≤0.010000AC负相关−0.2≤AC<023.0323.03
    0.01<P≤0.0523.03710.61−0.6≤AC<−0.234.5434.54
    P>0.052537.882233.33AC≤−0.62334.853045.46
    无关联χ2=01116.6723.03
    负相关P≤0.010000Pearson
    相关检验
    正相关P≤0.011319.7000
    0.01<P≤0.0569.0957.580.01<P≤0.050000
    P>0.052233.333045.45P>0.052537.883146.97
    无关联0<P<0.200000
    AC正相关AC≥0.6913.642030.30负相关P≤0.010000
    0.2≤AC<0.6812.1223.030.01<P≤0.050000
    0<AC<0.2812.12710.61P>0.052842.423553.03
    无关联AC =01319.7023.03
    下载: 导出CSV 
    | 显示表格

    上层木中总体显著率为19.70%(极显著正关联13个,P<0.01),不显著(P>0.05)正关联25个,占37.88%;不显著负关联28个,占比42.42%。细果秤锤树与其他树种为无联结关系,整个细果秤锤树群落处于优势发展趋势(表8)。下层木中总体显著率为0,不显著正关联31个,占46.97%;不显著负关联35个,占53.03%。细果秤锤树与水团花、毛花连蕊茶、杭州榆、短柄枹栎呈不显著正关联,与阔叶箬竹、茶、檵木、窄基红褐柃呈不显著负关联。

    建德市野生细果秤锤树群落动态监测样地内树种组成相对简单,细果秤锤树多生长在次生常绿阔叶林和针阔混交林中,群落优势树种主要为毛竹、柏木、板栗和细果秤锤树。这与秤锤属调查样地内的物种组成及数量相类似[13, 15-16]。调查发现:细果秤锤树群落中缺乏小径级个体或幼苗,这可能是因为秤锤属的种子萌发困难或遭受了人为的抚育等干扰,影响了幼苗的更新[13-14]。细果秤锤树是小流域生境群落中的优势种,早期生长喜较为荫蔽的环境,群落中高大上层木树种如樟树、毛竹、柏木等可在其幼苗更新时期起到遮光作用,以保护幼苗不受高温、强光照影响。在细果秤锤树生长后期,对光照需求增强,可间伐上层木,对高度接近细果秤锤树的树种进行一定程度的抚育,降低群落郁闭度[12, 16-17]

    生态位宽度作为植物群落的环境适应力和资源利用能力的衡量性指标,值越大,反映物种适应能力越强,在群落中更具优势[22, 27]。细果秤锤树在群落物种中重要值排在第4位,但生态位宽度却排在首位。可能是其喜光、耐贫瘠、喜微酸性土壤等生长特性有利于细果秤锤树在小溪流水域附近广泛分布。细果秤锤树的生态位宽度较大还可能与本研究的样地设置有关。本研究以细果秤锤树生长的位置为核心展开设置并调查,且呈聚集分布均匀的群落使得其占较大资源位或较大资源量,与极小群落植物圆叶玉兰Magnolia sinensis[28]、小花木兰Oyama sieboldii[29]、缙云秋海棠Begonia jinyunensis[30]在所处群落中生态位宽度均较大这一研究结果相同,表明在该分布点的研究区域生境条件下,生态位宽度大小与细果秤锤树致濒机制无必然联系。研究中有一些物种的生态位宽度大小排序与其重要值大小排序不同,如樟树、南酸枣等,这说明生态位宽度和重要值在物种之间的表现方式略有不同且并无显著关联性。

    生态位相似性特征反映种间资源利用的相似程度,重叠值特征衡量生态位相似的树种在特定空间环境下资源利用的差异性,两者结合衡量种间资源竞争程度[31-33]。细果秤锤树与上层优势树种樟树和黄檀的生态相似性与生态重叠性均最高。可能是因为樟树、黄檀是对环境适应性广泛的泛化种,也可能是适合调查区域环境的特化种,因此出现与细果秤锤树较高的生态位重叠值,也表明这些种对间生态学特性比较一致,或者对生境的要求比较相似[8]。一般来说,当多个物种同时具有较大的生态位宽度时,它们之间存在较高生态位重叠的可能性更大[21]。但是,具有较大生态位宽度的物种也可能与较小生态位宽度的物种间存在较大的生态位重叠[21, 31]。这是因为细果秤锤树与水团花、毛花连蕊茶为中生植物,在资源有限的条件下,它们对资源环境的竞争比较大,且对资源的利用和需求相近[32],因此,它们之间的联系也更为紧密,具有较高的生态位重叠[22, 26]。且细果秤锤树所在群落中物种之间的生态位重叠程度总体偏低,说明细果秤锤树群落中大多物种对资源利用的相似程度降低,物种之间竞争较弱,生态位可通过产生分化来降低种间竞争使得物种间在群落的结构与功能上互补且稳定[7, 22]。本研究发现:细果秤锤树群落大部分种对间的相关性比较弱,表明物种联结性较弱。种间负联结关系占主导,但大部分优势种种对间关联性比较低,说明样地中的不同物种间不存在紧密的相互关系,缺乏竞争或相互促进的趋势,物种间具有独立性,受外界的干扰较小[30]

    细果秤锤树群落中物种组成较为简单,群落结构相对单一,细果秤锤树群落幼树较少,更新相对较差。细果秤锤树生态位宽度最大,在时空上占据着优势地位,属于稍耐阴、耐贫瘠、适应力较强的植物,能更好利用资源和空间。调查样地中多数树种生态位重叠度较高,大部分物种间的竞争较强,对资源利用的相似程度高。树种间不存在较显著的种间相关联结,植物种间缺乏较强的相互依赖或竞争趋势。本研究明确了细果秤锤树生存的独特环境结构和群落间相互关系,对维持其野生群落的幼苗更新和群落规模增长具有重要作用。

  • 图  1  方法流程图

    Figure  1  Algorithm flow chart

    图  2  害虫收集装置示意图

    Figure  2  Schematic diagram of pest collection device

    图  3  白底样本分割效果对比

    Figure  3  Comparison of sample segmentation effects on white background

    图  4  OTSU算法与本研究算法对阴影分割的比较

    Figure  4  Comparison of shadow segmentation between OTSU and research algorithm

    图  5  HC算法与本研究算法对害虫细节分割的比较

    Figure  5  Comparison of HC algorithm and research algorithm on pest detail segmentation

    图  6  红底样本分割效果对比

    Figure  6  Comparison of sample segmentation effects on red background

    图  7  HC算法与本研究算法对深色害虫分割的比较

    Figure  7  Comparison of HC algorithm and research algorithm on dark pest segmentation

    表  1  准确率、召回率和分割速度对比

    Table  1.   Comparison of accuracy recall, and segmentation speed

    分割算法准确率召回率分割速度
    白底样本红底样本白底样本红底样本白底样本红底样本
    水平集算法0.720.740.620.641.684.35
    OTSU算法0.700.770.630.700.360.84
    阈值迭代算法0.830.830.810.690.560.87
    HC算法0.860.880.860.630.430.98
    本研究算法0.920.950.950.890.661.05
    下载: 导出CSV
  • [1] 张芳群. 基于Web技术的虫害预测系统的研究[D]. 杭州: 浙江理工大学, 2017.

    ZHANG Fangqun. Research on the Pest Forecasting System based on Web Technology[D]. Hangzhou: Zhejiang Sci-Tech University, 2017.
    [2] 陈梅香, 杨信廷, 石宝才, 等. 害虫自动识别与计数技术研究进展与展望[J]. 环境昆虫学报, 2015, 37(1): 176 − 183.

    CHEN Meixiang, YANG Xinting, SHI Baocai, et al. Research progress and prospect of technologies for automatic identifying and counting of pests [J]. J Environ Entomol, 2015, 37(1): 176 − 183.
    [3] 张书平, 余燕, 毕守东, 等. 灾变模型在马尾松毛虫幼虫发生量预报中的应用[J]. 浙江农林大学学报, 2020, 37(1): 93 − 99.

    ZHANG Shuping, YU Yan, BI Shoudong, et al. A catastrophe model to forecast larvae occurrence of Dendrolimus punctatus [J]. J Zhejiang A&F Univ, 2020, 37(1): 93 − 99.
    [4] 张鸣华, 丁鑫, 刘持标. 基于二值化图像识别的新型农业害虫诱杀监测系统研究[J]. 南昌航空大学学报(自然科学版), 2019, 33(1): 74 − 79.

    ZHANG Minghua, DING Xin, LIU Chibiao. Research on new agricultural pest trapping and killing monitoring system based on binarized image recognition [J]. J Nanchang Hangkong Univ Nat Sci, 2019, 33(1): 74 − 79.
    [5] 陈京, 刘德营, 谢堂胜, 等. 田间稻飞虱图像远程实时采集系统的研制[J]. 湖南农业大学学报(自然科学版), 2016, 42(6): 693 − 698.

    CHEN Jing, LIU Deying, XIE Tangsheng, et al. A remote real-time acquisition system for rice plant hopper images in the fields [J]. J Hunan Agric Univ Nat Sci, 2016, 42(6): 693 − 698.
    [6] SOLIS-SÁNCHEZ L O, GARCÍA-ESCALANTE J J, CASTAÑEDA-MIRANDA R, et al. Machine vision algorithm for whiteflies(Bemisia tabaci Genn)scouting under greenhouse environment [J]. J Appl Entomol, 2009, 133(7): 546 − 552.
    [7] NING Jifeng, ZHANG Lei, ZHANG David, et al. Interactive image segmentation by maximal similarity based region merging [J]. Pattern Recognition, 2009, 43(2): 445 − 456.
    [8] 吕金娜. 基于LAB空间和自适应聚类的害虫图像分割方法[J]. 河南科技学院学报(自然科学版), 2016, 44(1): 57 − 61.

    LÜ Jinna. Image segmentation based on LAB color space adaptive clustering algorithm for pest image segmentation [J]. J Henan Inst Sci Technol Nat Sci Ed, 2016, 44(1): 57 − 61.
    [9] 杨信廷, 刘蒙蒙, 许建平, 等. 自动监测装置用温室粉虱和蓟马成虫图像分割识别算法[J]. 农业工程学报, 2018, 34(1): 164 − 170.

    YANG Xinting, LIU Mengmeng, XU Jianping, et al. Image segmentation and recognition algorithm of greenhouse whitefly and thrip adults for automatic monitoring device [J]. Trans Chin Soc Agric Eng, 2018, 34(1): 164 − 170.
    [10] 陈树越, 吴正林, 朱军, 等. 基于凹点检测的粮仓粘连害虫图像分割算法[J]. 计算机工程, 2018, 44(6): 213 − 218.

    CHEN Shuyue, WU Zhenglin, ZHU Jun, et al. Image segmentation algorithm for grain overlapping pests based on pitting detection [J]. Comput Eng, 2018, 44(6): 213 − 218.
    [11] 李芝茹, 吴晓峰, 李全罡, 等. 自动追踪式可升降太阳能虫害监测装置的设计[J]. 森林工程, 2016, 32(6): 85 − 88.

    LI Zhiru, WU Xiaofeng, LI Quangang, et al. Design of the solar liftable device for pest monitoring with automatic trac-king [J]. For Eng, 2016, 32(6): 85 − 88.
    [12] 张红涛, 胡玉霞, 刘新宇, 等. 农田害虫实时检测装置的设计与实现[J]. 河南农业科学, 2007(12): 63 − 65.

    ZHANG Hongtao, HU Yuxia, LIU Xinyu, et al. Design and implementation of real-time detection device for farmland pests [J]. J Henan Agric Sci, 2007(12): 63 − 65.
    [13] 张昊辰, 申小艳, 张欣铃. 智能粮仓害虫监测系统的设计[J]. 通信电源技术, 2013, 30(4): 111 − 112, 126.

    ZHANG Haochen, SHEN Xiaoyan, ZHANG Xinling. Design of smart granary pests monitoring system [J]. Telecom Power Technol, 2013, 30(4): 111 − 112, 126.
    [14] CHENG Mingming, ZHANG Guoxin, MITRA N J, et al. Global contrast based salient region detection[C]// IEEE. International Conferene on Computer Vision and Pattern Recognition. Colorado: IEEE, 2011: 409 − 416.
    [15] 陈宁宁. 几种图像阈值分割算法的实现与比较[J]. 电脑知识与技术, 2011, 7(13): 3109 − 3111.

    CHEN Ningning. Achieve and comparison of image segmentation thresholding method [J]. Comput Knowl Technol, 2011, 7(13): 3109 − 3111.
    [16] 白雪冰, 许景涛, 郭景秋, 等. 基于局部二值拟合模型的板材表面节子与虫眼的图像分割[J]. 浙江农林大学学报, 2016, 33(2): 306 − 314.

    BAI Xuebing, XU Jingtao, GUO Jingqiu, et al. Segmentation of wood surface knots and wormholes based on an improvedLBF Model [J]. J Zhejiang A&F Univ, 2016, 33(2): 306 − 314.
    [17] BORJI J, CHENG Mingming, JIANG Huaizu, et al. Salient object detection: a benchmark [J]. IEEE Trans Image Process, 2012, 24(12): 5706 − 5722.
  • [1] 杨凡, 杨博凯, 李荣荣.  基于图像分割和深度学习的人造板表面缺陷检测 . 浙江农林大学学报, 2024, 41(1): 176-182. doi: 10.11833/j.issn.2095-0756.20230280
    [2] 程昱之, 钟丽辉, 何鑫, 王远, 李朝岚.  基于改进K-means聚类与分水岭的木材横截面管孔分割 . 浙江农林大学学报, 2022, 39(1): 173-179. doi: 10.11833/j.issn.2095-0756.20210219
    [3] 明浩, 苏喜友.  利用特征分割和病斑增强的杨树叶部病害识别 . 浙江农林大学学报, 2020, 37(6): 1159-1166. doi: 10.11833/j.issn.2095-0756.20190752
    [4] 叶晓明, 钱宇汀, 叶雯, 沈黄莹, 曾燕如, 喻卫武, 戴文圣.  香榧绿藻的生物学特性及物种鉴定 . 浙江农林大学学报, 2019, 36(4): 629-637. doi: 10.11833/j.issn.2095-0756.2019.04.001
    [5] 李洪滨, 朱诚棋, 周湘, 马良进, 苏秀.  红哺鸡竹异香柱菌的形态学和分子鉴定 . 浙江农林大学学报, 2016, 33(6): 1040-1044. doi: 10.11833/j.issn.2095-0756.2016.06.016
    [6] 李娟, 张小斌, 薛皎, 黄俊浩, 吴鸿.  基于Lucid多途径检索的竹林害虫智能诊断系统研建 . 浙江农林大学学报, 2016, 33(1): 122-129. doi: 10.11833/j.issn.2095-0756.2016.01.017
    [7] 詹伟君, 罗优波, 张立斌.  松墨天牛优良球孢白僵菌菌株生物学及毒力研究 . 浙江农林大学学报, 2013, 30(1): 83-89. doi: 10.11833/j.issn.2095-0756.2013.01.012
    [8] 梁香媚.  室带槌腹叶蜂的形态和生物学 . 浙江农林大学学报, 2012, 29(3): 431-434. doi: 10.11833/j.issn.2095-0756.2012.03.017
    [9] 闫争亮, 刘云彩, 胡光辉, 刘永刚, 孟梦.  思茅松林的树木多样性与主要害虫发生的关系 . 浙江农林大学学报, 2012, 29(2): 226-231. doi: 10.11833/j.issn.2095-0756.2012.02.012
    [10] 华正媛, 王井田, 刘剑, 王浩杰, 舒金平, 徐天森.  衢州市油茶害虫及天敌种类调查 . 浙江农林大学学报, 2012, 29(2): 232-243. doi: 10.11833/j.issn.2095-0756.2012.02.013
    [11] 张廷廷, 徐华潮, 江挺.  嘉兴市南湖区林业害虫调查及主要害虫风险分析 . 浙江农林大学学报, 2012, 29(4): 621-625. doi: 10.11833/j.issn.2095-0756.2012.04.022
    [12] 信晓颖, 江洪, 周国模, 余树全, 王永和.  加拿大森林火险气候指数系统(FWI)的原理及应用 . 浙江农林大学学报, 2011, 28(2): 314-318. doi: 10.11833/j.issn.2095-0756.2011.02.023
    [13] 石富超, 朱小兵, 龚宁, 王涛.  几种药剂对4种林业害虫的防治效果比较 . 浙江农林大学学报, 2009, 26(3): 358-362.
    [14] 王记祥, 马良进.  虫生真菌在农林害虫生物防治中的应用 . 浙江农林大学学报, 2009, 26(2): 286-291.
    [15] 沈幼莲, 劳冲, 冯林国.  慈溪市蔗扁蛾生物学特性及防治 . 浙江农林大学学报, 2008, 25(3): 367-370.
    [16] 和秋菊, 易传辉.  曲足侎缘蝽生物学特性及防治 . 浙江农林大学学报, 2008, 25(3): 363-366.
    [17] 楼雄伟, 方陆明, 徐爱俊, 叶祥根.  多源数据融合下的森林防火本底数据库建立 . 浙江农林大学学报, 2008, 25(4): 405-410.
    [18] 胡国良, 程益鹏, 楼君芳, 王国杜, 夏国华, 陈建刚, 章祖平.  山核桃花蕾蛆生物学特性及防治技术 . 浙江农林大学学报, 2007, 24(4): 463-467.
    [19] 王义平, 于振东, 吴鸿.  林木昆虫演变为重大害虫的主要环境因子 . 浙江农林大学学报, 2007, 24(6): 752-757.
    [20] 李克恩.  森得保和苏云金杆菌防治竹织叶野螟药效对比试验 . 浙江农林大学学报, 2006, 23(4): 445-448.
  • 加载中
  • 链接本文:

    https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.20200318

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2021/2/369

图(7) / 表(1)
计量
  • 文章访问数:  865
  • HTML全文浏览量:  168
  • PDF下载量:  39
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-05-13
  • 修回日期:  2020-11-21
  • 网络出版日期:  2021-04-01
  • 刊出日期:  2021-04-01

1种基于全局对比度的害虫收集方法

doi: 10.11833/j.issn.2095-0756.20200318
    基金项目:  浙江省基础公益研究计划项目(LGN19C140006);浙江省自然科学基金资助项目(LY16C160007);绍兴市科技计划项目(2018C20013)
    作者简介:

    柳懿祥(ORCID: 0000-0002-7561-7334),从事图像处理与模式识别研究。E-mail: 616731542@qq.com

    通信作者: 汪杭军(ORCID: 0000-0002-7698-1226),教授,博士,从事林业图像处理与模式识别研究。E-mail: whj@zafu.edu.cn
  • 中图分类号: S763

摘要:   目的  目前利用测报灯,通过灯光诱捕昆虫,并由计算机完成昆虫图像的采集、计数和识别已逐步成为害虫测报的重要方法。为了减少昆虫在采样盘上重叠造成的计数和识别误差,基于害虫图像,根据昆虫密度研究采样盘中昆虫的收集方法,从而提高采集效率和精度。  方法  根据昆虫在采样盘上姿态特点,提出基于全局对比度的图像分割方法,结合阈值迭代分割获得昆虫区域,计算昆虫比例,并控制采样盘翻转完成对昆虫的收集。  结果  通过对5种害虫的实际图像进行的试验表明:与水平集、大津法(OTSU)、阈值迭代法和基于直方图对比度的显著性检测(HC)4种算法相比,本研究方法在准确率和召回率上均提高10%以上,取得了较好的结果;同时,在分割速度上比水平集快3倍,与阈值和HC算法基本持平。  结论  基于全局对比度的分割方法简单、高效,在害虫自动测报中具有较高的实际应用价值。图7表1参17

English Abstract

方庆, 谭菊荣, 许惠春, 等. 珍稀濒危植物细果秤锤树群落物种组成与生态位分析[J]. 浙江农林大学学报, 2022, 39(5): 931-939. DOI: 10.11833/j.issn.2095-0756.20220197
引用本文: 柳懿祥, 汪杭军, 徐铁平. 1种基于全局对比度的害虫收集方法[J]. 浙江农林大学学报, 2021, 38(2): 369-376. DOI: 10.11833/j.issn.2095-0756.20200318
FANG Qing, TAN Jurong, XU Huichun, et al. Species composition and niche of Sinojackia microcarpa, a rare and endangered plant[J]. Journal of Zhejiang A&F University, 2022, 39(5): 931-939. DOI: 10.11833/j.issn.2095-0756.20220197
Citation: LIU Yixiang, WANG Hangjun, XU Tieping. A pest collection method based on global contrast[J]. Journal of Zhejiang A&F University, 2021, 38(2): 369-376. DOI: 10.11833/j.issn.2095-0756.20200318
  • 现代农业和林业害虫给粮食安全带来严重威胁,每年产生巨大的经济损失[1]。准确地监测害虫数量变化,预测害虫爆发趋势,可为害虫管理行动提供可靠依据和正确管理方法[2-3]。近年来,出现了许多用于虫害测报的新型害虫诱杀设备及监测系统,这些设备配有多种传感器,可以上传数据到用户手机端,且可以定时对诱杀的害虫拍摄然后识别计数[4-5],节省了大量人力。害虫诱杀设备及监测系统进行虫害测报的重要前提是如何准确分割害虫,针对此问题,许多国内外学者对害虫分割方法进行了研究。SOLIS-SÁNCHEZ等[6]利用目标的几何形态特征(偏心率、面积等)从诱虫板上分割粉虱Aleyrodidae。NING等[7]提出交互式分割方法从复杂背景中分割害虫。吕金娜[8]提出了基于LAB颜色空间的棉花Gossypium spp.害虫普适K聚类害虫图像分割方案,对典型棉花害虫图像进行分割。杨信廷等[9]提出了基于Prewitt、Canny边缘检测算子分割和支持向量机(SVM)的温室粉虱和蓟马Thripidae诱虫板的图像识别算法。陈树越等[10]提出改进的凹点检测和精确分割点定位的方法,实现对黏连害虫的分割。然而,上述方法大多用于单个害虫情况,对于野外诱虫装置获得图像、害虫产生轻微堆叠及颜色不一等问题还有待解决。害虫的分割效果与设备的收集功能密不可分。李芝茹等[11]采用了追踪式太阳能监测装置,在稳定供能的同时也可升降采集不同高度的害虫;张红涛等[12]、张昊辰等[13]针对捕虫诱集部分,设计了多种采集方法以及通道,来提升害虫图像质量。本研究针对害虫收集装置中的采集害虫图像分割精度的要求,根据害虫面积与采样盘面积的比例进行智能翻转,针对翻转功能需要克服的实际拍摄图像光照不均匀,阴影干扰大,及害虫种类颜色繁多等问题,提出了基于全局对比度的图像分割方法,对装置中的实际图像进行分割处理,使其满足采样盘的智能翻转要求。

    • 害虫的智能收集通过害虫面积占比控制采样盘翻转完成,为了准确获取此比例,需要对获取的害虫图像有较高精度的分割。本研究结合基于直方图对比度的显著性检测(HC)算法与阈值迭代算法,提出了满足害虫图像精度的分割方法。整体方法流程如图1所示。

      图  1  方法流程图

      Figure 1.  Algorithm flow chart

      为了减少多种颜色害虫同时出现的频率,首先将害虫图像分成150个同等大小的区域,每个区域大小约1只害虫的面积。在分割区域操作的同时,计算颜色的出现频率,高频颜色替换低频颜色,保证算法运行速度。在此基础上通过HC算法检测图像中的浅色害虫,结合阈值迭代算法分割背景颜色。得到的图像浅色害虫轮廓清晰,并避免了阴影的影响,但图像中深色害虫均未被分割,故将此前检测出的浅色害虫图像中的像素点更新为黑色,加强深色害虫与背景颜色的对比度。再根据直方图计算显著值完成对所有目标害虫的检测。显著值的计算公式为:

      $$S({{\mathit{\boldsymbol{a}}}}) = \frac{{{{\mathit{\boldsymbol{a}}}} - {a_{\min }}}}{{{a_{\max }} - {a_{\min }}}}\text{。}$$ (1)

      式(1)中:a为每个像素的颜色值构成的矩阵;amaxa矩阵中所有值的最大值;amina矩阵中所有值的最小值;S为显著值。

      最后通过sigmoid函数对结果的图像进行二值化计算得到最终的害虫分割图像,用此图像可对害虫面积比例进行计算,完成采样盘翻转的判定。公式为:

      $$\sigma (S) = \frac{1}{{1 + {{\rm{e}}^{ - S}}}}\text{。}$$ (2)

      式(2)中:S为显著值;$\sigma $为sigmoid函数。改进后的方法简易且运行快速,效果显著,可以更好地分割背景和害虫,同时也能分割出深色害虫清晰的轮廓。

    • 在显著性检测算法中,HC算法运行速度较快,且可以保留颜色特征、有效地体现不同颜色的边缘纹理。HC算法根据图像的颜色直方图分布,对颜色级数进行量化和平滑操作,最终计算每个像素的颜色与其他像素颜色的对比度定义为显著值[14]。其计算公式为:

      $$S({I_k}) = \sum\limits_{i = 0}^n {D({I_k},{I_i})} \text{。}$$ (3)

      式(3)中:Ikk点图像颜色转化为Lab颜色空间后的当前像素;Ii则是其余i点的每个像素;n是其余像素的个数;D为2个像素的欧式距离度量;S为显著值。

      在害虫图像中,会出现图像颜色过多影响处理速度的问题。为了解决这个问题,采用了量化图像颜色数据降低颜色数量的方法。对害虫图像压缩,将具有相同颜色像素的组合在一起。公式为:

      $$S({I_k}) = \sum\limits_{j = 0}^n {{f_j}D({c_i},{c_j})} \text{。}$$ (4)

      式(4)中:cii点当前像素Ik的颜色值;cj为其余j点的颜色值;n为颜色不同的像素数量;fjj点像素颜色出现的频率;D为颜色值间的欧式距离度量;S为显著值。

      由于自然图像中的颜色种类只占全颜色空间的小部分,采取忽略出现频率低的颜色值来进一步缩减颜色的数量,删除的颜色值用保留下来的相近颜色值所代替。为了防止缩减颜色数量导致害虫的轮廓细节被忽略,需保证选择留下的高频颜色在图像中覆盖的像素不少于95%。由于目的是完成对害虫轮廓的分割,从而计算害虫的面积,故此操作既不影响分割的精度,还成功提高了算法的运行效率。

      由于在缩减颜色数量时,相近的颜色也许会被量化为不同的值,为了减少这类颜色对检测害虫区域时产生的噪声,使用相近颜色的显著值加权平均值代替原本的显著值。其计算公式为:

      $$S'(c) = \frac{1}{{(m - 1)T}}\sum\limits_{i = 1}^m {[T - D(c,{c_i})]S({c_i})} \text{。}$$ (5)

      式(5)中:S′为平滑颜色后的显著值;S为原本的显著值;T为要改善颜色cm个最近的颜色ci之间的距离之和;D为要改善颜色c与最接近颜色ci的欧式距离度量。

      接着获取处理后图像的颜色直方图,计算1个像素与其余像素的欧式距离,求出所有距离之和,将此定义为显著值。最后将像素的颜色值更新为求出的显著值即处理完毕。

      同时存在多种颜色害虫时,单一使用HC算法存在无法提取深浅色害虫的缺陷。原因是图像中害虫的种类以及颜色多样,害虫间的对比度过大,且存在一定的阴影及轻微重叠干扰。故本研究决定结合阈值迭代对害虫图像进行处理,以此优化前景与背景之间的颜色差。

    • 利用阈值迭代分割算法解决了图像各区域亮度不同等因素,可更好地分割图像背景[15-16]。该算法可在不同区域采用不同的阈值进行分割,适合处理多种颜色害虫混合的图像,解决了害虫阴影的干扰问题,完成了分割单一颜色背景的目的。

      该算法主要通过迭代方式逼近最佳阈值,从而达到分割背景的目的,关键步骤在于找到最佳阈值。先将整幅图像分成多张互相之间有50%重叠的子图像,分别作出这些子图像的直方图。对这些子图像的直方图进行检测,查看是否为双峰,如果是双峰,则进行阈值迭代法确定阈值,如果不是双峰,就不进行处理。接着确定阈值,以此为依据划分前景和背景。然后分出阈值的2边数据为A和B部分。分别计算A和B部分的均值,然后再总的求A和B部分的均值即为新的阈值。设定迭代次数进行不断寻找最佳阈值,用得到的阈值进行图像的二值化。公式如下:

      $${R_0} = \frac{{\displaystyle\sum\limits_{R(i,j) = 0}^{{T_{\rm{K}}}} {R(i,j) \times N(i,j)} }}{{\displaystyle\sum\limits_{R(i.j) = 0}^{{T_{\rm{K}}}} {N(i,j)} }}\text{;}$$ (6)
      $${R_{\rm{G}}} = \frac{{\displaystyle\sum\limits_{R(i,j) = {T_{\rm{K}}}}^m {R(i,j) \times N(i.j)} }}{{\displaystyle\sum\limits_{R(i,j) = {T_{\rm{K}}}}^m {N(i,j)} }}\text{。}$$ (7)

      式(6)~(7)中:R(i, j)为图像(i, j)像素点的灰度值;N(i, j)为(i, j)像素点的权重系数;TK为确定阈值;m为像素个数;R0RG分别为目标和背景的平均灰度值。

    • 图像获取的步骤分为害虫的诱杀收集和拍摄2步。害虫的诱杀通过打开图2的黑光灯引诱害虫,害虫飞至黑光灯边上的电网时,被电死掉落到下方漏斗状装置内部,然后滑落到采样盘上。摄像头在LED灯点亮后进行拍摄,采集图像大小为2592像素×1944像素。最后将拍摄的图像传送至控制系统内部进行图像处理。

      图  2  害虫收集装置示意图

      Figure 2.  Schematic diagram of pest collection device

    • 本研究设备为Intel(R) Core(TM)i7-8750H CPU、16G 64位PC机和树莓派官方500万像素摄像头,系统和软件环境为Window10,Jupyter Notebook,Raspberry Pi3 B+。在实际装置托盘上方20 cm处使用装置内部的树莓派摄像头进行拍摄。研究样本为5种害虫在白底托盘和红底托盘上的实际图像。为了验证本研究的可靠性与有效性,共选用100张实际拍摄图像后运行算法,取平均值作为最后研究结果数据。因实际装置的硬件性能所限制,且实际装置对智能翻转的图像处理有一定的速度和精度需求,故排除与深度学习分割方法的对比。本研究选取了4种综合处理效果较好的经典算法作为比较。

      使用不同算法对实际摄像头拍摄的图片白底样本进行分割。由图3可见:本研究算法与实际比率的接近程度明显高于其他算法。且在白底托盘中,本研究算法可以更好避免阴影的影响,分割出目标害虫的轮廓。比较图4图5可知:大津算法(OTSU)出现了错误分割阴影的结果,它将较多的阴影区域分割,将会导致害虫比例严重误判。而本研究算法较好改善了该问题,至于仅剩的阴影噪声问题将通过改善拍摄条件弥补。

      图  3  白底样本分割效果对比

      Figure 3.  Comparison of sample segmentation effects on white background

      图  4  OTSU算法与本研究算法对阴影分割的比较

      Figure 4.  Comparison of shadow segmentation between OTSU and research algorithm

      图  5  HC算法与本研究算法对害虫细节分割的比较

      Figure 5.  Comparison of HC algorithm and research algorithm on pest detail segmentation

      白底样本的分割结果(图5)进一步表明:HC算法无法识别害虫颜色与背景颜色对比度不高的害虫,而本研究算法在改进其算法后,可以较清晰地分割出这些原本分割效果不佳的害虫。

      综合了白底样本采样盘易出现灯光照射以及阴影的干扰,且由于浅色害虫与白色采样盘背景对比度过于接近,更加不利于对害虫的分割,故改进装置的采样盘为红色。拍摄深色害虫与浅色害虫同时存在于采样盘上时的图像,作为红底样本,并使用不同算法对其进行分割处理。从图6图7可见:其他算法主要完成了浅色害虫的分割,但均无法分割出深色害虫并计算害虫面积。本研究算法则完成了深浅害虫的同时分割,从图7可清晰观察到深色害虫的翅膀及触角等细节,证明本研究算法在多种颜色害虫存在时分割结果更精确可靠。

      图  6  红底样本分割效果对比

      Figure 6.  Comparison of sample segmentation effects on red background

      图  7  HC算法与本研究算法对深色害虫分割的比较

      Figure 7.  Comparison of HC algorithm and research algorithm on dark pest segmentation

    • 为了评价算法的分割图像效果优劣,本研究将采用平均分割时间、准确率和召回率对算法的分割结果进行衡量[17]。为了比较各算法的效率,将使其对每幅图像进行多次分割,然后取平均处理时间来作最终评价。准确率(P)和召回率(R)的公式如下:

      $$P = \frac{{\displaystyle\sum\limits_{n = 1}^N {{G_n} \times {S_n}} }}{{\displaystyle\sum\limits_{n = 1}^N {{S_n}} }}\text{;}$$ (8)
      $$R = \frac{{\displaystyle\sum\limits_{n = 1}^N {{G_n} \times {S_n}} }}{{\displaystyle\sum\limits_{n = 1}^N {{G_n}} }}\text{。}$$ (9)

      式(8)~(9)中:N为像素个数;Gn表示基准图像第n个像素是否为分割目标像素;Sn表示分割图像第n个像素是否为分割目标像素,值为0或1。本研究取权重系数为0.3,来防止由于显著性检测导致召回率过高的问题。准确率表示目标分割的准确性,召回率为检测区域与基准图像区域的比值,表示算法分割的完整性。使用选取的5种算法对样本进行图像处理,并且将人工分割的图像作为基准图像,得到3个评价指标(表1)。

      表 1  准确率、召回率和分割速度对比

      Table 1.  Comparison of accuracy recall, and segmentation speed

      分割算法准确率召回率分割速度
      白底样本红底样本白底样本红底样本白底样本红底样本
      水平集算法0.720.740.620.641.684.35
      OTSU算法0.700.770.630.700.360.84
      阈值迭代算法0.830.830.810.690.560.87
      HC算法0.860.880.860.630.430.98
      本研究算法0.920.950.950.890.661.05

      表1可知:本研究算法的准确率及召回率最高,其他算法易出现分割阴影等非目标轮廓,或分割不出深色害虫。准确率相比其他4种算法均提高约13%。在召回率上,本研究算法提高了10%以上,最高达53%。准确率和召回率的提升,体现了本研究算法的精确性和可靠性。在分割速度方面,本研究算法高出水平集算法约3倍,且总耗时满足装置的实际需求。除水平集算法外,其他3种算法虽然分割时间短,但无法分割多种害虫同时存在的图像,故可以接受本研究算法小幅延长分割时间后,完成了更高要求的分割。

    • 本研究算法处理后得到的结果误差控制在5%以内。即使在光照不稳定和害虫种类以及姿态复杂的环境下采集的图像中,本研究算法也可以较精确避免轻微重叠和阴影的影响进行分割。在多种颜色害虫同时存在时,本研究算法也可分割出几乎接近实际害虫的轮廓。因此,本研究算法简单、高效且充分满足分割精度和速度的需求,适宜于害虫收集装置内部使用,达到对害虫比例进行准确计算的目的。

参考文献 (17)

目录

/

返回文章
返回