专利名称：一种自动纸浆纤维检测中分离交叉纤维的方法
一种自动纸浆纤维检测中分离交叉纤维的方法 技术领域 本发明涉及一种自动纸浆纤维检测中分离交叉纤维的方法。
背景技术：
在对纸浆纤维进行检测和识别时，纤维样品中经常存在纤维交叉重叠的现象，为 了对交叉纤维进行更好的识别，首先必须对其进行分离。而分离交叉纤维的关键在用于交 叉点的检测。传统的检测算法之一是通过判断图像非零像素8邻域内值的变化次数来判断 该像素点是否为交叉点，但通过对已有算法的实验、研究和分析，发现该方法并不能获得完 整的交叉点的集合，而且会出现直观认识中的一个点对应目标集合中多个相邻点的情况， 而且纤维交叉重叠的随机性比较大，重叠情况复杂，已有算法并不能涵盖大部分情况，准确 性不够，需要改进。另外一种传统的检测方法是首先按照8邻域中拥有一个相邻点即为顶 点的原则找到纤维中的第一个顶点a，再找到纤维中的一个交叉点O，接下来通过判断其他 顶点和交叉点ο的斜率与a和ο的斜率的比较确定它们的配对关系，已达到分离的目的。 但是实验证明，通过判断顶点和交叉点的斜率决定配对关系的方法也具有不稳定性，因为 纤维存在各种形态的弯曲程度，一旦弯曲程度过大就会导致判断失误，如图3所示，直观判 断X1X2应该是一根纤维，而经过判断后却是X1Y2是一根纤维，因此该方法具有严重的局限 性。同时纸浆纤维都具有一定的宽度，一旦纤维交叉在一起，其细化骨架的交叉点及其邻域 情况比较复杂，直接利用8邻域的个数判断是不是交叉点会漏掉许多交叉纤维，而认为是 一根纤维，如图4所示，虚线圆所圈区域为交叉段，但往往会被漏掉，而实线圆所圈区域不 是交叉点或段，却会被判断为交叉点。所以现有的方法在准确判断交叉点或交叉段位置及 并成功分离纤维都存在很大的误差。
发明内容
本发明的目的就是为了解决目前不能准确的分离交叉目标物体，影响对单一目标 物体的测量等操作，误差较大等问题，提供了一种在噪声干扰情况下，准确分离交叉纤维的 自动纸浆纤维检测中分离交叉纤维的方法。为实现上述目的，本发明采用了如下技术方案
一种自动纸浆纤维检测中分离交叉纤维的方法，它通过分析交叉点或交叉段的几何特 征，对交叉点或交叉段进行查找并标记，然后利用纤维边缘及骨架在交叉点或交叉段附近 的连续性和可微性，通过计算平均斜率，从而实现纤维的准确分离，它的具体步骤为，
(1)采集纤维图像，利用自动图像采集设备采集纸浆纤维的图像到计算机上，图像标准 是256位bmp格式；每采集一副图像就对其进行预处理并将测量得到的数据存入数据库；
(2)图像预处理，用阈值的方法对当前图像IpInput进行二值化，使背景变为白色，纤 维和噪声成为黑色；再二值化后的图像IpInput进行闭运算；然后进行骨骼化，使目标物成 为单一像素宽度；然后去除非纤维目标；(3)判断纤维种类，当图像存在N个连通域时，其中N l，分别对每个连通域Ipone进行判断是否是交叉纤维，若是则执行第(4)步，若不是则进行第(5)步；若交叉点或交叉段 太多(交叉点大于2个，交叉段大于6个)则进行第(6)步。(4)此连通域Ipone是交叉纤维，则将交叉部分及其每个点的邻域是目标点的像 素值改为中间值，剩下的是黑像素的点就被分离为四个子连通域lpl，1ρ2, 1ρ3, 1ρ4 ；对这 四个子连通域内的点分别与交叉部分中的一个确定的点计算平均斜率，斜率最相近的就是 属于同一根纤维；将配对好的纤维存放在两个图像空间Ipoutl和lpout2中，被设成中间 值的交叉部分同属于两根纤维，将其同时分配给Ipoutl和lpout2，并且将中间值改为黑像 素；对Ipoutl中的纤维去除小分枝然后测量小分枝即黑像素点的个数是五个以内的纤维； 再对lpout2中的纤维也去除小分枝后完成测量；
(5)此连通域Ipone是单根纤维，直接进行去除小分枝，进行测量；
(6)此连通域Ipone的交叉点多于所设定的阈值时，则说明是噪声，直接舍弃，不再对 其进行任何操作。所述步骤(3)中，对一个连通域中交叉点进行判断时，先判断当前点的8邻域，如 果有大于等于3个目标点即n8>=3，则继续判断16邻域内目标点的个数η16，若除去8邻 域中的目标点仍有4个以上的目标点即η16-η18>=4，则继续判断24邻域内目标点的个数 π24，若除去16邻域中的目标点仍有4个以上的目标点即η24-η16>=4，并且这4个以上的目 标点记个数为ncross，若ncross个目标点中有位置相连的，则只留一个，之后若剩下的目 标点有四个即代表四个不同的方向，则记录当前像素点位置，并统计这种交叉点的个数n， 判断完这一连通域后，将数组里的点和其8邻域的像素值都改为中间值，以使剩下的黑像 素点成为独立的四部分，即步骤(4)中的四个子连通域lpl，1ρ2, 1ρ3, 1ρ4 ；对这四个子连通 域进行标记，标记值为lpl, 1ρ2, 1ρ3, 1ρ4 ；
由于交叉点存入数组的顺序是从上到下从左到右的，取中间点为参考点m，在对上述四 个独立部分，它们的标记值分别为lpl，1ρ2, 1ρ3, 1ρ4，进行配对；配对原则具体为在标记 为Ipl的子连通域中从离m点最近的目标点开始依次取五个连续目标点，分别与m点计算 斜率；若目标点横坐标的值与m点的横坐标的值相同时，会导致所计算的斜率无穷大，这时 则将斜率自动赋予一个很大的数q (q=65535)，并将斜率的值依次按顺序存入数组1中，同 理标记为1ρ2，1ρ3，1ρ4的子连通域也分别依照上述方法计算斜率并存放到数组2，数组3， 数组4中；之后再对每个数组求平均斜率，比较这四个平均斜率，斜率值最相近的则为一根 纤维，将属于一根纤维的目标点存入空间lpoutl，另一部分存入lpout2，而标记为中间灰 度值gray (128)的交叉部分像素作为两根纤维的公共部分存入以上两个纤维空间的相应 位置中，并将中间值改回黑像素值。因为交叉点存入数组的顺序是从上到下从左到右的，取中间点为参考点m，在对 (1)中四个独立部分(标记值为lpl，1ρ2，1ρ3，1ρ4)进行配对，配对原则具体为在标记为 Ipl的子连通域中从离m点最近的目标点开始依次取五个连续目标点，分别与m点计算斜率 (若目标点横坐标的值与m点的横坐标的值相同时，会导致所计算的斜率无穷大，这时则将 斜率自动赋予一个很大的数q(q=65535))，并将斜率的值依次按顺序存入数组1中，同理标 记为1ρ2，1ρ3，1ρ4的子连通域也分别依照上述方法计算斜率并存放到数组2，数组3，数组 4中；之后再对每个数组求平均斜率，比较这四个平均斜率，值最相近的则为一根纤维，将属于一根纤维的目标点存入空间lpoutl，另一部分存入lpout2，而标记为中间灰度值gray (128)的交叉部分像素作为两根纤维的公共部分存入以上两个纤维空间Ipoutl和lpout2 的相应位置中，并将中间值改回黑像素值。本发明的有益效果是本发明通过精确判断交叉点，得到交叉纤维的公共部分，以 达到分离的目的的方法，并采用平均斜率变化最小的方法，提高了配对的准确性。通过VC++ 实现的交叉纤维分离的结果和大量数据显示了这一方法的准确性和鲁棒性。分离后的纤维 分别显示在两个图像空间中，结果显而易见。
图Ι-a交叉纤维；
图l_b交叉点或段查找图； 图1-c交叉点或段标记图Ι-d斜率最相近的两个纤维子段分配给一根纤维Ipoutl图； 图Ι-e两个纤维子段(1ρ2，1ρ4)属于另一根纤维lpout2图； 图ι-f被标记为中间灰度值的交叉点或交叉段属于两根纤维的公共部分； 图Ι-g分离后的两根纤维； 图Ι-h分离后的两根纤维； 图1-m为Ipoutl中的纤维去除小分叉后的图； 图1-n为lpout2中的纤维去除小分叉后的图； 图2两根交叉纤维分离算法流程图； 图3为现有方法误判的图； 图4为现有方法漏判的图5是采集到交叉纤维原图，格式是512 512大小256色位像；
图6是二值化后的结果；
图7是边界平滑和孔洞填充后的结果；
图8是骨骼化的结果；
图9是去除非纤维区域后的结果；
图10是分离后的结果；
图11小分叉去除后的结果。
具体实施例方式下面结合附图与实施实例对本发明作进一步说明。一种自动纸浆纤维检测中分离交叉纤维的方法，它的步骤为，
(1)采集纤维图像。利用自动图像采集设备采集纸浆纤维的图像到计算机上，图像标准 是标准灰度位像。每采集一幅图像IpInput就对其进行以下处理并将测量得到的数据 存入数据库，同时单次测量结果及统计结果。(2)图像预处理。预处理包括图像增强、二值化、纤维边界平滑、孔洞填充、骨骼 化、非纤维区域去除等。(3)判断纤维种类。当图像存在N(其中N l)个连通域时，分别对每个连通域Ipone进行判断是否是交叉纤维，若是则执行第(4)步，若不是则进行第(5)步，若交叉点或交叉段太多则进行第(6)步。(4)分离交叉在一起的两根纤维。连通域Ipone是交叉纤维，示意图如图；查 找交叉点或交叉段，如图l_b所示；如图1-c所示，将交叉点或交叉段的像素值标记为中间 灰度值(图1-c中视白色圆圈部分)，剩下的黑像素点集合就成为了四个子连通域，即两根 纤维的四个子段，标记为Ipl、1ρ2、1ρ3和1ρ4 ；计算Ipl、1ρ2、1ρ3和1ρ4在交叉区域附近 的斜率kl，k2，k3，k4，如图1-c所示，分别在lpl，1ρ2，1ρ3，1ρ4上的虚线矩形框就代表参 与斜率计算的五个连续目标点，根据配对原则将两个子连通域(lpl，1ρ3)存放到图像空间 lpoutl，如图Ι-d，另外两个子连通域(1ρ2，1ρ4)存放到图像空间lpout2，如图l_e ；被标记 为中间灰度值的交叉点或交叉段属于两根纤维的公共部分，如图ι-f，同时分配给两个图像 空间Ipoutl和lpout2，得到分离后的两根纤维，如图Ι-g，图Ι-h。对Ipoutl中的纤维去 除小分枝(黑像素点的个数是五个以内)，如图1-m，对lpout2中的纤维做同样的操作，如图 1-n。至此实现了两根交叉纤维的准确分离。然后分别对两根分离后的纤维进行特征参数 测量。(5)此连通域Ipone是单根纤维，直接进行去除小分枝，参数测量等。(6)此连通域Ipone的交叉点多于所设定的阈值时，则认为是絮聚纤维区域，说明 是噪声，直接舍弃，不再对其进行任何操作。所述步骤(4)中，
(1)对一个连通域中交叉点进行判断时，先判断当前点的8邻域，如果有大于等于3个 目标点即n8>=3，则继续判断16邻域内目标点的个数nl6，若除去8邻域中的目标点仍有 4个以上的目标点即nl6-nl8>=4，则继续判断24邻域内目标点的个数n24，若除去16邻 域中的目标点仍有4个以上的目标点即η24-η16>=4，并且这4个以上的目标点记个数为 ncross，若ncross个目标点中有位置相连的，则只留一个，之后若剩下的目标点有四个即 代表四个不同的方向，则记录当前像素点位置，并统计这种交叉点的个数n，判断完这一连 通域后，将数组里的点和其8邻域的像素值都改为中间值，以使剩下的黑像素点成为独立 的四个子连通域lpl、1ρ2、1ρ3和1ρ4。(2)因为交叉点存入数组的顺序是从上到下从左到右的，取中间点为参考点m，在 对(1)中四个子连通域(标记值为Ipl，lp2，lp3，lp4)进行配对，配对原则具体为在标记为 Ipl的子连通域中从离m点最近的目标点开始依次取五个连续目标点，分别与m点计算斜率 (若目标点横坐标的值与m点的横坐标的值相同时，会导致所计算的斜率无穷大，这时则将 斜率自动赋予一个很大的数q，q=65535)，并将斜率的值依次按顺序存入数组1中，同理标 记为1ρ2，1ρ3，1ρ4的子连通域也分别依照上述方法计算斜率并存放到数组2，数组3，数组 4中；之后再对每个数组求平均斜率，比较这四个平均斜率，值最相近的则为一根纤维，将 属于一根纤维的目标点存入空间lpoutl，另一部分存入lpout2，而标记为中间灰度值gray (128)的交叉部分像素作为两根纤维的公共部分存入以上两个纤维空间lpoutl和lpout2 的相应位置中，并将中间值改回黑像素值。两根交叉纤维分离算法流程图如图2所示。实施例
假定当前采集到的图像就是两根纤维交叉的图像如图5。用阈值的方法对当前图像lplnput进行二值化处理，使背景变为白色，纤维和一些噪声(干扰)成为黑色，结果如图 6所示；二值化后的图像lplnput会因为光线的原因纤维的某个部分的像素值未达到阈值， 从而变成背景色，需要进行闭运算进行纤维边界平滑和孔洞填充，如图7所示；然后进行骨 骼化，使目标物成为单一像素宽度，得到的结果是图8;然后去除非纤维目标，从而得到图 9 ；此连通域lpone是交叉纤维，则将交叉部分及其每个交叉点的邻域是目标点的像素值改 为中间灰度值，剩下的是黑像素的点就被分离为四个子连通域，对这四个子连通域内的点 分别与交叉部分中的一个确定的点计算平均斜率，斜率最相近的就是属于同一根纤维，将 配对好的纤维存放在两个图像空间lpoutl和lpout2中，被设成中间值的交叉部分属于两 根纤维的公共部分，将交叉部分存入两个图像空间，并且将中间值改为黑像素。对lpoutl 中的纤维去除小分枝，测量等，再对lpout2中的纤维做同样的操作。实验步骤
(1)采集纤维图像。利用自动图像采集设备采集纸浆纤维的图像到计算机上，图像标准 是标准灰度位像。每采集一幅图像lplnput就对其进行以下处理并将测量得到的数据 存入数据库，同时单次测量结果及统计结果。(2)图像预处理。预处理包括图像增强、二值化、纤维边界平滑、孔洞填充、细化、 非纤维区域去除等。(3)对预处理后的图像的各连通区域判断纤维种类。经过预处理后的纤维图像中 包括单根纤维、交叉纤维或絮聚纤维等多种目标区域。设各连通域用lpone表示，分别判 断每个连通域lpone是否为两根纤维交叉情况，若是则执行第(4)步，若不是则进行第(5) 步，若为絮聚纤维区域则进行第(6)步。(4)分离交叉在一起的两根纤维，如图9所示。连通域lpone是交叉纤维，查找交 叉点或交叉段；将交叉点或交叉段的像素值标记为中间灰度值，剩下的黑像素点集合就成 为了四个子连通域，即两根纤维的四个子段，标记为Ipl、lp2、lp3和lp4;计算Ipl、lp2、lp3 和lp4在交叉区域附近的平均斜率kl，k2, k3, k4，将斜率最相近的两个子连通域分配给图 像空间lpoutl，如图10左所示，另外两个子连通域(lp2，lp4)属于另一图像空间lpout2， 如图10右所示；被标记为中间灰度值的交叉点或交叉段属于两根纤维的公共部分，如图10 虚线圆所圈部，同时分配给两个图像空间lpoutl和lpout2，得到分离后的两根纤维，如图9 所示。分别对纤维lpoutl和lpout2去除小分枝(黑像素点的个数是五个以内)，结果如图 11。至此实现了两根交叉纤维的准确分离。然后分别对两根分离后的纤维进行特征参数测 量。(5)此连通域lpone是单根纤维，直接进行去除小分枝，参数测量等。(6)此连通域lpone的交叉点多于所设定的阈值时，则认为是絮聚纤维区域，说明 是噪声，直接舍弃，不再对其进行任何操作。实验证明，此方法分离效果明显，提高了交叉点或交叉段的定位精度和分离纤维 的准确度，具有很高的鲁棒性。此外，此方法还可以用来去除细化骨架上的小分枝。
权利要求
一种自动纸浆纤维检测中分离交叉纤维的方法，它通过分析交叉点或交叉段的几何特征，对交叉点或交叉段进行查找并标记，然后利用纤维边缘及骨架在交叉点或交叉段附近的连续性和可微性，通过计算平均斜率，从而实现纤维的准确分离，其特征是它的具体步骤为，(1)采集纤维图像，利用自动图像采集设备采集纸浆纤维的图像到计算机上，图像标准是256色bmp图像格式；(2) 图像预处理，用阈值的方法对当前图像lpInput进行二值化，使背景变为白色，纤维和噪声成为黑色；再二值化后的图像lpInput进行闭运算；然后进行骨骼化提取骨架，使目标物成为单一像素宽度；然后去除非纤维目标；(3) 判断纤维种类，当图像存在N个连通域时，其中N 1，分别对每个连通域lpone进行判断是否为交叉纤维，若是则执行第(4)步，若不是则进行第(5)步；若交叉点大于2个或交叉段大于6个则进行第(6)步；(4)此连通域lpone是交叉纤维，则将交叉部分及其每个点的邻域是目标点的像素值改为中间值，剩下的是0值像素点就被分离为四个子连通域lp1,lp2,lp3,lp4；对这四个子连通域内的点分别与交叉部分中的一个确定的点计算平均斜率，斜率最相近的则认为属于同一根纤维；将配对好的纤维存放在两个图像空间lpout1和lpout2中，被设成中间值的交叉部分同属于两根纤维，将其同时分配给lpout1和lpout2，并且将中间值改为0值；对lpout1中的纤维去除小分枝然后测量，小分枝即黑像素点的个数是五个以内的纤维；再对lpout2中的纤维也去除小分枝后完成测量；(5)此连通域lpone是单根纤维，直接进行去除小分枝，进行测量；(6)此连通域lpone的交叉点多于所设定的阈值时，则说明是噪声，直接舍弃，不再对其进行任何操作。
2.根据权利要求1所述的自动纸浆纤维检测中分离交叉纤维的方法，其特征是所述 步骤(3)中，对一个连通域中交叉点进行判断时，先判断当前点的8邻域，如果有大于等于 3个目标点即n8>=3，则继续判断16邻域内目标点的个数η16，若除去8邻域中的目标点 仍有4个以上的目标点即η16-η18>=4，则继续判断24邻域内目标点的个数η24，若除去16 邻域中的目标点仍有4个以上的目标点即η24-η16>=4，并且这4个以上的目标点记个数为 ncross，若ncross个目标点中有位置相连的，则只留一个，之后若剩下的目标点有四个即 代表四个不同的方向，则记录当前像素点位置，并统计这种交叉点的个数n，判断完这一连 通域后，将数组里的点和其8邻域的像素值都改为中间值，以使剩下的黑像素点成为独立 的四部分，即步骤(4)中的四个子连通域lpl，1ρ2, 1ρ3, 1ρ4 ；对这四个子连通域进行标记， 标记值为 lpl, 1ρ2, 1ρ3, 1ρ4 ；由于交叉点存入数组的顺序是从上到下从左到右的，取中间点为参考点m，在对上述四 个独立部分，它们的标记值分别为lpl，1ρ2, 1ρ3, 1ρ4，进行配对；配对原则具体为在标记 为Ipl的子连通域中从离m点最近的目标点开始依次取五个连续目标点，分别与m点计算 斜率；若目标点横坐标的值与m点的横坐标的值相同时，会导致所计算的斜率无穷大，这时 则将斜率自动赋予一个数q，并将斜率的值依次按顺序存入数组1中，同理标记为1ρ2，1ρ3， 1ρ4的子连通域也分别依照上述方法计算斜率并存放到数组2，数组3，数组4中；之后再对 每个数组求平均斜率，比较这四个平均斜率，斜率值最相近的则为一根纤维，将属于一根纤维的目标点存入空间lpoutl，另一部分存入lpout2，而标记为中间灰度值gray的交叉部分像素作为两根纤维的公共部分存入以上两个纤维空间的相应位置中，并将中间值改回黑像素值。
全文摘要
本发明涉及一种自动纸浆纤维检测中分离交叉纤维的方法。它解决了在测量单根纤维长度等参数时无法对两根交叉纤维进行精确测量，产生严重误差的问题，在存在小分枝等干扰因素下，该方法能成功地将交叉在一起的纤维准确分离开，并分别对分离后的纤维进行特征参数测量等操作。(1)采集纸浆纤维的图像到计算机；(2)对采集的图像进行预处理；(3)对预处理后的图像的各连通区域lpone判断纤维种类；(4)若连通域lpone为两根纤维交叉，则分离交叉在一起的两根纤维；(5)若连通域lpone非两根纤维交叉，则此lpone是单根纤维，直接进行去除小分枝，进行参数测量；(6)若连通域lpone为絮聚纤维区域，则认为是噪声，不再对其进行参数测量。
文档编号G06T5/00GK101866491SQ20101021717
公开日2010年10月20日 申请日期2010年7月5日 优先权日2010年7月5日
发明者孟荣爱, 张凯丽, 李庆华, 杨秀蔚, 邱书波 申请人:山东轻工业学院