JP3964155B2 - X-ray image correction method - Google Patents

X-ray image correction method Download PDF

Info

Publication number
JP3964155B2
JP3964155B2 JP2001150181A JP2001150181A JP3964155B2 JP 3964155 B2 JP3964155 B2 JP 3964155B2 JP 2001150181 A JP2001150181 A JP 2001150181A JP 2001150181 A JP2001150181 A JP 2001150181A JP 3964155 B2 JP3964155 B2 JP 3964155B2
Authority
JP
Japan
Prior art keywords
ray
image
center axis
imaging
image correction
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.)
Expired - Fee Related
Application number
JP2001150181A
Other languages
Japanese (ja)
Other versions
JP2002336237A (en
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.)
J Morita Manufaturing Corp
Original Assignee
J Morita Manufaturing 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 J Morita Manufaturing Corp filed Critical J Morita Manufaturing Corp
Priority to JP2001150181A priority Critical patent/JP3964155B2/en
Publication of JP2002336237A publication Critical patent/JP2002336237A/en
Application granted granted Critical
Publication of JP3964155B2 publication Critical patent/JP3964155B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/58Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
    • A61B6/582Calibration
    • A61B6/583Calibration using calibration phantoms

Description

【0001】
【発明の属する技術分野】
本発明は、コーンビームX線CT撮影、つまり、被写体を固定保持し、この被写体を中心としてX線発生器と二次元X線撮像手段とを対向配置し、前記X線発生器から前記被写体の一部である撮影すべき局所部位のみを包含するX線コーンビームを局所照射しながら、前記X線発生器と二次元X線撮像手段を撮影条件に応じた角度範囲で旋回させることによって、前記局所部位のX線画像を撮像し、このX線画像を演算処理して、前記局所部位の3次元的なX線吸収係数を算出するX線撮影において、被写体と前記被写体のX線原画像の位置的なずれを補正するX線画像補正方法に関する。
【0002】
【従来の技術】
X線CT(Computed Tomography)撮影においては、二次元X線撮像手段として、イメージインテンシファイア略してX線II(X線蛍光増倍管)とCCDカメラを組み合わせたもの(以後、「XII撮像管」という。)が用いられるが、このXII撮像管の場合、X線受光面が曲面であるのに対し、ここで可視化された光を電気信号に変換するCCDカメラの受光面は平面であるため、曲面を平面で検知することにより生じる管面歪みや、撮影中の地磁気による磁気歪みなどの歪み、つまり、被写体と前記被写体のX線原画像の位置的なずれが発生するが、これを補正するものとして種々の提案がなされていた。
【0003】
図10(a)はこの管面歪みの補正に用いる従来の画像補正ファントムの一例を示す図、(b)は、これを用いたX線CT撮影装置の基本構成図であり、特開平8−24248号公報に記載されたものである。
【0004】
図10(a)の画像補正ファントム121は、X線不透過性の金属平板161aに正方格子状にX線透過性の孔121bを形成したものであり、管面歪みの補正に好適に用いることができたが、金属平板121aへの孔121bの加工に多数の手間を要し、コストの高いものとなっていた。
【0005】
図10(b)に示すX線CT撮影装置は、被写体Oを挟んで対向配置させたX線発生器101と二次元X線撮像手段102が回転板103上に設置され、これを回転することで、X線発生器101から被写体OにX線コーンビームを旋回照射して、二次元X線撮像手段102で被写体OのX線画像を得、画像収集処理装置106において、得られたX線画像を位置対応テーブルなどを用いて補正などを行い、補正後のX線画像を再構成して、被写体Oの3次元的なX線吸収係数(再構成データ)を得て、これを画像表示装置Eによって表示するものであったが、この際、画像補正ファントム121は、二次元X線撮像手段102の管面に装着されて画像補正に用いられていた。
【0006】
また、被写体とこの被写体のX線原画像の位置的なずれの原因となるものとしては、対向配置されたX線発生器と二次元X線撮像手段とを被写体を中心として旋回させる回転中心軸の軸ぶれや、この回転中心軸と、X線発生器と二次元X線撮像手段の軸方向の位置ずれがあり、これを補正するための種々の提案もなされていた。
【0007】
図11(a)はこのような画像補正ファントムの他例を示す図、(b)は、これを用いたX線CT撮影装置の基本構成図であり、特開平9−173330号公報に記載されたものである。
【0008】
図11(a)の画像補正ファントム221は、X線透過性の支持体221aの内部にX線不透過性の複数の球体221bを上下に配列したものであり、回転中心軸の軸ぶれの補正に好適に用いることができたが、球体221bが多い分だけ、補正テーブルを作成するまでのデータ処理が煩雑なものであった。
【0009】
図11(b)に示すX線CT撮影装置は、図10(b)と同様にX線CT撮影に必須の、被写体Oを挟んで対向配置させたX線発生器201と二次元X線撮像手段202、回転板203、画像処理手段206などを備えると共に、画像補正ファントム221を用いたX線原画像の補正をするために、画像歪み補正手段206a、画像歪み補正テーブル作成手段261などを備えていた。
【0010】
図12(a)も同様の画像補正ファントムの他例を示す図、(b)は、これを用いたX線CT撮影装置の基本構成図であり、特開2000−201918号公報に記載されたものである。
【0011】
図12(a)の画像補正ファントム321は、X線透過性の支持体321aの内部にX線不透過性の棒体321bを上下方向に設置したものであり、回転中心軸と、X線発生器と二次元X線撮像手段の上下方向の位置ずれの補正に好適に用いることができたが、補正テーブルを作成するのに、このファントム321に全周照射したデータが必要であり、煩雑な処理となっていた。
【0012】
図12(b)に示すX線CT撮影装置は、図11(b)と同様にX線CT撮影に必須の、被写体Oを挟んで対向配置させたX線発生器301と二次元X線撮像手段302、回転板303、画像処理手段に相当するデータ処理部306などを備えると共に、画像補正ファントム361を用いたX線原画像の補正をするために、画像歪み補正手段に相当するジオメトリ推定部306aなどを備えていた。
【0013】
【発明が解決しようとする課題】
しかしながら、これら従来の画像補正ファントムは、所期の目的を達成するものの、なお上述したような解決すべき問題を抱えていた。
【0014】
本発明は、このような問題を解決しようとするもので、歪み補正の目的に対応して派生している種々の課題を解消する画像補正方法を提供することを目的とする。
【0015】
【課題を解決するための手段】
請求項1に記載のX線画像補正方法は、被写体を挟んで対向配置させたX線発生器と2次元撮像手段とを、回転中心軸を中心に旋回させながら、X線発生器からコーンビームX線を被写体に照射するコーンビームX線CT撮影における回転中心軸の軸方向の位置ずれを補正するX線画像補正方法であって、X線透過性の支持体内に、この支持体よりもX線吸収係数の高い2つの球体を収容して成り、この支持体を上記回転中心軸より所定位置離れた位置に設置し、上記X線発生器と上記2次元撮像手段とを上記回転中心軸を中心にして180度反転させ、それぞれの位置で上記2次元撮像手段上で得られた上記2つの球体のX線原画像に対して、所定の演算を行って上記2次元撮像手段上におけるX線原画像の基準位置を算出するようにしたことを特徴とする。
【0016】
請求項2に記載のX線画像補正方法は、請求項1において、X線CT撮影の回転中心軸の軸方向の位置ずれを補正するための中心軸ずれ補正テーブルを作成し、この中心軸ずれ補正テーブルを用いてX線CT撮影の各旋回角度毎に画像補正テーブルを作成し、該画像補正テーブルを用いて、それぞれのX線原画像を補正することを特徴とする。
【0017】
この補正方法は、請求項1に記載の補正方法を有効利用して、その効果を発揮させながら、かつ、X線CT撮影時に発生が考えられ得る歪みを補正する中心軸ずれ補正テーブルを用いて作成した画像補正テーブルによって、X線原画像を補正している。
請求項3に記載のX線画像補正方法は、請求項のX線画像補正方法において、更に、X線の支持体内にこの支持体よりもX線吸収係数の高い2つの球体を収容して成り、この支持体を、上記球体がX線CT撮影におけるほぼ回転中心軸方向に上下になるよう被写体保持位置に設置し、上記2つの球体の画像上の座標方向の動きの平均値から補正に用いる角度とシフト量とを算出することで、前記画像補正テーブルを作成することを特徴とする。
【0018】
【発明の実施の形態】
本発明の実施の形態について、添付図面を参照しながら、説明する。
【0019】
図1(a)は本発明で使用される画像補正ファントムの一例を示す正面図、(b)はこれを二次元X線撮像手段に装着した状態を示す縦断面図、(c)は(b)で撮像されたX線原画像の一例を示す図である。
【0020】
図1(a)に示す画像補正ファントム21は、X線透過性のプリント基板21aとX線吸収係数の高い小円体21bとから構成され、図示したように、プリント基板21a上に碁盤状の格子を想定し、この格子の交差点となる位置に小円体21bを配置したものである。
【0021】
プリント基板21aは、合成樹脂、例えば、ガラスエポキシ樹脂やフェノール樹脂の他、アクリル樹脂、塩化ビニール、ポリカーボネート等のプラスチック材料などX線吸収係数の低い素材で形成された板状体ものであれば良く、具体的には、電子機器の配列設置用に流布されているプリント基板であれば安価に入手可能で、好適である。
【0022】
小円体21bのX線不透過性材料としては、銅、タングステン、鉄−ニッケル−クロム合金、ステンレスなどの金属が好適であり、小円体21bの直径としては、1.5ミリ程度のものが好適であるが、これに限るものではなく、このファントムの全体寸法、つまり、ファントムを装着する二次元X線撮像手段の管面の大きさに対応して相対的に決定されるものである。
【0023】
この補正ファントム21は、補正の基準点となるX線吸収係数の高い小円体21bを格子状に配列させる板状物として、汎用のX線透過性のプリント基板21aを用いたので、基板が非常に安価なものとなる。また、格子配列は、このプリント基板21a上に、小円体21bをエッチング処理で形成するだけでよく、コストダウンとなる。
【0024】
二次元X線撮像手段2は、図1(b)に示すようにXII撮像管であり、公知のものであるが、簡単にその構成を説明すると、管面である入力窓2a、入力面基板2b、入力蛍光面2c、光電陰極2d、フォーカス電極2e、陽極2f、出力蛍光面2g、出力窓2h、光学レンズ2i、CCDカメラ2j、フォーカス電源E1、陽極電源E2などを備えたもので、X線の入射を受けて、これを可視光線化し、この可視光線をCCDカメラで撮像して、電気信号として出力する。
【0025】
この二次元X線撮像手段2には、その入力窓2aに画像補正ファントム21を装着するための取付枠25が設けられている。
【0026】
画像補正ファントム21は、図1(b)に示すように、この取付枠25によって二次元X線撮像手段2の管面に装着して用いられ、この状態で、X線照射を行い撮像したX線原画像IOが図1(c)に示したものである。
【0027】
このX線原画像IOを構成する小円体像IO(21b)と画像補正ファントム21の小円体21bの相対的な位置関係を示す補正テーブルによって、二次元X線撮像手段の管面歪みや、撮影中の地磁気による磁気歪みなどの歪み、つまり、被写体と前記被写体のX線原画像の位置的なずれを補正する。この補正テーブルを管面歪み補正テーブルという。
【0028】
図2(a)から(c)は、図1の画像補正ファントムを用いた補正方法の概念図であり、これらを用いて、この管面歪み補正テーブルの作成について説明する。
【0029】
1)X線原画像IOの2値化:これは、原画像IOにおける小円体像IO(21b)の位置を求めるためのものである。図2(a)には、図1(a)の1個の小円体21bとそれを囲む格子が示されている。ここで、格子の大きさが3ミリで、小円体21bの直径が1.5ミリであるので、小円体像IO(21b)の面積率は、
π*(1.5/2)*(1.5/2)/3*3≒0.196
となる。これがX線原画像IOのスレッシュホールドとなる。なお、ここでは、「*」は乗算、「/」は割り算を意味するものとする。
【0030】
ついで、図2(b)に示すように320画素*240画素のX線原画像IOを512画素*512画素の領域に視野境界を中心に境界内の画像データを対象に鏡像処理して外挿し、その結果にFFT(フーリエ変換)を掛け、ついで、周波数領域で図2(c)に示すようなバンドパスフィルタを掛け、更に、IFFT(逆フーリエ変換)を行い、この結果をスレッシュホールド(0.196)を用いて2値化する。
【0031】
2)補正テーブルの作成:これは原画像IOとファントム21の小円体21bの位置関係を示すものである。この問題は一般的には、図2(d)に示すように、xy座標系上の4点A,B,C,Dから、uv座標系上の4点A′,B′,C′,D′への変換問題と把握でき、ここでは、単なる座標系の回転、拡大に加え、拡大率が変化するものであるので、以下の拡張アフィン変換を用いる。
u = a11*xy + a10*x
+ a01*y + a00
v = b11*xy + b10*x
+ b01*y + b00
この場合、係数がu,vの各式毎に4個となるため、この問題のように4点相互間の関係を求めるのに適当である。
【0032】
3)小領域のアフィン変換の集合によって全体のアフィン変換を定義する。
【0033】
これは、図2(e)に示すように、左側部分の斜線部に対応する小領域のアフィン変換の集合として、全体のアフィン変換を定義するものである。
【0034】
このような手順によって、管面歪み補正テーブルが得られるが、この状態では、撮像系の傾きと上下左右の傾き、つまり、X線CT撮影の回転中心軸ぶれとX線CT撮影の回転中心軸の軸方向の位置ずれは補正されておらず、このためには、更に、以下に説明する二種類の画像補正ファントムを用いて、中心軸ぶれ補正テーブルと中心軸ずれ補正テーブルとを作成する必要がある。
【0035】
図3(a)、(b)は、画像補正ファントムの参考例を示す外観図、図4は、図3の画像補正ファントムを用いた補正方法の概念図である。
【0036】
この画像補正ファントム22は、X線透過性の支持体22aとX線吸収係数の高い2つの中球体22b(A),22b(B)とから構成され、支持体22aの内部に、2つの中球体22b(A),22b(B)を、この支持体22aを被写体保持位置に設置した際に、これらの中球体22bがX線CT撮影のほぼ回転中心軸3a方向に上下になるように配置した(図7参照)ものである。
【0037】
ここで支持体22aの素材としては、X線透過性のよい、透明な、また、入手容易で加工性の良いアクリル樹脂等がよい。中球体22bの素材は、X線吸収係数の高い鋼球などがよく、その直径は、3ミリから15ミリとするのが好適であるが、これに限られない。
【0038】
この補正ファントム22は、構成が簡単で、安いコストで作ることができると共に、このファントムを用いて補正テーブルを作成するのが容易になる。
【0039】
このファントム22は、X線CT撮影の回転中心軸ぶれを補正する中心軸ぶれ補正テーブルを作成するために用いられるが、そのために、ファントム22を図7に示すように、支持台26の上に載せ、その位置が図7のほぼ回転中心軸3a付近、つまり被写体保持位置になるようにする。この状態で、X線撮影を行い、得られたX線原画像から回転中心軸ぶれを算出する方法を図4に示す。
【0040】
また、補正ファントム22は、図3(a)のように直方体でも、(b)のように円柱体であってもよい。
【0041】
図4では、2個の中球体22b(A),22b(B)の画像を単にA、Bの符号で示している。
【0042】
まず、この図のA、Bは、管面歪み補正テーブルによって補正したものである必要がある(これは、図6についても同様である。)。ここで、角度θとシフト量XSを求めるのが目的である。
【0043】
旋回アーム3の360度の回転により、A,Bが図のように左右に動くので、簡単に、A,Bのx座標の平均値を求めることによって、点A′、B′を求めることができ、これによって、角度θとシフト量XSを算出する。より正確には、A,Bのx座標方向の動きをフーリエ変換し、それぞれの直流成分を点A′、B′のx座標とすることができる。
【0044】
角度θは小さいと仮定しているので、A,Bのy座標は一定(平均値)であり、A,Bは水平方向に動くだけと仮定しても問題ない。また、角度θとシフト量XSは、旋回アーム3の360度回転に対して1組のみ算出される。
【0045】
図5は画像補正ファントムの他例を示す外観図、図6は、図5の画像補正ファントムを用いた補正方法の概念図である。
【0046】
この画像補正ファントム23は、X線透過性の支持体23aとX線吸収係数の高い2つの小球体23b(A)、23b(B)とから構成され、支持体23aの内部に、これら2つの小球体23b(A)、23b(B)を、この支持体23aを被写体保持位置に設置した際にこれらの小球体23bがX線CT撮影のほぼ回転中心軸方向に上下になるように配置した(図7参照)ものである。
【0047】
ここで支持体23aの素材としては、上記支持体22aと同様のものでよい。小球体23bの素材としては、上記中球体22bと同様のものでよく、その直径は、1ミリから3ミリとするのが好適であるが、これに限られない。
【0048】
この補正ファントム23は、構成が簡単で、安いコストで作ることができると共に、このファントムを用いて補正テーブルを作成するのが容易になる。
【0049】
このファントム23は、X線CT撮影の回転中心軸の軸方向の位置ずれを補正するための中心軸ずれ補正テーブルを作成するために用いられるが、そのために、ファントム23を図7に示すように、支持台27の上に載せ、その位置が図7のほぼ回転中心軸3aから離れた位置、つまり被写体保持位置から所定距離離れた位置となるようにする。この状態で、X線撮影を行い、得られたX線原画像から回転中心軸の軸方向の位置ずれを算出する方法を図6に示す。
【0050】
図6では、2個の小球体23b(A),23b(B)を単にA、Bの符号で示している。また、符号1はX線発生器、符号2は二次元X線撮像手段、符号3aは旋回アーム3の回転中心軸3aである(図7参照)。
【0051】
符号a1,b1は、小球体A,Bに対してX線発生器1が最も離れた位置となる状態(一点鎖線で示す。)でこれらを照射した際の小球体A,BのX線原画像上の上下位置、符号a2,b2は、小球体A,Bに対してX線発生器1が最も近づいた位置(つまり、先の位置から回転中心軸3aを中心として180度回転させた位置)となる状態(実線で示す。)でこれらを照射した際の小球体A,BのX線原画像上の上下位置を示している。符号cはX線原画像の基準位置が回転中心軸3aに対して占める上下位置である。
【0052】
こうして得られた4つの位置a1,b1、a2,b2と、回転中心軸3aの上下位置cとの間には、
(a1 - c)/(a2 - c) =
(b1 - c)/(b2 - c)
の関係が成り立ち、これより、
c = (a2*b1 - b2*a1)/ ((a2 - b2) - (a1 - b1))
となる。
【0053】
これにより上下位置cが得られ、この上下位置cが機械上の回転中心軸3aに対するX線原画像の上下位置となるので、座標軸の上下移動式が生成でき、これが中心軸ずれ補正テーブルとなる。なお、この補正テーブルは、照射の全周について、一組用意すればよい。
【0054】
こうして、X線CT撮影の各旋回角度毎に管面歪み補正テーブル、また、全周についてそれぞれ一組の中心軸ぶれ補正テーブル、中心軸ずれ補正テーブルが得られると、管面歪み補正テーブルに対して、中心軸ぶれ補正テーブルによる座標軸回転変換、中心軸ずれ補正テーブルによる座標軸移動変換を行って、X線CT撮影の各旋回角度毎に画像補正テーブルが作成される。
【0055】
このような画像補正テーブルは、X線CT撮影時に発生が考え得る全ての歪みを補正する管面歪み補正テーブル、中心軸ぶれ補正テーブル、中心軸ずれ補正テーブルを併せ用いているので、すべての歪みを除去することができる。
【0056】
なお、中心軸ぶれ補正テーブル用の補正ファントム22の中球体22bと、中心軸ずれ補正テーブル用の補正ファントム23の小球体23bのように補正基準点として異なる大きさのものを用いたのは、それぞれの補正目的に対応したものである。
【0057】
つまり、中心軸ぶれ補正の場合は、補正ファントム22をほぼ回転中心軸付近、つまり、被写体保持位置に設置するが、この位置は、X線発生器1と二次元X線撮像手段2の真ん中であり、旋回中に、この補正ファントムの中球体22bのX線画像の大きさはあまり変化しない。そこで、二次元X線撮像手段2の撮像面にそれぞれの2個の中球体22bの画像が写る範囲で大きくすると、上述したような大きさとなる。
【0058】
一方、中心軸ずれ補正の場合は、補正ファントム23を、回転中心軸から離れた位置、被写体保持位置から所定距離離れた位置に設置する必要があるが、この場合旋回中に、この補正ファントム23の小球体23bのX線画像の大きさは大小に変化する。そこで、X線発生器1に補正ファントム23が最も近づく回転位置で、その小球体23bのX線画像の大きさが最大となった場合でも、二次元X線撮像手段2の撮像面にそれぞれの2個の小球体23bの画像が写るようにすると、この小球体23bは、中心軸ぶれの補正ファントムに用いられる中球体22bより小さくなるのである。
【0059】
図7は、画像補正ファントムを用いたX線CT撮影装置の一例の基本構成図、図8は、図7の画像補正手段の詳細な構成図である。
【0060】
このX線CT撮影装置20は、X線撮影手段A、X線ビーム調整手段B、旋回アーム駆動制御手段C、演算処理手段6、表示モニターE、被写体Oを固定保持する被写体保持手段4、被写体保持手段4を移動させる被写体位置移動手段5、主フレーム10、操作部11、操作パネル10eなどを備えている。
【0061】
X線撮影手段Aは旋回アーム3を有しており、この旋回アーム3は、X線発生器1と2次元X線撮像手段2とを対向した状態で吊り下げ配置している。
【0062】
X線発生器1に備えられたX線ビーム調整手段Bは、X線ビーム幅制限手段B1、X線ビームコントローラB2、出射制御スリットB3を有しており、X線管より発射するX線ビームをX線ビーム幅制限手段B1で調整して、所望のビーム幅のX線コーンビーム1aが放射できるようになっている。
【0063】
一方の2次元X線撮像手段2は、X線IIの表面に設けたシンチレータ層に当たったX線が可視光に変換され、この可視光を光電変換器により電子に変換し電子増倍してこの電子を蛍光体により可視光に変換しレンズを通して2次元配列されたCCD(固体撮像素子)カメラで撮影する構成である。
【0064】
撮像手段としては、これ以外にカドミウムテルル(CdTe)やMOSセンサなどのX線2次元撮像手段、シンチレータとグラスファイバとCCDとの組み合わせのCCD撮像手段などの公知の二次元X線撮像手段が使用できる。
【0065】
旋回アーム3には、XYテーブル31と昇降制御モータ32と回転制御モータ33とが設けられており、X軸制御モータ31a、Y軸制御モータ31bを制御することによって、その回転中心3aをXY方向に位置調整可能とし、昇降制御モータ32を駆動することによって上下に昇降するとともに、撮影時には回転制御モータ33を等速度で駆動させて旋回アーム3を被写体Oの周りに旋回できるようにしている。この昇降制御モータ32は、旋回アーム3のアーム上下位置調整手段を構成している。
【0066】
また、旋回アーム3の回転中心3a、つまり、旋回軸が鉛直に設けられ、旋回アーム3が水平に回転し、X線コーンビーム1aが水平に局所照射されるので、装置を占有床面積の少ない縦型として構成することができる。
【0067】
この回転制御モータ33は、旋回アーム3の旋回駆動手段を構成しており、サーボモータなどのように、その回転速度、回転位置を自由に制御することができるモータを用い、また、旋回アーム3の回転中心3aに軸直結で設置されている。
【0068】
したがって、旋回アーム3を等速度又は可変速で回転をさせることができるとともに、その回転位置も時間軸に沿って知ることができるので、タイミングを合わせて、2次元X線撮像手段2でX線透過画像を取り出すのに都合がよく、また、芯振れがなくコーンビームX線CT撮影を有効に実施することができる。
【0069】
旋回アーム3の回転中心3aには、中空部3bが設けられている。このような中空部3bを設けるためには、回転中心3a上に有る関連部品に全て、中空孔を設ける必要があるが、例えば、回転制御モータ33としては、そのために、中空軸を使用したサーボモータを使用することができる。
【0070】
この中空部3bは、旋回アーム3に吊り下げ配置されたX線発生器1と2次元X線撮像手段2と、主フレーム10側に設けた操作部11との間の接続線を配置するためのものである。
【0071】
回転部分に対して、電気配線を接続する場合、その接続線の配置方法が問題になるが、このように、旋回アーム3の回転中心3aを通して接続線を配置すると、回転による捻じれなどの影響を最小限にすることができるとともに、配線の美観上も好ましい効果を得ることができる。
【0072】
旋回アーム駆動制御手段Cは、この実施例ではXYテーブルからなる位置調整手段31と、昇降制御後モータ32と、回転制御モータ33とを組み合わせて構成されるが、このような構成に限られない。最も簡易な構造では、旋回アーム3の中心3aは、手回しハンドルを操作して、任意の位置に設定できるようにしてもよい。
【0073】
被写体保持手段4は、被写体(患者)Oを座位で保持する椅子4bと、この椅子4bの背部に設けられた頭部固定手段4aとを備えている。
【0074】
被写体位置移動手段5は、被写体保持手段4を左右に移動させるX軸制御モータ51、前後移動させるY軸制御モータ52、上下移動手段のZ軸制御モータ53を備えている。
【0075】
これらのモータ51〜54で駆動されるX軸、Y軸、Z軸直線移動テーブル(不図示)は、それぞれ周知のクロスローラガイドや、通常のベアリングとガイドを組み合わせたものなどで構成され、正確に直線移動ができるものである。また、モータ51〜54による、これらのX軸、Y軸、Z軸直線移動テーブルの移動は、ラックとピニオン方式や、ボールネジ方式や、通常のネジ軸を用いる方式などを適用できるが、正確に位置決めできるものが望ましい。
【0076】
こうして、被写体Oを椅子4bに座らせ、頭部固定手段4aで、被写体Oの頭部を固定保持し、被写体位置移動手段5を用いて、旋回アーム3の回転中心3aに、被写体Oの内部の局所部位の中心に合わせることができる。一方、被写体保持手段4を被写体位置移動手段5で移動させる替わりに、XYテーブル31と昇降制御モータ32を用いて、旋回アーム3側を移動させて、旋回アーム3の回転中心3aを、被写体Oの内部の局所部位の中心に合わせることも可能である。
【0077】
この装置20では、この回転中心3aの位置合わせ、つまり、X線撮影位置設置のために、被写体側を移動させる被写体位置移動手段5と、照射側である旋回アーム3を移動させるXYテーブル31と昇降制御モータ32の双方を備えているが、どちらか一方だけ設けてもよい。なお、コーンビームX線CT撮影の場合には、回転中心3aのブレのないことが重要なので、旋回アーム3側は、旋回だけとし、旋回中心3aは固定とするのが望ましい。
【0078】
演算処理手段6は、画像処理解析に高速で作動する演算プロセッサを含んでおり、2次元X線撮像手段2上に生成されたX線透過画像を前処理した後、所定の演算処理を実行することによって、X線を透過させた物体内部の3次元X線吸収係数データを算出し、また、このデータの投影面への投影などの演算を行って、外付けの表示手段Eに投影画像や、X線パノラマ画像を表示させ、また必要な記憶媒体に画像情報として記憶させる。
【0079】
この撮影装置では、この演算処理手段6に、上述した画像補正テーブルを作成し、この画像補正テーブルを用いて、X線原画像を補正するを行う画像補正手段6aを備えている。
【0080】
画像補正手段6aは、この目的を達成するため、種々の補正処理や演算処理を行う補正処理手段61、生成された画像補正テーブルをX線CT撮影の各旋回角度毎に記憶保存する画像補正テーブルメモリ62、管面歪み補正テーブルメモリ63、中心軸ぶれ補正テーブルメモリ64、中心軸ずれ補正テーブルメモリ65を備えている。
【0081】
図9(a)はX線CT撮影装置の一例の外観正面図、(b)は側面図(c)は部分側面図である。
【0082】
X線CT撮影装置20は、門型の非常に剛性の高い構造体である主フレーム10を全体の支持体として構成されている。
【0083】
この主フレーム10は、X線発生器1と2次元X線撮像手段2とを対向した状態で吊り下げ配置した旋回アーム3を回転可能に支持するアーム10a、このアーム10aの基端部を固定保持している1対の横ビーム10b、この横ビーム10bを支えている一対の縦ビーム10c、一対の縦ビーム10cが固定載置され、この装置20全体の基礎となっているベース10dから構成されている。
【0084】
この主フレーム10を構成する部材は、それぞれ、剛性の高い鋼鉄材が用いられ、また、適宜、筋交いや、角補強部材が設けられて変形に強いものとなっており、回転時に、旋回アーム3の回転中心3aが変動しないようになっている。
【0085】
このように主フレーム10は、旋回アーム3の旋回振れが生じないような構造体としているので、特に、旋回振れがないことが要求されるX線CT撮影装置として、ふさわしい。
【0086】
操作パネル10eは、主フレーム10の一方の縦ビーム10cの表面で、術者が、立位で操作がし易いような位置に設けられている。被写体保持手段4の椅子4bを載置しているのは、図7で説明した被写体位置移動手段5であり、椅子4bをX,Y,Z方向に、つまり、前後左右上下方向に移動させ、また、椅子4bの背板4baを傾動させて、被写体Oの頭部を傾動保持させることができる。
【0087】
本発明では、図9(b),(c)に示すように、この被写体保持手段4の頭部固定手段4aに着脱交換可能にファントム支持台26、27が取り付けられるようになっており、この支持台26、27に上述した画像補正ファントム22、23を載置することができる。
【0088】
ファントム支持台26は、図9(b)に示すように、画像補正ファントム22を旋回アーム3の回転中心軸3a付近に位置させるものである。
【0089】
この図9と図3から解るように、画像補正ファントム22をファントム支持台26に載せるだけで、このファントム22がいわゆる被写体保持位置付近に設置され、かつ、このファントム22に内蔵された中球体22b(A)、22b(B)がほぼ回転中心軸3a方向に所定間隔で配置される。
【0090】
ファントム支持台27は、図9(c)に示すように、支持台26に比べ、画像補正ファントム23を、より水平方向に旋回アーム3の回転中心軸3aから離れた位置に位置させるものである。
【0091】
この図9と図5から解るように、画像補正ファントム23をファントム支持台27に載せるだけで、このファントム23がいわゆる被写体保持位置から所定距離離れた位置に設置され、かつ、このファントム23に内蔵された小球体23b(A)、23b(B)がほぼ回転中心軸3a方向に所定間隔で配置される。
【0092】
【発明の効果】
請求項に記載のX線画像補正方法によれば、X線発生器からコーンビームX線を被写体に照射するコーンビームX線CT撮影における回転中心軸の軸方向の位置ずれを補正するX線画像補正方法であって、X線透過性の支持体内に、この支持体よりもX線吸収係数の高い2つの球体を収容して成り、この支持体を回転中心軸より所定位置離れた位置に設置し、X線発生器と2次元撮像手段とを回転中心軸を中心にして180度反転させ、それぞれの位置で2次元撮像手段上で得られた2つの球体のX線原画像に対して、所定の演算を行って2次元撮像手段上におけるX線原画像の基準位置を算出するようにしているので、補正テーブルを作成することが容易になる。
【0093】
請求項に記載のX線画像補正方法によれば、請求項1に記載の補正方法を有効利用して、その効果を発揮させながら、かつ、X線CT撮影時に発生が考えられ得る歪みを補正する中心軸ずれ補正テーブルを用いて作成した画像補正テーブルによって、X線原画像を補正することができ、X線CT撮影時に発生が考えうる全ての歪みを除去することができる。
【0094】
請求項に記載のX線画像補正方法は、請求項2のX線画像補正方法において、更に、X線の支持体内にこの支持体よりもX線吸収係数の高い2つの球体を収容して成り、この支持体を、球体がX線CT撮影におけるほぼ回転中心軸方向に上下になるよう被写体保持位置に設置し、2つの球体の画像上の座標方向の動きの平均値から補正に用いる角度とシフト量とを算出することで、画像補正テーブルを作成するので、補正テーブルを作成することが容易になる。
【図面の簡単な説明】
【図1】 (a)は本発明の画像補正ファントムの一例を示す正面図、(b)はこれを二次元X線撮像手段に装着した状態を示す縦断面図、(c)は(b)で撮像されたX線原画像の一例を示す図
【図2】 (a)から(e)は、図1の画像補正ファントムを用いた補正方法の概念図
【図3】 画像補正ファントムの参考例を示す外観図
【図4】 図3の画像補正ファントムを用いた補正方法の概念図
【図5】 本発明の画像補正ファントムの他例を示す外観図
【図6】 図5の画像補正ファントムを用いた補正方法の概念図
【図7】 本発明の画像補正ファントムを用いたX線CT撮影装置の一例の基本構成図
【図8】 図7の画像補正手段の詳細な構成図
【図9】 (a)は図7のX線CT撮影装置の外観正面図、(b)は側面図、(c)は部分側面図
【図10】 (a)は従来の画像補正ファントムの一例を示す図、(b)は、これを用いたX線CT撮影装置の基本構成図
【図11】 (a)は従来の画像補正ファントムの他例を示す図、(b)は、これを用いたX線CT撮影装置の基本構成図
【図12】 (a)は従来の画像補正ファントムの他例を示す図、(b)は、これを用いたX線CT撮影装置の基本構成図
【符号の説明】
1 X線発生器
1a X線コーンビーム
2 二次元X線撮像手段
21 X線画像補正ファントム
21a プリント基板
21b 小円体
22 X線画像補正ファントム
22a 支持体
22b 中球体
23 X線画像補正ファントム
23a 支持体
23b 小球体
3 旋回アーム
3a 回転中心軸
4 被写体保持手段
4a 椅子
4b 頭部固定手段
6 画像処理手段
6a 画像補正手段
20 X線CT撮影装置
[0001]
BACKGROUND OF THE INVENTION
  According to the present invention, cone beam X-ray CT imaging, that is, a subject is fixedly held, and an X-ray generator and a two-dimensional X-ray imaging unit are disposed opposite to each other with the subject as a center, and the X-ray generator is used to detect the subject. By rotating the X-ray generator and the two-dimensional X-ray imaging means within an angle range according to imaging conditions while locally irradiating an X-ray cone beam including only a local site to be imaged, In X-ray imaging in which an X-ray image of a local part is captured, and this X-ray image is processed to calculate a three-dimensional X-ray absorption coefficient of the local part, an X-ray image of the subject and the X-ray original image of the subject is obtained. Correcting positional deviationX-ray image correction methodAbout.
