发明名称 矿相成份的机器测量方法
摘要 本发明公开一种矿相成份的机器测量方法,采用矿相显微镜、摄像机、图像采集卡和计算机系统进行机器测量,包括如下步骤:1)由计算机建立矿物灰度分布特征数据库;2)制样;3)图像采集;4)计算机计算:将第3)步获取的显微图像输入计算机,利用第1)步建立的矿物灰度分布特征数据库中的各种矿物组合对待测矿物的灰度直方图曲线进行拟合计算,并使用遗传算法优化求解,最终得到待测矿物中不同组成的含量。本发明由于采用计算机系统进行智能测算,代替了实验员的工作,降低了劳动强度;对于含量较少的矿物敏感,测量准确、快速,操作简便;实现了测量过程的自动化,无需专业人员操作。
申请公布号 CN101162202B 申请公布日期 2010.05.26
申请号 CN200710093062.1 申请日期 2007.11.28
申请人 重庆大学 发明人 邱贵宝;吕学伟;白晨光;张生富;梁栋;扈玫珑;施瑞盟
分类号 G01N21/84(2006.01)I;G01N21/00(2006.01)I;G06F17/11(2006.01)I 主分类号 G01N21/84(2006.01)I
代理机构 重庆博凯知识产权代理有限公司 50212 代理人 张先芸
主权项 矿相成份的机器测量方法,采用矿相显微镜、摄像机、图像采集卡和计算机系统进行机器测量,其特征在于包括如下步骤:1)由计算机建立矿物灰度分布特征数据库:根据常见矿物标准样本的显微图像灰度直方图和高斯分布模型统计得到这些矿物在直方图中的灰度分布规律;2)制样:将待测矿物经过粗磨、细磨和抛光后备用;3)图像采集:由矿相显微镜和摄像机对上述第2)步获取的试样进行图像采集,获取矿物的显微图像;4)计算机计算:将第3)步获取的显微图像输入计算机,利用第1)步建立的矿物灰度分布特征数据库中的各种矿物组合对待测矿物的灰度直方图曲线进行拟合计算,并使用遗传算法优化求解,最终得到待测矿物中不同组成的含量;所述第4)步使用遗传算法计算矿物含量,包括:(1)灰度直方图:灰度直方图是图像的重要特征,是灰度分布密度函数的近似,它表示了在图像区域所有像素中,不同的灰度值出现的次数;其函数表达式如式(5)所示: <mrow> <msub> <mi>G</mi> <mi>i</mi> </msub> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <munderover> <mi>&Sigma;</mi> <mrow> <mi>n</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mi>P</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>其中,M、N分别表示图像的行数和列数,Gi为灰度为i的像素个数;P(i,m,n)的定义如式(6)所示: <mrow> <mi>P</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mi>P</mi> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>i</mi> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mi>P</mi> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>&NotEqual;</mo> <mi>i</mi> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>在灰度直方图中,每种矿物的灰度值在其固有的区间内呈正态分布,即如式(1)所不:Gm~N(μ,σ2) (1)其中,Gm为矿物m的灰度分布,μ、σ2分别是正态分布的均值和方差;对式(1)进行归一化处理,得到了式(7)所示的归一化直方图函数;其中,Gi′为归一化值,T为总像素数:Gi′=Gi/T (7)(2)几种矿物高斯分布迭加矿物服从如式(8)所示的正态分布规律,其中Gx为第x种矿物的灰度分布,μx是待测矿物x在灰度直方图中的灰度均值,σx2为第x种矿物在灰度直方图中的灰度分布方差,σx即为所述灰度分布方差的平方根;式(9)为第x种矿物灰度正态分布的密度函数,其中Cx为第x种矿物的含量;将不同矿物相加就得到了如式(10)所示的几种矿物迭加的密度分布函数,其中X为矿物的种类总数;Gx~Nx(μx,σx2) (8) <mrow> <msub> <mi>f</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <msub> <mi>C</mi> <mi>x</mi> </msub> <mrow> <msub> <mi>&sigma;</mi> <mi>x</mi> </msub> <msqrt> <mn>2</mn> <mi>&pi;</mi> </msqrt> </mrow> </mfrac> <msup> <mi>e</mi> <mfrac> <msup> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <msub> <mi>&mu;</mi> <mi>x</mi> </msub> <mo>)</mo> </mrow> </mrow> <mn>2</mn> </msup> <msup> <msub> <mrow> <mn>2</mn> <mi>&sigma;</mi> </mrow> <mi>x</mi> </msub> <mn>2</mn> </msup> </mfrac> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>x</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>X</mi> </munderover> <mfrac> <msub> <mi>C</mi> <mi>x</mi> </msub> <mrow> <msub> <mi>&sigma;</mi> <mi>x</mi> </msub> <msqrt> <mn>2</mn> <mi>&pi;</mi> </msqrt> </mrow> </mfrac> <msup> <mi>e</mi> <mfrac> <msup> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <msub> <mi>&mu;</mi> <mi>x</mi> </msub> <mo>)</mo> </mrow> </mrow> <mn>2</mn> </msup> <msup> <msub> <mrow> <mn>2</mn> <mi>&sigma;</mi> </mrow> <mi>x</mi> </msub> <mn>2</mn> </msup> </mfrac> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>(3)适应度函数采用高斯分布合成得到的曲线与实际灰度直方图曲线的差作为表征个体适应度的大小,其表达式如式11所示; <mrow> <mi>d</mi> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>i</mi> <mo>=</mo> <mn>256</mn> </mrow> </munderover> <mo>|</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>-</mo> <msup> <msub> <mi>G</mi> <mi>i</mi> </msub> <mo>&prime;</mo> </msup> <mo>|</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>变量为二进制编码,设置每个变量的编码长度;确定选择操作的方式,设置进化代数,选择概率和遗传概率参数;(4)求解使用遗传算法优化求解,最终得到使式(11)中的d最小时的各种矿物的含量多少。
地址 400030 重庆市沙坪坝区沙正街174号
您可能感兴趣的专利