(57)【特許請求の範囲】
【請求項１】
被検査液体物品の垂直透過画像をＤＲ結像により取得し、該垂直透過画像の少なくとも１行を指定し、前記指定された行に対応する位置において、前記被検査液体物品に対して平面断層ＣＴ走査を用いてダブルエネルギーＣＴ走査を行い、ダブルエネルギーＣＴの投影データを得るステップであって、ここで、前記被検査液体物品の垂直透過画像をＤＲ結像により取得することは、前記放射線源から放射線を放射させ、前記被検査液体物品を透過させて高エネルギー透過画像と低エネルギー透過画像を形成し、前記高エネルギー透過画像と前記低エネルギー透過画像とを結合して前記垂直透過画像を形成するか、又は、前記高エネルギー透過画像と前記低エネルギー透過画像の中から一つを選択して前記垂直透過画像とする、前記ステップと、
前記投影データに基づいてＣＴ再構成を行い、前記被検査液体物品の複数の物理属性値を表すＣＴ画像を得るステップと、
前記ＣＴ画像に基づいて前記被検査液体物品の複数の物理属性値を取り出すステップと、
前記抽出された前記被検査液体物品の前記複数の物理属性値を各成分とする多次元上の点が、当該多次元空間における所定の区域に含まれるか否かに基づいて前記被検査液体物品が危険物であるかどうかを判断するステップと、を含めることを特徴とするダブルエネルギーＣＴによる液体物品の検査方法。
【請求項２】
前記複数の物理属性値は被検査液体物品の密度と原子番号を含めることを特徴とする請求項１記載のダブルエネルギーＣＴによる液体物品の検査方法。
【請求項３】
前記垂直透過画像の少なくとも１行を取得することは、オペレータに入力装置により当該少なくとも一行を指定させるか、又は、前記垂直透過画像の垂直方向に対して画素差を計算することにより、少なくとも一つの液体層の位置を抽出して、前記少なくとも一つの液体層の垂直方向における中心行を指定することを含む、請求項１記載のダブルエネルギーＣＴによる液体物品の検査方法。
【請求項４】
ダブルエネルギーの投影データから被検査液体物品の複数の物理属性値を表すＣＴ画像を再構成するステップにおいて、
高エネルギーおよび低エネルギー放射線の各々を第１の基材および第２の基材の順で通過させることにより得られる高低エネルギー投影データと、前記高低エネルギー投影データを、前記第１および第２の基材の基材係数に対応付けた情報を記録するルックアップテーブルと、を生成するステップと、
前記ルックアップテーブルに基づいて再構成を行い、被検査液体に対応する二種類の基材係数のＣＴ画像を得るステップと、
前記基材係数を表すＣＴ画像に基づいて、被検査液体物品の物理属性値を表すＣＴ画像を生成するステップとを含むことを特徴とする請求項１記載のダブルエネルギーＣＴによる液体物品の検査方法。
【請求項５】
被検査液体物品の物理属性値を表すＣＴ画像に基づいて被検査液体物品の物理属性値を得るステップにおいて、
前記ＣＴ画像から液体部分に対応する画素を取り出すステップと、
液体部分の画素から密度平均値と原子番号平均値を算出して、前記被検査液体物品の密度と原子番号とするステップとを含むことを特徴とする請求項１記載のダブルエネルギーＣＴによる液体物品の検査方法。
【請求項６】
前記物理属性値に基づいて前記被検査液体物品が危険物であるかどうかを判断するステップにおいて、
前記密度と原子番号で確定した点が密度−原子番号の二次元座標空間における予定区域にあるかどうかを判断するステップと、
前記点は前記予定区域にあった場合、前記被検査液体物品が危険物であると判断するステップとを含むことを特徴とする請求項１記載のダブルエネルギーＣＴによる液体物品の検査方法。
【請求項７】
各位置に対してダブルエネルギーＣＴ走査を行った後、被検査液体物品のＣＴ画像を回転し、それを一回目のダブルエネルギーＣＴ走査で形成した画像と合わせることを特徴とする請求項１記載のダブルエネルギーＣＴによる液体物品の検査方法。
【請求項８】
各行に対してダブルエネルギーＣＴ走査を行った後、被検査液体物品を回転し、それを走査前の位置と同じようにすることを特徴とする請求項７記載のダブルエネルギーＣＴによる液体物品の検査方法。
【請求項９】
前記液体物品は複数の空間に分けられた容器に置かれていることを特徴とする請求項１記載のダブルエネルギーＣＴによる液体物品の検査方法。
【請求項１０】
所定のパターン画像とのマッチングによって、容器の存在を自動的に計測するステップと、
容器が存在する場合、ＣＴ画像の中の特定マークを計測するステップと、
前記特定マークに基づいて容器を予定の位置に回転するステップとを含むことを特徴とする請求項９記載のダブルエネルギーＣＴによる液体物品の検査方法。
【請求項１１】
さらに、被検査液体物品についての判断結果をスクリーンに表示するステップを含むことを特徴とする請求項８または１０記載のダブルエネルギーＣＴによる液体物品の検査方法。
【請求項１２】
さらに、各被検査液体物品についての判断結果をプリントアウトするステップを含むことを特徴とする請求項８または１０記載のダブルエネルギーＣＴによる液体物品の検査方法。
【請求項１３】
さらに、各被検査液体物品のＣＴ画像をカラーにするステップを含むことを特徴とする請求項８または１０記載のダブルエネルギーＣＴによる液体物品の検査方法。
【請求項１４】
放射線を出す放射線源と、
少なくとも一件の被検査液体物品を透過した放射線信号を計測採集する計測採集装置と、
前記計測採集装置を制御することによって、検査される液体物の垂直透過画像をＤＲ結像により取得し、該垂直透過画像の少なくとも１行を指定し、前記放射線源と計測採集装置を制御し、前記指定された行に対応する位置において、被検査液体物品に対して平面断層ＣＴ走査を用いてダブルエネルギーＣＴ走査を行って投影データを得るコントローラーであって、ここで、前記被検査液体物品の垂直透過画像をＤＲ結像により取得することは、前記放射線源から放射線を放射させ、前記被検査液体物品を透過させて高エネルギー透過画像と低エネルギー透過画像を形成し、前記高エネルギー透過画像と前記低エネルギー透過画像とを結合して前記垂直透過画像を形成するか、又は、前記高エネルギー透過画像と前記低エネルギー透過画像の中から一つを選択して前記垂直透過画像とする、コントローラーと、
投影データから被検査液体物品の複数の物理属性値のＣＴ画像を再構成する装置と、
前記抽出された前記被検査液体物品の前記複数の物理属性値を各成分とする多次元上の点が、当該多次元空間における所定の区域に含まれるか否かに基づいて前記被検査液体物品が危険物であるかどうかを判断する装置と、を備えることを特徴とするダブルエネルギーＣＴによる液体物品の検査設備。
【請求項１５】
前記複数の物理属性値は少なくとも被検査液体物品の密度と原子番号を含むことを特徴とする請求項１４記載のダブルエネルギーＣＴによる液体物品の検査設備。
【請求項１６】
前記垂直透過画像の中から少なくとも一行を指定する装置は、
オペレータにより当該少なくとも一行を指定させる入力装置と、
前記垂直透過画像の垂直方向に対して画素差を計算することにより、前記垂直透過画像を分層する装置と、
ダブルエネルギーＣＴ走査を行う行として各層の中の中間行を指定する装置とを含むことを特徴とする請求項１４記載のダブルエネルギーＣＴによる液体物品の検査設備。
【請求項１７】
投影データから被検査液体物品の複数の物理属性値のＣＴ画像を再構成する装置は、
被検査液体物品の密度を示す密度画像と被検査液体物品の原子番号を示す原子番号画像を結合して前記ＣＴ画像とする装置と、
前記ＣＴ画像から液体部分に対応する画素を取り出す装置と、
液体部分の画素から密度平均値と原子番号平均値を算出して前記被検査液体物品の密度と原子番号とする装置と、を含むことを特徴とする請求項１４記載のダブルエネルギーＣＴによる液体物品の検査設備。
【請求項１８】
前記物理属性値に基づいて前記被検査液体物品が危険物であるかどうかを判断する装置は、
前記密度と原子番号で確定した点が密度−原子番号の二次元座標空間における予定区域にあるかどうかを判断する装置を含め、前記点は前記予定区域にあった場合、前記被検査液体物品が危険物であると判断することを特徴とする請求項１４記載のダブルエネルギーＣＴによる液体物品の検査設備。
【請求項１９】
さらに、各行に対してダブルエネルギーＣＴ走査を行った後、被検査液体物品のＣＴ画像を回転し、それを一回目のダブルエネルギーＣＴ走査で形成した画像と合わせる装置を含むことを特徴とする請求項１４記載のダブルエネルギーＣＴによる液体物品の検査設備。
【請求項２０】
さらに、各行に対してダブルエネルギーＣＴ走査を行った後、被検査液体物品を回転し、それを走査前の位置と同じようにする装置を含むことを特徴とする請求項１４記載のダブルエネルギーＣＴによる液体物品の検査設備。
【請求項２１】
さらに、前記液体物品をそれぞれ置くため複数の空間に分けられた容器を含むことを特徴とする請求項１４記載のダブルエネルギーＣＴによる液体物品の検査設備。
【請求項２２】
さらに、
所定のパターン画像とのマッチングによって、容器の存在を自動的に計測する装置と、
容器が存在する場合、ＣＴ画像の中の特定マークを計測する装置と、
前記特定マークに基づいて容器を予定の位置に回転する装置と、を含むことを特徴とする請求項２１記載のダブルエネルギーＣＴによる液体物品の検査設備。
【請求項２３】
さらに、被検査液体物品についての判断結果を表示する表示装置を含むことを特徴とする請求項２０または２２記載のダブルエネルギーＣＴによる液体物品の検査設備。
【請求項２４】
さらに、各被検査液体物品についての判断結果をプリントアウトする装置を含むことを特徴とする請求項２０または２２記載のダブルエネルギーＣＴによる液体物品の検査設備。
【請求項２５】
さらに、各被検査液体物品のＣＴ画像をカラーにする装置を含むことを特徴とする請求項２０または２２記載のダブルエネルギーＣＴによる液体物品の検査設備。
【請求項２６】
さらに、前記被検査液体物品を支持する支持構造を含め、前記被検査液体物品を支持する前記支持構造の表面は、検査員で識別できる複数の区域に分けられることを特徴とする請求項１４記載のダブルエネルギーＣＴによる液体物品の検査設備。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、照射検査技術の分野に関するものであり、特に、液体物品に対してダブルエネルギーＣＴ結像を行って液体物品を快速かつ安全的に検査する方法及び設備に関するものである。
【背景技術】
【０００２】
米国の９１１事件後、航空領域において安全検査が益々重視されて来た。既存の手荷物安全検査のほかに、さらに手持ちの液体物品の安全検査が要求されている。従って、有効的な方法と設備によって、手荷物の中の液体物品を快速かつ安全的に検査することが急務となっている。
【０００３】
今では、液体物品の安全検査において、化学方法、電磁気方法、中性子方法と放射線方法の四種類の方法がある。その内容が下記の通りである。
【０００４】
１）化学方法は、さらに匂い識別、イオン走査チェックと物質分析に分けられる。実際、液体物品が密封されているから、匂い識別を行えないことがよくある。イオン走査チェックは高い感性で知られているが、背景環境の影響を受けやすく、誤報率が高い。物質分析は、高精度と的確性を有するが、一定の時間がかかるので、現場での快速検査には向かない。
【０００５】
２）電磁気方法は、主動的な測定方法であり、異なる液体物品の電磁波に対する誘電率の差を利用して液体物品を区別する。電磁気方法は金属包装や厚い材料包装から不利な影響を受けやすい。そして、包装材料が複雑となる今では、電磁気方法の応用が制限されている。
【０００６】
３）中性子検査方法の利用において、「中性子活性化」という現象が生じる。即ち、中性子の照射を受けた被検査物液体物品が放射能を残すことはある。それに、中性子の透過能力がもっと強いので、その遮蔽技術がさらに複雑で、遮蔽設備も大型化となり、従って、旅客航空の安全検査システムには向かない。
【０００７】
４）今、航空領域での多くの安全検査装置が放射線装置であり、その中でもっとも広く利用されている技術がＸ線二次元画像技術と三次元ＣＴ走査画像技術である。これらの技術は主に荷物物品の安全検査に使われ、特に荷物物品の中の液体物品について安全検査を行えない。
【０００８】
Ｘ線二次元画像技術では、被検査物体の三次元情報をＸ線入射方向で積分して、物体の二次元情報の投影画像を形成する。この画像は明暗段階や擬似カラーでの表現方法で物体内部の構造を表示し、設備操作者に対し直感的な表示を与える。しかしながら、Ｘ線二次元画像技術においては、どうしても被検査物体の一つの次元の情報を無くしており、この技術を利用して液体物品を検査するとき、被検査液体物品の外形や寸法からかなり影響を受ける。
【０００９】
三次元ＣＴ走査画像技術はＣＴ技術の延長と応用である。ＣＴ技術は、最初に医学に利用され、医者の診断を助けるツールとなっている。その技術は、被検査物の各断層を多角度に投影させることによって実現される。コンピュータを利用して前記各断層の多角度の投影データを再構成し、各断層の再構成画像を計算して、異なる明暗段階で再構成画像における異なる減衰率情報を表示して、被検査物の内部構造を表示する。ＣＴ技術の発展に伴って、非破壊検査の工業ＣＴや安全検査の荷物ＣＴなどが現れた。これらのＣＴ延長技術は、表示される物体の内部構造の断層画像を獲得することを目的とし、手作業での構造分析に依存して、被検査物が危険物品であるかどうかを判断する。だから、自動的、快速的検査には向かない。そのほか、前記ＣＴ型安全検査設備は、カバーする被検査物範囲が広いため、コストが高い、装置が重い、その応用を広げるのが困難である。
【００１０】
前述した通り、液体物品に対して快速検査を行うため、化学方法、電磁気方法、及び中性子方法では、快速安全検査に適合しない特徴があり、Ｘ線二次元画像技術と三次元ＣＴ走査画像技術では、物品内部構造を表現する画像のみを獲得し、液体物品に対する自動安全検査のための充分なデータを提供することができない。
【発明の開示】
【発明が解決しようとする課題】
【００１１】
本発明は、前述した既存技術の不足点を克服するため、放射線を利用して液体物品の梱包を破壊しない状態で快速検査を行い、被検査液体物品の定量的情報を得て、液体物品に対して安全検査を行う方法と設備を提供することを目的とする。
【課題を解決するための手段】
【００１２】
本発明の一つの側面は、被検査液体物品に対してダブルエネルギーＣＴ走査を行い、ダブルエネルギーＣＴの投影データを得るステップと、前記投影データに基づいてＣＴの再構成を行い、被検査液体物品の物理属性値を表すＣＴ画像を得るステップと、前記ＣＴ画像に基づいて被検査液体物品の物理属性値を取り出すステップと、前記の物理属性値に基づいて前記被検査液体物品が危険物であるかどうかを判断するステップとを含めることを特徴とするダブルエネルギーＣＴによる液体物品の検査方法にある。
【００１３】
本発明の実施例によれば、前記物理属性値は被検査液体物品の密度と原子番号を含める。
【００１４】
本発明の実施例によれば、前記ダブルエネルギーＣＴ走査は平面断層ＣＴ走査の方法を利用する。
【００１５】
本発明の実施例によれば、前記ダブルエネルギーＣＴ走査は慣用の螺旋ＣＴ走査の方法を利用する。
【００１６】
本発明の実施例によれば、前記ダブルエネルギーＣＴ走査は大ピッチ螺旋ＣＴ走査の方法を利用する。
【００１７】
本発明の実施例によれば、平面断層ＣＴ走査を行う前に、予め一組の走査位置を設定する。
【００１８】
本発明の実施例によれば、平面断層ＣＴ走査を行う前に、ＤＲ走査を行い、被検査部品の透過画像を得て、そして透過画像によって走査位置を確定する。
【００１９】
本発明の実施例によれば、透過画像を得た後、検査員が入力装置により透過画像の中から少なくとも一行を指定し、ＣＴ走査位置とする。
【００２０】
本発明の実施例によれば、透過画像を得た後、画像処理技術により透過画像の中から少なくとも一行を自動的に確定し、ＣＴ走査位置とする。
【００２１】
本発明の実施例によれば、透過画像の形成は、放射線源からの高エネルギー放射線と低エネルギー放射線が被検査物体を透過して、高エネルギー透過画像と低エネルギー透過画像を形成するステップと、高エネルギー透過画像と低エネルギー透過画像を結合して、前記透過画像を形成するステップとを含める。
【００２２】
本発明の実施例によれば、透過画像の形成は、放射線源からの高エネルギー放射線と低エネルギー放射線が被検査物体を透過して、高エネルギー透過画像と低エネルギー透過画像を形成するステップと、高エネルギー透過画像と低エネルギー透過画像の中から一つを選んで、前記透過画像とするステップとを含める。
【００２３】
本発明の実施例によれば、ダブルエネルギーの投影データから被検査液体物品の物理属性値を表すＣＴ画像を再構成するステップにおいて、高低エネルギー投影データに基づいて二種類の基材係数の投影データを生成ステップと、前記二種類の基材係数の投影データに基づいて再構成を行い、被検査液体に対応する二種類の基材係数のＣＴ画像を得るステップと、前記基材係数を表すＣＴ画像に基づいて、被検査液体物品の物理属性値を表すＣＴ画像を生成するステップとを含める。
【００２４】
本発明の実施例によれば、被検査液体物品の物理属性値を表すＣＴ画像に基づいて被検査液体物品の物理属性値を得るステップにおいて、前記ＣＴ画像から液体部分に対応する画素を取り出すステップと、液体部分の画素から密度平均値と原子番号平均値を算出して前記被検査液体物品の密度と原子番号とするステップとを含める。
【００２５】
本発明の実施例によれば、前記物理属性値に基づいて被検査液体物品が危険物であるかどうかを判断するステップにおいて、前記密度と原子番号で確定した点が密度−原子番号の二次元座標空間における予定区域にあるかどうかを判断するステップと、前記点は前記予定区域にあった場合、前記被検査液体物品が危険物であると判断するステップとを含める。
【００２６】
本発明の実施例によれば、各位置に対してダブルエネルギーＣＴ走査を行った後、被検査液体物品のＣＴ画像を回転し、それを一回目のダブルエネルギーＣＴ走査で形成した画像と合わせる。
【００２７】
本発明の実施例によれば、各行に対してダブルエネルギーＣＴ走査を行った後、被検査液体物品を回転し、それを走査前の位置と同じようにする。
【００２８】
本発明の実施例によれば、前記液体物品は複数の空間に分けられたバレルに置かれている。
【００２９】
本発明の実施例によれば、前記方法は、さらに予定のパターンによりバレルの存在を自動的に計測するステップと、バレルが存在する場合、ＣＴ画像の中の特定マークを計測するステップと、前記特定マークに基づいてバレルを予定の位置に回転するステップとを含める。
【００３０】
本発明の実施例によれば、前記方法は、さらに、被検査液体物品についての判断結果をスクリーンに表示するステップを含める。
【００３１】
本発明の実施例によれば、前記方法は、さらに、各被検査液体物品についての判断結果をプリントアウトするステップを含める。
【００３２】
本発明の実施例によれば、前記方法は、さらに、各被検査液体物品のＣＴ画像をカラーにするステップを含める。
【００３３】
本発明の実施例によれば、前記ダブルエネルギーＣＴ走査は予定の位置に基づいて行われる。
【００３４】
本発明のもう一つの側面は、放射線を出す放射線源と、少なくとも一件の被検査液体物品を透過した放射線信号を計測採集する計測採集装置と、前記放射線源と計測採集装置を制御し、被検査液体物品に対してダブルエネルギーＣＴ走査を行って投影データを得るコントローラーと、投影データから被検査液体物品の少なくとも一種類の物理属性値のＣＴ画像を再構成する装置と、前記物理属性値の相応する次元の空間での位置に基づいて前記被検査液体物品が危険物であるかどうかを判断する装置とを備えることを特徴とするダブルエネルギーＣＴによる液体物品の検査設備にある。
【００３５】
本発明の実施例によれば、前記ダブルエネルギーＣＴ走査は予定の位置に基づいて行われる。
【００３６】
本発明の実施例によれば、前記計測採集装置は少なくとも一件の被検査液体物品を透過した放射線信号を計測採集して透過画像を形成し、前記設備はさらに透過画像の中から少なくとも一行を指定する装置を含め、指定した行に基づいて前記ダブルエネルギーＣＴ走査を行う。
【００３７】
本発明の実施例によれば、前記物理属性値は少なくとも被検査液体物品の密度と原子番号を含める。
【００３８】
本発明の実施例によれば、放射線源からの高エネルギー放射線と低エネルギー放射線が被検査物体を透過して、高エネルギー透過画像と低エネルギー透過画像を形成し、前記設備はさらに高エネルギー透過画像と低エネルギー透過画像を結合して、前記透過画像を形成する装置を含める。
【００３９】
本発明の実施例によれば、放射線源からの高エネルギー放射線と低エネルギー放射線が被検査物体を透過して、高エネルギー透過画像と低エネルギー透過画像を形成し、前記設備はさらに高エネルギー透過画像と低エネルギー透過画像の中から一つを選んで、前記透過画像とする装置を含める。
【００４０】
本発明の実施例によれば、前記の透過画像の中から少なくとも一行を指定する装置は、検査員が入力装置により透過画像の中から少なくとも一行を選ぶ装置を含める。
【００４１】
本発明の実施例によれば、前記の透過画像の中から少なくとも一行を指定する装置は、前記透過画像を分層するため、前記透過画像の画素値を分析する装置と、ダブルエネルギーＣＴ走査を行う行として各層の中の中間行を指定する装置とを含める。
【００４２】
本発明の実施例によれば、投影データから被検査液体物品の物理属性値のＣＴ画像を再構成する装置は、被検査液体物品の密度でマークされた密度画像と被検査液体物品の原子番号でマークされた原子番号画像を結合して前記ＣＴ画像とする装置と、前記ＣＴ画像から液体部分に対応する画素を取り出す装置と、液体部分の画素から密度平均値と原子番号平均値を算出して前記被検査液体物品の密度と原子番号とする装置とを含める。
【００４３】
本発明の実施例によれば、前記物理属性値に基づいて前記被検査液体物品が危険物であるかどうかを判断する装置は、前記密度と原子番号で確定した点が密度−原子番号の二次元座標空間における予定区域にあるかどうかを判断する装置含め、前記点は前記予定区域にあった場合、前記被検査液体物品が危険物であると判断する。
【００４４】
本発明の実施例によれば、前記の設備は、さらに、各位置に対してダブルエネルギーＣＴ走査を行った後、被検査液体物品のＣＴ画像を回転し、それを一回目のダブルエネルギーＣＴ走査で形成した画像と合わせる装置を含める。
【００４５】
本発明の実施例によれば、前記の設備は、さらに、各行に対してダブルエネルギーＣＴ走査を行った後、被検査液体物品を回転し、それを走査前の位置と同じようにする装置を含める。
【００４６】
本発明の実施例によれば、前記の設備は、さらに、前記液体物品をそれぞれ置くため複数の空間に分けられたバレルを含める。
【００４７】
本発明の実施例によれば、前記の設備は、さらに、予定のパターンによりバレルの存在を自動的に計測する装置と、バレルが存在する場合、ＣＴ画像の中の特定マークを計測する装置と、前記特定マークに基づいてバレルを予定の位置に回転する装置とを含める。
【００４８】
本発明の実施例によれば、前記の設備は、さらに、被検査液体物品についての判断結果を表示する表示装置を含める。
【００４９】
本発明の実施例によれば、前記の設備は、さらに、各被検査液体物品についての判断結果プリントアウトする装置を含める。
【００５０】
本発明の実施例によれば、前記の設備は、さらに、各被検査液体物品のＣＴ画像をカラーにする装置を含める。
【００５１】
本発明の実施例によれば、前記の設備は、さらに、前記被検査液体物品を支持する支持構造を含め、前記被検査液体物品を支持する前記支持構造の表面は、検査員で識別できる複数の区域に分けられる。
【００５２】
本発明のさらにもう一つの側面は、放射線を出す放射線源と、少なくとも一件の被検査液体物品を透過した放射線信号を計測採集する計測採集装置と、前記放射線源と計測採集装置を制御し、被検査液体物品に対して螺旋ＣＴ走査を行ってそれぞれ被検査液体物品の少なくとも一種類の物理属性値を表す一組の螺旋ＣＴ画像を形成するコントローラーと、この一組の螺旋ＣＴ画像を分析して液体の螺旋ＣＴ画像部分を確定する装置と、液体の螺旋ＣＴ画像部分に含まれる物理属性値の相応する次元の空間での位置に基づいて被検査液体物品が危険物であるかどうかを判断する装置とを備えることを特徴とするダブルエネルギーＣＴによる液体物品の検査設備にある。
【発明の効果】
【００５３】
本発明の方法と設備によれば、透過画像を利用してダブルエネルギーＣＴ走査を誘導しているので、計測速度を高めると共に、計測精度を低下させないし、また、透過画像から液体物品にサンドイッチ構造があるかどうかを判断できる。
【００５４】
さらに、測定した密度と原子番号が二次元空間における予定区域にあるかどうかを判断することによって、被検査液体物品の中に爆発物があるかどうかを判断できる。
【００５５】
さらに、検査員は任意の指定位置でダブルエネルギーＣＴ走査を行うことができるので、再検査に便利である。
【００５６】
さらに、格子分けのバレルを利用するので、複数の被検査物がある場合、その被検査液体物品が爆発物であるかどうかを判断しやすい。
【発明を実施するための最良の形態】
【００５７】
以下、本発明の望ましい実施形態を、図面を参照しながら説明する。図面において、異なる図でも、同じまたは類似部品に対しては同じマークを使って表している。本発明の主旨を鮮明にするため、すでに既知の機能や構造について詳しい説明を省略する。
【００５８】
［第１実施形態］
図１は本発明の実施形態による検査設備の構造を示した図である。
【００５９】
図１に示すように、本実施形態による検査設備は、Ｘ線機などのように検査用ダブルエネルギーＸ線を出す放射線源１０と、被検査液体物品を上下昇降可能にまたはＺ軸に回転可能に支持して計測区域に進入させることができ、放射線源１０からの放射線を被検査液体物品に透過させる支持構造４０と、モジュール構造の計測器とデータ採集器を有し、被検査液体物品を透過したダブルエネルギー放射線を計測してアナログ信号を得て、デジタル信号に転換し、液体物品における高エネルギーＸ線と低エネルギーＸ線の走査データを出力する計測採集装置３０と、システムの各部分の同期的作動を制御するコントローラー５０と、データ採集器で採集したデータを処理し、検査結果を出力する計算機データ処理装置６０と、を備える。
【００６０】
図１に示すように、被検査液体物品のＤＲデータと多角度投影データを獲得するため、放射線源１０は被検査液体物品を置かれる支持構造４０の一側にあり、計測器とデータ採集器を含める計測採集装置３０は支持構造４０のほかの側にある。計測採集装置３０に含まれるデータ増幅成形回路が（電流）積分方式またはパルス（カウント）方式で作動する。計測採集装置３０のデータ出力ケーブルが計算機データ処理装置６０に接続し、採集したデータがトリガーコマンドによって計算機データ処理装置６０に記憶される。
【００６１】
また、検査設備は、さらに支持構造４０の上に設置され、外へのＸ線の漏洩を遮蔽する金属製の円筒状の物体通路２０を含める。被検査液体物品は被検査物体通路の中に置かれる。
【００６２】
図２は図１で示した計算機データ処理装置６０を説明するためのブロック図である。図２に示すように、データ採集器で採集したデータがインターフェースユニット６８とバス６４を通して記憶装置６１に記憶される。ＲＯＭ６２には計算機データ処理装置の配置情報とプログラムを記憶する。ＲＡＭ６３には処理装置６６の作動中において各種類のデータを一時的に記憶する。また、記憶装置６１にはデータ処理のための計算機プログラムを記憶する。バス６４は、前記の記憶装置６１、ＲＯＭ６２、ＲＡＭ６３、入力装置６５、処理装置６６、表示装置６７およびインターフェースユニット６８に接続している。
【００６３】
ユーザがキーボードやマウスなどの入力装置６５から作動コマンドを入力した後、処理装置６６は計算機プログラムのコードによって予定のデータ処理手順を作業し、データ処理の結果を得てから、ＬＣＤデイスプレーなどの表示装置６７に表示したり、プリンターなどのハードコピーの形で直接に出力する。
【００６４】
図３は本発明の実施形態によるコントローラの構造を示したブロック図である。図３に示すように、コントローラ５０は、計算機６０のコマンドによって放射線源１０、支持構造４０と計測採集装置３０を制御する制御ユニット５１と、制御ユニット５１の制御によってトリガーコマンドを生成して、放射線源１０、計測採集装置３０と支持構造４０の作動をトリガーするトリガー信号生成ユニット５２と、制御ユニット５１の制御によってトリガー信号生成ユニット５２が生成したトリガー信号に基づいて支持構造４０を上昇したり降下する第一駆動モーター５５と、支持構造４０の移動に伴って支持構造４０の高さに関する情報を制御ユニット５１へフィードバックする高さ情報獲得ユニット５３と、制御ユニット５１の制御によってトリガー信号生成ユニット５２が生成したトリガー信号に基づいて支持構造４０を回転する第二駆動モーター５６と、支持構造４０の回転において支持構造４０の回転角度を制御ユニット５１へフィードバックする角度情報獲得ユニット５４と、を備える。
【００６５】
本発明の実施形態によれば、前記高さ情報獲得ユニット５３と角度情報獲得ユニット５４は共に光コード式であり、干渉に強いという利点を有する。
【００６６】
本発明の実施例によれば、液体物品の中の液体の密度と原子番号を計測し、またその計測した密度と原子番号が予め設定した危険区域にあるかどうかを判断することによって、危険物を検査できる。密度と原子番号については、安全液体（たとえば、飲料、スキンケアー、化粧品など）と危険液体（たとえば、燃焼性液体、液体爆薬、強腐食性液体など）がかなり異なっている。たとえば、安全液体の水なら、その密度が約１であり（特徴密度が１．１１）、原子番号が約７．５１である。一方、危険液体のアルコールなら、その密度が約０．７９であり（特徴密度が０．８９）、原子番号が約６．４７である。したがって、モード識別方法を利用して、密度と原子番号から安全液体と危険液体を区別できる。
【００６７】
図４（Ａ）、（Ｂ）は本発明の液体物品の検査方法を示した図である。密度と原子番号で構成した二次元空間では、一つや複数の区域を作り、被検査物の密度と原子番号がこれらの区域にある場合、安全液体だと判断し、逆に危険液体だと判断する。安全液体区域の確定のし方として、沢山の安全液体と危険液体を計測して、安全液体が割り合い集中し、危険液体が疎らである区域を安全液体区域とする。図４（Ａ）に示すように、各種液体の密度と原子番号を計測して、密度と原子番号の二次元空間において、安全液体のある区域と危険液体のある区域を確定する。そして、図４（Ｂ）に示すように、安全液体のある区域を縮小して検査感度を向上する。したがって、前記の方法によって液体物品の密度と原子番号を獲得した後、その二次元空間においての位置から危険物であるかどうかを判断できる。
【００６８】
図５はＤＲ結像とＣＴ結像の関係を説明するための図である。本発明の実施形態によれば、まず液体物品に対してＤＲ結像を行い、液体物品の中の液体部分を判断する。その後、液体部分のある位置のみに対してＣＴ結像を行い、検査速度を向上する。
【００６９】
図６と図７はそれぞれＤＲ結像の例を示した。図６に示すように、ある液体物品に対してＤＲ結像を行った後、下記の通り、画素値について分析して液体物品の中の液体の位置を判断する。図６に示すように、この液体物品には一種類の液体しか含まれない。しかし、図７に示すように、異なる液体の液体に対する吸収係数が違うので、液体物品に二種類または二種類以上の液体が存在し、液体の中にいくつかの層に分かれる場合、ＤＲ結像によって得たＤＲ画像について分析して液体の間の境目の位置を判断する。そして、異なる液体に対してＣＴ結像を行う。
【００７０】
図８は本発明の第１実施形態による液体物品の検査方法の全体のフローチャートである。図８に示すように、たとえば、税関を通るとき、乗客が携帯している液体物品に対して安全検査をする。まず、ステップＳ１１１で、検査員が被検査液体物品を支持構造４０に載せる。
【００７１】
そして、図６と図７に示すように、ステップＳ１１２で、検査員がスタートボタンを押してＤＲ走査を行い、ＤＲ画像を生成する。
【００７２】
前述した通り、ＤＲ走査の目的は、被検査液体物品の透過画像を獲得し、被検査液体の内部構造をはっきりする事と、ソフトウェアによってＤＲ画像から液体位置を自動的に識別し、次のＣＴ結像を誘導する事と、ユーザがＤＲ画像においてＣＴ結像と検査を行う位置を指定し、次のＣＴ結像を誘導する事と、である。後にＤＲ結像の過程を詳しく説明する。
【００７３】
特に説明すべきは、必ずしもＤＲ走査をする必要がない。検査速度を向上するため、ＤＲ走査によってＣＴ走査を誘導することではなく、直接的に若干の位置を予定してＣＴ走査を行ってもよい。たとえば、調べた結果として、ボトルの底から５ｃｍのところをＣＴ走査すると、ほとんどの液体商品の液体位置を検査できるから、ボトルの底から５ｃｍの高さを予定走査高さとしてもよい。また、検査員が被検査物の寸法を目で見て、経験から適当な高さを設定してもよい。たとえば、缶詰コーラの走査高さが３ｃｍとし、底の厚い葡萄酒の走査高さが１０ｃｍとする。
【００７４】
ＤＲ画像を獲得した後、ステップＳ１１３Ａで、ＤＲ画像を分析し、ＣＴ結像の走査位置を自動的に確定する。あるいは、ステップＳ１１３Ｂで、検査員がマウスなどの入力装置６５でＣＴ走査の位置を指定する。これによって、液体物品における代表的な位置だけをＣＴ走査することによって、検査速度を向上できるし、検査品質を下げることはない。
【００７５】
ＣＴ走査の位置を確定した後、ステップＳ１１４でＣＴ走査過程に入る。即ち、前記確定した走査位置で液体物品に対してＣＴ走査を行って、ＣＴ走査データを獲得し、またそれによってＣＴ画像を再構成する。このＣＴ画像の各画素が液体物品の中の対応する部分の密度と原子番号を表している。
【００７６】
その後、ステップＳ１１５で、計算機がプログラムを実行してＣＴ画像を分析し、計測した密度と原子番号を獲得する。さらに、ステップＳ１１６で、計測した密度と原子番号で決めた点が密度と原子番号の二次元空間の予定区域にあるかどうかを判断する。ステップＳ１１７で、密度と原子番号で決めた点が予定の危険区域にあれば、この液体物品が疑わしい、または危険物である可能性があると見て、検査員に警告を出したり、検査結果をプリントアウトする。
【００７７】
次は、図９〜図１４を参照しながら、前記各ステップでの作業過程を詳しく説明する。図９はＤＲ結像過程のフローチャートであり、図１０はＤＲ結像過程における計測採集装置３０で採集したＤＲ画像データの行列方式を示した図である。
【００７８】
図９に示すように、ＤＲ結像過程において、ステップＳ２１０で、計算機６０からコントローラー５０へコマンドを伝送し、支持構造４０を物体通路２０に沿って垂直に移動する。支持構造４０が垂直に移動するとき、コントローラー５０は高さ情報獲得ユニット５３によって支持構造４０の高さを計測して制御する。
【００７９】
ステップＳ２１１で、コントローラー５０はある高さ（たとえば、１ｍｍ）ごとに計測採集装置３０へ一つのトリガー信号を伝送する。計測採集装置３０はこのトリガー信号を受信すると、各検出器からの出力信号を採集し、高エネルギー計測データと低エネルギー計測データを獲得して内蔵のバッファーに記憶する。
【００８０】
ステップＳ２１２で、支持構造４０が指定高さ、たとえば、５００ｍｍに到達したかどうかを判断する。到達しない場合、ステップＳ２１０へ進んで、前記の作業を続ける。
【００８１】
支持構造４０が前記指定高さに到達した場合、コントローラー５０はもう計測採集装置３０にトリガー信号を伝送しない。計算機６０は、計測採集装置３０から採集した高低エネルギー検出器の信号を読み取り、データ行列の形にしてＤＲ画像となる。ＤＲ画像の各画素は、低エネルギー放射線強度と高エネルギー放射線強度を含めて、物体を透過した放射線の残留強度を記録してある。
【００８２】
上述した通り、続けるＣＴ結像の走査位置がすべてＤＲ画像の上で決められる。自動識別方式または人工指定方式では、まずＤＲ画像において一つの行番号を取り、そして計算機がこの行番号を支持構造の高さに転換し、支持構造４０を指定位置へ移動させるようにコントローラー５０にコマンドを出して、ＣＴ結像を行う。
【００８３】
ＤＲ結像の過程から分かるように、ＤＲ画像の各行が支持構造４０の一つ特定の高さに対応する。ＤＲ結像を始めるときの支持構造４０の高さを０として、結像に伴って支持構造４０の高さは徐々に下がり、ｈ  ｍｍごとに一回のデータ採集をトリガーする場合、ＤＲ画像のｍ行目が支持構造４０の高さ−ｍ＊ｈに対応する。
【００８４】
図１１はＤＲ画像を処理してＣＴ走査位置を確定するフローチャートである。
【００８５】
ＤＲ画像において、被検査液体物品が普通ボトルの底、液体部分、ボトルネック、ボトル蓋など分けられる。画像分析によって、その中の液体部分（いくつかの層があるかもしれない）を取り出して、各層のＣＴ走査位置を確定することができる。
【００８６】
ステップＳ３１０で、ＤＲ画像の高低エネルギーデータを結合して平滑化し、雑音の少ない一意ＤＲ画像を得る。高低エネルギーデータを結合する方法としては、たとえば、高低エネルギーデータの中の一つを直接に選んで結合結果としてもよいし、高低エネルギーデータを重みつけて組み合わせてもよい。平滑化する方法としては、ガウスフィルターを利用して画像をフィルター処理する。
【００８７】
ステップＳ３１１で、平滑化されたＤＲ画像から被検査液体物品（フロント）を取り出し、空気（バック）を取り除く。具体的な方法としては、一つの閾値を設定し、この閾値より低い画素をフロント画素とし、ほかの画素をバック画素とする。閾値でバックを取り除く根拠は、被検査液体物品が放射線を遮断したので、それに対応するＤＲ画像の画素値が低くなることである（ＤＲ画像で記録したのが放射線の余剰強度である）。
【００８８】
ステップＳ３１２で、平滑化されたＤＲ画像の水平エッジ画素を取り出す。具体的な方法としては、垂直方法における隣り合う画素の差を算出し、この差がある閾値を越えると、その対応する画素が水平エッジであると認める。
【００８９】
ステップＳ３１３で、平滑化されたＤＲ画像の水平エッジ行を取り出す。水平エッジ行は、ボトル底と液体の界面、液体と空気の界面、ボトル蓋と空気の界面、または多層液体の間の界面に対応する。具体的な方法としては、ＤＲ画像における水平方向の各行に対して、水平エッジ画素数とフロント画素数の比例を統計する。この比例がある閾値（たとえば、５０％）を超えると、その行が水平エッジ行であると認める。
【００９０】
ステップＳ３１４で、ＤＲ画像を垂直に区画し、また、非液体区域を取り除く。ＤＲ画像の各水平エッジ行によって、ＤＲ画像を複数の区域に区画する。下から上へ順次に、ボトル底、液体（異なる密度によって多層があるかもしれない）、ボトル内の空気（もしあれば）、ボトル蓋などを含める。篩い分けルールを作ることによって、非液体区域を取り除ける。具体的なルールが次のとおりである。
【００９１】
（ａ）垂直方向においては、行数がある閾値より小さい区域を排除する。行数の割り合い小さい区域、即ち、厚さの小さい区域がボトル底、ボトル蓋、または容器内の液体と容器トップの間の小さい隙間（たとえば、缶詰の中のトップにある空気）であるかもしれない。各種の液体容器の底、蓋、及び容器内の空気層の厚さを調べて具体的な閾値を確定することができる。
【００９２】
（ｂ）水平方向においては、各行における平均フロント画素数がある閾値より小さい区域を排除する。これらの区域がよく細長いボトルネックに対応する。各種の液体容器のネックの幅を獲得して具体的な閾値を確定することができる。
【００９３】
ステップＳ３１５で、液体区域（複数の場合がある）のＣＴ走査位置を確定して、液体分層の位置決めを完成する。非液体区域を取り除いたら、残る区域（複数の場合がある）が液体区域に属する。高さ方向においてこれらの区域の中心行がＣＴ走査位置とする。
【００９４】
前述したのはＣＴ走査位置の自動的確定の過程である。しかし、手作業で走査位置を確定する場合、作業員が入力装置６５を通して表示されたＤＲ画像の上でＣＴ走査行を直接的に指定し、ＣＴ走査とする。
【００９５】
図１２はＣＴ結像過程を示した図である。図１３はＣＴ結像過程におけるＣＴ投影データの行列方式を示した図である。
【００９６】
図１２に示すように、ＣＴ走査の位置を確定した後、ＣＴ結像過程を行い、即ち、ＣＴ走査の位置でＣＴ結像を行い、被検査物の一つの断層における密度と原子番号の画像を生成して、液体の密度と原子番号を計測する。上述した通り、典型的な位置だけでＣＴ走査するので、税関を通る時間をかなり節約できる。
【００９７】
ステップＳ４１０で、計算機６０がコントローラ５０にコマンドを伝送し、支持構造４０を予定角度だけ、たとえば１度回転する。コントローラ５０は支持構造４０の回転過程において、角度情報獲得ユニット５４によって支持構造４０の角度を計測して制御する。
【００９８】
ステップＳ４１１で、１度回転すると、コントローラ５０は計測採集装置３０に一つのトリガー信号を伝送する。計測採集装置３０はこのトリガー信号を受信すると、各検出器からの出力信号を一度採集して、内蔵のバッファに記憶する。
【００９９】
そして、ステップＳ４１２で、累積回転角度が一周に達したかどうかを判断する。達しない場合、ステップＳ４１０へ進んで、前記の作業を続ける。
【０１００】
累積回転角度が指定角度（たとえば、３６０度）に達した場合、ステップＳ４１３で回転を止め、コントローラ５０はもう計測採集装置３０にトリガー信号を伝送しない。計算機６０は、計測採集装置３０から採集した高低エネルギー検出器の信号を読み取り、図１３に示すように、データ行列の形にして、ＣＴ投影データとなる。ＣＴ投影データにおける各画素が、低エネルギー放射線強度と高エネルギー放射線強度を含めて、物体を透過した放射線の残留強度を記録してある。
【０１０１】
ステップＳ４１４で、計算機６０はダブルエネルギー再構成方法を利用して、高低エネルギーＣＴ投影データに基づいて、断層における密度と原子番号画像、即ち、ＣＴ画像を再構成する。ＣＴ画像における各画素が、当画素の対応する位置での被検査物体の密度と原子番号を記録してある。
【０１０２】
次は、高低エネルギーＣＴ投影データから断層画像を再構成する過程を説明する。
【０１０３】
ＣＴ数学原理
二次元分布関数ｕ（ｘ，ｙ）をある方向θに沿って線積分して、一次元関数ｐθ（ｔ）を得る。この関数ｐθ（ｔ）はθ角度におけるｕ（ｘ，ｙ）の投影と呼ばれる。各方向の投影ｐθ（ｔ）を得たら、Ｒａｄｏｎ変換に従って二次元分布関数ｕ（ｘ，ｙ）を正確に計算できる。投影から二次元分布関数を得る過程は、再構成と言う。
【０１０４】
実際には、Ｘ線機と検出器が物体を回り一周回転して、物体のあるスライスにおける各方向の減衰係数分布の各方向の投影を計測でき、次第に、物体スライスの減衰係数の二次元分布をＣＴ原理によって再構成できる。
【０１０５】
基材分解モデル
小型Ｘ線安全検査システムのエネルギー範囲内（＜２００ＫｅＶ）においては、物質の線減衰係数が次の式１で近似的に表される。
【数１】
0005164777.tif 000002
【数２】
0005164777.tif 000003
【数３】
0005164777.tif 000004
【０１０６】
式１には、ｆp（Ｅ）  はエネルギーによる光電効果面積の変化を表し、ｆKN（Ｅ）はエネルギーによるコンプトン散乱面積の変化を表し、ｆp（Ｅ）  とｆKN（Ｅ）はすでに既存の解析式で表される。定数ａ1とａ2は物質の原子番号、質量数及び密度と関係し、式２と式３で表され、式の中で、Ｚは原子番号、Ｍは質量数、ρは密度（ｇ／ｃｍ3）、ｎは定数である。
【０１０７】
各物質の線減衰係数がすべて式１の二つ定数ａ1とａ2で唯一的に確定されるため、二種類の基材、たとえば、カーボンとアルミニウムを選び、式４のように、基材の線減衰係数の一次結合の組み合わせでその他のすべての材料を表すことができる。
【数４】
0005164777.tif 000005
【０１０８】
その中に、μ（Ｅ）  は任意材料の線減衰係数であり、μ1（Ｅ）  とμ2（Ｅ）は選んだ二種類の基材の線減衰係数である。ｂ1とｂ2を基材係数と呼ぶ。
【０１０９】
式５によって、原子番号の２倍と質量数の比を密度に乗じて特徴密度とする。
【数５】
0005164777.tif 000006
【０１１０】
二種類の基材の原子番号と特徴密度をそれぞれ（Ｚ1，ρ1*）、（Ｚ2，ρ2*）とし、前記式１〜４から、任意材料の原子番号と特徴密度が次の式で表される。
【数６】
0005164777.tif 000007
【数７】
0005164777.tif 000008
【０１１１】
基材投影モデル
Ｘ線源で発生したスペクトルが連続的なものであり、Ｘ線のエネルギーに対する検出器の反応関数も一定ではない。スペクトルＮ（Ｅ）とエネルギー反応関数Ｐd（Ｅ）の乗算の積をＳ（Ｅ）とし、また、Ｓ（Ｅ）を正規化すると
【数８】
0005164777.tif 000009
【０１１２】
そして、一つの投影線における投影値が次の積分式で表される。
【数９】
0005164777.tif 000010
【０１１３】
式９には、Ｉ0とＩは、それぞれ放射線が物体により減衰される前及び減衰された後の検出器でのカウント値を表す。Ｅmは放射線の最大エネルギーであり、ｌは放射線の透過ルートである。
【０１１４】
式９は実のシステムにおける計測値ｐと二次元分布μ（ｘ，ｙ）の関係を表す。この式から分かるように、Ｘ線の多色性のため、式９はμ（ｘ，ｙ）のある直線に沿う線積分を表していない、次第にＣＴ数学原理の要求を満たさない。慣用の再構成方法ではこの不一致性を無視し、再構成により得たμ（ｘ，ｙ）画像はカッピングアーチファクトを含み、いわゆる硬化アーチファクトと呼ぶ。
【０１１５】
既存のダブルエネルギーＣＴ方法では、まず慣用の再構成方法で２組のμ（ｘ，ｙ）を得て、原子番号と密度などの情報を算出する。しかし、これでは放射線の多色性の影響を取り除くことができない。本発明では基材分解のアイデアでこの問題を解決した。
【０１１６】
基材分解モデルを式９に入れて、基材係数に基づいて投影値を表す式が次の通りである。
【数１０】
0005164777.tif 000011
【０１１７】
前式の中のルートｌに沿う積分を下の式で表すと
【数１１】
0005164777.tif 000012
【数１２】
0005164777.tif 000013
【０１１８】
式１１、式１２の定義によれば、Ｂ1とＢ2は基材係数ｂ1（ｘ，ｙ）とｂ2（ｘ，ｙ）の投影を表す。各角度での完全な基材係数投影を獲得すれば、ＣＴ再構成理論で基材係数ｂ1とｂ2の分布を再構成でき、従って、基材分解モデルによって物体の原子番号と特徴密度分布、及び任意エネルギーでの線減衰係数値を算出する。
【０１１９】
基材係数投影を求める。
ダブルエネルギーＣＴにおいて、二つのエネルギーで投影データを採集し、下記のダブルエネルギー投影データを得る。
【数１３】
0005164777.tif 000014
【数１４】
0005164777.tif 000015
【０１２０】
計測で（ｐ1，ｐ2）を得た後、式１３と１４によって（Ｂ1，Ｂ2）を求めることができる。しかし、前記二つの式はすべて対数積分式であり、解析によって解を求められない。非線形繰り返し方法なら計算量が莫大であり、安定的な解を簡単に求められない。
【０１２１】
ただし、放射線が厚さｄ1とｄ2の基材１と２を透過した後、測定したダブルエネルギー投影が次の式で表される。
【数１５】
0005164777.tif 000016
【数１６】
0005164777.tif 000017
【０１２２】
式１３、式１４と式１５、式１６を比べて分かるように、（ｐ1、ｐ2）に対して測定した投影データが同じであれば、基材の投影データは（Ｂ1、Ｂ2）と基材の厚さ組み合わせ（ｄ1、ｄ2）に対してまったく同じである。
【０１２３】
異なる厚さの組み合わせにおけるダブルエネルギー投影を測定して、（ｐ1、ｐ2）へのダブルエネルギー投影データの対応関係、及び（Ｂ1、Ｂ2）への基材係数投影データの対応関係が得られ、参照表を形成する。（ｐ1、ｐ2）から（Ｂ1、Ｂ2）を計算する場合、参照表に基づき、線形補間方法を持って複雑な数学過程を換える。
【０１２４】
図１４は液体属性の測定過程を示した図である。
【０１２５】
図１４に示すように、ステップＳ５１０で、密度画像と原子番号画像を結合して平滑化し、雑音の少ない一意ＣＴ画像を得る。具体的な結合方法としては、密度画像と原子番号画像の中の一つを選んで結合結果としたり、両者を重みつけて組み合わせる。平滑化する方法としては、ガウスフィルターを利用して画像をフィルター処理する。
【０１２６】
ステップＳ５１１で、平滑化されたＣＴ画像から被検査物（フロント、液体及び包装を含める）を取り出し、空気（バック）を取り除く。具体的な方法としては、一つの閾値を設定し、この閾値より高い画素をフロント画素とし、ほかの画素をバック画素とする。これは、空気の密度と原子番号はほぼ０であり、被検査液体物品の密度と原子番号は比較的に大きいためである。
【０１２７】
ステップＳ５１２で、フロント画素から液体画素を取り出す。具体的な方法としては、まず、ＣＴ画像に対応するバイナリ画像を作り、フロント画素の値を１とし、バック画素の値を０とする。そして、このバイナリ画像に対して数回の腐食計算（erosion operation）を行う。毎回一つの画素を腐食する。最後に値１の画素が液体画素となる。これは、液体がいつも包装物の中にあるためである。数回の腐食によって徐々に包装物を取り除ける。包装物の厚さを確定してから、腐食の回数を設定する。
【０１２８】
ステップＳ５１３で、ＣＴ画像におけるすべての液体画素の平均密度と平均原子番号を統計して、今回測定の出力結果とする。
【０１２９】
また、ＤＲ画像分析の過程で多層の液体と確定すれば、各層の液体に対して前記のステップ作業を行い、進んで各層の液体が疑わしいかどうかを判断する。最後に、検査員に最終の検査結果を報告する。
【０１３０】
［第２実施形態］
前記の第１実施形態で紹介したのは毎回一件の液体物品を検査するが、次は、図１５〜１８参照しながら複数の液体物品を検査する過程を説明する。第１実施形態と第２実施形態の区別としては、複数の液体物品が存在するので、ＣＴ結像を完成したら、表示装置で表示されている結像の位置を支持構造の上の物体の位置に対応させて、どの液体物品が危険であるかを検査員に知らせる必要がある。図１５（Ａ）、（Ｂ）は本発明の第２実施形態による複数の液体物品を検査する場合、再構成したＣＴ画像を示した図である。
【０１３１】
たとえば、検査員が下向きの目視状態で支持構造の上の被検査物を観察しているので、検査後の各層のＣＴ画像における各物体の位置が支持構造の上面図と一致しなければならない。
【０１３２】
図１６（Ａ）〜（Ｋ）はＣＴ結像を完成した後、ＣＴ再構成画像および／または支持構造を回転して、それをＣＴ走査前の物体の相対位置と一致させる過程を示した図である。
【０１３３】
図１６（Ａ）は、一回目のＣＴ結像を完成した後の支持構造４０の上面図を示し、その中の矢印が支持構造の角度を示している。図１６（Ｂ）は、一回目のＣＴ画像を回転する前のＣＴ画像を示し、それが一回目の結像の始まりにおける支持構造の上面図と一致し、その中の点線が支持構造の回転範囲をマークする。図１６（Ｃ）は、一層目のＣＴ画像を回転した後の画像を示す。回転角度が０であるので、実際、この層のＣＴ画像を回転する必要がない。
【０１３４】
図１６（Ｄ）は、ｎ回目のＣＴ結像の始まりにおける支持構造の上面図を示す。図１６（Ｄ）から分かるように、ｎ回目のＣＴ結像の始まりにおいて、支持構造は一回目のＣＴ結像の前の状況と比べて一つの角度を有する。図１６（Ｅ）は、ｎ回目のＣＴ画像を回転する前の画像を示し、それがｎ回目の結像の始まりにおける支持構造の上面図と一致する。図１６（Ｆ）は、一層目のＣＴ画像に合わせるため、図１６（Ｄ）の画像を回転することを示している。
【０１３５】
図１６（Ｇ）は、最終回のＣＴ結像の始まりにおける支持構造の上面図を示す。このときの支持構造は一回目のＣＴ結像の始まりの状況と比べて一つの角度を有する。図１６（Ｈ）は、Ｎ回目のＣＴ画像を回転する前の画像を示し、それがＮ回目の結像の始まりにおける支持構造の上面図と一致する。図１６（Ｉ）は、一層目のＣＴ画像に合わせるため、図１６（Ｈ）の画像を回転することを示している。
【０１３６】
図１６（Ｊ）は、最終回のＣＴ結像の終わりにおける支持構造の上面図を示し、一回目のＣＴ結像の始まりの状況と比べて一つの回転角度を有する。図１６（Ｋ）は、一回目のＣＴ結像で得た画像と一致するため、結像終了後、支持構造を回転することを示している。
【０１３７】
すべてのＣＴ結像過程が終わると、各層のＣＴ画像と支持構造に対して角度調整を行う。まず、各回でＣＴ結像の始まりにおける支持構造の角度（この角度は角度情報獲得ユニット５４によって得られる）によって、各層のＣＴ画像を回転して、たとえば一層目のＣＴ画像に合わせるため、各層のＣＴ画像の中にある同一物体の位置を一致させる。そして、支持構造の角度を調整して、支持構造の正面図をＣＴ画像と一致する。
【０１３８】
たとえば、全部でＮ回のＣＴ結像を行ったとして、ｎ回目のＣＴ結像の始まりおける支持構造の角度がαnであり、終わりにおける支持構造の角度がβnであり、平面図において支持構造が逆時計方向に回転する。ｎ層目のＣＴ画像の中の物体位置を一層目での位置に一致させるため、ｎ層目のＣＴ画像を逆時計方向に角度αn−α1だけ回転すればよい。画像を回転した後、支持構造の正面図にある物体がＣＴ画像と一致するため、支持構造を逆時計方向に角度３６０−（βn−α1）だけ回転すればよい。
【０１３９】
図１７は複数の被検査物体を検査するフローチャートである。図１７に示すように、ステップＳ７１１で、検査員がシステムに電源を入れて、検査インターフェースに登録する。システムは自己検査を行った後、準備完了の状態に入る。そして、検査員は複数の被検査物、たとえば、物体Ａと物体Ｂを直接に支持構造４０の上に載せて、検査ボタンを押す。ここで、物体Ａは支持構造の右上の位置にあり、物体Ｂは支持構造の左下の位置にあるとする。
【０１４０】
次に、ステップＳ７１２で、検査員はスタートボタンを押して、上述のような方法でＤＲ走査を行い、ＤＲ画像を生成する。
【０１４１】
ＤＲ画像を獲得したら、ステップＳ７１３Ａで、ＤＲ画像を分析し、ＣＴ結像の走査位置を自動的に確定する。あるいは、ステップＳ７１３Ｂで、検査員がマウスなどの入力装置６５でＣＴ走査の位置を指定する。これによって、液体物品における代表的な位置だけをＣＴ走査することによって、検査速度を向上できるし、検査品質を下げることはない。
【０１４２】
ＣＴ走査の位置を確定した後、ステップＳ７１４でＣＴ走査過程に入る。即ち、前記確定した走査位置で液体物品に対してＣＴ走査を行って、ＣＴ走査データを獲得し、また再構成算法によってＣＴ画像を再構成する。このＣＴ画像の各画素が液体物品の中の対応する部分の密度と原子番号を表している。多層の液体が存在する場合、各層に対してＣＴ走査を行う。
【０１４３】
ステップＳ７１５とＳ７１６で、最終回のＣＴ結像が終わると、物体Ａと物体Ｂを区別するため、各層のＣＴ画像における各物体の位置が一致する（たとえば、一層目のＣＴ画像に合わせる）ように、同時に支持構造の上にある物体の実際の位置と一致させる（たとえば、平面図と一致させる）ように、上述の方法に従って、支持構造の角度と各層のＣＴ画像を調整する。
【０１４４】
ステップＳ７１７で、各層のＣＴ画像に対して（たとえば、分水嶺算法によって）画像分割処理を行い、各被検査物の液体区域を得る。ステップＳ７１８で、各液体区域内の各画素の密度と原子番号の平均値を統計する。ステップＳ７１９で、各液体区域内で測定した密度と原子番号で確定した点が密度−原子番号の座標空間における予定区域にあるかどうかを判断する。それによって、ステップＳ７２０で、当被検査物の当層の液体が危険物であるかどうかを判断する。
【０１４５】
ステップＳ７２１で、各層の検査結果をまとめて検査員に提示する。まとめの方法として、すべてのＣＴ画像のすべての液体区域が安全な危険品である場合、今回の検査の結論を「安全」とする。そうでなければ、今回の検査の結論を「疑わしい」とする。また、各層のＣＴ画像をカラーにしてユーザに表示する。ＣＴ画像において、疑わしい物品と判断する被検査物を一種類の色（たとえば、赤）で表示し、安全液体と判断するのをほかの色（たとえば、緑）で表示する。
【０１４６】
上述したのは二件の液体物品を同時に検査する状況である。もっと多くの液体物品を同時に検査する場合、各液体物品の位置をはっきり検査員に見せるため、図１８に示すように、液体物品を支持する支持構造の表面を、たとえば、図１８のＡ区、Ｂ区、Ｃ区とＤ区のように、複数の区域に区画する。これによって、検査員は各液体物品を対応する区域に入れて、各区域に各自の液体標識情報を入力する。
【０１４７】
［第３実施形態］
計測効率を向上するため、また、細長い被検査液体物品を直接に支持構造に乗せたら、計測において不安定になって計測効果に影響を及ぼす問題を避けるため、第３実施形態では、複数の物体を格子分けのバレル（容器）に入れる方法を提出した。
【０１４８】
従って、第３実施形態と第２実施形態の区別は、複数の物体を検査する場合、格子分けのバレルを利用していることにある。次は、図１９〜図２４を参照しながら、本実施形態の計測設備の具体的な操作過程を説明する。
【０１４９】
図１９は本実施形態における格子分けのバレルの側面図である。図２０は格子分けのバレルの平面図である。図２１は格子分けのバレルの底面図である。
【０１５０】
図１９に示すように、格子分けのバレルは、底部と、底部に結合する側壁を含める。底部の下表面には、複数の角錐形またはほかの形の突出物を設置する。この三つの角錐形の突出物は支持構造４０の対応する位置決め穴に挿入する。これによって、格子分けのバレル７０は支持構造４０と共に回転し、二者の間に相対的移動を防げる。
【０１５１】
さらに、図１９に示すように、検査員が掴んだり搬送しやすいように、側壁の上部の一周に沿って突出エッジを設置する。側壁は、円筒形または角錐形であり、弾性的な材料、たとえば、ポリエチレン（ＰＥ）やアルミニウムで作られる。
【０１５２】
図２０は三種類の格子分けのバレルの平面図である。図に示すように、底部と側壁で形成される空間の中に一つや複数の格子分け部品を設置し、当空間を液体物品を入れるための複数の空間に区画する。これで、複数の液体物品を同時に計測したい場合、これらの複数の液体物品を格子分け部品によって区画された空間に入れ、これらの液体物品を同時に計測する。この場合、側壁の外側表面にマークを設置し、計測の時、バレル内の物品の位置を決める。たとえば、四件の物品の中に一つ疑わしいものがあれば、側壁上のマークに対するこの物品の相対位置に基づいて、どの物品が疑わしいかをユーザに報告する。
【０１５３】
たとえば、図２０の左側図に示すように、格子分け部品の上側に予定寸法の円形断面のマークを設置している。あるいは、ある格子分け部品をほかの部品より短くすることによって、液体物品の位置決めマークとする。
【０１５４】
図２１は格子分けのバレルの底面図である。底部の下表面には三つの突出が均一的に分布している。ただし、この三つの突出の分布が均一でなくでもよい。
【０１５５】
図２２は検査過程における格子分けのバレルとマークを自動計測する過程を示した図である。上述した通り、格子分けのバレルは、特定の構造を有するので、画像処理でよく使われるパターンマッチング方法によって計測される。十字形の格子分けの第一種類のバレルを例とする。ステップＳ８１０で、まず十字形のパターン画像を作り、このパターンの中心を識別の対象となるＣＴ画像の中心に合わせて、マッチング値を得る。
【０１５６】
そして、ステップＳ８１１で、当マッチング値が予定の閾値より大きいかを判断する。もしそうでなければ、パターン画像とＣＴ画像のマッチング値が最大になるまで、パターン画像を回転する。マッチング値がある予定の閾値より大きくなれば、ＣＴ画像には格子分けのバレルが存在すると認め、そうでなければ、ＣＴ画像には格子分けのバレルが存在しない。
【０１５７】
ＣＴ画像に格子分けのバレルが存在する場合、さらに、ステップＳ８１２で、位置決めマークの特徴から位置決めマークを計測する。また、十字形の格子分けの第一種類のバレルを例とする。その位置決めマークは一つの格子分け線の端末にあり、この格子分け線をほかの三本の格子分け線より短くする。
【０１５８】
パターンマッチング方法によってＣＴ画像に格子分けのバレルの存在を確認した後、最大のマッチング値においては、パターン画像の十字線がちょうど格子分け線と重なる。四本の格子分け線の長さを比較して、その中の最も長いものを取り出して、位置決めマークを得る。
【０１５９】
多層ＣＴ結像後、まず各層のＣＴ画像について、格子分けのバレルと位置決めマークを計測する。各層のＣＴ画像について、格子分けのバレルを計測しなかった場合、検査員が格子分けのバレルを使っていないと認める。もし少なくでも一層のＣＴ画像について、格子分けのバレルを計測したら、格子分けのバレルを使っていると認める。各層で計測した位置決めマークの位置が一致しない場合、その中の信号強度最大のものを最終の位置決めマークとする。ある層の位置決めマークの信号強度を定義する方法として、この層の四本の分割線の長さの最大値から平均値を減らし、その差が大きいほど、位置決めマークの信号強度が大きい。
【０１６０】
図２３（Ａ）〜（Ｄ）は計測過程におけるバレルを回転することを示した図である。格子分けのバレルの位置決めマークを予定位置に調整すれば、各格子分けを指定位置に調整できる。第一種類のバレルを例とする。ＣＴ画像の中心を座標系の原点として、極座標系を作る。各格子分けの位置が位置決めマークの角度座標で唯一に決められる。図２４に示すように、各ＣＴ画像におけるシステムの最終的に確定された位置決めマークが極座標の角度γに対応し、システムにおける事前に設定した調整目標がθである場合、各層のＣＴ画像および支持構造を角度θ−γ回転すれば、調整を完成できる。
【０１６１】
図２３（Ａ）、（Ｂ）に示すように、格子分けのバレルを一つ角度回転して、事前に要求した位置に達する。
【０１６２】
図２４は第３実施形態による検査過程のフローチャートである。ステップＳ９１１ないしステップＳ９２１が前記第２実施状態のステップＳ７１１ないしステップＳ７２１とまったく同じであるので、ここではステップＳ９２２ないしステップＳ９２６だけを説明する。
【０１６３】
ステップＳ９２２で、ＣＴ画像を分析し、ステップＳ９２３で、格子分けの検査バレルがあるかどうかを判断する。検査バレルがない場合、今回の検査を終了する。
【０１６４】
検査バレルを使用した時、ステップＳ９２４で、検査バレルの位置決めマークを識別して、各格子分けの位置を確定する。そして、ステップＳ９２５で、支持構造の角度および各層のＣＴ画像の角度に対して二回目の調整を行い、支持構造の上の検査バレルにおける各格子分けの位置、および液体の各層のＣＴ画像における各格子分けの位置を共に予定の位置まで調整する。
【０１６５】
ステップＳ９２６で、各層のＣＴ画像における各格子分け部分内に落ちる各液体区域の検査結果をリストアウトし、各格子分けの検査結論をまとめてユーザに提示する。本実施形態による一つのまとめ方として、ある格子分けに対して、すべてのＣＴ画像における本格子内に落ちるすべての液体区域が安全液体と判断する場合のみ、本格子の検査結論が「安全」とする。そうでない場合、本格子の検査結論が「疑わしい」とする。
【０１６６】
［変形例１］
上述した通り、ＤＲ結像を行ってから、ＣＴ結像を行うことを例として説明したが、直接に螺旋ＣＴ結像を行って、液体物品を検査することもできる。
【０１６７】
液体物品に対して螺旋ＣＴ検査を行った結果は一組の螺旋ＣＴ画像であり、この組の螺旋ＣＴ画像を比較して分析して、液体物品の中の液体の位置、および液体の中に分層があるかを確定できる。そして、上述と同様な方法によって各層位置の液体の物理属性値、たとえば、密度と原子番号を獲得する。たとえば、図２５に示すような液体物品に対して螺旋ＣＴ検査を行うと、螺旋ピッチ２ｃｍで走査して、図２６（Ａ）〜（Ｍ）に示すような一組のＣＴ画像を得られる。このように、各螺旋ＣＴ画像の中の画素を分析したら、液体物品の中の液体の位置が分かる。ここでの螺旋ＣＴ結像は大螺旋ピッチＣＴ結像または慣用の螺旋ピッチＣＴ結像を利用することができる。
【０１６８】
本領域の技術者が分かる通り、上述した変形実施例も、第２実施形態や第３実施形態と同様に、複数のボトルを検査できる。
【０１６９】
また、前記では密度と原子番号を例として説明したが、本発明は密度または原子番号のみを使ったり、もっと多くのほかの物理属性を使って、危険物を識別することができる。前の状況では、検査過程において、一種類の物理属性値が予定の範囲内に落ちるかによって液体が危険物であるかどうかを判断する。後の状況では、検査過程において、複数の物理属性値で確定される点がこれらの属性に限られる多元空間内の予定区域に落ちるかによって、液体が危険物であるかどうかを判断する。
【０１７０】
［変形例２］
上述したとおり、ＤＲ結像したから、ダブルエネルギーＣＴ結像を行い、液体の密度と原子番号を獲得するが、ここでのＤＲ結像は必須ではない。異なる液体物品に対して、事前に指定した予定位置でダブルエネルギーＣＴ結像を行い、液体の密度と原子番号を獲得することもできる。
【０１７１】
上述したものが本発明の実施形態にのみ適合するが、本分野の技術者が分かるように、本発明の範囲を離れない前提で、本発明についていかなる修正、局部変更は、本発明の請求範囲内にある。従って、本発明の保護範囲は、特許請求の範囲によって限られるものとすべきである。
【図面の簡単な説明】
【０１７２】
【図１】本発明の実施形態による検査設備の構造を示した図である。
【図２】図１で示した計算機データ処理装置６０を説明するためのブロック図である。
【図３】本発明の第１実施形態によるコントローラの構造を示したブロック図である。
【図４】図４（Ａ）、（Ｂ）は、本発明の液体物品の検査方法を示した図である。
【図５】ＤＲ結像とＣＴ結像の関係を説明するための図である。
【図６】ＤＲ結像の結果の一例を示した図である。
【図７】ＣＴ結像の結果の一例を示した図である。
【図８】本発明の第１実施形態による液体物品の検査方法の全体のフローチャートである。
【図９】ＤＲ結像過程のフローチャートである。
【図１０】ＤＲ結像過程における計測採集装置３０で採集したＤＲ画像データの行列方式を示した図である。
【図１１】ＤＲ画像を処理してＣＴ走査位置を確定するフローチャートである。
【図１２】ＣＴ結像過程を示した図である。
【図１３】ＣＴ結像過程におけるＣＴ投影データの行列方式を示した図である。
【図１４】液体属性の測定過程を示した図である。
【図１５】図１５（Ａ）、（Ｂ）は、本発明の第２実施形態による複数の液体物品を検査する場合、再構成したＣＴ画像を示した図である。
【図１６】図１６（Ａ）〜（Ｋ）は、ＣＴ結像を完成した後、ＣＴ再構成画像および／または支持構造を回転して、それをＣＴ走査前の物体の相対位置と一致させる過程を示した図である。
【図１７】複数の被検査物体を検査するフローチャートである。
【図１８】本発明の第２実施形態による支持構造の平面図である。
【図１９】本実施形態における格子分けのバレルの側面図である。
【図２０】格子分けのバレルの平面図である。
【図２１】格子分けのバレルの底面図である。
【図２２】検査過程における格子分けのバレルとマークを自動計測する過程を示した図である。
【図２３】図２３（Ａ）〜（Ｄ）は、計測過程におけるバレルを回転することを示した図である。
【図２４】第３実施形態による検査過程のフローチャートである。
【図２５】液体物品に対して螺旋ＣＴ走査を行うことを示した図である。
【図２６】図２６（Ａ）〜（Ｍ）は、液体物品に対して螺旋ＣＴ走査を行って得た画像である。
【図１】
0005164777.tif 000018
【図２】
0005164777.tif 000019
【図３】
0005164777.tif 000020
【図４】
0005164777.tif 000021
【図５】
0005164777.tif 000022
【図６】
0005164777.tif 000023
【図７】
0005164777.tif 000024
【図８】
0005164777.tif 000025
【図９】
0005164777.tif 000026
【図１０】
0005164777.tif 000027
【図１１】
0005164777.tif 000028
【図１２】
0005164777.tif 000029
【図１３】
0005164777.tif 000030
【図１４】
0005164777.tif 000031
【図１５】
0005164777.tif 000032
【図１６】
0005164777.tif 000033
【図１７】
0005164777.tif 000034
【図１８】
0005164777.tif 000035
【図１９】
0005164777.tif 000036
【図２０】
0005164777.tif 000037
【図２１】
0005164777.tif 000038
【図２２】
0005164777.tif 000039
【図２３】
0005164777.tif 000040
【図２４】
0005164777.tif 000041
【図２５】
0005164777.tif 000042
【図２６】
0005164777.tif 000043