专利名称：一种多孔材料孔隙边界提取和数值化方法
技术领域：
本发明涉及一种孔隙边界的提取和描述方法，尤其是涉及一种多孔材 料孔隙边界提取和数值化方法。
背景技术：
多孔材料广泛存在于自然界和工程应用中，是一类重要的结构和功能 一体化材料，多孔材料内部所存在的大量孔隙使得其具有诸多不同于致密 材料的优异性能。多孔材料的应用性能主要取决于孔结构特征，孔结构参 数是多孔材料最重要的特征参量，对于孔结构的准确描述是多孔材料研 究、制造和应用的基础问题。但是由于孔隙大小不一、随机分布，孔隙空 间结构的杂乱无序，使得多孔材料孔结构的科学描述变得异常复杂。
常规多孔材料孔结构的描述方法是基于毛细管假设下的物理性能检 测，其将孔隙和孔通道分别简化为圆和直管，即所谓的"圆-直管模型"， 用不同的流体压差-流量数据计算得到等效孔径和孔径分布，这是一种近 似描述，而且由于检测方法的不同，结果差异较大，也难以反映孔隙形貌、 孔通道的弯曲程度等特征。

发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一 种多孔材料孔隙边界提取和数值化方法，其搡作简单且成本低，能对多孔 材料的孔隙形貌和孔通道的弯曲程度等性能参数进行较准确分析。
为解决上述技术问题，本发明采用的技术方案是 一种多孔材料孔隙边 界提取和数值化方法，其特征在于该方法包括以下步骤
步骤一、釆用图像釆集设备摄取被检测多孔材料的多个相互平行的断面的微观孔隙结构图像，其中，所述多孔材料的实体骨架材料和孔隙空间
分别对应的图像区域之间具有明显的色彩对比；
步骤二、将所述多个孔隙结构图像分别传输至计算机进行预处理，其 预处理过程包括以下步骤
(a) 将所述孔隙结构图像转化为多灰度级灰度图像；
(b) 对所述多灰度级灰度图像进行去噪和滤波处理；
(c) 将经去噪和滤波处理的所述多灰度级灰度图像转换为二灰度级 图像即黑白图像；
所述黑白图像由像素值为1的白色图像区域和像素值为0的黑色图像 区域组成，其中，所述实体骨架材料对应的图像区域为白色图像区域或黑 色图像区域，所述孔隙空间对应的图像区域相应为黑色图像区域或白色图 像区域；
步骤三、用计算机提取所述黑白图像中黑色图像区域和白色图像区域 的交界线即孔隙空间边界并用彩色线条进行标示，所标示出来的孔隙空间 边界为由多个像素点为顶点组成的多边形，确定所述平面多边形所有顶点 的二维坐标值并对其各顶点进行连续编号，即获得被检测多孔材料的多个 相互平行断面的二维孔隙边界图像；
步骤四、用计算机将相邻两个所述二维孔隙边界图像中所有平面多边 形的各顶点按编号顺序对应依次连接后，所组成的空间多边形即为被检测 多孔材料的三维孔隙边界。
步骤一中所述的图像釆集设备为金相显微镜或扫描电子显微镜。
步骤一中所述的孔隙结构图像为分辨率不低于500 x 500的清晰彩色 或多灰度级灰度图像。
步骤(c)中所述的滤波处理为最小值滤波和/或中值滤波和/或最大 值滤波和/或均值滤波和/或低通滤波和/或高通滤波处理。
步骤(c)中所述的实体骨架材料对应的图像区域为白色图像区域， 所述孔隙空间对应的图像区域相应为黑色图像区域。步骤三中还需对所述平面多边形进行简化处理，所述简化处理的方法 为减少所述平面多边形的边数并相应增加各边边长。
步骤三中确定所述平面多边形所有顶点的二维坐标值时，选取所述黑 白图像左下方的顶点为坐标原点。
步骤三中对所述平面多边形各顶点进行连续编号时，以二维坐标绝对 值最小的顶点为1 #顶点，然后按照逆时针顺序依次对剩余顶点进行连续编 号。
步骤一中所述相互平行断面的数量不少于三个。 步骤一中所述相互平行断面的数量为三个。
本发明与现有技术相比具有以下优点，1、不但可得到数字化的二维
孔隙边界，还可得到数字化的三维孔隙边界；2、所得到的孔隙边界即孔 隙空间边界反映了材料的真实孔隙结构，同时具有操作简单、成本低等优 点；3、通过数学分析方法和计算机数据处理技术，再结合孔隙边界数值 化后的大量结果数据，对多孔材料结构进行系统分析，便可获得多孔材料 的孔径、孔径分布、孔隙面积、孔隙周长等结构性能参数，可利用分形理 论对多孔材料进行数学分析，进而建立孔隙边界特征与材料应用性能间的 关系，以便于指导多孔材料的制造和应用；4、适用范围广，不仅可以对 已知的材质，例如金属、陶瓷、高分子、岩土等制成多孔材料的表面或剖 面孔隙结构进行分析，也可以对骨骼、木材等具有多孔结构的材料的表面 或剖面孔隙结构进行分析。
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。


图l为本发明所获取的二维孔隙边界多边形结构示意图。 图2为本发明所形成三个平行剖面的结构示意图。 图3为本发明重建形成的三维孔隙边界结构示意图。 附图标记说明1 —剖面一； 2 —剖面二； 3 —剖面三；
4 —简化多边形e
具体实施例方式
本发明所述的多孔材料孔隙边界的提取及其数值化方法包括以下步
骤
步骤一、采用图像釆集设备摄取被检测多孔材料的多个相互平行断面 的微观孔隙结构图像，其中，所述多孔材料的实体骨架材料和孔隙空间分 别对应的图像区域之间具有明显的色彩对比。其中，所述断面为多孔材料 的外表面和/或剖面。所述图像采集设备为金相显微镜或扫描电子显微镜， 并且所述相互平行断面的数量不少于三个。实践中，也可以釆用其它光学 或电子的材料结构图像釆集设备获取上述微观孔隙结构图像。所检测的多 孔材料可以是已知的材质，例如金属、陶瓷、高分子、岩土等制成多孔材 料的表面或剖面微观结构图像，也可以是骨骼、木材等具有多孔结构的材 料的表面或剖面微观结构图像。本实施例中，所采用的图像采集设备为金 相显微镜，被检测的多孔材料为烧结金属粉末多孔材料。并且，为保证所 提取孔隙边界的精度，所获取的多孔材料的孔隙结构图像应具有较高的清 晰度和分辨率，本实施例中，所述孔隙结构图像为分辨率不低于500 x 500 的清晰彩色图像或多灰度级灰度图像。
另外，结合图2,本实施例中，所述相互平行断面的数量为三个，并 且为被检测多孔材料的三个相互平行的剖面。具体的提取过程是首先准 备一块待分析多孔材料的样品，获取被检测多孔材料外表面即剖面一 1的 孔隙结构；再平行切除剖面一1上的一薄层材料，获取新形成的剖面二2 的孔隙结构图像；最后，再切除一薄层材料，获取新形成的剖面三3的孔 隙结构图像。
步骤二、将所述多个孔隙结构图像分别传输至计算机进行预处理，本 实施例中，即将所获取的三个相互平行剖面的孔隙结构图像分别传输至计算机，通过计算机依次对所述三个孔隙结构图像进行预处理，其预处理过
程包括以下步骤
U)将所述孔隙结构图像转化为多灰度级灰度图像即将RGB、 CMYK 等色彩模式转换为灰度的色彩模式。
(b )对所述多灰度级灰度图像进行去噪和滤波处理所述滤波处理为 最小值滤波和/或中值滤波和/或最大值滤波和/或均值滤波和/或低通滤 波和/或高通滤波处理。实践中，可以根据具体实际情况，相应选择最小 值滤波、中值滤波、最大值滤波、均值滤波、低通滤波和高通滤波中的一 种或几种滤波处理方式的组合方式进行滤波处理。
(c)将经去噪和滤波处理的所述多灰度级灰度图像转换为二灰度级 图像即黑白图像；
所述黑白图像由像素值为1的白色图像区域和像素值为0的黑色图像 区域组成，其中，所述实体骨架材料对应的图像区域为白色图像区域或黑 色图像区域，所述孔隙空间对应的图像区域相应为黑色图像区域或白色图 像区域。
也就是说，本步骤中，将所述孔隙结构图像中的所有像素值分别置为 0 (黑色图像区域)或1 (白色图像区域)，本实施例中，所述实体骨架材 料对应的图像区域为白色图像区域，所述孔隙空间对应的图像区域相应为 黑色图像区域。实际操作过程中，根据不同的多孔材料以及不同的孔隙结 构图像获取方式，黑色图像区域和白黑色图像区域分别与实体骨架材料和 孔隙空间的对应关系不同。
步骤三、用计算机提取所述黑白图像中黑色图像区域和白色图像区域 的交界线即孔隙空间边界并用彩色线条进行标示，所标示出来的孔隙空间 边界为由多个像素点为顶点组成的多边形，确定所述平面多边形所有顶点 的二维坐标值并对其各顶点进行连续编号，即获得被检测多孔材料的多个 相互平行断面的二维孔隙边界图像。另外，还需对所述平面多边形进行简 化处理，所述简化处理的方法为减少所述平面多边形的边数并相应增加各边边长。在确定所述平面多边形所有顶点的二维坐标值时，选取所述黑白 图像左下方的顶点为坐标原点，也就是说，以所述黑白图像左下方的顶点 为坐标原点，确定所述平面多边形各顶点的二维坐标值。对所述平面多边 形各顶点进行连续编号时，以二维坐标绝对值最小的顶点为1#顶点，然后
结合图1,本步骤中，从转换后的黑白图像中提取孔隙空间边界即孔 隙边界，为便于识别孔隙边界，利用计算机将黑色图像区域和白色图像区 域的交界线识别出来后，并用红色或蓝色等彩色线条在所述黑白图像中标 示出来，所标示出来的孔隙空间边界即孔隙边界是由多个像素点为顶点组 成的多边形。为降低后续计算处理的工作量，可对所标示出来的多边形进 行简化处理，其简化处理的方法是减少多边形的边数并增加多边形的边 长，优选的减少多边形的边数并增加多边形边长的实施方法是每隔若干个 像素点将首尾两个像素点连接起来构成简化的多边形，本实施例中，每隔
5个像素点将首尾两个像素点连接起来构成简化多边形4。
为方便对孔隙边界的多边形进行数学分析处理，则先对孔隙边界多边 形进行数值化处理，即利用计算机确定所述平面多边形所有顶点的二维坐 标值，并对所述平面多边形的所有顶点进行编号。本实施例中，确定所述 平面多边形所有顶点的二维坐标值时，选取所述黑白图像左下方的顶点为 坐标原点，而对所述平面多边形各顶点进行连续编号时，以二维坐标绝对 值最小的顶点为1#顶点，然后按照逆时针顺序依次对剩余的12个顶点进 行连续编号。
分别按照步骤三提取所述三个孔隙结构图像的孔隙边界并分别进行 数值化，便获得被检测多孔材料的三个相互平行剖面的二维孔隙边界图
步骤四、用计算机将相邻两个所述二维孔隙边界图像中所有平面多边 形的各顶点按编号顺序对应依次连接后，所组成的空间多边形即为被检测 多孔材料的三维孔隙边界图像。本实施例中，通过三维重建形成被检测多孔材料孔结构的三维孔隙边
界。如图3所示，按照步骤三获取剖面一l、剖面二2和剖面三3三个相 互平行剖面的二维孔隙边界图像后，将相邻两个二维孔隙边界图像中所述 平面多边形的各顶点按照编号顺序依次进行连接，即将对应的坐标值最接 近的顶点连接起来，例如，将剖面一1中的1#、 2#……13#顶点分别对应 与剖面二2中的1，#、 2，#……13，#顶点连接起来构成一个空间多边形，再 将剖面二 2中的顶点1，#、 2，#……13，#顶点分别对应与剖面三3中的1"#、 2"#...…13"#顶点连接起来构成另外一个空间多边形。将所有相邻剖面的 二维孔隙边界图像中的各顶点连接后所组成的空间多边形即三维面，便是 被检测多孔材料的三维孔隙边界。
以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制，凡是 根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构 变化，均仍属于本发明技术方案的保护范围内。
权利要求
1. 一种多孔材料孔隙边界提取和数值化方法，其特征在于该方法包括以下步骤步骤一、采用图像采集设备摄取被检测多孔材料的多个相互平行的断面的微观孔隙结构图像，其中，所述多孔材料的实体骨架材料和孔隙空间分别对应的图像区域之间具有明显的色彩对比；步骤二、将所述多个孔隙结构图像分别传输至计算机进行预处理，其预处理过程包括以下步骤(a)将所述孔隙结构图像转化为多灰度级灰度图像；(b)对所述多灰度级灰度图像进行去噪和滤波处理；(c)将经去噪和滤波处理的所述多灰度级灰度图像转换为二灰度级图像即黑白图像；所述黑白图像由像素值为1的白色图像区域和像素值为0的黑色图像区域组成，其中，所述实体骨架材料对应的图像区域为白色图像区域或黑色图像区域，所述孔隙空间对应的图像区域相应为黑色图像区域或白色图像区域；步骤三、用计算机提取所述黑白图像中黑色图像区域和白色图像区域的交界线即孔隙空间边界并用彩色线条进行标示，所标示出来的孔隙空间边界为由多个像素点为顶点组成的平面多边形，确定所述平面多边形所有顶点的二维坐标值并对其各顶点进行连续编号，即获得被检测多孔材料的多个相互平行断面的二维孔隙边界图像；步骤四、用计算机将相邻两个所述二维孔隙边界图像中所有平面多边形的各顶点按编号顺序对应依次连接后，所组成的空间多边形即为被检测多孔材料的三维孔隙边界。
2. 按照权利要求1所述的 一 种多孔材料孔隙边界提取和数值化方法， 其特征在于步骤一中所述的图像采集设备为金相显微镜或扫描电子显微 镜。
3. 按照权利要求1或2所述的 一 种多孔材料孔隙边界提取和数值化方 法，其特征在于步骤一中所述的孔隙结构图像为分辨率不低于500 x 500的清晰彩色或多灰度级灰度图像。
4. 按照权利要求l或2所述的一种多孔材料孔隙边界提取和数值化方 法，其特征在于步骤(c)中所述的滤波处理为最小值滤波和/或中值滤波和/或最大值滤波和/或均值滤波和/或低通滤波和/或高通滤波处理。
5. 按照权利要求l或2所述的一种多孔材料孔隙边界提取和数值化方 法，其特征在于步骤(c)中所述的实体骨架材料对应的图像区域为白色图像区域，所述孔隙空间对应的图像区域相应为黑色图像区域。
6. 按照权利要求1或2所述的 一 种多孔材料孔隙边界提取和数值化方法，其特征在于步骤三中还需对所述平面多边形进行简化处理，所述简 化处理的方法为减少所述平面多边形的边数并相应增加各边边长。
7. 按照权利要求1或2所述的 一 种多孔材料孔隙边界提取和数值化方 法，其特征在于步骤三中确定所述平面多边形所有顶点的二维坐标值时， 选取所述黑白图像左下方的顶点为坐标原点。
8. 按照权利要求1或2所述的 一种多孔材料孔隙边界提取和数值化方 法，其特征在于步骤三中对所述平面多边形各顶点进行连续编号时，以 二维坐标绝对值最小的顶点为1#顶点，然后按照逆时针顺序依次对剩余顶 点进行连续编号。
9. 按照权利要求1或2所述的一种多孔材料孔隙边界提取和数值化方 法，其特征在于步骤一中所述相互平行断面的数量不少于三个。
10. 按照权利要求1或2所述的一种多孔材料孔隙边界提取和数值化 方法，其特征在于步骤一中所述相互平行断面的数量为三个。
全文摘要
本发明公开了一种多孔材料孔隙边界提取和数值化方法，该方法包括以下步骤摄取被检测多孔材料的多个相互平行断面的微观孔隙结构图像；将多个孔隙结构图像分别传输至计算机进行预处理；提取黑白孔隙结构图像中黑色图像区域和白色图像区域的交界线即孔隙空间边界并进行标示，标示出来的孔隙空间边界为由多个像素点为顶点组成的多边形，确定多边形所有顶点的二维坐标值并对其各顶点进行连续编号，即获得二维孔隙边界；将相邻两个二维孔隙边界图像中所有平面多边形的各顶点按编号顺序对应依次连接后，所组成的空间多边形即为三维孔隙边界。本发明操作简单且成本低，能对多孔材料的孔隙形貌和孔通道的弯曲程度等性能参数进行数值化处理，分析准确。
文档编号G01B15/00GK101413792SQ200810232608
公开日2009年4月22日 申请日期2008年12月5日 优先权日2008年12月5日
发明者奚正平, 敖庆波, 朱纪磊, 汤慧萍, 王建永 申请人:西北有色金属研究院