(57)【要約】【課題】検出精度を低下させることなく、飛跡を検出するための演算量を削減することができる飛跡検出装置、飛跡検出方法、及び飛跡検出プログラムを提供する。【解決手段】飛跡検出装置１２は、乳剤層５１における深さ方向が異なる複数の断面の各々を撮像して得られ、かつμ粒子の入射側に最も近い乳剤層５１の断面によって得られた断層画像５７１と、μ粒子の入射側に最も遠い乳剤層５１の断面によって得られた断層画像５７ｎとを含み、各々複数の画素５８が配置された断層画像群５６を取得する断層画像取得部２０を備える。また、飛跡検出装置１２は、断層画像５７１の基準点６０を通る方向が各々異なり、かつ各々が断層画像５７ｎの異なる画素５８を通る複数の飛跡候補線Ｃを特定する飛跡候補線特定部２４と、乳剤層５１における飛跡Ｔを、複数の飛跡候補線Ｃから検出する検出部２６と、を備える。【選択図】図５
2020165918.tif 000002
【特許請求の範囲】
【請求項１】
入射された荷電粒子が透過することにより、前記荷電粒子の飛跡を潜像として記録する荷電粒子飛跡記録層における、前記飛跡を検出する飛跡検出装置であって、
前記荷電粒子飛跡記録層における深さ方向が異なる複数の断面の各々を撮像して得られ、かつ前記荷電粒子の入射側に最も近い前記荷電粒子飛跡記録層の断面によって得られた第１断層画像と、前記荷電粒子の入射側に最も遠い前記荷電粒子飛跡記録層の断面によって得られた第２断層画像とを含み、各々複数の画素が配置された複数の断層画像を取得する断層画像取得部と、
前記第１断層画像の基準点を通る方向が各々異なり、かつ各々が前記第２断層画像の異なる画素を通る複数の飛跡候補線を特定する飛跡候補線特定部と、
前記荷電粒子飛跡記録層における前記飛跡を、前記複数の飛跡候補線から検出する検出部と、
を備えた飛跡検出装置。
【請求項２】
前記飛跡候補線特定部は、前記基準点を通る方向が各々異なり、各々が前記第２断層画像を通る複数の仮想飛跡候補線から、前記飛跡候補線を選択する、
請求項１に記載の飛跡検出装置。
【請求項３】
前記飛跡候補線特定部は、前記複数の仮想飛跡候補線のうち、前記第２断層画像における同一の画素を通る仮想飛跡候補線が複数ある場合、いずれか一つの仮想飛跡候補線のみを前記複数の飛跡候補線に含め、
前記第２断層画像における同一の画素を通る仮想飛跡候補線が一つの場合、前記同一の画素を通る仮想飛跡候補線を前記複数の飛跡候補線に含める、
請求項２に記載の飛跡検出装置。
【請求項４】
前記飛跡候補線特定部は、前記第２断層画像における同一の画素を通る仮想飛跡候補線が複数ある場合、前記同一の画素を通る複数の仮想飛跡候補線のうち、前記同一の画素の中央に最も近い位置を通る仮想飛跡候補線を、前記一つの仮想飛跡候補線とする、
請求項３に記載の飛跡検出装置。
【請求項５】
前記飛跡候補線特定部は、前記複数の飛跡候補線同士が成す角度が、予め定められた角度以下の場合、少なくとも１つの飛跡候補線を、特定した複数の飛跡候補線から除外する、
請求項１から請求項４のいずれか１項に記載の飛跡検出装置。
【請求項６】
前記荷電粒子飛跡記録層の荷電粒子が入射される面に平行な面内角が大きいとした条件、及び小さいとした条件のいずれか一方の条件が削除条件として予め定められており、
前記飛跡候補線特定部は、前記複数の飛跡候補線同士のうち、前記削除条件を満たす飛跡候補線を、特定した複数の飛跡候補線から除外する、
請求項５に記載の飛跡検出装置。
【請求項７】
前記飛跡候補線特定部は、天頂角が同一の複数の飛跡候補線について、互いの成す角度が前記予め定められた角度以下となる複数の飛跡候補線同士が存在する場合、特定した複数の飛跡候補線のうちの一部を除外することにより、残存する複数の飛跡候補線の分布を、天頂角が０度の軸に対して、回転対称に最も近くなる状態とする、
請求項５に記載の飛跡検出装置。
【請求項８】
前記検出部は、予め定められた数の隣接する断層画像同士を組み合わせて得られた複数の合成断層画像に含まれる前記潜像に応じた飛跡像に基づいて、前記飛跡候補線から前記飛跡を検出する、
請求項１から請求項７のいずれか１項に記載の飛跡検出装置。
【請求項９】
前記合成断層画像は、前記飛跡候補線が前記荷電粒子飛跡記録層に入射する方向に応じて、前記隣接する断層画像同士の位置をずらして合成した画像である、
請求項８に記載の飛跡検出装置。
【請求項１０】
入射された荷電粒子が透過することにより、前記荷電粒子の飛跡を潜像として記録する荷電粒子飛跡記録層における、前記飛跡を検出する飛跡検出方法であって、
前記荷電粒子飛跡記録層における深さ方向が異なる複数の断面の各々を撮像して得られ、かつ前記荷電粒子の入射側に最も近い前記荷電粒子飛跡記録層の断面によって得られた第１断層画像と、前記荷電粒子の入射側に最も遠い前記荷電粒子飛跡記録層の断面によって得られた第２断層画像とを含み、各々複数の画素が配置された複数の断層画像を取得し、
前記第１断層画像の基準点を通る方向が各々異なり、かつ各々が前記第２断層画像の異なる画素を通る複数の飛跡候補線を特定し、
前記荷電粒子飛跡記録層における前記飛跡を、前記複数の飛跡候補線から検出する、
処理を含む飛跡検出方法。
【請求項１１】
入射された荷電粒子が透過することにより、前記荷電粒子の飛跡を潜像として記録する荷電粒子飛跡記録層における、前記飛跡を検出するコンピュータに、
前記荷電粒子飛跡記録層における深さ方向が異なる複数の断面の各々を撮像して得られ、かつ前記荷電粒子の入射側に最も近い前記荷電粒子飛跡記録層の断面によって得られた第１断層画像と、前記荷電粒子の入射側に最も遠い前記荷電粒子飛跡記録層の断面によって得られた第２断層画像とを含み、各々複数の画素が配置された複数の断層画像を取得し、
前記第１断層画像の基準点を通る方向が各々異なり、かつ各々が前記第２断層画像の異なる画素を通る複数の飛跡候補線を特定し、
前記荷電粒子飛跡記録層における前記飛跡を、前記複数の飛跡候補線から検出する、
処理を実行させるための飛跡検出プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、飛跡検出装置、飛跡検出方法、及び飛跡検出プログラムに関する。
【背景技術】
【０００２】
従来、非破壊検査等において、荷電粒子の飛跡を潜像として記録する荷電粒子飛跡記録層を含む原子核乾板を用いて、宇宙線に含まれるミュー粒子等の荷電粒子の飛跡を検出することが行われている（特許文献１参照）。
【０００３】
特許文献１に記載の技術では、原子核乾板を現像後、焦点面の深さを変えながら原子核乾板の荷電粒子飛跡記録層を連続的に撮像して、複数の断層画像を取得する。そして、特許文献１に記載の技術では、複数の断層画像を飛跡の角度に応じたシフト量でシフトさせるいわゆるシフト法により、断層画像を重ね合わせて得られた飛跡の存在を示す指標の分布に基づいて、飛跡を検出する。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１０−１０１６７６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
上述のシフト法では、荷電粒子の飛跡として想定される飛跡の飛跡角度に応じたシフト量で各断層画像をずらし、該当する飛跡角度を有する飛跡が、各断層画像上のほぼ同じｘ、ｙ座標に位置する状態として、各断層画像の同一座標における黒点の有無に基づいて、飛跡を検出する。
【０００６】
しかしながら、想定される飛跡の飛跡角度全てに対して、上述のシフト法を適用した場合、演算量が膨大となる場合があった。一方、単純に想定される飛跡の飛跡角度を削減した場合、飛跡検出精度が低下する懸念があった。
【０００７】
本開示は、以上の事情を鑑みてなされたものであり、検出精度を低下させることなく、飛跡を検出するための演算量を削減することができる飛跡検出装置、飛跡検出方法、及び飛跡検出プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記目的を達成するために、本開示の第１の態様の飛跡検出装置は、入射された荷電粒子が透過することにより、荷電粒子の飛跡を潜像として記録する荷電粒子飛跡記録層における、飛跡を検出する飛跡検出装置であって、荷電粒子飛跡記録層における深さ方向が異なる複数の断面の各々を撮像して得られ、かつ荷電粒子の入射側に最も近い荷電粒子飛跡記録層の断面によって得られた第１断層画像と、荷電粒子の入射側に最も遠い荷電粒子飛跡記録層の断面によって得られた第２断層画像とを含み、各々複数の画素が配置された複数の断層画像を取得する断層画像取得部と、第１断層画像の基準点を通る方向が各々異なり、かつ各々が第２断層画像の異なる画素を通る複数の飛跡候補線を特定する飛跡候補線特定部と、荷電粒子飛跡記録層における飛跡を、複数の飛跡候補線から検出する検出部と、を備える。
【０００９】
第１の態様の飛跡検出装置によれば、検出精度を低下させることなく、飛跡を検出するための演算量を削減することができる。
【００１０】
なお、本開示の第２の態様の飛跡検出装置は、第１の態様の飛跡検出装置において、飛跡候補線特定部は、基準点を通る方向が各々異なり、各々が第２断層画像を通る複数の仮想飛跡候補線から、飛跡候補線を選択する。
【００１１】
第２の態様の飛跡検出装置によれば、飛跡候補線の特定を容易に行うことがきる。
【００１２】
また、本開示の第３の態様の飛跡検出装置は、第２の態様の飛跡検出装置において、飛跡候補線特定部は、複数の仮想飛跡候補線のうち、第２断層画像における同一の画素を通る仮想飛跡候補線が複数ある場合、いずれか一つの仮想飛跡候補線のみを複数の飛跡候補線に含め、第２断層画像における同一の画素を通る仮想飛跡候補線が一つの場合、同一の画素を通る仮想飛跡候補線を複数の飛跡候補線に含める。
【００１３】
第３の態様の飛跡検出装置によれば、検出精度を低下させることなく、飛跡を検出するための演算量をより適切に削減することができる。
【００１４】
また、本開示の第４の態様の飛跡検出装置は、第３の態様の飛跡検出装置において、飛跡候補線特定部は、第２断層画像における同一の画素を通る仮想飛跡候補線が複数ある場合、同一の画素を通る複数の仮想飛跡候補線のうち、同一の画素の中央に最も近い位置を通る仮想飛跡候補線を、一つの仮想飛跡候補線とする。
【００１５】
第４の態様の飛跡検出装置によれば、検出精度を向上させることができる。
【００１６】
また、本開示の第５の態様の飛跡検出装置は、第１の態様から第４の態様のいずれか１態様の飛跡検出装置において、飛跡候補線特定部は、複数の飛跡候補線同士が成す角度が、予め定められた角度以下の場合、少なくとも１つの飛跡候補線を、特定した複数の飛跡候補線から除外する。
【００１７】
第５の態様の飛跡検出装置によれば、検出精度を低下させることなく、飛跡を検出するための演算量をより削減することができる。
【００１８】
また、本開示の第６の態様の飛跡検出装置は、第５の態様の飛跡検出装置において、荷電粒子飛跡記録層の荷電粒子が入射される面に平行な面内角が大きいとした条件、及び小さいとした条件のいずれか一方の条件が削除条件として予め定められており、飛跡候補線特定部は、複数の飛跡候補線同士のうち、削除条件を満たす飛跡候補線を、特定した複数の飛跡候補線から除外する。
【００１９】
第６の態様の飛跡検出装置によれば、特定された飛跡候補線を間引くことにより、規則的に並ぶ飛跡候補線を特定することができる。
【００２０】
また、本開示の第７の態様の飛跡検出装置は、第５の態様の飛跡検出装置において、飛跡候補線特定部は、天頂角が同一の複数の飛跡候補線について、互いの成す角度が予め定められた角度以下となる複数の飛跡候補線同士が存在する場合、特定した複数の仮想飛跡候補線のうちの一部を除外することにより、残存する複数の飛跡候補線の分布を、天頂角が０度の軸に対して、回転対称に最も近くなる状態とする。
【００２１】
第７の態様の飛跡検出装置によれば、特定された飛跡候補線を間引くことにより、より規則的に並ぶ飛跡候補線を特定することができる。
【００２２】
また、本開示の第８の態様の飛跡検出装置は、第１の態様から第７の態様のいずれか１態様の飛跡検出装置において、検出部は、予め定められた数の隣接する断層画像同士を組み合わせて得られた複数の合成断層画像に含まれる潜像に応じた飛跡像に基づいて、飛跡候補線から飛跡を検出する。
【００２３】
第８の態様の飛跡検出装置によれば、検出精度を低下させることなく、飛跡を検出するための演算量をさらに削減することができる。
【００２４】
また、本開示の第９の態様の飛跡検出装置は、第８の態様の飛跡検出装置において、合成断層画像は、飛跡候補線が荷電粒子飛跡記録層に入射する方向に応じて、隣接する断層画像同士の位置をずらして合成した画像である。
【００２５】
第９の態様の飛跡検出装置によれば、検出精度をより向上させることができる。
【００２６】
また、本開示の第１０の態様の飛跡検出方法は、入射された荷電粒子が透過することにより、荷電粒子の飛跡を潜像として記録する荷電粒子飛跡記録層における、飛跡を検出する飛跡検出方法であって、荷電粒子飛跡記録層における深さ方向が異なる複数の断面の各々を撮像して得られ、かつ荷電粒子の入射側に最も近い荷電粒子飛跡記録層の断面によって得られた第１断層画像と、荷電粒子の入射側に最も遠い荷電粒子飛跡記録層の断面によって得られた第２断層画像とを含み、各々複数の画素が配置された複数の断層画像を取得し、第１断層画像の基準点を通る方向が各々異なり、かつ各々が第２断層画像の異なる画素を通る複数の飛跡候補線を特定し、荷電粒子飛跡記録層における飛跡を、複数の飛跡候補線から検出する、処理を含む。
【００２７】
また、本開示の第１１の態様の飛跡検出プログラムは、入射された荷電粒子が透過することにより、荷電粒子の飛跡を潜像として記録する荷電粒子飛跡記録層における、飛跡を検出するコンピュータに、荷電粒子飛跡記録層における深さ方向が異なる複数の断面の各々を撮像して得られ、かつ荷電粒子の入射側に最も近い荷電粒子飛跡記録層の断面によって得られた第１断層画像と、荷電粒子の入射側に最も遠い荷電粒子飛跡記録層の断面によって得られた第２断層画像とを含み、各々複数の画素が配置された複数の断層画像を取得し、第１断層画像の基準点を通る方向が各々異なり、かつ各々が第２断層画像の異なる画素を通る複数の飛跡候補線を特定し、荷電粒子飛跡記録層における飛跡を、複数の飛跡候補線から検出する、処理を実行させるためのものである。
【００２８】
また、本開示の飛跡検出装置は、プロセッサ及びメモリを有し、入射された荷電粒子が透過することにより、荷電粒子の飛跡を潜像として記録する荷電粒子飛跡記録層における、飛跡を検出する飛跡検出装置であって、プロセッサが、荷電粒子飛跡記録層における深さ方向が異なる複数の断面の各々を撮像して得られ、かつ荷電粒子の入射側に最も近い荷電粒子飛跡記録層の断面によって得られた第１断層画像と、荷電粒子の入射側に最も遠い荷電粒子飛跡記録層の断面によって得られた第２断層画像とを含み、各々複数の画素が配置された複数の断層画像を取得し、第１断層画像の基準点を通る方向が各々異なり、かつ各々が第２断層画像の異なる画素を通る複数の飛跡候補線を特定し、荷電粒子飛跡記録層における飛跡を、複数の飛跡候補線から検出する。
【発明の効果】
【００２９】
本開示によれば、検出精度を低下させることなく、飛跡を検出するための演算量を削減することができる。
【図面の簡単な説明】
【００３０】
【図１】実施形態の原子核乾板の一例の側面断面図である。
【図２】実施形態の非破壊検査システムの構成の一例を表すブロック図である。
【図３】実施形態の原子核乾板から生成される断層画像を説明する図である。
【図４】実施形態の飛跡検出装置のハードウェア構成の一例を示すブロック図である。
【図５】実施形態の飛跡検出装置の機能的な構成の一例を示すブロック図である。
【図６】飛跡の飛跡角度を説明するための図である。
【図７】乳剤層によって得られる断層画像群に含まれる、複数の断層画像を通る飛跡を説明するための図である。
【図８】基準点を通る仮想飛跡候補線の例を説明するための図である。
【図９】基準点を通る飛跡候補線の例を説明するための図である。
【図１０】第１実施形態の飛跡検出装置で実行される飛跡検出処理の一例を示すフローチャートである。
【図１１】図１０に示した飛跡検出処理にて実行される飛跡候補線特定処理の一例を示すフローチャートである。
【図１２】図１１に示した飛跡候補線特定処理にて実行される飛跡候補角度導出処理の一例を示すフローチャートである。
【図１３】図１２に示した飛跡候補角度導出処理にて用いられる一覧表の一例を表す図である。
【図１４】図１２に示した飛跡候補角度導出処理にて用いられるＩＤ配列表の一例を表す図である。
【図１５】図１０に示した飛跡検出処理にて実行される検出処理の一例を示すフローチャートである。
【図１６Ａ】第１実施形態の検出処理におけるシフト法による飛跡の検出について説明するための図である。
【図１６Ｂ】第１実施形態の検出処理におけるシフト法による飛跡の検出について説明するための図である。
【図１７】第２実施形態の飛跡検出装置において実行される飛跡候補角度導出処理の一例を示すフローチャートである。
【図１８Ａ】第２実施形態の飛跡候補角度処理における飛跡候補線の除外（間引き）の例を説明するための図である。
【図１８Ｂ】第２実施形態の飛跡候補角度処理における飛跡候補線の除外（間引き）の例を説明するための図である。
【図１９】断層画像の合成について説明するための図である。
【図２０】断層画像の合成について説明するための図である。
【図２１】第３実施形態の飛跡検出処理において実行される検出処理の一例を示すフローチャートである。
【図２２】断層画像を合成して得られる合成断層画像を説明するための図である。
【図２３Ａ】第３実施形態の検出処理におけるシフト法による飛跡の検出について説明するための図である。
【図２３Ｂ】第３実施形態の検出処理におけるシフト法による飛跡の検出について説明するための図である。
【発明を実施するための形態】
【００３１】
以下、図面を参照して、本開示の技術を実施するための形態例を詳細に説明する。
【００３２】
［第１実施形態］
原子核乾板５０を用いた非破壊検査では、検査対象物の周辺に原子核乾板５０を設置し、設置中に原子核乾板５０に記録されたμ粒子の飛跡Ｔを検出し、検出した飛跡Ｔを解析することによって検査対象物の密度等が得られる。
【００３３】
図１には、本実施形態の原子核乾板５０の一例の側面断面図を示す。図１に示すように、本実施形態の原子核乾板５０は、２つの乳剤層５１Ａ及び５１Ｂ、及び１つのフィルムベース５２を備え、フィルムベース５２の両面に、乳剤層５１Ａ及び５１Ｂが塗布されている。フィルムベース５２は、例えば、プラスチック製であり、厚さが３００μｍとされる。乳剤層５１Ａ及び５１Ｂの各々は、例えば、ハロゲン化銀を含み、厚さが６０μｍとされる。本実施形態の乳剤層５１Ａ及び５１Ｂが、本開示における、入射された荷電粒子が透過することにより、荷電粒子の飛跡を潜像として記録する荷電粒子飛跡記録層の一例である。なお、飛跡記録層に含まれるハロゲン化銀以外の飛跡記録材料として、例えば、シュウ酸第二鉄塩を含む鉄塩感光材料、又は芳香族ジアゾニウム塩を含むジアゾ感光材を用いることも可能である。なお、本実施形態では、乳剤層５１Ａ及び５１Ｂの各々を区別せずに総称する場合、「乳剤層５１」という。
【００３４】
宇宙線に含まれるμ粒子が原子核乾板５０を透過すると、乳剤層５１内での電離作用によって発生した電子により、乳剤層５１内にμ粒子の飛跡Ｔに沿った潜像が記録される。図１の例では、乳剤層５１Ａ側からμ粒子が入射され、入射されたμ粒子によって飛跡Ｔが形成される。なお、図１の例では、飛跡Ｔのうち、乳剤層５１内に潜像として記録される部分を実線で示し、乳剤層５１内に潜像として記録される部分以外の部分を破線で示している。本実施形態のμ粒子が開示の技術に係る荷電粒子の一例である。なお、荷電粒子は、μ粒子に限定されない。例えば、荷電粒子は、γ線又は中性子線等から２次的に発生した電子又は陽子でもよい。また、荷電粒子は、原子核乾板５０を直接感光させるものに限定されず、例えば、原子核乾板５０の飛跡記録材料中の陽子と衝突し、散乱させることにより陽子が検出されることで飛跡が形成されるものでもよい。
【００３５】
まず、原子核乾板５０を用いた、本実施形態の非破壊検査システム１の構成について説明する。図２には、本実施形態の非破壊検査システム１の構成の一例を表すブロック図を示す。図２に示すように非破壊検査システム１は、断層画像生成装置１０、飛跡検出装置１２、及び解析装置１４を備える。断層画像生成装置１０、飛跡検出装置１２、及び解析装置１４は、それぞれネットワークＮに接続され、ネットワークＮを介して互いに通信が可能とされる。なお、本実施形態では、断層画像生成装置１０、飛跡検出装置１２、及び解析装置１４の各々を別個の装置とした形態について示したが、これらのうち、２つまたは全部を組み合わせて１つの装置とした形態であってもよいことはいうまでもない。
【００３６】
断層画像生成装置１０は、μ粒子によって、潜像が形成された原子核乾板５０の乳剤層５１の断層画像を生成する装置である。上述のようにして潜像が形成された原子核乾板５０の現像処理を行うことで、潜像が、１μｍ程度の大きさの銀粒子（黒点）として原子核乾板５０中に残留する。原子核乾板５０を通過したμ粒子の飛跡Ｔは、銀粒子（黒点）の立体的な並びとして記録される。この銀粒子（黒点）の立体的な並びは、光学顕微鏡を用いた観察により抽出される。本実施形態の黒点が、本開示の飛跡像の一例である。
【００３７】
断層画像生成装置１０は、光学顕微鏡及び撮像装置（いずれも図示省略）を備える。光学顕微鏡の対物レンズは、乳剤層５１Ａのμ粒子が入射される面側に配置される。断層画像生成装置１０は、乳剤層５１の面内方向に、原子核乾板５０を移動させながら連続的に撮像を行うことで、乳剤層５１のμ粒子が入射される面全体の画像を撮像する。また、原子核乾板５０は、対物レンズを原子核乾板５０における、乳剤層５１及びフィルムベース５２の積層方向（図１では縦方向、以下、「深さ方向」という）に移動させることで、対物レンズの焦点位置を変化させて、乳剤層５１の深さが異なる複数の画像を撮像する。
【００３８】
このようにして断層画像生成装置１０により、図３に示すように、乳剤層５１Ａからは、深さ方向が異なり、複数の画素５８（図７参照）を各々含む、ｎ（ｎは２以上の整数）枚の断層画像５７１〜５７ｎを含む断層画像群５６Ａが生成される。また、乳剤層５１Ｂからは、深さ方向が異なる、ｎ枚の断層画像５７１〜５７ｎを含む断層画像群５６Ｂが生成される。断層画像５７１及び断層画像５７ｎのうち断層画像５７１側が、飛跡Ｔの入射側に最も近い断面であり、一例として、断層画像群５６の最上層に対応する。また、断層画像５７１及び断層画像５７ｎのうち断層画像５７ｎ側が、飛跡Ｔの入射側に最も遠い断面であり、一例として、断層画像群５６の最下層に対応する。本実施形態の断層画像５７１が本開示の第１断層画像の一例に対応し、本実施形態の断層画像５７ｎが本開示の第２断層画像の一例に対応する。なお、以下では、断層画像群５６Ａ及び断層画像群５６Ｂ各々の断層画像５７１〜５７ｎを区別せずに総称する場合、「断層画像５７」という。また、断層画像群５６Ａ及び断層画像群５６Ｂを区別せずに総称する場合、「断層画像群５６」という。断層画像生成装置１０により生成された複数の断層画像５７の画像データは、飛跡検出装置１２に送信される。
【００３９】
飛跡検出装置１２は、断層画像生成装置１０により生成された断層画像５７の画像データを取得し、取得した断層画像５７に基づいて、原子核乾板５０の乳剤層５１における、飛跡Ｔを検出する装置である。飛跡検出装置１２による飛跡Ｔの検出結果は、解析装置１４に送信される。
【００４０】
解析装置１４は、飛跡検出装置１２による飛跡Ｔの検出結果を取得し、飛跡Ｔを解析することによって検査対象物の密度等を導出する装置である。
【００４１】
次に、図４を参照して、本実施形態の飛跡検出装置１２のハードウェア構成を説明する。図４に示すように、飛跡検出装置１２は、ＣＰＵ（Central Processing Unit）３０、一時記憶領域としてのメモリ３１、及び不揮発性の記憶部３２を含む。また、飛跡検出装置１２は、液晶ディスプレイ等の表示部３３、キーボードとマウス等の入力部３４、及びネットワークＮに接続されるネットワークＩ／Ｆ（InterFace）３６を含む。ＣＰＵ３０、メモリ３１、記憶部３２、表示部３３、入力部３４、及びネットワークＩ／Ｆ３６は、バス３９に接続される。解析装置１４の例としては、サーバコンピュータが挙げられる。
【００４２】
記憶部３２は、ＨＤＤ（Hard Disk Drive）、ＳＳＤ（Solid State Drive）、及びフラッシュメモリ等によって実現される。記憶媒体としての記憶部３２には、飛跡検出プログラム３７が記憶される。ＣＰＵ３０は、記憶部３２から飛跡検出プログラム３７を読み出してからメモリ３１に展開し、展開した飛跡検出プログラム３７を実行する。また、記憶部３２には、詳細を後述する飛跡角度データ３８が記憶される。
【００４３】
上述のように、ＣＰＵ３０が飛跡検出プログラム３７を実行することにより、図５に示す、断層画像取得部２０、前処理部２２、飛跡候補線特定部２４、検出部２６、及び出力部２８の各々として機能する。図５には、本実施形態の飛跡検出装置１２の機能的な構成の一例を表すブロック図が示されている。
【００４４】
断層画像取得部２０は、断層画像生成装置１０から送信された断層画像５７の画像データを取得し、取得した断層画像５７の画像データを前処理部２２に出力する。
【００４５】
前処理部２２は、入力された断層画像５７の画像データに対して、前処理となる画像処理を行う。前処理とは、後段の飛跡候補線特定部２４において断層画像５７から飛跡Ｔを検出し易くするための処理であり、断層画像５７から飛跡Ｔを表す黒点を検出し易くするための処理である。前処理としては、例えば、オフセット補正、ゲイン補正、及び画像ムラの等の各種補正が挙げられる。前処理後の断層画像５７の画像データは、飛跡候補線特定部２４に出力される。
【００４６】
飛跡候補線特定部２４は、入力された画像データが表す断層画像５７に基づいて、想定される飛跡Ｔの飛跡角度に応じて、断層画像５７を通る飛跡Ｔを表す飛跡線の候補となる、飛跡線候補を特定し、特定結果を表す情報を検出部２６に出力する。なお、「飛跡角度」とは、乳剤層５１Ａにμ粒子が入射する入射角度であり、本実施形態では、図６に示すように、天頂角θ、及び乳剤層５１（断層画像５７）の面に平行な面内角ωの組み合わせで表される角度であり、μ粒子が入射する方向ともいう。なお、図６に示すように、面内角ωは、面内方向のうち予め定められた方向をω＝０度とし、天頂角θ＝０度の軸に対して、半時計回りの方向を正とした角度で表され、ω＝０度以上、ω＝３６０度未満で定義される。
【００４７】
飛跡Ｔの飛跡角度（θ，ω）が（０，０）以外の場合、換言すると、飛跡Ｔが断層画像５７に斜めに入射した場合、一例として図７に示すように、断層画像群５６の各断層画像５７を通る飛跡Ｔの各断層画像５７内における位置が異なる。一例として図７では、断層画像群５６が４枚の断層画像５７（ｎ＝４）を含む場合を示している。図７に示す例では、断層画像５７１の画素５８（ｘ１，ｙ１）の基準点６０を通る飛跡Ｔは、断層画像５７２、断層画像５７３、及び断層画像５７２の各々における基準点６０に対応する位置と異なる位置であり、かつ各断層画像５７で異なる位置を通る。なお、本実施形態では、一例として基準点６０を、画素５８の中心点としている。
【００４８】
想定される飛跡Ｔの飛跡角度として、天頂角θ、及び面内角ωを網羅的に定めた場合の仮想的な飛跡Ｔによる飛跡線を、仮想飛跡候補線とする。この場合、図８に示す例のように、最上層である断層画像５７１の画素５８（ｘ１，ｙ１）の基準点６０を通り、最下層である断層画像５７ｎの同一の画素５８（ｘｎ，ｙｎ）を通る、複数（図８では５本）の仮想飛跡候補線ＶＣが存在する。
【００４９】
後段の検出部２６では、断層画像５７に形成された飛跡Ｔに応じた黒点に基づいて、断層画像５７の画素単位で演算を行い、飛跡候補線特定部２４で特定した飛跡候補線から、飛跡Ｔによる飛跡線を検出する。そのため、最下層である断層画像５７ｎの同一の画素５８を通る、複数の仮想飛跡候補線ＶＣの各々を用いて演算を行うと、同一の黒点に対応する演算を行う、換言すると同等の演算を繰り返し行うことになる。そこで、本実施形態の飛跡候補線特定部２４は、図９に示すように、最上層である断層画像５７１の画素５８（ｘ１，ｙ１）の基準点６０を通る複数の仮想飛跡候補線ＶＣのうちから、最下層である断層画像５７ｎの同一の画素５８を通る１本の飛跡候補線Ｃを特定することで、後段の検出部２６における、同等の演算を繰り返すことによる演算量の増加を抑制する。
【００５０】
一例として、本実施形態の飛跡候補線特定部２４では、断層画像５７ｎを通る複数の仮想飛跡候補線ＶＣのうち、各画素５８の中心点６１に近い１本の仮想飛跡候補線ＶＣを、飛跡候補線Ｃとして特定する。図９に示した例では、飛跡候補線特定部２４が、断層画像５７ｎの９つの画素５８の各々を通る飛跡候補線Ｃを、仮想飛跡候補線ＶＣから特定した例を示している。
【００５１】
従って、本実施形態では、飛跡候補線の飛跡角度（以下、「飛跡候補角度」という）（θ，ω）は、断層画像５７における画素５８の大きさと、断層画像５７１から断層画像５７ｎまでの距離と、天頂角θ及び面内角ωの刻みと、に応じて定まる。本実施形態の飛跡候補線特定部２４は、飛跡角度（θ，ω）を、断層画像５７における画素５８の大きさと、断層画像５７１から断層画像５７ｎまでの距離と、天頂角θ及び面内角ωの刻みとの３つの条件に対応付けた情報を飛跡角度データ３８として記憶部３２に記憶している。
【００５２】
検出部２６は、飛跡候補線特定部２４が特定した複数の飛跡候補線から、飛跡Ｔを検出する。一例として、本実施形態の検出部２６は、いわゆる、シフト法を用いて飛跡Ｔを検出する。具体的には検出部２６は、飛跡Ｔの飛跡角度（θ，ω）に応じたシフト量で各断層画像５７をずらし、該当する飛跡角度（θ，ω）を有する飛跡Ｔが、各断層画像５７上のほぼ同じｘ、ｙ座標に位置する状態として、各断層画像５７の同一座標における黒点の有無に基づいて、飛跡Ｔを検出し、検出結果を表す情報を出力部２８に出力する。
【００５３】
出力部２８は、検出部２６の検出結果を表す情報、すなわち、飛跡Ｔを表す情報を、ネットワークＩ／Ｆ３６により、ネットワークＮを介して解析装置１４へ出力する。
【００５４】
次に、本実施形態の飛跡検出装置１２の作用を説明する。まず、ＣＰＵ３０が飛跡検出プログラム３７を実行することによって、図１０に示す飛跡検出処理が実行される。図１０に示す飛跡検出処理は、例えば、飛跡検出装置１２のユーザによる飛跡検出の指示が入力部３４により入力されたタイミング等で実行される。
【００５５】
図１０のステップＳ１００で断層画像取得部２０は、上述したように、断層画像生成装置１０から断層画像５７の画像データを取得する。断層画像取得部２０により取得された断層画像５７の画像データは、上述したように、前処理部２２により前処理が行われる。
【００５６】
次のステップＳ１０２で飛跡候補線特定部２４は、詳細を後述する、図１１に一例を示した飛跡候補線特定処理を行い、上述したように、前処理が行われた断層画像５７に基づいて、飛跡Ｔを表す飛跡線の候補となる、飛跡候補線を特定する。
【００５７】
次のステップＳ１０４で検出部２６は、詳細を後述する、図１５に一例を示した検出処理を行い、上述したように、上記ステップＳ１０２で特定された飛跡候補線Ｃから飛跡Ｔを検出した後、本飛跡検出処理を終了する。
【００５８】
上述した飛跡検出処理（図１０参照）のステップＳ１０２で実行される、飛跡候補線特定処理について図１１を参照して詳細に説明する。
【００５９】
図１１のステップＳ１２０で飛跡候補線特定部２４は、飛跡候補角度（θ，ω）を導出するか否か判定する。本実施形態では、飛跡候補線Ｃの特定に、記憶部３２に記憶されている飛跡角度データ３８による飛跡候補角度（θ，ω）を用いる場合と、新たに導出した飛跡候補角度（θ，ω）を用いる場合とがある。一例として本実施形態の飛跡候補線特定部２４では、ユーザから入力部３４により、飛跡候補角度（θ，ω）の導出を指示された場合、及び飛跡角度データ３８に記憶されている情報と、上記の条件が一致しない場合、飛跡候補角度（θ，ω）を導出すると判定する。上記の条件が一致しない場合の例としては、飛跡角度データ３８において、飛跡候補角度（θ，ω）が対応付けられている画素５８の大きさと、上記ステップＳ１００で取得した断層画像５７の画素５８の大きさとが異なり、その画素５８の大きさの差が許容範囲を越えている場合が挙げられる。
【００６０】
飛跡候補角度（θ，ω）を導出する場合、ステップＳ１２０の判定が肯定判定となり、ステップＳ１２２へ移行する。ステップＳ１２２で飛跡候補線特定部２４は、詳細を後述する、図１２に一例を示した飛跡候補角度導出処理を実行した後、本飛跡候補線特定処理を終了し、図１０に示した飛跡検出処理のステップＳ１０４へ移行する。
【００６１】
一方、飛跡候補角度（θ，ω）を導出しない場合、換言すると、飛跡角度データ３８を用いる場合、飛跡候補線特定処理のステップＳ１２０の判定が否定判定となり、ステップＳ１２４へ移行する。
【００６２】
ステップＳ１２４で飛跡候補線特定部２４は、飛跡角度データ３８から、飛跡候補角度（θ，ω）を特定する。具体的には、飛跡候補線特定部２４は、断層画像取得部２０が取得した断層画像５７に応じた、断層画像５７における画素５８の大きさと、断層画像５７１から断層画像５７ｎまでの距離と、天頂角θ及び面内角ωの刻みとに対応付けられている飛跡候補角度（θ，ω）を飛跡角度データ３８から取得する。ステップＳ１２４を終了することにより、本飛跡候補線特定処理が終了し、図１０に示した飛跡検出処理のステップＳ１０４へ移行する。
【００６３】
上述した飛跡候補線特定処理（図１１参照）のステップＳ１２０で実行される、飛跡候補角度導出処理について図１２を参照して詳細に説明する。
【００６４】
図１２のステップＳ２００で飛跡候補線特定部２４は、飛跡候補角度（θ，ω）として想定される、天頂角θ、及び面内角ωの全ての組み合わせに応じた一覧表を生成する。換言すると、全ての仮想飛跡候補線ＶＣの飛跡角度（以下、「仮想飛跡候補角度」という）（θ，ω）の天頂角θ、及び面内角ωに応じた一覧表を生成する。図１３には、一覧表７０の一例を示す。ステップＳ２００では、図１３に示した一覧表７０における、「線ＩＤ」、「θ」、及び「ω」の欄に、データが入力された状態の一覧表７０が生成される。
【００６５】
なお、一般的にμ粒子は、天頂から降り注ぐため、原子核乾板５０を、法線方向が天頂角θ＝０度となる状態に配置した場合、原子核乾板５０に到達するμ粒子の量は、天頂角θ＝９０度未満の範囲内が多く、天頂角θ＝４５度以下の範囲内が特に多い。そこで、本実施形態の飛跡候補線特定部２４は、天頂角θが０度以上、４５度以下、かつ面内角ωが０度以上、３６０度未満の範囲で１度ごとに網羅的に仮想飛跡候補角度（θ，ω）を定めている。なお、天頂角θが０度の場合、面内角ωの角度にかかわらず、得られる飛跡候補線Ｃが、同一となる。そのため、本実施形態の飛跡候補線特定部２４では、天頂角θが０度の場合、面内角ωが０度の場合のみを仮想飛跡候補角度（θ，ω）に含める。従って、仮想飛跡候補角度（θ，ω）は、１６２０１（４５×３６０＋１＝１６２０１）通りとなり、組み合わせ毎に、図１３に示すように、線ＩＤ１〜ＩＤ１６２０１が順次、対応付けされている。
【００６６】
次のステップＳ２０２で飛跡候補線特定部２４は、最上段の断層画像５７１の基準点６０を通り、仮想飛跡候補角度（θ，ω）を有する全ての仮想飛跡候補線ＶＣが通る、最下段の断層画像５７ｎの画素５８の、断層画像５７ｎにおける位置を導出し、一覧表７０の「画素座標ｘ」及び「画素座標ｙ」の欄の各々に入力する。なお、本実施形態では、断層画像５７における画素５８の位置として、画素単位で表される座標を用いている。図１３に示した一覧表７０では、線ＩＤが「３」で仮想飛跡候補角度（１，１）である仮想飛跡候補線ＶＣは、最下段の断層画像５７ｎにおける画素座標（５，０）の位置にある画素５８を通ることを示している。
【００６７】
次のステップＳ２０４で飛跡候補線特定部２４は、全ての仮想飛跡候補線ＶＣについて、仮想飛跡候補線ＶＣが通る、最下段の断層画像５７ｎの画素５８における中心点６１と、通過位置との距離を導出し、一覧表７０の「Δｄ」の欄に入力する。Δｄは、仮想飛跡候補角度（θ，ω）、断層画像５７における画素５８の大きさ、及び断層画像５７１から断層画像５７ｎまでの距離から導出される。
【００６８】
次のステップＳ２０６で飛跡候補線特定部２４は、図１４に一例を示したＩＤ配列表７２の各値を「０」にリセットする。図１４に一例を示したＩＤ配列表７２は、最下層の断層画像５７ｎにおける画素５８の座標（位置）を項目に採用した表である。本実施形態のＩＤ配列表７２は、画素５８を通る仮想飛跡候補線ＶＣが存在しない場合、「０」が入力されている。また、詳細は後述するが、ＩＤ配列表７２は、画素５８を通る１本の仮想飛跡候補線ＶＣが存在する場合、その画素５８を通る仮想飛跡候補線ＶＣの線ＩＤが入力される。さらに、ＩＤ配列表７２には、画素５８を通る複数の仮想飛跡候補線ＶＣがある場合、その画素５８を通る複数の仮想飛跡候補線ＶＣのうち、中心点６１に最も近い位置を通る１本の仮想飛跡候補線ＶＣの線ＩＤが入力される。
【００６９】
次のステップＳ２０８で飛跡候補線特定部２４は、注目ＩＤを１（注目ＩＤ＝１）とする。
【００７０】
次のステップＳ２１０で飛跡候補線特定部２４は、一覧表７０から、注目ＩＤに対応する画素座標（ｘ，ｙ）を取得する。図１３に示した例では、注目ＩＤが「１」の場合、一覧表７０から画素座標（０，０）を取得し、注目ＩＤが「３」の場合、画素座標（５，０）を一覧表７０から取得する。
【００７１】
次のステップＳ２１２で飛跡候補線特定部２４は、注目ＩＤに対応するＩＤ配列表７２の画素座標（ｘ，ｙ）の値が０（（ｘ，ｙ）＝０）であるか否かを判定する。画素座標（ｘ，ｙ）の値が０の場合、ステップＳ２１２の判定が肯定判定となり、ステップＳ２１４へ移行する。すなわち、注目ＩＤに対応する仮想飛跡候補線ＶＣの他に、上記ステップＳ２１０で取得した画素座標（ｘ，ｙ）の画素５８を通る仮想飛跡候補線ＶＣが存在しない場合、ステップＳ２１２の判定が肯定判定となり、ステップＳ２１４へ移行する。
【００７２】
ステップＳ２１４で飛跡候補線特定部２４は、一覧表７０の注目ＩＤに対応する「有効／無効」の欄に、有効を表す「１」を入力した後、ステップＳ２２８へ移行する。また、飛跡候補線特定部２４は、上記ステップＳ２１０で取得した画素座標（ｘ，ｙ）に対応する、ＩＤ配列表７２の（ｘ，ｙ）の値に、注目ＩＤの番号を入力する。なお、一覧表７０に表される「有効／無効」とは、その仮想飛跡候補線ＶＣが、飛跡候補線Ｃとして有効であるか無効であるかを表しており、「有効」の場合、上述のように「１」が入力され、無効の場合「０」が入力される。
【００７３】
一方、画素座標（ｘ、ｙ）の値が０ではない場合、ステップＳ２１２の判定が否定判定となり、ステップＳ２１６へ移行する。すなわち、注目ＩＤに対応する仮想飛跡候補線ＶＣの他に、上記ステップＳ２１０で取得した画素座標（ｘ，ｙ）の画素５８を通る仮想飛跡候補線ＶＣが存在する場合、ステップＳ２１２の判定が否定判定となり、ステップＳ２１６へ移行する。
【００７４】
ステップＳ２１６で飛跡候補線特定部２４は、一覧表７０の注目ＩＤに対応するΔｄの値と、ＩＤ配列表７２における注目ＩＤの画素座標（ｘ，ｙ）に入力されている線ＩＤを有する仮想飛跡候補線ＶＣのΔｄとを比較する。次のステップＳ２１８で飛跡候補線特定部２４は、注目ＩＤのΔｄが、上記ステップＳ２１６で比較した線ＩＤのΔｄよりも小さいか否かを判定する。
【００７５】
注目ＩＤのΔｄが、上記ステップＳ２１６で比較した線ＩＤのΔｄよりも小さい場合、ステップＳ２１８の判定が肯定判定となり、ステップＳ２２２へ移行する。このようにステップＳ２１８が肯定判定となる場合とは、注目ＩＤに対応する仮想飛跡候補線ＶＣの方が、上記ステップＳ２１６で比較した線ＩＤに対応する仮想飛跡候補線ＶＣよりも、画素５８の中心点６１に近い位置を通る場合である。
【００７６】
ステップＳ２２２で飛跡候補線特定部２４は、一覧表７０の注目ＩＤに対応する「有効／無効」の欄に、有効を表す「１」を入力する。次のステップＳ２２４で飛跡候補線特定部２４は、一覧表７０の、上記ステップＳ２１６で比較した線ＩＤに対応する「有効／無効」の欄に入力されている値を、有効を表す「１」から、無効を表す「０」にする。
【００７７】
次のステップＳ２２６で飛跡候補線特定部２４は、注目ＩＤの画素座標（ｘ，ｙ）に対応する、ＩＤ配列表７２の（ｘ，ｙ）の値に、注目ＩＤの番号を入力した後、ステップＳ２２８へ移行する。
【００７８】
一方、注目ＩＤのΔｄが、上記ステップＳ２１６で比較した線ＩＤのΔｄよりも小さくない場合、換言すると、注目ＩＤのΔｄが、上記ステップＳ２１６で比較した線ＩＤのΔｄ以上の場合、ステップＳ２１８の判定が否定判定となり、ステップＳ２２０へ移行する。このようにステップＳ２１８が否定判定となる場合とは、注目ＩＤに対応する仮想飛跡候補線ＶＣの方が、上記ステップＳ２１６で比較した線ＩＤに対応する仮想飛跡候補線ＶＣよりも、画素５８の中心点６１から遠い位置を通る場合、または注目ＩＤに対応する仮想飛跡候補線ＶＣと、上記ステップＳ２１６で比較した線ＩＤに対応する仮想飛跡候補線ＶＣとが同じ位置を通る場合である。
【００７９】
ステップＳ２２０で飛跡候補線特定部２４は、一覧表７０の注目ＩＤに対応する「有効／無効」の欄に、無効を表す「０」を入力した後、ステップＳ２２８へ移行する。
【００８０】
例えば、注目ＩＤが「３」の場合、一覧表７０によれば、注目ＩＤの画素座標（ｘ，ｙ）は、「（５，０）」であり、Δｄは、「０．０８０」である。ＩＤ配列表７２によれば、画素座標（５，０）に入力されている線ＩＤは、「２」である。一覧表７０によれば、線ＩＤが「２」に対応する、Δｄの値は「０．０７８」であり、上記注目ＩＤに対応する仮想飛跡候補線ＶＣよりも中心点６１に近い距離を通るため、ステップＳ２１８の判定が否定判定となる。次のステップＳ２２０で飛跡候補線特定部２４は、注目ＩＤ＝線ＩＤが「３」に対応する「有効／無効」の欄に、無効を表す「０」を入力する。
【００８１】
なお、本実施形態では、上記ステップＳ２１８において、注目ＩＤのΔｄと、上記ステップＳ２１６で比較した線ＩＤのΔｄとが等しい場合、否定判定となりステップＳ２２０へ移行することで、先に入力されている線ＩＤに対応する仮想飛跡候補線ＶＣを飛跡候補線Ｃとして採用する形態について説明した。しかしながら、注目ＩＤのΔｄと、上記ステップＳ２１６で比較した線ＩＤのΔｄとが等しい場合、いずれの仮想飛跡候補線ＶＣを、飛跡候補線Ｃとして採用するかは、本実施形態に限定されない。例えば、天頂角θの大小、または面内角ωの大小に応じて何れを採用するかを定めておいてもよい。
【００８２】
ステップＳ２２８で飛跡候補線特定部２４は、注目ＩＤが線ＩＤの最大値未満（注目ＩＤ＜最大値）か否かを判定する。上述したように、本実施形態では、線ＩＤの最大値が「１６２０１」であるため、注目ＩＤの値が「１６２０１」未満の場合、ステップＳ２２８の判定が肯定判定となり、ステップＳ２３０へ移行する。ステップＳ２３０で飛跡候補線特定部２４は、注目ＩＤに１を加算した値を新たな注目ＩＤとした（注目ＩＤ＝注目ＩＤ＋１）後、ステップＳ２１０に戻り、上記ステップＳ２１０〜Ｓ２２８の処理を繰り返す。
【００８３】
一方、注目ＩＤが線ＩＤの最大値未満ではない場合、換言すると、注目ＩＤが線ＩＤの最大値に達した場合、ステップＳ２２８の判定が否定判定となり、ステップＳ２３２へ移行する。ステップＳ２３２で飛跡候補線特定部２４は、一覧表７０の、「有効／無効」の欄の値が有効を表す「１」である、線ＩＤに対応する仮想飛跡候補角度（θ，ω）を、飛跡候補線Ｃの飛跡候補角度（θ，ω）として特定した後、本飛跡候補角度導出処理を終了する。本飛跡候補角度導出処理を終了することにより、図１１に示した飛跡候補線特定処理も終了し、図１０に示した飛跡検出処理のステップＳ１０４へ移行する。
【００８４】
なお、上記飛跡候補角度導出処理により導出された飛跡候補角度（θ，ω）のデータを、断層画像５７における画素５８の大きさと、断層画像５７１から断層画像５７ｎまでの距離と、天頂角θ及び面内角ωの刻みとの３つの条件に対応付けて飛跡角度データ３８として記憶部３２に記憶させてもよい。このように飛跡角度データ３８を記憶させることにより、次回以降、飛跡検出処理を行う場合、上記３つの条件を満たせば、飛跡候補線特定部２４は、飛跡角度データ３８を利用して飛跡候補線Ｃを特定することができるため、飛跡候補線Ｃを特性するための演算量を削減することができる。
【００８５】
上述した飛跡検出処理（図１０参照）のステップＳ１０４で実行される、検出処理について図１５を参照して詳細に説明する。
【００８６】
図１５のステップＳ１４０で検出部２６は、上記ステップＳ１０２で特定された複数の飛跡候補線Ｃのうちから、ある飛跡候補角度（θ，ω）に応じたシフト量で、各断層画像５７をシフトさせる。図１６Ａに示す、飛跡候補角度（θ，ω）を有する飛跡候補線Ｃが、５枚の断層画像５７（ｎ＝５）を含む断層画像群５６を通る場合を具体例として説明する。この場合、検出部２６は、図１６Ｂに示すように、各断層画像５７を飛跡候補角度（θ，ω）に応じたシフト量でずらして各画素５８の座標を変換することより、飛跡候補線Ｃを、断層画像５７に対して垂直、換言すると飛跡候補線Ｃを断層画像５７における法線に一致させた状態とする。
【００８７】
次のステップＳ１４２で検出部２６は、各断層画像５７の、ある座標における画素５８の画素値を加算する。本実施形態では、図１６Ｂに示すように、各断層画像５７を、飛跡候補角度（θ，ω）に応じたシフト量でずらした状態において、全ての断層画像５７が重なる領域を飛跡Ｔの検出対象とし、一部の断層画像５７が重ならない各断層画像５７の端部の領域は、飛跡Ｔの検出対象外とする。なお、このように飛跡Ｔの検出において除外される各断層画像５７の端部の領域は、断層画像５７１から断層画像５７ｎまでの厚みに対し、合成断層画像５９１から断層画像５７ｎの大きさが比較的大きいため、相対的に面積が小さくなるため、端部の領域を検出対象から除外することによる影響は抑制される。
【００８８】
そのため、検出部２６は、飛跡Ｔの検出対象内のある座標の画素５８について、各断層画像５７の画素値を加算する。換言すると、検出部２６は、飛跡候補線Ｃが通る各断層画像５７の画素５８の画素値を加算する。
【００８９】
次のステップＳ１４４で検出部２６は、上記ステップＳ１４２で得られた加算値が、予め定められた閾値以下（加算値≦閾値）であるか否かを判定する。本実施形態では、画素値が小さいほど、画像が暗く、画素値が大きいほど、画像が明るくなる。換言すると、画素値が小さいほど画像が黒く、黒点となり、画素値が大きいほど画像が白く、白点となる。従って、飛跡Ｔにより断層画像５７に黒点が形成されている場合、黒点となった画素５８の画素値は小さくなる。本実施形態では、全ての断層画像５７が黒点であった場合の画素の加算値に基づいて、飛跡Ｔにより断層画像群５６に黒点が形成されていることを判別するための閾値を予め定めておく。全ての断層画像５７が黒点であった場合の画素の加算値と閾値の関係は、使用する乳剤の特性により変わる。これは、μ粒子が、通過した乳剤をどれだけの割合で感光させるか、つまりは何個の黒点が発生するかが、乳剤の感度に依存するためである。さらに、天頂角θが大きくなるとμ粒子の通過距離が長くなるので発生する黒点の数が増える。実験では、天頂角θが０度の場合に、全ての断層画像５７が黒点であった場合の画素の加算値に対して、８０％の断層画像５７が黒点であった場合の画素の加算値を閾値にすることで良好な結果を得た。
【００９０】
加算値が閾値以下ではない場合、換言すると、加算値が閾値を超える場合、ステップＳ１４４の判定が否定判定となり、ステップＳ１４８へ移行する。一方、加算値が閾値以下の場合、ステップＳ１４４の判定が肯定判定となり、ステップＳ１４６へ移行する。
【００９１】
ステップＳ１４６で検出部２６は、上記ステップＳ１４０における飛跡候補角度（θ，ω）を有し、上記ステップＳ１４２における座標を有する飛跡候補線Ｃを飛跡Ｔとして特定する。
【００９２】
次のステップＳ１４８で検出部２６は、飛跡Ｔの検出対象内の全ての座標（画素５８）について、上記ステップＳ１４２〜Ｓ１４６の処理を行ったか否かを判定する。未だ、上記ステップＳ１４２〜Ｓ１４６の処理を行っていない座標（画素５８）が存在する場合、ステップＳ１４８の判定が否定判定となり、ステップＳ１５０へ移行する。ステップＳ１５０で検出部２６が、飛跡Ｔの検出対象となる座標（画素５８）を変更した後、ステップＳ１４２に戻り、上記ステップＳ１４２〜Ｓ１４６の処理を繰り返す。
【００９３】
一方、飛跡Ｔの検出対象となる全ての座標（画素５８）について、上記ステップＳ１４２〜Ｓ１４６の処理を行った場合、ステップＳ１４８の判定が肯定判定となり、ステップＳ１５２へ移行する。
【００９４】
ステップＳ１５２で検出部２６は、上記ステップＳ１０２で特定された複数の飛跡候補線Ｃの全てについて上記ステップＳ１４０〜Ｓ１５０の処理を行ったか否かを判定する。未だ、上記ステップＳ１４０〜Ｓ１５０の処理を行っていない、飛跡候補線Ｃがある場合、ステップＳ１５２の判定が否定判定となり、ステップＳ１５４へ移行する。次のステップＳ１５４で検出部２６が、飛跡Ｔの検出対象となる飛跡候補線Ｃを変更した後、ステップＳ１４０に戻り、変更した飛跡候補線Ｃの飛跡候補角度（θ，ω）に基づき、上記ステップＳ１４０〜Ｓ１５２の処理を繰り返す。
【００９５】
一方、上記ステップＳ１０２で特定された複数の飛跡候補線Ｃの全てについて上記ステップＳ１４０〜Ｓ１５０の処理を行った場合、ステップＳ１５２の判定が肯定判定となり、ステップＳ１５６へ移行する。
【００９６】
ステップＳ１５６で検出部２６は、上記ステップＳ１４６で飛跡Ｔとして特定した全ての飛跡Ｔを表す情報を、検出結果として出力部２８に出力した後、本検出処理を終了する。本検出処理が終了することにより、飛跡検出処理（図１０参照）のステップＳ１０４が終了し、本飛跡検出処理が終了する。本願では図示していないが、この処理の後に、重複飛跡のまとめ処理を行う。黒点は、合成断層画像５９や断層画像５７上で、複数の画素５８を占める大きさを持っている場合、原子核乾板５０の一つの飛跡Ｔが、複数の飛跡Ｔとして検出される場合がある。ただし、これらの検出された複数の飛跡Ｔは、隣り合った画素５８など極めて近い位置を通過しており、かつ、飛跡Ｔの方向がほぼ同じであるという特徴をもっているので、重複して検出された飛跡Ｔとして抽出できる。重複して検出された飛跡Ｔは、それらの飛跡Ｔの平均的な通過位置、飛翔方向を持った一つの飛跡Ｔとしてまとめる処理を行っている。
【００９７】
なお、飛跡検出装置１２の検出部２６が実行する飛跡検出処理（図１５参照）における飛跡候補線Ｃが飛跡Ｔであるか否かの判定方法について、飛跡候補線Ｃが通る各断層画像５７の画素５８の画素値の加算値と、予め定められた閾値とを比較した比較結果に基づく方法について説明したが、上記方法に限定されない。例えば、以下の方法を採用してもよい。まず、検出部２６は、飛跡候補線Ｃが通る各断層画像５７の画素５８毎に、画素５８に黒点が形成されているか否かを判定するための閾値と比較し、黒点が形成されているか否かを判定する。さらに、検出部２６は、黒点であると判定した画素５８（断層画像５７）の数が飛跡候補線Ｃが飛跡Ｔであるか否かを判定するための閾値以上の場合に、飛跡候補線Ｃが飛跡Ｔであると判定する方法を用いてもよい。
【００９８】
このように本実施形態の飛跡検出装置１２では、飛跡候補線特定部２４が、断層画像群５６の最上層である断層画像５７１の基準点６０を通る方向が各々異なり、かつ各々が断層画像群５６の最下層である断層画像５７ｎの異なる画素５８を通る複数の飛跡候補線Ｃを特定する。本実施形態によれば、想定される天頂角θ、及び面内角ωを網羅的に定めて飛跡候補角度（θ，ω）として定めた飛跡候補線Ｃを用いて飛跡Ｔを検出する場合に比べて、飛跡Ｔの検出に用いる飛跡候補線Ｃの数を減少させることができる。
【００９９】
従って、本実施形態の飛跡検出装置１２によれば、検出精度を低下させることなく、飛跡を検出するための演算量を削減することができる。
【０１００】
［第２実施形態］
以下、第２実施形態について詳細に説明する。なお、本実施形態では、上記第１実施形態で説明した構成及び作用と同一の構成及び作用については同一の符号を付し、詳細な説明を省略する。
【０１０１】
本実施形態の非破壊検査システム１の全体の構成、及び飛跡検出装置１２の構成は、第１実施形態の非破壊検査システム１及び飛跡検出装置１２の構成（図２、図４、及び図５参照）と同様であるため、説明を省略する。
【０１０２】
本実施形態は、飛跡検出装置１２の作用の一部が異なるため、本実施形態の飛跡検出装置１２の作用について説明する。本実施形態の飛跡検出装置１２では、第１実施形態の飛跡検出処理（図１０参照）のステップＳ１０２で実行される飛跡候補線特定処理（図１１参照）のステップＳ１２２において実行される飛跡候補角度導出処理の一部が異なるため、本実施形態の飛跡検出装置１２が実行する飛跡候補角度導出処理について説明する。
【０１０３】
図１７には、本実施形態の飛跡検出装置１２が実行する飛跡候補角度導出処理の一例を表すフローチャートが示されている。本実施形態の飛跡候補角度導出処理は、ステップＳ２３２の後に、ステップＳ２３４〜Ｓ２５０の処理を行う点で第１実施形態の飛跡候補角度導出処理（図１２参照）と異なっている。
【０１０４】
図１７に示すように、ステップＳ２３２の後のステップＳ２３４で飛跡候補線特定部２４は、上記ステップＳ２３２で特定した複数の飛跡候補線Ｃのうちから、注目する飛跡候補線Ｃ（以下、「注目飛跡候補線」という）を１つ決定する。具体的には、飛跡候補線特定部２４は、一覧表７０で、「有効／無効」の欄に入力されている値が、有効を表す「１」である線ＩＤに対応する飛跡候補線Ｃから、注目飛跡候補線Ｃを決定する。
【０１０５】
次のステップＳ２３６で飛跡候補線特定部２４は、特定した複数の飛跡候補線Ｃのうちから、注目飛跡候補線Ｃの他の飛跡候補線Ｃを１つ選択する。
【０１０６】
次のステップＳ２３８で飛跡候補線特定部２４は、上記ステップＳ２３４で決定した注目飛跡候補線Ｃと、上記ステップＳ２３６で選択した飛跡候補線Ｃとの成す角度が所定値以下であるか否かを判定する。
【０１０７】
注目飛跡候補線Ｃと、選択した飛跡候補線Ｃとの成す角度が所定値を越える場合、ステップＳ２３８の判定が否定判定となり、ステップＳ２４４へ移行する。一方、注目飛跡候補線Ｃと、選択した飛跡候補線Ｃとの成す角度が所定値以下の場合、ステップＳ２３８の判定が肯定判定となり、ステップＳ２４０へ移行する。
【０１０８】
ステップＳ２４０で飛跡候補線特定部２４は、選択した飛跡候補線Ｃが、特定された飛跡候補線Ｃから除外するための除外条件を満たすか否かを判定する。飛跡候補線Ｃ同士がなす角度が所定値以下の場合、この飛跡候補線Ｃ同士の密度は、他の複数の飛跡候補線Ｃの密度に比べて高い状態にある。特に、飛跡候補線Ｃの飛跡候補角度（θ，ω）における天頂角θが比較的小さい場合、図１８Ａに示すように飛跡候補角度（θ，ω）における天頂角θが比較的大きい場合に比べて、飛跡候補線Ｃの密度が高くなる傾向がある。なお、図１８Ａでは、上記ステップＳ２３２で特定された飛跡候補線Ｃの線分の長さを同一とした場合に、飛跡Ｔの入射側から断層画像５７をみた場合の飛跡候補線Ｃの端部の位置の一例が示されている。
【０１０９】
そこで、本実施形態の飛跡候補線特定部２４は、密度が高い状態の飛跡候補線Ｃを間引くことにより、特定した複数の飛跡候補線Ｃ全体の密度を、天頂角θ及び面内角ωによらず、均等な状態に近づける。図１８Ｂには、本飛跡候補角度導出処理により間引き後の飛跡候補線Ｃについて、線分の長さを同一とした場合に、飛跡Ｔの入射側から断層画像５７をみた場合の飛跡候補線Ｃの端部の位置の一例が示されている。図１８Ａ及び図１８Ｂを比較すればわかるように、間引き後の飛跡候補線Ｃ群では、飛跡候補線Ｃが規則的に設けられており、密度が均等な状態に近付いている。
【０１１０】
本実施形態では、このように密度が高い状態にある飛跡候補線Ｃを間引く条件として、除外条件を設けている。除外条件は、飛跡Ｔの検出精度を損なうことなく、かつ演算量を適切に削減できる条件であればよい。
【０１１１】
例えば、注目飛跡候補線Ｃ及び選択した飛跡候補線Ｃのうち、面内角ωが大きい方を除外するとした条件、及び飛跡候補線Ｃ及び選択した飛跡候補線Ｃのうち、面内角ωが小さい方を除外するとした条件のいずれか一方を除外条件としてもよい。この場合、飛跡候補線Ｃをより規則的に設けることができる。
【０１１２】
また例えば、天頂角θが同一の複数の飛跡候補線Ｃについて、互いの成す角度が所定値以下となる複数の飛跡候補線Ｃ同士が存在する場合、特定した複数の飛跡候補線Ｃのうちの一部を除外することにより、残存する複数の飛跡候補線Ｃの分布を、天頂角が０度の軸に対して、回転対称に最も近くなる状態とするとした、除外条件としてもよい。換言すると、互いの成す角度が所定値以下となる複数の飛跡候補線Ｃ同士が存在する場合、天頂角θが０度の軸に対して、回転対称に最も近くなる状態の飛跡候補線Ｃを残し、他の飛跡候補線Ｃを除外するとした除外条件であってもよい。この場合、面内角ωに対してより均等に間引くことができ、飛跡候補線Ｃの偏りを抑制することができる。
【０１１３】
選択した飛跡候補線Ｃが除外条件を満たさない場合、ステップＳ２４０の判定が否定判定となり、ステップＳ２４４へ移行する。一方、選択した飛跡候補線Ｃが除外条件を満たす場合、ステップＳ２４０の判定が肯定判定となり、ステップＳ２４２へ移行する。
【０１１４】
ステップＳ２４２で飛跡候補線特定部２４は、選択した飛跡候補線Ｃを、特定した飛跡候補線Ｃから除外ずる。一例として本実施形態の飛跡候補線特定部２４は、除外する飛跡候補線Ｃに対応する、一覧表７０の線ＩＤに対応付けられた「有効／無効」の欄に入力されている値を、有効を表す「１」から、無効を表す「０」にすることで、特定した飛跡候補線Ｃから除外する。
【０１１５】
次のステップＳ２４４で検出部２６は、特定した複数の飛跡候補線Ｃのうちの、注目飛跡候補線Ｃの他の飛跡候補線Ｃを、上記ステップＳ２３６において、全て選択したか否かを判定する。未だ選択していない飛跡候補線Ｃがある場合、ステップＳ２４４の判定が否定判定となり、ステップＳ２４６へ移行する。ステップＳ２４６で飛跡候補線特定部２４は、特定した複数の飛跡候補線Ｃのうちの、注目飛跡候補線Ｃの他の飛跡候補線Ｃで、未だ上記ステップＳ２３６において選択されていない飛跡候補線Ｃを１つ選択した後、ステップＳ２３８に戻り、上記ステップＳ２３８〜Ｓ２４４の処理を繰り返す。
【０１１６】
一方、特定した複数の飛跡候補線Ｃのうちの、注目飛跡候補線Ｃの他の飛跡候補線Ｃを、上記ステップＳ２３６において、全て選択した場合、ステップＳ２４４の判定が肯定判定となり、ステップＳ２４８へ移行する。
【０１１７】
ステップＳ２４８で飛跡候補線特定部２４は、特定した飛跡候補線Ｃを、全て注目飛跡候補線Ｃとしたか否かを判定する。具体的には、飛跡候補線特定部２４は、一覧表７０で、「有効／無効」の欄に入力されている値が、有効を表す「１」である線ＩＤに対応する飛跡候補線Ｃのうち、注目飛跡候補線Ｃとしていない飛跡候補線Ｃがあるか否かを判定する。未だ、注目飛跡候補線Ｃとしていない飛跡候補線Ｃがある場合、ステップＳ２４８の判定が否定判定となり、ステップＳ２５０へ移行する。ステップＳ２５０で飛跡候補線特定部２４は、注目飛跡候補線Ｃを変更した後、ステップＳ２３６に戻り、上記ステップＳ２３６〜Ｓ２４８の処理を繰り返す。
【０１１８】
一方、特定した飛跡候補線Ｃを、全て注目飛跡候補線Ｃとした場合、ステップＳ２４８の判定が肯定判定となり、本飛跡候補角度導出処理を終了する。
【０１１９】
このように本実施形態の飛跡検出装置１２では、飛跡候補線特定部２４が、さらに飛跡候補線Ｃを間引くため、飛跡Ｔの検出に用いる飛跡候補線Ｃの数をより少なくすることができるとともに、飛跡候補線Ｃ同士の偏りを抑制することができる。従って、本実施形態の飛跡検出装置１２によれば、検出精度を低下させることなく、飛跡を検出するための演算量をより削減することができる。
【０１２０】
［第３実施形態］
以下、第３実施形態について詳細に説明する。なお、本実施形態では、上記第１実施形態で説明した構成及び作用と同一の構成及び作用については同一の符号を付し、詳細な説明を省略する。
【０１２１】
本実施形態の非破壊検査システム１の全体の構成、及び飛跡検出装置１２の構成は、第１実施形態の非破壊検査システム１及び飛跡検出装置１２の構成（図２、図４、及び図５参照）と同様であるため、説明を省略する。
【０１２２】
本実施形態は、飛跡検出装置１２の作用の一部が異なるため、本実施形態の飛跡検出装置１２の作用について説明する。本実施形態の飛跡検出装置１２では、第１実施形態の飛跡検出処理（図１０参照）のステップＳ１０４で実行される検出処理（図１５参照）の一部が異なるため、本実施形態の飛跡検出装置１２が実行する検出処理について説明する。
【０１２３】
上述したように、飛跡Ｔの飛跡角度（θ，ω）が（０，０）以外の場合、換言すると、飛跡Ｔが断層画像５７に斜めに入射した場合、一例として図１９に示すように、断層画像群５６の各断層画像５７を通る飛跡Ｔの各断層画像５７内における位置が異なる。一例として図１９では、断層画像群５６が５枚の断層画像５７（ｎ＝５）を含み、断層画像５７１の基準点６３を飛跡Ｔが通る場合を示している。図１９では、以降の説明のために、飛跡Ｔが断層画像５７3と断層画像５７5を通る位置が、断層画像５７3と断層画像５７5の中心になるように図示しているので、断層画像群５６の断層画像５７１と断層画像５７２、断層画像５７3、断層画像５７4、断層画像５７5の互いの間隔が、一見等しくないように見えるが、これらは、図１６Ａと同様に等間隔である。
【０１２４】
図２０には、図１９の断層画像群５６における、断層画像５７１及び断層画像５７２を矢印Ａの方向からみた場合の、基準点６３を中心にした９つの画素５８の一例が示されている。図２０には、飛跡Ｔの飛跡角度（θ，ω）において、天頂角θを最大値とし、面内角ωを０度〜３６０度とした場合の、飛跡Ｔが断層画像５７２を通る位置が点線６４で示されている。図２０に示すように、従って、断層画像５７１の基準点６３を通る飛跡Ｔが、断層画像５７２を通る位置は、飛跡角度（θ，ω）によらず、点線６４で示される円の内側となる。これと同様のことが、図１９の断層画像群５６における、断層画像５７３及び断層画像５７４についても成り立つ。
【０１２５】
このように、隣接する断層画像５７同士の間隔が比較的短い場合、隣接する断層画像５７同士において、飛跡Ｔが通過する位置のずれ量は、±１画素の範囲内となり、断層画像５７全体としてみれば、非常に小さなずれ量となる。
【０１２６】
そこで、本実施形態の検出部２６は、隣接する２つの断層画像５７を組み合わせて得られた合成断層画像に含まれる黒点、より具体的には、合成断層画像の画素５８の画素値に基づいて飛跡Ｔを検出する。
【０１２７】
図２１には、本実施形態の飛跡検出装置１２が実行する検出処理の一例を表すフローチャートが示されている。本実施形態の検出処理は、ステップＳ１３９が設けられている点、及びステップＳ１４０及びＳ１４２に代わり、ステップＳ１４１及びＳ１４３が設けられている点で第１実施形態の検出処理（図１５参照）と異なっている。
【０１２８】
図２１に示すように、本実施形態では、検出処理がスタートすると、ステップＳ１３９で検出部２６は、隣接する２つの断層画像５７を組み合わせて合成した、上述の合成断層画像を生成する。上述したように、本実施形態では、隣接する２つの断層画像５７同士における、飛跡Ｔが通過する位置のずれ量は、±１画素の範囲内である。そのため、本実施形態では、図２２に示すように、隣接する２つの断層画像５７のうち、上側の断層画像５７ｍ（ｍは１以上、ｎ以下の奇数）に対して下側の断層画像５７ｍ＋１をずらすことなくそのまま重畳した合成断層画像５９Ｅ、及び±１画素、断層画像５７の面内におけるｘ方向及びｙ方向にずらして重畳した合成断層画像５９Ａ〜５９Ｄ、５９Ｆ〜５９Ｉを生成する。以下では、合成断層画像５９Ａ〜５９Ｉの個々を区別せずに総称する場合、単に合成断層画像５９という。なお、検出部２６は、断層画像５７ｍ及び断層画像５７ｍ＋１が重畳する領域の画像のみを合成断層画像５９としており、「重畳」とは、重畳する断層画像５７同士の画素５８の画素値を加算することをいう。
【０１２９】
このように、本実施形態では、隣接する２つの断層画像５７から、９つの合成断層画像５９が生成される。
【０１３０】
なお、断層画像群５６に含まれる断層画像５７の数（ｎ）が奇数の場合、最下層の断層画像５７ｎは、他の断層画像５７と合成せずに、そのまま用いる。
【０１３１】
次のステップＳ１４１で検出部２６は、上記ステップＳ１０２で特定された複数の飛跡候補線Ｃのうちから、ある飛跡候補角度（θ，ω）に応じたシフト量で、各合成断層画像５９、または合成断層画像５９及び最下層の断層画像５７ｎをシフトさせる。すなわち、検出部２６は、第１実施形態の検出処理（図１５参照）のステップＳ１４０の断層画像５７を、合成断層画像５９に置き換えた処理を実行する。なお、合成断層画像５９Ａ〜５９Ｉのうち、いずれの合成断層画像５９を用いるかは、飛跡候補角度（θ，ω）に応じて定められる。例えば、飛跡候補角度（θ，ω）における天頂角θが比較的小さい場合、合成断層画像５９Ｅが用いられる。
【０１３２】
具体的には、検出部２６は、図２３Ａに示す、飛跡候補角度（θ，ω）を有する飛跡候補線Ｃが、２枚の合成断層画像５９（５９１及び５９２）と断層画像５７５（ｎ＝５）を通る場合を具体例として説明する。合成断層画像５９１は断層画像５７１と断層画像５７２を重畳した断層画像であり、合成断層画像５９２は断層画像５７３と断層画像５７４を重畳した断層画像である。この場合、検出部２６は、図２３Ｂに示すように、合成断層画像５９及び断層画像５７ｎの各々を飛跡候補角度（θ，ω）に応じたシフト量でずらして各画素５８の座標を変換することより、飛跡候補線Ｃを、合成断層画像５９及び断層画像５７ｎに対して垂直、換言すると飛跡候補線Ｃを合成断層画像５９及び断層画像５７ｎにおける法線に一致させた状態とする。
【０１３３】
次のステップＳ１４３で検出部２６は、各合成断層画像５９、または合成断層画像５９及び断層画像５７ｎの各々の、ある座標における画素５８の画素値を加算する。すなわち、検出部２６は、第１実施形態の検出処理（図１５参照）のステップＳ１４２の断層画像５７を、合成断層画像５９に置き換えた処理を実行する。
【０１３４】
具体的には、検出部２６は、図２３Ｂに示すように、合成断層画像５９及び断層画像５７ｎの各々を、飛跡候補角度（θ，ω）に応じたシフト量でずらした状態において、合成断層画像５９及び断層画像５７ｎの全てが重なる領域を飛跡Ｔの検出対象とする。また、検出部２６は、検出対象の領域における、合成断層画像５９及び断層画像５７ｎの画素５８の画素値を加算する。
【０１３５】
ステップＳ１４３以降の、ステップＳ１４４〜Ｓ１５６の処理は第１実施形態の検出処理と同様であるため、説明を省略する。
【０１３６】
このように本実施形態の飛跡検出装置１２では、検出部２６が、予め定められた数の隣接する断層画像５７を組み合わせて得られた複数の合成断層画像５９に含まれる黒点に基づいて、飛跡Ｔを検出する。飛跡Ｔの検出に合成断層画像５９を用いるため、上記各実施形態のように断層画像５７を用いる場合に比べて、画像の数を約半数にすることができる。
【０１３７】
従って、本実施形態の飛跡検出装置１２によれば、検出精度を低下させることなく、飛跡を検出するための演算量をより削減することができる
【０１３８】
なお、本実施形態では、隣接する２つの断層画像５７を合成して合成断層画像５９とする形態について説明したが、合成する断層画像５７の数は、２つに限定されない。なお、合成する断層画像５７の数が多くなるほど、合成する断層画像５７のうちの最上層の断層画像５７１と最下層の断層画像５７ｎにおける、飛跡Ｔが通過する位置のずれ量が大きくなる。飛跡Ｔの通過する位置のずれ量が大きくなるほど、合成断層画像５９の数は、多くなり、例えば、演算のためのメモリに必要な量が多くなる。また、隣接する断層画像５７同士の間隔に応じても、合成する断層画像５７のうちの最上層の断層画像５７１と最下層の断層画像５７ｎにおける、飛跡Ｔが通過する位置のずれ量は変化し、隣接する断層画像５７の間隔が広くなるほど、ずれ量は大きくなる。
【０１３９】
従って、いくつの断層画像５７を合成して合成断層画像５９とするかは、断層画像群５６に含まれる断層画像５７の数、隣接する断層画像５７同士の間隔、及びメモリの容量等に応じて選択することが好ましい。
【０１４０】
以上説明したように、各実施形態の飛跡検出装置１２は、入射されたμ粒子が透過することにより、μ粒子の飛跡Ｔを潜像として記録する乳剤層５１における、飛跡Ｔを検出する飛跡検出装置１２である。飛跡検出装置１２は、乳剤層５１における深さ方向が異なる複数の断面の各々を撮像して得られ、かつμ粒子の入射側に最も近い乳剤層５１の断面によって得られた断層画像５７１と、μ粒子の入射側に最も遠い乳剤層５１の断面によって得られた断層画像５７ｎとを含み、各々複数の画素５８が配置された断層画像群５６を取得する断層画像取得部２０を備える。また、飛跡検出装置１２は、断層画像５７１の基準点６０を通る方向が各々異なり、かつ各々が断層画像５７ｎの異なる画素５８を通る複数の飛跡候補線Ｃを特定する飛跡候補線特定部２４と、乳剤層５１における飛跡Ｔを、複数の飛跡候補線Ｃから検出する検出部２６と、を備える。
【０１４１】
従って、上記構成を有する各実施形態の飛跡検出装置１２によれば、検出精度を低下させることなく、飛跡を検出するための演算量を削減することができる。
【０１４２】
また、上記実施形態において、例えば、断層画像取得部２０、前処理部２２、飛跡候補線特定部２４、検出部２６、及び出力部２８といった各種の処理を実行する処理部（processing unit）のハードウェア的な構造としては、次に示す各種のプロセッサ（processor）を用いることができる。上記各種のプロセッサには、前述したように、ソフトウェア（プログラム）を実行して各種の処理部として機能する汎用的なプロセッサであるＣＰＵに加えて、ＦＰＧＡ（Fieid Programmable Gate Array）等の製造後に回路構成を変更可能なプロセッサであるプログラマブルロジックデバイス（Programmable Logic Device：ＰＬＤ）、ＡＳＩＣ（Application Specific Integrated Circuit）等の特定の処理を実行させるために専用に設計された回路構成を有するプロセッサである専用電気回路等が含まれる。
【０１４３】
１つの処理部は、これらの各種のプロセッサのうちの１つで構成されてもよいし、同種又は異種の２つ以上のプロセッサの組み合わせ（例えば、複数のＦＰＧＡの組み合わせや、ＣＰＵとＦＰＧＡとの組み合わせ）で構成されてもよい。また、複数の処理部を１つのプロセッサで構成してもよい。
【０１４４】
複数の処理部を１つのプロセッサで構成する例としては、第１に、クライアント及びサーバ等のコンピュータに代表されるように、１つ以上のＣＰＵとソフトウェアの組み合わせで１つのプロセッサを構成し、このプロセッサが複数の処理部として機能する形態がある。第２に、システムオンチップ（System On Chip：ＳｏＣ）等に代表されるように、複数の処理部を含むシステム全体の機能を１つのＩＣ（Integrated Circuit）チップで実現するプロセッサを使用する形態がある。このように、各種の処理部は、ハードウェア的な構造として、上記各種のプロセッサの１つ以上を用いて構成される。
【０１４５】
更に、これらの各種のプロセッサのハードウェア的な構造としては、より具体的には、半導体素子などの回路素子を組み合わせた電気回路（circuitry）を用いることができる。
【０１４６】
また、上記実施形態では、飛跡検出プログラム３７が記憶部３２に予め記憶（インストール）されている態様を説明したが、これに限定されない。飛跡検出プログラム３７は、ＣＤ−ＲＯＭ（Compact Disc Read Only Memory）、ＤＶＤ−ＲＯＭ（Digital Versatile Disc Read Only Memory）、及びＵＳＢ（Universal Serial Bus）メモリ等の記録媒体に記録された形態で提供されてもよい。また、飛跡検出プログラム３７は、ネットワークを介して外部装置からダウンロードされる形態としてもよい。
【符号の説明】
【０１４７】
１  非破壊検査システム
１０  断層画像生成装置
１２  飛跡検出装置
１４  解析装置
２０  断層画像取得部
２２  前処理部
２４  飛跡候補線特定部
２６  検出部
２８  出力部
３０  ＣＰＵ
３１  メモリ
３２  記憶部
３３  表示部
３４  入力部
３６  ネットワークＩ／Ｆ
３７  飛跡検出プログラム
３８  飛跡角度データ
３９  バス
５０  原子核乾板
５１、５１Ａ、５１Ｂ  乳剤層
５２  フィルムベース
５６、５６Ａ、５６Ｂ   断層画像群
５７、５７１〜５７ｎ、５７ｍ、５７ｍ＋１  断層画像
５８、５８（ｘ１，ｙ１）、５８（ｘｎ，ｙｎ）  画素
５９、５９Ａ〜５９Ｉ、５９１、５９２  合成断層画像
６０、６３  基準点
６１  中心点
６４  点線
７０  一覧表
７２  ＩＤ配列表
Ｃ  飛跡候補線
Ｎ  ネットワーク
Ｔ  飛跡
ＶＣ  仮想飛跡候補線
θ  天頂角
ω  面内角
【図１】
2020165918.tif 000003
【図２】
2020165918.tif 000004
【図３】
2020165918.tif 000005
【図４】
2020165918.tif 000006
【図５】
2020165918.tif 000007
【図６】
2020165918.tif 000008
【図７】
2020165918.tif 000009
【図８】
2020165918.tif 000010
【図９】
2020165918.tif 000011
【図１０】
2020165918.tif 000012
【図１１】
2020165918.tif 000013
【図１２】
2020165918.tif 000014
【図１３】
2020165918.tif 000015
【図１４】
2020165918.tif 000016
【図１５】
2020165918.tif 000017
【図１６Ａ】
2020165918.tif 000018
【図１６Ｂ】
2020165918.tif 000019
【図１７】
2020165918.tif 000020
【図１８Ａ】
2020165918.tif 000021
【図１８Ｂ】
2020165918.tif 000022
【図１９】
2020165918.tif 000023
【図２０】
2020165918.tif 000024
【図２１】
2020165918.tif 000025
【図２２】
2020165918.tif 000026
【図２３Ａ】
2020165918.tif 000027
【図２３Ｂ】
2020165918.tif 000028
【手続補正書】
【提出日】平成31年4月22日(2019.4.22)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】００４７
【補正方法】変更
【補正の内容】
【００４７】
飛跡Ｔの飛跡角度（θ，ω）が（０，０）以外の場合、換言すると、飛跡Ｔが断層画像５７に斜めに入射した場合、一例として図７に示すように、断層画像群５６の各断層画像５７を通る飛跡Ｔの各断層画像５７内における位置が異なる。一例として図７では、断層画像群５６が４枚の断層画像５７（ｎ＝４）を含む場合を示している。図７に示す例では、断層画像５７１の画素５８（ｘ１，ｙ１）の基準点６０を通る飛跡Ｔは、断層画像５７２、断層画像５７３、及び断層画像５７４の各々における基準点６０に対応する位置と異なる位置であり、かつ各断層画像５７で異なる位置を通る。なお、本実施形態では、一例として基準点６０を、画素５８の中心点としている。