CN102509017A - 一种用计算机预测刨花板强度的方法 - Google Patents
一种用计算机预测刨花板强度的方法 Download PDFInfo
- Publication number
- CN102509017A CN102509017A CN2011103552931A CN201110355293A CN102509017A CN 102509017 A CN102509017 A CN 102509017A CN 2011103552931 A CN2011103552931 A CN 2011103552931A CN 201110355293 A CN201110355293 A CN 201110355293A CN 102509017 A CN102509017 A CN 102509017A
- Authority
- CN
- China
- Prior art keywords
- particieboard
- representes
- point
- shaving boards
- intensity
- 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.)
- Pending
Links
Images
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开了一种用计算机预测刨花板强度的方法。现有技术估计刨花板强度的方法是由经验丰富的工人通过观察板材表面碎木块的走向对板材进行检测,但这是个繁琐而又主观的过程。本发明通过对刨花板表面图像进行分析估计其静力弯曲强度,提供了一种计算机辅助的刨花板强度估计方法。本发明的方法由三个有效步骤组成:首先对刨花板的表面图像进行边缘检测,得到刨花碎料的边缘;其次对边缘图像进行霍夫线检测,得到直线边缘;最后对刨花板表面直线的长度和走向进行分析,从而估计刨花板的强度。本发明与通常的人工检测相比,能有效地提高生产效率,并提高对刨花板强度估计的精度,有效地降低了人工估计的主观不稳定性,同时大大减少了人工工作量。
Description
技术领域
本发明属于数字图像处理领域,涉及一种用计算机预测刨花板强度的方法。
背景技术
数字图像处理起源于二十世纪二十年代,经过近一个世纪的发展,目前已经成为信息科学、生物学、医学等诸多领域的基础研究对象。随着信息时代的到来,数字图像处理在计算机视觉、机器学习、人工智能等领域都得到广泛的应用,其重要性日益凸显。
刨花板(article board)是以木材、木质碎料、木质纤维或其他植物纤维为原料,加胶粘剂和其他添加剂制成的板材的统称。刨花板比天然木材的稳定性好,广泛用于建筑、车船、家具、包装等方面。
由于刨花板的制作过程采用刨花碎料进行压制合成,因此刨花板的力学性能与刨花碎料的形状和分布密切相关,尤其是静力弯曲强度,一般情况下形状细长且大多沿水平分布的碎料组成的刨花板强度较大。通常估计刨花板强度的方法是由经验丰富的工人通过观察板材表面碎木块的走向对板材进行检测,但这是个繁琐而又主观的过程。
发明内容
本发明针对现有技术的不足,提供了一种用计算机通过数字图像处理技术自动分析刨花板表面碎木块的长度和基本纹路朝向,以此估计出刨花板强度的方法。该方法包括以下步骤:
步骤一、对刨花板表面图像进行边缘检测
首先使用高斯卷积模板对图像进行平滑滤波,有效去除图像中的噪声,减少图像噪声对边缘检测的影响。
高斯卷积模板各个位置的权重与相应位置到模板中心点的距离满足如下关系:
式(1)中,r表示当前点到模板中心的距离,f表示当前点的权值,σ表示高斯模板的标准方差,e表示自然对数;
使用Canny边缘检测算法检测出边缘图像,该边缘检测算法计算图像中每一点的梯度幅值及角度,然后使用非极大值抑制的方法找出梯度幅值的局部最大值点作为边缘的候选点,最后使用边缘跟踪算法将边缘线勾勒出来;
步骤二、对检测出来的边缘图像进行霍夫线检测,图像中的直线满足如下关系:
xcosθ+ysinθ=ρ (2)
式(2)中,ρ表示直线与原点的距离,θ表示直线的倾斜角,x表示像素点的水平坐标,y表示像素点的垂直坐标;
检测的结果是得到关于ρ和θ的一个二维矩阵,矩阵中每一个元素存储了属于某一特定ρ和θ的直线的点序列;
步骤三、对检测出来的线段进行分析,估计刨花板强度。
刨花板的力学性能尤其是静力弯曲强度与刨花碎料的形状和分布密切相关,通常形状细长且大多沿水平分布的碎料组成的刨花板强度较大。基于以上规律,刨花板的强度与表面直线的长度和水平分布成正比,三者的关系如下所示:
式(3)中,F表示刨花板的估计强度,l表示线段的长度,θ表示线段与水平方向的夹角,θlow表示对刨花板的强度有贡献的角度阈值。
有益的效果:
本发明与通常的人工检测相比,能有效地提高生产效率,并提高对刨花板强度估计的精度,有效地降低了人工估计的主观不稳定性,同时大大减少了人工工作量。
(1)精确性。本发明方法能够对产品表面的碎料纹路做全面的分析,客观地根据纹路的统计情况进行强度估计,而人工估计则会因人因地因时而异,因此,使用计算机方法能够提高估计精度。
(2)高效性。本发明方法可以对刨花板强度进行实时的估计,因此可以将其加入生产流水线,对生产出的刨花板按强度自动分类。大大降低人工成本。
附图说明
图1是以100倍放大了的9乘9的高斯滤波模板;
图2直观地描述了某一边缘上各点梯度的方向;
图3是某一幅刨花板表面图像的梯度幅值图;
图4描述了梯度方向与梯度幅值比较点的映射关系;
图5显示了霍夫线检测完并经过断线后产生的线段图。
具体实施方式
一种用计算机预测刨花板强度的方法具体流程如下:
1)对刨花板原始图像进行边缘检测
使用高斯卷积模板对图像进行平滑滤波,然后使用Canny边缘检测算法检测出图像的边缘信息;
第一步,使用高斯卷积模板对图像进行滤波,基于的高斯分布的模板产生函数为:
其中,f(r)为离模板中心点距离为r的点的权值,σ为高斯模板的标准方差;
使用该高斯函数产生一个9乘9的滤波模板,使用时,由于做了归一化处理,模板中所有的值相加为1,如图1所示,该模板中数值放大100倍的结果,使用该滤波模板与图像进行卷积,能够减少图像中的噪声,其原理是,图像中噪声点往往与周围的点差异很大,与高斯滤波模板做卷积能够以将周围的像素以一定的权重加到操作点上,这样,与周围像素点差异过大的噪声点就会被掩盖;
第二步,使用Canny边缘检测算法检测出图像的边缘信息,该边缘检测算法计算图像中每一点的梯度幅值及其角度,某个点梯度的幅值表示该点在某一方向上的最大变化值,这个方向就是梯度的角度方向,可以用图2直观地表示某一边缘上各点的梯度方向,梯度幅值和角度的计算公式为:
θx,y=arctan(Gy/Gx)
其中,Gx为某点水平方向上的梯度分量,Gy为某点垂直方向上的梯度分量,G为梯度的幅值,θxy代表梯度方向的角度。
如图3所示,用灰度图的形式直观地表示出一幅刨花板表面图像的梯度幅值,计算出梯度幅值及角度后,使用非极大值抑制的方法找出梯度幅值的局部最大值点,即要判断一点是否为梯度幅值最大值点,只要判断其梯度方向上的相邻两个点的梯度幅值是否比自身的小,若是,则该点是梯度幅值的局部最大值点,这样点作为边缘的候选点,图4画出了用于角度映射的图,因为图像是一个离散函数,其相邻的只有8个点,而梯度的方向角却是一个连续的值,需要将角度映射成相应的用于比较的点的位置,每两条相邻的虚线间的连续角度映射成两虚线间实线对应方向上的两个点的位置。候选点中有一部分梯度值特别大的点一开始就被标记成边缘点,这里涉及到一个阈值,即有多少比例的高梯度值点成为边缘点。最后使用边缘跟踪算法将边缘线勾勒出来,之所以需要进行边缘跟踪,是因为边缘上的某些点梯度值可能并不大,不一定能够在一开始就被标记成边缘点,如果不进行连接,将造成边缘的不连续性,边缘跟踪是一个递归的深度优先的计算过程,以梯度值特别大的那部分点作为跟踪起点,遍历与其连接的边缘候选点,并标记为边缘点,该过程以找不到相邻的候选点为终止条件。
2)对检测出来的边缘图像进行霍夫线检测,霍夫线检测算法使用的直线标准式为:
xcosθ+ysinθ=ρ
其中,ρ为直线与原点的距离,θ为直线的倾斜角。
使用边缘点匹配上述直线方程,这个过程对每一个边缘点(x,y)计算其在0至360度时符合的直线方程,即使用(x,y,θ)求出ρ,将其结果记录进矩阵,结果是得到关于ρ和θ的一个二维矩阵,矩阵中每一个元素存储了属于某一特定ρ和θ的直线的点序列。
3)对检测出来的线段进行统计
第一步,断线,由于霍夫线检测输出的矩阵元素是属于某一直线的点集,这些点集需要通过分割,产生一条或多条连续的线段,以便进行线条的有效统计,断线过程是扫描每一个直线点集,对于每一个点集,按下标顺序检查,如果一个点与其邻近点的距离大于某个阈值,就认为这根线段断开了,并记录找到的长度大于某阈值的有效线段,图5显示了霍夫线检测完并经过断线后产生的线段图。
第二步,计算木质人造板强度的预测值,由于刨花板的力学性能尤其是静力弯曲强度与刨花碎料的形状和分布密切相关,通常形状细长且大多沿水平分布的碎料组成的刨花板强度较大。基于以上规律,刨花板的强度与表面直线的长度和水平分布成正比,三者的关系如下所示:
其中,F表示刨花板的估计强度,l表示线段的长度,θ表示线段与水平方向的夹角,θlow表示对刨花板的强度有贡献的角度阈值。
公式的分子部分表示碎料边缘在水平方向上的分量,即对强度的贡献,分母在分子的基础上添加了一个直线边缘在垂直方向上的分量,二者的比值反映了刨花碎料的分布情况,在水平方向上分布越多,这个比值就越大。然后再乘上线段的长度,表示越长的碎料组成的刨花板强度越大。因此F值越大,则刨花板的强度越大。这个公式反映了“形状细长且大多沿水平分布的碎料组成的刨花板强度较大”的规律。
Claims (1)
1.一种用计算机预测刨花板强度的方法,其特征在于该方法包括以下步骤:
步骤一、对刨花板表面图像进行边缘检测;
首先使用高斯卷积模板对图像进行平滑滤波,有效去除图像中的噪声,减少图像噪声对边缘检测的影响;
高斯卷积模板各个位置的权重与相应位置到模板中心点的距离满足如下关系:
式(1)中,r表示当前点到模板中心的距离,f表示当前点的权值,σ表示高斯模板的标准方差,e表示自然对数,
使用Canny边缘检测算法检测出边缘图像,该边缘检测算法计算图像中每一点的梯度幅值及角度,然后使用非极大值抑制的方法找出梯度幅值的局部最大值点作为边缘的候选点,最后使用边缘跟踪算法将边缘线勾勒出来;
步骤二、对检测出来的边缘图像进行霍夫线检测;
图像中的直线满足如下关系:
xcosθ+ysinθ=ρ (2)
式(2)中,ρ表示直线与原点的距离,θ表示直线的倾斜角,x表示像素点的水平坐标,y表示像素点的垂直坐标,
检测的结果是得到关于ρ和θ的一个二维矩阵,矩阵中每一个元素存储了属于某一特定ρ和θ的直线的点序列;
步骤三、对检测出来的线段进行分析,估计刨花板强度;
刨花板的力学性能尤其是静力弯曲强度与刨花碎料的形状和分布密切相关,通常形状细长且大多沿水平分布的碎料组成的刨花板强度较大;基于以上规律,刨花板的强度与表面直线的长度和水平分布成正比,三者的关系如下所示:
式(3)中,F表示刨花板的估计强度,l表示线段的长度,θ表示线段与水平方向的夹角,θlow表示对刨花板的强度有贡献的角度阈值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011103552931A CN102509017A (zh) | 2011-11-10 | 2011-11-10 | 一种用计算机预测刨花板强度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011103552931A CN102509017A (zh) | 2011-11-10 | 2011-11-10 | 一种用计算机预测刨花板强度的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102509017A true CN102509017A (zh) | 2012-06-20 |
Family
ID=46221102
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011103552931A Pending CN102509017A (zh) | 2011-11-10 | 2011-11-10 | 一种用计算机预测刨花板强度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102509017A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107710229A (zh) * | 2016-10-20 | 2018-02-16 | 深圳配天智能技术研究院有限公司 | 图像中的形状识别方法、装置、设备及计算机存储介质 |
CN116883993A (zh) * | 2023-09-06 | 2023-10-13 | 临沂大学 | 基于视觉的玫瑰花茶干花分选方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5011886A (en) * | 1983-12-21 | 1991-04-30 | Rwe-Dea Aktiengesellschaft Fur Mineraloel Und Chemie | Process for producing modified phenolic resin bonding agents and use thereof for producing particle board |
CN101013126A (zh) * | 2007-01-15 | 2007-08-08 | 浙江大学 | 用计算机预测木质人造板强度的方法 |
-
2011
- 2011-11-10 CN CN2011103552931A patent/CN102509017A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5011886A (en) * | 1983-12-21 | 1991-04-30 | Rwe-Dea Aktiengesellschaft Fur Mineraloel Und Chemie | Process for producing modified phenolic resin bonding agents and use thereof for producing particle board |
CN101013126A (zh) * | 2007-01-15 | 2007-08-08 | 浙江大学 | 用计算机预测木质人造板强度的方法 |
Non-Patent Citations (2)
Title |
---|
JOSEPH KHEDARI等: "New insulating particleboards from durian peel and coconut coir", 《BUILDING AND ENVIRONMENT》 * |
沈哲红等: "如何提高刨花板静曲强度和降低甲醛释放量", 《林产工业》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107710229A (zh) * | 2016-10-20 | 2018-02-16 | 深圳配天智能技术研究院有限公司 | 图像中的形状识别方法、装置、设备及计算机存储介质 |
CN107710229B (zh) * | 2016-10-20 | 2021-02-26 | 深圳配天智能技术研究院有限公司 | 图像中的形状识别方法、装置、设备及计算机存储介质 |
CN116883993A (zh) * | 2023-09-06 | 2023-10-13 | 临沂大学 | 基于视觉的玫瑰花茶干花分选方法 |
CN116883993B (zh) * | 2023-09-06 | 2023-12-01 | 临沂大学 | 基于视觉的玫瑰花茶干花分选方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109272524B (zh) | 一种基于阈值分割的小尺度点云噪声去噪方法 | |
CN106683093B (zh) | 板材外观质量综合量化评价方法 | |
CN103390280A (zh) | 基于灰度-梯度二维对称Tsallis交叉熵的快速阈值分割方法 | |
CN101995412B (zh) | 一种鲁棒的玻璃划伤缺陷检测方法及其装置 | |
CN107656902B (zh) | 不同采样间距下的结构面粗糙度系数统计方法 | |
CN108662983B (zh) | 钢轨磨耗检测计算相关度系数的方法 | |
CN101013126A (zh) | 用计算机预测木质人造板强度的方法 | |
CN109612412A (zh) | 一种预制混凝土构件结合面粗糙度计算方法及评价系统 | |
CN107133565A (zh) | 基于线激光的激光刻型线特征提取方法 | |
CN105571511A (zh) | 一种船舶外板成型精度在线检测方法 | |
CN106482674A (zh) | 基于中智数函数的结构面粗糙度尺寸效应的近似表达方法 | |
CN104463845A (zh) | 一种水火弯板特征配准点的选取方法及系统 | |
CN106340010A (zh) | 一种基于二阶轮廓差分的角点检测方法 | |
CN110263422A (zh) | 基于简单平面滑动法的环形凸坡稳定性评价方法 | |
CN105627935A (zh) | 产品料厚三维分析系统及方法 | |
CN106662438B (zh) | 线位移评价方法、线位移评价装置以及记录介质 | |
CN107036561A (zh) | 基于中智数函数的结构面粗糙度各向异性的近似表达方法 | |
CN102509017A (zh) | 一种用计算机预测刨花板强度的方法 | |
CN104517295A (zh) | 一种三维太赫兹图像的图像分割方法及系统 | |
CN106680299B (zh) | 板材节疤外观质量评价方法 | |
CN204555910U (zh) | 快速平面度检测治具 | |
CN104517296A (zh) | 一种三维毫米波图像的分割方法及系统 | |
CN104180789B (zh) | 基于图形匹配算法的叶片检测方法 | |
CN201081763Y (zh) | 一种纺织品缺陷检测装置 | |
CN102980536B (zh) | 刹车片钢背平面度检测系统及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20120620 |