[0002]
[Prior art]
  In X-ray CT (Computed Tomography) imaging, as a two-dimensional X-ray imaging means, an image intensifier is abbreviated to an X-ray II (X-ray fluorescence multiplier) and a CCD camera (hereinafter referred to as “XII imaging tube”). In the case of this XII imaging tube, the X-ray light receiving surface is a curved surface, whereas the light receiving surface of a CCD camera that converts the light visualized here into an electric signal is a flat surface. , Tube distortion caused by detecting a curved surface in a plane, and distortion such as magnetic distortion due to geomagnetism during imaging, that is, positional deviation between the subject and the X-ray original image of the subject is corrected. Various proposals have been made.
[0003]
  FIG. 10A is a view showing an example of a conventional image correction phantom used for correcting the tube surface distortion, and FIG. 10B is a basic configuration diagram of an X-ray CT imaging apparatus using the image correction phantom. No. 24248.
[0004]
  The image correction phantom 121 shown in FIG. 10A is obtained by forming X-ray transparent holes 121b in a square lattice shape on a metal flat plate 161a that is X-ray opaque, and is preferably used for correcting the tube surface distortion. However, it took a lot of work to process the hole 121b in the metal flat plate 121a, and the cost was high.
[0005]
  In the X-ray CT imaging apparatus shown in FIG. 10B, an X-ray generator 101 and a two-dimensional X-ray imaging unit 102 arranged opposite to each other with a subject O interposed therebetween are installed on a rotating plate 103 and rotated. Then, the X-ray generator 101 turns and irradiates the subject O with the X-ray cone beam, obtains an X-ray image of the subject O by the two-dimensional X-ray imaging means 102, and the image acquisition processing device 106 obtains the obtained X-ray. The image is corrected using a position correspondence table, etc., and the corrected X-ray image is reconstructed to obtain a three-dimensional X-ray absorption coefficient (reconstruction data) of the subject O, which is displayed as an image. In this case, the image correction phantom 121 was mounted on the tube surface of the two-dimensional X-ray imaging unit 102 and used for image correction.
