[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

JP3825492B2 - Image reconstruction processing apparatus and X-ray CT apparatus - Google Patents

Image reconstruction processing apparatus and X-ray CT apparatus Download PDF

Info

Publication number
JP3825492B2
JP3825492B2 JP00101596A JP101596A JP3825492B2 JP 3825492 B2 JP3825492 B2 JP 3825492B2 JP 00101596 A JP00101596 A JP 00101596A JP 101596 A JP101596 A JP 101596A JP 3825492 B2 JP3825492 B2 JP 3825492B2
Authority
JP
Japan
Prior art keywords
data
ray
voxel
centering
detector
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 - Lifetime
Application number
JP00101596A
Other languages
Japanese (ja)
Other versions
JPH09187449A (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.)
Toshiba Corp
Original Assignee
Toshiba 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 Toshiba Corp filed Critical Toshiba Corp
Priority to JP00101596A priority Critical patent/JP3825492B2/en
Priority to US08/778,820 priority patent/US5838756A/en
Publication of JPH09187449A publication Critical patent/JPH09187449A/en
Application granted granted Critical
Publication of JP3825492B2 publication Critical patent/JP3825492B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)

Description

【0001】
【発明の属する技術分野】
この発明は、X線源から円錐状に放射されたX線を対象物に照射し、この対象物を透過したX線をX線検出器により検出し、このX線検出器から得られた検出データに基づいて、前記対象物の透過画像を再構成する画像再構成処理装置及びX線CT装置に関する。
【0002】
【従来の技術】
従来のX線CT(computed tomography )装置では、図23に示すように、X線源101からX線ビームがファン状( 扇形状 )に放射されるファンビームを使用するものが知られている。
このようなX線CT装置は、X線源101から放射されたX線ビームを被写体に照射し、この被写体を通過したX線を扇状に1列に約1000チャンネル配列したX線検出器102で検出してデータ収集を行い、X線源101及びX線検出器102を被写体の周囲を回転させながら、1回転する間に1000回程度データ収集し( 1回のデータ収集を1ビューと称する )、その収集されたデータに基づいて被写体のX線の透過画像を再構成する。なお、FOV103は、有効視野を示すものである。
このファンビームを使用したときの画像再構成式は、( 式1 )により算出される。
【0003】
【数1】

