Page 41 - 理化检验-物理分册2021年第十二期
P. 41

刘吉华, 等: MATLAB软件在铁路车轮钢组织定量分析中的应用



                                                                                       与波形的交点所对应的
                                                                   ( 3 )取d 的最大值d m
                                                               阈值h 就是所选的阈值。

                                                                   ( 4 )在最大波峰对应位置不在直方图亮侧时,
                                                               需要进行直方图翻转, 翻转后求出阈值 T , 再用255
                                                               减去阈值T , 即为所求的阈值h 。
                                                                   结合上文直方图分析可得, 含量较高的珠光体
                                                               相为 SEM 图中较亮的区域, 而含量较少的先共析
                                                               铁素体较暗, 因此直方图波峰均出现在图像亮侧, 所
                                                               以在计算中不用考虑原三角形算法中直方图翻转的
                                                               情况, 可以将算法进行简化, 简化后的算法为

                                                                                               , )( 波峰) 和
                                                                   ( 1 )连接单峰直方图最大点( x m ym
                                                                          , ), 形成线段L 。
                                                               最小值点( x sy s

                                                                   ( 2 )作L 的垂线与波形相交, 形成垂线段d 。
                                                                                       与波形的交点所对应的

                                                                   ( 3 )取d 的最大值d m
                                                               阈值h 就是所选的阈值。
                 图7 对图1c ) 进行直方图规定化前后的图像分割对比                       简化后的算法更加简洁, 同时可避免因为材料

                Fi g  7 Com p aretheima g ese g mentationbeforeandafter
                                                               中杂质造成计算过程中直方图翻转计算的情况, 更

                       histo g rams p ecificationforFi g  1c
                                                               适合珠光体 - 铁素体材料 SEM 图的图像分割。图9
             a  beforehisto g rams p ecification b   afterhisto g rams p ecification
                                                               为图1中各图二值化分割后的图像。
            时, 小的波峰干扰导致图像分割失败; 使用直方图规
                                                              3.4 特征提取
            定化对图像增强后, 图像直方图符合三角形算法要
                                                                   使用三角形算法计算出合适的阈值对图像进行
            求, 使用三角形算法能成功对图像进行分割。对图
                                                               分割后, 图像不再涉及其他灰度值像素, 只有 0 或
            像进行预处理后, 能有效地减少金相图像中存在的
            噪声, 同时还能对图像进行有目的地增强, 为后续图                         255 。图像中的目标和背景即被分离, 小于阈值的黑
                                                               色部分即为所需计算的铁素体相。
            像分割奠定了基础。
            3.3 图像分割                                               计算出图像中黑色像素点数量后, 与所计算图
                 使用三角形算法对图像进行分割, 三角形算法                         像总像素点之比即为该图像中先共析铁素体相所占
            是 PAUL于2001年提出的, 基于纯几何方法寻找                         比例。基于三角形算法对图像二值化分割, 进行特
            最佳阈值对图像进行分割的方法                 [ 12 ] 。图 8 为三角     征提取后, 计算结果见表2 。
            形算法说明示意图, 三角形算法的原理为
                                                              4 结果对比

                                                                   传统金相分析与基于数字图像处理的定量金相
                                                               分析结果对比如图 10 所示。可以发现, 1 号材料
                                                              3000×图像计算结果相差较大, 分析认为, 数字图

                                                               像处理过程中, 材料中的类珠光体碳化物因为颜色
                                                               较深, 在二值化时被识别为先共析铁素体, 导致计算

                                                               结果偏大。 1号材料除3000×图像外, 考虑计算中
                                                               的误差, 本算法计算的先共析铁素体结果与传统定
                           图8 三角形算法示意图                         量金相分析结果接近。 2号材料传统金相方法测算

                  Fi g  8 Schematicdia g ramofthetrian g leal g orithm  结果均大于数字图像处理计算得到的结果, 造成该

                                                  , )( 波       结果的原因是相分布的不均匀性。金相分析测定相
                 ( 1 )连接单峰直方图的最大值点( x m ym
                             , ), 波谷在图像的左侧还是右                  含量会因为放大倍数和测量视场造成误差                    [ 14 ] , 在传
            峰) 和最小值点( x s y s
            侧, 视实际情况而定, 形成线段L 。                                统金相方法测量相含量过程中, 会选择视场内有先

                 ( 2 )作L 的垂线与波形相交, 形成垂线段d 。                    共析铁素体的区域进行测算, 当先共析铁素体在视
                                                                                                         2 5
   36   37   38   39   40   41   42   43   44   45   46