[0006]
  Further, as a cause of the positional deviation between the subject and the X-ray original image of the subject, a rotation center axis for rotating the X-ray generator and the two-dimensional X-ray imaging means arranged opposite to each other about the subject is used. There have been various proposals for correcting this, as well as the rotational center axis, the axial displacement of the X-ray generator and the two-dimensional X-ray imaging means.
[0007]
  FIG. 11A is a view showing another example of such an image correction phantom, and FIG. 11B is a basic configuration diagram of an X-ray CT imaging apparatus using the image correction phantom, which is described in JP-A-9-173330. It is a thing.
[0008]
  An image correction phantom 221 shown in FIG. 11A is obtained by vertically arranging a plurality of X-ray opaque spheres 221b inside an X-ray transparent support 221a, and corrects the shake of the rotation center axis. However, data processing up to the creation of the correction table is complicated by the amount of the sphere 221b.
[0009]
  The X-ray CT imaging apparatus shown in FIG. 11 (b) is indispensable for X-ray CT imaging as in FIG. 10 (b), and the X-ray generator 201 and two-dimensional X-ray imaging arranged opposite to each other with the subject O interposed therebetween. In addition to means 202, rotating plate 203, image processing means 206, etc., in order to correct an X-ray original image using image correction phantom 221, image distortion correction means 206 a, image distortion correction table creation means 261, etc. are provided. It was.
