CN103808646A - 测量颗粒材料线分形维数的图像分析仪及其方法 - Google Patents
测量颗粒材料线分形维数的图像分析仪及其方法 Download PDFInfo
- Publication number
- CN103808646A CN103808646A CN201410044578.7A CN201410044578A CN103808646A CN 103808646 A CN103808646 A CN 103808646A CN 201410044578 A CN201410044578 A CN 201410044578A CN 103808646 A CN103808646 A CN 103808646A
- Authority
- CN
- China
- Prior art keywords
- fractal dimension
- image
- line
- particle
- particle image
- 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.)
- Granted
Links
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明提供一种测量颗粒材料线分形维数的图像分析仪及其方法,所要解决的问题是:现有的图像分析仪仅能获得颗粒材料二维图像数据,在对二维图像数据进行处理之后,应用上述幂指数与分形维数关系回归分析才能得到的颗粒材料的线分维数,而且回归分析系数相对较低、无法直接读取分形维数结果;本发明的技术要点是:数据处理系统包括测控模块、颗粒图像处理模块和线分形维数计算模块;以及配合图像分析仪使用的方法;本发明通过在颗粒图像分析系统的测控软件中增加颗粒图像处理模块和线分形维数计算模块,使颗粒图像分析系统在现有功能基础上,增加了颗粒二维图像轮廓线分形维数测量新功能,并能够直接输出颗粒材料二维图像线分形维数的测量结果。
Description
技术领域
本发明涉及颗粒材料检测技术领域,具体说是一种测量颗粒材料二维图像轮廓线分形维数的测量颗粒材料线分形维数的图像分析仪及其方法。
背景技术
分形是部分与整体以某种方式相似的形体。分形维数(又称分数维)在表征不规则形状和传统方法难以表征的材料性质(如颗粒的团聚、针状材料和片状材料)方面正在发挥越来越重要的作用。大量实验结果表明:颗粒材料二维图像轮廓满足统计自相似,属于无规分形。
目前,颗粒材料(特别是粉体材料)二维图像轮廓分形维数定义方法大多采用改变粗视化程度求维数的方法如:码尺法、盒维数法。其共同点是:1)分形满足幂律,且幂指数与分形维数的关系是线性关系,其中码尺法为: ,盒维数法为:;2)通过线性回归分析计算得到相近材料分形维数的结果波动很大;3)图像分析仪只能获得颗粒材料的二维图像及其他表征量结果,不能对二维图像轮廓进行计算与分析从而得到线分形维数。
现有的图像分析仪仅能获得颗粒材料二维图像数据,在对二维图像数据进行处理之后,应用上述幂指数与分形维数关系回归分析才能得到的颗粒材料的线分维数,而且回归分析系数相对较低、无法直接读取分形维数结果。
发明内容
本发明的目的是提供一种测量颗粒材料二维图像轮廓线的测量颗粒材料线分形维数的图像分析仪及其方法。
本发明的技术方案是:光学显微镜与摄像机和数字照相机连接,摄像机上设置有图像采集卡,摄像机和数字照相机与数据处理系统连接;所述摄像机镜头内设置有光源;
所述数据处理系统包括原测控软件模块、颗粒图像处理模块和线分形维数计算模块;所述颗粒图像处理模块包括:
颗粒图像数据接收单元,用于接收一次图像处理结果;
颗粒图像二次处理单元,用于对一次图像处理结果进行筛选及二次处理,按顺序储存;
所述线分形维数计算模块包括:
单颗粒图像最大弦直径计算单元,用于计算每个二次处理后颗粒图像的最大弦直径,按顺序储存;
单颗粒图像线分形维数计算单元,用于计算每个颗粒图像的线分形维数,按顺序储存;
颗粒群线分形维数计算单元,用于计算颗粒群的线分形维数;
线分形维数储存和/或发送单元,用于储存和/或发送颗粒群线分形维数值。
测量颗粒材料线分形维数的方法,采用下列步骤:
1)、将颗粒材料图像二次处理程序及线分形维数计算程序按顺序储存导入数据处理系统的程序储存器中;
2)、启动该计算机系统,CPU对通过光学显微镜获取并经过测控模块一次处理的颗粒材料二维图像数据进行读取,并选择指令、执行操作;
在选择分析指令时,执行颗粒图像处理程序,对一次处理后的颗粒材料二值化图片进行数学形态化运算,获得单颗粒图像轮廓坐标;
在选择储存指令时,对单颗粒图像轮廓坐标按顺序储存。
3)、使用经过二次处理的单颗粒图像轮廓坐标数据进行建模:建模计算公式如下:
4)、使用经过归一化处理的数据进行建模,建模计算公式如下:
6)、对每个经过二次处理后的颗粒图像轮廓坐标,重复步骤3)、步骤4)和步骤5),逐一计算单颗粒材料二维图像轮廓的线分形维数,按顺序储存;
光学显微镜带有至少两组成像系统,一组设置在竖直与台座的位置,另一组设置在其水平偏转+60o的位置上.
本发明的有益效果是:通过分形理论建立单颗粒材料二维图像轮廓分形数学模型,得到了描述轮廓分形自相似规律的幂指数与颗粒材料线分形维数之间的非线性关系为:,解决了由大量颗粒组成的颗粒群材料二维图像轮廓线分形维数自动测量的难题,克服了相近材料线分形维数有波动的不足;颗粒材料二维图像数据的二次处理有效提高了回归分析的相关系数;本发明适用于颗粒材料二维图像轮廓线分形维数表征,同样也适用于其它领域满足自相似规律二维轮廓曲线的线分形维数表征;本发明通过在颗粒图像分析系统的测控软件中增加颗粒图像处理模块和线分形维数计算模块,使颗粒图像分析系统在现有功能基础上,增加了颗粒二维图像轮廓线分形维数测量新功能,并能够直接输出颗粒材料二维图像线分形维数的测量结果。
附图说明
图1是本发明模块结构示意图;
图2是颗粒处理前示意图;
图3是颗粒经一次处理后示意图;
图4是颗粒轮廓点坐标示意图;
图5是图4的A部放大图。
下面将通过实例对本发明作进一步详细说明,但下述的实例仅仅是本发明其中的例子而已,并不代表本发明所限定的权利保护范围。
具体实施方式
参见图1,光学显微镜1与摄像机2和数字照相机3连接,摄像机与图像采集卡4连接,摄像机和数字照相机与数据处理系统5连接;所述摄像机镜头内设置有光源;
所述数据处理系统包括原测控软件模块51、颗粒图像处理模块52和线分形维数计算模块53;所述颗粒图像处理模块包括:
颗粒图像数据接收单元,用于接收一次图像处理结果;
颗粒图像二次处理单元,用于对一次图像处理结果进行筛选及二次处理,按顺序储存;
所述线分形维数计算模块包括:
单颗粒图像最大弦直径计算单元,用于计算每个二次处理后颗粒图像的最大弦直径,按顺序储存;
单颗粒图像线分形维数计算单元,用于计算每个颗粒图像的线分形维数,按顺序储存;
颗粒群线分形维数计算单元,用于计算颗粒群的线分形维数;
线分形维数储存和/或发送单元,用于储存和/或发送颗粒群线分形维数值。
光学显微镜带有至少两组成像系统,一组设置在竖直与台座的位置,另一组设置在其水平偏转+60o的位置上.
测量颗粒材料线分形维数的方法采用下列步骤:
1)、将颗粒材料图像二次处理程序及线分形维数计算程序导入数据处理系统的程序储存器中;
2)、启动该计算机系统,CPU对通过光学显微镜获取并经过原测控软件模块一次处理的颗粒材料二维图像数据进行读取,并选择指令、执行操作;
在选择分析指令时,执行颗粒图像处理程序,对一次处理后的颗粒材料二值化图片进行数学形态化运算,获得单颗粒图像轮廓坐标;
在选择储存指令时,对单颗粒图像轮廓坐标按顺序储存。
3)、使用经过二次处理的单颗粒图像轮廓坐标数据进行建模:建模计算公式如下:
建立颗粒图像轮廓步长和周长归一化计算模型,具体计算过程为:首先计算颗粒图像最大弦直径 ,其中是任意两个不同的颗粒图像轮廓坐标点之间的距离,即,所述是一个维的颗粒图像轮廓坐标矩阵,且,;将结果储存到数据处理系统的存储器中;其次是选取步长和计算,步长分别选取5、10、15、20、25、30、35、40、45、50个像素点为步长,并记为,与不同步长相对应的轮廓度量周长的步数为,其中步数分别选取23、19、9、7、6、5、4、3、1、1,相对应的度量周长为;最后是步长的归一化处理,用最大弦直径对步长归一化处理的结果是0.0937、0.1875、0.2812、0.375、0.4687、0.5624、0.6562、0.7499、0.8437、0.9374;
4)、使用经过归一化处理的数据进行建模,建模计算公式如下:
建立单颗粒材料二维图像轮廓分形理论模型,并得到线分形维数与归一化后的步长和步数的关系为:,其中:是单颗粒材料二维图像轮廓线分形维数,为幂指数;该模型中幂指数与线分形维数的关系为非线性关系,对上述关系两边取对数得:;
采用最小二乘法回归分析计算幂指数,即双对数坐标曲线的斜率为-1.2655;
6)、对每个经过二次处理后的颗粒图像轮廓坐标,重复步骤3)、步骤4)和步骤5),逐一计算单颗粒材料二维图像轮廓的线分形维数,按顺序储存;
本发明的图像处理方法如下:经原系统图像处理得到的单颗粒图像,颗粒边缘存在模糊不清的像素点,参见图2;
图像一次处理,消除模糊不清的像素,参见图3;图像二次处理,即数学形态化运算,计算单个颗粒图像边缘轮廓点的坐标,参见图4;
读取步骤3处理所得单个颗粒图像轮廓边缘点的坐标,“x”号代表颗粒边缘轮廓坐标点,分别选取任一轮廓边缘坐标点为起始点,计算与其他各轮廓边缘坐标点之间距离,即:,然后依次选取剩余坐标点为起始点,进行所述计算,得到距离数组阵,选取大弦直径,即:;
Claims (2)
1.一种测量颗粒材料线分形维数的图像分析仪,它包括摄像机,其特征是:光学显微镜与摄像机和数字照相机连接,摄像机上设置有图像采集卡,摄像机和数字照相机与数据处理系统连接;所述摄像机镜头内设置有光源;
所述数据处理系统包括测控模块、颗粒图像处理模块和线分形维数计算模块;所述颗粒图像处理模块包括:
颗粒图像数据接收单元,用于接收一次图像处理结果;
颗粒图像二次处理单元,用于对一次图像处理结果进行筛选及二次处理,按顺序储存;
所述线分形维数计算模块包括:
单颗粒图像最大弦直径计算单元,用于计算每个二次处理后颗粒图像的最大弦直径,按顺序储存;
单颗粒图像线分形维数计算单元,用于计算每个颗粒图像的线分形维数,按顺序储存;
颗粒群线分形维数计算单元,用于计算颗粒群的线分形维数;
线分形维数储存和/或发送单元,用于储存和/或发送颗粒群线分形维数值。
2.一种测量颗粒材料线分形维数的方法,其特征在于采用下列步骤:
1)、将颗粒材料图像二次处理程序及线分形维数计算程序按顺序储存导入数据处理系统的程序储存器中;
2)、启动该计算机系统,CPU对通过光学显微镜获取并经过测控模块一次处理的颗粒材料二维图像数据进行读取,并选择指令、执行操作;
在选择分析指令时,执行颗粒图像处理程序,对一次处理后的颗粒材料二值化图片进行数学形态化运算,获得单颗粒图像轮廓坐标;
在选择储存指令时,对单颗粒图像轮廓坐标按顺序储存;
3)、使用经过二次处理的单颗粒图像轮廓坐标数据进行建模:建模计算公式如下:
4)、使用经过归一化处理的数据进行建模,建模计算公式如下:
6)、对每个经过二次处理后的颗粒图像轮廓坐标,重复步骤3)、步骤4)和步骤5),逐一计算单颗粒材料二维图像轮廓的线分形维数,按顺序储存;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410044578.7A CN103808646B (zh) | 2014-02-07 | 2014-02-07 | 测量颗粒材料线分形维数的图像分析仪及其方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410044578.7A CN103808646B (zh) | 2014-02-07 | 2014-02-07 | 测量颗粒材料线分形维数的图像分析仪及其方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103808646A true CN103808646A (zh) | 2014-05-21 |
CN103808646B CN103808646B (zh) | 2016-04-06 |
Family
ID=50705737
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410044578.7A Active CN103808646B (zh) | 2014-02-07 | 2014-02-07 | 测量颗粒材料线分形维数的图像分析仪及其方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103808646B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2574173C1 (ru) * | 2014-10-29 | 2016-02-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Южно-Уральский государственный университет" (национальный исследовательский университет) (ФГБОУ ВПО "ЮУрГУ" (НИУ)) | Способ определения фрактальной размерности границ зерен формовочного песка |
CN105894001A (zh) * | 2016-03-30 | 2016-08-24 | 中国矿业大学 | 一种磨屑轮廓自动提取与分形识别系统及应用方法 |
CN106596357A (zh) * | 2016-11-28 | 2017-04-26 | 江苏大学 | 一种柴油机颗粒物形貌的表征方法 |
CN113776469A (zh) * | 2021-08-10 | 2021-12-10 | 同济大学 | 一种粉体颗粒表面粗糙度的检测方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5848177A (en) * | 1994-12-29 | 1998-12-08 | Board Of Trustees Operating Michigan State University | Method and system for detection of biological materials using fractal dimensions |
JPH11306362A (ja) * | 1998-04-15 | 1999-11-05 | Sysmex Corp | 画像解析装置 |
CN101655913A (zh) * | 2009-09-17 | 2010-02-24 | 上海交通大学 | 分形维数的计算机生成图像被动检测方法 |
CN102800113A (zh) * | 2012-07-18 | 2012-11-28 | 合肥工业大学 | 一种基于分形维数的数字图像分析方法 |
-
2014
- 2014-02-07 CN CN201410044578.7A patent/CN103808646B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5848177A (en) * | 1994-12-29 | 1998-12-08 | Board Of Trustees Operating Michigan State University | Method and system for detection of biological materials using fractal dimensions |
JPH11306362A (ja) * | 1998-04-15 | 1999-11-05 | Sysmex Corp | 画像解析装置 |
CN101655913A (zh) * | 2009-09-17 | 2010-02-24 | 上海交通大学 | 分形维数的计算机生成图像被动检测方法 |
CN102800113A (zh) * | 2012-07-18 | 2012-11-28 | 合肥工业大学 | 一种基于分形维数的数字图像分析方法 |
Non-Patent Citations (1)
Title |
---|
宋延勇 等: "颗粒显微图像分析程序设计", 《仪器仪表学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2574173C1 (ru) * | 2014-10-29 | 2016-02-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Южно-Уральский государственный университет" (национальный исследовательский университет) (ФГБОУ ВПО "ЮУрГУ" (НИУ)) | Способ определения фрактальной размерности границ зерен формовочного песка |
CN105894001A (zh) * | 2016-03-30 | 2016-08-24 | 中国矿业大学 | 一种磨屑轮廓自动提取与分形识别系统及应用方法 |
CN106596357A (zh) * | 2016-11-28 | 2017-04-26 | 江苏大学 | 一种柴油机颗粒物形貌的表征方法 |
CN113776469A (zh) * | 2021-08-10 | 2021-12-10 | 同济大学 | 一种粉体颗粒表面粗糙度的检测方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103808646B (zh) | 2016-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109545072B (zh) | 地图构建的位姿计算方法、装置、存储介质和系统 | |
CN104408055B (zh) | 一种激光雷达点云数据的存储方法及装置 | |
CN102831427A (zh) | 一种融合视觉显著性和灰度共生矩的纹理特征提取方法 | |
CN115541030B (zh) | 高炉炉顶料面温度分布的识别方法、装置及存储介质 | |
CN101599181A (zh) | 一种代数b样条曲面的实时绘制方法 | |
CN111275821B (zh) | 一种电力线拟合方法、系统及终端 | |
Zhou et al. | A streaming framework for seamless building reconstruction from large-scale aerial lidar data | |
CN103808646B (zh) | 测量颗粒材料线分形维数的图像分析仪及其方法 | |
CN111709999A (zh) | 标定板、相机标定方法、装置、电子设备及相机系统 | |
CN104299241A (zh) | 基于 Hadoop 的遥感图像显著性目标检测方法及系统 | |
CN108596032B (zh) | 一种视频中打架行为检测方法、装置、设备及介质 | |
CN111161348B (zh) | 一种基于单目相机的物体位姿估计方法、装置及设备 | |
Zhang et al. | Fully automatic system for fish biomass estimation based on deep neural network | |
CN116266365A (zh) | 一种点云数据提取方法、装置、电子设备及存储介质 | |
CN109238264A (zh) | 一种家畜位置姿势归一化方法及装置 | |
CN117058022A (zh) | 深度图像去噪方法、装置、计算机设备及存储介质 | |
CN114723808B (zh) | 一种基于fpga的光条亚像素中心提取方法 | |
CN116824371A (zh) | 生态保护红线植被异常快速监测方法和装置 | |
CN113920269B (zh) | 工程进度获取方法、装置、电子设备及介质 | |
Mo et al. | Automatic power transmission towers detection based on the deep learning algorithm | |
CN110322415A (zh) | 基于点云的高精度表面三维重构方法 | |
CN103954536A (zh) | 测量颗粒材料体分形维数的激光粒度仪及其方法 | |
CN108278988A (zh) | 测量叶片面积的方法以及装置 | |
CN109658489B (zh) | 一种基于神经网络的立体网格数据处理方法和系统 | |
CN103530647A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |