基于图像处理的甘蔗茎节识别与定位_张卫正

理、Sobel 边 缘 提 取，结 合 数 学 形 态 学 中 的 膨 胀 和 腐 蚀，获得甘蔗 茎 节 的 边 缘 直 线，达 到 识 别 甘 蔗 茎 节 的 目的［11］； 但只定位了茎节的中间位置，没有考虑采集 的甘蔗图像倾斜时定位的范围。
本文提出基于数字图像处理技术的甘蔗茎节识 别定位方法，依据沿包含茎节的甘蔗的长轴方向在灰
1 甘蔗图像采集及茎节形态分析
以浙江省义乌市甘蔗优质生产基地主产的“义红 一号”果蔗为研究对象。试验台采用 OEM 品牌的型 号为 C1745B 的 摄 像 头，该 摄 像 头 采 用 CMOS 传 感 器，可进行 自 动 曝 光 控 制、白 平 衡 和 颜 色 校 正。 考 虑 到相较于其它颜色，黑色背景可以消除甘蔗边缘阴影
研究中，陆尚平等以甘蔗图像 HSV 颜色空间的 S 分 角的茎节快速、精确定位的问题。
量经阈值分割、数学形态滤波处理作为模板，和 H 分 量经阈值分割的反图像进行与运算得到合成图，采用 支持向量机对合成图各列块识别茎节类与茎间类，对 茎节列块聚类识别得到茎节数与茎节位置［10］； 但平均 识别率为 93． 359%，且 采 用 的 算 法 需 要 较 多 的 输 入 量。黄亦其提出一种甘蔗茎节识别方法，通过灰度处
在实际的甘蔗图 像 采 集 过 程 中 ，不 能 保 证 甘 蔗 区 域处于完全水平状态，会导致茎节与垂直方向有一定 的夹角。为了提高茎节的定位精度，此时的定位范围 应如图 1（ b） 两条直线所示。
如图 4（ b） 所示。
第4期
图 2 茎节识别与定位流程图 Fig．2 Flowcharts of internode identification and location
阈值处理效果好于黑色背景，原因是增大了灰度图像
背景与待识别目标的灰度差。因此，将图 1（ a） 中的
黑色背景转为白色，如图 3 所示。
在 MatLab2010b 中，可以使用图像处理工具箱中
的 fgb2gray 函数将 ＲGB 图像转换为灰度图像，则有
Gray（ x，y） = 0．299Ｒ（ x，y） + 0．587G（ x，y） +
图 4 图像预处理 Fig．4 mage preprocessing
中值滤波法对消 除 椒 盐 噪 声 非 常 有 效 ，对 甘 蔗 表 面的白色果粉小斑点也有较好的消除效果。中值滤 波采用 5×5 的模板，滤波程序语句为： K = medfilt2（ I， ［5，5］） 。
自动阈值函数 graythresh 会 根 据 背 景 的 不 同，使 用 Otsu 算法自动调整阈值。Otsu 算法也称最大类间 差法，被认为 是 图 像 分 割 中 阈 值 选 取 的 最 佳 算 法，不 受图像亮度和对比度的影响。设灰度图像灰度级是 L，则灰度范围为［0，L－1］，利用 Ostu 算法计算图像的 最佳阈值为
种自动切 断 过 程 中 防 伤 芽 的 功 能，影 响 了 产 业 的 发 度值上会有明显的上升来识别茎节，并以最大值所在
展。现有的方法为人工切断和种植机的定长切断，效 位置进行茎节两端的定位。为了提高定位精度，分割
率有待提高。计算机图像处理技术广泛的应用在缺 出甘蔗区域，并计算该区域的质心、倾角和等效长短
通 过 regionprops 函 数 测 定 该 区 域 的 质 心 、倾 角 和 等 效 长 短 轴 长 度 ； 以 分 割 后 的 甘 蔗 图 像 对 二 值 化 图 像 掩 膜 ，得 到
含有 茎 节 和 干 扰 信 息 的 图 像 ； 对 该 图 像 进 行 旋 转 ，计 算 每 列 像 素 值 之 和 ，统 计 分 析 最 大 值 所 在 列 ，并 结 合 质 心 、等
程序语句为： stats = regionprops （ TIANC，' Orientation'，' BoundingBox'，'Centroid'，'MajorAxisLength'，'MinorAxisLength'） ； 可以采集到倾角 α 为 stats．Orientation = 11．5579，长 轴等效长度为 stats．MajorAxisLength = 370．090，短轴等 效长度为 stats．MinorAxisLength = 180．1494。 2．3 掩膜与旋转图像 将经孔洞填充后分割出来的甘蔗图像的两边宽 度各减少 5 列，然后对甘蔗二值化图像（ 见图 5） 做掩 膜，即可提取 出 来 甘 蔗 图 像 的 带 有 茎 节 信 息 的 区 域 ， 如图 9（ a） 所示。从图 9（ a） 中也可看出： 除了茎节信 息外，还有其他信息的干扰。
图 1 甘蔗图像及茎节分析 Fig．1 Sugarcane image and internode analyze image
2 图像分析处理
为了快速识别出 茎 节 ，采 用 目 前 常 用 的 灰 度 级 转
换、中值滤波 及 阈 值 处 理 等 对 原 始 图 像 进 行 预 处 理 ，
获取二值图像； 然后分割出甘蔗区域，测定该区域的
·217·
2016 年 4 月
农机化研究
像直 接 存 入 计 算 机。 计 算 机 操 作 系 统 是 Microsoft Windows 7，图像处理软件采用 MatLab2010b。
蔗芽环绕茎节分 布 ，有 时 所 拍 摄 图 像 中 蔗 芽 位 于 茎节的后侧，因 而 在 图 像 中 没 有 显 示，因 此 以 茎 节 的 识别间接得到蔗芽的位置。图 1（ b） 中矩形框所围部 分为茎节，宽度在 7 ～ 16 mm 之间变化； 茎节左端的白 色部分是甘蔗叶与甘蔗茎节的连接部分，宽度在 1．5 ～ 5 mm 。 由于茎节白色 部 分 特 征 相 对 明 显 ，本 文 先 提 取 茎节白色部分然后进行定位，进而将定位范围扩展到 整个茎节。
质心、倾斜角和等效长短轴等相关参数； 利用分割出
的甘蔗区域图像对二值化图像进行掩膜，以倾角反向
旋转掩膜后 的 图 像，使 茎 节 处 于 垂 直 方 向，以 便 于 后
续的茎节特征的计算和分析，流程如图 2 所示。
2．1 图像预处理
经过试验验证，在 MatLab 中的 图 像 自 动 阈 值 处
理函数 graythresh 对白色背景中的甘蔗感兴趣区域的
最佳阈值。用阈值 t 分割图像对原图像 f（ x，y） 使用上
面的准则找到阈值 t，进行二值化。公式为
{0 f（ x，y） ＜ T
g（ x，y） = 1 f（ x，y） ≥ T
（ 3）
二值化的图像如图 5 所示。
相关参数进行提取 （ 如甘蔗区域的质心，相对于水平 的 X 轴 的 倾 角 α ，等 效 的 长 轴 ST、短 轴 UV 的 长 度 等） ，其示意图如图 8 所示。
·218·
2016 年 4 月
农机化研究
第4期
t = max［w0（ t） * （ u0（ t） － u） 2 +
w1（ t） * （ u1（ t） － u） 2］
（ 2）
其中，当分割的阈值为 t 时，w0为背景比例： u0 为背
景均值； w1为前景比例； u1为前景均值； u 为整幅图像的
均值［12］。使以上表达式值最大的 t，即为分割图像的
为了测定图像中 甘 蔗 的 倾 角 ，需 要 确 定 甘 蔗 的 实 际区域，将采集的原始图像的两边宽度各增加 5 列， 并填充黑色使之与背景相同，如图 6 所示。
图 8 甘蔗的长短轴及倾角示意图 Fig．8 Schematic diagram of sugarcane major axis，minor axis and tilt angle
图 6 将图像宽度扩展并填充黑色 Fig．6 Expand the image width and fill with black
陷检测、品质 检 测、水 果 分 级、蔬 菜、农 业 和 食 品 工 程 等领域［4－8］。运 用 计 算 机 图 像 处 理 技 术 识 别 甘 蔗 茎
轴长度； 以分割后的甘蔗图像对二值化图像掩膜，得 到含有茎节和干扰信息的图像； 将该图像以倾角反向
节，防止甘蔗种植机在切割过程中的伤芽，实现甘蔗 旋转，计算每列像素值之和，统计分析最大值所在列，
的精细化种植，具有重要意义。
并结合质心、等效长短轴，得到茎节上下端点坐标，以
Moshashai 等利用灰度图像阈值分割的方法对甘 蔗茎节识别做了初步研究［9］。目前，国内在此领域的
倾角的度数 再 次 进 行 反 向 旋 转，最 终 得 到 茎 节 位 置。 该研究创新性地解决了实际生产过程中具有一定倾
mail） zjqiu@ zju．edu．cn。
的影响，所以在普通光源下采集黑色背景上的甘蔗彩 色图像。拍摄前先将甘蔗叶剥干净，保持镜头主光轴 与工作台垂直。采集的图像如图 1（ a） 所示。
其中，图像大小为 320 × 240 像素，输出的图像格 式是 JPEG。摄像头通过 USB 2．0 接口将采集到的图
以 倾 角 的 度 数 进 行 反 向 旋 转 ，最 终 得 到 茎 节 位 置 。 试 验 结 果 表 明 ： 甘 蔗 茎 节 识
别 与 定 位 方 法 处 理 速 度 快 ，茎 节 识 别 率 高 ，左 右 端 的 定 位 误 差 分 别 小 于 0 ． 9 mm 和 2 ． 4 mm 。
收稿日期： 2015－03－18 基金项目： 国家高技术研究发展计划项目（ 2012AA101903） 作者简介： 张卫正（ 1982 －） ，男，河南安阳人，博士研究生，（ E － mail ）
weizheng008@ 126．com。 通讯作者： 裘正军（ 1971 －） ，男，浙江嵊州人，教授，博士生导 师，（ E －
关键词： 甘蔗茎节； 图像处理； MatLab； 掩膜； 定位