Figure 0003825492
【0004】
この( 式1 )から判るように、ファンビームの再構成では、X線検出器から得られたデータに、再構成すべきピクセルの位置に依存した重み付けを乗算して逆投影する必要があるので複雑な処理になる。
すなわち、図24に示すように、有効視野FOV103に対して再構成すべき画像を構成するピクセルが設定されており、X線検出器102の各チャンネルで得られたデータを、重み付けとして焦点( X線源101のX線ビームの放射点 )−ピクセル間距離 FpixelD(X) の2乗の逆数を乗算して、該当するピクセルに逆投影する。なお、FpixelD は、Focus-Pixel-Distanceである。
また、直接逆投影する方法もあるが、この場合には極座標変換が必要となり複雑な計算になる。
【0005】
そこで、ファンビームを使用したX線CT装置では、現在のところ2種類の画像再構成法が考案されている。
1つの方法は、ファン−パラ変換法と呼ばれるものであり、これは、図25に示すように、ファンビームによるX線検出器から得られた投影データを並び替えかつ補間してパラレルビーム投影データを作成( 変換を含む )し、これにより得られたデータを、従来のパラレルビームを使用したX線CT装置で行われるように逆投影する方法である。
データ変換の計算と補間処理などが必要になる反面、逆投影時には、ファンビームのときの再構成ピクセル毎に異なった重み付け処理などが不要で、1つのデータ( パラレルビーム投影データ )をビーム路( パラレルビームとなるときの放射点とX線検出器のチャンネルとを結ぶ直線 )の全てのピクセルに逆投影すれば良いので処理が単純になる。
【0006】
他の1つの方法は、センタリング軸を使用したファンビーム再構成法であり、その詳細は特開昭55−99240号に開示している。
図26に示すように、一度ピクセル列に平行な所定の基準軸( センタリング軸 )104にX線検出器から得られたデータを射影( 逆投影 )し、それを再度再構成のピクセル列毎に逆投影するものである。このように一度センタリング軸に逆投影することで、ピクセル毎に異なる重み付け処理がピクセル列毎には同一となるので、高速かつ単純な処理が可能になる。
【0007】
このセンタリング軸を使用したファンビーム再構成法のステップを以下に説明する。
1.投影データData-Proj のX 線強度補正などの種々の補正とcos 項の乗算して、生データData-Rawを得る。この生データData-Rawと再構成関数とのコンボリューション演算をしてData-Conv を得る。
2.Data-Conv をある基準軸 (センタリング軸、例えばピクセルが配列された基準となるX軸及びY軸)上の点に( 式2 )の重み付け処理して射影し、Data- Centerを得る。この( 式2 )式において、FcpD(X) は、焦点−センタリング軸点間距離である。
【数2】
Figure 0003825492
【0008】
3.2の処理を繰り返し、全ビューのファンビーム投影データを対応する基準軸に射影する。
4.基準軸に射影されたあるビューのデータData-Center を、再構成する画像の全ピクセル列に対しピクセル列毎に同一の( 式3 )の重み付けしてData-Back を得る。この得た射影データData-Back を、逆投影( 画像メモリのピクセルに相当するアドレスに加算) する。なおA、Bは、図26に示す。
【数3】
Figure 0003825492
【0009】
5.4の処理を全ビュー繰り返し全ビューの射影データを逆投影する。
【0010】
なおここで、図26を参照して( 式3 )を変形すると、
【数4】
Figure 0003825492
【0011】
となり、これはファンビームの再構成式、( 式1 )の積分の中身と一致する。
【0012】
センタリング軸に一度射影した後に各pixel 列に逆投影することで、本来ピクセル毎に異なっていた逆投影時の重み( 式1 )を、 pixel列単位では等しい重み( 式2 )にして,重みの発生の計算自体の簡便化と計算回数の削減を達成している。
また、ここでは詳細には述べないが、円弧上に等角度で配列された検出器の素子と直線上に等ピッチで配列されたpixel との複雑な対応関係をも簡略化している。
従来のファンビームを使用したX線CT装置では、以上説明した2つの方法のうちいずれかにより、高速な画像再構成処理を実現している。
【0013】
さらに一方、図27に示すように、X線源201からX線ビームが円錐状に放射されるコーンビームと、ファンビーム用検出器列をZ軸方向にN列積み重ねたような、円筒面上に検出器の素子(Mチャンネル×N列)を配列した2次元X線検出器202とを使用して、X線透視画像を撮影するX線CT装置が考案されている。
このようなコーンビームを使用したX線CT装置における代表的なコーンビーム再構成(Feldkamp再構成)は、下記の文献に開示されている。
【0014】
"Practical cone-beam algorithm"
L.A.Feldkamp, L.C.Davis, and J.W.Kress
J. Opt. Soc. Am. A/Vol.1, No.6, pp.612-619/June 1984
これは、数学的に厳密な再構成法であるファンビーム(2次元平面内) 再構成アルゴリズム[ Filtered-Backprojection(フィルタ補正逆投影法) ] を、Z軸方向に拡張することによって得られた近似的な3次元再構成アルゴリズムである。
【0015】
このコーンビーム再構成法では、コーンビームによるコンベンショナルスキャンを対象としており,以下のステップからなる。なお、このコーンビームでは、2次元的な画素としてのピクセルの代わりに、図28に示すように、3次元的な画素としてのボクセルが使用される。
【0016】
1.投影データの重み付け
投影データに、Z座標に依存した項とcos 項を乗算する。
【0017】
2.コンボリューション演算
1の処理により得たデータと、ファンビームと同じ再構成関数とのコンボリューション演算を行う。
【0018】
3.BackProjection(逆投影)
2の処理により得たデータを、X線が通過した( 焦点から検出器のチャンネルまでの) パス上に逆投影する。すなわち、焦点から逆投影するボクセルを通る直線が検出器面と交差する点を計算し、その点の周囲の2の処理のデータから逆投影するデータを補間などで作成し、それをFvoxelD(X)の2乗の逆数で重み付けして逆投影する。この逆投影は360°( 1回転 )にわたって行なう。
【0019】
ファンビーム再構成式と類似な式で表現すると、下記となる。
なお、 FvoxelD(X)=Focus-Voxel-Distanceは、焦点−ボクセル間距離である。
【数5】
Figure 0003825492
【0020】
この3次元再構成式( コーンビーム再構成式 )( 式4 )について、式上ではファンビーム再構成と非常に似ているが、Data-Back の逆投影方法が大きく異なることを説明する。
2次元的なファンビーム再構成においては、図29に示すように、再構成面内の全画素(ピクセル)に対して1 次元に配列された検出器のデータから逆投影するのに対し、コーンビーム(Feldkamp)再構成においては、図30に示すように、焦点と再構成するボクセル(voxel) を結んだ直線が2次元のX線検出器面と交差する点を求め、その交差点に関与する検出器素子から得られるデータをその直線上に位置する全てのボクセルに逆投影する。
【0021】
従って、コーンビーム再構成で、ファンビーム再構成のようにある面を再構成する場合には、特定の検出器列かつチャンネルのデータが再構成面の一部のボクセルにのみ逆投影されるため、各ボクセルに対して逆投影するデータ(検出器列と検出器チャンネル)を選択する必要があるので、再構成ボクセルと焦点を結んだ直線とX線検出器面の3次元的な位置関係が重要になる。
しかも、Z座標が同じ検出器列を考え、その検出器素子と焦点を結んだ直線を考えた場合、ある面( 再構成面 )においてそれらの直線が通過するボクセルは、焦点を中心とした検出器面の相似図形(円筒検出器の場合、同心円)上に並ぶため、この位置関係の計算は非常に複雑になる。
【0022】
【発明が解決しようとする課題】
上述したように、従来のX線CT装置において、コーンビーム及び2次元的なX線検出器を使用した場合、画像の再構成ではその計算が複雑で膨大な量になり、一般的に普及しているコンピュータ等では処理時間が長くかかり過ぎて実現できないという問題があった。
そこでこの発明は、コーンビームを使用して撮影された画像の正確な再構成を実現することができる画像再構成処理装置及びX線CT装置を提供することを目的とする。
【0023】
【課題を解決するための手段】
本発明の第1局面は、X線源から円錐状に放射されたX線を対象物に照射し、この対象物を透過したX線を、複数の検出素子を2次元状に配列してなるX線検出器により検出し、このX線検出器から得られた検出データに基づいて、前記対象物の画像を再構成する画像再構成処理装置において、前記複数の検出素子で検出された検出データを、ボクセルの特定面と略平行にボクセル集合体内に設定されたセンタリング面上のデータに置き換える手段と、前記センタリング面に置き換えたデータを、前記ボクセル集合体の各ボクセルに対して逆投影する手段とを有することを特徴とする。
本発明の第2局面は、X線を発生するX線源と、被検体を透過したX線を検出する複数の検出素子が前記被検体の体軸方向に複数配列されたX線検出器と、前記X線検出器からの出力信号に基づいて逆投影するデータを得る手段と、前記被検体の画像を再構成する画像再構成処理手段を有するX線CT装置において、前記画像再構成処理手段は、3次元的に空間配置された画像であるボクセルと、前記X線源との3次元的な位置関係に基づいて、当該ボクセル列を前記X線検出器に投影した投影曲線を計算し、前記逆投影するデータを特定する手段と、前記特定した逆投影するデータを、前記ボクセルに逆投影する手段とを備えたことを特徴とする。
【0025】
【発明の実施の形態】
以下、この発明の第1の実施の形態を図1〜図16を参照して説明する。
図1は第1の実施の形態によるX線CT装置の構成図である。図2は、図1のガントリの外観図である。
投影データ測定系としてのガントリ(架台ともいう)1は、円錐に近似したコーンビーム状のX線束を発生するX線源3と、複数の検出素子を2次元状に配列してなる2次元アレイ型のX線検出器5とを収容する。前記X線源3と前記X線検出器5とは、寝台6のスライド天板に載置された被検体を挟んで対向した状態で回転リング2に装備される。
前記X線検出器5としては、複数( 1000チャンネル )の検出素子が、ファンビーム用の1次元的に配列された1次元アレイ型検出器を複数列( 10列 )積み重ねられたように配列されて構成されたもの( 図27参照 )で、前記回転リング2に実装される。ここで、1つの検出素子は1チャンネルに相当するものと定義する。
【0026】
前記X線源3からのX線はX線フィルタ4を介して被検体に曝射される。被検体を通過したX線は前記X線検出器5で電気信号として検出される。
X線制御器8は高圧発生器7にトリガ信号を供給する。この高圧発生器7はトリガ信号を受けたタイミングで前記X線源3に高電圧を印加する。これによりX線源3からはX線が曝射される。
架台寝台制御器9は、前記ガントリ1の前記回転リング2の回転と、前記寝台6のスライド天板のスライドとを同期して制御する。システム全体の制御中枢としてのシステム制御器10は、被検体から見て前記X線源3が螺旋軌道を移動するいわゆるヘリカルスキャンを実行するように、前記X線制御器8と前記架台寝台制御器9を制御する。
具体的には、前記回転リング2が一定の角速度で連続回転し、前記寝台6のスライド天板が一定の速度で移動し、前記X線源3から連続的又は一定角度毎に間欠的にX線が曝射される。
【0027】
前記X線検出器5からの出力信号は、チャンネル毎にデータ収集部11で増幅され、ディジタル信号に変換される。このデータ収集部11から出力される投影データは、再構成処理部12に取り込まれる。
この再構成処理部12は、投影データに基づいてボクセル毎にX線吸収率を反映した逆投影データを求める。
コーンビームを使用したヘリカルスキャン方式のX線CT装置において、有効視野( FOV、撮影領域 )は、ヘリカルスキャンの回転中心軸を中心として円筒形状となり、再構成処理部12は、この有効視野に複数のボクセル( 3次元的に配置された画素 )を規定し( 図28参照 )、X線検出器5からの投影データから各ボクセルの逆投影データを求める。この逆投影データに基づいて作成された3次元画像データ又は断層像データは表示装置14に送られ3次元画像又は断層像としてビジュアルに表示される。
【0028】
図3(a)及び図3(b)に示すように、このX線CT装置のジオメトリは、
検出器列数 M=20、
各列のZ軸方向の高さ Dseg =2mm、
X線検出器の厚み M×Dseg =40mm、
チャンネル数 N=1000、
焦点−回転中心間距離 FCD(Focus-center-Distance )=600mm、
焦点−検出器間距離 FDD(Focus-Detector-Distance )=1200mm、
有効視野直径 FOV(Field of View )=500mm、
有効視野角(ファン角)θ=50°
となっている。
【0029】
図4は、前記再構成処理部12の要部構成を示すブロック図である。
21は、コンボリューション処理及び逆投影処理におけるデータ選択、重み付け、センタリング処理、逆投影などの計算及び3次元再構成処理全体を制御するこの再構成処理部12の制御部本体を構成する再構成処理制御部である。
コンボリューション演算部22は、前記データ収集部11で収集された投影データをコンボリューション処理し、このコンボリューション処理により得られたコンボリューションデータは第1のデータメモリ23に記憶される。
【0030】
第1の逆投影部24は、前記第1のデータメモリ23に記憶されたコンボリューションデータを予め設定されたセンタリング面に逆投影( 射影 )処理し、この逆投影されたセンタリングデータは第2のデータメモリ25に記憶される。
【0031】
第2の逆投影部26は、前記第2のデータメモリ23-2記憶されたセンタリングデータをボクセルに逆投影( 3次元逆投影 )処理し、この逆投影された再構成データ( 画像 )は画像メモリ27に記憶される。
【0032】
ここで再構成についてコーンビームによる3次元再構成について、幾つかの点について考察する。
再構成ボクセル列( 直線 )、センタリング面( 平面 )、検出器列( 円弧 )、検出器面( 円筒、但しデータメモリ上のように展開して考えるときは平面である) などの座標変換を伴う関係について説明する。
まず、再構成ボクセル列、検出器面、センタリング面を考える。
図5に示すように、1列のボクセル列( 直線 )の検出器面への投影が図6に示すような曲線になることを説明する。
【0033】
再構成ボクセル列、検出器面、センタリング面をZ軸方向( 上方 )から観察した図を、図7( a )に示す。X線検出器5のチャンネル( 検出器列 )は焦点FからZ軸方向において等角度の円弧上に配列されており、図のように焦点回転( ビュー) による回転角をφ、チャンネル方向の角度をθで示す。
また、この焦点位置においてセンタリング面はX軸上にあり( X軸Z軸平面に含まれている )、X軸はセンタリング面のX軸であるXcp軸に一致する。再構成ボクセル列の座標を( Xv ,Yv )…( Xv はボクセルと共に変化するが、Yv は一定 )で定義する。
図7( b )に示すように、線分FCの長さ( 焦点−回転中心間距離 )をFCD、線分FV0 の長さ( 焦点−ボクセル間距離 )をFCD´、あるボクセルVを考えるときに線分FVをXY平面へ射影したときの( 図の点線における )長さを FCD”、焦点- 検出器間距離をFDDとする。また、焦点とボクセルを通る直線を引き、XY平面上での焦点からの距離がFCDのときのZ座標をZ0 、ボクセルVのZ座標をZv とする。さらに、展開した検出器面上での座標を図7( c )のように横軸( Xdet=θ・FDD )と縦軸( Zdet)とする。
【0034】
まず、再構成ボクセル列を検出器面に投影することを考える。なお、Xcpは線分FVのセンタリング面上でのX座標を示し、次の( 式5 )により求められる。
【数6】
Figure 0003825492
【0035】
そこで、△FCXcpと△FVo Xv は相似なので、次に示す関係式( 式6 )を使用して、
【数7】
Figure 0003825492
【0036】
となる。この( 式7 )がXv →θ変換、再構成ボクセル列の検出器面への投影におけるチャンネル方向の式である。
【0037】
さて、線分FXcpと線分FCDとの間の関係は次に示す( 式8 )であるから、( 式9 )を得る。
【数8】
Figure 0003825492
【0038】
従って、再構成ボクセル列の検出器面への投影における列方向の式は、次に示す( 式10 )となる。
【数9】
Figure 0003825492
【0039】
( 式7 )から等ピッチに並んだ再構成ボクセル列を検出器面へ投影すると非線形な配列になることがわかる。従って、図中左方例えば第1ボクセルと第2ボクセルを検出器面上に投影したときの間隔(例えば2.4チャンネル分)と、図中右方例えば第511ボクセルと第512ボクセルの間隔(例えば3.5チャンネル分)とは異なり、角度θに依存して非線形な配列になる。
また、( 式10 )からZ座標がZv で固定された直線の再構成ボクセル列を検出器面へ投影すると角度θに依存して図6のような非線形な曲線になることがわかる。ただし、再構成ボクセル列のZ座標Zv とZdet とは比例関係にある。
【0040】
これを示すのが図8及び図9であり、図8には、ボクセル列V、センタリング面C、検出器面Dとそれぞれの変数および端点と中心点の定義を示す。なお、図10には、ボクセル列Vとセンタリング面Cとを示す。このボクセル列Vに対してセンタリング面Cは平行に配置されている。
図9( b )に示した直線であるボクセル列のセンタリング面への投影も、図9( c )に示すように直線になる。ボクセルのピッチとセンタリング面上の点のピッチは一定の比になっており、歪みは全く生じない。すなわち、逆投影すべきデータが、センタリング面上で直線かつ等ピッチであれば、当然そのボクセルへの逆投影も直線かつ等ピッチとなる。すなわち、ボクセルとセンタリング面との関係は単純な拡大縮小関係となる。
【0041】
しかし、ボクセル列を検出器面に投影すると、図9( a )に示すように、チャンネル方向、列方向ともに非線形な歪みが発生する。しかし、Z座標Zv の異なる2本のボクセル列の投影が示すように、検出器面上の2本のボクセル列の投影像どうしの間には、上述の拡大縮小関係的な比例関係が成立している。
【0042】
以上の考察をふまえて、最も単純なコーンビームの再構成法は、以下に説明するステップである。
1.X線検出器からの投影データをX線強度補正等の補正処理後、Feldkamp重みづけ処理し、データメモリに記憶させる。
2.データメモリに記憶された補正投影データを再構成関数とコンボリューション処理し、データメモリに記憶させる。
3.( 式7 )および( 式10 )に従って逆投影するボクセル列を検出器面に投影した投影曲線を計算し、逆投影するデータを選択してそのアドレスを発生させる。
4.該当するデータを読み出し、所定の重みづけ処理後、画像メモリ該当するボクセルの位置に加算する。
【0043】
これを“直接逆投影法”と称する。3の処理は、予め投影曲線を計算しておきテーブルとしてデータメモリ等に記憶しておいても良い。これを“テーブル法”と称する。
いずれの方法でも投影曲線は近似曲線でも良い。しかしこの方法ではファンビーム投影データの2次元逆投影でも存在した、逆投影する (チャンネル方向の) データ選択の計算と、重み計算の他に、列方向のデータ選択計算が加わるため、膨大な計算量になってしまう。または“テーブル法”による場合には巨大なテーブルを記憶するメモリが必要になる。
【0044】
そこで、この発明の第1の実施の形態においては、X線検出器5から得られた逆投影するデータ( 投影データ )を一度センタリング面へ逆投影後、各ボクセルへ逆投影する方法を行う。
図9における考察とは逆に、検出器面上(あるいはデータメモリ上)のピッチおよび直線がセンタリング面に投影される場合を検討する。前述の( 式5 )と、( 式10の変形 )を( 式8 )により解いた( 式11 )を次に示す。
【数10】
Figure 0003825492
【0045】
また、センタリング面上の直線が検出器面上に投影される場合は、次の( 式7の変形 )及び( 式11の変形 )となる。
【数11】
Figure 0003825492
【0046】
従って( 式7 )と( 式10 )に従った非線形歪みの代わりに、( 式5 )と( 式11 )に従った非線形歪みが発生する。
【0047】
これを示すのが図11である。図11( b )に直線で示された検出器列の全チャンネルのデータをセンタリング面に投影すると、図11( c )に示すように、上述の非線形( 式4及び式8 )によってXcp、Zcp方向共に歪む。
Z方向の歪み及びX方向の等ピッチの配列の歪みは、共に図9の場合とは逆になる。前述のようにセンタリング面と再構成ボクセル列との間に歪みはないが、図11( a )に示すように、検出器列の全チャンネルのデータをボクセルに投影すると同様の歪みが発生する。すなわち、図11( c )に示す前述のセンタリング面への投影のZ方向の歪み及びX方向の等ピッチの配列の歪みと図11( a )に示すボクセルへの投影のZ方向の歪み及びX方向の等ピッチの配列の歪みとは、同様のもので拡大・縮小関係になっている。
【0048】
そこで、センタリング面上で補間処理などによってデータのリサンプリングを行い、直線上に等ピッチでデータが並ぶように処理する。その結果が図12( a )である。
センタリング面上の座標系Xcp、Zcpと再構成ボクセルの座標系Xv 、Zv の関係は、( 式5 )、( 式6 )、( 式9 )、( 式10 )を応用して、次の( 式12 )及び( 式13 )を得る。
【数12】
Figure 0003825492
【0049】
従って、再構成ボクセルにおいてスライス位置Z=Zcpのアキシャル断面を再構成するとき、そのアキシャル断面を図12( b )に示すような正方形に内接する円形FOVと考えると、( 式12 )、( 式13 )でZv を断面内で常に定数、ボクセル列内でYv を定数とし、ボクセル単位でXv を変化させることになる。従って、( 破線で示す )正方形とその内接する円形FOVの投影は、( 破線で示す )台形とそれに内接する円の変形( 図示せず )になる。
また、ボクセル列をセンタリング面に投影した直線は、ボクセル列の位置に対応してセンタリング面内を上下に平行移動する。そこで、逆投影するボクセル列に対応するZ座標Zcpを( 式13 )で求め、逆投影するボクセルに対応するX座標Xcpを( 式12 )で求め、対応するデータを目的のボクセルに逆投影する。これを全ボクセルに全ビュー繰り返して逆投影を行う。
【0050】
式で表現すると下のようになる。なお、FdpDは焦点−検出器素子間距離であり、FcpD( X,Z )は焦点−センタリング点間距離である。
【数13】
Figure 0003825492
【0051】
ここで、( 式14 )は、コンボリューション処理したデータを示す式であり、( 式15 )は、センタリング処理したデータを示す式であり、( 式16 )は、ボクセルに逆投影したデータを示す式である。
【0052】
この( 式16 )は、さらに
【数14】
Figure 0003825492
【0053】
となり、この( 式17 )は3次元再構成式( 式3 )と一致する。
【0054】
すなわち、この発明のセンタリング面を使用した3次元再構成法( コーンビーム再構成法 )のステップを以下に説明する。
1.再構成処理制御部21は、データ収集部11からの投影データData-Proj をX線強度補正などの補正処理して生データData-Rawを得て、第1のデータメモリ23に記憶させる。
2.再構成処理制御部21はコンボリューション演算部22により、第1のデータメモリ23のデータを読み出し、Feldkamp重みづけ( ( 式14 )の第1項 )処理後、再構成関数Conv-Function とコンボリューションし( ( 式14 )の第2項) 、第1のデータメモリ23に記憶させる。
【0055】
3.再構成処理制御部21は第1の逆投影部24により、第1のデータメモリ23に記憶されたコンボリューションデータに基づいて、次の(i) あるいは(ii)のいずれか一方により( 式15 )のセンタリング処理を行う。
(i) ( 式5 )及び( 式11 )に従って、検出器面のデータをセンタリング面に投影した投影曲線を計算し、コンボリューションデータに重み付けを行った後にセンタリング面に射影してセンタリングデータを計算し、更にセンタリング面上のセンタリングデータを図12( a )のように格子状にリサンプリングして、第2のデータメモリ25に記憶する。
(ii) ( 式7変形 )及び( 式11変形 )に従って、センタリング面上で(i) でリサンプリングしたような格子状のデータの位置を検出器面に投影した投影点の位置を計算し、投影点周囲の4個 (2列×2CH) の検出器素子のコンボリューションデータを重みづけ加算後に( 式15 )の3次元逆投影における重み付けを行ってセンタリング面に射影してセンタリングデータData-Center を計算し、第2のデータメモリ25に記憶する。
【0056】
4.再構成処理制御部21は、( 式12 )及び( 式13 )に従って、再構成するボクセル (列) をセンタリング面に投影した投影点 (直線) を計算し、逆投影するデータを例えば4個 (2列×2CH) 選択してそのアドレスを発生させる。 5.再構成処理制御部21は第2の逆投影部26により、第2のデータメモリ25から該当するセンタリングデータData-Center を読み出し、データ数が複数の場合は重みづけ加算し、2次元ファンビーム再構成におけるセンタリングデータの逆投影時と同様に、A/Bの2乗の重みづけ処理後(この正当性は( 式17 )にて証明済み)、画像メモリ27の該当するボクセルの位置に加算する。
以上で3次元再構成ができる。
この3次元再構成処理の流れを図13に示す。
【0057】
なお、上述した3の処理及び5の処理における検出器面上の点又はセンタリング面上の点からセンタリング面上の点又はボクセルへの逆投影( 射影 )データの計算では、4点Bi-Linear 補間等の線形補間やSpline補間等の非線形補間、あるいはその他の補間を使用してリサンプリングしても良いものである。さらに、補間を行わずに、例えばNearest Neighborとして、該当する投影点に最も近い点を選択して逆投影データを計算しても良い。
【0058】
なお、4点Bi-Linear 補間においては、1つの投影点に対して例えばj列及びj+1列とnチャンネル及びn+1チャンネルの4点Data(j,n) ,Data(j,n+1) ,Data(j+1,n) ,Data(j+1,n+1) が補間するデータとして計算対象となる。
このとき、4点Bi-Linear 補間においては、以下の計算が行われる。
Data(j,n) ×wch+Data(j,n+1) ×( 1−wch )…CH(j)
Data(j+1,n) ×wch+Data(j+1,n+1) ×( 1−wch )…CH(j+1)
CH(j) ×wro+CH(j+1) ×( 1−wro )…SEG(j,j+1)
以上の処理を複数の投影点について順番に繰り返して処理することになる。すなわちCH(j) を計算し、次にCH(j+1) を計算し、それらの計算結果によりSEG(j,j+1) を計算する。そして次の投影点について、CH(j) を計算し、次にCH(j+1) を計算し、それらの計算結果によりSEG(j,j+1) を計算する。以下同様にして各投影点について、繰り返して補間処理が行われる。この繰り返し補間処理は時間がかかる。
そこで、CH(j) とCH(j+1) とをそれぞれ並列処理し、この並列処理を連結するようにSEG(j,j+1) をパイプライン処理することにより、処理時間を短縮して、2点補間とほぼ同じ処理時間で補間処理を行うことができる。
【0059】
また、センタリング列数がボクセル列数以上の場合、Nearest Neighborとして該当する投影点に最も近い点を選択して逆投影データを計算することにより、画像劣化が少なく、計算回数が4点Bi-Linear 補間に比べて約1/3になり処理時間が短縮されるという効果が得られる。
【0060】
ここでは省略したが、シングルスライスCT( ファンビームCT )と同様にチャンネル方向の非線形と逆投影時の重み発生の簡便化だけを図る目的で、3の処理でのセンタリング面上でのリサンプリング処理を省略し、4の処理で逆投影するボクセルに対応するセンタリングデータに該当するような投影曲線を発生させても良い。リサンプリング処理を省略する代償として逆投影時の投影曲線が複雑になるが、トータルの補間回数が減るというメリットがある。
【0061】
このようにこの第1の実施の形態によれば、3次元逆投影において、一度センタリング面に射影することで、円弧上に等角度で配列された検出器CH列と直線上に等ピッチで配列されたボクセル列の複雑な対応関係を簡略化している。
すなわち、一度センタリング面に射影することで逆投影のときのボクセル単位に異なる重みをボクセル列単位に同じにでき、重みの計算回数削減、計算自体を簡便化した。
コーンビームCT特有の逆投影時のZ方向の歪み(非線形)を、センタリング面上で格子状にデータを配列することで是正し、再構成ボクセルへの逆投影を簡便にできた。このような是正は、格子状センタリング点への投影又はリサンプリング処理で実現した。
【0062】
以上により、第2の逆投影部26の処理を簡便化し、現存の2次元逆投影用( ファンビーム再構成 )の逆投影部の構成と比べても、大幅な拡張をせずに単純な構成で3次元逆投影における画像再構成が実現できる。
【0063】
なお、この第1の実施の形態においては、円筒型X線検出器を持つコーンビームX線CT装置について記述したが、この発明はこれに限定されるものではなく、例えば、図14に示すような平面型X線検出器を持つコーンビームCT装置でも同様の効果が得られるものである。
すなわち、図15に示すように、有効視野FOV中の直線矢印及び破線矢印で示されたボクセル列について、平面型X線検出器面への投影では、図16に示すように、傾いた直線矢印及び破線矢印になるのが、チャンネル方向( 横方向 )の歪み( 非線形な歪みを含めて )は発生しない。さらに、列方向( 縦方向 )についても歪みは発生しないが、その投影線が傾斜した直線となるため、円筒型X線検出器の場合に比べて軽減されるものの、データ選択及び逆投影時の重み発生は複雑である。
【0064】
そこで、仲介的にセンタリング面への逆投影( 射影 )するセンタリング処理を行うことによって、図17に示すように、ボクセル列はセンタリング面へ平行に投影されるので、有効視野FOV中の直線矢印及び破線矢印で示されたボクセル列の平面型X線検出器面への投影すると、ボクセル列単位での逆投影データ計算に使うセンタリング列を固定できるので、データ選択、重み発生を共に単純にすることができる。
なお、逆投影時の重みについて、ボクセル単位の計算からボクセル列単位の計算に削減できることは円筒型X線検出器の場合と同様である。また、センタリング列数を増やせば、補間精度を向上する等の効果も円筒型X線検出器の場合と同様である。
【0065】
この発明の第2の実施の形態を図18〜図22を参照して説明する。なお、この第2の実施の形態におけるX線CT装置( コーンビームCT装置 )の基本的な構成は上述した第1の実施の形態とほとんど同一なので、ここではその説明は省略し、異なる点としてのが画像再構成処理部における処理について説明する。
【0066】
補間を繰り返すと画像に生じるボケが増幅してしまう。しかし投影曲線と重みの簡単化のために、センタリング面に一度逆投影後にボクセルへ逆投影する場合、2回の補間は避けられない。
すなわち、図18に示すように1回補間の場合に、ボクセルデータDBが、検出器列をボクセルへ逆投影した( オリジナルの )データ位置( D1 ,D2 ,D3 )に対して、D1とD2との間を8:1に分ける位置にあるとき、
DB=( 1/9 )×D1+( 8/9 )×D2
となり、データ位置D1とD2とにより確定し、D3等の他のデータの干渉が完全に排除されている。
一方、図19に示すようにセンタリング面に一度逆投影する2回補間の場合には、まず、センタリング点Dc1及びDc2が、検出器列をセンタリング面へ逆投影( 射影 )した( オリジナルの )データ位置( D1 ,D2 ,D3 )に対して、D1とD2との間を5:4に分ける位置及びD2とD3との間を2:1に分ける位置にあるとき、
Dc1=( 4/9 )×D1+( 5/9 )×D2
Dc2=( 1/3 )×D2+( 2/3 )×D3
となり、さらに、ボクセルデータ( 点 )DBが、Dc1とDc2との間を3:5に分ける位置にあるとき、
Figure 0003825492
となる。すなわち、実際にはボクセルデータDBは、D1とD2との間にあるにもかかわらず、2回の補間によりD3の干渉項が加わりその分だけD1及びD2の項も誤差を含むようになり、実際の画像においてボケが増幅することになる。そこで、補間の精度を向上することにより画像に生じるボケを減少させることが必要になる。
【0067】
補間精度を向上する方法としては、次に示す2つの例がある。
図20は高精度な2回補間の第1の例である。
検出器列をセンタリング面へ逆投影( 射影 )した( オリジナルの )データ位置( D1 ,D2 ,D3 )と一致するように、センタリング面上でセンタリング点( Dc1,Dc4,Dc6 )を設定すると共に、それらの各センタリング点の間に等ピッチで補間点( Dc2,Dc3,Dc5 )を設定したことである。これにより2回補間のボケを減少させることができる。
例えば、ボクセルデータDBが、Dc3とDc4との間を2:1に分ける位置にあるとき、
Figure 0003825492
となり、オリジナルのデータ位置においてD3の干渉項を完全に排除している。しかし、センタリング面上における検出器列の投影曲線は非線形な歪みを持ち、隣接するチャンネル間すなわちXcp方向に隣接したセンタリング点におけるオリジナルの検出器データが存在するZcp座標が微妙に食い違ってくるため、全部のセンタリング列でオリジナルのデータ位置に補間点( センタリング点 )を設定することは困難である。
【0068】
図21は、高精度の2回補間の第2の例である。
等ピッチの補間点( センタリング点 )の個数を計算時間とのバランスから可能なだけ多数にしたもので、オリジナルのデータ位置には必ずしも補間点は存在しなくても良い。しかし、少なくとも検出器列の個数よりセンタリング列の個数を多くしたことにより、大部分の位置において補間精度は高く、補間によるボケが生じるのは第2回目の補間位置( ボクセル )がオリジナルのデータ位置を挟む第1回目の補間点2点の間にあるとき(例えばDc9とDc10 の間)だけである。
【0069】
例えば、ボクセルデータDBが、Dc8とDc9との間を1:1に分ける位置にあるとき、
Figure 0003825492
となり、オリジナルのデータ位置においてD3の干渉項を完全に排除している。また、ボクセルデータDBが、Dc9とDc10 との間を1:1に分ける位置( 丁度D2に対応する位置 )にあるとき、
Figure 0003825492
となり、D3( この場合ではD1も該当する )の干渉項が排除されていないが、D1,D3の重み( 係数 )がD2に比べて小さいので、従来のボケよりは減少されており、しかもこのボケが発生する範囲は、上述したように、Dc9とDc10 との間の距離に限定される。
【0070】
これをセンタリング面に応用した例を図22に示す。図22( a )は検出器列5列のセンタリング面への投影曲線であり、オリジナルのデータ位置を示す。前述の第1の実施の形態において説明したセンタリング面を使用した3次元再構成法( コーンビーム再構成法 )のステップに記載されているように、このオリジナルデータからセンタリング面において格子状に配列されたセンタリングデータを作成する。
格子の列(センタリング列)の数を検出器列数と同じ数の5列とした場合を、図22( b )に示す。
あるボクセル列のセンタリング面への投影直線Sに対して、2列目のセンタリング列C2は軽い重みで補間に使用され、3列目のセンタリング列C3は重い重みで補間に使用される。
この補間に使用されるデータの幅がZcp方向に広くなり、ボケが生じることが分かる。これは図19の2回補間の例に相当する。
【0071】
さて、センタリング列の数を検出器列数より多くした例を図22( c )に示す。この図22( c )では10列としている。
あるボクセル列に逆投影するデータ位置を上と同様に直線Sで示した。4列目のセンタリング列C4は軽い重みで補間に使用され、5列目のセンタリング列C5は重い重みで補間に使用される。
この補間に使われるデータの幅がZcp方向に狭くなり、ボケを抑制できることが分かる。これは図22の高精度な2回補間の第2の例に相当する。
なお、図22において補間方法は同様に距離の逆比による線形1次補間である。しかし例に示した矢印の点( Dc8とDc9との間にDBが位置する場合 )の補間精度は非常に高い。
【0072】
また、検出器列数に対してセンタリング列数を格段に多くして(例えば検出器列5列に対してセンタリング列50列あるいは500列)、第2回目の補間を0次補間すなわち最も近いセンタリング列を選択するNearest Neighborにしても良いものである。
【0073】
このようにこの第2の実施の形態によれば、センタリング面を使用した逆投影処理を行うときに避けられない補間によるデータのボケを、補間精度を高めることにより最小に抑制することができる。
【0074】
【発明の効果】
以上詳述したようにこの発明によれば、コーンビームを使用して撮影された画像の正確な再構成を実現することができる画像再構成処理装置及びX線CT装置を提供できる。
【図面の簡単な説明】
【図1】この発明の第1の実施の形態のX線CT装置を示すの構成図。
【図2】同実施の形態のX線CT装置のガントリを示す外観図
【図3】同実施の形態のX線CT装置のジオメトリを説明するための図。
【図4】同実施の形態のX線CT装置の再構成処理部の要部構成を示すブロック図。
【図5】同実施の形態のX線CT装置の再構成処理部におけるボクセル列の検出器面への投影を示す図。
【図6】同実施の形態のX線CT装置の再構成処理部におけるボクセル列の検出器面上の投影曲線を示す図。
【図7】同実施の形態のX線CT装置の再構成処理部におけるボクセル列の検出器面上の投影曲線を説明するための図。
【図8】同実施の形態のX線CT装置の再構成処理部におけるボクセル列、センタリング面、検出器面とそれぞれの変数および端点と中心点の定義を示すセンタリング面を示す図。
【図9】同実施の形態のX線CT装置の再構成処理部におけるボクセル列の検出器面への投影曲線、ボクセル列及びボクセル列のセンタリング面への投影曲線を示す図。
【図10】同実施の形態のX線CT装置の再構成処理部におけるボクセル列及びセンタリング面を示す図。
【図11】同実施の形態のX線CT装置の再構成処理部における検出器列のボクセルへの投影曲線、検出器列及び検出器列のセンタリング面への投影曲線を示す図。
【図12】同実施の形態のX線CT装置の再構成処理部におけるセンタリング面上での補間処理などのデータのリサンプリングを説明するための図。
【図13】同実施の形態のX線CT装置の再構成処理部における3次元再構成処理の流れを示す図。
【図14】同実施の形態のX線CT装置で使用されるX線検出器の他の例としての平面型X線検出器を示す図。
【図15】同実施の形態の平面型X線検出器を使用したX線CT装置におけるボクセル列と平面型X線検出器との位置関係を示す図。
【図16】同実施の形態の平面型X線検出器を使用したX線CT装置におけるボクセル列の平面型X線検出器の検出器面への投影曲線を示す図。
【図17】同実施の形態の平面型X線検出器を使用したX線CT装置におけるボクセル列のセンタリング面への投影曲線を示す図。
【図18】検出器列をボクセルへの逆投影における1回補間の例を説明するための図。
【図19】検出器列をボクセルへの逆投影における2回補間の例を説明するための図。
【図20】この発明の第2の実施の形態のX線CT装置の再構成処理部における検出器列をボクセルへの逆投影における高精度な2回補間の第1の例を説明するための図。
【図21】同実施の形態のX線CT装置の再構成処理部における検出器列をボクセルへの逆投影における高精度な2回補間の第2の例を説明するための図。
【図22】同実施の形態のX線CT装置の再構成処理部における5列の検出器列のセンタリング面への投影曲線、センタリング列の数を5列としたときのボクセル列のセンタリング面への投影直線と重み付けされるセンタリング列及びセンタリング列の数を10列としたときのボクセル列のセンタリング面への投影直線と重み付けされるセンタリング列を示す図。
【図23】従来例のX線CT装置におけるファンビーム構成及びFOVを示す図。
【図24】同従来例のX線CT装置におけるピクセルを説明するための図。
【図25】従来例のX線CT装置におけるファン−パラ変換法を説明するための図。
【図26】従来例のX線CT装置におけるセンタリング軸を使用したファンビーム再構成法を説明するための図。
【図27】従来例のX線CT装置におけるコーンビームを示す図。
【図28】従来例のX線CT装置におけるコーンビームに対するボクセルを説明するための図。
【図29】従来例のファンビームを使用したX線CT装置における検出器データのピクセルへの逆投影を説明するための図。
【図30】従来例のコーンビームを使用したX線CT装置における検出器データのボクセルへの逆投影を説明するための図。
【符号の説明】
3…X線源、
5…X線検出器、
12…再構成処理部、
21…再構成処理制御部、
22…コンボリューション演算部、
23…第1のデータメモリ、
24…第1の逆投影部、
25…第2のデータメモリ、
26…第2の逆投影部、
27…画像メモリ。[0001]
BACKGROUND OF THE INVENTION
The present invention irradiates a target object with X-rays radiated in a conical shape from an X-ray source, detects X-rays transmitted through the target object with an X-ray detector, and obtains the detection obtained from the X-ray detector. The present invention relates to an image reconstruction processing device and an X-ray CT device for reconstructing a transmission image of the object based on data.
[0002]
[Prior art]
As shown in FIG. 23, a conventional X-ray CT (computed tomography) apparatus uses a fan beam in which an X-ray beam is emitted from a X-ray source 101 in a fan shape (fan shape).
Such an X-ray CT apparatus irradiates a subject with an X-ray beam emitted from the X-ray source 101, and an X-ray detector 102 in which the X-rays passing through the subject are arranged in a row in a row with about 1000 channels. Data is collected by detection, and the X-ray source 101 and the X-ray detector 102 are rotated around the subject while collecting data about 1000 times during one rotation (one data collection is referred to as one view). Then, an X-ray transmission image of the subject is reconstructed based on the collected data. Note that the FOV 103 indicates an effective visual field.
An image reconstruction formula when this fan beam is used is calculated by (Formula 1).
[0003]
[Expression 1]
Figure 0003825492
[0004]
As can be seen from this (Equation 1), in the reconstruction of the fan beam, it is necessary to multiply the data obtained from the X-ray detector by weighting depending on the position of the pixel to be reconstructed and backproject it. Complicated processing.
That is, as shown in FIG. 24, pixels constituting an image to be reconstructed are set with respect to the effective field of view FOV 103, and data obtained in each channel of the X-ray detector 102 is used as a focus (X The radiation point of the X-ray beam of the radiation source 101) -pixel distance FpixelD (X) is multiplied by the reciprocal of the square and backprojected to the corresponding pixel. Note that FpixelD is Focus-Pixel-Distance.
There is also a method of direct backprojection, but in this case, polar coordinate conversion is required, resulting in a complicated calculation.
[0005]
Therefore, at present, two types of image reconstruction methods have been devised for X-ray CT apparatuses using a fan beam.
One method is called a fan-para conversion method, which is a parallel beam projection data obtained by rearranging and interpolating projection data obtained from a fan beam X-ray detector as shown in FIG. Is generated (including conversion), and data obtained thereby is back-projected so as to be performed by an X-ray CT apparatus using a conventional parallel beam.
While calculation of data conversion and interpolation processing are required, at the time of back projection, different weighting processing is not required for each reconstructed pixel at the time of fan beam, and one data (parallel beam projection data) is transferred to the beam path ( Processing can be simplified because it is sufficient to perform back projection on all the pixels on the straight line connecting the radiation point and the channel of the X-ray detector when parallel beams are formed.
[0006]
Another method is a fan beam reconstruction method using a centering axis, the details of which are disclosed in Japanese Patent Laid-Open No. 55-99240.
As shown in FIG. 26, once the data obtained from the X-ray detector is projected (backprojected) onto a predetermined reference axis (centering axis) 104 parallel to the pixel row, it is again reconstructed for each pixel row to be reconstructed. Back projection. By performing back projection once on the centering axis in this way, different weighting processing for each pixel becomes the same for each pixel column, so that high-speed and simple processing is possible.
[0007]
The steps of the fan beam reconstruction method using this centering axis will be described below.
1. The raw data Data-Raw is obtained by multiplying various corrections such as X-ray intensity correction of the projection data Data-Proj by the cos term. Data-Conv is obtained by performing a convolution operation between the raw data Data-Raw and the reconstruction function.
2. Data-Conv is projected to a point on a certain reference axis (centering axis, for example, X and Y axes serving as a reference on which pixels are arranged) by weighting processing of (Equation 2) to obtain Data-Center. In this equation (Equation 2), FcpD (X) is the distance between the focal point and the centering axis point.
[Expression 2]
Figure 0003825492
[0008]
The process of 3.2 is repeated to project the fan beam projection data of all views onto the corresponding reference axis.
4). Data-Back is obtained by weighting the data Data-Center of a certain view projected on the reference axis with the same weighting (Equation 3) for every pixel column for all the pixel columns of the image to be reconstructed. The obtained projection data Data-Back is back-projected (added to an address corresponding to a pixel in the image memory). A and B are shown in FIG.
[Equation 3]
Figure 0003825492
[0009]
The process of 5.4 is repeated for all views, and projection data for all views is backprojected.
[0010]
Here, when (Equation 3) is transformed with reference to FIG.
[Expression 4]
Figure 0003825492
[0011]
This agrees with the content of the integral of the fan beam reconstruction formula (Equation 1).
[0012]
By projecting once onto the centering axis and then projecting back to each pixel column, the weight at the time of back projection (Equation 1), which was originally different for each pixel, is set to the same weight (Equation 2) for each pixel column. Simplification of generation calculation itself and reduction of the number of calculations have been achieved.
Although not described in detail here, a complicated correspondence between detector elements arranged at equal angles on an arc and pixels arranged at equal pitch on a straight line is also simplified.
A conventional X-ray CT apparatus using a fan beam realizes high-speed image reconstruction processing by one of the two methods described above.
[0013]
On the other hand, as shown in FIG. 27, on a cylindrical surface, such as a cone beam in which an X-ray beam is emitted from an X-ray source 201 in a conical shape, and N rows of fan beam detector rows are stacked in the Z-axis direction. An X-ray CT apparatus has been devised that captures X-ray fluoroscopic images using a two-dimensional X-ray detector 202 in which detector elements (M channels × N columns) are arranged.
A typical cone beam reconstruction (Feldkamp reconstruction) in such an X-ray CT apparatus using a cone beam is disclosed in the following document.
[0014]
"Practical cone-beam algorithm"
LAFeldkamp, LCDavis, and JWKress
J. Opt. Soc. Am. A / Vol.1, No.6, pp.612-619 / June 1984
This is an approximation obtained by extending the fan beam (in a two-dimensional plane) reconstruction algorithm [Filtered-Backprojection], which is a mathematically exact reconstruction method, in the Z-axis direction. This is a typical three-dimensional reconstruction algorithm.
[0015]
This cone beam reconstruction method is intended for conventional scanning using cone beams, and includes the following steps. In this cone beam, voxels as three-dimensional pixels are used as shown in FIG. 28 instead of pixels as two-dimensional pixels.
[0016]
1. Projection data weighting
The projection data is multiplied by a term depending on the Z coordinate and a cos term.
[0017]
2. Convolution calculation
A convolution operation between the data obtained by the process 1 and the same reconstruction function as the fan beam is performed.
[0018]
3. BackProjection
Data obtained by the process 2 is projected back onto the path (from the focal point to the detector channel) through which the X-ray passes. That is, the point where the straight line passing through the voxel back-projected from the focal point intersects the detector plane is calculated, and the data to be back-projected is created by interpolation etc. from the data of the two processes around that point, and it is generated by FvoxelD (X ) Is weighted by the reciprocal of the square of) and backprojected. This backprojection is performed over 360 ° (one rotation).
[0019]
When expressed by a formula similar to the fan beam reconstruction formula, it is as follows.
Note that FvoxelD (X) = Focus-Voxel-Distance is the focus-voxel distance.
[Equation 5]
Figure 0003825492
[0020]
The three-dimensional reconstruction formula (cone beam reconstruction formula) (formula 4) is very similar to the fan beam reconstruction on the formula, but the data-back backprojection method is greatly different.
In the two-dimensional fan beam reconstruction, as shown in FIG. 29, the projection data is back-projected from the data of the detectors arranged in one dimension for all the pixels (pixels) in the reconstruction plane. In the reconstruction of the beam (Feldkamp), as shown in FIG. 30, a point where the straight line connecting the focal point and the voxel to be reconstructed intersects the two-dimensional X-ray detector surface is obtained and involved in the intersection. Data obtained from the detector elements is backprojected to all voxels located on the line.
[0021]
Therefore, when a certain surface is reconstructed by cone beam reconstruction, such as fan beam reconstruction, data of a specific detector array and channel is backprojected only to a part of the voxels of the reconstruction surface. Since it is necessary to select data to be back-projected for each voxel (detector row and detector channel), the three-dimensional positional relationship between the straight line connecting the reconstructed voxel and the focus and the X-ray detector plane is Become important.
In addition, when considering a detector array with the same Z coordinate and a straight line connecting the detector element and the focal point, the voxel through which the straight line passes on a certain plane (reconstructed plane) is detected with the focus at the center. Since this is arranged on a similar figure on the instrument surface (concentric circle in the case of a cylindrical detector), the calculation of this positional relationship becomes very complicated.
[0022]
[Problems to be solved by the invention]
As described above, when a cone beam and a two-dimensional X-ray detector are used in a conventional X-ray CT apparatus, the calculation of the image reconstruction is complicated and enormous, and is generally popular. However, there is a problem that the processing time is too long for a computer or the like that can be realized.
Accordingly, an object of the present invention is to provide an image reconstruction processing apparatus and an X-ray CT apparatus capable of realizing accurate reconstruction of an image photographed using a cone beam.
[0023]
[Means for Solving the Problems]
The first aspect of the present invention irradiates an object with X-rays radiated in a conical shape from an X-ray source, and transmits X-rays transmitted through the object. Two or more detection elements are arranged in two dimensions In an image reconstruction processing apparatus for reconstructing an image of the object based on detection data detected by an X-ray detector and obtained from the X-ray detector, Means for replacing detection data detected by the plurality of detection elements with data on a centering surface set in a voxel assembly substantially parallel to a specific surface of the voxel; and the data replaced with the centering surface, Means for backprojecting each voxel of the body It is characterized by that.
A second aspect of the present invention includes an X-ray source that generates X-rays, an X-ray detector in which a plurality of detection elements that detect X-rays transmitted through the subject are arranged in the body axis direction of the subject; In the X-ray CT apparatus having means for obtaining backprojected data based on an output signal from the X-ray detector and image reconstruction processing means for reconstructing the image of the subject, the image reconstruction processing means Is based on a three-dimensional positional relationship between a voxel that is an image spatially arranged in three dimensions and the X-ray source. And calculating a projection curve obtained by projecting the voxel row onto the X-ray detector, and specifying the data to be back-projected. Means for backprojecting the specified backprojected data onto the voxels.
[0025]
DETAILED DESCRIPTION OF THE INVENTION
A first embodiment of the present invention will be described below with reference to FIGS.
FIG. 1 is a configuration diagram of an X-ray CT apparatus according to the first embodiment. FIG. 2 is an external view of the gantry of FIG.
A gantry (also referred to as a gantry) 1 as a projection data measuring system is a two-dimensional array in which an X-ray source 3 that generates a cone-beam X-ray bundle approximate to a cone and a plurality of detection elements are arranged in a two-dimensional manner. A type X-ray detector 5 is accommodated. The X-ray source 3 and the X-ray detector 5 are mounted on the rotating ring 2 in a state where they face each other with a subject placed on a slide top plate of a bed 6.
As the X-ray detector 5, a plurality (1000 channels) of detection elements are arranged in such a manner that a plurality of rows (10 rows) of one-dimensional array detectors arranged in a one-dimensional manner for fan beams are stacked. (See FIG. 27) and mounted on the rotating ring 2. Here, one detection element is defined to correspond to one channel.
[0026]
X-rays from the X-ray source 3 are exposed to the subject via the X-ray filter 4. X-rays that have passed through the subject are detected as electrical signals by the X-ray detector 5.
The X-ray controller 8 supplies a trigger signal to the high voltage generator 7. The high voltage generator 7 applies a high voltage to the X-ray source 3 at the timing of receiving the trigger signal. As a result, X-rays are emitted from the X-ray source 3.
The gantry bed controller 9 controls the rotation of the rotating ring 2 of the gantry 1 and the slide of the slide top plate of the bed 6 in synchronization. The system controller 10 as a control center of the entire system is configured so that the X-ray controller 8 and the pedestal bed controller are configured so as to execute a so-called helical scan in which the X-ray source 3 moves along a spiral trajectory as viewed from the subject. 9 is controlled.
Specifically, the rotating ring 2 continuously rotates at a constant angular velocity, the slide top plate of the bed 6 moves at a constant speed, and continuously or intermittently from the X-ray source 3 at every constant angle. The line is exposed.
[0027]
The output signal from the X-ray detector 5 is amplified by the data acquisition unit 11 for each channel and converted into a digital signal. Projection data output from the data collection unit 11 is taken into the reconstruction processing unit 12.
The reconstruction processing unit 12 obtains back projection data reflecting the X-ray absorption rate for each voxel based on the projection data.
In a helical scan type X-ray CT apparatus using a cone beam, an effective field of view (FOV, imaging region) has a cylindrical shape centering on the rotational center axis of the helical scan. Voxels (pixels arranged three-dimensionally) are defined (see FIG. 28), and back projection data of each voxel is obtained from the projection data from the X-ray detector 5. The three-dimensional image data or tomographic image data created based on the backprojection data is sent to the display device 14 and visually displayed as a three-dimensional image or tomographic image.
[0028]
As shown in FIGS. 3A and 3B, the geometry of this X-ray CT apparatus is
Number of detector rows M = 20,
Height of each row in the Z-axis direction Dseg = 2mm,
X-ray detector thickness M × Dseg = 40 mm,
Number of channels N = 1000,
Focus-rotation center distance FCD (Focus-center-Distance) = 600mm,
Focus-detector distance FDD (Focus-Detector-Distance) = 1200 mm,
Effective field diameter FOV (Field of View) = 500mm,
Effective viewing angle (fan angle) θ = 50 °
It has become.
[0029]
FIG. 4 is a block diagram showing a main configuration of the reconfiguration processing unit 12.
21 is a reconstruction process that constitutes the control unit main body of the reconstruction processing unit 12 that controls the calculation, such as data selection, weighting, centering process, backprojection, and the entire three-dimensional reconstruction process in the convolution process and the backprojection process. It is a control unit.
The convolution operation unit 22 performs a convolution process on the projection data collected by the data collection unit 11, and the convolution data obtained by the convolution process is stored in the first data memory 23.
[0030]
The first back projection unit 24 performs back projection (projection) processing on the convolution data stored in the first data memory 23 onto a preset centering surface, and the back projected centering data Stored in the data memory 25.
[0031]
The second backprojection unit 26 performs backprojection (three-dimensional backprojection) processing on the centering data stored in the second data memory 23-2 on the voxel, and the backprojected reconstruction data (image) is an image. Stored in the memory 27.
[0032]
Here, regarding the reconstruction, several points will be considered regarding the three-dimensional reconstruction by the cone beam.
With coordinate transformations such as reconstructed voxel row (straight line), centering surface (plane), detector row (arc), detector surface (cylindrical, but flat when considered as expanded on data memory) The relationship will be described.
First, consider the reconstructed voxel array, detector plane, and centering plane.
As shown in FIG. 5, it will be explained that the projection of one voxel row (straight line) onto the detector surface becomes a curve as shown in FIG.
[0033]
FIG. 7A shows a view of the reconstructed voxel array, detector surface, and centering surface observed from the Z-axis direction (upward). The channels (detector rows) of the X-ray detector 5 are arranged on an arc having an equal angle in the Z-axis direction from the focus F, and the rotation angle by the focus rotation (view) is φ and the angle in the channel direction as shown in the figure. Is represented by θ.
At this focal position, the centering surface is on the X-axis (included in the X-axis Z-axis plane), and the X-axis coincides with the Xcp-axis that is the X-axis of the centering surface. The coordinates of the reconstructed voxel sequence are defined as (Xv, Yv)... (Xv varies with the voxel, but Yv is constant).
As shown in FIG. 7 (b), when the length of the line segment FC (focus-rotation center distance) is FCD, the length of the line segment FV0 (focus-voxel distance) is FCD ', and a certain voxel V is considered. When the line segment FV is projected onto the XY plane, the length (in the dotted line in the figure) is FCD ", and the focus-detector distance is FDD. Also, a straight line passing through the focus and voxel is drawn on the XY plane. When the distance from the focal point is FCD, the Z coordinate is Z0, the Z coordinate of the voxel V is Zv, and the coordinate on the developed detector surface is the horizontal axis (Xdet = θ · FDD) and the vertical axis (Zdet).
[0034]
First, consider projecting the reconstructed voxel train onto the detector plane. Xcp represents the X coordinate on the centering plane of the line segment FV, and is obtained by the following (Equation 5).
[Formula 6]
Figure 0003825492
[0035]
Therefore, ΔFCXcp and ΔFVo Xv are similar, so using the following relational expression (Equation 6),
[Expression 7]
Figure 0003825492
[0036]
It becomes. This (Expression 7) is an expression of the channel direction in the projection of Xv → θ conversion and reconstruction voxel array onto the detector surface.
[0037]
Now, since the relationship between the line segment FXcp and the line segment FCD is (Expression 8) shown below, (Expression 9) is obtained.
[Equation 8]
Figure 0003825492
[0038]
Therefore, the expression of the column direction in the projection of the reconstructed voxel array onto the detector surface is as follows (Expression 10).
[Equation 9]
Figure 0003825492
[0039]
From (Equation 7), it can be seen that a reconstructed voxel array arranged at equal pitches is projected onto the detector surface to form a non-linear array. Accordingly, the left side in the figure, for example, the interval when the first voxel and the second voxel are projected on the detector surface (for example, 2.4 channels), and the right side in the figure, for example, the interval between the 511 voxel and the 512th voxel ( Unlike, for example, 3.5 channels), the arrangement becomes nonlinear depending on the angle θ.
Further, it can be seen from (Equation 10) that when a linearly reconstructed voxel array whose Z coordinate is fixed at Zv is projected onto the detector surface, a nonlinear curve as shown in FIG. 6 is formed depending on the angle θ. However, the Z coordinates Zv and Zdet of the reconstructed voxel row are in a proportional relationship.
[0040]
FIG. 8 and FIG. 9 show this, and FIG. 8 shows the definition of the voxel row V, the centering surface C, the detector surface D, the variables, and the end points and the center points. FIG. 10 shows the voxel row V and the centering surface C. The centering surface C is arranged in parallel to the voxel row V.
The projection of the voxel row, which is a straight line shown in FIG. 9B, onto the centering plane is also a straight line as shown in FIG. 9C. The pitch of the voxels and the pitch of the points on the centering surface is a fixed ratio, and no distortion occurs. That is, if the data to be backprojected is a straight line and an equal pitch on the centering surface, the backprojection onto the voxel is naturally a straight line and an equal pitch. That is, the relationship between the voxel and the centering surface is a simple enlargement / reduction relationship.
[0041]
However, when a voxel row is projected onto the detector surface, nonlinear distortion occurs in both the channel direction and the row direction, as shown in FIG. However, as shown by the projection of two voxel rows having different Z coordinates Zv, the above-described scaling relationship is established between the projected images of the two voxel rows on the detector surface. ing.
[0042]
Based on the above considerations, the simplest cone beam reconstruction method is the step described below.
1. The projection data from the X-ray detector is subjected to correction processing such as X-ray intensity correction, and then subjected to Feldkamp weighting processing and stored in a data memory.
2. The corrected projection data stored in the data memory is convolved with the reconstruction function and stored in the data memory.
3. According to (Equation 7) and (Equation 10), a projection curve obtained by projecting the back-projected voxel train on the detector surface is calculated, and the data to be back-projected is selected and its address is generated.
4). The corresponding data is read out, and after a predetermined weighting process, is added to the position of the voxel corresponding to the image memory.
[0043]
This is referred to as “direct backprojection”. In the process 3, a projection curve may be calculated in advance and stored in a data memory or the like as a table. This is called “table method”.
In either method, the projection curve may be an approximate curve. However, in this method, the calculation of data selection in the column direction is added in addition to the data selection calculation for the back projection (channel direction) and the weight calculation that existed even in the two-dimensional back projection of the fan beam projection data. It becomes quantity. In the case of the “table method”, a memory for storing a huge table is required.
[0044]
Therefore, in the first embodiment of the present invention, a method of back-projecting data (projection data) obtained from the X-ray detector 5 once onto the centering plane and then back-projecting it onto each voxel is performed.
Contrary to the consideration in FIG. 9, the case where the pitch and straight line on the detector surface (or data memory) are projected onto the centering surface will be considered. (Equation 11) obtained by solving (Equation 5) and (Modification of Equation 10) by (Equation 8) is shown below.
[Expression 10]
Figure 0003825492
[0045]
Further, when a straight line on the centering surface is projected onto the detector surface, the following (formula 7) and (formula 11) are obtained.
[Expression 11]
Figure 0003825492
[0046]
Therefore, nonlinear distortion according to (Expression 5) and (Expression 11) is generated instead of the nonlinear distortion according to (Expression 7) and (Expression 10).
[0047]
This is shown in FIG. When the data of all the channels of the detector array indicated by the straight line in FIG. 11 (b) is projected onto the centering plane, as shown in FIG. 11 (c), Xcp, Zcp are obtained by the above-described nonlinearity (Equation 4 and Equation 8). Both directions are distorted.
Both the distortion in the Z direction and the distortion of the array with the same pitch in the X direction are opposite to those in FIG. As described above, there is no distortion between the centering plane and the reconstructed voxel array. However, as shown in FIG. 11A, when the data of all the channels of the detector array are projected onto the voxels, the same distortion occurs. That is, the Z-direction distortion of the projection onto the centering plane shown in FIG. 11 (c) and the distortion of the X-direction equal pitch array and the distortion of the projection onto the voxel shown in FIG. The distortion of the array with the equal pitch in the direction is the same and has an enlargement / reduction relationship.
[0048]
Therefore, data resampling is performed on the centering surface by interpolation processing or the like, and processing is performed so that the data is arranged at a constant pitch on a straight line. The result is FIG. 12 (a).
The relationship between the coordinate system Xcp, Zcp on the centering plane and the coordinate system Xv, Zv of the reconstructed voxel is expressed by applying (Equation 5), (Equation 6), (Equation 9), (Equation 10) as follows: Equations 12) and 13 are obtained.
[Expression 12]
Figure 0003825492
[0049]
Therefore, when reconstructing an axial section at the slice position Z = Zcp in the reconstructed voxel, if the axial section is considered as a circular FOV inscribed in a square as shown in FIG. In 13), Zv is always a constant in the cross section, Yv is a constant in the voxel row, and Xv is changed in units of voxels. Therefore, the projection of a square (shown by a broken line) and its inscribed circular FOV is a deformation (not shown) of a trapezoid (shown by a broken line) and its inscribed circle.
Further, the straight line obtained by projecting the voxel row onto the centering plane translates up and down in the centering plane corresponding to the position of the voxel row. Therefore, the Z coordinate Zcp corresponding to the voxel string to be backprojected is obtained by (Equation 13), the X coordinate Xcp corresponding to the voxel to be backprojected is obtained by (Equation 12), and the corresponding data is backprojected to the target voxel. . This is repeated for all voxels for all views to perform back projection.
[0050]
It is as follows when expressed by an expression. FdpD is a focus-detector element distance, and FcpD (X, Z) is a focus-centering point distance.
[Formula 13]
Figure 0003825492
[0051]
Here, (Expression 14) is an expression indicating the data subjected to the convolution process, (Expression 15) is an expression indicating the data subjected to the centering process, and (Expression 16) is the data back-projected to the voxel. It is a formula.
[0052]
This (Equation 16)
[Expression 14]
Figure 0003825492
[0053]
This (Equation 17) agrees with the three-dimensional reconstruction equation (Equation 3).
[0054]
That is, the steps of the three-dimensional reconstruction method (cone beam reconstruction method) using the centering surface of the present invention will be described below.
1. The reconstruction processing control unit 21 obtains raw data Data-Raw by performing correction processing such as X-ray intensity correction on the projection data Data-Proj from the data collection unit 11 and stores the raw data Data-Raw in the first data memory 23.
2. The reconstruction processing control unit 21 reads the data in the first data memory 23 by the convolution operation unit 22, and after the Feldkamp weighting (the first term of (Equation 14)) processing, the reconstruction function Conv-Function and the convolution (The second term of (Equation 14)) and stored in the first data memory 23.
[0055]
3. Based on the convolution data stored in the first data memory 23 by the first back projection unit 24, the reconstruction processing control unit 21 performs either of the following (i) or (ii) (Equation 15 ) Is performed.
(i) According to (Equation 5) and (Equation 11), a projection curve obtained by projecting the detector surface data onto the centering surface is calculated, and the convolution data is weighted and then projected onto the centering surface to calculate the centering data. Further, the centering data on the centering surface is resampled in a lattice shape as shown in FIG. 12A and stored in the second data memory 25.
(ii) According to (Equation 7 variant) and (Equation 11 variant), the position of the projection point obtained by projecting the grid-like data position resampled in (i) on the centering plane onto the detector plane is calculated. After the weighted addition of the convolution data of 4 detector elements (2 rows x 2 CH) around the projection point and weighting in the three-dimensional backprojection of (Equation 15), it is projected onto the centering surface and centered data Data-Center Is calculated and stored in the second data memory 25.
[0056]
4). The reconstruction processing control unit 21 calculates a projection point (straight line) obtained by projecting the voxel (column) to be reconstructed onto the centering plane according to (Equation 12) and (Equation 13), and, for example, four pieces of backprojected data ( 2 columns x 2CH) Select and generate its address. 5). The reconstruction processing control unit 21 reads out the corresponding centering data Data-Center from the second data memory 25 by the second back projection unit 26, and adds the weighting when the number of data is plural, and reconstructs the two-dimensional fan beam. Similar to the back projection of the centering data in the configuration, after the weighting process of the square of A / B (this validity has been proved by (Equation 17)), it is added to the position of the corresponding voxel in the image memory 27. .
As described above, three-dimensional reconstruction is possible.
The flow of this three-dimensional reconstruction process is shown in FIG.
[0057]
In the calculation of back projection (projection) data from a point on the detector surface or a point on the centering surface to a point or voxel on the centering surface in the processes 3 and 5 described above, four-point Bi-Linear interpolation is used. Re-sampling may be performed using non-linear interpolation such as linear interpolation, Spline interpolation, etc., or other interpolation. Further, back projection data may be calculated by selecting a point closest to the corresponding projection point, for example, as Nearest Neighbor without performing interpolation.
[0058]
In the four-point Bi-Linear interpolation, for example, four points Data (j, n), Data (j, n + 1), Data of j column and j + 1 column, n channel and n + 1 channel with respect to one projection point. (j + 1, n) and Data (j + 1, n + 1) are subject to calculation as data to be interpolated.
At this time, the following calculation is performed in the four-point Bi-Linear interpolation.
Data (j, n) × wch + Data (j, n + 1) × (1-wch) ... CH (j)
Data (j + 1, n) × wch + Data (j + 1, n + 1) × (1-wch) ... CH (j + 1)
CH (j) * wro + CH (j + 1) * (1-wro) ... SEG (j, j + 1)
The above processing is repeated for a plurality of projection points in order. That is, CH (j) is calculated, then CH (j + 1) is calculated, and SEG (j, j + 1) is calculated from those calculation results. Then, CH (j) is calculated for the next projection point, then CH (j + 1) is calculated, and SEG (j, j + 1) is calculated from the calculation results. In the same manner, interpolation processing is repeatedly performed for each projection point. This repeated interpolation process takes time.
Therefore, CH (j) and CH (j + 1) are processed in parallel, and the processing time is shortened by pipeline processing of SEG (j, j + 1) so as to connect the parallel processing. Interpolation processing can be performed in substantially the same processing time as two-point interpolation.
[0059]
If the number of centering columns is greater than or equal to the number of voxel columns, the nearest projection point is selected as the Nearest Neighbor and the backprojection data is calculated to reduce image degradation and the number of calculations to be 4 points Bi-Linear. Compared to interpolation, the effect is about 1/3 and the processing time is shortened.
[0060]
Although omitted here, resampling processing on the centering surface in step 3 is performed only for the purpose of simplifying non-linearity in the channel direction and weight generation at the time of back projection as in the case of single slice CT (fan beam CT). May be omitted, and a projection curve corresponding to the centering data corresponding to the voxel to be back-projected in the process of 4 may be generated. As a price for omitting the resampling process, the projection curve at the time of back projection becomes complicated, but there is an advantage that the total number of interpolations is reduced.
[0061]
As described above, according to the first embodiment, in the three-dimensional backprojection, once the projection is performed on the centering plane, the detector CH array arranged at an equal angle on the arc and the same pitch on the straight line are arranged. This simplifies the complicated correspondence between the voxel strings.
That is, once the projection is performed on the centering plane, different weights for each voxel unit at the time of backprojection can be made the same for each voxel column unit, thereby reducing the number of weight calculations and simplifying the calculation itself.
The Z-direction distortion (non-linearity) at the time of back projection peculiar to cone beam CT was corrected by arranging the data in a grid pattern on the centering plane, and back projection onto the reconstructed voxel could be simplified. Such correction was realized by projection onto a grid centering point or resampling processing.
[0062]
As described above, the processing of the second backprojection unit 26 is simplified, and a simple configuration without significant expansion compared to the configuration of the existing two-dimensional backprojection (fan beam reconstruction) backprojection unit. Thus, image reconstruction in three-dimensional backprojection can be realized.
[0063]
In the first embodiment, a cone beam X-ray CT apparatus having a cylindrical X-ray detector has been described. However, the present invention is not limited to this, and for example, as shown in FIG. A cone beam CT apparatus having a flat planar X-ray detector can provide the same effect.
That is, as shown in FIG. 15, with respect to the voxel row indicated by the straight and broken arrows in the effective field of view FOV, in the projection onto the planar X-ray detector surface, as shown in FIG. In addition, the broken line arrow does not cause distortion (including non-linear distortion) in the channel direction (lateral direction). Furthermore, although distortion does not occur in the column direction (vertical direction), the projection line becomes an inclined straight line, which is reduced compared to the case of a cylindrical X-ray detector, but at the time of data selection and back projection. Weight generation is complex.
[0064]
Therefore, by performing a centering process that performs back projection (projection) on the centering surface in an intermediary manner, as shown in FIG. 17, the voxel row is projected in parallel on the centering surface, so that the straight arrows in the effective field of view FOV and By projecting the voxel sequence indicated by the dashed arrow onto the plane X-ray detector surface, the centering sequence used for backprojection data calculation in units of voxel sequence can be fixed, thus simplifying both data selection and weight generation. Can do.
As in the case of the cylindrical X-ray detector, the weight at the time of back projection can be reduced from the calculation for each voxel column to the calculation for each voxel column. Further, if the number of centering rows is increased, the effect of improving the interpolation accuracy is the same as that of the cylindrical X-ray detector.
[0065]
A second embodiment of the present invention will be described with reference to FIGS. The basic configuration of the X-ray CT apparatus (cone beam CT apparatus) in the second embodiment is almost the same as that of the first embodiment described above. The processing in the image reconstruction processing unit will be described.
[0066]
When interpolation is repeated, blurring generated in the image is amplified. However, in order to simplify the projection curve and the weight, when back projection is performed once on the centering surface and then back projection onto the voxel, two interpolations are inevitable.
That is, in the case of one-time interpolation as shown in FIG. 18, the voxel data DB is D1 and D2 with respect to the (original) data positions (D1, D2, D3) obtained by back-projecting the detector row onto the voxels. Is in the position of dividing 8: 1
DB = (1/9) × D1 + (8/9) × D2
Thus, it is determined by the data positions D1 and D2, and interference of other data such as D3 is completely eliminated.
On the other hand, as shown in FIG. 19, in the case of two-time interpolation in which back projection is once performed on the centering plane, first, centering points Dc1 and Dc2 are backprojected (projected) to the centering plane (original) data. With respect to the position (D1, D2, D3), the position between D1 and D2 is divided into 5: 4 and the position between D2 and D3 is divided into 2: 1.
Dc1 = (4/9) × D1 + (5/9) × D2
Dc2 = (1/3) × D2 + (2/3) × D3
In addition, when the voxel data (point) DB is in a position that divides Dc1 and Dc2 into 3: 5,
Figure 0003825492
It becomes. That is, although the voxel data DB is actually between D1 and D2, the interference term of D3 is added by two interpolations, and the terms of D1 and D2 also include errors accordingly. The blur will be amplified in the actual image. Therefore, it is necessary to reduce blurring generated in an image by improving the accuracy of interpolation.
[0067]
There are the following two examples of methods for improving the interpolation accuracy.
FIG. 20 is a first example of high-precision two-time interpolation.
Set the centering points (Dc1, Dc4, Dc6) on the centering surface so that they coincide with the (original) data positions (D1, D2, D3) back-projected (projected) onto the centering surface; Interpolation points (Dc2, Dc3, Dc5) are set at equal pitches between these centering points. As a result, it is possible to reduce the blur of the interpolation twice.
For example, when the voxel data DB is in a position that divides Dc3 and Dc4 into 2: 1,
Figure 0003825492
Thus, the interference term of D3 is completely eliminated at the original data position. However, the projection curve of the detector array on the centering plane has non-linear distortion, and the Zcp coordinates where the original detector data exists between the adjacent channels, that is, the centering points adjacent in the Xcp direction, slightly differ. It is difficult to set interpolation points (centering points) at the original data positions in all centering rows.
[0068]
FIG. 21 is a second example of high-precision two-time interpolation.
The number of equal-pitch interpolation points (centering points) is made as large as possible from the balance with the calculation time, and the original data positions do not necessarily have to have interpolation points. However, by increasing the number of centering columns at least than the number of detector columns, the interpolation accuracy is high in most positions, and blurring due to interpolation occurs because the second interpolation position (voxel) is the original data position. This is only when it is between the two interpolation points for the first time across (for example, between Dc9 and Dc10).
[0069]
For example, when the voxel data DB is in a position that divides Dc8 and Dc9 into 1: 1,
Figure 0003825492
Thus, the interference term of D3 is completely eliminated at the original data position. In addition, when the voxel data DB is at a position that divides Dc9 and Dc10 into 1: 1 (a position corresponding to D2),
Figure 0003825492
Although the interference term of D3 (in this case, D1 also applies) is not eliminated, the weights (coefficients) of D1 and D3 are smaller than D2, so they are reduced compared to the conventional blur, and this As described above, the range in which the blur occurs is limited to the distance between Dc9 and Dc10.
[0070]
An example in which this is applied to the centering surface is shown in FIG. FIG. 22 (a) is a projection curve onto the centering plane of five detector rows, and shows the original data position. As described in the step of the three-dimensional reconstruction method (cone beam reconstruction method) using the centering surface described in the first embodiment, the original data is arranged in a grid pattern on the centering surface. Create centering data.
FIG. 22B shows a case where the number of grid rows (centering rows) is five, which is the same as the number of detector rows.
The second centering column C2 is used for interpolation with a light weight, and the third centering column C3 is used for interpolation with a heavy weight with respect to a projection straight line S on a centering surface of a certain voxel column.
It can be seen that the width of data used for this interpolation becomes wider in the Zcp direction and blurring occurs. This corresponds to the double interpolation example of FIG.
[0071]
FIG. 22 (c) shows an example in which the number of centering rows is larger than the number of detector rows. In FIG. 22 (c), there are 10 columns.
The data position to be back-projected onto a certain voxel row is indicated by a straight line S as above. The fourth centering column C4 is used for interpolation with a light weight, and the fifth centering column C5 is used for interpolation with a heavy weight.
It can be seen that the width of the data used for this interpolation becomes narrower in the Zcp direction and blur can be suppressed. This corresponds to the second example of the high-precision double interpolation shown in FIG.
In FIG. 22, the interpolation method is linear linear interpolation based on the inverse distance ratio. However, the interpolation accuracy at the point of the arrow shown in the example (when DB is located between Dc8 and Dc9) is very high.
[0072]
Further, the number of centering columns is significantly increased with respect to the number of detector columns (for example, 50 or 500 centering columns with respect to 5 detector columns), and the second interpolation is performed by zero-order interpolation, that is, nearest centering. You can also use Nearest Neighbor to select a column.
[0073]
As described above, according to the second embodiment, it is possible to minimize the data blur due to the interpolation, which is unavoidable when performing the back projection process using the centering surface, by increasing the interpolation accuracy.
[0074]
【The invention's effect】
As described above in detail, according to the present invention, it is possible to provide an image reconstruction processing apparatus and an X-ray CT apparatus capable of realizing accurate reconstruction of an image captured using a cone beam.
[Brief description of the drawings]
FIG. 1 is a configuration diagram showing an X-ray CT apparatus according to a first embodiment of the present invention.
FIG. 2 is an external view showing a gantry of the X-ray CT apparatus according to the embodiment.
FIG. 3 is a view for explaining the geometry of the X-ray CT apparatus according to the embodiment;
FIG. 4 is a block diagram showing a main configuration of a reconstruction processing unit of the X-ray CT apparatus according to the embodiment.
FIG. 5 is a view showing projection of a voxel array onto a detector surface in the reconstruction processing unit of the X-ray CT apparatus according to the embodiment;
FIG. 6 is a diagram showing a projection curve on a detector surface of a voxel array in the reconstruction processing unit of the X-ray CT apparatus according to the embodiment;
FIG. 7 is a diagram for explaining a projection curve on a detector surface of a voxel array in the reconstruction processing unit of the X-ray CT apparatus according to the embodiment;
FIG. 8 is a diagram showing a centering surface showing definitions of voxel rows, centering surfaces, detector surfaces, variables, and end points and center points in the reconstruction processing unit of the X-ray CT apparatus according to the embodiment;
FIG. 9 is a diagram showing a projection curve of a voxel array on a detector surface, a projection curve of a voxel array and a centering surface of the voxel array in the reconstruction processing unit of the X-ray CT apparatus according to the embodiment;
FIG. 10 is a view showing a voxel row and a centering surface in the reconstruction processing unit of the X-ray CT apparatus according to the embodiment.
FIG. 11 is a diagram showing a projection curve onto a voxel of a detector array, a detector array, and a projection curve onto a centering plane of the detector array in the reconstruction processing unit of the X-ray CT apparatus according to the embodiment;
FIG. 12 is a diagram for explaining data resampling such as interpolation processing on the centering plane in the reconstruction processing unit of the X-ray CT apparatus according to the embodiment;
FIG. 13 is a view showing a flow of three-dimensional reconstruction processing in a reconstruction processing unit of the X-ray CT apparatus according to the embodiment;
FIG. 14 is a view showing a planar X-ray detector as another example of the X-ray detector used in the X-ray CT apparatus according to the embodiment;
FIG. 15 is a diagram showing a positional relationship between a voxel array and a planar X-ray detector in the X-ray CT apparatus using the planar X-ray detector according to the embodiment;
FIG. 16 is a view showing a projection curve of a voxel array on the detector surface of the planar X-ray detector in the X-ray CT apparatus using the planar X-ray detector of the embodiment;
FIG. 17 is a diagram showing a projection curve on a centering surface of a voxel array in the X-ray CT apparatus using the planar X-ray detector according to the embodiment;
FIG. 18 is a diagram for explaining an example of one-time interpolation in back projection of detector rows onto voxels;
FIG. 19 is a diagram for explaining an example of twice interpolation in back projection of detector rows onto voxels;
FIG. 20 is a diagram for explaining a first example of high-accuracy two-time interpolation in back projection of a detector array in a reconstruction processing unit of an X-ray CT apparatus according to a second embodiment of the present invention onto voxels; Figure.
FIG. 21 is a diagram for explaining a second example of high-accuracy two-time interpolation in backprojecting the detector row to the voxel in the reconstruction processing unit of the X-ray CT apparatus according to the embodiment;
FIG. 22 shows the projection curves on the centering surfaces of the five detector rows and the centering surfaces of the voxel rows when the number of centering rows is five in the reconstruction processing unit of the X-ray CT apparatus of the embodiment. The centering row weighted with the projected straight line to the centering surface of the voxel row when the number of the centering row and the number of centering rows to be weighted with the projected straight line is 10.
FIG. 23 is a diagram showing a fan beam configuration and FOV in a conventional X-ray CT apparatus.
FIG. 24 is a view for explaining pixels in the conventional X-ray CT apparatus;
FIG. 25 is a diagram for explaining a fan-para conversion method in a conventional X-ray CT apparatus;
FIG. 26 is a diagram for explaining a fan beam reconstruction method using a centering axis in a conventional X-ray CT apparatus.
FIG. 27 is a diagram showing a cone beam in a conventional X-ray CT apparatus.
FIG. 28 is a diagram for explaining voxels with respect to a cone beam in a conventional X-ray CT apparatus.
FIG. 29 is a diagram for explaining back projection of detector data onto pixels in an X-ray CT apparatus using a conventional fan beam.
FIG. 30 is a diagram for explaining back projection of detector data onto voxels in an X-ray CT apparatus using a conventional cone beam;
[Explanation of symbols]
3 ... X-ray source,
5 ... X-ray detector,
12 ... Reconfiguration processing unit,
21 ... Reconfiguration processing control unit,
22: Convolution calculation unit,
23. First data memory,
24. First back projection unit,
25 ... second data memory,
26 ... second back projection unit,
27: Image memory.

