专利名称：一种异形纤维的尺度形状特征的提取方法
技术领域：
本发明涉及一种异形纤维的尺度形状特征的提取方法，属于异形 纤维形状特征检测技术领域。
背景技术：
纤维截面的异形化是物理改性的一项重要手段，它通过模拟天 然纤维的各种形态改善合成纤维的服用性能。各种异形截面纤维在国 外已大量生产，国内也有研制和批量生产，其中中空、中空异形、三
角形、三叶形等使用已非常广泛，多角形、多叶形、T形、C形、H形、 十字形、扁平形等也已有生产。异形截面纤维不规则的截面形状赋予 织物独特的服用和外观性能。
截面特征参数是表征截面形状的主要指标，也是研究服用性能的 基础之一。目前对异形纤维截面的识别主要还限于人工法或半自动 法。在异形纤维检测领域广泛采用视频显微镜、普通生物显微镜测试 异形截面纤维特征参数。视频显微镜通过与之配套兼容的 PinnacleStudio8. 0图像采集软件，采集在计算机屏幕上观察到的 一定放大倍数下的异形纤维截面图像，再利用与视频显微镜配套兼容 的Imager- Proplus图像分析软件，测量异形纤维截面的内切圆半 径、外接圆半径、异形纤维截面周长和面积、异形纤维截面外接圆周 长和面积、中空纤维截面的中空面积、"V"字形纤维的角度等特征
5参数，由此计算出异形纤维的异形度指标。生物显微镜法用装在显微 镜上的显微描绘仪在描图纸上描绘异形纤维截面图像，再用复印机放 大2 4倍，剪纸、称重或用软尺测量，计算得到异形截面纤维异形 度指标。
这两种测试方法，需要检测人员手工操作，方法费时，费力，准 确性差，对检测人员的主观技术经验依赖性强，且只适用于已知异形 纤维形态测量异形度，不能进行异形纤维的自动识别、分类和统计。
发明内容
本发明的目的在于提出一种异形纤维的尺度形状特征的提取方 法，该方法在实现异形纤维横截面显微图像全自动拍摄采集基础上， 利用异形纤维边沿的凹凸特性自动识别异形纤维的几何形状，准确地 对异形纤维自动分类，实现了异形纤维自动分类和混合纤维比例自动 计算。
为达到上述目的，本发明的技术方案是利用三维自动载物台、 光学显微镜、工业控制计算机和数字摄像机组成图像采集系统，由控 制软件经串行口分别对三维自动载物台XYZ方向的步进电机控制， 实现载物平台在水平和垂直方向的移动控制，使光学显微镜工作在最 佳聚焦状态进行图像采集。
本发明的提取方法的步骤是首先给定'异形纤维横截向样片'， '装配三维自动平载物台的显微镜'经'连续迂回图像采集，、'自动 聚焦'到对图像进行'是图像目标？'判断，如果否，返回到'连续 迂回图像采集'如果是，到'分图像采集'进行分图像采集，然后可以返回到^连续迂回图像采集'继续采集，直到分图像采集结束进到 '大图拼接'，得到'清晰的纤维图像全景图'、经对采集到的图像进
行'图像分割'、'链码表示''TAR计算'、'提取TARN'、'提取TARP'、 '纤维几何形状确定'最后'纤维分类、统计'。
通过上述步骤，拍摄清晰的异形纤维横截面图像，在试样片的水 平面上找到各个纤维的位置，其次是针对各个纤维进行垂直方向的聚 焦。利用三维自动载物台的光学显微镜，通过以下的步骤在水平(X, Y)和垂直(Z)方向自动搜索纤维并采集图像，实现异形纤维横截 面显微图像全自动拍摄采集和自动分类统计。
本发明的有益效果采用本发明的提取方法，省时、省力，准确 性高，客观，可以自动获得试样片的全景图和纤维的分类统计结果， 取代人工检测异形纤维，具有一定的使用价值。
以下结合附图和实施例对本发明的方案作比较详细的说明。
图1为本发明的异形纤维提取方法的步骤流程图2为本发明的异形纤维提取方法的连续迂回采集示意图3为本发明的异形纤维提取方法的采集时毗连的图像间留有
重叠区域示意图4为本发明的异形纤维提取方法的拼接预处理后的全景图； 图5为本发明的异形纤维提取方法的异形纤维几何形状自动确
定示意图。
具体实施方式
参照图1，这是本发明的异形纤维提取方法的步骤流程图。 如图所示，利用图像提取设备的控制软件经串行口分别对三维自 动载物台XYZ方向的步进电机控制，实现载物平台在水平和垂直方向
的移动控制，使光学显微镜工作在最佳聚焦状态对图像进行采集和分 类统计。
首先要对异形纤维自动分类，包括图像分割，链码表示，三角形
面积(TAR)计算，三角形负面积(TARN)表示，三角形正面积(TARP) 表示，异形纤维分类。然后进行图像采集，图像采集是采集异形纤维 样片目标范围的图像内容，、实现步骤是先收索目标区域，全自动扫描 目标定位，对寻找到目标自动聚焦采集分图像。然后大图拼接，得到 异形纤维横截面的全景图。
本发明提取方法的步骤如下首先给定'异形纤维横截向样片'，
'装配三维自动平载物台的显微镜'经^连续迂回图像采集，、'自动 聚焦'到对图像进行'是图像目标？'判断，如果否，返回到'连续 迂回图像采集'如果是，到'分图像采集'进行分图像采集，然后返 回到'连续迂回图像采集'继续采集，直到分图像采集结束进到'大 图拼接'，得到"清晰的纤维图像全景图'、经对采集到的图像进行'图 像分割'、'链码表示，'TAR计算'、'提取TARN'、'提取TARP'、
'纤维几何形状确定'最后'纤维分类、统计'。
参照图2，这是本发明的异形纤维提取方法的连续迂回采集示意图。
如图所示，连续迂回的采集，通过三维自动载物台光学显微镜的控制软件，控制显微镜的三维自动载物台通过迂回的连续采集方式，获 取试样片上纤维目标区域。
参照图3，这是本发明的异形纤维提取方法的采集时毗连的图像 间留有重叠区域示意图。
如图所示，为了保证在对分图像进行拼接时不遗漏纤维截面的信 息，各个分图像在采集时都需要和毗连的图像间有一公有的重叠区
域。系统中每一个分图像的大小为640*480像素，左右的重叠部分取 128个像素，上下的重叠部分取98个像素。这样目标区域被划分为 多个分图像子图，计算机控制移动载物台聚焦自动采集分图像信息数 据，同时纪录下分图像的坐标和序号，这些信息用于确定分图像在全 景图缝合中的位置，以便后续拼图。
参照图4，这是本发明的异形纤维提取方法的拼接预处理后的全 景图。
如图所示是拼接预处理后的全景图。全图拼接是对每一张分图像 经裁剪后根据原先的编号顺序进行拼接，得到全景图。
参照图5，这是本发明的异形纤维提取方法的异形纤维几何形状 自动确定示意图。
如图所示，异形纤维几何形状自动确定的实现方法 通过对异形纤维图像边沿链码的TAR计算，得到TARN，确定边 沿的凹点个数和在边沿上的位置；得到TARP，确定边沿的凸点个数 和在边沿上的位置，计算机根据边沿点上凹凸点的个数和在边沿上的 分布位置自动识别异形纤维的形态。下面将有关问题作比较详细的说明。
所述的全自动扫描目标定位，是在载物台运动的过程中，计算机 不断计算每一个区域图像灰度的方差，对纤维目标进行锁定。载玻片 的图像灰度比较浅，而且变化不大，其方差近似相等，数值较小；纤 维目标区域的图像灰度较深，且变化较大，其方差增大。根据方差的 变化区分载波片和目标区域，保存下纤维目标区域的起点坐标和终点 坐标，锁定纤维目标区域，实现扫描目标的全自动定位。
分图像采集，为了保证纤维目标图像的清晰度，显微镜的放大倍 数取100倍，因此整个目标纤维的横截面不可能一次全部采集到，必 须采用分图像采集的方法。为了保证在对分图像进行拼接时不遗漏纤 维截面的信息，各个分图像在采集时都需要和毗连的图像间有一公有
的重叠区域。系统中每一个分图像的大小为640x480像素，左右的重 叠部分取128个像素，上下的重叠部分取98个像素。这样目标区域 被划分为多个分图像子图，计算机控制移动载物台聚焦自动采集分图 像信息数据，同时纪录下分图像的坐标和序号，这些信息用于确定分 图像在全景图缝合中的位置，以便后续拼图。
自动聚焦，系统中自动聚焦和全景图扫描通过计算机串行口对各 个基本控制合理的结合实现。自动聚焦由算法中定义的两个参数确 定FP参数(图像的聚焦参数)，用来确定自动平台的Z轴的正确聚 焦位置，通过图像中心区域计算获得;参数CFP(正确聚焦位置参数)， 用于确定采集到最大的FP图像时自动平台的Z轴的位置。系统的聚 焦实现方法从一系列不同Z轴位置采集的图像中，找具有最大FP像采集时的Z轴位置作为最佳聚焦点。
大图拼接，对于采集到的纤维分图像，必须进行拼接，将一个载 波片上的纤维横截面信息集中在一个总图上，以便进行纤维成份分 析。要求拼接后的大图尽可能的保证每一个纤维的完整性和不重叠 性，失真尽可能小，没有明显的缝合线，拼接的方法必须科学合理。 分析纤维图像发现，在纤维目标区域中，背景和纤维的灰度又有差异， 纤维图像的背景点图像值较小，方差变化不大；纤维目标图像的灰度 值较大，方差变化较大。可以根据重叠区域方差的数值确定图像的拼 接线。从第一个分图像开始逐个对相毗邻的两个分图像分别计算每一 个像素八邻域的灰度方差，在相应的位置保存得到的计算结果。
横向裁剪，对左右重叠区域寻找每一行方差的最小值，将得到的 方差最小值连接得到一条最小方差曲线，这条曲线一定在图像的背景 点上，以这条曲线作为两张左右毗邻分图像的裁剪线，对左边的分图 像裁剪掉裁剪线右边的重叠部分，对右边的图像裁剪掉裁剪线左边的 重叠部分。
纵向裁剪，对上下重叠区域寻找每一列方差的最小值，将得到的 方差最小值连接得到一条最小方差曲线，这条曲线一定在图像的背景 点上，以这条曲线作为两张上下毗邻分图像的裁剪线，对上面的分图 像裁剪掉裁剪线下面的重叠部分，对下面的图像裁剪掉裁剪线上面的 重叠部分。，
全图拼接，对经裁剪后的每一张分图像根据原先的编号顺序进行 拼接，就得到全景图。上述的图像分割，分以下几个方面
1) .原始图像的预处理，包括灰度处理、去噪、图像平滑；
2) .纤维轮廓提取，对预处理后的图像进行边缘检测，提取纤维 对象的连续轮廓；
3) .纤维图像分离，对经图像拼接预处理得到的纤维图中的单个
纤维对象进行提取，获取横截面中所有的单个纤维对象。
所述的链码表示，分以下几个方面
1. 通过细化的算法去除单个纤维边沿的冗余元素和分岔，使每 个边界点在序列中仅有两个相邻边界点，得到最短闭合边沿链码。
2. 单个纤维连续轮廓信息用N个链码保存，设第一个链码为
(xo，w)，最后一个链码为(A.n。，轮廓链码的保存格式-
，i
■^n-liVn-l
异形纤维轮廓三角形面积表示(TAR)的实现方法 用行列式计算由顶点P。"，h), A"，力)和A"，h)组成的三角形面积 S， S的绝对值和符号与顶点的空间位置有关。改变行列式中任意两 行的顺序，行列式结果符号将发生变化。所以，符号与顶点序列的弯 曲方向有关。
在由^个点构成的异形纤维闭合边界路径P上，沿逆时针方向绕行对 每一个元素/ ,.(0&'<^)依次计算三角形面积&，下标W人0增加到 7V-1 。
定义步长Z为介于BKWC之间的一个整数。当w为偶数时，元素
P, (0< iV )表示三角形的中间顶点，前一个顶点的下标表示为
12卜r, /a
jV + / - elsewhere (1)
后一个顶点的下标表示为
z' + /， / + " jV
/ + /1 - iV， elsewhere (2) 所构成的三角形面积记为S,.，，，表示顶点为"步长为纟的三角形面积。 定义
,r (4) 式(3)的TAR称为闭合2D序列的多尺度描述算子。
TARN提取以及对凹点特征描述的实现方法
按逆时针旋转方向沿闭合边界顺序选取三个顶点组成三角形，当 边界出现凹点时，其局部弯曲方发生改变，三角形面积值变为负值。 因此边沿序列凹点组成的三角形面积值TAR为负值，用 窗# = {咖,,}表示。定义
jVw,,, to*,, < o ,、
tor", f = "j , ,。 (5)
' elsewhere
对不同角度的凹点形态，其对应的TARN输出值不同，TARN的 形状也不同，若凹点的两边不对称，多尺度TARN值的分布向凹点 的较长边倾斜。TARN图形的底部随着角度的增大而增宽，可以用底 部的宽度判断凹点的形态。
TARP提取以及对凸点特征描述的实现方法 边界上凸点可以通过下列步骤确定1. 对于不同步长hl,2,A,7V/2，重复下面计算。
2. 通过公式(4)计算当/ = 0， 1,2, A ， iV -1时的 一,,值。
3. 对每个步长,计算to^的平均值和/^,与平均值的差值
，'1&气' (6)△ , 二 一 , — avg, 。 (7)
4. a,"。,&。，-OTg,大于零表示位置o是凸点。用r乂及/^(似iy,、j表
不o
「△— ,， △ m。 ' >0
" = 3 (8)' ^0， elsewhere
不同步长^对应的TARP输出称为TARP的多尺度结果。
异形纤维几何形状自动确定的实现方法
通过对异形纤维边沿链码的TAR计算，得到TARN，确定边沿的凹点个数及位置；得到TARP，确定边沿的凸点个数及位置，计算机根据边沿点上凹凸点的个数和在边沿上的分布位置自动识别异形纤维的形态。
自动分类、统计的实现方法
根据上述方法可以将一个横截面图像中拥有的异形纤维种类分开，对每一种类型的异形纤维统计个数，计算出横截面上各种形状纤维的总数，用每一种异形纤维的数量和各种形状纤维的总数进行计算，就能统计出每一种异形纤维在横截面中的比例，从而算出每一种异形纤维的含量。 像采集(1) 连续迂回的采集通过三维自动载物台光学显微镜的控制软件，控制显微镜的三维自动载物台通过迂回的连续采集方式，获取试样片上纤维目标区域。
(2) 自动聚焦在搜索到纤维目标的区域进行自动聚焦和分图像采集。
(3)大图拼接
对采集到的分图像集合进行拼接，得到异形纤维横截面试样片指定区域的全景图。
2.纤维分类
首先对纤维图像进行预处理，然后对图像的凹凸性进行提取，得到异形纤维的几何形状，给出纤维的分类结果。
(1) 图像分割从采集到的纤维图像中提取纤维的轮廓信息。
(2) 链码表示用链码形式表示纤维图像的边沿信息。
(3) 三角形面积(TAR)计算沿纤维边沿按逆时针方向计算不同尺度对应的三角形面积。
(4) 三角形负面积(TARN)表示TARN反映纤维边沿的凹陷特征。
(5) 三角形正面积(TARP)表示TARP表征纤维边沿的凸起特征。
(6) 异形纤维分类根据TARN和TARP的特征得出异形纤维的边沿几何形状，实现异形纤维种类的分类。
(7) 计算得到异形纤维的混合比例
权利要求
1、一种异形纤维的尺度形状特征的提取方法，其特征在于利用提取设备，由控制软件经串行口分别对三维自动载物台XYZ方向的步进电机控制，实现载物平台在水平和垂直方向的移动控制，使光学显微镜工作在最佳聚焦状态；然后对异形纤维自动分类，包括图像分割，链码表示，三角形面积(TAR)计算，三角形负面积(TARN)表示，三角形正面积(TARP)表示，异形纤维分类，然后进行图像采集，图像采集是采集异形纤维样片目标范围的图像内容，实现步骤是先收索目标区域，全自动扫描目标定位，对寻找到目标自动聚焦采集分图像，然后大图拼接，得到异形纤维横截面的全景图；其提取方法的步骤如下给定‘异形纤维横截向样片’，‘装配三维自动平载物台的显微镜’经‘连续迂回图像采集’、‘自动聚焦’到对图像进行‘是图像目标？’判断，如果否，返回到‘连续迂回图像采集’如果是，到‘分图像采集’进行分图像采集，然后返回到‘连续迂回图像采集’继续采集，直到分图像采集结束进到‘大图拼接’，得到‘清晰的纤维图像全景图’、经对采集到的图像进行‘图像分割’、‘链码表示’‘TAR计算’、‘提取TARN’、‘提取TARP’、‘纤维几何形状确定’最后‘纤维分类、统计’。
2、如权利要求l所述的异形纤维的尺度形状特征的提取方法，其特征在于所述的异形纤维自动分类，包括得到异形纤维的几何形状给出纤维的分类结果；从采集到的纤维图像中提取纤维的轮廓信息；用链码形式表示纤维图像的边沿信息；三角形面积TAR计算；沿纤维边沿按逆时针方向计算不同尺度对应的三角形面积。
3、如权利要求1或2所述的异形纤维的尺度形状特征的提取方法，其特征在于所述的计算不同尺度对应的三角形面积1)、三角形负面积TARN表示反映纤维边沿的凹陷特征；2)、三角形正面积TARP表示表征纤维边沿的凸起特征；3)、异形纤维分类根据TARN和TARP的特征得出异形纤维的边沿几何形状，实现异形纤维种类的分类；4)、计算得到异形纤维的混合比例。
4、 如权利要求1所述的—异形纤维的尺度形状特征的提取方法，其特征在于所述的图像采集，是采集异形纤维样片目标范围的图像内容，先收索目标区域，全自动扫描目标定位，对寻找到目标自动聚焦采集分图像，然后大图拼接，得到异形纤维横截面的全景图。
5、 如权利要求1或4所述的异形纤维的尺度形状特征的提取方法，其特征在于所述的全自动扫描目标定位，在载物台运动的过程中，计算机不断计算每一个区域图像灰度的方差，对纤维目标进行锁定；载玻片的图像灰度比较浅，而且变化不大，其方差近似相等，数值较小；纤维目标区域的图像灰度较深，且变化较大，其方差增大，根据方差的变化区分载波片和目标区域，保存下纤维目标区域的起点坐标和终点坐标，锁定纤维目标区域，实现扫描目标的全自动定位。
6、 如权利要求1或4所述的异形纤维的尺度形状特征的提取方法，其特征在于所述的分图像采集，为保证纤维目标图像的清晰度，整个目标纤维的横截面不可能一次全部采集到，而采用分图像采集的方法，在对分图像进行拼接时不遗漏纤维截面的信息。
7、如权利要求1或4所述的异形纤维的尺度形状特征的提取方法，其特征在于所述的各个分图像采集，需要和毗连的图像间有一公有的重叠区域，系统中每一个分图像的大小为640*480像素，左右的重叠部分取128个像素，上下的重叠部分取98个像素，这样目标区域被划分为多个分图像子图，计算机控制移动载物台聚焦自动采集分图像信息数据，同时纪录下分图像的坐标和序号，这些信息用于确定分图像在全景图缝合中的位置，以便后续拼图。
全文摘要
本发明涉及异形纤维的尺度形状特征的提取方法，在经串行口对自动载物台XYZ方向的步进电机控制，实现在水平和垂直方向的移动控制，使光学显微镜工作在最佳聚焦状态。提取方法的步骤如下给定“异形纤维横截向样片”，“装配三维自动平载物台的显微镜”经“连续迂回图像采集”、“自动聚焦”到“是图像目标？”判断，如果否，返回到“连续迂回图像采集”如果是，到“分图像采集”，然后返回到“连续迂回图像采集”直到分图像采集结束进到“大图拼接”，得到“清晰的纤维图像全景图”、经对采集到的图像进行“图像分割”、“链码表示”“TAR计算”、“提取TARN”、“提取TARP”、“纤维几何形状确定”最后“纤维分类、统”。
文档编号G01N33/36GK101487838SQ20081020440
公开日2009年7月22日 申请日期2008年12月11日 优先权日2008年12月11日
发明者吴雄英, 唐莉萍, 曾培峰 申请人:东华大学