[0010]
  FIG. 12A is a diagram showing another example of the same image correction phantom, and FIG. 12B is a basic configuration diagram of an X-ray CT imaging apparatus using the same, which is described in Japanese Patent Laid-Open No. 2000-201918. Is.
[0011]
  An image correction phantom 321 shown in FIG. 12A has an X-ray opaque rod 321b installed in the vertical direction inside an X-ray transparent support 321a. Can be suitably used to correct the vertical misalignment of the detector and the two-dimensional X-ray imaging means. However, in order to create a correction table, the data irradiated to the entire phantom 321 is necessary, which is complicated. It was processing.
[0012]
  The X-ray CT imaging apparatus shown in FIG. 12 (b) is indispensable for X-ray CT imaging as in FIG. 11 (b), and an X-ray generator 301 and two-dimensional X-ray imaging arranged opposite to each other with the subject O interposed therebetween. Means 302, a rotating plate 303, a data processing unit 306 corresponding to an image processing unit, and the like, and a geometry estimation unit corresponding to an image distortion correction unit in order to correct an X-ray original image using the image correction phantom 361. 306a and the like.
[0013]
[Problems to be solved by the invention]
  However, although these conventional image correction phantoms achieve the intended purpose, they still have the problems to be solved as described above.
[0014]
  The present invention is intended to solve such problems and solves various problems derived from the purpose of distortion correction.Image correction methodThe purpose is to provide.
