JP6131606B2 - 放射線撮影装置およびそれにおける画像処理方法 - Google Patents

放射線撮影装置およびそれにおける画像処理方法 Download PDF

Info

Publication number
JP6131606B2
JP6131606B2 JP2013008437A JP2013008437A JP6131606B2 JP 6131606 B2 JP6131606 B2 JP 6131606B2 JP 2013008437 A JP2013008437 A JP 2013008437A JP 2013008437 A JP2013008437 A JP 2013008437A JP 6131606 B2 JP6131606 B2 JP 6131606B2
Authority
JP
Japan
Prior art keywords
parameter
parameters
radiation
stage
calculation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2013008437A
Other languages
English (en)
Other versions
JP2014139533A (ja
Inventor
雄介 田川
雄介 田川
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shimadzu Corp
Original Assignee
Shimadzu Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shimadzu Corp filed Critical Shimadzu Corp
Priority to JP2013008437A priority Critical patent/JP6131606B2/ja
Priority to US14/154,922 priority patent/US9476844B2/en
Priority to CN201410028530.7A priority patent/CN103940835B/zh
Publication of JP2014139533A publication Critical patent/JP2014139533A/ja
Application granted granted Critical
Publication of JP6131606B2 publication Critical patent/JP6131606B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/30—Accessories, mechanical or electrical features
    • G01N2223/33—Accessories, mechanical or electrical features scanning, i.e. relative motion for measurement of successive object-parts
    • G01N2223/3304—Accessories, mechanical or electrical features scanning, i.e. relative motion for measurement of successive object-parts helicoidal scan
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/40—Imaging
    • G01N2223/419—Imaging computed tomograph
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/60—Specific applications or type of materials
    • G01N2223/611—Specific applications or type of materials patterned objects; electronic devices
    • G01N2223/6113—Specific applications or type of materials patterned objects; electronic devices printed circuit board [PCB]

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Pulmonology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Description

この発明は、放射線照射手段から照射されて対象物を透過した放射線を放射線検出手段で検出することにより得られた投影画像に基づいて放射線撮影を行う放射線撮影装置およびそれにおける画像処理方法に関する。
放射線として、Ｘ線を例に採って説明する。なお、対象物としては、実装基板、多層基板のスルーホール／パターン／はんだ接合部、パレット上に配置された集積回路（IC: Integrated Circuit）のような実装前の電子部品、金属などの鋳物、ビデオデッキのような成型品などがある。具体的には、電子部品の検査（例えば、基板の配線検査やBGA(Ball Grid Array)やはんだ接合部、ボイド（空隙）などの検査）やこれら対象物の内部欠陥の検査などに用いられる。
従来知られているＣＴ(Computed Tomography)装置では、図１０に示すように、互いに対向するように配置されたＸ線管Ｔ（放射線照射手段）とＸ線検出器Ｄ（放射線検出手段）との間に対象物Ｏを回転ステージＳ上に載置する。回転ステージＳの面に直交する軸を回転軸Ａｘとして回転ステージＳを回転駆動することにより、回転ステージＳ上に載置された対象物Ｏを回転させる。このように回転軸Ａｘの周りに対象物Ｏを回転させることで様々な角度からの対象物Ｏの投影画像を撮影し、それらの投影画像を再構成することで３次元画像を算出している。
Ball Grid Array(BGA)や配線などの非常に微細な構造を有する対象物に対する断層撮影によりＸ線検査を行う際には、拡大率を大きくして撮影する必要がある。しかし、拡大率を大きくするには、Ｘ線管に代表される放射線源と対象物とを近づけて撮影する必要があるので、対象物が平面に広い形状の場合にはＸ線管と対象物とが互いに干渉してしまう恐れがある。その結果、干渉を避けるために拡大率をあまり上げることができない。
そこで、図１１に示すように、回転軸Ａｘからラミノ角傾いた斜め方向にＸ線管ＴおよびＸ線検出器Ｄを配置して、斜め方向でＸ線管ＴからＸ線を照射することで、対象物Ｏが干渉しないように断層撮影を行うプレナーＣＴ(PCT: Planar Computed Tomography)が知られている。また、ステージに回転機構のないＸ線透視装置でも、図１２に示すように、対象物Ｏを載置したステージＳを回転軸Ａｘに垂直な平面（図１２では水平面）で円軌道を描くように平行移動させ、ステージＳの移動に同期させて同一の回転軸Ａｘの軸心周りにＸ線検出器Ｄを回転駆動することで、ＰＣＴを実現する方法がある（例えば、特許文献１、２参照）。
図１２の撮影方法では、ステージとＸ線検出器とが互いに独立した駆動機構となっているので、理想的な断層撮影の走査軌道を得るために、高精度な位置決め、同期が可能な機構と制御とが必要になり、高価になる。そこで、実際の走査軌道と理想走査軌道とのズレがある場合でも、校正用ファントムでキャリブレーション（校正）することで実際の走査軌道を高精度に算出し、再構成用の補正パラメータとして保持する方法がある（例えば、特許文献３参照）。このときの補正パラメータ算出は、コンピュータビジョンの分野で知られているカメラキャリブレーションと同一の問題設定であり、算出アルゴリズムとしてBundle Adjustmentなどの方法が知られている。
Bundle Adjustment法は、画像から抽出した特徴点から、非線形最適化演算により特徴点の３次元座標および撮影時の幾何学的なモデルのパラメータを算出する手法である。このとき繰り返し演算を行うために、初期値の設定方法や特徴点数やフレーム数や推定パラメータ数などの種々の条件によって計算時間が長くなることがある。そこで、計算時間を減らす方法が提案されている（例えば、特許文献４、５参照）。
特許文献４の特開２００７−４８０６８号公報では、従来のヤコビ行列が“０”を多く含む疎行列であることを利用した演算の高速化が、画像から抽出する特徴量の種類が２つ以上（点パターンと正方形パターンなど）ある場合に未知パラメータが増えることによって、直接に適用することができなくなるという問題に対して、以下のように対処している。すなわち、未知パラメータを３個以上に分割して高速化演算可能な行列の形にそれぞれ落とし込み、段階的に計算していくことで計算時間の短縮を図っている。
特許文献５の特開２００９−１４６２９号公報では、画像から抽出した特徴点のフレーム間での対応付けに誤りがある場合や、不適切な特徴点を抽出してしまった場合に、非線形最適化演算の収束性が悪くなり、パラメータ算出精度および計算時間に悪影響を及ぼすという問題に対して、以下のように対処している。すなわち、評価関数計算時に外れ値を除くようなロバスト推定法（LmedS法、RANSAC法など）を適用し外れ値の影響を抑えることで、非線形最適化演算の高精度化や高速化を図っている。
特開２０１０−２２２１号公報 特許第３６９４８３３号 特許第４４１５７６２号 特開２００７−４８０６８号公報 特開２００９−１４６２９号公報
しかしながら、特許文献４の特開２００７−４８０６８号公報では、特徴量の種類が２つ以上あるという特定の条件での高速化手法であり、画像から特徴点のみを抽出してBundle Adjustment法を実施する場合には無関係である。また、特許文献５の特開２００９−１４６２９号公報では、不適切な特徴点抽出や、フレーム間での特徴点の誤対応を前提としているが、校正用ファントムを用いてコントロールされた環境の下でキャリブレーションを実施する場合にはこのような影響（パラメータ算出精度および計算時間への悪影響）は少なく、高速化に結びつくとは言い難い。また、従来手法ではパラメータ数が増加するにつれて非線形最適化演算での計算時間が長くなるという問題に対する高速化手法については提案されていない。
この発明は、このような事情に鑑みてなされたものであって、非線形最適化演算の収束精度を確保しつつ、計算を高速化することができる放射線撮影装置およびそれにおける画像処理方法を提供することを目的とする。
この発明は、このような目的を達成するために、次のような構成をとる。
すなわち、この発明の放射線撮影装置は、対象物を載置するステージと、そのステージを間に挟んで互いに対向するように配置された放射線照射手段および放射線検出手段とを備え、前記放射線照射手段から照射されて前記対象物を透過した放射線を前記放射線検出手段で検出することにより得られた投影画像に基づいて放射線撮影を行う放射線撮影装置であって、校正用ファントムもしくは前記対象物を前記ステージに載置した状態で、前記放射線照射手段、前記ステージの少なくともいずれかを駆動する駆動手段と、複数の前記校正用ファントムの投影画像に基づいて、放射線照射手段、ステージ、放射線検出手段の幾何関係を表すパラメータを非線形最適化演算により算出するパラメータ算出手段と、複数の前記対象物の投影画像および前記パラメータに基づいて、対象物の断層画像を算出し、当該算出時に前記対象物を撮影したときの撮影条件に適した前記パラメータを利用して再構成する断層画像算出手段とを備え、前記パラメータ算出時に、前記校正用ファントムを撮影したときの撮影条件に基づいて、非線形最適化演算によって算出するパラメータの数を減らして限定するパラメータ数限定手段と、過去に算出されたパラメータを更新するパラメータ更新手段とを備え、当該パラメータ更新手段は、前記パラメータ数限定手段で限定されたパラメータのみを更新することを特徴とするものである。
［作用・効果］この発明に係る放射線撮影装置によれば、パラメータ算出手段は、複数の校正用ファントムの投影画像に基づいて、放射線照射手段、ステージ、放射線検出手段の幾何関係を表すパラメータを非線形最適化演算により算出し、複数の対象物の投影画像および上述のパラメータ（補正パラメータ）に基づいて、断層画像算出手段は、対象物の断層画像を算出し、当該算出時に対象物を撮影したときの撮影条件に適したパラメータを利用して再構成する。このパラメータ算出時に、校正用ファントムを撮影したときの撮影条件に基づいて、パラメータ数限定手段は、非線形最適化演算によって算出するパラメータの数を減らして限定する。このように既知の撮影条件に基づいて、放射線照射手段、ステージ、放射線検出手段の幾何関係を表すパラメータの初期値を推定し、装置の機構や駆動特性から考えて、パラメータの初期値と実際に撮影したときのパラメータとの誤差が大きいと考えられるパラメータに限定して、非線形最適化演算を行う。その結果、撮影条件という断層撮影の情報を利用することで、非線形最適化演算の収束精度を確保しつつ、計算を高速化することができる。さらに、パラメータの更新時には、限定されたパラメータのみを更新しているので、パラメータ更新手段で更新するパラメータの数が少なくて済み、非線形最適化演算の収束精度を確保しつつ、計算を高速化することができる。
上述したこの発明に係る放射線撮影装置において、上述したパラメータ数限定手段は、撮影条件を変更したときに変更された当該撮影条件に基づいて、非線形最適化演算によって算出するパラメータの数を減らして限定するのが好ましい。例えば、装置の出荷後や据え付け後に一度だけ非線形最適化演算を行い、撮影条件を変更したときに変更された当該撮影条件に基づいてパラメータの数を減らして限定する。駆動系の駆動誤差だけでなく、放射線撮影装置の各構成要素を出荷して据え付けるときの据え付け誤差も存在し、パラメータの初期値との誤差となる。そこで、装置の出荷後や据え付け後に一度だけ非線形最適化演算を行うことで据え付け誤差を含んだパラメータを算出し、撮影条件を変更したときには駆動系の駆動誤差が関係する部分のパラメータのみを推定する。
なお、パラメータの数を減らして限定するタイミングは、撮影条件を実際に変更したときのみに限定されない。例えば、装置の電源投入時に駆動系の駆動誤差が生じる可能性があるので、装置の電源投入後の度に当該電源投入時の撮影条件に基づいてパラメータの数を減らして限定してもよい。なお、パラメータの数を減らして限定する設定については、撮影条件が変更されたのを中央演算処理装置（ＣＰＵ）などが判断して自動に行ってもよいし、オペレータが手動で行ってもよいし、自動と手動とを組み合わせて行ってもよい。
上述したこれらの発明に係る放射線撮影装置において、上述したパラメータ数限定手段は、非線形最適化演算によって算出するパラメータを、放射線照射手段とステージとの座標位置関係を表す３つのパラメータに限定するのが好ましい。例えば、３軸駆動において、拡大率を変更した場合や、放射線照射手段から放射線検出手段へ垂線を下ろした距離（SID: Source Image Distance）を変更する場合には、３つのパラメータに限定するのが可能となる。なお、放射線撮影は透視投影であるので、放射線検出手段では必ず拡大された像が投影される。したがって、放射線検出手段の駆動誤差が投影画像に及ぼす影響よりも、放射線照射手段やステージの駆動誤差が投影画像に及ぼす影響の方が大きい。そこで、影響の大きい放射線照射手段とステージとの座標位置関係を表す３つのパラメータに限定して算出することで、非線形最適化演算の収束精度を確保しつつ、演算の高速化が可能となる。
特に、駆動手段が、ステージの面に直交する軸を回転軸としてステージを回転駆動し、回転軸からラミノ角傾いた斜め方向で放射線照射手段から放射線を照射するプレナーＣＴ(PCT: Planar Computed Tomography)を行う場合に適用することができる。プレナーＣＴ（ＰＣＴ）では通常のＣＴよりも高拡大率での放射線撮影が可能となる。高拡大率にすればするほど、放射線照射手段やステージの駆動誤差が投影画像に及ぼす影響が大きくなるので、３つのパラメータに限定することによる評価関数の収束精度への悪影響は小さくなる。
また、この発明の画像処理方法は、対象物を載置するステージと、そのステージを間に挟んで互いに対向するように配置された放射線照射手段および放射線検出手段とを備え、前記放射線照射手段から照射されて前記対象物を透過した放射線を前記放射線検出手段で検出することにより得られた投影画像に基づいて放射線撮影を行う放射線撮影装置における画像処理方法であって、複数の校正用ファントムの投影画像に基づいて、放射線照射手段、ステージ、放射線検出手段の幾何関係を表すパラメータを非線形最適化演算により算出する際に、前記校正用ファントムを撮影したときの撮影条件に基づいて、非線形最適化演算によって算出するパラメータの数を減らして限定するパラメータ数限定過程と、非線形最適化演算によって前記パラメータを算出するパラメータ算出過程とを備え、前記パラメータ算出過程は、過去に算出されたパラメータを更新するパラメータ更新過程を備え、当該パラメータ更新では、前記パラメータ数限定過程で限定されたパラメータのみを更新することを特徴とするものである。
［作用・効果］この発明に係る画像処理方法によれば、パラメータ算出時に、校正用ファントムを撮影したときの撮影条件に基づいて、非線形最適化演算によって算出するパラメータの数を減らして限定する。さらには、非線形最適化演算によって上述のパラメータを算出するパラメータ算出過程を備え、パラメータ算出過程は、過去に算出されたパラメータを更新するパラメータ更新過程を備え、当該パラメータ更新では、上述のパラメータ数限定過程で限定されたパラメータのみを更新する。このように既知の撮影条件に基づいて、放射線照射手段、ステージ、放射線検出手段の幾何関係を表すパラメータの初期値を推定し、装置の機構や駆動特性から考えて、パラメータの初期値と実際に撮影したときのパラメータとの誤差が大きいと考えられるパラメータに限定して、非線形最適化演算を行う。その結果、撮影条件という断層撮影の情報を利用することで、非線形最適化演算の収束精度を確保しつつ、計算を高速化することができる。さらに、パラメータの更新時には、限定されたパラメータのみを更新しているので、パラメータ更新過程で更新するパラメータの数が少なくて済み、非線形最適化演算の収束精度を確保しつつ、計算を高速化することができる。
具体的には、上述のパラメータ算出過程は、特徴点３次元座標推定過程，再投影座標計算過程，再投影二乗誤差計算過程および収束判定過程を、さらに備え、上述のパラメータ更新過程では、上述のパラメータ数限定過程で限定されたパラメータのみ更新する。特徴点３次元座標推定過程では、校正用ファントムの投影画像から抽出された特徴点、およびパラメータの初期値に基づいて、特徴点の３次元座標を推定する。再投影座標計算過程では、特徴点３次元座標推定過程で推定した特徴点３次元座標に基づいて、再投影座標を計算する。再投影二乗誤差計算過程では、再投影座標計算過程で算出した再投影座標、および上述した特徴点に基づいて、再投影二乗誤差を算出する。収束判定過程では、再投影二乗誤差計算過程で算出した再投影二乗誤差の収束判定を行う。そして、上述のパラメータ更新過程では、再投影二乗誤差の値が収束していない場合には、再投影二乗誤差の値が小さくなるようにパラメータを更新するが、その更新の際には上述のパラメータ数限定過程で限定されたパラメータのみ更新する。このようにパラメータの更新時には、算出すべきパラメータの数を撮影条件に基づいて予め限定しているので、パラメータ更新過程で更新するパラメータの数が少なくて済み、非線形最適化演算の収束精度を確保しつつ、計算を高速化することができる。
この発明に係る放射線撮影装置およびそれにおける画像処理方法によれば、パラメータ算出時に、校正用ファントムを撮影したときの撮影条件に基づいて、パラメータ数限定手段は、非線形最適化演算によって算出するパラメータの数を減らして限定する。このように既知の撮影条件に基づいて、放射線照射手段、ステージ、放射線検出手段の幾何関係を表すパラメータの初期値を推定し、装置の機構や駆動特性から考えて、パラメータの初期値と実際に撮影したときのパラメータとの誤差が大きいと考えられるパラメータに限定して、非線形最適化演算を行う。その結果、撮影条件という断層撮影の情報を利用することで、非線形最適化演算の収束精度を確保しつつ、計算を高速化することができる。さらに、パラメータの更新時には、限定されたパラメータのみを更新しているので、更新するパラメータの数が少なくて済み、非線形最適化演算の収束精度を確保しつつ、計算を高速化することができる。
実施例に係るＸ線検査装置の概略構成図である。 実施例に係るＸ線検査装置の概略斜視図である。 校正用ファントムによるキャリブレーション法を併記した実施例に係るＸ線検査装置の概略斜視図である。 各座標系を併記した実施例に係るＸ線検査装置の概略斜視図である。 実施例に係るＸ線検査装置のブロック図である。 放射線撮影時に校正用ファントムを搭載可能にステージを構成する一例の概略図であり、（ａ）はステージにマーカを埋め込んで校正用ファントムとする場合の概略斜視図、（ｂ）はその概略断面図である。 放射線撮影時に校正用ファントムを搭載可能にステージを構成する他の一例の概略図であり、（ａ）はステージの表面（上面）に、パターンをマーカとして設置したときの校正用ファントムの概略斜視図、（ｂ）はステージの裏面（下面）に、パターンをマーカとして設置したときの概略斜視図である。 校正用ファントムの投影画像からパラメータを算出するまでの一連の処理を示すフローチャートである。 図８のパラメータ算出の具体的なフローチャートである。 従来の撮影の概略図である。 従来の斜め撮影の概略図である。 ステージを平行移動させ、ステージの移動に同期させてＸ線検出器を回転駆動させたときの従来の斜め撮影の概略図である。
以下、図面を参照してこの発明の実施例を説明する。
図１は、実施例に係るＸ線検査装置の概略構成図であり、図２は、実施例に係るＸ線検査装置の概略斜視図であり、図３は、校正用ファントムによるキャリブレーション法を併記した実施例に係るＸ線検査装置の概略斜視図であり、図４は、各座標系を併記した実施例に係るＸ線検査装置の概略斜視図であり、図５は、実施例に係るＸ線検査装置のブロック図である。本実施例では、放射線撮影装置として、Ｘ線検査装置を例に採って説明する。
図１に示すように、Ｘ線検査装置１は、対象物Ｏを載置するステージ２と、そのステージ２を間に挟んで互いに対向するように配置されたＸ線管３およびＸ線検出器４とを備えている。Ｘ線検出器４については、イメージインテンシファイア（Ｉ．Ｉ）やフラットパネル型Ｘ線検出器（FPD: Flat Panel Detector）などに例示されるように、特に限定されない。本実施例では、Ｘ線検出器４としてフラットパネル型Ｘ線検出器（ＦＰＤ）を例に採って説明する。ステージ２は、この発明におけるステージに相当し、Ｘ線管３は、この発明における放射線照射手段に相当し、Ｘ線検出器４は、この発明における放射線検出手段に相当する。
ＦＰＤは、画素に対応して縦横に並べられた複数の検出素子からなり、Ｘ線を検出素子が検出して、検出されたＸ線のデータ（電荷信号）をＸ線検出信号として出力する。このようにして、Ｘ線管３から照射されて対象物Ｏを透過したＸ線をＦＰＤからなるＸ線検出器４が検出してＸ線検出信号を出力し、Ｘ線検出信号に基づく画素値を画素に対応してそれぞれ並べることで、Ｘ線検出器４の検出面に投影された投影画像を取得する。
その他に、Ｘ線検査装置１は、図１に示すように、Ｘ線検出器４を矢印Ｒ１周りに回転駆動する検出器回転機構５と、Ｘ線検出器４を矢印Ｒ２方向に傾動させる検出器傾動機構６とを備えている。検出器傾動機構６は、Ｘ線検出器４を支持する円弧状のガイド部６ａ、および回転モータ（図示省略）からなり、回転モータが回転駆動することで、ガイド部６ａに沿ってＸ線検出器４が矢印Ｒ２方向に傾動する。
検出器回転機構５は、回転モータ（図示省略）からなり、回転モータが検出器傾動機構６のガイド部６ａを矢印Ｒ１周りに回転駆動することで、ガイド部６ａに支持されたＸ線検出器４も矢印Ｒ１周りに回転駆動する。また、本実施例では、検出器回転機構５は、ステージ２の駆動に同期させてＸ線検出器４を矢印Ｒ１周りに回転駆動する。特に、Ｘ線管３から照射されたＸ線が対象物Ｏの注目点を透過してＸ線検出器４の中心部分で検出されるように、検出器回転機構５はＸ線検出器４を矢印Ｒ１周りに回転駆動する。
本実施例では、Ｘ線管３は装置筐体（図示省略）に固定されており、Ｘ線管３から広範囲に照射されたＸ線が対象物Ｏの注目点を通り、Ｘ線検出器４で検出される。
その他に、Ｘ線検査装置１は、図５に示すように、ステージ２を直交座標系であるワールド座標系ＸＷ，ＹＷ，ＺＷ（図１、図４を参照）でそれぞれ矢印Ｒ３，Ｒ４，Ｒ５（図１、図２を参照）方向に直進駆動するステージ駆動機構７と、校正用ファントムＰｈ（図３、図５〜図７を参照）の投影画像に基づいて補正用のパラメータを算出するパラメータ算出部８と、パラメータ算出部８で算出されたパラメータを書き込んで記憶するパラメータ記憶部９と、複数の対象物Ｏの投影画像およびパラメータに基づいて、対象物Ｏの断層画像を算出して演算する断層画像算出部１０と、これらを統括制御するコントローラ１１と、断層画像算出部１０で得られた断層画像を出力（モニタに表示出力あるいはプリンタに印刷出力）する画像出力部１２とを備えている。上述の検出器回転機構５や検出器傾動機構６を含め、ステージ駆動機構７は、この発明における駆動手段に相当し、パラメータ算出部８は、この発明におけるパラメータ算出手段に相当し、パラメータ算出部８は、この発明におけるパラメータ数限定手段にも相当し、断層画像算出部１０は、この発明における断層画像算出手段に相当する。
ステージ駆動機構７は、ステージ２をＸＷ（図１、図４を参照）に平行な矢印Ｒ３（図１、図２を参照）方向に直進駆動（ここでは水平駆動）するＸ軸直進モータ（図示省略）、ステージ２をＹＷ（図１、図４を参照）に平行な矢印Ｒ４（図２を参照）方向に直進駆動（ここでは水平駆動）するＹ軸直進モータ（図示省略）、およびステージ２をＺＷ（図１、図４を参照）に平行な矢印Ｒ５（図１、図２を参照）方向に直進駆動（ここでは昇降駆動）するＺ軸直進モータ（図示省略）からなる。本実施例では、各々のＸ軸直進モータ，Ｙ軸直進モータによる軌道の合成が円軌道で、かつ検出器回転機構５によるＸ線検出器４の回転駆動に同期するようにコントローラ１１はステージ２を駆動する。
このようにステージ２を円軌道上に駆動しつつ、Ｘ線検出器４を回転駆動することで、ステージ２の向きを一定に固定することができる。また、図３に示すように、校正用ファントムＰｈをステージ２に載置し、実施する断層撮影条件に基づいて、実線で示されている回転軸Ａｘを中心とした円軌道のようにステージ２およびＸ線検出器４を同期して駆動させて放射線撮影を行うことで、複数の校正用ファントムＰｈの投影画像を取得する。
なお、本実施例では、図３に示すような円軌道（従来の図１２に示す円軌道と同じ）で動作させて断層撮影を行っているが、従来の図１０〜図１２のような断層撮影が可能な装置構成であればよく、図１〜図４以外の軌道でも様々な装置構成を採用することが可能である。また、ステージ２とＸ線管３とＸ線検出器４とが相対的に同じ軌道で動作すればよいので、例えば図１０の場合では、ステージを回転駆動する替りにＸ線管およびＸ線検出器を駆動してもよい。また、図１０〜図１２に限らず直線断層撮影に例示されるように断層撮影であれば特に限定されない。このように、Ｘ線管に代表される放射線照射手段、ステージの少なくともいずれかを駆動する態様であれば、放射線照射手段およびステージの両方を駆動してもよいし、放射線照射手段のみを駆動してもよいし、ステージのみを駆動してもよい。
なお、図４に示すように、ＷＣをワールド座標系とし、ＣＣをカメラ座標系とし、ＩＣを画像面座標系とし、ＰＣをピクセル座標系とする。ワールド座標系は任意に設定した基準座標系であり、原点としては例えば対象物Ｏ（図１、図２を参照）や校正用ファントムＰｈ（図３、図５〜図７を参照）の位置に設定してもよい。カメラ座標系はＸ線管３（いわゆる光源）を原点として、ＺＣ軸をＸ線管３からの照射軸（いわゆる光軸）方向とし、ＸＣ軸およびＹＣ軸をＸ線検出器４の検出面の水平方向・垂直方向（画像面座標系のｘ軸・ｙ軸に平行）に合わせた座標系である。画像面座標系は光軸と検出面とが交わる位置を原点として、ｘ軸・ｙ軸を検出面の水平方向・垂直方向に合わせた座標系である。さらに、ピクセル座標系は画像の左上を原点とした画素単位（ピクセル単位）の座標系である。
図５の説明に戻り、パラメータ算出部８は、校正用ファントムＰｈ（図３、図６および図７も参照）の投影画像に基づいて補正用のパラメータを算出し、算出されたパラメータを撮影条件毎にパラメータ記憶部９に書き込んで記憶する。対象物Ｏの断層画像を補正する際に、パラメータ記憶部９に記憶されたパラメータを読み出して、その読み出されたパラメータを利用する。パラメータ記憶部９は、ＲＡＭ（Random Access Memory）などに代表される記憶媒体で構成されている。
さらに、パラメータの初期値をパラメータ記憶部９に書き込んで記憶し、パラメータ記憶部９に記憶されたパラメータを読み出して、パラメータ算出部８はパラメータを更新して、更新されたパラメータをパラメータ記憶部９に書き込んで記憶する。パラメータの更新の際には限定されたパラメータのみ更新する。よって、パラメータ算出部８は、この発明におけるパラメータ数限定手段の機能を有する。パラメータ数限定手段の具体的な機能については後述する。
断層画像算出部１０は、複数の対象物Ｏの投影画像および上述のパラメータに基づいて、対象物Ｏの断層画像を算出し、当該算出時に対象物Ｏを撮影したときの撮影条件に適したパラメータを利用して再構成する。コントローラ１１は、Ｘ線検査装置１を構成する各部分を統括制御し、特に、検出器回転機構５の回転モータ（図示省略）、検出器傾動機構６の回転モータ（図示省略）、ステージ駆動機構７のＸ軸／Ｙ軸／Ｚ軸直進モータ（図示省略）をそれぞれ制御する。上述のパラメータ算出部８や断層画像算出部１０やコントローラ１１は、中央演算処理装置（ＣＰＵ）などで構成されている。
図１に示すように、Ｘ線管３、対象物ＯおよびＸ線検出器４を配置することで、図１２と同様に、ラミノ角傾いた斜め方向にＸ線管３とＸ線検出器４とを配置して斜め方向から撮影（プレナーＣＴ）することができる。そして、Ｘ線管３と、ステージ２ひいては対象物Ｏとを近づけることができ、Ｘ線管３と対象物Ｏとが互いに干渉することなく高拡大率の投影画像を得ることができる。ステージ２を駆動する度に投影画像を取得することで複数の角度からの投影画像を取得し、図５に示す断層画像算出部１０は、複数の投影画像に基づいて断層画像を算出して演算する。
次に、校正用ファントムＰｈについて、図６および図７を参照して説明する。図６は、放射線撮影時に校正用ファントムを搭載可能にステージを構成する一例の概略図であり、図６（ａ）は、ステージにマーカを埋め込んで校正用ファントムとする場合の概略斜視図であり、図６（ｂ）は、その概略断面図であり、図７は、放射線撮影時に校正用ファントムを搭載可能にステージを構成する他の一例の概略図であり、図７（ａ）は、ステージの表面（上面）に、パターンをマーカとして設置したときの校正用ファントムの概略斜視図であり、図７（ｂ）は、ステージの裏面（下面）に、パターンをマーカとして設置したときの校正用ファントムの概略斜視図である。
図６（ａ）の概略斜視図や図６（ｂ）の概略断面図に示すように、ステージ２に校正用ファントムＰｈを埋め込んでもよいし、図７の概略斜視図に示すように、ステージ２の表面（上面）（図７（ａ）を参照）あるいは裏面（下面）（図７（ｂ）を参照）にパターンからなるマーカを設置してもよい。また、図７（ａ）および図７（ｂ）を組み合わせて、ステージ２の表裏面ともにマーカを設置してもよいし、その他の面（側面、正面あるいは背面）にマーカを設置してもよい。
一般に、ステージ２は放射線透過率が高い物質で形成されているので、マーカを放射線透過率が低い物質（例えば鉛）で形成して、ステージ２と区別する。図６の場合には、マーカＭを放射線透過率が低い球状物質で形成し、その球状物質をステージ２に埋め込んだものを校正用ファントムＰｈとする。
また、図７の場合には、マーカを設置することで、マーカからなる校正用ファントムＰｈを搭載可能にステージ２を構成する。図７のマーカも、パターンの交点を特徴点とすることができるので、マーカの概念に包括されることに留意されたい。
また、対象物Ｏが校正用ファントムＰｈを兼用してもよい。例えば対象物ＯがＢＧＡの場合には、対象物Ｏ中の球状物質（Ｂａｌｌ）のいずれかにマーカを設置する。マーカを設置することで、マーカからなる校正用ファントムＰｈを対象物Ｏが兼用することができ、対象物Ｏとは別の校正用ファントムを用意する必要がなくなる。
以上のような校正用ファントムＰｈを予め搭載可能にステージ２を構成する構造に限定されない。図４に示すように、放射線撮影時（必要なとき）のみにステージ２に校正用ファントムＰｈを載置してもよい。
いずれの場合においても、校正用ファントムＰｈは、投影画像から画像処理により特徴点が抽出可能であり、さらに各投影画像間で特徴点の対応付けが可能であるようなパターンを有している必要がある。特徴点としては、例えば図６に示すような球形状の中心点や、図７に示すような格子形状の交点などが考えられる。また、非線形最適化演算（非線形的な最適化方法）によりパラメータを算出する場合には、「特徴点数×２×フレーム数＞特徴点数×３＋パラメータ数×フレーム数」が成り立つ数だけ、特徴点が３次元的に配置されている必要がある。
次に、校正用ファントムによるキャリブレーション（校正）について、図８および図９を参照して説明する。図８は、校正用ファントムの投影画像からパラメータを算出するまでの一連の処理を示すフローチャートであり、図９は、図８のパラメータ算出の具体的なフローチャートである。
（ステップＳ１）断層撮影
図３に示すように、校正用ファントムＰｈをステージ２上に載置し、放射線撮影を実施し、複数の校正用ファントムＰｈの断層画像を取得する。このときの走査軌道は、対象物Ｏ（図１、図２を参照）を断層撮影するときと同じ走査軌道で駆動すればよい。撮影フレーム数については、対象物Ｏを断層撮影するときのフレーム数以上取得するのが好ましい。
（ステップＳ２）特徴点算出
ステップＳ１で取得した校正用ファントムＰｈの投影画像から特徴点を抽出し、特徴点座標を算出する。特徴点抽出の方法については様々あり、特に限定されない。例えば、図６に示すように球形状のマーカを有する校正用ファントムＰｈを撮影した場合には、球の中心を特徴点として抽出することが考えられる。このときの球の中心としては、オープニング処理やエロージョン（erosion）処理などに代表されるモルフォロジー（morphology）演算を利用して球形状の中心位置を算出してもよいし、画素値（輝度値）の情報を利用して、例えば距離変換などにより距離画像の極大値座標を中心位置として算出してもよい。
（ステップＳ３）同一特徴点識別
ステップＳ２で算出された特徴点のフレーム間での対応付けを行う。対応付けの方法としては、時間的に連続するフレームでは対応する特徴点は画像上で近い位置に存在するという特徴を利用して、特徴点座標が互いに近い点を対応する特徴点とする方法などが考えられる。
（ステップＳ４）パラメータ算出
ステップＳ３で対応付けされた特徴点座標からパラメータを算出する。このときのパラメータ算出は、上述したように、コンピュータビジョンの分野で知られているカメラキャリブレーションと同一の問題設定であり、算出アルゴリズムとしてBundle Adjustmentなどの方法が知られている。Bundle Adjustment法は、上述したように、非線形最適化演算により特徴点の３次元座標および撮影時の幾何学的なモデルのパラメータを算出する手法である。
具体的には、透視投影モデルを考えた場合に、図４に示すような座標系をそれぞれ定義すると、特徴点に関するワールド座標系からピクセル座標系への変換は、斉次座標で下記（１）式のように表される。
Figure 0006131606
行列Ｍ中の行列Ｒは３行３列の回転行列であり、回転の３自由度（θｘ，θｙ，θｚ）を有する。また、ベクトルｔは３行１列の並進ベクトルであり、並進３自由度（ｔｘ，ｔｙ，ｔｚ）を有する。つまり、行列Ｍは６つのパラメータからなる。
行列Ｐは、Ｘ線管３（図１を参照）、いわゆる光源からＸ線検出器４（図１、図２を参照）へ垂線を下ろした距離（すなわちＳＩＤ）を表すパラメータｆからなる。
行列Ａは、画素単位（ピクセル単位）へ変換するスケール変数ｋｘ，ｋｙ、画像面座標系とピクセル座標系との原点位置関係を表すσｘ，σｙ、ピクセル座標系のｘ軸・ｙ軸がなる角度を表すｓｋｅｗの合計５つのパラメータからなる。
つまり、透視投影モデルでは、ワールド座標系からピクセル座標系への変換が１２個の幾何関係を表すパラメータ（自由度は１１）で表され、Bundle Adjustment法ではこれらのパラメータを算出する。
図９は、図８のパラメータ算出の具体的なフローチャートであり、Bundle Adjustment法によるパラメータ算出処理フローを示す。また、図９のステップＴ１〜Ｔ５は、この発明におけるパラメータ算出過程に相当する。パラメータの初期値については、断層撮影条件に基づいて設定してもよいし、以前のキャリブレーションで算出された、ある適切なパラメータをパラメータの初期値として設定してもよい。
（ステップＴ１）特徴点３次元座標推定
対応付けされた特徴点（図９では「対応付け特徴点座標群」で表記）（ｐｉｍ）１，１…（ｐｉｍ）ｍ，ｎ、パラメータの初期値Ａ＾１Ｐ＾１Ｍ＾１…Ａ＾ｎＰ＾ｎＭ＾ｎから、特徴点の３次元座標を推定する。このステップＴ１では、全てのパラメータの初期値Ａ＾１Ｐ＾１Ｍ＾１…Ａ＾ｎＰ＾ｎＭ＾ｎが入力される。このステップＴ１は、特徴点３次元座標推定過程に相当する。上記（１）式中の変換行列をｚ＾ｊ＝Ａ＾ｊＰ＾ｊＭ＾ｊで書き直して、ｎフレームの場合に拡張すると、下記（２）式のようになる。
Figure 0006131606
上記（２）式より、擬似逆行列を計算することで、特徴点３次元座標の最小二乗解（ｐ＾Ｗ）ｉが下記（３）式により求まる。
Figure 0006131606
（ステップＴ２）再投影座標計算
ステップＴ１で推定した特徴点３次元座標ｐ＾Ｗを上記（２）式に代入することで、再投影座標（ｐ＾ｉｍ）１，１…（ｐ＾ｉｍ）ｍ，ｎが求まる。このステップＴ２は、再投影座標計算過程に相当する。
（ステップＴ３）再投影二乗誤差計算
ステップＴ２で算出した再投影座標（ｐ＾ｉｍ）１，１…（ｐ＾ｉｍ）ｍ，ｎ、対応付けされた特徴点（ｐｉｍ）１，１…（ｐｉｍ）ｍ，ｎから、再投影二乗誤差Ｃを下記（４）式で算出する。このステップＴ３は、再投影二乗誤差計算過程に相当する。
Figure 0006131606
（ステップＴ４）収束判定
ステップＴ３で算出した再投影二乗誤差Ｃの値が収束したか否かを判定する。このステップＴ４は、収束判定過程に相当する。再投影二乗誤差Ｃの値が十分に収束している場合には、そのときのパラメータを、最終的に数が限定されたパラメータであるとして一連のフローを終了する。
（ステップＴ５）パラメータ更新
再投影二乗誤差Ｃの値が十分に収束していない場合には、Ｃの値が小さくなるようにパラメータを更新する。このステップＴ５は、パラメータ更新過程に相当する。パラメータ更新方法は、非線形最適化の分野で一般的に扱われている問題であり、レーベンバーグマーカート法など種々の方法がある。
全てのパラメータが未知であると考えた場合、Bundle Adjustment法によって推定する変数の数は、上述したように（特徴点数×３＋パラメータ数×フレーム数）個である。「発明が解決しようとする課題」の欄でも述べたように、非線形最適化演算では、変数の数が多いほど計算時間が長くなるという問題がある。そこで、断層撮影の情報（撮影条件）を利用することで、収束精度を確保しつつ、推定する変数の数を減らして限定して計算を高速化する。したがって、ステップＴ５においてパラメータを更新する際には、撮影条件に基づいて予め限定されたパラメータのみ更新するので、更新するパラメータの数が少なくて済む。
非線形最適化演算において断層撮影条件から初期値を設定した場合、初期値と真値との誤差は、装置の設計図面と実際の装置との組み立て誤差および駆動系の駆動誤差が合わさったものであるといえる。組み立て誤差に起因して誤差が生じるパラメータについては、一度キャリブレーションを行えば、断層撮影条件を変更したときにも使用することができる。よって、一度だけ非線形最適化演算により全パラメータを推定し、断層撮影条件を変更してキャリブレーションを再び行う場合には、駆動誤差に起因して誤差が生じるパラメータに限定して非線形最適化演算を実施することで、収束精度を確保しつつ、変数の数を減らすことができる。
また、各駆動系で駆動誤差が生じている場合、それらの駆動誤差が投影画像に及ぼす影響が収束性に影響することが上記（４）式からわかる。放射線撮影は透視撮影であるので、拡大された像がＸ線検出器４に投影される。よって、ステージ２の駆動誤差の影響は拡大率の分だけ大きくなってＸ線検出器４上で現れる。その一方でＸ線検出器４の駆動誤差の影響は拡大されないので、収束性を向上させるためには、ステージ２の駆動精度がより重要である。そこで、例えば推定するパラメータを、ステージ２の駆動精度に関係するパラメータθｘ，θｙ，θｚ，ｔｘ，ｔｙ，ｔｚに限定することで、収束精度の悪化を抑えて、変数の数を減らすことができる。
また、ステージ２が直線駆動要素（直進駆動要素）で構成されている場合には、駆動軸が装置に強固に固定されていれば、座標系が傾くことはさほどにないと考えられる。そこで、推定するパラメータを並進パラメータであるｔｘ，ｔｙ，ｔｚの３つに限定することにより、収束精度の悪化を抑えて、変数の数を減らすことができる。
上述したように、高拡大率で放射線撮影を行うほどステージ２の駆動精度が重要となるので、ＰＣＴで高拡大率で断層撮影を行う場合などは、Ｘ線検出器４の駆動精度に対するステージ２の駆動精度の重要性が増す。つまり、パラメータを上述の並進パラメータｔｘ，ｔｙ，ｔｚに限定する場合と限定しない場合との収束性の差がより変わらなくなってくるので、さらに効果的に収束精度の悪化を抑えて、変数の数を減らすことができる。
図８のステップＳ２〜Ｓ４については、パラメータ算出部８（図５を参照）が行う。また、図９のステップＴ１〜Ｔ５についても、パラメータ算出部８が行う。パラメータ算出部８がパラメータを更新する際には、限定されたパラメータのみ更新する。また、パラメータの数を減らして限定する過程は、この発明におけるパラメータ数限定過程に相当する。
断層画像算出部１０（図５を参照）は対象物Ｏの投影画像の逆投影を実施して、断層画像を算出する。断層画像算出アルゴリズムは公知の方法があり、例えばＦＢＰ法（Filtered Back Projection法）などが知られている。
上述の構成を備えた本実施例に係るＸ線検査装置およびそれにおける画像処理方法によれば、パラメータ算出部８は、複数の校正用ファントムＰｈの投影画像に基づいて、放射線照射手段（本実施例ではＸ線管３）、ステージ２、放射線検出手段（本実施例ではＸ線検出器４）の幾何関係を表すパラメータを非線形最適化演算により算出し、複数の対象物Ｏ投影画像および上述のパラメータ（補正パラメータ）に基づいて、断層画像算出部１０は、対象物Ｏの断層画像を算出し、当該算出時に対象物Ｏを撮影したときの撮影条件に適したパラメータを利用して再構成する。このパラメータ算出時に、校正用ファントムＰｈを撮影したときの撮影条件に基づいて、パラメータ数限定手段（本実施例ではパラメータ算出部８）は、非線形最適化演算によって算出するパラメータの数を減らして限定する。このように既知の撮影条件に基づいて、放射線照射手段（Ｘ線管３）、ステージ２、放射線検出手段（Ｘ線検出器４）の幾何関係を表すパラメータの初期値を推定し、装置の機構や駆動特性から考えて、パラメータの初期値と実際に撮影したときのパラメータとの誤差が大きいと考えられるパラメータに限定して、非線形最適化演算を行う。その結果、撮影条件という断層撮影の情報を利用することで、非線形最適化演算の収束精度を確保しつつ、計算を高速化することができる。
好ましくは、本実施例において、上述したパラメータ数限定手段（パラメータ算出部８）は、撮影条件を変更したときに変更された当該撮影条件に基づいて、非線形最適化演算によって算出するパラメータの数を減らして限定している。例えば、装置の出荷後や据え付け後に一度だけ非線形最適化演算を行い、撮影条件を変更したときに変更された当該撮影条件に基づいてパラメータの数を減らして限定する。駆動系の駆動誤差だけでなく、放射線撮影装置（本実施例ではＸ線検査装置１）の各構成要素を出荷して据え付けるときの据え付け誤差も存在し、パラメータの初期値との誤差となる。そこで、装置の出荷後や据え付け後に一度だけ非線形最適化演算を行うことで据え付け誤差を含んだパラメータを算出し、撮影条件を変更したときには駆動系の駆動誤差が関係する部分のパラメータのみを推定する。
「課題を解決するための手段」の欄でも述べたように、パラメータの数を減らして限定するタイミングは、撮影条件を実際に変更したときのみに限定されない。例えば、装置の電源投入時に駆動系の駆動誤差が生じる可能性があるので、装置の電源投入後の度に当該電源投入時の撮影条件に基づいてパラメータの数を減らして限定してもよい。なお、パラメータの数を減らして限定する設定については、撮影条件が変更されたのを中央演算処理装置（ＣＰＵ）など（本実施例ではコントローラ１１）が判断して自動に行ってもよいし、オペレータが手動で行ってもよいし、自動と手動とを組み合わせて行ってもよい。
好ましくは、本実施例において、上述したパラメータ数限定手段（パラメータ算出部８）は、非線形最適化演算によって算出するパラメータを、放射線照射手段（Ｘ線管３）とステージ２との座標位置関係を表す３つのパラメータ（本実施例では並進パラメータｔｘ，ｔｙ，ｔｚ）に限定している。例えば、３軸駆動（本実施例ではＸ線管３に関しては固定、ステージ２に関しては水平および昇降駆動のみ）において、拡大率を変更した場合や、放射線照射手段（Ｘ線管３）から放射線検出手段（Ｘ線検出器４）へ垂線を下ろした距離（SID: Source Image Distance）を変更する場合には、３つのパラメータ（並進パラメータｔｘ，ｔｙ，ｔｚ）に限定するのが可能となる。なお、放射線撮影は透視投影であるので、放射線検出手段（Ｘ線検出器４）では必ず拡大された像が投影される。したがって、放射線検出手段（Ｘ線検出器４）の駆動誤差が投影画像に及ぼす影響よりも、放射線照射手段（Ｘ線管３）やステージ２の駆動誤差が投影画像に及ぼす影響の方が大きい。そこで、影響の大きい放射線照射手段（Ｘ線管３）やステージ２との座標位置関係を表す３つのパラメータ（並進パラメータｔｘ，ｔｙ，ｔｚ）に限定して算出することで、非線形最適化演算の収束精度を確保しつつ、演算の高速化が可能となる。
特に、本実施例のように駆動手段（本実施例ではステージ駆動機構７）が、ステージ２の面に直交する軸を回転軸Ａｘとしてステージ２を回転駆動し、回転軸Ａｘからラミノ角傾いた斜め方向で放射線照射手段（Ｘ線検出器４）から放射線（本実施例ではＸ線）を照射するプレナーＣＴ(PCT: Planar Computed Tomography)を行う場合に適用することができる。プレナーＣＴ（ＰＣＴ）では通常のＣＴよりも高拡大率での放射線撮影が可能となる。高拡大率にすればするほど、放射線照射手段（Ｘ線検出器４）やステージ２の駆動誤差が投影画像に及ぼす影響が大きくなるので、３つのパラメータ（並進パラメータｔｘ，ｔｙ，ｔｚ）に限定することによる評価関数の収束精度への悪影響は小さくなる。
本実施例に係る画像処理方法の具体的は以下の通りである。具体的には、非線形最適化演算によって上述のパラメータを算出するパラメータ算出過程（図９のステップＴ１〜Ｔ５）を備え、パラメータ算出過程（ステップＴ１〜Ｔ５）は、特徴点３次元座標推定過程（ステップＴ１），再投影座標計算過程（ステップＴ２），再投影二乗誤差計算過程（ステップＴ３），収束判定過程（ステップＴ４）およびパラメータ更新過程（ステップＴ５）を備え、パラメータ更新過程（ステップＴ５）では、限定されたパラメータのみ更新する。特徴点３次元座標推定過程（ステップＴ１）では、校正用ファントムの投影画像から抽出された特徴点、およびパラメータの初期値に基づいて、特徴点の３次元座標を推定する。再投影座標計算過程（ステップＴ２）では、特徴点３次元座標推定過程（ステップＴ１）で推定した特徴点３次元座標に基づいて、再投影座標を計算する。再投影二乗誤差計算過程（ステップＴ３）では、再投影座標計算過程（ステップＴ２）で算出した再投影座標、および上述した特徴点に基づいて、再投影二乗誤差を算出する。収束判定過程（ステップＴ４）では、再投影二乗誤差計算過程（ステップＴ３）で算出した再投影二乗誤差の収束判定を行う。そして、パラメータ更新過程（ステップＴ５）では、再投影二乗誤差の値が収束していない場合には、再投影二乗誤差の値が小さくなるようにパラメータを更新するが、その更新の際には限定されたパラメータのみ更新する。このようにパラメータの更新時には、算出すべきパラメータの数を撮影条件に基づいて予め限定しているので、パラメータ更新過程（ステップＴ５）で更新するパラメータの数が少なくて済み、非線形最適化演算の収束精度を確保しつつ、計算を高速化することができる。
この発明は、上記実施形態に限られることはなく、下記のように変形実施することができる。
（１）上述した実施例では、放射線撮影装置として、Ｘ線検査装置を例に採って説明したが、放射線照射手段から照射されて対象物を透過した放射線を放射線検出手段で検出することにより得られた投影画像に基づいて放射線撮影を行う装置であれば、放射線についてはＸ線に限定されず、Ｘ線以外の放射線（α線、β線、γ線など）であってもよい。
（２）上述した実施例では、上述したように実装基板、多層基板のスルーホール／パターン／はんだ接合部、パレット上に配置された集積回路（ＩＣ）のような実装前の電子部品、金属などの鋳物、ビデオデッキのような成型品などに例示されるように、対象物に対する放射線撮影を行うのであればよい。
２ … ステージ
３ … Ｘ線管
４ … Ｘ線検出器
７ … ステージ駆動機構
８ … パラメータ算出部
１０ … 断層画像算出部
ｔｘ，ｔｙ，ｔｚ … 並進パラメータ
Ｏ … 対象物

Claims (6)

  1. 対象物を載置するステージと、
    そのステージを間に挟んで互いに対向するように配置された放射線照射手段および放射線検出手段と
    を備え、
    前記放射線照射手段から照射されて前記対象物を透過した放射線を前記放射線検出手段で検出することにより得られた投影画像に基づいて放射線撮影を行う放射線撮影装置であって、
    校正用ファントムもしくは前記対象物を前記ステージに載置した状態で、前記放射線照射手段、前記ステージの少なくともいずれかを駆動する駆動手段と、
    複数の前記校正用ファントムの投影画像に基づいて、放射線照射手段、ステージ、放射線検出手段の幾何関係を表すパラメータを非線形最適化演算により算出するパラメータ算出手段と、
    複数の前記対象物の投影画像および前記パラメータに基づいて、対象物の断層画像を算出し、当該算出時に前記対象物を撮影したときの撮影条件に適した前記パラメータを利用して再構成する断層画像算出手段と
    を備え、
    前記パラメータ算出時に、前記校正用ファントムを撮影したときの撮影条件に基づいて、非線形最適化演算によって算出するパラメータの数を減らして限定するパラメータ数限定手段と、
    過去に算出されたパラメータを更新するパラメータ更新手段と
    を備え、
    当該パラメータ更新手段は、前記パラメータ数限定手段で限定されたパラメータのみを更新することを特徴とする放射線撮影装置。
  2. 請求項１に記載の放射線撮影装置において、
    前記パラメータ数限定手段は、撮影条件を変更したときに変更された当該撮影条件に基づいて、非線形最適化演算によって算出するパラメータの数を減らして限定することを特徴とする放射線撮影装置。
  3. 請求項１または請求項２に記載の放射線撮影装置において、
    前記パラメータ数限定手段は、非線形最適化演算によって算出するパラメータを、前記放射線照射手段と前記ステージとの座標位置関係を表す３つのパラメータに限定することを特徴とする放射線撮影装置。
  4. 請求項３に記載の放射線撮影装置において、
    前記駆動手段は、前記ステージの面に直交する軸を回転軸としてステージを回転駆動し、
    前記回転軸からラミノ角傾いた斜め方向で前記放射線照射手段から放射線を照射するプレナーＣＴを行うことを特徴とする放射線撮影装置。
  5. 対象物を載置するステージと、
    そのステージを間に挟んで互いに対向するように配置された放射線照射手段および放射線検出手段と
    を備え、
    前記放射線照射手段から照射されて前記対象物を透過した放射線を前記放射線検出手段で検出することにより得られた投影画像に基づいて放射線撮影を行う放射線撮影装置における画像処理方法であって、
    複数の校正用ファントムの投影画像に基づいて、放射線照射手段、ステージ、放射線検出手段の幾何関係を表すパラメータを非線形最適化演算により算出する際に、前記校正用ファントムを撮影したときの撮影条件に基づいて、非線形最適化演算によって算出するパラメータの数を減らして限定するパラメータ数限定過程と、
    非線形最適化演算によって前記パラメータを算出するパラメータ算出過程と
    を備え、
    前記パラメータ算出過程は、
    過去に算出されたパラメータを更新するパラメータ更新過程を備え、
    当該パラメータ更新では、前記パラメータ数限定過程で限定されたパラメータのみを更新することを特徴とする放射線撮影装置における画像処理方法。
  6. 請求項５に記載の放射線撮影装置における画像処理方法において、
    前記パラメータ算出過程は、
    前記校正用ファントムの投影画像から抽出された特徴点、およびパラメータの初期値に基づいて、特徴点の３次元座標を推定する特徴点３次元座標推定過程と、
    前記特徴点３次元座標推定過程で推定した特徴点３次元座標に基づいて、再投影座標を計算する再投影座標計算過程と、
    前記再投影座標計算過程で算出した再投影座標、および前記特徴点に基づいて、再投影二乗誤差を算出する再投影二乗誤差計算過程と、
    前記再投影二乗誤差計算過程で算出した再投影二乗誤差の収束判定を行う収束判定過程と、
    を、さらに備え、
    前記パラメータ更新過程では、前記再投影二乗誤差の値が収束していない場合には、前記再投影二乗誤差の値が小さくなるようにパラメータを更新することを特徴とする放射線撮影装置における画像処理方法。
JP2013008437A 2013-01-21 2013-01-21 放射線撮影装置およびそれにおける画像処理方法 Active JP6131606B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2013008437A JP6131606B2 (ja) 2013-01-21 2013-01-21 放射線撮影装置およびそれにおける画像処理方法
US14/154,922 US9476844B2 (en) 2013-01-21 2014-01-14 Radiographic apparatus and an image processing method therefore
CN201410028530.7A CN103940835B (zh) 2013-01-21 2014-01-21 放射线摄影装置及其图像处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2013008437A JP6131606B2 (ja) 2013-01-21 2013-01-21 放射線撮影装置およびそれにおける画像処理方法

Publications (2)

Publication Number Publication Date
JP2014139533A JP2014139533A (ja) 2014-07-31
JP6131606B2 true JP6131606B2 (ja) 2017-05-24

Family

ID=51188593

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013008437A Active JP6131606B2 (ja) 2013-01-21 2013-01-21 放射線撮影装置およびそれにおける画像処理方法

Country Status (3)

Country Link
US (1) US9476844B2 (ja)
JP (1) JP6131606B2 (ja)
CN (1) CN103940835B (ja)

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107076684B (zh) 2014-09-02 2021-04-02 株式会社尼康 测量处理装置、测量处理方法和测量处理程序
EP3190403B1 (en) 2014-09-02 2023-03-15 Nikon Corporation Measurement processing device, measurement processing method, measurement processing program, and structure production method
FR3030753A1 (fr) * 2014-12-23 2016-06-24 Safran Fantome destine a etre utilise pour le controle de la qualite d'images tomographiques
WO2016124554A1 (en) * 2015-02-03 2016-08-11 Koninklijke Philips N.V. Object localization in projective x-ray images by geometric considerations
CN104665862B (zh) * 2015-02-16 2017-05-03 清华大学 在cbct中消除几何伪影的方法以及使用该方法的cbct***
US9606072B2 (en) * 2015-06-09 2017-03-28 Shimadzu Corporation Radiation inspecting apparatus
WO2017100475A1 (en) * 2015-12-09 2017-06-15 Integrated-X, Inc. Systems and methods for inspection using electromagnetic radiation
CN108738341B (zh) * 2015-12-15 2022-01-28 株式会社岛津制作所 螺旋ct装置
JP6861470B2 (ja) * 2016-03-04 2021-04-21 三星電子株式会社Ｓａｍｓｕｎｇ Ｅｌｅｃｔｒｏｎｉｃｓ Ｃｏ．，Ｌｔｄ． Ｘ線検査装置
US11145431B2 (en) * 2016-08-16 2021-10-12 Massachusetts Institute Of Technology System and method for nanoscale X-ray imaging of biological specimen
WO2018035171A1 (en) * 2016-08-16 2018-02-22 Massachusetts Institute Of Technology Nanoscale x-ray tomosynthesis for rapid analysis of integrated circuit (ic) dies
CN112313502A (zh) * 2018-03-07 2021-02-02 Ge传感与检测技术有限公司 计算机断层摄影应用的摆动补偿
JP2019158534A (ja) * 2018-03-12 2019-09-19 株式会社ミツトヨ 計測用ｘ線ｃｔ装置、及び、その断層画像生成方法
GB2574188B (en) * 2018-04-11 2022-03-09 E M & I Maritime Ltd Inspection of hazardous area equipment
TWI708217B (zh) * 2019-10-30 2020-10-21 行政院原子能委員會核能研究所 用於雙軸數位斷層合成造影系統的幾何校正方法及其系統
US11437218B2 (en) 2019-11-14 2022-09-06 Massachusetts Institute Of Technology Apparatus and method for nanoscale X-ray imaging
CN112683933B (zh) * 2020-11-30 2022-07-19 北京星航机电装备有限公司 一种增材制造多层结构检测射线灵敏度的测定方法
WO2022213274A1 (zh) * 2021-04-06 2022-10-13 中国科学技术大学 通用的ct对轴方法
CN115588011B (zh) * 2022-12-12 2023-03-10 佗道医疗科技有限公司 一种x射线机投影模型参数校准方法

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6837892B2 (en) * 2000-07-24 2005-01-04 Mazor Surgical Technologies Ltd. Miniature bone-mounted surgical robot
JP3694833B2 (ja) * 2002-05-10 2005-09-14 独立行政法人産業技術総合研究所 ユーセントリック型傾斜三次元ｘ線ｃｔ及びそれによる三次元画像の撮影方法
JP4415762B2 (ja) * 2003-06-10 2010-02-17 株式会社島津製作所 断層撮影装置
US7186023B2 (en) * 2003-06-10 2007-03-06 Shimadzu Corporation Slice image and/or dimensional image creating method
JP4537090B2 (ja) * 2004-02-19 2010-09-01 東芝Ｉｔコントロールシステム株式会社 トモシンセシス装置
FR2881941B1 (fr) * 2005-02-16 2007-04-13 Gen Electric Procede de determination des parametres geometriques d'un dispositif d'imagerie par rayon x
EP1861734A2 (de) * 2005-03-09 2007-12-05 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Korrektur von nichtlinearitäten eines abbildungssytems durch ein a-priori wissen bei durchstrahlungsaufnahmen
JP2007048068A (ja) 2005-08-10 2007-02-22 Canon Inc 情報処理方法および装置
CN100536770C (zh) * 2007-03-29 2009-09-09 新奥博为技术有限公司 一种磁共振图像引导下的手术***及手术导航方法
JP5184830B2 (ja) 2007-07-06 2013-04-17 株式会社トプコン 位置計測装置及び位置計測方法
US7787689B2 (en) 2007-07-06 2010-08-31 Topcon Corporation Location measuring device and method
CN101361660B (zh) * 2008-05-16 2011-03-09 深圳先进技术研究院 一种多磁性目标的定位方法及定位***
JP5396751B2 (ja) 2008-06-18 2014-01-22 オムロン株式会社 Ｘ線利用の基板検査装置
CN102652674B (zh) 2011-03-04 2014-02-19 首都师范大学 一种消除ct图像中的几何伪影的方法和***

Also Published As

Publication number Publication date
US9476844B2 (en) 2016-10-25
CN103940835B (zh) 2017-01-11
US20140205058A1 (en) 2014-07-24
JP2014139533A (ja) 2014-07-31
CN103940835A (zh) 2014-07-23

Similar Documents

Publication Publication Date Title
JP6131606B2 (ja) 放射線撮影装置およびそれにおける画像処理方法
US9442080B2 (en) Method and apparatus for generating a three-dimensional model of a region of interest using an imaging system
KR102493849B1 (ko) 화상 취득 장치 및 화상 취득 방법, 그리고 화상 보정 프로그램
JP2006175213A (ja) ３次元画像処理装置
US9606072B2 (en) Radiation inspecting apparatus
JP6383707B2 (ja) 撮影画像のずれ補正装置、撮影画像のずれ補正方法および撮影画像のずれ補正プログラム
CN117871557A (zh) 检查装置、以及检查区域的确定方法
JP2012112790A (ja) Ｘ線ｃｔ装置
JP5938860B2 (ja) 放射線撮影装置および断層画像補正方法
JP6676023B2 (ja) 検査位置の特定方法及び検査装置
US20080275666A1 (en) Region-In-Object Measuring System, Computing Device for Measuring a Region-In-Object, Program for Measuring Region-In-Object and Computer Readable Recording Medium on Which the Program Is Recorded
WO2011089658A1 (ja) 補正用治具を用いた検査装置の補正方法、補正用治具を搭載した検査装置
JP2009237704A (ja) ３次元画像化方法およびｘ線断層像撮影装置
JP2013160710A (ja) 放射線撮影装置および投影画像処理方法
JP6015404B2 (ja) 放射線検査装置
JP4926734B2 (ja) 放射線検査装置、放射線検査方法および放射線検査プログラム
JP2005127886A (ja) 傾斜型ｘ線ｃｔ装置における回転中心軸の較正方法および傾斜型ｘ線ｃｔ装置
JP4400737B2 (ja) 画像処理方法および放射線断層撮像装置
JP2023018659A (ja) Ｘ線撮像において散乱防止グリッド・アーティファクトを除去する方法及びシステム
JP2022141454A (ja) Ｘ線検査装置およびｘ線検査方法
WO2020012520A1 (ja) 医用ｘ線画像処理装置およびｘ線画像撮影装置
JP5999516B2 (ja) 放射線撮影装置
JP2022140949A (ja) Ｘ線検査装置およびｘ線検査方法
JP5851200B2 (ja) 基板検査装置、基板検査方法および基板検査プログラム
JP2018179657A (ja) Ｘ線検査装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20151005

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20160729

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20160805

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20161003

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20170321

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20170403

R151 Written notification of patent or utility model registration

Ref document number: 6131606

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151