Claims (4)

X線源から円錐状に放射されたX線を対象物に照射し、この対象物を透過したX線を、複数の検出素子を2次元状に配列してなるX線検出器により検出し、このX線検出器から得られた検出データに基づいて、前記対象物の画像を再構成する画像再構成処理装置において、
前記複数の検出素子で検出された検出データを、ボクセルの特定面と略平行にボクセル集合体内に設定されたセンタリング面上のデータに置き換える手段と、
前記センタリング面に置き換えたデータを、前記ボクセル集合体の各ボクセルに対して逆投影する手段とを有することを特徴とする画像再構成処理装置。
X-rays radiated in a conical shape from an X-ray source are irradiated onto an object, and X-rays transmitted through the object are detected by an X-ray detector in which a plurality of detection elements are arranged in a two-dimensional manner , In the image reconstruction processing device for reconstructing the image of the object based on the detection data obtained from the X-ray detector,
Means for replacing detection data detected by the plurality of detection elements with data on a centering surface set in a voxel assembly substantially parallel to a specific surface of the voxel;
An image reconstruction processing apparatus comprising: means for back-projecting the data replaced with the centering plane onto each voxel of the voxel aggregate .
前記X線検出器の列数より前記センタリング面のデータ列数を多くしたことを特徴とする請求項1記載の画像再構成処理装置。 The image reconstruction processing apparatus according to claim 1 , wherein the number of data columns on the centering surface is larger than the number of columns of the X-ray detector . X線を発生するX線源と、被検体を透過したX線を検出する複数の検出素子が前記被検体の体軸方向に複数配列されたX線検出器と、前記X線検出器からの出力信号に基づいて逆投影するデータを得る手段と、前記被検体の画像を再構成する画像再構成処理手段を有するX線CT装置において、
前記画像再構成処理手段は、
3次元的に空間配置された画像であるボクセルと、前記X線源との3次元的な位置関係に基づいて、当該ボクセル列を前記X線検出器に投影した投影曲線を計算し、前記逆投影するデータを特定する手段と、
前記特定した逆投影するデータを、前記ボクセルに逆投影する手段とを備えたことを特徴とするX線CT装置。
An X-ray source that generates X-rays, an X-ray detector in which a plurality of detection elements that detect X-rays transmitted through the subject are arranged in the body axis direction of the subject, and an X-ray detector In an X-ray CT apparatus having means for obtaining backprojected data based on an output signal and image reconstruction processing means for reconstructing the image of the subject,
The image reconstruction processing means includes
Based on the three-dimensional positional relationship between the voxel, which is a three-dimensional spatially arranged image, and the X-ray source, a projection curve in which the voxel row is projected onto the X-ray detector is calculated, and the inverse Means for identifying the data to be projected ;
An X-ray CT apparatus comprising means for back-projecting the specified back-projected data onto the voxel.
前記特定手段は、前記ボクセルと前記X線源との3次元的な位置関係に基づいて逆投影するデータを特定するテーブルを有し、前記位置関係と前記テーブルに基づいて前記逆投影するデータを特定することを特徴とする請求項3記載のX線CT装置。  The specifying unit includes a table for specifying data to be backprojected based on a three-dimensional positional relationship between the voxel and the X-ray source, and the backprojected data is determined based on the positional relationship and the table. The X-ray CT apparatus according to claim 3, wherein the X-ray CT apparatus is specified.
JP00101596A 1996-01-08 1996-01-08 Image reconstruction processing apparatus and X-ray CT apparatus Expired - Lifetime JP3825492B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP00101596A JP3825492B2 (en) 1996-01-08 1996-01-08 Image reconstruction processing apparatus and X-ray CT apparatus
US08/778,820 US5838756A (en) 1996-01-08 1997-01-03 Radiation computed tomography apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP00101596A JP3825492B2 (en) 1996-01-08 1996-01-08 Image reconstruction processing apparatus and X-ray CT apparatus