[0015]
[Means for Solving the Problems]
  The X-ray image correction method according to claim 1 comprises:The cone beam X-ray CT that irradiates the subject with cone beam X-rays from the X-ray generator while rotating the X-ray generator and the two-dimensional imaging means arranged opposite to each other with the subject interposed therebetween. The axial displacement of the rotation center axis during shootingX-ray image correction method for correctionThe X-ray transparent support body accommodates two spheres having a higher X-ray absorption coefficient than the support body, and the support body is installed at a position away from the rotation center axis by a predetermined position. The X-ray generator and the two-dimensional imaging means are inverted 180 degrees around the rotation center axis, and the original X-ray images of the two spheres obtained on the two-dimensional imaging means at the respective positions. On the other hand, a predetermined calculation is performed to calculate the reference position of the X-ray original image on the two-dimensional imaging means.
[0016]
  Claim 2The X-ray image correction method described inIn claim 1,A center axis deviation correction table is created to correct the axial position deviation of the rotation center axis of X-ray CT imaging.,thisAn image correction table is created for each rotation angle of X-ray CT imaging using the center axis deviation correction table, and each X-ray original image is corrected using the image correction table.
[0017]
  The correction method is the correction according to claim 1.MethodThe X-ray original image is corrected by an image correction table created using a center axis deviation correction table that corrects distortion that may occur during X-ray CT imaging while effectively utilizing is doing.
  Claim 3The X-ray image correction method described in claim 12In the X-ray image correction method, two spheres having an X-ray absorption coefficient higher than that of the support are accommodated in the X-ray support, and the support is substantially the same as the sphere in X-ray CT imaging. The angle and shift amount used for correction from the average value of the movement in the coordinate direction on the image of the two spheres installed at the subject holding position so as to be up and down in the rotation center axis directionAnd calculatingThe image correction table is created.
[0018]
DETAILED DESCRIPTION OF THE INVENTION
  Of the present inventionEmbodiments will be described with reference to the accompanying drawings.
[0019]
  FIG. 1 (a) shows the present invention.Used inA front view showing an example of an image correction phantom, (b) is a longitudinal sectional view showing a state in which the image correction phantom is mounted on a two-dimensional X-ray imaging means, and (c) is an example of an original X-ray image taken in (b). FIG.
[0020]
  An image correction phantom 21 shown in FIG. 1A is composed of an X-ray transmissive printed board 21a and a small circle 21b having a high X-ray absorption coefficient, and as shown in the figure, a grid-like shape is formed on the printed board 21a. Assuming a lattice, small circles 21b are arranged at positions that are intersections of the lattice.
[0021]
  The printed circuit board 21a may be a plate-like body formed of a material having a low X-ray absorption coefficient such as a plastic material such as an acrylic resin, vinyl chloride, or polycarbonate in addition to a synthetic resin, for example, a glass epoxy resin or a phenol resin. Specifically, a printed circuit board that is widely used for arranging and arranging electronic devices can be obtained at low cost and is suitable.
[0022]
  As the X-ray opaque material of the small circular body 21b, metals such as copper, tungsten, iron-nickel-chromium alloy, and stainless steel are preferable, and the small circular body 21b has a diameter of about 1.5 mm. However, the present invention is not limited to this, and is relatively determined according to the overall size of the phantom, that is, the size of the tube surface of the two-dimensional X-ray imaging means on which the phantom is mounted. .
[0023]
  The correction phantom 21 uses a general-purpose X-ray transmissive printed board 21a as a plate-like object for arranging small circles 21b having a high X-ray absorption coefficient serving as a reference point for correction in a grid pattern. It will be very cheap. In addition, the lattice arrangement can be achieved by simply forming the small circles 21b on the printed circuit board 21a by etching, thereby reducing the cost.
[0024]
  The two-dimensional X-ray image pickup means 2 is an XII image pickup tube as shown in FIG. 1B, and is a well-known one, but its configuration will be briefly described. An input window 2a that is a tube surface, an input surface substrate. 2b, input phosphor screen 2c, photocathode 2d, focus electrode 2e, anode 2f, output phosphor screen 2g, output window 2h, optical lens 2i, CCD camera 2j, focus power supply E1, anode power supply E2, etc. Upon receiving the incident line, this is converted into visible light, and the visible light is imaged by a CCD camera and output as an electrical signal.
[0025]
  The two-dimensional X-ray imaging means 2 is provided with an attachment frame 25 for attaching the image correction phantom 21 to the input window 2a.
[0026]
  As shown in FIG. 1B, the image correction phantom 21 is used by being attached to the tube surface of the two-dimensional X-ray imaging means 2 by this attachment frame 25. In this state, X The line original image IO is shown in FIG.
[0027]
  The correction table indicating the relative positional relationship between the small-circle image IO (21b) constituting the X-ray original image IO and the small-circle body 21b of the image correction phantom 21 is used to correct the tube surface distortion of the two-dimensional X-ray imaging means. Then, distortion such as magnetostriction due to geomagnetism during imaging, that is, positional displacement between the subject and the X-ray original image of the subject is corrected. This correction table is called a tube surface distortion correction table.
[0028]
  FIGS. 2A to 2C are conceptual diagrams of a correction method using the image correction phantom of FIG. 1, and the creation of the tube surface distortion correction table will be described using these.
[0029]
  1) Binarization of X-ray original image IO: This is for obtaining the position of the small circle image IO (21b) in the original image IO. FIG. 2A shows one small circular body 21b of FIG. 1A and a grid surrounding it. Here, since the size of the lattice is 3 mm and the diameter of the small circle 21b is 1.5 mm, the area ratio of the small circle image IO (21b) is
    π * (1.5 / 2) * (1.5 / 2) /3*3≈0.196
It becomes. This is a threshold for the X-ray original image IO. Here, “*” means multiplication, and “/” means division.
[0030]
  Next, as shown in FIG. 2B, an X-ray original image IO of 320 pixels * 240 pixels is extrapolated by applying a mirror image process to the image data within the boundary around the field boundary in the region of 512 pixels * 512 pixels. The result is subjected to FFT (Fourier transform), then a band pass filter as shown in FIG. 2C is applied in the frequency domain, and further IFFT (inverse Fourier transform) is performed. .196)).
[0031]
  2) Creation of correction table: This indicates the positional relationship between the original image IO and the small circle 21b of the phantom 21. As shown in FIG. 2 (d), this problem generally occurs from four points A, B, C, D on the xy coordinate system to four points A ′, B ′, C ′, on the uv coordinate system. It can be understood that this is a conversion problem to D ′. Here, in addition to simple rotation and enlargement of the coordinate system, the enlargement ratio changes, so the following extended affine transformation is used.
   u = a11 * xy + a10 * x
+ a01 * y + a00
   v = b11 * xy + b10 * x
+ b01 * y + b00
  In this case, there are four coefficients for each equation of u and v, which is suitable for obtaining the relationship between the four points as in this problem.
[0032]
  3) Define the whole affine transformation by a set of small area affine transformations.
[0033]
  As shown in FIG. 2 (e), the entire affine transformation is defined as a set of small area affine transformations corresponding to the hatched portion on the left side.
[0034]
  By such a procedure, a tube surface distortion correction table is obtained. In this state, the tilt of the imaging system and the vertical and horizontal tilts, that is, the rotation center axis shake of X-ray CT imaging and the rotation center axis of X-ray CT imaging. In order to achieve this, it is necessary to create a center axis shake correction table and a center axis deviation correction table using two types of image correction phantoms described below. There is.
