(57)【要約】【課題】高感度で安定した異物検出が容易に行えるよう被測定物のデュアルエネルギーＸ線画像に合わせて最適な差分処理用の重みパラメータを自動設定する。【解決手段】画像処理手段（６）は、デュアルエネルギーＸ線画像からローエネルギーとハイエネルギーの等価厚画像ペアを生成する等価厚画像生成手段（６１）と、前記等価厚画像ペアに対し混入異物の成分と被測定物の成分とに分離するように独立成分分析を適用して分離行列を求める分離行列算出手段（６３）と、前記分離行列の要素である２つの分離ベクトルに基いて前記差分処理の重みパラメータを求めるパラメータ算出手段（６４）と、前記重みパラメータを用いた前記差分処理によって前記等価厚画像ペアから前記混入異物の成分画像を分離する異成分画像分離手段（６７）と、前記混入異物の成分画像を閾値処理して異物検出する異物検出手段（６８）と、を具備している。 【選択図】図１
2010091483.tif 000002
【特許請求の範囲】
【請求項１】
被測定物のデュアルエネルギーＸ線画像を差分処理して混入異物の有無検出を行う異物検出方法であって、前記デュアルエネルギーＸ線画像としてローエネルギー等価厚画像とハイエネルギー等価厚画像の等価厚画像ペアを生成する段階と、前記等価厚画像ペアに対し前記混入異物の成分と前記被測定物の成分とに分離するように独立成分分析を適用して分離行列を求める段階と、前記分離行列の要素である２つの分離ベクトルに基いて前記差分処理の重みパラメータを求める段階と、前記重みパラメータを用いた前記差分処理によって前記等価厚画像ペアから前記混入異物の成分画像を分離する段階と、前記混入異物の成分画像を閾値処理して異物検出する段階と、を含む構成を有していることを特徴とする異物検出方法。
【請求項２】
前記等価厚画像ペアに対し前記混入異物の成分と前記被測定物の成分とに分離するように独立成分分析を適用して分離行列を求める段階の前に、前記等価厚画像ペアの画像ノイズを低減する段階を付加した請求項１に記載の異物検出方法。
【請求項３】
被測定物（３）にＸ線を照射するＸ線照射手段（２）と、前記被測定物を透過したＸ線強度をローエネルギーとハイエネルギーのＸ線画像ペアに変換するＸ線検出手段（４）と、前記Ｘ線画像ペアをデュアルエネルギーＸ線画像として記憶する画像入力手段（５）と、前記デュアルエネルギーＸ線画像を差分処理して被測定物に混入している異物の有無を検出する画像処理手段（６）と、を具備した異物検出装置において、前記画像処理手段（６）は、前記デュアルエネルギーＸ線画像からローエネルギー等価厚画像とハイエネルギー等価厚画像の等価厚画像ペアを生成する等価厚画像生成手段（６１）と、前記等価厚画像ペアに対し前記混入異物の成分と前記被測定物の成分とに分離するように独立成分分析を適用して分離行列を求める分離行列算出手段（６３）と、前記分離行列の要素である２つの分離ベクトルに基いて前記差分処理の重みパラメータを求めるパラメータ算出手段（６４）と、前記重みパラメータを用いた前記差分処理によって前記等価厚画像ペアから前記混入異物の成分画像を分離する異成分画像分離手段（６７）と、前記混入異物の成分画像を閾値処理して異物検出する異物検出手段（６８）と、を備えたことを特徴とする異物検出装置。
【請求項４】
前記等価厚画像ペアの画像ノイズを低減する画像ノイズ低減手段（６２）を備え、前記分離行列を求める手段は、ノイズ低減された等価厚画像ペアに対し分離行列を求めることを特徴とする請求項３に記載の異物検出装置。
【請求項５】
前記画像処理手段（６）に、前記パラメータ算出手段で求められた重みパラメータを記憶させるパラメータ記憶手段（６９）と、前記記憶させた重みパラメータを参照するパラメータ参照手段（６０）と、を付加した請求項３または４に記載の異物検出装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被測定物のデュアルエネルギーＸ線画像を差分処理して混入異物の有無検出を行う、異物検出方法および異物検出装置に関する。
【背景技術】
【０００２】
食品中の外来性異物や食肉中の残骨といった混入異物のＸ線検査を高コントラストな異成分画像を用い高感度に実施したい場合に有効な標準的方法として、ローエネルギーとハイエネルギーのデュアルエネルギーＸ線画像をローエネルギーとハイエネルギーの等価厚画像ペアに対数変換してから重み付き差分処理することにより異成分画像を分離するエネルギーサブトラクション（特許文献１参照）が知られている。
【０００３】
食品や食肉といった被測定物を透過したローエネルギーＸ線とハイエネルギーＸ線の強度は、間接変換型Ｘ線検出器を複数用いたエネルギー別フォトインテグレーションＸ線光学系（特許文献２参照）によって、または、Ｘ線光子到着毎に当該光子エネルギーを弁別する直接変換型Ｘ線検出器を単数用いたエネルギー別フォトカウンティングＸ線光学系によって、ローエネルギーとハイエネルギーのＸ線画像ペアすなわちデュアルエネルギーＸ線画像に変換できる。
【０００４】
前記エネルギー別フォトインテグレーションＸ線光学系の場合、当該Ｘ線光学系に依存して発生する前記デュアルエネルギーＸ線画像の相対位置ずれを無くすためイメージレジストレーション処理が必要になる。前記エネルギー別フォトカウンティングＸ線光学系の場合、当該相対位置ずれは発生しない。
【０００５】
混入異物を異成分として高コントラストに検出できるようにエネルギー帯域を適切に設定しておくとともに、前記重み付き差分処理の重みパラメータの設定が最適であれば、相対位置ずれの無いデュアルエネルギーＸ線画像に前記エネルギーサブトラクションを適用することによって直ちに高コントラストな異成分画像が分離できて、混入異物の有無判定が当該異成分画像に対する簡単な閾値処理で実現できる。
【０００６】
【特許文献１】特開平０４−１１４７３号公報
【特許文献２】特開２００２−３６５３６８号公報
【発明の開示】
【発明が解決しようとする課題】
【０００７】
前記エネルギーサブトラクションは食品中の外来性異物や食肉中の残骨といった混入異物のＸ線検査を高コントラストな異成分画像を用いて高感度に実施したい場合に効果を発揮する方法ではあるが，デュアルエネルギーＸ線画像の撮像条件や食品や食肉といった被測定物に依存する前記重み付き差分処理の重みパラメータを最適に設定する作業が微妙で煩わしいという欠点の克服が課題であった。
【０００８】
本発明の目的は、前記重みパラメータを自動設定することによって当該欠点を克服した異物検出方法および異物検出装置の提供にある。
【課題を解決するための手段】
【０００９】
本発明の発明者は、食品中の外来性異物や食肉中の残骨といった混入異物のＸ線吸収スペクトルと食品や食肉といった被測定物のＸ線吸収スペクトルが異なっている点に着眼して独立成分分析を行った結果、混入異物が被測定物の異成分画像として分離できるという事実を実験により見出していたが、この事実を応用すれば設定が煩わしいエネルギーサブトラクションにおける重みパラメータを最適に自動設定することができると考え、本発明に到達したものである。
【００１０】
本発明の請求項１に係わる異物検出方法は、被測定物のデュアルエネルギーＸ線画像を差分処理して混入異物の有無検出を行う異物検出方法であって、前記デュアルエネルギーＸ線画像としてローエネルギー等価厚画像とハイエネルギー等価厚画像の等価厚画像ペアを生成する段階と、前記等価厚画像ペアに対し前記混入異物の成分と前記被測定物の成分とに分離するように独立成分分析を適用して分離行列を求める段階と、前記分離行列の要素である２つの分離ベクトルに基いて前記差分処理の重みパラメータを求める段階と、前記重みパラメータを用いた前記差分処理によって前記等価厚画像ペアから前記混入異物の成分画像を分離する段階と、前記混入異物の成分画像を閾値処理して異物検出する段階と、を含む構成を有している。
【００１１】
この構成により、前記差分処理に用いる最適な重みパラメータを被測定物のデュアルエネルギーＸ線画像に合わせて自動設定することが可能となり高感度な異物検出ができる。
【００１２】
本発明の請求項２に係わる異物検出方法は、本発明の請求項１に係わる異物検出方法において、前記等価厚画像ペアに独立成分分析を適用して分離行列を求める段階の前に、前記等価厚画像ペアの画像ノイズを低減する段階を付加したものである。
すなわち被測定物のデュアルエネルギーＸ線画像を差分処理して混入異物の有無検出を行う異物検出方法であって、前記デュアルエネルギーＸ線画像としてローエネルギー等価厚画像とハイエネルギー等価厚画像の等価厚画像ペアを生成する段階と、前記等価厚画像ペアの画像ノイズを低減する段階と、前記画像ノイズを低減された等価厚画像ペアに対し前記混入異物の成分と前記被測定物の成分とに分離するように独立成分分析を適用して分離行列を求める段階と、前記分離行列の要素である２つの分離ベクトルに基いて前記差分処理の重みパラメータを求める段階と、前記重みパラメータを用いた前記差分処理によって前記等価厚画像ペアから前記混入異物の成分画像を分離する段階と、前記混入異物の成分画像を閾値処理して異物検出する段階と、を含む構成を有している。
【００１３】
この構成により、前記差分処理に用いる最適な重みパラメータを被測定物のノイズ低減されたデュアルエネルギーＸ線画像に合わせて自動設定することが可能となり高感度な異物検出ができる。
【００１４】
本発明の請求項３に係わる異物検出装置は、被測定物３にＸ線を照射するＸ線照射手段２と、前記被測定物を透過したＸ線強度をローエネルギーとハイエネルギーのＸ線画像ペアに変換するＸ線検出手段４と、前記Ｘ線画像ペアをデュアルエネルギーＸ線画像として記憶する画像入力手段５と、前記デュアルエネルギーＸ線画像を差分処理して被測定物に混入している異物の有無を検出する画像処理手段６と、を具備した異物検出装置において、画像処理手段６に請求項１の異物検出方法を用いている。
すなわち請求項１の異物検出方法を用いた前記画像処理手段６は、前記デュアルエネルギーＸ線画像からローエネルギー等価厚画像とハイエネルギー等価厚画像の等価厚画像ペアを生成する等価厚画像生成手段６１と、前記等価厚画像ペアに対し前記混入異物の成分と前記被測定物の成分とに分離するように独立成分分析を適用して分離行列を求める分離行列算出手段６３と、前記分離行列の要素である２つの分離ベクトルに基いて前記差分処理の重みパラメータを求めるパラメータ算出手段６４と、前記重みパラメータを用いた前記差分処理によって前記等価厚画像ペアから前記混入異物の成分画像を分離する異成分画像分離手段６７と、前記混入異物の成分画像を閾値処理して異物検出する異物検出手段６８と、を具備している。
【００１５】
この構成により、前記差分処理に用いる最適な重みパラメータを被測定物に合わせて自動設定することが可能となり高感度な異物検出装置が実現できる。
【００１６】
本発明の請求項４に係わる異物検出装置は、本発明の請求項３に係わる異物検出装置において、前記等価厚画像ペアの画像ノイズを低減する画像ノイズ低減手段６２を備え、前記分離行列を求める手段は、ノイズ低減された等価厚画像ペアに対し分離行列を求めるようにしたものである。
【００１７】
この構成により、前記差分処理に用いる最適な重みパラメータを被測定物のノイズ低減されたデュアルエネルギーＸ線画像に合わせて自動設定することが可能となり高感度な異物検出装置が実現できる。
【００１８】
本発明の請求項５に係わる異物検出装置は、本発明の請求項３または４に係わる異物検出装置において、前記画像処理手段６の中に、パラメータ算出手段６４で求められた重みパラメータを記憶させるパラメータ記憶手段６９と、前記記憶させた重みパラメータを参照するパラメータ参照手段６０と、を付加したものである。
【００１９】
この構成により、特定の被測定物に合わせて求めらた前記差分処理に用いる最適な重みパラメータを記憶させておき、同種の被測定物に対しても毎回行っていた分離行列算出から重みパラメータ算出にいたる計算処理に替えて、当該重みパラメータを参照させることによって、最適な重みパラメータを求めた特定の被測定物と同種の被測定物を高速に検査できるようにした高感度な異物検出装置が実現できる。
【発明の効果】
【００２０】
本発明の請求項１に係わる異物検出方法は、前記差分処理に用いる最適な重みパラメータを被測定物のデュアルエネルギーＸ線画像に合わせて自動設定することが可能となり、高感度で安定した異物検出を行うことができる。
【００２１】
本発明の請求項２に係わる異物検出方法は、前記差分処理に用いる最適な重みパラメータを被測定物のノイズ低減されたデュアルエネルギーＸ線画像に合わせて自動設定することが可能となり、格段に高感度で安定した異物検出を行うことができる。
【００２２】
本発明の請求項３または４に係わる異物検出装置は、前記差分処理に用いる最適な重みパラメータを被測定物に合わせて自動設定することが可能となるため、高感度で安定した異物検出を行うことができる。
【００２３】
本発明の請求項５に係わる異物検出装置は、特定の被測定物に合わせて求められた前記差分処理に用いる最適な重みパラメータを記憶させておき、当該重みパラメータを参照させることによって、同種の被測定物を高速に検査できるため、高感度で高速な異物検出を行うことができる。
【発明を実施するための最良の形態】
【００２４】
以下、本発明の好ましい実施の形態について図面を参照しつつ説明する。
【００２５】
１．第１実施形態（図１〜図２）
本実施形態の異物検出装置は図１に示すように、被測定物３を搬送する搬送手段１と、被測定物３にＸ線を照射するＸ線照射手段（Ｘ線源）２と、前記被測定物を透過したＸ線強度をローエネルギーとハイエネルギーのＸ線画像ペアに変換するＸ線検出手段４と、前記Ｘ線画像ペアをデュアルエネルギーＸ線画像として記憶する画像入力手段５と、前記デュアルエネルギーＸ線画像を差分処理して被測定物に混入している異物の有無を検出する画像処理手段６とを備え、必要に応じて画像表示手段７を備えている。
【００２６】
搬送手段１は、例えばＸ線を良く透過するベルトコンベアで実現され、対向配置されたＸ線源２とＸ線検出手段４の間を通して被測定物３を搬送する。Ｘ線源２から照射されたＸ線は、被測定物３による吸収と前記ベルトコンベアによる僅かな吸収を受けてこれらを透過した後、Ｘ線検出手段４に到達する。
【００２７】
Ｘ線検出手段４は、例えばＸ線光子到着毎に当該光子エネルギーを弁別する直接変換型Ｘ線ラインセンサで実現され、エネルギー別フォトカウンティングＸ線光学系を成すことによって、前記Ｘ線ラインセンサによる１ライン上のサンプリングピッチと略等しいサンプリングピッチで搬送方向にサンプリングすることによって被測定物３と前記ベルトコンベアを透過したＸ線強度をエネルギー帯域別の画像データとしてメモリー上に取り込み、ローエネルギーとハイエネルギーのＸ線画像ペアに変換する。
【００２８】
なお図示しないが、Ｘ線検出手段４として前記直接変換型Ｘ線ラインセンサに替えてローエネルギー用とハイエネルギー用に２台の間接変換型Ｘ線ラインセンサを用いることで、エネルギー別フォトインテグレーションＸ線光学系を成しても良い。
【００２９】
画像入力手段５は、被測定物３の前記Ｘ線画像ペアを相対位置ずれの無いデュアルエネルギーＸ線画像として記憶し、画像処理手段６に出力する。なお前記エネルギー別フォトインテグレーションＸ線光学系を成した場合は、メモリー上に取り込まれたＸ線画像ペアの相対位置ずれを無くすため、ここでイメージレジストレーション処理を行う。
【００３０】
画像処理手段６は画像処理プログラム等が実装されたＣＰＵ等で、図２に示す本発明に係わる異物検出方法が実装されており、画像入力手段５から出力された相対位置ずれの無いデュアルエネルギーＸ線画像の入力を受けて（ステップＳ１０１）、等価厚画像生成手段６１で対数変換を行うことによってローエネルギー等価厚画像とハイエネルギー等価厚画像の等価厚画像ペアを生成し（ステップＳ１０２）、必要に応じて画像ノイズ低減手段６２で前記等価厚画像ペアにエッジ保存平滑化処理等を適用して画像ノイズを低減し（ステップＳ１０３）、分離行列計算手段６３で当該等価厚画像ペアに独立成分分析を適用して分離行列を求め（ステップＳ１０４）、分離ベクトル判別手段６４１でローエネルギーのハイエネルギーからの差分に近い方の分離ベクトルを異成分の分離ベクトルとして判別し（ステップＳ１０５）、分離ベクトル正規化手段６４２で当該異成分の分離ベクトルをローエネルギーからスケーリングされたハイエネルギーが差し引かれる形に正規化し（ステップＳ１０６）、パラメータ設定手段６４３で前記正規化された分離ベクトルを前記差分処理の重みパラメータとして自動設定する（ステップＳ１０７）。すなわち手段６４１〜６４３（ステップＳ１０５〜Ｓ１０７）から構成されるパラメータ算出手段６４によって前記分離行列の要素である２つの分離ベクトルに基いて前記差分処理の重みパラメータを求めて自動設定する。さらに、前記自動設定された重みパラメータを用いて異成分画像分離手段６７で前記重みパラメータを用いた前記差分処理によって前記等価厚画像ペアから前記混入異物の成分画像を分離し（ステップＳ１０８）、異物検出手段６８で前記混入異物の成分画像を閾値処理して異物検出する（ステップＳ１０９）ことによって混入異物の有無検出を行う。
【００３１】
前記独立成分分析としてはFastICA（Fast
Independent Component Analysis）やSOBI（Second Order Blind Identification）といった一般的で公知な信号分離アルゴリズムを適用することが可能である。独立成分分析の概略を図５に示す。図５（a）はローエネルギー等価厚画像を示し、図５（ｂ）はハイエネルギー等価厚画像を示し、図５（ｃ）は混入異物の成分画像を示し、図５（ｄ）は被測定物の成分画像を示しており、独立成分分析は前記等価厚画像ペアから前記混入異物の成分と前記被測定物の成分とに分離するように分離行列を求めるものである。なお、前記独立成分分析の入力となる等価厚画像ペアには少なくとも被測定物３の主要部分が含まれていることが必要であり、分離対象となる異成分異物は含まれていても含まれていなくても良い。なぜなら、例えばソーセージの異物検査のように、被測定物を略単一成分とみなすことが可能であれば、前記独立成分分析により独立成分として分離できる２種類の信号は被測定物（ソーセージ）の信号とそれ以外（残骨、金属、ガラス、石などの硬質異物）の信号であるとみなすことができるからである。
【００３２】
ステップＳ１０４で計算された分離行列の要素である２つの分離ベクトルを（ａ，ｂ）、（ｃ，ｄ）とすれば、前記等価厚画像ペアの対応画素における濃度ペアを成分としたベクトル（Ｌ，Ｈ）に対して、２つの独立成分は  Ｘ＝ａＬ＋ｂＨ、Ｙ＝ｃＬ＋ｄＨ  で与えられるが、このままでは、ＸとＹのどちらが異成分画像の濃度に相当するかは不定であり、かつ、ＸとＹの濃度倍率も反転状態を含み不定であり、最適な重みパラメータによるエネルギーサブトラクションに結びつけることができない。
【００３３】
そこで、ステップＳ１０５では、Ｌをローエネルギー成分、Ｈをハイエネルギー成分、正値の重みパラメータをｋとして、前記差分処理に該当するところのＺ＝Ｌ−ｋＨ  なるエネルギーサブトラクションが分離ベクトル（１，−ｋ）によって分離できる成分に等しいことに着目し、（ａ，ｂ）と（ｃ，ｄ）のうち（１，−ｋ）に近い方を異成分の分離ベクトルとして判別する。具体的には、ローエネルギー等価厚画像はハイエネルギー等価厚画像よりも高コントラストかつ鮮明になるため前記等価厚画像ペアの濃度が同程度になるように前記対数変換しておけばｋは１に近い値になるという知見に基き、例えば  ｂ／ａ  と  ｄ／ｃ  の計算値を評価して−1に近い値を算出した方の分離ベクトルを異成分の分離ベクトルと判別すれば良い。
【００３４】
そして、ステップＳ１０６では、前記判別された分離ベクトルを（ａ，ｂ）として、（ａ，ｂ）に倍率  １／ａ  を掛けて、（１，ｂ／ａ）という形に変形することによって当該分離ベクトルを正規化する。
【００３５】
そして、ステップＳ１０７では、前記正規化された分離ベクトルを（１，ｂ／ａ）として、Ｚ＝Ｌ−ｋＨ  なる前記差分処理に対して、ｋ＝−ｂ／ａ  と自動設定する。
【００３６】
そして、ステップＳ１０８では、Ｚ＝Ｌ−ｋＨ  なる前記差分処理を実行する。
【００３７】
画像表示手段７は、異物検出手段６８の異物検出結果に対応する検査結果情報をフラットパネルディスプレイ等に表示する。
【００３８】
２．第２実施形態（図３〜図４）
本実施の形態は、第１実施形態の画像処理手段６において、自動設定された重みパラメータを記憶させるためのパラメータ記憶手段６９と、前記記憶させた重みパラメータを参照させるためのパラメータ参照手段６０と、を付加したものであり、その他の構成に関しては第１実施形態と同様である。それゆえ、パラメータ記憶手段６９とパラメータ参照手段６０について説明し、第１実施形態と同様の部分は図中に第１実施形態と同じ符号を付して第１実施形態の説明を援用することによって、再度の説明は省略する。
【００３９】
パラメータ記憶手段６９は、パラメータ算出手段６４で自動設定した前記差分処理の重みパラメータを受けて異成分画像分離手段６７に出力するとともに、図示しない運転モード信号がパラメータ記憶を示しているとき、前記パラメータを図示しない参照用メモリーに記憶する（ステップＳ２０２）。
【００４０】
パラメータ参照手段６０は、前記運転モード信号がパラメータ参照を示していないとき、等価厚画像生成手段６１が出力した、画像ノイズ低減手段６２が存在する場合はさらに画像ノイズ低減手段６２を経由した、等価厚画像ペアを受けて分離行列算出手段６３に出力する。そして前記運転モード信号がパラメータ参照を示しているとき、前述したステップＳ１０４からステップＳ２０２までの流れを切り替えて、記憶していた前記差分処理の重みパラメータを前記参照用メモリーから読み出して異成分画像分離手段６７に直接出力する（ステップＳ２０１）。
【産業上の利用可能性】
【００４１】
以上説明したように本発明は、高感度で安定した異物検出が容易に行えるよう被測定物のデュアルエネルギーＸ線画像に合わせて最適なエネルギーサブトラクション用の重みパラメータを自動設定するため、デュアルエネルギーＸ線データを出力する検出手段を備えた異物検出装置全般に有用である。
【図面の簡単な説明】
【００４２】
【図１】本発明の第１実施形態に係わる異物検出装置の概略構成を示す図である。
【図２】本発明の第１実施形態に係わる異物検出方法の概略の流れを示すフローチャートである。
【図３】本発明の第２実施形態に係わる異物検出装置の概略構成を示す図である。
【図４】本発明の第２実施形態に係わる異物検出方法の概略の流れを示すフローチャートである。
【図５】一般的で公知な独立成分分析の概略を示す図である。
【符号の説明】
【００４３】
１…搬送手段
２…Ｘ線源
３…被測定物
４…Ｘ線検出手段
５…画像入力手段
６…画像処理手段
７…画像表示手段 
６０…パラメータ参照手段
６１…等価厚画像生成手段
６２…画像ノイズ低減手段
６３…分離行列算出手段
６４…パラメータ算出手段
６４１…分離ベクトル判別手段
６４２…分離ベクトル正規化手段
６４３…パラメータ設定手段
６７…異成分画像分離手段
６８…異物検出手段
６９…パラメータ記憶手段
【図１】
2010091483.tif 000003
【図２】
2010091483.tif 000004
【図３】
2010091483.tif 000005
【図４】
2010091483.tif 000006
【図５】
2010091483.tif 000007