发明名称 铁矿石粗糙度自动数字量化测量方法
摘要 本发明为一种铁矿石粗糙度自动数字量化测量方法,采用扫描电子显微镜和计算机图像处理系统进行铁矿石粗糙度自动数字化测量。本发明的方法通过计算机系统对图像进行处理得到的数据完全符合人的量化特征,表明粗糙度的量化完全可以用数字代替,避免了人为的将铁矿石粗糙度粗分为有限的几个量化等级;并且避免了人工判断的主观性,测量准确、快速,操作简便,为烧结制粒提供配料依据;实现了测量过程的自动化,无须专业人员操作。
申请公布号 CN103267498B 申请公布日期 2016.04.27
申请号 CN201310167762.6 申请日期 2013.05.09
申请人 北京科技大学 发明人 国宏伟;苏步新;张建良;白真龙;李新宇
分类号 G01B11/30(2006.01)I 主分类号 G01B11/30(2006.01)I
代理机构 北京金智普华知识产权代理有限公司 11401 代理人 巴晓艳
主权项 一种铁矿石粗糙度自动数字量化测量方法,其特征在于,所述方法采用图像采集系统和计算机图像处理系统,具体方法包括以下步骤:1)不同种类铁矿石图像采集:利用扫描电子显微镜对不同种类的铁矿石进行图像采集;2)图像预处理:对不同种类的铁矿石分别进行图像预处理,首先对得到的图像进行分段灰度拉升,在灰度拉升的基础上,对图像背景进行滤除;3)局部灰度面平滑特征处理:对预处理后的图像进行局部灰度面平滑特征处理,得到的灰度平滑像素点用红色标记,计算总红色像素点占总有效像素点的比值作为量化铁矿石的粗糙度的指标,比值越小,说明铁矿石的粗糙度越大;所述步骤2)具体是:<img file="FDA0000869725260000013.GIF" wi="54" he="52" />分段灰度拉升在灰度拉升操作中,通过对原始图像中每个像素点赋予一个新的灰度值来增强图像的视觉效果,如式(1)所示:t=f(s)        (1)其中,s为原始图像某一位置的灰度值,t为增强后输出图像同一位置的灰度值,对于目标区域灰度动态范围偏窄的图像,则采用分段线性变换如式(2)所示:<maths num="0001"><math><![CDATA[<mrow><mi>t</mi><mo>=</mo><mo>{</mo><mrow><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mn>0</mn><mo>&lt;</mo><mi>s</mi><mo>&le;</mo><mi>b</mi></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mn>255</mn><mrow><mi>c</mi><mo>-</mo><mi>b</mi></mrow></mfrac><mrow><mo>(</mo><mi>s</mi><mo>-</mo><mi>b</mi><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mi>b</mi><mo>&lt;</mo><mi>s</mi><mo>&le;</mo><mi>c</mi></mrow></mtd></mtr><mtr><mtd><mn>255</mn></mtd><mtd><mrow><mi>c</mi><mo>&lt;</mo><mi>s</mi><mo>&le;</mo><mn>255</mn></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo></mrow><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000869725260000011.GIF" wi="1501" he="319" /></maths>其中b和c是整个映射规则的分段点,通过灰度归一化算法进行处理取值b=0,c=255;<img file="FDA0000869725260000014.GIF" wi="54" he="52" />背景滤除利用式(3)对图像中低灰度的黑色部分滤除:<img file="FDA0000869725260000012.GIF" wi="1158" he="158" />其中g(x,y)为某点的灰度值,T为阈值,小于T的黑色背景区域将被滤除;所述步骤3)具体为:对于图像中有效灰度像素点g(x,y),扫描其周围的以正方形区(x‑r,y‑r)到(x+r,y+r)的所有点,并把其周围所有点的灰度像素值与g(x,y)相减,求差值的绝对平均,如公式(4)所示:<maths num="0002"><math><![CDATA[<mrow><mi>s</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>k</mi></mfrac><mo>&CenterDot;</mo><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mo>-</mo><mi>r</mi></mrow><mi>r</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mo>-</mo><mi>r</mi></mrow><mi>r</mi></munderover><mo>|</mo><mi>g</mi><mrow><mo>(</mo><mi>x</mi><mo>+</mo><mi>m</mi><mo>,</mo><mi>y</mi><mo>+</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mi>g</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000869725260000021.GIF" wi="1175" he="135" /></maths>其中g(x+m,y+n)为正方形区域(x‑r,y‑r)到(x+r,y+r)的一个有效点,m、n分别为行数和列数,k为区域(x‑r,y‑r)到(x+r,y+r)有效像素点的总数,s(x,y)为差值的绝对平均,平滑点的识别则由公式(5)给出:<img file="FDA0000869725260000022.GIF" wi="1206" he="159" />其中T<sub>0</sub>为给定的特定灰度值,在图像上用红色进行标记平滑点,计算总红色像素点占总有效像素点的比值,比值越大表示平滑度越高,粗糙度越小。
地址 100083 北京市海淀区学院路30号