Publications (2)

Publication Number Publication Date
JPH09187449A JPH09187449A (en) 1997-07-22
JP3825492B2 true JP3825492B2 (en) 2006-09-27

Family

ID=11489758

Family Applications (1)

Application Number Title Priority Date Filing Date
JP00101596A Expired - Lifetime JP3825492B2 (en) 1996-01-08 1996-01-08 Image reconstruction processing apparatus and X-ray CT apparatus

Country Status (1)

Country Link
JP (1) JP3825492B2 (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1277175B1 (en) * 2000-04-14 2010-05-12 General Electric Company Reconstruction of computed tomographic images using interpolation between projection views
JP4612290B2 (en) * 2002-09-09 2011-01-12 株式会社東芝 Computer tomography apparatus and program
JP4642331B2 (en) * 2003-08-21 2011-03-02 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー X-ray CT apparatus and control method thereof
JP5179007B2 (en) * 2005-12-01 2013-04-10 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー X-ray CT apparatus and X-ray CT image reconstruction method thereof
DE112009005019B4 (en) * 2009-06-30 2022-02-03 Analogic Corp. Efficient quasi-exact 3D image reconstruction algorithm for CT scanners
JP6085824B2 (en) * 2013-08-09 2017-03-01 朝日レントゲン工業株式会社 X-ray imaging apparatus and image processing method
CN106934763B (en) * 2017-04-17 2023-08-22 北京灵起科技有限公司 Panoramic camera, automobile data recorder, image processing method and device

Also Published As

Publication number Publication date
JPH09187449A (en) 1997-07-22

Similar Documents

Publication Publication Date Title
US6292530B1 (en) Method and apparatus for reconstructing image data acquired by a tomosynthesis x-ray imaging system
US5598453A (en) Method for X-ray fluoroscopy or radiography, and X-ray apparatus
US6452996B1 (en) Methods and apparatus utilizing generalized helical interpolation algorithm
US6256370B1 (en) Method and apparatus for performing tomosynthesis
JP5221394B2 (en) How to reconstruct image functions from radon data
US5838756A (en) Radiation computed tomography apparatus
JP2914891B2 (en) X-ray computed tomography apparatus
JP4360817B2 (en) Radiation tomography equipment
EP0520778B1 (en) Tomographic image reconstruction using cross-plane rays
JPH08509408A (en) Image reconstruction from cone beam data
JP4342164B2 (en) Computed tomography equipment
JPH06214031A (en) Apparatus and method for display of image
CN102781327A (en) Phase contrast imaging
JP2001095794A (en) Apparatus and procedure for reconstruction of image at computer tomography using gradient slice
JP2008012319A (en) Method and system for reducing artifact in tomosynthesis/imaging/system
JP3980696B2 (en) Image reconstruction processing device
JP3913798B2 (en) Image reconstruction processing device
JP3825492B2 (en) Image reconstruction processing apparatus and X-ray CT apparatus
JP3917684B2 (en) Method and apparatus for creating a tomographic image of an object
JPH10201751A (en) Method and device for eliminating boundary error in conical beam image
CN1759811B (en) Method and apparatus for reconstruction of tilted cone beam data
JPH09192126A (en) Image reconstitution processor
JP3878788B2 (en) Method and apparatus for performing tomosynthesis
JP4582997B2 (en) High speed computed tomography method
WO1999038439A1 (en) X-ray image pickup device

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20040507

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20050906

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20051107

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20060630

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100707

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100707

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110707

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120707

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130707

Year of fee payment: 7

EXPY Cancellation because of completion of term