[0035]
3A and 3B are external views showing a reference example of the image correction phantom, and FIG. 4 is a conceptual diagram of a correction method using the image correction phantom of FIG.
[0036]
  The image correction phantom 22 is composed of an X-ray transparent support 22a and two medium spheres 22b (A) and 22b (B) having a high X-ray absorption coefficient. The spheres 22b (A) and 22b (B) are arranged so that when the support 22a is placed at the subject holding position, these spheres 22b are vertically moved substantially in the direction of the rotation center axis 3a of X-ray CT imaging. (See FIG. 7).
[0037]
  Here, the material of the support 22a is preferably an acrylic resin having a good X-ray transparency, transparent, easily available and easy to process. The material of the intermediate sphere 22b is preferably a steel ball having a high X-ray absorption coefficient, and its diameter is preferably 3 mm to 15 mm, but is not limited thereto.
[0038]
  The correction phantom 22 has a simple configuration and can be manufactured at a low cost, and it becomes easy to create a correction table using the phantom.
[0039]
  The phantom 22 is used to create a center axis shake correction table for correcting the rotation center axis shake in X-ray CT imaging. For this purpose, the phantom 22 is placed on a support base 26 as shown in FIG. It is set so that its position is almost the vicinity of the rotation center axis 3a in FIG. FIG. 4 shows a method of performing X-ray imaging in this state and calculating the rotation center axis shake from the obtained X-ray original image.
[0040]
  Further, the correction phantom 22 may be a rectangular parallelepiped as shown in FIG. 3A or a cylindrical body as shown in FIG.
[0041]
  In FIG. 4, the images of the two spheres 22 b (A) and 22 b (B) are simply indicated by the symbols A and B.
[0042]
  First, A and B in this figure need to be corrected by the tube surface distortion correction table (the same applies to FIG. 6). Here, the purpose is to obtain the angle θ and the shift amount XS.
[0043]
  Since A and B move to the left and right as shown in the figure by the rotation of the swivel arm 360, the points A 'and B' can be easily obtained by obtaining the average value of the x coordinates of A and B. Thus, the angle θ and the shift amount XS are calculated. More precisely, the movement of A and B in the x-coordinate direction can be subjected to Fourier transform, and the respective DC components can be used as the x-coordinates of the points A ′ and B ′.
[0044]
  Since it is assumed that the angle θ is small, the y-coordinates of A and B are constant (average value), and it can be assumed that A and B only move in the horizontal direction. Further, only one set of the angle θ and the shift amount XS is calculated for the 360 degree rotation of the swing arm 3.
[0045]
  FIG.Other examples of image correction phantomFIG. 6 is a conceptual diagram of a correction method using the image correction phantom of FIG.
[0046]
  The image correction phantom 23 is composed of an X-ray transparent support 23a and two small spheres 23b (A) and 23b (B) having a high X-ray absorption coefficient. The small spheres 23b (A) and 23b (B) are arranged so that the small spheres 23b are vertically arranged in the direction of the rotation center axis of X-ray CT imaging when the support 23a is placed at the subject holding position. (See FIG. 7).
[0047]
  Here, the material of the support 23a may be the same as that of the support 22a. The material of the small sphere 23b may be the same as that of the intermediate sphere 22b, and the diameter is preferably 1 to 3 mm, but is not limited thereto.
[0048]
  The correction phantom 23 has a simple configuration and can be manufactured at a low cost, and it becomes easy to create a correction table using the phantom.
[0049]
  The phantom 23 is used to create a center axis deviation correction table for correcting the axial position deviation of the rotation center axis of X-ray CT imaging. For this purpose, the phantom 23 is shown in FIG. Then, it is placed on the support base 27 so that the position thereof is substantially away from the rotation center axis 3a in FIG. 7, that is, a position away from the subject holding position by a predetermined distance. FIG. 6 shows a method of performing X-ray imaging in this state and calculating the axial displacement of the rotation center axis from the obtained X-ray original image.
[0050]
  In FIG. 6, the two small spheres 23b (A) and 23b (B) are simply indicated by the symbols A and B. Reference numeral 1 denotes an X-ray generator, reference numeral 2 denotes a two-dimensional X-ray imaging means, and reference numeral 3a denotes a rotation center axis 3a of the swivel arm 3 (see FIG. 7).
[0051]
  Reference numerals a1 and b1 denote the X-ray sources of the small spheres A and B when the X-ray generator 1 is irradiated with the X-ray generator 1 in the position farthest away from the small spheres A and B (indicated by a one-dot chain line). The vertical positions on the image, the symbols a2 and b2, are the positions where the X-ray generator 1 is closest to the small spheres A and B (that is, positions rotated 180 degrees around the rotation center axis 3a from the previous position). The vertical positions of the small spheres A and B on the X-ray original image when these are irradiated in the state (shown by a solid line). Reference symbol c is the vertical position occupied by the reference position of the X-ray original image with respect to the rotation center axis 3a.
[0052]
  Between the four positions a1, b1, a2, b2 thus obtained and the vertical position c of the rotation center shaft 3a,
   (a1-c) / (a2-c) =
(b1-c) / (b2-c)
The relationship of
    c = (a2 * b1-b2 * a1) / ((a2-b2)-(a1-b1))
It becomes.
[0053]
  Thus, the vertical position c is obtained, and this vertical position c becomes the vertical position of the X-ray original image with respect to the rotation center axis 3a on the machine, so that the vertical movement type of the coordinate axis can be generated, and this becomes the central axis deviation correction table. . In addition, what is necessary is just to prepare one set of this correction table about the perimeter of irradiation.
[0054]
  Thus, when a tube surface distortion correction table is obtained for each turning angle of X-ray CT imaging, and a set of center axis shake correction table and center axis deviation correction table is obtained for the entire circumference, the tube surface distortion correction table is obtained. Then, coordinate axis rotation conversion by the center axis shake correction table and coordinate axis movement conversion by the center axis deviation correction table are performed, and an image correction table is created for each turning angle of X-ray CT imaging.
[0055]
  Such an image correction table uses a tube surface distortion correction table, a center axis shake correction table, and a center axis deviation correction table for correcting all possible distortions that may occur during X-ray CT imaging. Can be removed.
[0056]
  It should be noted that those having different sizes as the correction reference points, such as the middle sphere 22b of the correction phantom 22 for the center axis shake correction table and the small sphere 23b of the correction phantom 23 for the center axis deviation correction table, It corresponds to each correction purpose.
[0057]
  In other words, in the case of the center axis blur correction, the correction phantom 22 is placed almost in the vicinity of the rotation center axis, that is, in the subject holding position, but this position is in the middle of the X-ray generator 1 and the two-dimensional X-ray imaging means 2. Yes, during turning, the size of the X-ray image of the middle sphere 22b of the correction phantom does not change much. Therefore, if the size of each image of the two spheres 22b is reflected on the imaging surface of the two-dimensional X-ray imaging means 2, the size is as described above.
[0058]
On the other hand, in the case of center axis deviation correction, it is necessary to install the correction phantom 23 at a position away from the rotation center axis and at a predetermined distance from the subject holding position. The size of the X-ray image of the small sphere 23b changes from large to small. Therefore, even when the size of the X-ray image of the small sphere 23b is maximized at the rotational position where the correction phantom 23 is closest to the X-ray generator 1, the respective surfaces of the two-dimensional X-ray imaging means 2 are set on the imaging surface. If the images of the two small spheres 23b are captured, the small sphere 23b is smaller than the middle sphere 22b used in the correction phantom of the center axis shake.
[0059]
  FIG.Image correction phantomFIG. 8 is a detailed configuration diagram of the image correction unit in FIG. 7.
[0060]
  The X-ray CT imaging apparatus 20 includes an X-ray imaging unit A, an X-ray beam adjusting unit B, a turning arm drive control unit C, an arithmetic processing unit 6, a display monitor E, a subject holding unit 4 that holds the subject O fixedly, a subject A subject position moving unit 5 for moving the holding unit 4, a main frame 10, an operation unit 11, an operation panel 10 e and the like are provided.
[0061]
  The X-ray imaging means A has a swivel arm 3, and the swivel arm 3 is suspended from the X-ray generator 1 and the two-dimensional X-ray imaging means 2.
[0062]
  The X-ray beam adjusting means B provided in the X-ray generator 1 has an X-ray beam width limiting means B1, an X-ray beam controller B2, and an emission control slit B3, and the X-ray beam emitted from the X-ray tube Is adjusted by the X-ray beam width limiting means B1, so that an X-ray cone beam 1a having a desired beam width can be emitted.
[0063]
  One two-dimensional X-ray imaging means 2 converts X-rays hitting a scintillator layer provided on the surface of X-ray II into visible light, converts the visible light into electrons by a photoelectric converter, and multiplies the electrons. This electron is converted into visible light by a phosphor and photographed with a CCD (solid-state imaging device) camera arranged two-dimensionally through a lens.
[0064]
  In addition to this, known two-dimensional X-ray imaging means such as an X-ray two-dimensional imaging means such as cadmium tellurium (CdTe) or a MOS sensor, or a CCD imaging means using a combination of a scintillator, a glass fiber, and a CCD are used. it can.
[0065]
  The revolving arm 3 is provided with an XY table 31, an elevation control motor 32, and a rotation control motor 33. By controlling the X-axis control motor 31a and the Y-axis control motor 31b, the rotation center 3a is set in the XY direction. The position of the revolving arm 3 can be moved up and down by driving the elevating control motor 32 and the revolving arm 3 can be rotated around the subject O by driving the rotation control motor 33 at a constant speed during photographing. The lift control motor 32 constitutes an arm vertical position adjusting means for the swing arm 3.
[0066]
  Further, the rotation center 3a of the swivel arm 3, that is, the swivel axis is provided vertically, the swivel arm 3 rotates horizontally, and the X-ray cone beam 1a is locally irradiated horizontally, so that the apparatus occupies a small floor space. It can be configured as a vertical type.
[0067]
  This rotation control motor 33 constitutes a turning drive means for the turning arm 3, and uses a motor that can freely control its rotation speed and rotation position, such as a servo motor. The shaft is directly connected to the rotation center 3a.
[0068]
  Therefore, the swivel arm 3 can be rotated at a constant speed or a variable speed, and the rotation position can be known along the time axis. It is convenient for taking out a transmission image, and there is no center shake, and cone beam X-ray CT imaging can be effectively performed.
[0069]
  A hollow portion 3 b is provided at the rotation center 3 a of the turning arm 3. In order to provide such a hollow portion 3b, it is necessary to provide a hollow hole in all the related parts on the rotation center 3a. For example, the rotation control motor 33 is a servo using a hollow shaft. A motor can be used.
[0070]
  The hollow portion 3b is used for arranging a connection line between the X-ray generator 1 and the two-dimensional X-ray imaging means 2 suspended from the swivel arm 3 and the operation unit 11 provided on the main frame 10 side. belongs to.
[0071]
  When electrical wiring is connected to the rotating portion, the arrangement method of the connecting wire becomes a problem. However, if the connecting wire is arranged through the rotation center 3a of the swivel arm 3 as described above, the influence of twisting due to rotation, etc. Can be minimized, and a favorable effect can be obtained from the aesthetics of the wiring.
[0072]
  In this embodiment, the swing arm drive control means C is configured by combining the position adjusting means 31 composed of an XY table, the motor 32 after lift control, and the rotation control motor 33, but is not limited to such a configuration. . In the simplest structure, the center 3a of the swivel arm 3 may be set at an arbitrary position by operating a handwheel handle.
[0073]
  The subject holding means 4 includes a chair 4b that holds the subject (patient) O in a sitting position, and a head fixing means 4a provided on the back of the chair 4b.
[0074]
  The subject position moving unit 5 includes an X-axis control motor 51 that moves the subject holding unit 4 left and right, a Y-axis control motor 52 that moves back and forth, and a Z-axis control motor 53 that moves up and down.
[0075]
  The X-axis, Y-axis, and Z-axis linear movement tables (not shown) driven by these motors 51 to 54 are each composed of a well-known cross roller guide or a combination of ordinary bearings and guides. Can be moved linearly. Moreover, the movement of these X-axis, Y-axis, and Z-axis linear movement tables by the motors 51 to 54 can be applied with a rack and pinion method, a ball screw method, a method using a normal screw shaft, or the like. What can be positioned is desirable.
[0076]
  In this way, the subject O is seated on the chair 4 b, the head of the subject O is fixedly held by the head fixing means 4 a, and the subject position moving means 5 is used to bring the inside of the subject O into the rotation center 3 a of the swivel arm 3. It can be adjusted to the center of the local site. On the other hand, instead of moving the subject holding means 4 by the subject position moving means 5, the rotation arm 3 side is moved using the XY table 31 and the lifting control motor 32, and the rotation center 3 a of the turning arm 3 is moved to the subject O It is also possible to match with the center of the local site inside.
[0077]
  In this apparatus 20, in order to align the rotation center 3a, that is, to set the X-ray imaging position, the subject position moving means 5 for moving the subject side, and the XY table 31 for moving the swivel arm 3 on the irradiation side, Although both the lift control motors 32 are provided, only one of them may be provided. In the case of cone beam X-ray CT imaging, it is important that there is no blurring of the rotation center 3a. Therefore, it is desirable that only the turning arm 3 is turned and the turning center 3a is fixed.
[0078]
  The arithmetic processing means 6 includes an arithmetic processor that operates at high speed for image processing analysis, and after pre-processing the X-ray transmission image generated on the two-dimensional X-ray imaging means 2, executes predetermined arithmetic processing. Thus, the three-dimensional X-ray absorption coefficient data inside the object through which X-rays are transmitted are calculated, and calculation such as projection of this data on the projection surface is performed, and the projection image or An X-ray panoramic image is displayed and stored as image information in a necessary storage medium.
[0079]
  In this imaging apparatus, the arithmetic processing unit 6 includes an image correction unit 6a that creates the above-described image correction table and corrects an X-ray original image using the image correction table.
[0080]
  In order to achieve this object, the image correction unit 6a performs correction processing unit 61 that performs various correction processes and arithmetic processes, and an image correction table that stores and saves the generated image correction table for each rotation angle of X-ray CT imaging. A memory 62, a tube surface distortion correction table memory 63, a center axis shake correction table memory 64, and a center axis deviation correction table memory 65 are provided.
[0081]
  FIG. 9 (a)X-ray CT imaging systemFIG. 2B is a side view of the outer appearance of the example, FIG.
[0082]
  The X-ray CT imaging apparatus 20 is configured with a main frame 10 that is a portal-type extremely rigid structure as an overall support.
[0083]
  The main frame 10 has an arm 10a for rotatably supporting a swing arm 3 suspended from the X-ray generator 1 and the two-dimensional X-ray imaging means 2, and a base end portion of the arm 10a is fixed. A pair of horizontal beams 10b that are held, a pair of vertical beams 10c that support the horizontal beams 10b, and a pair of vertical beams 10c are fixedly mounted, and are composed of a base 10d that forms the basis of the entire apparatus 20. Has been.
[0084]
  The members constituting the main frame 10 are each made of a steel material having high rigidity, and are appropriately provided with braces or corner reinforcing members to be resistant to deformation. The rotation center 3a is not changed.
[0085]
  As described above, the main frame 10 has a structure that prevents the swing arm 3 from swinging. Therefore, the main frame 10 is particularly suitable as an X-ray CT imaging apparatus that requires no swing swing.
[0086]
  The operation panel 10e is provided on the surface of one vertical beam 10c of the main frame 10 at a position where the operator can easily operate in a standing position. The chair 4b of the subject holding means 4 is placed on the subject position moving means 5 described in FIG. 7, and the chair 4b is moved in the X, Y, and Z directions, that is, in the front-rear, left-right, up-down directions, Further, the back plate 4ba of the chair 4b can be tilted to hold the head of the subject O tilted.
[0087]
  In the present invention, as shown in FIGS. 9B and 9C, phantom support bases 26 and 27 are attached to the head fixing means 4a of the subject holding means 4 so as to be detachable and replaceable. The above-described image correction phantoms 22 and 23 can be placed on the support bases 26 and 27.
[0088]
  The phantom support 26 is for positioning the image correction phantom 22 in the vicinity of the rotation center axis 3a of the swivel arm 3, as shown in FIG.
[0089]
  As can be seen from FIGS. 9 and 3, the phantom 22 is installed in the vicinity of a so-called subject holding position by simply placing the image correction phantom 22 on the phantom support base 26, and the sphere 22 b incorporated in the phantom 22. (A) and 22b (B) are arranged at predetermined intervals substantially in the direction of the rotation center axis 3a.
[0090]
  As shown in FIG. 9C, the phantom support base 27 is for positioning the image correction phantom 23 at a position further away from the rotation center axis 3 a of the swivel arm 3 in the horizontal direction than the support base 26. .
[0091]
  As can be seen from FIGS. 9 and 5, the image correction phantom 23 is installed at a position away from a so-called subject holding position by a predetermined distance only by placing the image correction phantom 23 on the phantom support 27, and is built in the phantom 23. The small spheres 23b (A) and 23b (B) thus arranged are arranged at predetermined intervals substantially in the direction of the rotation center axis 3a.
[0092]
【The invention's effect】
  Claim1X-ray image correction described inMethodAccording toAn X-ray image correction method for correcting an axial displacement of a rotation center axis in cone beam X-ray CT imaging in which a subject is irradiated with cone beam X-rays from an X-ray generator. , Which contains two spheres having a higher X-ray absorption coefficient than this support, this support is placed at a predetermined position away from the rotation center axis, and the X-ray generator and the two-dimensional imaging means are rotated. X-rays on the two-dimensional imaging means are obtained by performing a predetermined operation on the two spherical X-ray images obtained on the two-dimensional imaging means at respective positions, inverted by 180 degrees around the central axis. Since the reference position of the original image is calculated, it is easy to create a correction table.
[0093]
  Claim2According to the X-ray image correction method according to claim 1, the correction according to claim 1.MethodEffectively corrects distortion that may occur during X-ray CT imaging while demonstrating its effectsAn X-ray original image can be corrected by an image correction table created using the center axis deviation correction table, and all possible distortions that may occur during X-ray CT imaging can be removed.
[0094]
  Claim3The X-ray image correction method described in3. The X-ray image correction method according to claim 2, wherein two spheres having an X-ray absorption coefficient higher than that of the support are accommodated in an X-ray support, and the support is formed by X-ray CT imaging. The image correction table is calculated by calculating the angle and the shift amount used for correction from the average value of the movement in the coordinate direction on the image of the two spheres. Therefore, it becomes easy to create a correction table.
