(57)【要約】【課題】より正確に火花の破裂を認識することが可能な破裂火花認識装置、鋼材識別装置、破裂火花認識方法及び破裂火花認識プログラムを提供する。【解決手段】鋼材識別装置１は、鋼材Ｂを研削して発生する火花の画像である火花画像に含まれる複数の流線の各々の長さに基づいて、主幹流線Ｒｍと、枝流線Ｒｂとを分類する流線分類部１５と、流線分類部１５により主幹流線Ｒｍに分類された線上の任意の点を中心にした探索領域Ｅに少なくとも一部を含む枝流線Ｒｂの数が所定数以上の場合に、破裂として認識する破裂認識部１６と、を備える。【選択図】図１
2020134204.tif 000002
【特許請求の範囲】
【請求項１】
鋼材を研削して発生する火花の画像である火花画像に含まれる複数の流線の各々の長さに基づいて、主幹流線と、枝流線とを分類する流線分類部と、
前記流線分類部により前記主幹流線に分類された線上の任意の点を中心にした規定領域に少なくとも一部を含む前記枝流線の数が所定数以上の場合に、破裂として認識する破裂認識部と、
を備える破裂火花認識装置。
【請求項２】
請求項１に記載の破裂火花認識装置において、
前記破裂認識部は、既に破裂として認識された前記枝流線を除き、前記流線分類部により前記主幹流線に分類された線上の任意の点を中心にした規定領域に少なくとも一部を含む前記枝流線の数が所定数以上の場合に、破裂として認識する破裂火花認識装置。
【請求項３】
請求項１又は請求項２に記載の破裂火花認識装置において、
前記破裂認識部は、前記規定領域を、前記火花の根本から前記火花の先端に向かうほど大きな領域にする破裂火花認識装置。
【請求項４】
請求項１から請求項３までのいずれか一項に記載の破裂火花認識装置において、
前記破裂認識部は、前記主幹流線に平行する前記枝流線を、前記枝流線の数から除外する破裂火花認識装置。
【請求項５】
請求項１から請求項４までのいずれか一項に記載の破裂火花認識装置において、
前記流線について、複数の角度に相当する短直線を表す複数のテンプレートを用いて前記火花の先端から前記火花の根本に向かって順番にマッチングしながら短直線を連結し、各マッチングした前記短直線が前記流線に収まるテンプレートの種類、位置及び連結関係を記憶するマッチング処理部を備え、
前記流線分類部は、前記マッチング処理部によりマッチングした前記短直線の数に基づいて前記分類を行う破裂火花認識装置。
【請求項６】
請求項５に記載の破裂火花認識装置において、
前記流線分類部は、前記マッチング処理部によりマッチングした前記短直線の数が第１の値以上の場合に、前記主幹流線に分類し、前記マッチング処理部によりマッチングした前記短直線の数が前記第１の値未満であって、前記第１の値より小さい第２の値以上の場合に、前記枝流線に分類する破裂火花認識装置。
【請求項７】
請求項１から請求項６までのいずれか一項に記載の破裂火花認識装置を含む鋼材識別装置であって、
前記破裂認識部により認識された前記破裂の数と、前記流線分類部により分類された前記主幹流線に係る数とに基づいて前記火花画像全体の破裂密度を算出する破裂密度算出部と、
前記破裂密度算出部により算出した前記破裂密度に基づいて前記鋼材内の成分含有量を推定し、前記鋼材を識別する鋼材識別処理部と、
を備える鋼材識別装置。
【請求項８】
鋼材を研削して発生する火花の画像である火花画像から破裂を認識する破裂火花認識方法であって、
コンピュータが、
前記火花画像に含まれる複数の流線の各々の長さに基づいて、主幹流線と、枝流線とを分類する流線分類ステップと、
前記流線分類ステップにより前記主幹流線に分類された線上の任意の点を中心にした規定領域に少なくとも一部を含む前記枝流線の数が所定数以上の場合に、破裂として認識する破裂認識ステップと、
を含む破裂火花認識方法。
【請求項９】
コンピュータである破裂火花認識装置で実行する破裂火花認識プログラムであって、
前記破裂火花認識装置を、
鋼材を研削して発生する火花の画像である火花画像に含まれる複数の流線の各々の長さに基づいて、主幹流線と、枝流線とを分類する流線分類手段、及び
前記流線分類手段により前記主幹流線に分類された線上の任意の点を中心にした規定領域に少なくとも一部を含む前記枝流線の数が所定数以上の場合に、破裂として認識する破裂認識手段
として機能させるための破裂火花認識プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、破裂火花認識装置、鋼材識別装置、破裂火花認識方法及び破裂火花認識プログラムに関する。
【背景技術】
【０００２】
鋼の製造工程では、異材を鑑別して排除するために火花試験が広く用いられている。火花試験とは、グラインダを使用して鋼材を研削し、発生する火花の特徴を観察することによって、鋼種の推定又は異材の鑑別を行う試験のことをいい、ＪＩＳ  Ｇ  ０５６６に規定されている（例えば、非特許文献１参照）。
【０００３】
従来、火花試験は、鉄鋼材の検査工程等で熟練した経験をもった検査員が目視観察にて官能検査として行われてきた。しかし、個人差や環境の変動によって判定結果がばらついて適正検査結果を得ることが困難であった。また、官能検査の必然性として検査結果が記録できないため、検査技術はもっぱら経験もしくは伝承によるところが大きく、技術改善を評価することが困難であった。
【０００４】
この点、判定結果が個人差や環境の変動に影響される度合いを抑制する技術として、火花試験を装置により自動的に行う技術が開示されている（例えば、特許文献１参照）。
具体的には、特許文献１は、鋼材を研削して発生する火花の画像をグレースケールの火花画像に変換する輝度変換処理部と、グレースケールの火花画像の画素毎に所定の閾値で二値化を行う二値化処理部と、二値化された火花画像の所定幅の芯線に対し、複数の角度に相当する複数のテンプレートを、複数のテンプレートの短直線を延長して各々マッチングし、延長した短直線が芯線に収まるテンプレートの種類と位置を記憶する短直線マッチング処理部と、マッチングされたテンプレートが任意の火花画像の範囲において所定の数以上である場合に火花の破裂部とみなして破裂部を抽出し、火花の破裂数及び火花画像内のマッチングされたテンプレートの総数をカウントする破裂部抽出処理部と、火花の破裂数及びマッチングされたテンプレートの総数に基づいて鋼材を識別する鋼材識別処理部とを有し、とりわけ鋼材に含まれる炭素を識別することが可能な鋼材成分識別装置を開示している。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１６−０９９１２６号公報
【非特許文献】
【０００６】
【非特許文献１】日本工業規格、ＪＩＳ  Ｇ  ０５６６  鋼の火花試験方法、日本規格協会
【発明の概要】
【発明が解決しようとする課題】
【０００７】
上述した特許文献１で開示される技術は、テンプレートが任意の火花画像の範囲において所定の数以上である場合に、火花の破裂部とみなすものである。この技術は、鋼材中の炭素量の確認に有効なものではあるものの、実際の破裂と比較すると、若干の誤差が生じるものであった。
【０００８】
そこで、本発明は、より正確に火花の破裂を認識することが可能な破裂火花認識装置、鋼材識別装置、破裂火花認識方法及び破裂火花認識プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００９】
（１）本開示の一態様は、破裂火花認識装置が、鋼材を研削して発生する火花の画像である火花画像に含まれる複数の流線の各々の長さに基づいて、主幹流線と、枝流線とを分類する流線分類部と、前記流線分類部により前記主幹流線に分類された線上の任意の点を中心にした規定領域に少なくとも一部を含む前記枝流線の数が所定数以上の場合に、破裂として認識する破裂認識部と、を備える。
【００１０】
（２）  （１）に記載の破裂火花認識装置において、前記破裂認識部は、既に破裂として認識された前記枝流線を除き、前記流線分類部により前記主幹流線に分類された線上の任意の点を中心にした規定領域に少なくとも一部を含む前記枝流線の数が所定数以上の場合に、破裂として認識するものであってもよい。
【００１１】
（３）  （１）又は（２）に記載の破裂火花認識装置において、前記破裂認識部は、前記規定領域を、前記火花の根本から前記火花の先端に向かうほど大きな領域にするものであってもよい。
【００１２】
（４）  （１）から（３）までのいずれかに記載の破裂火花認識装置において、前記破裂認識部は、前記主幹流線に平行する前記枝流線を、前記枝流線の数から除外するものであってもよい。
【００１３】
（５）  （１）から（４）までのいずれかに記載の破裂火花認識装置において、前記流線について、複数の角度に相当する短直線を表す複数のテンプレートを用いて前記火花の先端から前記火花の根本に向かって順番にマッチングしながら短直線を連結し、各マッチングした前記短直線が前記流線に収まるテンプレートの種類、位置及び連結関係を記憶するマッチング処理部を備え、前記流線分類部は、前記マッチング処理部によりマッチングした前記短直線の数に基づいて前記分類を行うものであってもよい。
【００１４】
（６）  （５）に記載の破裂火花認識装置において、前記流線分類部は、前記マッチング処理部によりマッチングした前記短直線の数が第１の値以上の場合に、前記主幹流線に分類し、前記マッチング処理部によりマッチングした前記短直線の数が前記第１の値未満であって、前記第１の値より小さい第２の値以上の場合に、前記枝流線に分類するものであってもよい。
【００１５】
（７）  （１）から（６）までのいずれかに記載の破裂火花認識装置を含む鋼材識別装置が、前記破裂認識部により認識された前記破裂の数と、前記流線分類部により分類された前記主幹流線に係る数とに基づいて前記火花画像全体の破裂密度を算出する破裂密度算出部と、前記破裂密度算出部により算出した前記破裂密度に基づいて前記鋼材内の成分含有量を推定し、前記鋼材を識別する鋼材識別処理部と、を備える。
【００１６】
（８）  本開示の一態様は、鋼材を研削して発生する火花の画像である火花画像から破裂を認識する破裂火花認識方法であって、コンピュータが、前記火花画像に含まれる複数の流線の各々の長さに基づいて、主幹流線と、枝流線とを分類する流線分類ステップと、前記流線分類ステップにより前記主幹流線に分類された線上の任意の点を中心にした規定領域に少なくとも一部を含む前記枝流線の数が所定数以上の場合に、破裂として認識する破裂認識ステップと、を含む。
【００１７】
（９）  本開示の一態様は、コンピュータである破裂火花認識装置で実行する破裂火花認識プログラムであって、前記破裂火花認識装置を、鋼材を研削して発生する火花の画像である火花画像に含まれる複数の流線の各々の長さに基づいて、主幹流線と、枝流線とを分類する流線分類手段、及び前記流線分類手段により前記主幹流線に分類された線上の任意の点を中心にした規定領域に少なくとも一部を含む前記枝流線の数が所定数以上の場合に、破裂として認識する破裂認識手段として機能させる。
【発明の効果】
【００１８】
本発明によれば、より正確に火花の破裂を認識することが可能な破裂火花認識装置、鋼材識別装置、破裂火花認識方法及び破裂火花認識プログラムを提供することができる。
【図面の簡単な説明】
【００１９】
【図１】本実施形態に係る鋼材識別システムの全体構成図である。
【図２】本実施形態に係る鋼材識別装置の機能ブロックを示す図である。
【図３】本実施形態に係る鋼材識別装置の画像分析処理を示すフローチャートである。
【図４】本実施形態に係る鋼材識別装置で生成される火花画像をグレースケール化した画像の例を示す図である。
【図５Ａ】本実施形態に係る鋼材識別装置で用いる十字二値化法を説明するための図である。
【図５Ｂ】本実施形態に係る鋼材識別装置で用いる十字二値化法を説明するための図である。
【図６】本実施形態に係る鋼材識別装置で生成される火花画像を二値化処理した画像の例を示す図である。
【図７】本実施形態に係る鋼材識別装置の短直線マッチング処理を示すフローチャートである。
【図８】本実施形態に係る鋼材識別装置での基準短直線探索処理を説明するための図である。
【図９】本実施形態に係る鋼材識別装置での基準短直線探索処理を説明するための図である。
【図１０】本実施形態に係る鋼材識別装置での基準短直線探索処理の例を示す図である。
【図１１Ａ】本実施形態に係る鋼材識別装置での流線分類処理を説明するための図である。
【図１１Ｂ】本実施形態に係る鋼材識別装置での流線分類処理を説明するための図である。
【図１１Ｃ】本実施形態に係る鋼材識別装置での流線分類処理を説明するための図である。
【図１２】本実施形態に係る鋼材識別装置の破裂認識処理を示すフローチャートである。
【図１３】本実施形態に係る鋼材識別装置での破裂認識処理を説明するための図である。
【図１４】本実施形態に係る鋼材識別装置の探索領域を説明するための図である。
【図１５】本実施形態に係る鋼材識別装置の実施例における撮影環境を示す図である。
【図１６】本実施形態に係る鋼材識別装置を用いて各方法にて算出した破裂密度の相関係数を比較した図である。
【図１７】鋼材の炭素含有量と新旧の手法によって実測した破裂密度とを比較したグラフである。
【発明を実施するための形態】
【００２０】
以下、本発明を実施するための形態について、図を参照しながら説明する。なお、これは、あくまでも一例であって、本発明の技術的範囲はこれに限られるものではない。
（実施形態）
〔鋼材識別システム１００〕
図１は、本実施形態に係る鋼材識別システム１００の全体構成図である。
図２は、本実施形態に係る鋼材識別装置１の機能ブロックを示す図である。
【００２１】
図１に示す鋼材識別システム１００は、検査対象の鋼材ＢとグラインダＣとが接触して発生した火花を撮像するカメラ５と、カメラ５で撮像された火花画像に対する画像処理を行う鋼材識別装置１とを備える。
【００２２】
まず、鋼材識別システム１００の前提事項について説明する。
鋼材識別システム１００において、定量的な鋼材分類を行うためには、同種の鋼材の火花に差異が発生しないよう、常に同条件で火花を発生させることが好ましい。そこで、鋼材識別システム１００では、モータ（不図示）によって発生したトルク回転力を直線運動に変換し、グラインダＣに鋼材Ｂを押し当てる機構を採用している。モータは、出力トルクを制御できるトルクモータを採用し、一定の力で、一定の位置で、かつ、一定の角度での鋼材Ｂの押し当てを可能とした。
【００２３】
鋼材識別システム１００での撮影環境は、周囲が薄暗い環境の中、明るい火花を撮影するという特殊な環境である。そのため、カメラ５は、シャッタースピードや絞り等をマニュアルで設定できるものが好ましい。また、鋼材ＢがグラインダＣに削られることにより、鋼材ＢとグラインダＣの接触面積が増加して、押し当て圧が下がり火花の発生が時間変化するという問題がある。そこで、鋼材識別システム１００では、火花の発生の時間変化の影響を低減するために、短時間で高速の連続撮影が可能なカメラ５を用いる。
【００２４】
同様に、時間経過による火花の変化を避けるため、カメラ５は、撮影時間を短くすることが好ましい（例えば１．０秒以下）。また、画像で確認できる火花の線（以降の説明において、「流線」ともいう。）の長さや破裂の数は、カメラ５のシャッタースピードに依存し変化する。このため、鋼材識別システム１００のカメラ５では、解析に適するシャッタースピードの選定（例えば、１／１６０ｓ等）をした上で、撮影を行うことが好ましい。
【００２５】
なお、鋼材識別システム１００の詳細については、後述の実施例で記載する。
以下の説明においては、例として、ＤＩＴＥＣＴ社製のカメラＨＡＳ−Ｌ２を用いて３０ｆｐｓで撮像し、解像度が２５９２×２０４８ｐｉｘｅｌの画像を５０枚用いて解析したケースにおける数値を用いることがある。しかし、本発明は、これに限定されるものではない。
【００２６】
〔鋼材識別装置１〕
鋼材識別装置１は、カメラ５が撮影した火花画像から破裂を認識し、認識した破裂の数を用いて鋼材Ｂの成分含有量を推定する。
鋼材識別装置１は、例えば、サーバである。鋼材識別装置１は、１台のサーバに限定されるものではなく、複数台のサーバによって構成されるものであってもよい。また、鋼材識別装置１は、その他、パーソナルコンピュータ（ＰＣ）等であってもよい。
鋼材識別装置１は、以下に説明する処理を行うことにより、火花画像に対する画像処理を行って、鋼材の火花の定量化を行う。
【００２７】
図２に示すように、鋼材識別装置１は、演算処理を行うＣＰＵ１０と、データのワークエリアであるＲＡＭ２１と、ＣＰＵ１０の制御プログラム等を記憶するＲＯＭ２２とを備える。
ＣＰＵ１０は、ＲＯＭ２２に格納された制御プログラムを読み出して実行することにより、上述したハードウェアと協働し、制御プログラムにしたがって各種機能を実行する。
【００２８】
ＣＰＵ１０は、画像取込処理部１１と、輝度変換処理部１２と、二値化処理部１３と、短直線マッチング処理部１４と、流線分類部１５と、破裂認識部１６と、破裂密度算出部１７と、鋼材識別処理部１８とを備える。
画像取込処理部１１は、カメラ５により撮像された火花画像を、カラー画像としてＲＡＭ２１に取り込む。
輝度変換処理部１２は、画像取込処理部１１が取り込んだ火花画像の各画素の輝度を変換処理して、画像を輝度値のみのグレースケール画像に変換し、変換後の火花画像をＲＡＭ２１に格納する。
二値化処理部１３は、輝度変換処理部１２による変換後の火花画像に対して二値化処理を実行し、二値化処理後の火花画像をＲＡＭ２１に格納する。
【００２９】
短直線マッチング処理部１４は、二値化処理部１３による二値化処理後の火花画像に対して、短直線マッチング処理を行う。
流線分類部１５は、短直線マッチング処理部１４による短直線マッチング処理を行った後の短直線の数により表される流線の長さに基づいて、主幹流線と、枝流線とに分類する。
破裂認識部１６は、流線分類部１５により主幹流線に分類された線上の任意の点を中心にした規定領域に少なくとも一部を含む枝流線の数が所定数以上の場合に、破裂として認識する。
【００３０】
なお、破裂火花認識装置として、短直線マッチング処理部１４と、流線分類部１５と、破裂認識部１６との機能を実現することも可能である。そして、破裂火花認識装置は、本実施形態のように、鋼材識別装置１に含まれてもよいし、鋼材識別装置１とは別のコンピュータによって実現されてもよい。
【００３１】
破裂密度算出部１７は、破裂認識部１６により認識された破裂の数と、流線分類部１５により分類された主幹流線に係る数とに基づいて、火花画像全体の破裂密度を算出する。
鋼材識別処理部１８は、破裂密度算出部１７により算出した破裂密度に基づいて鋼材Ｂ内の成分含有量を推定し、鋼材Ｂを識別する。
【００３２】
ＲＯＭ２２は、画像取込処理部１１と、輝度変換処理部１２と、二値化処理部１３と、短直線マッチング処理部１４と、流線分類部１５と、破裂認識部１６と、破裂密度算出部１７と、鋼材識別処理部１８とを実現する制御プログラムとして、鋼材識別プログラムを含む。なお、ＣＰＵ１０が実行する機能部ごとのプログラムを有してもよいし、例えば、短直線マッチング処理部１４と、流線分類部１５と、破裂認識部１６とを実現する制御プログラムとして、破裂火花認識プログラムを有してもよい。
【００３３】
〔処理の説明〕
次に、鋼材識別装置１における処理について、フローチャートに基づき説明する。
図３は、本実施形態に係る鋼材識別装置１の画像分析処理を示すフローチャートである。
図３のステップＳ（以下、「ステップＳ」を単に「Ｓ」という。）１において、鋼材識別装置１の画像取込処理部１１は、カメラ５が撮影した火花画像を、カラー画像としてＲＡＭ２１に取り込む。ここで、鋼材識別装置１は、無線ＬＡＮ等の通信ネットワークを介してカメラ５から火花画像を受信してもよい。また、鋼材識別装置１は、記憶媒体を介してカメラ５が撮影した火花画像を受け付けてもよい。
【００３４】
Ｓ２において、輝度変換処理部１２は、ＲＡＭ２１に取り込んだ火花画像の各画素の輝度を変換処理して、画像を輝度値のみのグレースケール画像に変換し、変換後の火花画像をＲＡＭ２１に格納する。
【００３５】
図４は、火花画像をグレースケール化した画像３１の例である。画像３１により、火花画像は、複数の火花の線（流線）を含むことがわかる。本実施形態では、火花画像中の各画素のＲＧＢ値を、所定の変換式によってＹＩＱ表示系の輝度値Ｙへ変換している。所定の変換式として、例えば、以下の式を用いることが可能である。
【数１】
2020134204.tif 000003
なお、後述するが、カラー画像からグレースケール画像に変換するための方法は、これに限定されるものではない。
【００３６】
図３のＳ３において、二値化処理部１３は、ＲＡＭ２１に格納したグレースケール化した火花画像に対して二値化処理を実行し、変換後の火花画像をＲＡＭ２１に格納する。
本実施形態では、二値化処理の対象となる画素毎に、所定の閾値で二値化処理を行う。二値化を行うことにより流線がより明瞭になる。本実施形態では、抽出対象が火花のような線状の場合に適している十字二値化を用いる。
【００３７】
図５Ａ及び図５Ｂは、二値化処理部１３が実行する十字二値化法を説明する図である。十字二値化法は、十字内の中心画素以外の画素の平均輝度値Ｘから十字の中心画素の輝度値Ｙを引き、その値が設定した閾値よりも大きい場合に、十字の中心画素を黒とし（図５Ａ）、そうでなければ白とする（図５Ｂ）手法である。さらに、今回は火花を黒画素とするため、十字二値化画像の白黒の出力を反転させている。図６に、画像３１に対して十字二値化して白黒反転させた画像３２の例を示す。
【００３８】
図３のＳ４において、短直線マッチング処理部１４は、火花画像の短直線マッチング処理を行う。
ここで、短直線マッチング処理部１４による処理について、図７に示すフローチャートに基づき説明する。
【００３９】
〔短直線マッチング処理〕
図７は、本実施形態に係る鋼材識別装置１の短直線マッチング処理を示すフローチャートである。
以下に説明する短直線マッチング処理は、火花画像に含まれる流線を認識して記録するための処理である。
図７のＳ１１において、短直線マッチング処理部１４は、基準短直線探索処理を行う。
図８に示すように、短直線マッチング処理部１４は、二値化処理された画像に含まれる１つの流線Ｒについて、流線Ｒ内で、座標を（ｘ，ｙ）とする任意の黒画素を、基準短直線Ｌｓ（短直線）の開始点（始点）とする。この座標（ｘ，ｙ）の黒画素は、既に求めた基準短直線Ｌｓから与えられた黒画素でもよく、流線Ｒ内の画素を順番に探すことにより検知された黒画素でもよい。図８に示す例では、火花の線の先端に座標（ｘ，ｙ）が設定されている。
【００４０】
ここで、短直線マッチング処理部１４は、座標を（ｘ’，ｙ’）＝（ｘ＋Ｌ×ｓｉｎθ，ｙ＋Ｌ×ｃｏｓθ）（ただし、１≦Ｌ≦１５、Ｌは整数）とする領域を、θを５度刻みとしながら探索する。そして、短直線マッチング処理部１４は、（ｘ，ｙ）から（ｘ’，ｙ’）までの間に、１０ｐｉｘ／１５ｐｉｘ（６６％）以上、黒画素が存在すれば、（ｘ，ｙ）から（ｘ’，ｙ’）までを基準短直線Ｌｓとして設定する。
そして、図９に示すように、短直線マッチング処理部１４は、上記の処理で求めた基準短直線Ｌｓの先端を新たな開始点として、流線Ｒの終点（火花の根本方向）に向かって基準短直線探索処理を繰り返す。
【００４１】
この短直線マッチング処理部１４による処理は、具体的には、予め角度の異なる短直線を表す複数のテンプレートを用意して行うことができる。
特許文献１に開示されているものは、２８方向の角度（６．５度毎）に相当し、長さが１０ｐｉｘｅｌの短直線を表すテンプレートを用いている。今回、本実施形態では、精度を向上させるべく、上述した条件に合致する、５度毎に相当し、長さを１５ｐｉｘｅｌとする３７方向のテンプレートを用いた。
【００４２】
なお、短直線マッチング処理部１４は、短直線が複数存在する場合には、最も多く黒画素がマッチングした短直線を基準短直線Ｌｓとする。
また、黒画素の数が同数の短直線が複数存在する場合には、短直線マッチング処理部１４は、各テンプレートの短直線を延長し、各テンプレートに含まれる画素のうち一致する流線Ｒの画素数が最多のテンプレートに対応する短直線を基準短直線Ｌｓとする。
さらに、各テンプレートの短直線を延長し、各テンプレートに含まれる画素のうち一致する流線Ｒの画素数が最多のテンプレートが複数ある場合には、短直線マッチング処理部１４は、角度θが最も大きな直線を、基準短直線Ｌｓとする。
【００４３】
そして、短直線マッチング処理部１４は、マッチングしたテンプレートの種類（パターン番号）と位置とを、ＲＡＭ２１に格納する。また、短直線マッチング処理部１４は、流線Ｒについての連結関係を示すものとして、１つの流線Ｒを構成する連結された１以上の基準短直線Ｌｓに一意となる同一のＩＤを付与して、ＲＡＭ２１に格納する。なお、テンプレートの種類には、角度情報が含まれる。
【００４４】
このように、テンプレートを用いた短直線マッチング処理では、角度の異なる複数のテンプレートが合致した場合、その短直線を仮想的に延長し、より合致するものを基準短直線Ｌｓとして選択することができる。こうすることで、ＣＰＵ１０は、より正確な流線Ｒの角度を得ることができる。
【００４５】
図１０は、流線Ｒについて、短直線マッチング処理部１４によって基準短直線Ｌｓをマッチングさせていく態様を示す。
図１０（Ａ）は、短直線マッチング処理部１４による処理を行う前の流線Ｒである。
図１０（Ｂ）は、短直線マッチング処理部１４による処理を繰り返し行っている途中の状態を示し、３つの基準短直線Ｌｓを設定した流線Ｒである。ここで、短直線マッチング処理部１４は、マッチングした基準短直線Ｌｓを順番に連結させることで、基準短直線Ｌｓをつなげていく。
図１０（Ｃ）は、短直線マッチング処理部１４による処理後の流線Ｒである。この図１０に示す流線Ｒは、５つの基準短直線Ｌｓによって構成されていることが分かる。
【００４６】
この図７のＳ１１の処理によって、各流線Ｒの長さを、基準短直線Ｌｓの数によって算出できる。
【００４７】
〔流線分類処理〕
次に、流線Ｒについて、短直線マッチング処理部１４によるマッチング後の処理について説明する。
図７のＳ１２からＳ１５までの処理は、流線Ｒを分類するための流線分類処理である。
流線分類部１５は、流線Ｒを、マッチングした基準短直線Ｌｓの数（以降、マッチングした基準短直線Ｌｓの数を、単に「マッチング数」ともいう。）に応じて、主幹流線と、枝流線と、それ以外との３つに分類する。本実施形態では、例えば、マッチング数が７以上の場合に、主幹流線とする。図１１Ａに示す流線Ｒ１は、マッチング数の数が８であるため、主幹流線に分類される。また、マッチング数が２以上７未満の場合に、枝流線とする。図１１Ｂに示す流線Ｒ２は、マッチング数が５であるため、枝流線に分類される。さらに、マッチング数が２未満の場合には、主幹流線でも枝流線でもないとする。図１１Ｃに示す流線Ｒ３は、マッチング数が１であるため、主幹流線にも枝流線にも分類されない。
【００４８】
図７のＳ１２において、流線分類部１５は、基準短直線探索処理によって探索をした流線Ｒに含まれるマッチング数がＮ（第１の値；例えば、７）以上であるか否かを判断する。マッチング数がＮ以上である場合（Ｓ１２：ＹＥＳ）には、流線分類部１５は、処理をＳ１３に移す。他方、マッチング数がＮ以上ではない場合（Ｓ１２：ＮＯ）には、流線分類部１５は、処理をＳ１４に移す。
Ｓ１３において、流線分類部１５は、当該流線Ｒを主幹流線に分類して、ＲＡＭ２１に格納する。その後、流線分類部１５は、処理をＳ１６に移す。
【００４９】
Ｓ１４において、流線分類部１５は、マッチング数がＭ（第２の値；例えば、２）以上であるか否かを判断する。マッチング数がＭ以上である場合（Ｓ１３：ＹＥＳ）には、流線分類部１５は、処理をＳ１５に移す。他方、マッチング数がＭ以上ではない場合（Ｓ１３：ＮＯ）には、流線分類部１５は、処理をＳ１６に移す。
Ｓ１５において、流線分類部１５は、当該流線Ｒを枝流線に分類して、ＲＡＭ２１に格納する。
なお、上述した例では、Ｎを７とし、Ｍを２としているが、ＮとＭとには、ＮがＭよりも大きい適宜の整数を設定してよい。
【００５０】
Ｓ１６において、短直線マッチング処理部１４は、画像中の全ての黒画素に対して処理をしたか否かを判断する。全ての黒画素に対して処理をした場合（Ｓ１６：ＹＥＳ）には、短直線マッチング処理部１４は、処理を図３のＳ５に移す。他方、全ての黒画素に対して処理をしていない場合（Ｓ１６：ＮＯ）には、短直線マッチング処理部１４は、処理をＳ１１に移す。
【００５１】
図３のＳ５において、破裂認識部１６は、破裂認識処理を行う。
ここで、破裂認識部１６による処理について、図１２に示すフローチャートに基づき説明する。
【００５２】
〔破裂認識処理〕
図１２は、本実施形態に係る鋼材識別装置１の破裂認識処理を示すフローチャートである。
以下に説明する破裂認識処理は、主幹流線に基づいて破裂を認識して記録するための破裂認識処理である。
図１２のＳ２１において、破裂認識部１６は、画像の左上の座標を、例えば、（０，０）として横方向Ｘに順番に画像を走査し、右側に到達したら高さ方向Ｙ（下方向）に走査して、主幹流線を探索する。
図１３（Ａ）は、画像を左上から順番に走査して、主幹流線Ｒｍを探索している態様を示す。
【００５３】
図１２のＳ２２において、破裂認識部１６は、主幹流線Ｒｍを見つけたか否かを判断する。主幹流線Ｒｍを見つけた場合（Ｓ２２：ＹＥＳ）には、破裂認識部１６は、処理をＳ２３に移す。主幹流線Ｒｍを見つけていない場合（Ｓ２２：ＮＯ）には、破裂認識部１６は、処理をＳ２６に移す。
【００５４】
Ｓ２３において、破裂認識部１６は、主幹流線Ｒｍを見つけた位置を基準点とし、基準点を中心とした探索領域Ｅ（規定領域）を探索範囲として、探索範囲に少なくとも部分的に含む枝流線ＲｂがＫ（Ｋは、所定数であり、任意の整数）本以上あるか否かを判断する。
図１３（Ｂ）は、主幹流線Ｒｍを見つけた際に、探索領域Ｅを設定した態様を示す。
ここで、探索領域Ｅは、図１４に示すように、基準点（ｘ，ｙ）を中心とした矩形であり、火花の先端部分である方が、火花の根本部分よりも大きさが大きいものとする。
探索領域Ｅの一辺の長さＬｅは、以下の式により設定する。
【００５５】
【数２】
2020134204.tif 000004
但し、ｗｉｄｔｈは、画像横幅（ｘ方向）であり、ｍａｘは、探索範囲最大値であり、ｍｉｎは、探索範囲最小値である。そして、例えば、ｍａｘ＝４０ｐｉｘ、ｍｉｎ＝２０ｐｉｘである。
上記の式に当てはめれば、ｘ＝０の場合には、長さＬｅがｍａｘになり、ｘ＝ｗｉｄｔｈの場合には、長さＬｅがｍｉｎになる。
このように、火花の先端部分と、火花の根本部分とで探索範囲を異なるものにするのは、火花の先端部ほど、火花の広がりが大きくなるため、探索範囲を広げる必要があるためである。
【００５６】
ここで、破裂認識部１６は、探索範囲に含まれる枝流線Ｒｂのうち、座標（ｘ，ｙ）における主幹流線Ｒｍと平行な枝流線Ｒｂは、カウントから除外する。これは、破裂として認識される枝流線Ｒｂは、主幹流線Ｒｍと平行にはならないという原則による。主幹流線Ｒｍと平行な枝流線Ｒｂがある場合には、これは、他の主幹流線Ｒｍの枝流線Ｒｂである可能性が高い。
主幹流線Ｒｍと平行な枝流線Ｒｂであるとは、例えば、枝流線Ｒｂを構成する短直線のテンプレートが、主幹流線Ｒｍを構成する短直線のテンプレートと同じテンプレートである場合をいう。また、主幹流線Ｒｍを構成する短直線のテンプレートに対して前後１つのテンプレートのいずれかを、枝流線Ｒｂを構成する短直線のテンプレートが有している場合にも、主幹流線Ｒｍと平行な枝流線Ｒｂであるとする。
【００５７】
また、破裂認識部１６は、既に破裂として認識された枝流線Ｒｂは、カウントから除外する。これは、既に破裂として認識された枝流線Ｒｂをカウントすると、１つの枝流線Ｒｂが複数の主幹流線Ｒｍでカウントされることになり、破裂を正確に認識できないためである。
【００５８】
探索範囲に枝流線ＲｂがＫ本以上ある場合（Ｓ２３：ＹＥＳ）には、破裂認識部１６は、処理をＳ２４に移す。他方、探索範囲に枝流線ＲｂがＫ本以上ない場合（Ｓ２３：ＮＯ）には、破裂認識部１６は、処理をＳ２６に移す。
例えば、Ｋを３とすると、図１３に示す例では、図１３（Ｂ）において、探索領域Ｅには、枝流線Ｒｂが１本もない。よって、破裂認識部１６は、探索範囲に枝流線Ｒｂが３本以上ないとして、処理をＳ２６に移す。他方、図１３（Ｃ）において、探索領域Ｅには、枝流線Ｒｂが３本ある。よって、破裂認識部１６は、探索範囲に枝流線Ｒｂが３本以上あるとして、処理をＳ２４に移し、破裂を認識することになる。
【００５９】
Ｓ２４において、破裂認識部１６は、破裂と認識し、ＲＡＭ２１に破裂番号を含む破裂情報を格納する。
Ｓ２５において、破裂認識部１６は、画像から探索範囲に少なくとも一部が含まれる枝流線Ｒｂを消去する。この処理は、上述したように、枝流線Ｒｂを１つの破裂の要素としてのみ用い、複数回カウントされないようにするためである。
Ｓ２６において、破裂認識部１６は、画像内を全て探索したか否かを判断する。画像内を全て探索した場合（Ｓ２６：ＹＥＳ）には、ＣＰＵ１０は、処理を図３のＳ６に移す。他方、画像内を全て探索していない場合（Ｓ２６：ＮＯ）には、ＣＰＵ１０は、処理をＳ２１に移す。
【００６０】
図３のＳ６において、破裂密度算出部１７は、火花の破裂密度の算出を行い、算出した破裂密度の数値により、鋼材中の炭素量を確認する。
ここで、破裂密度は、以下の式により算出される。
【数３】
2020134204.tif 000005
【００６１】
なお、本発明者らは、分母とした主幹流線Ｒｍを構成する基準短直線Ｌｓの総数（以下、主幹短直線数ともいう。）を、他の主幹流線Ｒｍに係る数に代えて破裂密度を算出した。しかし、上記の式で表される主幹短直線数を用いるのが、最も適切であり、つまり、破裂密度をより正確に表すものであった。
【００６２】
図３のＳ７において、鋼材識別処理部１８は、火花の破裂密度から鋼材中の炭素量を推定し、鋼材種類の認識を行う。
本実施形態では、テンプレートでマッチングされた主幹短直線数と、破裂の総数とに基づいて、鋼材の種類を定量的に分類することが可能となる。これが、破裂の総数／主幹短直線数で定義した破裂密度の技術的意味となる。
【００６３】
また、火花の形態は、立体的かつ「発生」、「成長」、「破裂」及び「消滅」の動的かつ不安定な形態変化を伴うものであるため、火花の根本部、中央部、及び、先端部の各部において火花の明度や密度等の特性には相違が生じる。そこで、本実施形態では、複数の火花画像から各々破裂密度を算出し、それらを平均した平均破裂密度を鋼材中の炭素量の確認に使用する。これにより、１枚の火花画像から求めた破裂密度のみで炭素量を確認した場合と比較して、予測の誤差を小さくすることが可能となる。
【００６４】
また、ＲＡＭ２１には、予め実測した、破裂密度、炭素量、及び、鋼種の関係がデータベースとして格納されている。鋼材識別処理部１８は、このデータベースを参照して、上記算出された破裂密度から測定対象の鋼材の鋼種を識別する。
【００６５】
〔実施例：撮影環境〕
図１５は、実施例における撮影環境を示す図である。カメラ５としては、ＤＩＴＥＣＴ社製のＨＡＳ−Ｌ２を用い、鋼材ＢとグラインダＣとの接触点から６００ｍｍの距離に設置した。
また、カメラ５によって、鋼材Ｂの軸線から側方に１１０ｍｍ〜４８０（＝１１０＋３７０）ｍｍの範囲における、鋼材Ｂを研削して発生する火花の画像を、３０ｆｐｓで撮像する。撮像した火花画像の解析度は、上下方向が２０４８ｐｉｘｅｌで左右方向が２５９２ｐｉｘｅｌである。解析時には、画像を５０枚用いると共に、各火花画像の上側及び下側を各１／４と、左右方向の火花の根本側を１／１０と、火花の先端側を１／４とをそれぞれカットして、火花以外の不要な部分を取り除いた。
【００６６】
画像の取得方法としては、複数枚の火花画像を取得する場合には所定の時間間隔で取得したり、火花発生後から所定の時間を経過した後に連続して複数枚の火花画像を取得したりする等の手法を用いればよい。
【００６７】
〔実施例：グレースケール化及び破裂密度の算出〕
この実施例では、ＲＧＢからグレースケールへの変換方法として、以下の３つの方法を行った。
１つ目は、火花画像中の各画素のＲＧＢ値を用いて、所定の変換式によってＹＩＱ表色系の輝度値Ｙへ変換したものである（図１６（Ａ））。この方法は、上述したとおり、一般的なグレースケールの手法である。
【００６８】
２つ目は、火花画像中の各画素のＲＧＢ値を、所定の変換式によってＸＹＺ表色系の輝度値Ｙへ変換したものである（図１６（Ｂ））。ここで、所定の変換式としては、下記の式を用いた。
【数４】
2020134204.tif 000006
３つ目は、火花画像中の各画素のＲＧＢ値を、ＲＧＢ値の最大値にしてＨＳＶ表色系のＶ（明度）値に変換したものである（図１６（Ｃ））。
【００６９】
また、各方法でグレースケール化した画像から二値化処理を行った画像を用いて、火花の破裂を認識し、破裂密度を算出した。その際、破裂の総数に対して除算する分母を、図１６（１）に示す全ての短直線数と、（２）に示す主幹短直線数と、（３）に示す主幹流線Ｒｍの総本数にして、破裂密度を算出した。
図１６に示す表は、火花画像を上述の３つの方法でグレースケール化したものについて、各々上述の３つのパターンによる破裂密度を算出した場合の破裂密度の相関係数を比較したものである。
この結果によれば、グレースケール化については、（Ｂ）に示すＸＹＺのＹに関するものが若干良い値が算出されたが、（Ａ）に示す輝度値（ＹＩＱのＹ）を用いるものとの差は殆どない。また、破裂密度の算出式の分母としては、（２）の主幹短直線数が良い値が算出された。
【００７０】
図１７は、鋼材の炭素含有量と新旧の手法によって実測した破裂密度とを比較したグラフである。ここで、特許文献１に記載の手法によって算出した破裂密度を、旧プログラムといい、上述した手法のうち、ＹＩＱの輝度値によりグレースケール化し、破裂密度を算出するための分母を主幹短直線数にしたもの（図１６の（Ａ）（２）の場合）を新プログラムという。
グラフには旧プログラムと新プログラムそれぞれについて、鋼材別の平均破裂密度を算出した。図１７からは、炭素量と破裂密度には相関があり、どれも炭素量にほぼ比例して破裂密度が増加していることから、鋼材の種類を定量的に分類することが可能であることがわかる。
また、旧プログラムの相関係数は、０．９１８４であった。この結果より、今回の手法の方が、相関係数が大きいものになり、より精度が高くなったことが分かった。
【００７１】
以上説明したように、本実施形態の鋼材識別装置１及びその方法によれば、以下のような効果がある。
（１）二値化後の火花画像に含む各流線Ｒを、流線Ｒの長さに応じて主幹流線Ｒｍと、枝流線Ｒｂとに分類し、主幹流線Ｒｍに分類された線上の任意の点を中心にした探索領域Ｅに少なくとも一部を含む枝流線Ｒｂの数がＫ以上の場合に、破裂として認識するようにした。よって、主幹流線Ｒｍと、枝流線Ｒｂとに基づいて、破裂をより正確に認識できる。
（２）探索領域Ｅに少なくとも一部を含む枝流線Ｒｂから、既に破裂として認識された枝流線Ｒｂを除いて破裂を認識するので、枝流線Ｒｂを重複して用いることがない。よって、実際の破裂に合った方法で破裂を認識できる。
（３）探索領域Ｅを、火花の根本から火花の先端に向かうほど大きな領域にする。よって、先端に行くほど、火花の広がりが大きくなる火花の特性に合った探索領域Ｅにすることができ、効率良く破裂を認識できる。
（４）探索領域Ｅに少なくとも一部を含む枝流線Ｒｂが、主幹流線Ｒｍに平行する場合に、その枝流線Ｒｂを除いて破裂を認識するので、主幹流線Ｒｍから生じないと思われる枝流線Ｒｂを除外できる。よって、実際の破裂に合った方法で破裂を認識できる。
（５）複数の角度に相当する短直線を表す複数のテンプレートを用いて、流線Ｒをマッチングするので、テンプレートを用いたマッチングによって、処理の効率化を図ることができる。また、マッチング数に基づいて、主幹流線Ｒｍと枝流線Ｒｂとを分類するので、主幹流線Ｒｍと枝流線Ｒｂとの分類を、簡単に行うことができる。
（６）認識した破裂の数と、主幹短直線数とによって、破裂密度を算出し、算出した破裂密度に基づいて鋼材Ｂ内の成分含有量を推定することで、鋼材Ｂを識別する。この方法を用いて認識した破裂数を処理した「破裂密度」が、鋼材の成分、特に炭素量との相関が高いことから、鋼材Ｂの成分評価を、定量的かつ容易に行うことができる。
【００７２】
以上、本発明の実施形態について説明したが、本発明は上述した実施形態に限定されるものではない。また、実施形態に記載した効果は、本発明から生じる最も好適な効果を列挙したに過ぎず、本発明による効果は、実施形態に記載したものに限定されない。なお、上述した実施形態及び後述する変形形態は、適宜組み合わせて用いることもできるが、詳細な説明は省略する。
【００７３】
（変形形態）
（１）本実施形態では、火花画像についてカット処理を行ったものを用いたが、これに限定されない。処理能力の高いコンピュータを用いれば、カット処理を行う必要がない。
（２）本実施形態においては、火花先端から根本に向かってマッチングをしたが、これに限定されるものではなく、火花根本から先端に向かってマッチングしてもよい。
【符号の説明】
【００７４】
１    鋼材識別装置
５    カメラ
１０  ＣＰＵ
１４  短直線マッチング処理部
１５  流線分類部
１６  破裂認識部
１７  破裂密度算出部
１８  鋼材識別処理部
２１  ＲＡＭ
２２  ＲＯＭ
１００  鋼材識別システム
Ｂ    鋼材
Ｃ    グラインダ
Ｅ    探索領域
Ｒ，Ｒ１，Ｒ２  流線
Ｒｍ  主幹流線
Ｒｂ  枝流線
【図１】
2020134204.tif 000007
【図２】
2020134204.tif 000008
【図３】
2020134204.tif 000009
【図４】
2020134204.tif 000010
【図５Ａ】
2020134204.tif 000011
【図５Ｂ】
2020134204.tif 000012
【図６】
2020134204.tif 000013
【図７】
2020134204.tif 000014
【図８】
2020134204.tif 000015
【図９】
2020134204.tif 000016
【図１０】
2020134204.tif 000017
【図１１Ａ】
2020134204.tif 000018
【図１１Ｂ】
2020134204.tif 000019
【図１１Ｃ】
2020134204.tif 000020
【図１２】
2020134204.tif 000021
【図１３】
2020134204.tif 000022
【図１４】
2020134204.tif 000023
【図１５】
2020134204.tif 000024
【図１６】
2020134204.tif 000025
【図１７】
2020134204.tif 000026