CN105403482B - 一种相对面密度的光学算法 - Google Patents
一种相对面密度的光学算法 Download PDFInfo
- Publication number
- CN105403482B CN105403482B CN201510703493.XA CN201510703493A CN105403482B CN 105403482 B CN105403482 B CN 105403482B CN 201510703493 A CN201510703493 A CN 201510703493A CN 105403482 B CN105403482 B CN 105403482B
- Authority
- CN
- China
- Prior art keywords
- sample
- fiber
- surface density
- region
- stratiform
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N9/00—Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity
- G01N9/24—Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity by observing the transmission of wave or particle radiation through the material
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明涉及一种由纤维或细薄物体组成的层状集合体相对面密度的光学测量计算方法,其特征在于,包括以下步骤:第一步、用均匀光源照射层状纤维集合体材料,利用数码成像技术获取试样的透射光数字灰度图像;第二步、计算层状集合体上任一点的相对面密度;或者计算层状集合体上任一个区域的平均面密度的相对值。运用本发明算法获得的纤维层相对面密度更加精确,从而为精确测量分析纤维长度分布、纤维或细薄物体堆砌成的层状材料结构和性能打下了基础。
Description
技术领域
本发明涉及纤维集合体结构分布特性的光学测量技术,具体来说,涉及由纤维或其他高分子物体堆砌而成的层状材料透光信号计算其面密度分布的高精度方法。
背景技术
层状纺织纤维集合体的面密度分布,是表征多种层状纤维集合体结构和质量分布等特性的中间变量,是研究纤维、纤维网、织物等纤维集合体结构和性能的重要依据,被广泛应用于纺织材料的光电法检测等领域,照影仪法测量棉纤维长度就是一例,但是,早期利用透光强度表征纤维层厚度的照影仪法存在很大误差。
2008年,严漂利用摄像头拍摄透射光照射的纤维束,获得数字图像,根据图像中像素点的亮度值与透光强度成正比的物理概念做出自称的照影仪曲线(东华大学,硕士论文:棉纤维色度和长度测量方法的研究与设计)。因透光量与纤维层厚度(相对面密度)不是线性关系,将透光量作为纤维层厚度计算的纤维层面密度存在较大的偏差。
2012年,王府梅、吴红艳(专利号:ZL201210328089.5)利用透射式扫描仪获取双端须丛的灰度值图像后,将灰度值作为透射光强代入考虑材料吸收与厚度关系的Lambert-Beer定律计算须丛的厚度分布信息,以下称此法为Lambert-Beer定律的衍生算法。虽然较前述方法可以在很大程度提高了须丛厚度的计算精度,但是与真实厚度相比,仍然存在一定偏差。后来吴红艳(东华大学,博士论文:一种纤维长度快速低成本测量方法——双须测量法)基于一组实验总结出了棉纤维须丛曲线的修正方法,但是,此经验公式也仅适用于棉纤维。
在理论上,由于纤维是半透明物质,当光照射纤维集合体时,一部分光线透过纤维集合体成为透射光,另一部分光线因纤维集合体对光线的作用而损失。损失的光线,除了包括被纤维吸收的光能量以外,还包括纤维集合体内部众多纤维表面产生的四面八方的反射光和纤维内部结构不均或消光剂等嵌入物产生的各个方向的散射光。因为基于单向吸收的Lambert-Beer定律没有考虑光被介质的多向散射作用和散射光进一步被吸收,所以计算出的纤维层面密度仍然存在明显偏差。为了高精度测量和分析各种纺织纤维集合体的结构和性质,需要建立考虑光的多向散射和吸收的光学理论模型。
许多研究者对纺织纤维的光学性质进行过实验研究与理论研究,但是综合考虑吸收、散射和反射的纤维层光学面密度研究却未见报道。
当光入射到均质片状材料时,Kubelka和Munk将材料对光的散射和吸收作用在各个方向上造成的光衰减近似折算到向上和向下两个通道光的变化,通过光衰减量推导出与材料的厚度、散射系数、透射系数、透射信号和反射信号相关的物理方程,此方程被称为Kubelka-Munk双通道模型。Kubelka-Munk理论模型已被广泛应用于油墨纸张、颜色混色、医药等均质材料以及膜层原位定量光密度分析领域,而在固体与气体混合构成的纤维集合体面密度定量计算方面未见研究或应用报道。
因此,为了准确测量和分析各种纺织纤维集合体的结构和性质,迫切需要能够综合考虑由纤维和空气混合构成的层状纤维集合体对光线的吸收和散射作用,采用纤维集合体的透射光信号计算其相对面密度分布信息的精准算法。
发明内容
本发明的目的是提供一种由纤维层透光信号计算其相对面密度分布的高精度计算方法。
为了达到上述目的,本发明的技术方案是提供了一种相对面密度的光学算法,其特征在于,包括以下步骤:
第一步、获取由纤维或其他高分子物体堆砌而成的层状材料试样,用均匀光源照射层状材料试样,利用数码成像技术获取层状材料试样透射数字灰度图像;
第二步、或者用公式(1)计算层状材料试样任一点的相对面密度Wr,即m行n列的透射光数字灰度图像矩阵中第i行第j列的像素点上材料层的面密度Wij与该层状材料试样最厚点的面密度Wmax的比值:
式中,R∞为所测纤维或其他高分子物体堆砌到无穷厚时的反射率,是纤维或其他高分子物体自身的光学参数;Tij=Iij/A为第i行第j列像素点处纤维层的透射率,A为该数字灰度图像中无纤维处像素点的灰度值(等于入射光强度),Iij为第i行第j列的像素点的灰度值即透射光强度;Tmax为该试样最厚处的透射率;
或者计算层状材料试样任一个区域平均面密度的相对值Wr,即计算第x个区域平均面密度Wij与该试样最厚区域平均面密度Wmax的比值,包括以下步骤:
(1)试样的数字灰度图像划分为N个不同区域,N≥2,每个区域的数字灰度图像均为m行n列的灰度值矩阵;
(2)计算任意一个区域即第x个区域的相对面密度Wr,x=1,…,N;设第x个区域的数字灰度图像上第i行第j列的像素点的透射率为Tij,厚度最大区域的数字灰度图像上第i行第j列的像素点的透射率为Tmaxij,则用公式(2)计算第x个区域的相对平均面密度Wr:
优选地,所述无穷厚材料的反射率R∞为常用的紫外可见分光光度仪在380nm~720nm波长范围内测量的层状材料厚度为0.5~5cm范围的反射率。
优选地,所述层状材料试样为层状纤维集合体试样或由高分子膜或碎屑堆砌而成的层状试样。
优选地,所述第一步包括:用均匀光源照射由纤维或其他高分子物体堆砌而成的层状试样,同时,利用数码成像技术获取该试样的透射数字灰度图像。
运用本发明算法获得的纤维层及高分子物体堆砌成的层状试样的相对面密度更加精确,从而为精确测量分析纤维长度分布、纤维集合体及高分子物体堆砌的层状材料的结构和性能打下了基础。
附图说明
图1为羊毛须丛透射光灰度图
图2为须从的相对质量随横截面位置的变化曲线
图3为纤维层排列及正方形测试孔;
图4A至图4J为实施例2中的第1个层状棉纤维集合体试样至第10个层状棉纤维集合体试样的透光图像;
图5为本发明计算的棉纤维层相对面密度与Lambert衍生算法比较;
图6A至图6J为实施例2中的第1个层状羊毛纤维集合体试样至第10个层状羊毛纤维集合体试样的透光图像;
图7为本发明计算的毛纤维层相对面密度与Lambert衍生算法比较;
图8A至图8Q为0层至16层叠加的聚酯膜透光图像;
图9为本发明计算的薄膜相对面密度与Lambert衍生算法、光强直接算法比较。
具体实施方式
为使本发明更明显易懂,兹以优选实施例,并配合附图作详细说明如下。
实施例1:羊毛须丛曲线的测量计算与效果考查
须丛制作:随机抽取出一段长度大于最长纤维两倍的毛条试样,用专用夹持器垂直于毛条长度方向夹持住毛条,然后用梳子对毛条一端进行梳理,去除浮游纤维,制得待测羊毛须丛,参见图1。
须丛曲线的物理意义:图1所示的须丛实际是一个层状纤维集合体,纺纱牵伸机构设计和纺纱工艺设计等许多领域都需要知道该须丛任一横截面上的纤维量与横截面位置的关系曲线,俗称须丛曲线。
须丛透光图像测量:将制作好的羊毛须丛放入光学检测装置,获取须丛的透光数字灰度图像,如图1所示。灰度值范围为0~255,分辨率为1000,精度为0.0254mm。平行和垂直于纤维轴建立图像的坐标系。
用本发明专利计算须丛曲线:根据公式(2)可以计算须丛任一横截面上即第j列像素点上纤维的相对面密度Wjr,Wjr=m行n列的数字灰度图像矩阵中第j列像素点上总计面密度与该试样最厚列(须丛起始横截面)的总面密度的比值,即:
式中,该纤维无穷厚时的反射率R∞取为0.6;Tij=Iij/A为须丛数字灰度图像上第i行第j列像素点处的透射率,A为数字灰度图像中无纤维处像素点的灰度值,Iij为须丛数字灰度图像上第i行第j列像素点的灰度值;Timax=Iimax/A为须丛数字灰度图中第i行最厚列(起始列)的透射率,Iimam为须丛数字灰度图像上第i行最厚列(起始列)像素点的灰度值。
以上面求得的第j列像素点上纤维的相对面密度Wjr为纵坐标,以横截面所在位置为横坐标,绘制须丛曲线,如图2中的黑点所示。
效果考查:另一方面,利用Almeter100仪器的电容式传感器测得该须丛任一横截面的相对质量随横截面位置的变化曲线,如图2中的“△”所示。图2中的2曲线高度一致,说明本发明专利计算的须丛曲线精确度很高。由于本发明适合采用目前高度发达的数码成像技术,比传统的电容式检测等方法更具应用价值。
实施例2:棉纤维层相对面密度的测量计算及其效果
为验证本发明的正确性,利用棉纤维制作了10个厚度依次增加的层状纤维集合体试样,每个试样的厚度均匀,参见图3中部2cm2的窗口。按厚度依次增加的顺序,分别被称作第1、第2、…第x个…、第10个层状纤维集合体试样。
利用数码成像设备在相同光照等测试条件下分别获取第1个、第2个、…第x个…、第10个层状纤维集合体试样中部1cm2的透射光灰度图像,如图4A至图4J所示。灰度值越高表明透射光越强,纤维层越薄。第x个试样的透光图像中第i行第j列像素点的透射率以Tij表示,则第10个即最厚试样第i行第j列象素点的透射率以Tmaxij表示。
采用公式(2),分别计算第1个、第2个、…第x个…、第10个层状纤维集合体试样的相对面密度,第x个试样的相对面密度Wr=第x个试样的平均面密度/第10个试样的平均面密度即
式中,R∞为棉纤维层无穷厚时的反射率,在本实施例中取为0.9;Tij=Iij/A为第x个试样的数字灰度图像上第i行第j列像素点处的透射率,A为数字灰度图像中无纤维处像素点的灰度值,Iij为第x个试样的数字灰度图像上第i行第j列像素点的灰度值;Tmaxij=Imaxij/A为第10个试样的数字灰度图像上第i行第j列像素点处的透射率,Imaxij为第10个试样数字灰度图像上第i行第j列像素点的灰度值。
根据Lambert定律的衍生算法,将第x个试样的平均面密度与最厚试样的平均面密度相比,得到该试样的相对面密度Wr,则有:
式中A为数字灰度图像中无纤维处像素点的灰度值即入射光强度,Iij即为第x个试样透光图像中第i行第j列像素点的灰度值;Imaxij为第10个试样透光图像的第i行第j列像素点的灰度值。
每个试样的基准相对面密度的测量方法如下:完成光学测量后的每个棉纤维层状试样,利用手术剪剪下其中部2×2cm2正方形测试孔内的纤维,利用精密天平对剪下的纤维进行称重,计算该试样的面密度,结果如表1所示。并将该面密度与最重纤维层即第10个试样的面密度相比,获得该试样的基准相对面密度。
本发明方法的效果,可以从图5看得很清楚,利用本发明方法计算的光学相对面密度与基准相对面密度非常一致,几乎是线性关系,而纤维领域目前使用的Lambert定律衍生算法的效果明显比较差,误差较大,说明本发明算法可以更精确地计算纤维层的相对面密度。
表1 棉纤维层的面密度
实施例3:羊毛纤维层相对面密度的测量计算及其效果
利用直径为19.7μm澳洲细毛纤维制作10个厚度依次增加的层状纤维集合体试样,每个试样的厚度均匀,试样具体情况参见实施例2。
获取每个试样透光图像的方法同实施例2,结果如图6A至图6J所示。
根据本发明算法及Lambert定律衍生算法计算每个毛纤维试样的相对面密度方法同实施例2。R∞取为0.6。
每个试样的基准相对面密度的测量方法同实施例2,结果如表2所示。
本发明方法的效果,可以从图7所示的与基准相对面密度关系看得非常清楚,同时也可看到至今使用的Lambert定律衍生算法的效果。显然,本发明算法得到的羊毛纤维层相对光学面密度与其基准面密度的一致性比Lambert定律衍生算法高得多。似乎在羊毛纤维领域更能凸显本发明算法的优越性,分析原因有二:(1)羊毛纤维的表面散射光更强,不可忽略;(2)实施例3的10个试样厚度差异大,如表1和表2所示。试样越厚纤维集合体内部的散射光越强,试样厚度差异越大越能显示考虑散射的本发明方法的优越性。
表2 羊毛纤维层的面密度
实施例4:多层叠加的聚酯薄膜累计厚度测量计算及其效果
即便不是纤维而是薄膜,遇到多层叠加情况时,若用光电方法测量累积厚度也存在各层表面对光线反射和内部非均匀结构等因素的散射问题,层间还可能存在空隙问题,都可视作集合体内部散射,用本发明方法计算其叠加厚度。
选用均匀的聚酯薄膜,将若干层薄膜平铺叠加,分别获取16层叠加、15层叠加…1层薄膜、0层薄膜的透光图像,如图8A至图8Q所示。
首先分析薄膜透光图的像素信息与相对面密度的定性关系。图8A为0层叠加,以此类推最末图像为16层叠加,层数越多透过的光强越低即灰度值越小,因此16层叠加的灰度值最小,0层叠加的灰度值最大。
根据本发明算法及Lambert定律衍生算法计算每个叠加薄膜试样相对面密度的方法同实施例2。
每个试样的基准相对面密度获取方法:由于每层薄膜是均质材料,将最厚16层叠加的聚酯薄膜的基准相对厚度看作1,则x层叠加聚酯薄膜试样的基准相对厚度为x/16。
早期曾经用透射光强度表征纤维集合体材料的厚度,按照该算法x层叠加的薄膜相对面密度式中为x层叠加薄膜的平均厚度;为16层叠加的薄膜厚度;透光强度Iij为x层叠加薄膜图像的第i行第j列像素点的灰度值;Imaxij为16层叠加薄膜透光图像的第i行第j列像素点的灰度值。
图9显示本发明计算的光学相对面密度与Lambert衍生算法、透光强度算法的效果比较,显然,本发明算法得到的薄膜相对光学面密度与其基准相对面密度几乎相等,是很直的线性关系,而早期使用的透光强度算法和Lambert衍生算法的误差明显较大,说明本发明算法可以更精确地反映多层叠加高分子材料的相对面密度。
Claims (4)
1.一种相对面密度的光学算法,其特征在于,包括以下步骤:
第一步、获取由纤维或其他高分子物体堆砌而成的层状材料试样,用均匀光源照射层状材料试样,利用数码成像技术获取层状材料试样的透射光数字灰度图像;
第二步、或者用公式(1)计算层状材料试样任一点的相对面密度Wr,即m行n列的透射光数字灰度图像矩阵中第i行第j列的像素点上材料层的面密度Wij与该层状材料试样最厚点的面密度Wmax的比值:
式中,R∞为所测纤维或其他高分子物体堆砌到无穷厚时的反射率,是纤维或其他高分子物体自身的光学参数;Tij=Iij/A为第i行第j列像素点处纤维层的透射率,A为该数字灰度图像中无纤维处像素点的灰度值,等于入射光强度,Iij为第i行第j列的像素点的灰度值即透射光强度;Tmax为该试样最厚处的透射率;
或者计算层状材料试样任一个区域平均面密度的相对值Wr,即计算第x个区域平均面密度Wij与该试样最厚区域平均面密度Wmax的比值,包括以下步骤:
(1)试样的数字灰度图像划分为N个不同区域,N≥2,每个区域的数字灰度图像均为m行n列的灰度值矩阵;
(2)计算任意一个区域即第x个区域的相对面密度Wr,x=1,…,N;设第x个区域的数字灰度图像上第i行第j列的像素点的透射率为Tij,厚度最大区域的数字灰度图像上第i行第j列的像素点的透射率为Tmaxij,则用公式(2)计算第x个区域的相对平均面密度Wr:
2.如权利要求1所述的一种相对面密度的光学算法,其特征在于,所述无穷厚材料的反射率R∞为常用的紫外可见分光光度仪在380nm~720nm波长范围内测量的层状材料厚度为0.5~5cm范围的反射率。
3.如权利要求1所述的一种相对面密度的光学算法,其特征在于,所述第一步包括:用均匀光源照射由纤维或其他高分子物体堆砌而成的层状试样,同时,利用数码成像技术获取该试样的透射光数字灰度图像。
4.如权利要求1所述的一种面密度的光学算法,其特征在于,所述层状材料试样为层状纤维集合体试样或由高分子膜或碎屑堆砌而成的层状试样。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510703493.XA CN105403482B (zh) | 2015-10-26 | 2015-10-26 | 一种相对面密度的光学算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510703493.XA CN105403482B (zh) | 2015-10-26 | 2015-10-26 | 一种相对面密度的光学算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105403482A CN105403482A (zh) | 2016-03-16 |
CN105403482B true CN105403482B (zh) | 2019-01-11 |
Family
ID=55469107
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510703493.XA Active CN105403482B (zh) | 2015-10-26 | 2015-10-26 | 一种相对面密度的光学算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105403482B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106407524B (zh) * | 2016-09-05 | 2019-11-08 | 东华大学 | 一种纺纱质量的定量预测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1442669A (zh) * | 2003-03-18 | 2003-09-17 | 周胜 | 实时检测纤维集合体平面产品厚度和纤维取向度的方法 |
CN102853775A (zh) * | 2012-09-06 | 2013-01-02 | 东华大学 | 须丛曲线的获取方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE19608930A1 (de) * | 1996-03-07 | 1997-09-11 | Frieder A Grieshaber | Verfahren zur Bestimmung der Dichte von Stoffen |
JP2000275193A (ja) * | 1999-03-25 | 2000-10-06 | Canon Inc | 表面密度の測定方法および測定装置 |
CN1156708C (zh) * | 2002-01-14 | 2004-07-07 | 清华大学 | 一种羊毛纤维细度测量方法及系统 |
CN1282869C (zh) * | 2004-08-10 | 2006-11-01 | 东华大学 | 纤维集合体变密度力学性能和密度分布测量装置与用途 |
CN101241082B (zh) * | 2008-03-14 | 2010-06-02 | 东华大学 | 一种非织造布纤维取向分布测量方法 |
CN202583071U (zh) * | 2012-01-06 | 2012-12-05 | 青岛检验检疫技术发展中心 | 织物经纬密度自动分析仪 |
CN103364397B (zh) * | 2012-03-31 | 2015-10-28 | 佛山市南海天富科技有限公司 | 一种织物纬密测量方法及装置 |
-
2015
- 2015-10-26 CN CN201510703493.XA patent/CN105403482B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1442669A (zh) * | 2003-03-18 | 2003-09-17 | 周胜 | 实时检测纤维集合体平面产品厚度和纤维取向度的方法 |
CN102853775A (zh) * | 2012-09-06 | 2013-01-02 | 东华大学 | 须丛曲线的获取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105403482A (zh) | 2016-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102253043B (zh) | 一种单目ccd多角度纱线外观数字化分析方法 | |
CN102105781B (zh) | 卷材检查校准系统及相关方法 | |
CN104596989B (zh) | 基于干涉条纹图处理测量透明介质折射率分布的方法 | |
CN202362250U (zh) | 一种单目ccd多角度纱线外观数字化分析装置 | |
US20040246493A1 (en) | Apparatus for measuring thickness profile and refractive index distribution of multiple layers of thin films by means of two-dimensional refrectometry and method of measuring the same | |
CN106517086A (zh) | 一种大面积高分辨率宽视场在线测量装置及其测量方法 | |
CN106323748A (zh) | 一种用于碳纤维复丝应变测量和力学性能表征的非接触视频测量系统 | |
Valla et al. | Comparison of two optical methods for contactless, full field and highly sensitive in-plane deformation measurements using the example of plywood | |
CN206990338U (zh) | 一种利用数字激光散斑法测量杨氏弹性模量的实验系统 | |
CN104713489B (zh) | 一种三维云纹干涉仪及材料表面测量方法 | |
CN105403482B (zh) | 一种相对面密度的光学算法 | |
CN106918597A (zh) | 薄膜质量检测方法和薄膜质量检测系统 | |
Safavizadeh et al. | Utilizing digital image correlation (DIC) in asphalt pavement testing | |
CN106226314B (zh) | 基于纱线截面周长的纱线条干均匀度测量方法 | |
Li et al. | Visualization of non-uniform soil deformation during triaxial testing | |
CN104122209B (zh) | 微晶硅薄膜生长过程的可视化观测系统及测量方法 | |
CN106404715B (zh) | 测量折射率的方法 | |
Chen et al. | Exploring digital image correlation technique for the analysis of the tensile properties of all-cellulose composites | |
CN102853775B (zh) | 须丛曲线的获取方法 | |
Ezazshahabi et al. | A Review on the Poisson’s Ratio of Fabrics | |
CN104238002B (zh) | 一种侧视成像保偏光纤定轴方法 | |
CN105806231A (zh) | 一种基于光纤白光干涉原理的测厚仪 | |
Aslan et al. | Characterization of individual cotton fibers via light-scattering experiments | |
CN104792733B (zh) | 一种快速定标模块及应用 | |
CN106403829B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |