JPWO2019008620A1 - Ｘ線ｃｔ装置 - Google Patents

Ｘ線ｃｔ装置 Download PDF

Info

Publication number
JPWO2019008620A1
JPWO2019008620A1 JP2019528189A JP2019528189A JPWO2019008620A1 JP WO2019008620 A1 JPWO2019008620 A1 JP WO2019008620A1 JP 2019528189 A JP2019528189 A JP 2019528189A JP 2019528189 A JP2019528189 A JP 2019528189A JP WO2019008620 A1 JPWO2019008620 A1 JP WO2019008620A1
Authority
JP
Japan
Prior art keywords
ray
marker
flat plate
subject
image
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.)
Granted
Application number
JP2019528189A
Other languages
English (en)
Other versions
JP7164524B2 (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
Publication of JPWO2019008620A1 publication Critical patent/JPWO2019008620A1/ja
Priority to JP2021187680A priority Critical patent/JP7251602B2/ja
Application granted granted Critical
Publication of JP7164524B2 publication Critical patent/JP7164524B2/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
    • 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/401—Imaging image processing

Landscapes

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

Abstract

Ｘ線源と回転ステージ１３との間には、マーカＭが形成された平板２１と、平板２１を支持する支持部２２から成るマーカ部材が配置されている。平板２１におけるマーカＭの形成位置は、Ｘ線検出器１２の検出範囲内において、マーカＭのいずれもが被写体の投影像と重畳しない領域、かつ、Ｘ線焦点の移動があっても常に検出範囲内に含まれる領域の中で、マーカＭ同士の距離が最も離れる位置としている。また、支持部２２の長さは、投影像において、平板２１およびマーカＭが被写体Ｗに重畳することがないＸ線検出器の検出範囲の側端の位置となる長さに調整される。

Description

この発明は、各種製品の内部構造の観察および３次元形状の測定を非破壊で行うＸ線ＣＴ装置に関する。
産業用のＸ線ＣＴ装置は、互いに対向配置させたＸ線源とＸ線検出器との間に、工業製品等の被写体を載置する回転ステージを配置し、被写体の周囲の各方向からのＸ線投影データを収集して断層画像を再構成することにより、その被写体の内部構造を３次元的に観察するものである（特許文献１参照）。このような被写体の３次元構造を捉えることができる性質から、近年では、観察用途だけでなく、３次元形状測定にＸ線ＣＴ装置が用いられている。
断層画像は、ＣＴスキャン時のＸ線焦点、被写体、検出器の幾何的な情報を基に投影像を逆投影して再構成される。このため、３次元形状測定用のＸ線ＣＴ装置においては、幾何情報に生じた誤差が、再構成像に空間的な歪みを生じさせ、寸法測定精度を低下させる原因となる。特にＸ線焦点は、Ｘ線照射中にＸ線管の熱膨張やターゲットの劣化等により常に変動する。例えば、Ｘ線焦点が１０μｍ変動したと仮定すると、Ｘ線検出器で検出される１００倍に拡大された投影像では、１ｍｍの誤差となる。３次元形状測定用としてＸ線ＣＴ装置が満たすべき寸法測定精度を実現できるように被写体の断層画像を空間歪みなく再構成するためには、ＣＴスキャン中のＸ線焦点の位置を３次元的に検出して、断層画像の再構成に必要な情報を補正する、あるいは、Ｘ線焦点の変動を抑制する必要がある。
特許文献２、特許文献３および非特許文献１には、ＣＴスキャン中のＸ線焦点の位置を検出する手法として、マーカの投影像から２次元の焦点移動量を求める技術が記載されている。また、特許文献４には、装置の温度変化の影響による検出精度の低下を抑制するために、Ｘ線源をハウジングに収容し、ハウジング内に温度調節された気体を供給することでＸ線源を含むハウジングの熱を冷却する構造が記載されている。
特開２００５−３５１８７９号公報 国際公開２００９／０３６９８３号 特許第３７４３５９４号 特許第５８５００５９号
Ｆｒｅｄｅｒｉｋ Ｖｏｇｅｌｅｒ、Ｗｅｓｌｅｙ Ｖｅｒｈｅｅｃｋｅ、Ａｎｄｒｅ Ｖｏｅｔ、Ｊｅａｎ−Ｐｉｅｒｒｅ Ｋｒｕｔｈ、Ｗｉｍ Ｄｅｗｕｌｆ、Ｐｏｓｉｔｉｎａｌ Ｓｔａｂｉｌｉｔｙ ｏｆ ２Ｄ Ｘ−ｒａｙ Ｉｍａｇｅｓ ｆｏｒ Ｃｏｍｐｕｔｅｒ Ｔｏｍｏｇｒａｐｈｙ、Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｓｙｍｐｏｓｉｕｍ ｏｎ Ｄｉｇｉｔａｌ Ｉｎｄｕｓｔｒｉａｌ Ｒａｄｉｏｌｏｇｙ ａｎｄ Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ−Ｍｏ．３．３
特許文献２に記載のＸ線焦点の位置検出では、検出器を介して得られるマーカの投影像が大きく、被写体の撮像視野がマーカにより阻害されるため、被写体とマーカを同時に撮像することができない。このため、被写体をスキャンするときには、マーカを退避させている。また、マーカを撮像するときには、マーカと検出器との位置関係が同じになるようにマーカの基準位置への配置を再現する必要がある。このため、精度の高いマーカの位置決め機構が必要となり、製作面での負担が大きくなる。
非特許文献１に記載のＸ線焦点の位置検出では、被写体を載置する回転ステージの近傍に被写体の視野を阻害しないようマーカを設置し、マーカを被写体と同時に撮像することでＣＴスキャン中のＸ線焦点の位置を常時検出できる。しかしながら、検出するＸ線焦点の位置は検出器面上の２方向のみである。
特許文献３に記載のＸ線焦点の位置検出では、マーカを被写体に貼設することでマーカと被写体の投影像を同時に得ることができ、ＣＴスキャン中のＸ線焦点の位置を常時検出することができる。しかしながら、検出器を介して得られる投影像は、被写体にマーカが重なったものとなるため、マーカの検出精度が低下するとともに、被写体の断層画像にマーカの影響が及ぶ。
また、Ｘ線焦点の変動を抑制するために、Ｘ線源側の温度を制御するとしても、特許文献４に記載のように、Ｘ線源をハウジングに収容して固定し、そのハウジングをさらに冷却装置に接続する構造を採用すると、装置の制作コストが高くなる。
この発明は、上記課題を解決するためになされたものであり、被写体のＣＴスキャン中に変動するＸ線の焦点位置を３次元的に検出し、断層画像を空間歪なく再構成することが可能なＸ線ＣＴ装置を提供することを目的とする。
請求項１に記載の発明は、被写体に複数の角度からＸ線を照射して取得した投影データに基づいて前記被写体の断層画像を再構成するＸ線ＣＴ装置であって、Ｘ線を発生するＸ線源と、前記Ｘ線源に対向配置されるＸ線を検出するためのＸ線検出器と、前記Ｘ線源と前記Ｘ線検出器との間に配置される前記被写体を載置するための回転ステージと、前記Ｘ線検出器が検出した前記被写体の投影データに基づいて演算処理を実行する制御装置と、マーカが少なくとも２か所設けられた平板と、前記平板を支持する支持部とから成り、前記回転ステージと前記Ｘ線源との間の位置であって、ＣＴスキャンの実行中に前記マーカが前記Ｘ線検出器の検出範囲内に含まれる位置、かつ、前記平板が前記被写体の投影像に重畳しない位置に配置されるマーカ部材と、を備え、前記制御装置は、前記Ｘ線源の焦点と前記マーカとを通る直線が前記Ｘ線検出器と交差する点である特徴点を、前記マーカの投影像を画像処理することにより検出し、２次元Ｘ線画像上の前記特徴点の座標を求めるマーカ特徴点検出部と、前記マーカ特徴点検出部により求めた２つの異なるフレームの前記特徴点の座標を用いて、前記Ｘ線源の焦点が前記２つの異なるフレームが得られる間に移動した量を３次元的に算出する焦点移動量算出部と、を備え、前記Ｘ線源の焦点移動量に基づいて前記被写体の断層画像を再構成するときの座標系を修正することを特徴とする。
請求項２に記載の発明は、請求項１に記載のＸ線ＣＴ装置において、前記マーカ部材は、前記Ｘ線検出器に対する前記平板の設置角度が変化しないように固定され、前記焦点移動量算出部は、前記平板の設置角度情報を利用して、前記Ｘ線源の焦点の移動量を３次元的に算出する。
請求項３に記載の発明は、請求項１に記載のＸ線ＣＴ装置において、前記マーカ部材は、前記平板を前記Ｘ線源のＸ線照射口に近接させて配置される。
請求項４に記載の発明は、請求項１に記載のＸ線ＣＴ装置において、前記マーカ特徴点検出部は、前記マーカの投影像を画像処理して輝度重心点を求め、前記輝度重心点を前記特徴点の座標とする。
請求項５に記載の発明は、請求項１に記載のＸ線ＣＴ装置において、前記マーカ部材における前記平板は、Ｘ線を減衰させる材料から成る。
請求項６に記載の発明は、請求項５に記載のＸ線ＣＴ装置において、前記マーカは、前記平板内に設けられた空洞、あるいは、前記平板に形成された貫通孔または凹部である。
請求項７に記載の発明は、請求項５に記載のＸ線ＣＴ装置において、前記マーカは、前記平板の表面に配設されたＸ線を減衰させる材料から成る円柱または円錐台形状の部材である。
請求項１から請求項７に記載の発明によれば、Ｘ線焦点位置を検出するためのマーカを少なくとも２カ所設けた平板を有するマーカ部材を備え、少なくとも２つのマーカを被写体の投影像に重畳させることなく同時にスキャンして、ＣＴスキャン中のＸ線焦点位置を３次元的に検出し、Ｘ線源の焦点移動量に基づいて断層画像を再構成する際の座標系を補正することから、断層画像を空間歪なく再構成することが可能となる。このため、３次元形状測定用として高い寸法測定精度を実現することができる。
請求項２に記載の発明によれば、Ｘ線検出器に対する平板の設置角度を固定し、設置角度をＸ線の焦点移動量の算出パラメータとして加えることにより、ＣＴスキャン中のＸ線焦点の位置を、より正確に求めることが可能となる。
請求項３に記載の発明によれば、マーカ部材は、平板をＸ線源のＸ線照射口に近接させて配置されることから、投影像におけるマーカの拡大率を大きくすることができる。これにより、平板に設けるマーカの大きさを小さくすることが可能となる。
請求項４に記載の発明によれば、マーカの投影像を画像処理することにより、輝度重心点を特徴点の座標とすることから、マーカの投影像が変形している場合でも、特徴点の座標を正確に決めることができる。
請求項５から請求項７に記載の発明によれば、平板の厚みと異なる厚みを持つマーカを設けることで、投影像において平板の板部分とマーカ部分とでＸ線透過強度に明確な差を生じさせ、画像処理による特徴点の検出を容易に行うことができる。
請求項６に記載の発明によれば、平板にマーカを容易に形成することができる。
この発明に係るＸ線ＣＴ装置の概要図である。 Ｘ線検出器１２の検出面におけるマーカＭの位置を説明する概要図である。 この発明に係るＸ線ＣＴ装置の主要な制御系を説明するブロック図である。 マーカＭの投影像の取得からＣＴ画像再構成に至る処理を示すフローチャートである。
この発明の実施の形態を、図面に基づいて説明する。図１は、この発明に係るＸ線ＣＴ装置の概要図である。図２は、Ｘ線検出器１２の検出面におけるマーカＭの位置を説明する概要図である。図３は、この発明に係るＸ線ＣＴ装置の主要な制御系を説明するブロック図である。
Ｘ線ＣＴ装置は、Ｘ線を発生させるＸ線源１１と、Ｘ線源１１に対向配置されたＸ線検出器１２と、被写体Ｗを載置するための回転ステージ１３と、Ｘ線ＣＴ装置全体の制御を行う制御装置３０を備える。制御装置３０は、パーソナルコンピュータ４０と、Ｘ線コントローラ３１と、ステージコントローラ３２とから構成される。
Ｘ線源１１とＸ線検出器１２は、図示を省略した支持機構により定盤１６に固定されている。Ｘ線源１１とＸ線検出器１２との間には、回転軸Ｒを中心に回転する回転ステージ１３が配置されている。回転ステージ１３は、ステージ移動機構１５により定盤１６上をＸ線光軸に沿う方向に移動する。回転ステージ１３とＸ線源１１との距離をステージ移動機構１５により変更することで、Ｘ線検出器１２で検出される被写体Ｗの投影像の拡大率が変更される。なおここで、本明細書では、図２に記載されているように、Ｘ軸はＸ線光軸に沿う水平方向であり、Ｚ軸は上下（鉛直）方向であり、Ｙ軸はＸ軸とＺ軸に直交する方向である。
Ｘ線源１１は、Ｘ線検出器１２に向けてＸ線を円錐状に照射するＸ線管を有する。被写体Ｗの材質やＸ線透過特性に応じて、Ｘ線管に供給する管電圧、管電流をＸ線コントローラ３１により制御する。Ｘ線コントローラ３１は、パーソナルコンピュータ４０の制御下に置かれている。Ｘ線検出器１２により検出された透過Ｘ線像は、パーソナルコンピュータ４０に取り込まれる。パーソナルコンピュータ４０では、投影データに基づいてＣＴ画像再構成部５３（図３参照）により回転軸Ｒに直交するＸ−Ｙ平面に沿った面でスライスした断層画像が構築される。
回転ステージ１３およびステージ移動機構１５は、それぞれ独立の駆動モータを有し、それらの駆動モータは、ステージコントローラ３２から供給される駆動信号によって制御される。ステージコントローラ３２は、パーソナルコンピュータ４０の制御下に置かれている。
また、Ｘ線源１１と回転ステージ１３との間には、マーカＭが形成された平板２１と、平板２１を支持する支持部２２から成るマーカ部材２０が配置されている。マーカ部材２０は、回転ステージ１３よりもＸ線源１１に近い位置に配置される。
平板２１は、支持部２２で支持した状態で変形することがない厚みを持つ薄板により作成される。そして、平板２１には、２カ所にマーカＭが設けられている。この実施形態では、マーカＭは、矩形状の平板２１の左右対称な位置に穿設された貫通孔である。なお、貫通孔に代えて、平板２１の材料と同様のＸ線を減衰させる材料により作成された円柱や円錐台形状の部材を平板２１の表面に配設してマーカＭとしてもよい。また、貫通孔の形状も、円柱孔に限らず、孔の内面が傾斜しているテーパ孔であってもよい。さらに、貫通孔ではなく、平板２１に凹部を形成することでマーカＭとしてもよく、平板２１の材質によっては、平板内に空洞（ボイド）を設けることでマーカＭとしてもよい。すなわち、マーカＭは、平板２１の厚みと異なる厚みを持つことで、投影像において平板２１とのコントラストを得ることができるものであればよく、マーカＭの形状は、平板２１の板部分とマーカ部分とでＸ線透過強度に明確な差を生じさせるものであればよい。また、この実施形態では、マーカＭを２カ所設けているが、マーカＭの数をさらに増やしてもよい。
平板２１は、支持部２２に固定されており、設置角度が変化しないようになっている。Ｘ線源１１の焦点からＸ線を照射して被写体Ｗの投影データを取得する際には、図２に示すように、マーカＭがＸ線検出器１２の検出範囲内に含まれるように、マーカ部材２０の位置をＸ線源１１に近い側に固定する。これにより、Ｘ線源１１のＸ線照射口の前の近接した位置に平板２１が固定される。なお、平板２１の位置を規定する支持部２２の長さは、図２に示すように、マーカＭを含む平板２１の投影像が被写体Ｗの投影像と重畳することがないＸ線検出器１２の検出範囲の側端の位置となる長さとなっている。このように、支持部２２を介しマーカ部材２０を定盤１６に配置することで、マーカＭが形成された平板２１を、Ｘ線源１１の筐体とは独立して定位置に固定することができる。
平板２１におけるマーカＭの形成位置は、マーカＭのいずれもが被写体の投影像と重畳しない範囲でマーカＭ同士の離間距離が長くなるようにすることで、後述するＸ線焦点の３次元的な移動量の算出精度がよくなる。
Ｘ線焦点から円錐状に照射され、平板２１および回転ステージ１３上の被写体Ｗを透過したＸ線は、Ｘ線検出器１２により検出され、検出データはパーソナルコンピュータ４０に入力される。
パーソナルコンピュータ４０は、ＲＡＭ、ＲＯＭなどのメモリ４１と、各種の演算処理を実行するＣＰＵなどの演算装置４２と、ＣＴスキャンにより取得された投影データを保存するＨＤＤなどの記憶装置４３と、Ｘ線コントローラ３１やステージコントローラ３２に制御信号を送信するための通信部４４を備える。メモリ４１、演算装置４２、記憶装置４３および通信部４４は、相互のデータ通信を可能とする内部バス６１により接続されている。メモリ４１には、演算装置４２を動作させて機能を実現するプログラムが格納されている。図３においては、パーソナルコンピュータ４０にインストールされているプログラムを機能ブロックとして示している。この実施形態では、機能ブロックとして、マーカ特徴点検出部５１と、焦点移動量算出部５２と、ＣＴ画像再構成部５３とを備える。また、パーソナルコンピュータ４０には、演算装置４２の動作により構築された断層画像などを表示する表示装置４８と、各種指令を装置に与えるためのマウスやキーボードなどの入力装置４９が、通信部４４を介して接続されている。
次に、上述した構成のＸ線ＣＴ装置によりＣＴスキャンを実行する際に、スキャン中のＸ線の焦点位置を検出する手順について説明する。断層画像は、ＣＴスキャン時のＸ線焦点、被写体Ｗ、検出器の位置関係を含む幾何情報を基に、投影像を逆投影して再構成される。したがって、この発明では、断層画像の再構成に必要な幾何情報の一つであるＸ線の焦点位置の移動を追跡し、移動量を断層画像の再構成に反映するようにしている。図４は、マーカＭの投影像の取得からＣＴ画像再構成に至る処理を示すフローチャートである。
まず、マーカＭが形成された平板２１がＸ線源１１のＸ線照射口の直前に配置されるように、マーカ部材２０を定盤１６に固定する（ステップＳ１１）。このとき、Ｘ線検出器１２の検出面と平板２１の板面とが平行になるように、マーカ部材２０を固定する。これにより、マーカＭの位置、平板２１のＸ線検出器１２の検出面に対する設置角度が固定される。そして、Ｘ線焦点の基準となるＸ線焦点位置（基準焦点位置）を決める（ステップＳ１２）。なお、Ｘ線焦点位置は、Ｘ線源１１におけるＸ線管において、熱電子がターゲットに衝突してＸ線を発生させる位置であって、この明細書では一つの点であると見なしている。
基準焦点位置でのマーカＭの投影像を取得し（ステップＳ１３）、得られた投影像からマーカＭの特徴点を検出する（ステップＳ１４）。すなわち、被写体ＷへのＣＴスキャンを開始する前に、Ｘ線管の熱変形が起こる前の基準となるＸ線焦点位置でのマーカＭの特徴点（基準マーカ特徴点）の位置を取得する。
マーカＭの特徴点は、Ｘ線焦点とマーカＭを通る直線がＸ線検出器１２の検出面と交差する点である。Ｘ線検出器１２より得られた投影像を画像処理して決定したマーカＭの投影像の輝度重心点を、特徴点の座標としている。なお、マーカＭの投影像の輝度重心点は、２次元Ｘ線画像上の点であり、マーカＭの特徴点の移動は、このＸ線検出器１２の検出面であるｕ−ｖ平面座標におけるベクトルで表される。また、マーカＭの特徴点の座標を輝度重心点とすることで、マーカＭへのＸ線の入射角との関係で、平板２１に形成された円形のマーカＭが投影像では楕円形に変形しても、マーカＭの特徴点の座標を正確に求めることができる。
続いて、被写体ＷとマーカＭとを同時に撮像しながらスキャンする（ステップＳ１５）。ＣＴスキャンに際しては、Ｘ線源１１からＸ線を照射しながら、被写体Ｗを載置した回転ステージ１３を回転させて、所定の回転角度ごとにＸ線検出器１２からパーソナルコンピュータ４０に入力された投影像が、フレームデータとして記憶装置４３に記憶される。なお、フレームデータとは投影像の１画面分のデータのことであり、得られた時刻の順番にしたがって記憶装置４３に記憶されるものである。
得られた投影像を画像処理することにより、マーカＭの特徴点を検出する（ステップＳ１６）。ＣＴスキャンでは、設定されたスキャンモードで設定されたビュー数に応じたサンプリングピッチに対応する投影データが収集される。スキャン中は、被写体Ｗには複数の角度からＸ線が照射されるが、Ｘ線検出器１２に対する平板２１の位置は固定されていることから、各フレームから、マーカＭの投影像が得られる。スキャンの経過に伴い、Ｘ線管の熱変形などによりＸ線焦点位置が移動すると、それは、Ｘ線検出器１２の検出面におけるｕ−ｖ平面座標上でのマーカＭの特徴点の移動として観察できる。フレームデータに対して画像処理を実行することで、異なる２つのフレームのマーカＭの特徴点を検出する（ステップＳ１６）。なお、マーカＭの特徴点の検出は、全てのフレームデータに対して実行してもよく、所定の間隔で抽出したフレームデータに対して行ってもよい。すなわち、後述するステップＳ１７で、２つの異なるフレーム間で、どれだけマーカＭの特徴点が移動したかを計算することができればよい。
なお、ステップＳ１４およびステップＳ１６におけるマーカＭの特徴点の検出は、演算装置４２がメモリ４１におけるマーカ特徴点検出部５１から読み込んだプログラムを実行することにより実現される。
ステップＳ１４で得られたＸ線焦点が基準位置にあるときのマーカＭの特徴点のベクトルと、ステップＳ１６で得られたＣＴスキャン中のあるフレームのマーカＭの特徴点のベクトルから、特徴点がどれだけどの方向に移動したかを求める（ステップＳ１７）。すなわち、ある時間経過後のマーカＭの特徴点の座標位置の変化から、その時間内に、Ｘ線焦点がどれだけの量どの方向に移動したかが算出される。このＸ線焦点移動量の算出は、演算装置４２がメモリ４１における焦点移動量算出部５２から読み込んだプログラムを実行することにより実現される。
なお、この実施形態では、２つのマーカＭの特徴点の位置を求めており、２つのマーカＭの位置関係と各マーカＭの特徴点の移動量から、Ｘ線焦点のＸ線光軸に沿う方向（Ｘ方向）の移動も計算することができる。このように、少なくとも２つのマーカＭの２つの異なるフレームで検出した特徴点の座標を用いて、Ｘ線の焦点が当該時間内で移動した量を、移動した方向も含めて３次元的に算出することができる。
このＸ線ＣＴ装置では、ＣＴスキャン中のＸ線検出器１２の検出面に対して平板２１の板面をほぼ平行に設置している。このため、Ｘ線焦点の移動量を算出する際に、特に設置角度情報をパラメータとして設定しなくても、ＣＴスキャン中のＸ線焦点の移動量を算出することができる。
一方、Ｘ線焦点の移動量の３次元的算出において、Ｘ線検出器１２の検出面に対する平板２１の設置角度をパラメータとして加えることもできる。Ｘ線検出器１２の検出面に対する平板２１の設置角度が平行からずれるほど、２つのマーカＭの各々とＸ線検出器１２の検出面との距離に差が生じ、投射像における両者の拡大率が異なることになる。Ｘ線焦点の空間移動の検出に、既知パラメータとして予め取得しておいた設置角度情報を利用することで、Ｘ線焦点移動の検出精度をさらに向上させることができる。
Ｘ線検出器１２の検出面に対する平板２１の設置角度は、例えば、マーカ部材２０を定盤１６に固定したときに、レーザ変位計などを使用して取得しておけばよい。なお、このＸ線ＣＴ装置のようにＸ線源１１とＸ線検出器１２の位置が固定されており、ＣＴスキャン中のマーカ部材２０の位置が定盤１６に固定されているに等しい状態の装置では、Ｘ線検出器１２とマーカ部材２０を配置して以降、Ｘ線ＣＴ装置の使用時においてこの設置角度が変わることがないため、レーザ変位計などを使用した設置角度情報の取得は一度で済むことになる。
Ｘ線の焦点移動量が求まると、その焦点移動量に基づいて、ＣＴ画像再構成の座標系を修正する（ステップＳ１８）。Ｘ線の焦点位置は、断層画像の再構成に必要な幾何情報の一つであり、焦点移動量に基づいてＣＴ画像再構成の座標系を修正することで、ＣＴスキャン中の幾何情報の誤差が修正される。しかる後、ＣＴ画像再構成を実行し（ステップＳ１９）、被写体Ｗの断層画像を再構成する。なお、ＣＴ画像再構成の座標系の修正（ステップＳ１８）およびＣＴ画像再構成の実行（ステップＳ１９）は、演算装置４２がメモリ４１におけるＣＴ画像再構成部５３から読み込んだプログラムを実行することにより実現される。
上述したように、この発明のＸ線ＣＴ装置では、少なくとも２つのマーカＭを被写体Ｗと同時に撮像しながらＣＴスキャンし、被写体Ｗの投影データと同じ時間軸でマーカＭの特徴点を検出することで、ＣＴスキャン中のＸ線焦点位置を３次元的に検出している。これにより、ＣＴスキャン中のＸ線焦点の移動を断層画像の再構成に反映できることから、断層画像を空間歪みなく再構成すること可能となる。したがって、この発明のＸ線ＣＴ装置を高い寸法測定精度が要求される３次元形状測定用に使用することも可能となる。
１１ Ｘ線源
１２ Ｘ線検出器
１３ 回転ステージ
１５ ステージ移動機構
１６ 定盤
２０ マーカ部材
２１ 平板
２２ 支持部
３０ 制御装置
３１ Ｘ線コントローラ
３２ ステージコントローラ
４０ パーソナルコンピュータ
４１ メモリ
４２ 演算装置
４３ 記憶装置
４４ 通信部
４８ 表示装置
４９ 入力装置
５１ マーカ特徴点検出部
５２ 焦点移動量算出部
５３ ＣＴ画像再構成部
Ｍ マーカ
Ｗ 被写体
この発明は、各種製品の内部構造の観察および３次元形状の測定を非破壊で行うＸ線ＣＴ装置に関する。
産業用のＸ線ＣＴ装置は、互いに対向配置させたＸ線源とＸ線検出器との間に、工業製品等の被写体を載置する回転ステージを配置し、被写体の周囲の各方向からのＸ線投影データを収集して断層画像を再構成することにより、その被写体の内部構造を３次元的に観察するものである（特許文献１参照）。このような被写体の３次元構造を捉えることができる性質から、近年では、観察用途だけでなく、３次元形状測定にＸ線ＣＴ装置が用いられている。
断層画像は、ＣＴスキャン時のＸ線焦点、被写体、検出器の幾何的な情報を基に投影像を逆投影して再構成される。このため、３次元形状測定用のＸ線ＣＴ装置においては、幾何情報に生じた誤差が、再構成像に空間的な歪みを生じさせ、寸法測定精度を低下させる原因となる。特にＸ線焦点は、Ｘ線照射中にＸ線管の熱膨張やターゲットの劣化等により常に変動する。例えば、Ｘ線焦点が１０μｍ変動したと仮定すると、Ｘ線検出器で検出される１００倍に拡大された投影像では、１ｍｍの誤差となる。３次元形状測定用としてＸ線ＣＴ装置が満たすべき寸法測定精度を実現できるように被写体の断層画像を空間歪みなく再構成するためには、ＣＴスキャン中のＸ線焦点の位置を３次元的に検出して、断層画像の再構成に必要な情報を補正する、あるいは、Ｘ線焦点の変動を抑制する必要がある。
特許文献２、特許文献３および非特許文献１には、ＣＴスキャン中のＸ線焦点の位置を検出する手法として、マーカの投影像から２次元の焦点移動量を求める技術が記載されている。また、特許文献４には、装置の温度変化の影響による検出精度の低下を抑制するために、Ｘ線源をハウジングに収容し、ハウジング内に温度調節された気体を供給することでＸ線源を含むハウジングの熱を冷却する構造が記載されている。
特開２００５−３５１８７９号公報 国際公開２００９／０３６９８３号 特許第３７４３５９４号 特許第５８５００５９号
Ｆｒｅｄｅｒｉｋ Ｖｏｇｅｌｅｒ、Ｗｅｓｌｅｙ Ｖｅｒｈｅｅｃｋｅ、Ａｎｄｒｅ Ｖｏｅｔ、Ｊｅａｎ−Ｐｉｅｒｒｅ Ｋｒｕｔｈ、Ｗｉｍ Ｄｅｗｕｌｆ、Ｐｏｓｉｔｉｎａｌ Ｓｔａｂｉｌｉｔｙ ｏｆ ２Ｄ Ｘ−ｒａｙ Ｉｍａｇｅｓ ｆｏｒ Ｃｏｍｐｕｔｅｒ Ｔｏｍｏｇｒａｐｈｙ、Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｓｙｍｐｏｓｉｕｍ ｏｎ Ｄｉｇｉｔａｌ Ｉｎｄｕｓｔｒｉａｌ Ｒａｄｉｏｌｏｇｙ ａｎｄ Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ−Ｍｏ．３．３
特許文献２に記載のＸ線焦点の位置検出では、検出器を介して得られるマーカの投影像が大きく、被写体の撮像視野がマーカにより阻害されるため、被写体とマーカを同時に撮像することができない。このため、被写体をスキャンするときには、マーカを退避させている。また、マーカを撮像するときには、マーカと検出器との位置関係が同じになるようにマーカの基準位置への配置を再現する必要がある。このため、精度の高いマーカの位置決め機構が必要となり、製作面での負担が大きくなる。
非特許文献１に記載のＸ線焦点の位置検出では、被写体を載置する回転ステージの近傍に被写体の視野を阻害しないようマーカを設置し、マーカを被写体と同時に撮像することでＣＴスキャン中のＸ線焦点の位置を常時検出できる。しかしながら、検出するＸ線焦点の位置は検出器面上の２方向のみである。
特許文献３に記載のＸ線焦点の位置検出では、マーカを被写体に貼設することでマーカと被写体の投影像を同時に得ることができ、ＣＴスキャン中のＸ線焦点の位置を常時検出することができる。しかしながら、検出器を介して得られる投影像は、被写体にマーカが重なったものとなるため、マーカの検出精度が低下するとともに、被写体の断層画像にマーカの影響が及ぶ。
また、Ｘ線焦点の変動を抑制するために、Ｘ線源側の温度を制御するとしても、特許文献４に記載のように、Ｘ線源をハウジングに収容して固定し、そのハウジングをさらに冷却装置に接続する構造を採用すると、装置の制作コストが高くなる。
この発明は、上記課題を解決するためになされたものであり、被写体のＣＴスキャン中に変動するＸ線の焦点位置を３次元的に検出し、断層画像を空間歪なく再構成することが可能なＸ線ＣＴ装置を提供することを目的とする。
請求項１に記載の発明は、被写体に複数の角度からＸ線を照射して取得した投影データに基づいて前記被写体の断層画像を再構成するＸ線ＣＴ装置であって、Ｘ線を発生するＸ線源と、前記Ｘ線源に対向配置されるＸ線を検出するためのＸ線検出器と、前記Ｘ線源と前記Ｘ線検出器との間に配置される前記被写体を載置するための回転ステージと、前記Ｘ線検出器が検出した前記被写体の投影データに基づいて演算処理を実行する制御装置と、マーカが少なくとも２か所設けられた平板と、前記平板を支持する支持部とから成り、前記回転ステージと前記Ｘ線源との間の位置であって、ＣＴスキャンの実行中に前記マーカが前記Ｘ線検出器の検出範囲内に含まれる位置、かつ、前記平板が前記被写体の投影像に重畳しない位置に配置されるマーカ部材と、を備え、前記制御装置は、前記Ｘ線源の焦点と前記マーカとを通る直線が前記Ｘ線検出器と交差する点である特徴点を、前記マーカの投影像を画像処理することにより検出し、２次元Ｘ線画像上の前記特徴点の座標を求めるマーカ特徴点検出部と、前記マーカ特徴点検出部により求めた２つの異なるフレームの前記特徴点の座標を用いて、前記Ｘ線源の焦点が前記２つの異なるフレームが得られる間に移動した量を３次元的に算出する焦点移動量算出部と、を備え、前記Ｘ線源の焦点移動量に基づいて前記被写体の断層画像を再構成するときの座標系を修正することを特徴とする。
請求項２に記載の発明は、請求項１に記載のＸ線ＣＴ装置において、前記マーカ部材は、前記Ｘ線検出器に対する前記平板の設置角度が変化しないように固定され、前記焦点移動量算出部は、前記平板の設置角度情報を利用して、前記Ｘ線源の焦点の移動量を３次元的に算出する。
請求項３に記載の発明は、請求項１に記載のＸ線ＣＴ装置において、前記マーカ部材は、前記平板を前記Ｘ線源のＸ線照射口に近接させて配置される。
請求項４に記載の発明は、請求項１に記載のＸ線ＣＴ装置において、前記マーカ特徴点検出部は、前記マーカの投影像を画像処理して輝度重心点を求め、前記輝度重心点を前記特徴点の座標とする。
請求項５に記載の発明は、請求項１に記載のＸ線ＣＴ装置において、前記マーカ部材における前記平板は、Ｘ線を減衰させる材料から成る。
請求項６に記載の発明は、請求項５に記載のＸ線ＣＴ装置において、前記マーカは、前記平板内に設けられた空洞、あるいは、前記平板に形成された貫通孔または凹部である。
請求項７に記載の発明は、請求項５に記載のＸ線ＣＴ装置において、前記マーカは、前記平板の表面に配設されたＸ線を減衰させる材料から成る円柱または円錐台形状の部材である。
請求項８に記載の発明は、請求項１に記載のＸ線ＣＴ装置において、前記焦点移動量算出部は、Ｘ線焦点が基準位置にあるときの前記特徴点の座標と、ＣＴスキャンの実行中の前記特徴点の座標とを用いて、前記Ｘ線源の焦点がＣＴスキャンの実行中に移動した量を３次元的に算出する。
請求項９に記載の発明は、請求項１に記載のＸ線ＣＴ装置において、前記焦点移動量算出部は、各マーカの位置関係と各マーカの特徴点の移動量から、前記Ｘ線源の焦点のＸ線光軸に沿う方向の移動量を算出する。
請求項１０に記載の発明は、請求項１に記載のＸ線ＣＴ装置において、前記平板が、前記Ｘ線源の筐体とは独立して定位置に固定されている。
請求項１から請求項１０に記載の発明によれば、Ｘ線焦点位置を検出するためのマーカを少なくとも２カ所設けた平板を有するマーカ部材を備え、少なくとも２つのマーカを被写体の投影像に重畳させることなく同時にスキャンして、ＣＴスキャン中のＸ線焦点位置を３次元的に検出し、Ｘ線源の焦点移動量に基づいて断層画像を再構成する際の座標系を補正することから、断層画像を空間歪なく再構成することが可能となる。このため、３次元形状測定用として高い寸法測定精度を実現することができる。
請求項２に記載の発明によれば、Ｘ線検出器に対する平板の設置角度を固定し、設置角度をＸ線の焦点移動量の算出パラメータとして加えることにより、ＣＴスキャン中のＸ線焦点の位置を、より正確に求めることが可能となる。
請求項３に記載の発明によれば、マーカ部材は、平板をＸ線源のＸ線照射口に近接させて配置されることから、投影像におけるマーカの拡大率を大きくすることができる。これにより、平板に設けるマーカの大きさを小さくすることが可能となる。
請求項４に記載の発明によれば、マーカの投影像を画像処理することにより、輝度重心点を特徴点の座標とすることから、マーカの投影像が変形している場合でも、特徴点の座標を正確に決めることができる。
請求項５から請求項７に記載の発明によれば、平板の厚みと異なる厚みを持つマーカを設けることで、投影像において平板の板部分とマーカ部分とでＸ線透過強度に明確な差を生じさせ、画像処理による特徴点の検出を容易に行うことができる。
請求項６に記載の発明によれば、平板にマーカを容易に形成することができる。
この発明に係るＸ線ＣＴ装置の概要図である。 Ｘ線検出器１２の検出面におけるマーカＭの位置を説明する概要図である。 この発明に係るＸ線ＣＴ装置の主要な制御系を説明するブロック図である。 マーカＭの投影像の取得からＣＴ画像再構成に至る処理を示すフローチャートである。
この発明の実施の形態を、図面に基づいて説明する。図１は、この発明に係るＸ線ＣＴ装置の概要図である。図２は、Ｘ線検出器１２の検出面におけるマーカＭの位置を説明する概要図である。図３は、この発明に係るＸ線ＣＴ装置の主要な制御系を説明するブロック図である。
Ｘ線ＣＴ装置は、Ｘ線を発生させるＸ線源１１と、Ｘ線源１１に対向配置されたＸ線検出器１２と、被写体Ｗを載置するための回転ステージ１３と、Ｘ線ＣＴ装置全体の制御を行う制御装置３０を備える。制御装置３０は、パーソナルコンピュータ４０と、Ｘ線コントローラ３１と、ステージコントローラ３２とから構成される。
Ｘ線源１１とＸ線検出器１２は、図示を省略した支持機構により定盤１６に固定されている。Ｘ線源１１とＸ線検出器１２との間には、回転軸Ｒを中心に回転する回転ステージ１３が配置されている。回転ステージ１３は、ステージ移動機構１５により定盤１６上をＸ線光軸に沿う方向に移動する。回転ステージ１３とＸ線源１１との距離をステージ移動機構１５により変更することで、Ｘ線検出器１２で検出される被写体Ｗの投影像の拡大率が変更される。なおここで、本明細書では、図２に記載されているように、Ｘ軸はＸ線光軸に沿う水平方向であり、Ｚ軸は上下（鉛直）方向であり、Ｙ軸はＸ軸とＺ軸に直交する方向である。
Ｘ線源１１は、Ｘ線検出器１２に向けてＸ線を円錐状に照射するＸ線管を有する。被写体Ｗの材質やＸ線透過特性に応じて、Ｘ線管に供給する管電圧、管電流をＸ線コントローラ３１により制御する。Ｘ線コントローラ３１は、パーソナルコンピュータ４０の制御下に置かれている。Ｘ線検出器１２により検出された透過Ｘ線像は、パーソナルコンピュータ４０に取り込まれる。パーソナルコンピュータ４０では、投影データに基づいてＣＴ画像再構成部５３（図３参照）により回転軸Ｒに直交するＸ−Ｙ平面に沿った面でスライスした断層画像が構築される。
回転ステージ１３およびステージ移動機構１５は、それぞれ独立の駆動モータを有し、それらの駆動モータは、ステージコントローラ３２から供給される駆動信号によって制御される。ステージコントローラ３２は、パーソナルコンピュータ４０の制御下に置かれている。
また、Ｘ線源１１と回転ステージ１３との間には、マーカＭが形成された平板２１と、平板２１を支持する支持部２２から成るマーカ部材２０が配置されている。マーカ部材２０は、回転ステージ１３よりもＸ線源１１に近い位置に配置される。
平板２１は、支持部２２で支持した状態で変形することがない厚みを持つ薄板により作成される。そして、平板２１には、２カ所にマーカＭが設けられている。この実施形態では、マーカＭは、矩形状の平板２１の左右対称な位置に穿設された貫通孔である。なお、貫通孔に代えて、平板２１の材料と同様のＸ線を減衰させる材料により作成された円柱や円錐台形状の部材を平板２１の表面に配設してマーカＭとしてもよい。また、貫通孔の形状も、円柱孔に限らず、孔の内面が傾斜しているテーパ孔であってもよい。さらに、貫通孔ではなく、平板２１に凹部を形成することでマーカＭとしてもよく、平板２１の材質によっては、平板内に空洞（ボイド）を設けることでマーカＭとしてもよい。すなわち、マーカＭは、平板２１の厚みと異なる厚みを持つことで、投影像において平板２１とのコントラストを得ることができるものであればよく、マーカＭの形状は、平板２１の板部分とマーカ部分とでＸ線透過強度に明確な差を生じさせるものであればよい。また、この実施形態では、マーカＭを２カ所設けているが、マーカＭの数をさらに増やしてもよい。
平板２１は、支持部２２に固定されており、設置角度が変化しないようになっている。Ｘ線源１１の焦点からＸ線を照射して被写体Ｗの投影データを取得する際には、図２に示すように、マーカＭがＸ線検出器１２の検出範囲内に含まれるように、マーカ部材２０の位置をＸ線源１１に近い側に固定する。これにより、Ｘ線源１１のＸ線照射口の前の近接した位置に平板２１が固定される。なお、平板２１の位置を規定する支持部２２の長さは、図２に示すように、マーカＭを含む平板２１の投影像が被写体Ｗの投影像と重畳することがないＸ線検出器１２の検出範囲の側端の位置となる長さとなっている。このように、支持部２２を介しマーカ部材２０を定盤１６に配置することで、マーカＭが形成された平板２１を、Ｘ線源１１の筐体とは独立して定位置に固定することができる。
平板２１におけるマーカＭの形成位置は、マーカＭのいずれもが被写体の投影像と重畳しない範囲でマーカＭ同士の離間距離が長くなるようにすることで、後述するＸ線焦点の３次元的な移動量の算出精度がよくなる。
Ｘ線焦点から円錐状に照射され、平板２１および回転ステージ１３上の被写体Ｗを透過したＸ線は、Ｘ線検出器１２により検出され、検出データはパーソナルコンピュータ４０に入力される。
パーソナルコンピュータ４０は、ＲＡＭ、ＲＯＭなどのメモリ４１と、各種の演算処理を実行するＣＰＵなどの演算装置４２と、ＣＴスキャンにより取得された投影データを保存するＨＤＤなどの記憶装置４３と、Ｘ線コントローラ３１やステージコントローラ３２に制御信号を送信するための通信部４４を備える。メモリ４１、演算装置４２、記憶装置４３および通信部４４は、相互のデータ通信を可能とする内部バス６１により接続されている。メモリ４１には、演算装置４２を動作させて機能を実現するプログラムが格納されている。図３においては、パーソナルコンピュータ４０にインストールされているプログラムを機能ブロックとして示している。この実施形態では、機能ブロックとして、マーカ特徴点検出部５１と、焦点移動量算出部５２と、ＣＴ画像再構成部５３とを備える。また、パーソナルコンピュータ４０には、演算装置４２の動作により構築された断層画像などを表示する表示装置４８と、各種指令を装置に与えるためのマウスやキーボードなどの入力装置４９が、通信部４４を介して接続されている。
次に、上述した構成のＸ線ＣＴ装置によりＣＴスキャンを実行する際に、スキャン中のＸ線の焦点位置を検出する手順について説明する。断層画像は、ＣＴスキャン時のＸ線焦点、被写体Ｗ、検出器の位置関係を含む幾何情報を基に、投影像を逆投影して再構成される。したがって、この発明では、断層画像の再構成に必要な幾何情報の一つであるＸ線の焦点位置の移動を追跡し、移動量を断層画像の再構成に反映するようにしている。図４は、マーカＭの投影像の取得からＣＴ画像再構成に至る処理を示すフローチャートである。
まず、マーカＭが形成された平板２１がＸ線源１１のＸ線照射口の直前に配置されるように、マーカ部材２０を定盤１６に固定する（ステップＳ１１）。このとき、Ｘ線検出器１２の検出面と平板２１の板面とが平行になるように、マーカ部材２０を固定する。これにより、マーカＭの位置、平板２１のＸ線検出器１２の検出面に対する設置角度が固定される。そして、Ｘ線焦点の基準となるＸ線焦点位置（基準焦点位置）を決める（ステップＳ１２）。なお、Ｘ線焦点位置は、Ｘ線源１１におけるＸ線管において、熱電子がターゲットに衝突してＸ線を発生させる位置であって、この明細書では一つの点であると見なしている。
基準焦点位置でのマーカＭの投影像を取得し（ステップＳ１３）、得られた投影像からマーカＭの特徴点を検出する（ステップＳ１４）。すなわち、被写体ＷへのＣＴスキャンを開始する前に、Ｘ線管の熱変形が起こる前の基準となるＸ線焦点位置でのマーカＭの特徴点（基準マーカ特徴点）の位置を取得する。
マーカＭの特徴点は、Ｘ線焦点とマーカＭを通る直線がＸ線検出器１２の検出面と交差する点である。Ｘ線検出器１２より得られた投影像を画像処理して決定したマーカＭの投影像の輝度重心点を、特徴点の座標としている。なお、マーカＭの投影像の輝度重心点は、２次元Ｘ線画像上の点であり、マーカＭの特徴点の移動は、このＸ線検出器１２の検出面であるｕ−ｖ平面座標におけるベクトルで表される。また、マーカＭの特徴点の座標を輝度重心点とすることで、マーカＭへのＸ線の入射角との関係で、平板２１に形成された円形のマーカＭが投影像では楕円形に変形しても、マーカＭの特徴点の座標を正確に求めることができる。
続いて、被写体ＷとマーカＭとを同時に撮像しながらスキャンする（ステップＳ１５）。ＣＴスキャンに際しては、Ｘ線源１１からＸ線を照射しながら、被写体Ｗを載置した回転ステージ１３を回転させて、所定の回転角度ごとにＸ線検出器１２からパーソナルコンピュータ４０に入力された投影像が、フレームデータとして記憶装置４３に記憶される。なお、フレームデータとは投影像の１画面分のデータのことであり、得られた時刻の順番にしたがって記憶装置４３に記憶されるものである。
得られた投影像を画像処理することにより、マーカＭの特徴点を検出する（ステップＳ１６）。ＣＴスキャンでは、設定されたスキャンモードで設定されたビュー数に応じたサンプリングピッチに対応する投影データが収集される。スキャン中は、被写体Ｗには複数の角度からＸ線が照射されるが、Ｘ線検出器１２に対する平板２１の位置は固定されていることから、各フレームから、マーカＭの投影像が得られる。スキャンの経過に伴い、Ｘ線管の熱変形などによりＸ線焦点位置が移動すると、それは、Ｘ線検出器１２の検出面におけるｕ−ｖ平面座標上でのマーカＭの特徴点の移動として観察できる。フレームデータに対して画像処理を実行することで、異なる２つのフレームのマーカＭの特徴点を検出する（ステップＳ１６）。なお、マーカＭの特徴点の検出は、全てのフレームデータに対して実行してもよく、所定の間隔で抽出したフレームデータに対して行ってもよい。すなわち、後述するステップＳ１７で、２つの異なるフレーム間で、どれだけマーカＭの特徴点が移動したかを計算することができればよい。
なお、ステップＳ１４およびステップＳ１６におけるマーカＭの特徴点の検出は、演算装置４２がメモリ４１におけるマーカ特徴点検出部５１から読み込んだプログラムを実行することにより実現される。
ステップＳ１４で得られたＸ線焦点が基準位置にあるときのマーカＭの特徴点のベクトルと、ステップＳ１６で得られたＣＴスキャン中のあるフレームのマーカＭの特徴点のベクトルから、特徴点がどれだけどの方向に移動したかを求める（ステップＳ１７）。すなわち、ある時間経過後のマーカＭの特徴点の座標位置の変化から、その時間内に、Ｘ線焦点がどれだけの量どの方向に移動したかが算出される。このＸ線焦点移動量の算出は、演算装置４２がメモリ４１における焦点移動量算出部５２から読み込んだプログラムを実行することにより実現される。
なお、この実施形態では、２つのマーカＭの特徴点の位置を求めており、２つのマーカＭの位置関係と各マーカＭの特徴点の移動量から、Ｘ線焦点のＸ線光軸に沿う方向（Ｘ方向）の移動も計算することができる。このように、少なくとも２つのマーカＭの２つの異なるフレームで検出した特徴点の座標を用いて、Ｘ線の焦点が当該時間内で移動した量を、移動した方向も含めて３次元的に算出することができる。
このＸ線ＣＴ装置では、ＣＴスキャン中のＸ線検出器１２の検出面に対して平板２１の板面をほぼ平行に設置している。このため、Ｘ線焦点の移動量を算出する際に、特に設置角度情報をパラメータとして設定しなくても、ＣＴスキャン中のＸ線焦点の移動量を算出することができる。
一方、Ｘ線焦点の移動量の３次元的算出において、Ｘ線検出器１２の検出面に対する平板２１の設置角度をパラメータとして加えることもできる。Ｘ線検出器１２の検出面に対する平板２１の設置角度が平行からずれるほど、２つのマーカＭの各々とＸ線検出器１２の検出面との距離に差が生じ、投射像における両者の拡大率が異なることになる。Ｘ線焦点の空間移動の検出に、既知パラメータとして予め取得しておいた設置角度情報を利用することで、Ｘ線焦点移動の検出精度をさらに向上させることができる。
Ｘ線検出器１２の検出面に対する平板２１の設置角度は、例えば、マーカ部材２０を定盤１６に固定したときに、レーザ変位計などを使用して取得しておけばよい。なお、このＸ線ＣＴ装置のようにＸ線源１１とＸ線検出器１２の位置が固定されており、ＣＴスキャン中のマーカ部材２０の位置が定盤１６に固定されているに等しい状態の装置では、Ｘ線検出器１２とマーカ部材２０を配置して以降、Ｘ線ＣＴ装置の使用時においてこの設置角度が変わることがないため、レーザ変位計などを使用した設置角度情報の取得は一度で済むことになる。
Ｘ線の焦点移動量が求まると、その焦点移動量に基づいて、ＣＴ画像再構成の座標系を修正する（ステップＳ１８）。Ｘ線の焦点位置は、断層画像の再構成に必要な幾何情報の一つであり、焦点移動量に基づいてＣＴ画像再構成の座標系を修正することで、ＣＴスキャン中の幾何情報の誤差が修正される。しかる後、ＣＴ画像再構成を実行し（ステップＳ１９）、被写体Ｗの断層画像を再構成する。なお、ＣＴ画像再構成の座標系の修正（ステップＳ１８）およびＣＴ画像再構成の実行（ステップＳ１９）は、演算装置４２がメモリ４１におけるＣＴ画像再構成部５３から読み込んだプログラムを実行することにより実現される。
上述したように、この発明のＸ線ＣＴ装置では、少なくとも２つのマーカＭを被写体Ｗと同時に撮像しながらＣＴスキャンし、被写体Ｗの投影データと同じ時間軸でマーカＭの特徴点を検出することで、ＣＴスキャン中のＸ線焦点位置を３次元的に検出している。これにより、ＣＴスキャン中のＸ線焦点の移動を断層画像の再構成に反映できることから、断層画像を空間歪みなく再構成すること可能となる。したがって、この発明のＸ線ＣＴ装置を高い寸法測定精度が要求される３次元形状測定用に使用することも可能となる。
１１ Ｘ線源
１２ Ｘ線検出器
１３ 回転ステージ
１５ ステージ移動機構
１６ 定盤
２０ マーカ部材
２１ 平板
２２ 支持部
３０ 制御装置
３１ Ｘ線コントローラ
３２ ステージコントローラ
４０ パーソナルコンピュータ
４１ メモリ
４２ 演算装置
４３ 記憶装置
４４ 通信部
４８ 表示装置
４９ 入力装置
５１ マーカ特徴点検出部
５２ 焦点移動量算出部
５３ ＣＴ画像再構成部
Ｍ マーカ
Ｗ 被写体

Claims (7)

  1. 被写体に複数の角度からＸ線を照射して取得した投影データに基づいて前記被写体の断層画像を再構成するＸ線ＣＴ装置であって、
    Ｘ線を発生するＸ線源と、
    前記Ｘ線源に対向配置されるＸ線を検出するためのＸ線検出器と、
    前記Ｘ線源と前記Ｘ線検出器との間に配置される前記被写体を載置するための回転ステージと、
    前記Ｘ線検出器が検出した前記被写体の投影データに基づいて演算処理を実行する制御装置と、
    マーカが少なくとも２か所設けられた平板と、前記平板を支持する支持部とから成り、前記回転ステージと前記Ｘ線源との間の位置であって、ＣＴスキャンの実行中に前記マーカが前記Ｘ線検出器の検出範囲内に含まれる位置、かつ、前記平板が前記被写体の投影像に重畳しない位置に配置されるマーカ部材と、
    を備え、
    前記制御装置は、
    前記Ｘ線源の焦点と前記マーカとを通る直線が前記Ｘ線検出器と交差する点である特徴点を、前記マーカの投影像を画像処理することにより検出し、２次元Ｘ線画像上の前記特徴点の座標を求めるマーカ特徴点検出部と、
    前記マーカ特徴点検出部により求めた２つの異なるフレームの前記特徴点の座標を用いて、前記Ｘ線源の焦点が前記２つの異なるフレームが得られる間に移動した量を３次元的に算出する焦点移動量算出部と、
    を備え、
    前記Ｘ線源の焦点移動量に基づいて前記被写体の断層画像を再構成するときの座標系を修正することを特徴とするＸ線ＣＴ装置。
  2. 請求項１に記載のＸ線ＣＴ装置において、
    前記マーカ部材は、前記Ｘ線検出器に対する前記平板の設置角度が変化しないように固定され、
    前記焦点移動量算出部は、前記平板の設置角度情報を利用して、前記Ｘ線源の焦点の移動量を３次元的に算出するＸ線ＣＴ装置。
  3. 請求項１に記載のＸ線ＣＴ装置において、
    前記マーカ部材は、前記平板を前記Ｘ線源のＸ線照射口に近接させて配置されるＸ線ＣＴ装置。
  4. 請求項１に記載のＸ線ＣＴ装置において、
    前記マーカ特徴点検出部は、前記マーカの投影像を画像処理して輝度重心点を求め、前記輝度重心点を前記特徴点の座標とするＸ線ＣＴ装置。
  5. 請求項１に記載のＸ線ＣＴ装置において、
    前記マーカ部材における前記平板は、Ｘ線を減衰させる材料から成るＸ線ＣＴ装置。
  6. 請求項５に記載のＸ線ＣＴ装置において、
    前記マーカは、前記平板内に設けられた空洞、あるいは、前記平板に形成された貫通孔または凹部であるＸ線ＣＴ装置。
  7. 請求項５に記載のＸ線ＣＴ装置において、
    前記マーカは、前記平板の表面に配設されたＸ線を減衰させる材料から成る円柱または円錐台形状の部材であるＸ線ＣＴ装置。
JP2019528189A 2017-07-03 2017-07-03 Ｘ線ｃｔ装置 Active JP7164524B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2021187680A JP7251602B2 (ja) 2017-07-03 2021-11-18 Ｘ線ｃｔ装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2017/024319 WO2019008620A1 (ja) 2017-07-03 2017-07-03 Ｘ線ｃｔ装置

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2021187680A Division JP7251602B2 (ja) 2017-07-03 2021-11-18 Ｘ線ｃｔ装置

Publications (2)

Publication Number Publication Date
JPWO2019008620A1 true JPWO2019008620A1 (ja) 2020-02-27
JP7164524B2 JP7164524B2 (ja) 2022-11-01

Family

ID=64949871

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019528189A Active JP7164524B2 (ja) 2017-07-03 2017-07-03 Ｘ線ｃｔ装置

Country Status (6)

Country Link
US (1) US11002690B2 (ja)
EP (1) EP3620778A4 (ja)
JP (1) JP7164524B2 (ja)
CN (1) CN110621985B (ja)
TW (1) TWI680293B (ja)
WO (1) WO2019008620A1 (ja)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3690428A4 (en) * 2017-09-28 2021-06-09 Saki Corporation INSPECTION POSITION SPECIFICATION PROCESS, THREE-DIMENSIONAL IMAGE GENERATION PROCESS AND INSPECTION DEVICE
JP7143567B2 (ja) * 2018-09-14 2022-09-29 株式会社島津テクノリサーチ 材料試験機および放射線ｃｔ装置
EP3932317A4 (en) 2019-02-26 2022-10-05 National University Corporation Shizuoka University X-RAY IMAGING DEVICE
JP2020186913A (ja) * 2019-05-09 2020-11-19 株式会社リガク Ｘ線ｃｔ装置、及びｃｔ画像再構成方法
CN111449668B (zh) * 2020-04-23 2023-06-23 深圳市安健科技股份有限公司 三维扫描重建中实时几何校正的标记装置、方法及系统
CN111812130A (zh) * 2020-07-15 2020-10-23 深圳市金园智能科技有限公司 一种基于x射线的材料内部3d成像方法及装置
CN114280086B (zh) * 2021-11-16 2024-01-23 中国电子科技集团公司第三十八研究所 一种ct成像装置
CN114199905B (zh) * 2021-12-13 2024-02-20 中国航发南方工业有限公司 一种机匣内部缺陷的空间定位方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000298105A (ja) * 1999-04-14 2000-10-24 Toshiba Fa Syst Eng Corp コンピュータ断層撮影装置
JP2002333408A (ja) * 2001-05-08 2002-11-22 Hitachi Medical Corp 産業用ｘ線ｃｔ装置
JP2004340630A (ja) * 2003-05-13 2004-12-02 Sony Corp コンピュータ断層撮像方法及び装置
JP2005270297A (ja) * 2004-03-24 2005-10-06 Canon Inc 放射線ｃｔ撮影装置及び放射線ｃｔ撮影システム及びそれを用いた放射線ｃｔ撮影方法
JP2007255951A (ja) * 2006-03-20 2007-10-04 Ratoc System Engineering Co Ltd Ｘ線検査装置
JP2017086816A (ja) * 2015-11-17 2017-05-25 東芝メディカルシステムズ株式会社 Ｘ線コンピュータ断層撮影装置

Family Cites Families (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4054402B2 (ja) 1997-04-25 2008-02-27 株式会社東芝 Ｘ線断層撮影装置
JPS5850059B2 (ja) 1977-02-04 1983-11-08 日本電気株式会社 受信信号制御回路
JPS62284250A (ja) * 1986-05-31 1987-12-10 Toshiba Corp 産業用ｃｔスキヤナ
JP3743594B2 (ja) 1998-03-11 2006-02-08 株式会社モリタ製作所 Ｃｔ撮影装置
US6960020B2 (en) * 2001-08-31 2005-11-01 Analogic Corporation Image positioning method and system for tomosynthesis in a digital X-ray radiography system
DE10206190A1 (de) * 2002-02-14 2003-09-04 Siemens Ag Verfahren und Vorrichtung zur Erzeugung eines Volumendatensatzes
US7310404B2 (en) 2004-03-24 2007-12-18 Canon Kabushiki Kaisha Radiation CT radiographing device, radiation CT radiographing system, and radiation CT radiographing method using the same
JP4375555B2 (ja) 2004-05-14 2009-12-02 株式会社島津製作所 Ｘ線ｃｔ装置
US7344307B2 (en) * 2004-11-12 2008-03-18 General Electric Company System and method for integration of a calibration target into a C-arm
DE102005021020A1 (de) * 2005-05-06 2006-11-16 Siemens Ag Verfahren zur Berechnung einer Orthogonal-Röntgenabschwächung eines Objekts anhand einer gemessenen Referenz-Röntgenabschwächung sowie zugehörige Vorrichtung
WO2008038283A2 (en) * 2006-09-25 2008-04-03 Mazor Surgical Technologies Ltd. C-arm computerized tomography system
US7835811B2 (en) * 2006-10-07 2010-11-16 Voxelogix Corporation Surgical guides and methods for positioning artificial teeth and dental implants
EP2091459B1 (en) * 2006-10-23 2013-09-11 Hirdesh Sahni An image guided whole body stereotactic needle placement device with falling arc
DE102007043820B4 (de) 2007-09-13 2020-06-04 Carl Zeiss Industrielle Messtechnik Gmbh Verfahren zur Bestimmung eines Korrekturwerts einer Bremsfleckposition einer Röntgenquelle einer Messanordnung sowie eine Messanordnung zum Erzeugen von Durchstrahlungsbildern
CN101398397B (zh) * 2007-09-30 2012-05-30 首都师范大学 多次扫描模式的ct成像方法
CN102123664B (zh) * 2008-08-13 2014-05-07 皇家飞利浦电子股份有限公司 使用旋转中心寻找算法进行环形伪影校正的校准方法
WO2010064287A1 (ja) * 2008-12-01 2010-06-10 株式会社 島津製作所 放射線撮像装置
JP5400546B2 (ja) * 2009-09-28 2014-01-29 株式会社日立メディコ Ｘ線ｃｔ装置
US9408579B2 (en) * 2010-05-11 2016-08-09 Takara Telesystems Corp. Radiation imaging apparatus and phantom used for the same
CN102652674B (zh) * 2011-03-04 2014-02-19 首都师范大学 一种消除ct图像中的几何伪影的方法和系统
US9693748B2 (en) * 2011-07-23 2017-07-04 Broncus Medical Inc. System and method for automatically determining calibration parameters of a fluoroscope
DE102011083416A1 (de) * 2011-09-26 2013-03-28 CT Imaging GmbH Röntgengerät
EP2765408B1 (en) 2011-10-04 2018-07-25 Nikon Corporation X-ray device, x-ray irradiation method, and manufacturing method for structure
JP6032729B2 (ja) 2012-05-08 2016-11-30 国立研究開発法人理化学研究所 イメージングマーカーおよびその利用
CN106030293B (zh) * 2014-01-23 2019-11-26 株式会社蛟簿 X射线检查装置以及x射线检查方法
JP2015215218A (ja) 2014-05-09 2015-12-03 株式会社島津製作所 放射線検査装置
JP6394082B2 (ja) 2014-06-09 2018-09-26 株式会社島津製作所 Ｘ線検査装置
KR102144273B1 (ko) * 2014-10-14 2020-08-13 가부시키가이샤 리가쿠 X선 박막 검사 장치
CN106361367B (zh) * 2016-12-01 2019-10-08 上海联影医疗科技有限公司 一种检测器的校正方法和使用该校正方法的装置及设备

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000298105A (ja) * 1999-04-14 2000-10-24 Toshiba Fa Syst Eng Corp コンピュータ断層撮影装置
JP2002333408A (ja) * 2001-05-08 2002-11-22 Hitachi Medical Corp 産業用ｘ線ｃｔ装置
JP2004340630A (ja) * 2003-05-13 2004-12-02 Sony Corp コンピュータ断層撮像方法及び装置
JP2005270297A (ja) * 2004-03-24 2005-10-06 Canon Inc 放射線ｃｔ撮影装置及び放射線ｃｔ撮影システム及びそれを用いた放射線ｃｔ撮影方法
JP2007255951A (ja) * 2006-03-20 2007-10-04 Ratoc System Engineering Co Ltd Ｘ線検査装置
JP2017086816A (ja) * 2015-11-17 2017-05-25 東芝メディカルシステムズ株式会社 Ｘ線コンピュータ断層撮影装置

Also Published As

Publication number Publication date
JP7164524B2 (ja) 2022-11-01
TWI680293B (zh) 2019-12-21
WO2019008620A1 (ja) 2019-01-10
CN110621985A (zh) 2019-12-27
US11002690B2 (en) 2021-05-11
EP3620778A4 (en) 2020-08-19
EP3620778A1 (en) 2020-03-11
CN110621985B (zh) 2022-03-11
US20200124545A1 (en) 2020-04-23
TW201907154A (zh) 2019-02-16

Similar Documents

Publication Publication Date Title
TWI680293B (zh) X線電腦斷層裝置
JP2019128163A (ja) 計測用ｘ線ｃｔ装置、及び、その校正方法
US20030043962A1 (en) Image positioning method and system for tomosynthesis in a digital X-ray radiography system
JP5125423B2 (ja) Ｘ線断層画像によるはんだ電極の検査方法およびこの方法を用いた基板検査装置
JP2009505083A (ja) コンピュータ断層撮影用測定装置および方法
JP4561990B2 (ja) Ｘ線撮影装置
US20070114428A1 (en) Radiation Image Capturing Condition Correction Apparatus and Radiation Image Capturing Condition Correction Method
US10729390B2 (en) Rotary collimator for determining the position of an element provided with sensors in an X-ray imaging system
JP2012112790A (ja) Ｘ線ｃｔ装置
JP6394082B2 (ja) Ｘ線検査装置
JPWO2018083930A1 (ja) 放射線断層撮影装置の撮像倍率校正方法
CN117871557A (zh) 检查装置、以及检查区域的确定方法
JP5011085B2 (ja) 検査装置，検査装置のアライメント方法及び校正用ファントムの製作方法
JP5487519B2 (ja) 産業用ｘ線ｃｔ装置および撮像方法
US20160183892A1 (en) Method for imaging by means of an x-ray device and x-ray device
JP2012042340A (ja) Ｘ線ｃｔ装置
JP7251602B2 (ja) Ｘ線ｃｔ装置
JP2021050937A (ja) 計測用ｘ線ｃｔ装置の校正方法、測定方法、及び、計測用ｘ線ｃｔ装置
JP2009276142A (ja) 放射線検査システム及び放射線検査の撮像方法
JP4479503B2 (ja) 断層撮影装置
JP2013061257A (ja) Ｘ線検査装置、ｘ線検査装置の制御方法、ｘ線検査装置を制御するためのプログラム、および、当該プログラムを格納した記録媒体
JP2005134213A (ja) Ｘ線断層撮像方法及び装置
JP2006292462A (ja) コンピュータ断層撮影装置
JP2005195530A (ja) 校正器、及び、校正方法
JP2012120564A (ja) Ｘ線撮影装置

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20191015

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20191015

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20200220

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20201013

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20201202

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210608

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20211019

C60 Trial request (containing other claim documents, opposition documents)

Free format text: JAPANESE INTERMEDIATE CODE: C60

Effective date: 20211118

C22 Notice of designation (change) of administrative judge

Free format text: JAPANESE INTERMEDIATE CODE: C22

Effective date: 20220111

C13 Notice of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: C13

Effective date: 20220628

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220720

C23 Notice of termination of proceedings

Free format text: JAPANESE INTERMEDIATE CODE: C23

Effective date: 20220906

C03 Trial/appeal decision taken

Free format text: JAPANESE INTERMEDIATE CODE: C03

Effective date: 20221004

C30A Notification sent

Free format text: JAPANESE INTERMEDIATE CODE: C3012

Effective date: 20221004

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20221020

R151 Written notification of patent or utility model registration

Ref document number: 7164524

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151