[Brief description of the drawings]
1A is a front view showing an example of an image correction phantom of the present invention, FIG. 1B is a longitudinal sectional view showing a state in which the image correction phantom is mounted on a two-dimensional X-ray imaging means, and FIG. Showing an example of an original X-ray image captured by
2A to 2E are conceptual diagrams of a correction method using the image correction phantom of FIG.
FIG. 3 is an external view showing a reference example of an image correction phantom.
4 is a conceptual diagram of a correction method using the image correction phantom of FIG.
FIG. 5 is an external view showing another example of the image correction phantom of the present invention.
6 is a conceptual diagram of a correction method using the image correction phantom of FIG.
FIG. 7 is a basic configuration diagram of an example of an X-ray CT imaging apparatus using the image correction phantom of the present invention.
FIG. 8 is a detailed configuration diagram of the image correction unit in FIG. 7;
9A is an external front view of the X-ray CT imaging apparatus of FIG. 7, FIG. 9B is a side view, and FIG. 9C is a partial side view.
10A is a diagram showing an example of a conventional image correction phantom, and FIG. 10B is a basic configuration diagram of an X-ray CT imaging apparatus using the image correction phantom.
11A is a diagram showing another example of a conventional image correction phantom, and FIG. 11B is a basic configuration diagram of an X-ray CT imaging apparatus using the image correction phantom.
12A is a diagram showing another example of a conventional image correction phantom, and FIG. 12B is a basic configuration diagram of an X-ray CT imaging apparatus using the image correction phantom.
[Explanation of symbols]
      1 X-ray generator
      1a X-ray cone beam
      2 Two-dimensional X-ray imaging means
      21 X-ray image correction phantom
      21a Printed circuit board
      21b small circle
      22 X-ray image correction phantom
      22a Support
      22b Medium sphere
      23 X-ray image correction phantom
      23a Support
      23b small sphere
      3 swivel arm
      3a Rotation center axis
      4 Subject holding means
      4a chair
      4b Head fixing means
      6 Image processing means
      6a Image correction means
    20 X-ray CT imaging system

Claims (3)

被写体を挟んで対向配置させたX線発生器と2次元撮像手段とを、回転中心軸を中心に旋回させながら、X線発生器からコーンビームX線を被写体に照射するコーンビームX線CT撮影における回転中心軸の軸方向の位置ずれを補正するX線画像補正方法であって、
X線透過性の支持体内に、この支持体よりもX線吸収係数の高い2つの球体を収容して成り、
この支持体を上記回転中心軸より所定位置離れた位置に設置し、
上記X線発生器と上記2次元撮像手段とを上記回転中心軸を中心にして180度反転させ、それぞれの位置で上記2次元撮像手段上で得られた上記2つの球体のX線原画像に対して、所定の演算を行って上記2次元撮像手段上におけるX線原画像の基準位置を算出するようにしたことを特徴とするX線画像補正方法
Cone beam X-ray CT imaging that irradiates a subject with cone beam X-rays from an X-ray generator while rotating an X-ray generator and a two-dimensional imaging unit opposed to each other with the subject interposed therebetween. An X-ray image correction method for correcting a positional deviation in the axial direction of the rotation center axis in
An X-ray transparent support body accommodates two spheres having an X-ray absorption coefficient higher than that of the support body.
This support is installed at a position away from the rotation center axis by a predetermined position,
The X-ray generator and the two-dimensional imaging means are inverted 180 degrees around the rotation center axis, and the X-ray original images of the two spheres obtained on the two-dimensional imaging means at the respective positions are obtained. An X-ray image correction method characterized in that a predetermined calculation is performed to calculate a reference position of an X-ray original image on the two-dimensional imaging means.
請求項1において、
X線CT撮影の回転中心軸の軸方向の位置ずれを補正するための中心軸ずれ補正テーブルを作成し、この中心軸ずれ補正テーブルを用いてX線CT撮影の各旋回角度毎に画像補正テーブルを作成し、該画像補正テーブルを用いて、それぞれのX線原画像を補正することを特徴とするX線画像補正方法。
In claim 1,
A center axis deviation correction table for correcting the axial position deviation of the rotation center axis of the X-ray CT imaging is created , and an image correction table for each turning angle of the X-ray CT imaging is created using this center axis deviation correction table. And correcting each original X-ray image using the image correction table.
請求項2において、
更に、X線の支持体内にこの支持体よりもX線吸収係数の高い2つの球体を収容して成り、この支持体を、上記球体がX線CT撮影におけるほぼ回転中心軸方向に上下になるよう被写体保持位置に設置し、上記2つの球体の画像上の座標方向の動きの平均値から補正に用いる角度とシフト量とを算出することで、前記画像補正テーブルを作成することを特徴とするX線画像補正方法。
In claim 2,
Further, two spheres having an X-ray absorption coefficient higher than that of the support are accommodated in the X-ray support, and the support is vertically moved in the direction of the rotation center axis in X-ray CT imaging. The image correction table is created by calculating the angle and shift amount used for correction from the average value of the movement in the coordinate direction on the image of the two spheres. X-ray image correction method.
JP2001150181A 2001-05-18 2001-05-18 X-ray image correction method Expired - Fee Related JP3964155B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2001150181A JP3964155B2 (en) 2001-05-18 2001-05-18 X-ray image correction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2001150181A JP3964155B2 (en) 2001-05-18 2001-05-18 X-ray image correction method

Publications (2)

Publication Number Publication Date
JP2002336237A JP2002336237A (en) 2002-11-26
JP3964155B2 true JP3964155B2 (en) 2007-08-22

Family

ID=18995235

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001150181A Expired - Fee Related JP3964155B2 (en) 2001-05-18 2001-05-18 X-ray image correction method

Country Status (1)

Country Link
JP (1) JP3964155B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019099544A1 (en) * 2017-11-14 2019-05-23 Imatrex, Inc. Geometric calibration of x-ray imaging systems
EP4063838A1 (en) * 2021-03-26 2022-09-28 Rigaku Corporation Tolerance error estimating apparatus, method, program, reconstruction apparatus and control apparatus

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3918932B2 (en) * 2002-11-29 2007-05-23 株式会社島津製作所 Center axis calibration jig in X-ray CT apparatus and calibration method using this jig
US7570734B2 (en) 2003-07-25 2009-08-04 J. Morita Manufacturing Corporation Method and apparatus for X-ray image correction
JP2006000225A (en) 2004-06-15 2006-01-05 Canon Inc X-ray ct apparatus
JP4610304B2 (en) * 2004-11-05 2011-01-12 株式会社日立メディコ X-ray CT system
JP4881796B2 (en) * 2007-06-08 2012-02-22 日立アロカメディカル株式会社 X-ray CT system
JP5415885B2 (en) * 2009-09-30 2014-02-12 富士フイルム株式会社 Radiation CT apparatus and image processing apparatus
KR101098417B1 (en) 2011-08-11 2011-12-23 테크밸리 주식회사 Position modification jig for voxels of ct apparatus
KR101214377B1 (en) 2011-09-23 2012-12-21 한국전기연구원 Method for calibrating the focal spot position of x-ray using calibration phantom

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019099544A1 (en) * 2017-11-14 2019-05-23 Imatrex, Inc. Geometric calibration of x-ray imaging systems
US11402524B2 (en) 2017-11-14 2022-08-02 Imatrex, Inc. Geometric calibration of X-ray imaging systems
EP4063838A1 (en) * 2021-03-26 2022-09-28 Rigaku Corporation Tolerance error estimating apparatus, method, program, reconstruction apparatus and control apparatus
US11885754B2 (en) 2021-03-26 2024-01-30 Rigaku Corporation Tolerance error estimating apparatus, method, program, reconstruction apparatus and control apparatus

Also Published As

Publication number Publication date
JP2002336237A (en) 2002-11-26

Similar Documents

Publication Publication Date Title
JP4163991B2 (en) X-ray CT imaging apparatus and imaging method
JP4415762B2 (en) Tomography equipment
JP3964155B2 (en) X-ray image correction method
US9724051B2 (en) Medical X-ray CT photography apparatus using positional shift based on rotation angle
JPWO2010029862A1 (en) X-ray inspection apparatus and X-ray inspection method
WO2018159820A1 (en) X-ray tomography apparatus and x-ray tomography method
JP6837400B2 (en) X-ray imaging device and X-ray imaging method
JP6187959B2 (en) System and method for orienting an x-ray detector
JP7071410B2 (en) X-ray tomography equipment with added scanner function
KR20140044158A (en) Radiography system for dental model and jig device used thereof
JP6666283B2 (en) X-ray tomography apparatus and X-ray tomography method
JP4268996B2 (en) Local X-ray CT imaging apparatus and image display method thereof
US8213565B2 (en) Method for correcting truncated projection data
JP4573593B2 (en) X-ray image correction method and apparatus
JP4169292B2 (en) X-ray diagnostic equipment
JP3907158B2 (en) X-ray CT imaging system
JP3614214B2 (en) Tomography equipment
KR20140044174A (en) X-ray imaging apparatus and method
JP2006288719A (en) X-ray ct photographing method and device
JP2013047644A (en) Radiographing apparatus and tomographic image correction method
JP4504442B2 (en) X-ray diagnostic equipment
JP2008125909A (en) X-ray ct apparatus
JP2015100361A (en) X-ray diagnostic apparatus
CN102283661A (en) Mammary machine with three-degree-of-freedom-of-motion C-shaped arm
JP2009118950A (en) X-ray ct apparatus

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20040126

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20051205

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20060207

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060410

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20061219

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070216

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: 20070424

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20070523

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees