JP3964155B2 - X-ray image correction method - Google Patents
X-ray image correction method Download PDFInfo
- 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
Links
- 238000003702 image correction Methods 0.000 title claims description 87
- 238000000034 method Methods 0.000 title claims description 37
- 238000003384 imaging method Methods 0.000 claims description 58
- 238000013170 computed tomography imaging Methods 0.000 claims description 43
- 238000010521 absorption reaction Methods 0.000 claims description 15
- 238000010586 diagram Methods 0.000 description 19
- 239000000463 material Substances 0.000 description 8
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 5
- 238000006073 displacement reaction Methods 0.000 description 5
- 230000009466 transformation Effects 0.000 description 5
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000002591 computed tomography Methods 0.000 description 3
- 229910052751 metal Inorganic materials 0.000 description 3
- 239000002184 metal Substances 0.000 description 3
- 239000004925 Acrylic resin Substances 0.000 description 2
- 229920000178 Acrylic resin Polymers 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- RPPBZEBXAAZZJH-UHFFFAOYSA-N cadmium telluride Chemical compound [Te]=[Cd] RPPBZEBXAAZZJH-UHFFFAOYSA-N 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- 238000000844 transformation Methods 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- BZHJMEDXRYGGRV-UHFFFAOYSA-N Vinyl chloride Chemical compound ClC=C BZHJMEDXRYGGRV-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- BIJOYKCOMBZXAE-UHFFFAOYSA-N chromium iron nickel Chemical compound [Cr].[Fe].[Ni] BIJOYKCOMBZXAE-UHFFFAOYSA-N 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009429 electrical wiring Methods 0.000 description 1
- 230000003028 elevating effect Effects 0.000 description 1
- 239000003822 epoxy resin Substances 0.000 description 1
- 238000005530 etching Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 239000003365 glass fiber Substances 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 150000002739 metals Chemical class 0.000 description 1
- 229910000623 nickel–chromium alloy Inorganic materials 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000005011 phenolic resin Substances 0.000 description 1
- 239000004033 plastic Substances 0.000 description 1
- 229920003023 plastic Polymers 0.000 description 1
- 239000004417 polycarbonate Substances 0.000 description 1
- 229920000515 polycarbonate Polymers 0.000 description 1
- 229920000647 polyepoxide Polymers 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000003014 reinforcing effect Effects 0.000 description 1
- 239000010935 stainless steel Substances 0.000 description 1
- 229910001220 stainless steel Inorganic materials 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 239000000057 synthetic resin Substances 0.000 description 1
- 229920003002 synthetic resin Polymers 0.000 description 1
- WFKWXMTUELFFGS-UHFFFAOYSA-N tungsten Chemical compound [W] WFKWXMTUELFFGS-UHFFFAOYSA-N 0.000 description 1
- 229910052721 tungsten Inorganic materials 0.000 description 1
- 239000010937 tungsten Substances 0.000 description 1
- 238000004876 x-ray fluorescence Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
- A61B6/582—Calibration
- A61B6/583—Calibration 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線画像補正方法は、請求項2の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】
【発明の効果】
請求項1に記載のX線画像補正方法によれば、X線発生器からコーンビームX線を被写体に照射するコーンビームX線CT撮影における回転中心軸の軸方向の位置ずれを補正するX線画像補正方法であって、X線透過性の支持体内に、この支持体よりもX線吸収係数の高い2つの球体を収容して成り、この支持体を回転中心軸より所定位置離れた位置に設置し、X線発生器と2次元撮像手段とを回転中心軸を中心にして180度反転させ、それぞれの位置で2次元撮像手段上で得られた2つの球体のX線原画像に対して、所定の演算を行って2次元撮像手段上におけるX線原画像の基準位置を算出するようにしているので、補正テーブルを作成することが容易になる。
【0093】
請求項2に記載のX線画像補正方法によれば、請求項1に記載の補正方法を有効利用して、その効果を発揮させながら、かつ、X線CT撮影時に発生が考えられ得る歪みを補正する中心軸ずれ補正テーブルを用いて作成した画像補正テーブルによって、X線原画像を補正することができ、X線CT撮影時に発生が考えうる全ての歪みを除去することができる。
【0094】
請求項3に記載の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
[0005]
In the X-ray CT imaging apparatus shown in FIG. 10B, an
[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
[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
[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
[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
[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
[0016]
Claim 2The X-ray image correction method described
[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
[0021]
The printed
[0022]
As the X-ray opaque material of the small
[0023]
The
[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
[0025]
The two-dimensional X-ray imaging means 2 is provided with an
[0026]
As shown in FIG. 1B, the
[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-
[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
π * (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
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
[0037]
Here, the material of the
[0038]
The
[0039]
The
[0040]
Further, the
[0041]
In FIG. 4, the images of the two
[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
[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
[0047]
Here, the material of the
[0048]
The
[0049]
The
[0050]
In FIG. 6, the two
[0051]
Reference numerals a1 and b1 denote the X-ray sources of the small spheres A and B when the
[0052]
Between the four positions a1, b1, a2, b2 thus obtained and the vertical position c of the
(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
[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
[0057]
In other words, in the case of the center axis blur correction, the
[0058]
On the other hand, in the case of center axis deviation correction, it is necessary to install the
[0059]
FIG.Image correction phantomFIG. 8 is a detailed configuration diagram of the image correction unit in FIG. 7.
[0060]
The X-ray
[0061]
The X-ray imaging means A has a
[0062]
The X-ray beam adjusting means B provided in the
[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
[0066]
Further, the
[0067]
This
[0068]
Therefore, the
[0069]
A hollow portion 3 b is provided at the
[0070]
The hollow portion 3b is used for arranging a connection line between the
[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
[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
[0073]
The subject holding means 4 includes a
[0074]
The subject position moving unit 5 includes an
[0075]
The X-axis, Y-axis, and Z-axis linear movement tables (not shown) driven by these
[0076]
In this way, the subject O is seated on the
[0077]
In this
[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
[0080]
In order to achieve this object, the
[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
[0083]
The
[0084]
The members constituting the
[0085]
As described above, the
[0086]
The
[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
[0088]
The
[0089]
As can be seen from FIGS. 9 and 3, the
[0090]
As shown in FIG. 9C, the
[0091]
As can be seen from FIGS. 9 and 5, the
[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
[0094]
Claim3The X-ray image correction method described in3. The X-ray image correction method according to
[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線透過性の支持体内に、この支持体よりも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.
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.
更に、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.
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)
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)
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 |
-
2001
- 2001-05-18 JP JP2001150181A patent/JP3964155B2/en not_active Expired - Fee Related
Cited By (4)
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 |