发明名称 一种二值化的纳米多孔材料表征方法
摘要 一种基于二值化方法的、结合了扫描云纹技术和灰度共生矩阵处理方法的纳米多孔材料表征方法,通过扫描电子显微镜拍摄被测纳米多孔材料的云纹条纹图像,并采用灰度共生矩阵方法处理该云纹条纹图,通过对获得的云纹条纹图进行二值化处理,可将图像的灰度级别降低为原来的1/128,相应地,由二值化云纹条纹图所获得的灰度共生矩阵的维数也减少到原来的1/128。本发明方法具有快速高效等优点,同时兼具灵敏度高、操作简单,以及可批量处理等优点。
申请公布号 CN104359929A 申请公布日期 2015.02.18
申请号 CN201410673713.4 申请日期 2014.11.21
申请人 清华大学 发明人 谢惠民;李传崴;刘战伟
分类号 G01N23/22(2006.01)I 主分类号 G01N23/22(2006.01)I
代理机构 代理人
主权项 一种基于二值法的用于表征纳米多孔材料的结构主周期和结构主方向的方法,其特征在于,该方法包括以下步骤:a).将待测纳米多孔材料加工成厚度为0.1~0.3mm的矩形薄膜试样并将该试样置于扫描电镜的样品台上,调节扫描电镜的放大倍数和工作距离,以得到所述纳米多孔材料的云纹条纹,记录所述扫描电镜的放大倍数并拍摄所述云纹条纹图像;b).将所拍摄的云纹条纹图像进行二值化处理,求出所述扫描电镜二值化云纹条纹图的灰度共生矩阵;c).根据所述扫描电镜的放大倍数,确定所述扫描电镜的扫描线周期;d).将扫描电镜的扫描参数设置成:搜索区间为[1,z],搜索角度分别为0°、45°、90°和135°,搜索初始值设为1像素,分别在0°、45°、90°和135°这4个方向上以搜索初始值1个像素为搜索起点,以1个像素为单次增量逐渐增加,直至搜索位移达到z像素,停止搜索;e).根据下述公式计算所述灰度共生矩阵的对比度F<sub>1</sub>:F<sub>1</sub>=g<sub>00</sub>+g<sub>11</sub>其中,g<sub>00</sub>是所述的灰度共生矩阵中第1行、第1列的元素;g<sub>11</sub>是所述的灰度共生矩阵中第2行、第2列的元素;根据下述公式计算所述灰度共生矩阵的相关性F<sub>2</sub>:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>F</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><mo>(</mo><msub><mi>g</mi><mn>11</mn></msub><mo>-</mo><mrow><mo>(</mo><msub><mi>g</mi><mn>01</mn></msub><mo>+</mo><msub><mi>g</mi><mn>11</mn></msub><mo>)</mo></mrow><mo>&times;</mo><mrow><mo>(</mo><msub><mi>g</mi><mn>10</mn></msub><mo>+</mo><msub><mi>g</mi><mn>11</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow><mrow><mn>2</mn><msqrt><mrow><mo>(</mo><msub><mi>g</mi><mn>01</mn></msub><mo>+</mo><msub><mi>g</mi><mn>11</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>g</mi><mn>00</mn></msub><mo>+</mo><msub><mi>g</mi><mn>10</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>g</mi><mn>10</mn></msub><mo>+</mo><msub><mi>g</mi><mn>11</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>g</mi><mn>00</mn></msub><mo>+</mo><msub><mi>g</mi><mn>01</mn></msub><mo>)</mo></mrow></msqrt></mrow></mfrac><mo>+</mo><mn>0.5</mn><mo>;</mo></mrow>]]></math><img file="FSA0000110600260000011.GIF" wi="1136" he="140" /></maths>其中,g<sub>01</sub>是所述的灰度共生矩阵中第1行、第2列的元素;g<sub>10</sub>是所述的灰度共生矩阵中第2行、第1列的元素;根据下述公式计算所述灰度共生矩阵的总变化量F<sub>3</sub>:F<sub>3</sub>=(g<sub>01</sub>+g<sub>10</sub>+2g<sub>11</sub>)‑(g<sub>01</sub>+g<sub>10</sub>+2g<sub>11</sub>)<sup>2</sup>+2g<sub>11</sub>;根据下述公式计算所述灰度共生矩阵的差熵F<sub>4</sub>:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>F</mi><mtext>4</mtext></msub><mo>=</mo><mfrac><mrow><mn>0.75</mn><mo>-</mo><mn>2</mn><mrow><mo>(</mo><msub><mi>g</mi><mn>00</mn></msub><mo>+</mo><msub><mi>g</mi><mn>11</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>g</mi><mn>10</mn></msub><mo>+</mo><msub><mi>g</mi><mn>01</mn></msub><mo>)</mo></mrow></mrow><mn>0.75</mn></mfrac><mo>;</mo></mrow>]]></math><img file="FSA0000110600260000012.GIF" wi="726" he="115" /></maths>根据下述公式计算所述灰度共生矩阵的κ统计量F<sub>5</sub>:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>F</mi><mn>5</mn></msub><mo>=</mo><mfrac><mrow><mn>2</mn><mo>-</mo><mrow><mo>(</mo><msub><mi>g</mi><mn>01</mn></msub><mo>+</mo><msub><mi>g</mi><mn>10</mn></msub><mo>)</mo></mrow><mo>/</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mrow><mo>(</mo><msub><mi>g</mi><mn>00</mn></msub><mo>+</mo><msub><mi>g</mi><mn>01</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>g</mi><mn>10</mn></msub><mo>+</mo><msub><mi>g</mi><mn>11</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>g</mi><mn>01</mn></msub><mo>+</mo><msub><mi>g</mi><mn>11</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>g</mi><mn>00</mn></msub><mo>+</mo><msub><mi>g</mi><mn>10</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><mo>;</mo></mrow>]]></math><img file="FSA0000110600260000013.GIF" wi="1331" he="112" /></maths>根据下述公式计算所述灰度共生矩阵的对比度F<sub>1</sub>、相关性F<sub>2</sub>、总变化量F<sub>3</sub>、差熵F<sub>4</sub>和κ统计量F<sub>5</sub>的算数平均值<img file="FSA0000110600260000021.GIF" wi="65" he="70" /><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mover><mi>F</mi><mo>~</mo></mover><mo>=</mo><mfrac><mn>1</mn><mn>5</mn></mfrac><mi>&Sigma;</mi><msub><mi>F</mi><mi>i</mi></msub><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2,3,4,5</mn><mo>;</mo></mrow>]]></math><img file="FSA0000110600260000022.GIF" wi="520" he="110" /></maths>f).分别计算前述d)中提到的搜索区间[1,z]内的所有搜索位移的<img file="FSA0000110600260000023.GIF" wi="38" he="61" />值,当<img file="FSA0000110600260000024.GIF" wi="37" he="61" />取最大值时的搜索位移即为所述扫描电镜云纹条纹二值化图的结构主周期d,此时的方向即为结构主方向θ;g).根据下述公式计算被测纳米多孔材料试样的结构主周期p′:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msup><mi>p</mi><mo>&prime;</mo></msup><mo>=</mo><mfrac><mi>d</mi><msqrt><mn>1</mn><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mi>d</mi><mi>p</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>2</mn><mfrac><mi>d</mi><mi>p</mi></mfrac><mi>cos</mi><mi>&phi;</mi></msqrt></mfrac><mo>;</mo></mrow>]]></math><img file="FSA0000110600260000025.GIF" wi="517" he="181" /></maths>根据下述公式计算被测纳米多孔材料试样的结构主方向sinθ:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>sin</mi><mi>&theta;</mi><mo>=</mo><mfrac><mrow><mi>sin</mi><mi>&phi;</mi></mrow><msqrt><msup><mi>sin</mi><mn>2</mn></msup><mi>&phi;</mi><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mi>d</mi><mi>p</mi></mfrac><mo>+</mo><mi>cos</mi><mi>&phi;</mi><mo>)</mo></mrow><mn>2</mn></msup><mo></mo></msqrt></mfrac><mo>.</mo></mrow>]]></math><img file="FSA0000110600260000026.GIF" wi="592" he="181" /></maths>
地址 100084 北京市海淀区清华大学