(57)【特許請求の範囲】
【請求項１】
被検体に対し多方向からＸ線を照射することによって所定の投影角度間隔における複数の投影データを収集する収集手段と、
前記被検体内部のＸ線吸収係数分布を表す画像を、前記投影データに基づいて再構成する再構成手段と、
前記再構成手段において行われる逆投影演算に際し、前記投影角度間隔の不均一性に応じたパワーの差異を補正するための補正係数を各々の逆投影データに掛けることにより前記再構成画像の劣化を抑制する抑制手段と、
を具備することを特徴とする画像再構成装置。
【請求項２】
前記補正係数は、基準となる投影角度間隔に対する所定の投影角度における投影角度間隔の比に基づいて算出されることを特徴とする請求項１に記載の画像再構成装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、被検体内部のＸ線の吸収係数に応じた分布画像を、異なる角度からの複数の投影データに基づいて再構成する画像再構成装置に関する。
【０００２】
【従来の技術】
被検体を透過した放射線の線量分布を複数方向で計測した投影データから、被検体のＸ線吸収係数の分布画像を再構成する方法としては、参考文献（「医用画像工学ハンドブック」日本医用画像工学会監修、日本医用画像工学会・医用画像工学ハンドブック編集委員会編集、篠原出版、Ｉ．基礎論理、第５章、画像再構成）に幾つか紹介されている。その中でも直接的（解析的）再構成法の１つであるフィルター補正逆投影法（Filtered Backprojection ）がよく知られている。
【０００３】
現在のＣＴ（コンピュータ断層撮影）装置において適用される画像再構成法は、ハードウェア化を容易にするために多少の近似が行われていたり、空間フィルターを補正フィルターに加味している等の相違点を有しているが、基本的にはこののような方法を採用している。
【０００４】
この画像再構成法は、投影データが等角度間隔、すなわち各々が等しいピッチで収集されることを前提としている。しかしながら、Ｘ線を用いて画像再構成を行う装置として最もメジャーなＣＴ装置においても、Ｘ線管球の回転速度あるいはＸ線の照射角度にはＸ線管球の重み等に起因する微小なブレが生じるため、現実には不等角度間隔で投影データが収集される場合が起こり得る。さらに、近年話題になっている三次元ＣＴを、回転アーム機構を有する循環器Ｘ線撮影装置によって実現する場合においては、このような影響は顕著となる。具体的に言えば、回転開始から加速している期間及び減速から回転終了までの期間においては、回転速度が動的に（ダイナミックに）変化する。ところがＸ線曝射間隔は常に等時間間隔である。このため投影データ収集の角度ピッチがこの間にダイナミックに変化してしまう。
【０００５】
また、たとえ回転速度がほぼ一定であっても、Ｘ線管球や検出器の重さによる回転速度への直接的な影響や、回転アームのたわみ等による間接的な影響によって、投影データ収集の角度ピッチが変化することも考えられる。
【０００６】
したがって、このような投影データ収集における角度ピッチの不均一性により、フィルター補正逆投影法による再構成処理により得られた再構成像が著しく劣化してしまうという問題点がある。
【０００７】
上記した三次元ＣＴに係る画像再構成法に参考文献（L.A.Feldkamp他著、"Practical cone-beam algorithm" 、J. Opt. Soc. Am. A, No.6, 612-619(1984) ）あるいは（大石  悟他著、「コーンビーム投影系を用いた三次元再構成」、Med.Imag.Tech.,Vol.8, No.2, 127-130(1990) ）に記載のFeldkampのフィルター補正逆投影法がある。この方法は、Feldkamp等によって提案され、三次元再構成法において最も有名で且つ実用的な方法の一つとされており、投影データに対してコンボリューションフィルターを掛けた後、コーンビームの広がり補正を行いながら三次元的な逆投影演算を行うものである。
【０００８】
しかしながらこのFeldkampのフィルター補正逆投影法においても上述した問題は解消されることなく、投影データ収集における不均一な角度ピッチが再構成像を劣化させるように働いてしまう。
【０００９】
【発明が解決しようとする課題】
本発明は上述した事情を考慮してなされたものであり、その目的とするところは、投影データ収集時における角度ピッチの不均一性を起因として発生する偽像（アーチファクト）による再構成画像の劣化を、処理時間を大幅に増やすことなく効果的に抑制し得る画像再構成装置を提供することにある。
【００１０】
【課題を解決するための手段】
上記課題を解決し目的を達成するために本発明は以下の手段を備えている。
（１）本発明の画像再構成装置は、被検体に対し多方向からＸ線を照射することによって所定の投影角度間隔における複数の投影データを収集する収集手段と、前記被検体内部のＸ線吸収係数分布を表す画像を、前記投影データに基づいて再構成する再構成手段と、前記再構成手段において行われる逆投影演算に際し、前記投影角度間隔の不均一性に応じたパワーの差異を補正するための補正係数を各々の逆投影データに掛けることにより前記再構成画像の劣化を抑制する抑制手段と、を具備する。
このため投影角度間隔の不均一性に応じたパワーの差異を補正して逆投影演算を行うことができる。
（２）本発明の画像再構成装置は、上記（１）に記載の装置であって、且つ前記補正係数は、基準となる投影角度間隔に対する所定の投影角度における投影角度間隔の比に基づいて算出されることを特徴とする。
このため、各々の投影角度毎の補正係数を算出できる。
【００１１】
このため投影角度間隔の不均一性により生じる再構成画像の劣化を抑制できる。
（２）本発明の画像再構成装置は、上記（１）に記載の装置であって、且つ前記抑制手段は、前記再構成手段において行われる逆投影演算に際し、前記投影角度間隔の不均一性に応じたパワーの差異を補正するための補正係数を各々の逆投影データに掛ける手段を具備することを特徴とする。
【００１２】
このため投影角度間隔の不均一性に応じたパワーの差異を補正して逆投影演算を行うことができる。
（３）本発明の画像再構成装置は、上記（２）に記載の装置であって、且つ前記補正係数は、基準となる投影角度間隔に対する所定の投影角度における投影角度間隔の比に基づいて算出されることを特徴とする。
【００１３】
このため、各々の投影角度毎の補正係数を算出できる。
（４）本発明の画像再構成装置は、上記（１）に記載の装置であって、且つ前記抑制手段は、前記収集手段により収集された不等角度間隔の投影データを補間処理することにより等角度間隔の投影データを作成する補間手段を具備し、前記再構成手段は、前記補間手段により作成された投影データに基づいて前記画像を再構成することを特徴とする。
このため、投影角度間隔の不均一性を再構成演算前に改善できる。
【００１４】
【発明の実施の形態】
以下、図面を参照しながら本発明の画像再構成装置の実施形態を説明する。
（第１実施形態）
図１及び図２は、本発明による画像再構成装置の第１実施形態に係る循環器用Ｘ線診断システムの構成を示す外観図及びブロック図である。図１に示すＸ線管球１からコーンビーム状に照射されたＸ線は図示しない被検体を透過し、二次元検出器２に入射する。この二次元検出器２は被検体を透過した透過Ｘ線の強度分布を検出する。Ｘ線管球１と二次元検出器２とは、Ｃ型の回転アーム３に接続されており、この回転アーム３を回転させながらＸ線を複数回曝射（又は連続的に曝射）することにより、被検体の周囲多方向からの投影データを収集するように構成されている。
【００１５】
回転アーム３をスライド回転可能に支持するホルダ３ａは、駆動機構３ｂ（不図示）を備えており、これにより回転アーム３は矢印Ｒ１の方向にスライド回転される。また、ホルダ３ａは天井に対し回転可能に支持されており、駆動機構３ｃ（不図示）によって矢印Ｒ２の方向に回転駆動される。なお、これら駆動機構３ｂ、３ｃは別途設けられた制御装置（不図示）によって制御されるものとなっている。
【００１６】
図２の（ａ）に示すように、本実施形態のＸ線診断システムはデータ収集部４、Ａ／Ｄ変換器５、投影データ用メモリ６、再構成部７、画像メモリ８、Ｄ／Ａ変換器９、そしてモニタ１０を具備している。なおデータ収集部４は回転アーム３の二次元検出器２を備えている。また、図２の（ｂ）は再構成部７の内部構成を示すブロック図であって、この再構成部７は歪補正部７１、対数（ｌｏｇ）変換部７２、フィルタリング部７３、係数補正部７４、そして逆投影部７５から構成されている。
【００１７】
図３は、アーム３の回転速度に伴い投影データの収集角度間隔が変化している様子を示す図である。
アーム３の回転速度は厳密には一定ではなく、このためＸ線の照射角度間隔も一定ではない。したがって、投影データの収集角度間隔はアーム３のスライド回転の始まりから終わりまでの間でダイナミックに変化する。例えば０度〜１５度においてはアーム３の回転速度は比較的遅く、このときの収集角度間隔はΔθ1 であり、１４５度〜１６０度においてはアーム３の回転速度は比較的速く、このときの収集角度間隔はΔθ2   である。なお、図３は理解のためにΔθ1   とΔθ2   の差を誇張して示したものであって、実際にはΔθ1   とΔθ2   の差はわずかである。
【００１８】
投影データを収集する場合、制御装置は駆動機構３ｂに対し、回転アーム３をスライド回転させるように制御信号を出力する。二次元検出器２は、このスライド回転中に一定の時間間隔で周期的にＸ線像を撮影してこれを出力する。
【００１９】
さて、以上述べたような回転アーム機構によるＸ線曝射により、例えば二次元検出器２に対する１２０パターン分の投影に応じた１８０度分の投影データが収集された場合を考える。実際には１８０度分の投影データだけでは完全な再構成像を得ることは出来ない。しかし回転アームを回転機構として有するこのような装置では、１８０度以上回転することは難しいことが多い。一方、このような再構成処理の目標は血管等の形状を観察したり計測したりすることであり、その内部のＣＴ値分布までを正確に求めることではない。従ってこのような不十分な投影データの影響によるアーチファクト（偽像）の発生は、例えば閾値処理等によって無視することができる。
【００２０】
二次元検出器２により検出された１２０パターン分の投影データは、イメージインテンシファイアを含むＴＶ系でアナログ信号に変換され、データ収集部４に送られた後、Ａ／Ｄ変換器５によってディジタル信号に変換され、投影データ用メモリ６に格納される。
【００２１】
再構成部７は、投影データ用メモリ６に格納された投影データを読み出して再構成処理を行う。ここでは、フィルター補正逆投影法による画像再構成を行うこととする。なお、再構成領域は全投影方向のＸ線束に内接する円筒若しくは球として定義される。
【００２２】
再構成部７は、Ｘ線管球１の回転軸と垂直をなす平面内において、例えば二次元検出器２の１つの検出素子の幅に投影される再構成領域中心部における長さｄによる離散化を行い、これにより得られる離散点のデータの集合によって表される再構成像を得る。なお、離散間隔は特定の長さｄに限定されるものではなく、装置やメーカーによって相違する。したがって、原則として装置毎で定義された離散間隔を用いる。さらに、二次元検出器の１素子の（回転軸方向の）幅を考慮し、回転軸方向について断層像内と同様の離散化を行う。かくして再構成領域は三次元の格子状に離散化される。この格子状に離散化されたデータの各素子は、一般にボクセルと称される。
【００２３】
ここで、再構成部７における再構成処理をより詳細に説明する。再構成部７は、二次元検出器２から得られた１２０フレームの投影データに対しフィルター補正逆投影法のための例えばShepp & Logan やRamachandranのような適当なフィルターを掛けてコンボリューションし、かかるコンボリューション結果に対し逆投影演算を行うことにより再構成データを得る。より具体的には、歪補正部７１によるイメージインテンシファイアの歪補正および対数（ｌｏｇ）変換部７２による投影データの対数変換を行い、フィルタリング部７３により、投影データに対しフィルター補正逆投影法に係る補正フィルタを用いてコンボリューションし、係数補正部７４により、コンボリューション結果に対し後述するパワー補正のための補正係数を掛ける。逆投影部７５は、補正係数が掛け合わされた投影データを逆投影して画像を再構成する。なお、ＴＶ系としてイメージインテンシファイアを具備しない場合、歪補正部７１は不要である。
【００２４】
特に本実施形態の再構成部７は、離散化された再構成領域の再構成演算を実施する際に、投影データの収集角度間隔の不均一性により生じる再構成画像の劣化を抑制する抑制手段を備えている。
【００２５】
この抑制手段は、投影データに対しフィルター補正逆投影法に係るコンボリューションフィルタを掛けた後、投影角度間隔の不均一性に伴うパワーの差異を補正するための補正係数を掛けて逆投影を行う。これによりパワーの差異に起因するアーチファクトを低減することができ、再構成画像の劣化を防止できる。
【００２６】
投影角度によっては、逆投影されるデータのパワーは異なる。例えば説明を簡単にするために、図４に示すようにΔθ1   ＝Δθ（Δθ  ０）間隔の投影データとΔθ2   （＝Δθ／２）間隔の投影データとが混在する場合を考える。３６０度すべてがΔθ間隔で収集されたデータであればアーチファクトは発生しない。しかしながら、Δθ／２間隔で収集された任意の方向Ｎ・Δθ（Ｎは整数）近傍においては、本来Δθ間隔で収集された場合の次の方向のデータ（Ｎ＋１）・Δθ以外に、（Ｎ＋１／２）・Δθが存在する。ここでΔθ  ０であるから、（Ｎ＋１／２）・ΔθのデータはＮ・Δθのデータとほぼ等しいと考えられる。このことからＮ・Δθと（Ｎ＋１／２）・Δθのデータを逆投影することは、Ｎ・Δθ投影データを２倍して逆投影することに等しく、言い替えれば投影角度Ｎ・Δθにおいて本来の投影データの２倍のデータが収集されたことになる。このように不等角度間隔で収集された投影方向では、その投影角度間隔によってパワーの異なる投影データが収集されていることと等しくなっているので、アーチファクトが生じる。
【００２７】
そこで本実施形態は、逆投影演算に際し、投影角度間隔ｄθの差異によるパワーの変化を補正するために、以下のような補正係数を逆投影データに掛けた後に逆投影演算を行うものとなっている。
【００２８】
Ｗ（ｄθ）＝ｄθ／ｄθ0（ｄθ0は基準となる投影角度間隔）…（１）
上記（１）式を適用することにより、投影角度間隔がΔθ0／２で収集された投影データには、逆投影演算に際し基準投影角度との比に応じた補正係数すなわち１／２が掛けられることになる。したがって、不均一な投影収集角度に起因して２倍となったパワーを適切な値に補正できる。
【００２９】
なお、投影角度間隔ｄθは、投影データ収集のたびに回転アーム３のスライド回転量から収集角度を検出するなどして計測により求めてもよいし、事前に計測された収集角度に基づいて各投影データの間の値を固定値としてテーブル等に記憶しておき既知の値として用いても良い。
【００３０】
次に、逆投影演算により得られた再構成データをコントラストスケールで正規化する。一般には逆投影により得られた再構成データを投影方向数Ｎ（一例として挙げると１２０）で除算するが、特に本実施形態においては、次式（２）で示されるようなＮdで除算する。
【００３１】
Ｎd＝θa／ｄθ0  …（２）
ここで、θaは投影データを収集した角度であって、本実施形態では例えば１８０度である。
【００３２】
再構成部７における画像再構成処理によって得られた画像、すなわち被検体の断層像は画像メモリ８に送られて記憶保持される。この断層像はモニタ１０に表示させることもできるし、また図示しない他の画像処理装置において閾値処理等により血管領域を抽出し、血管の表面画像に対する陰影付けのためのシェーディング処理等を行うことにより、任意の方向から観察した血管の表面画像を表示できる。
【００３３】
以上説明したように、本実施形態によれば、投影角度間隔の不均一性に応じたパワーの差異を補正しながらフィルター補正逆投影法に係る再構成処理を行うので、例えば代数的再構成法（ＡＲＴ：algebraic reconstruction technique）のような逐次近似法に比べて処理時間を大幅に増加させることなく短時間で、且つ効果的に再構成画像の劣化を抑制できる。
【００３４】
（第２実施形態）
次に本発明の第２実施形態を説明する。第２実施形態は上述した第１実施形態の変形に係り、循環器用Ｘ線診断システムにおいてディジタル・サブトラクション・アンギオグラフィ（ＤＳＡ：Digital Subtraction Angiography ）と称される撮影手法を実施するものである。図４は本実施形態の構成を示すブロック図である。
【００３５】
本実施形態では、被検体に対し造影剤を注入しながらＸ線撮影を行う。例えば血管のみを表す画像の再構成を行う場合は、投影データの収集を造影剤の注入前と注入後の２回行い、それぞれを第１（マスク像）の投影データ用メモリ１１１、第２（ライブ像）の投影データ用メモリ１１２に記憶する。双方のメモリに所要のデータが蓄積されると、サブトラクション処理部１２においてサブトラクション（差分化）処理が行われる。すなわち、マスク像の投影データとライブ像の投影データとをそれぞれｌｏｇ（対数）変換した後、ｌｏｇ変換結果同士を減算処理する。
【００３６】
サブトラクション処理結果は再構成部７に送られる。再構成部７における画像の再構成処理は第１実施形態のものとほぼ同様であるが若干異なっている。すなわち、サブトラクション処理部１２において既にｌｏｇ変換が施されているので、本実施形態の再構成部７はｌｏｇ変換部７２を備えない。
【００３７】
そして本実施形態は、第１実施形態と同様に、離散化された再構成領域の再構成演算を実施する際に、投影データの収集角度間隔の不均一性により生じる再構成画像の劣化を抑制する抑制手段を備えており、フィルタリング部７３により投影データに対しフィルター補正逆投影法に係る補正フィルタを用いてコンボリューションし、係数補正部７４により、コンボリューション結果に対し投影データの収集角度間隔の不均一性に伴うパワーの差異を補正するための補正係数を掛ける。逆投影部７５は補正係数が掛け合わされた投影データを逆投影して画像を再構成する。
【００３８】
したがって第２実施形態によれば、血管のみが抽出された画像（サブトラクション像）を得ることができるとともに、パワーの差異に起因するアーチファクトを低減することができ、再構成画像の劣化を防止できる。
【００３９】
（第３実施形態）
次に本発明の第３実施形態を説明する。上記第１及び第２実施形態では、投影角度間隔の不均一性に起因する再構成画像の劣化を抑制するために、コンボリューション結果に対し投影角度間隔の差異に伴うパワーの変化を補正するための補正係数を掛けて逆投影するという再構成処理の構成について説明した。これに対し第３実施形態の画像再構成装置は、投影角度間隔が均等な投影データを補間演算によって求め、この補間により得られた投影データをコンボリューションして逆投影する。その他の構成については、第１実施形態のものと同様でありその詳細な説明は省略する。
【００４０】
図６は本発明の第３実施形態の概略構成を示すブロック図である。図６（ａ）に示すように、投影データ２０と再構成部２１との間に、補間処理部２０が設けられている点において第１実施形態のものとは構成が異なっている。補間処理部２０の具体的な構成及び動作については後述する。なお、本実施形態の再構成部２１は、図６（ｂ）に示すように、第１実施形態に係る再構成部７の係数補正部７４を具備していない。
【００４１】
フィルター補正逆投影法において用いられる逆投影フィルターは、高周波強調フィルターであって、ディラックのδ関数を再構成する場合を考えると、当該δ関数が存在する位置以外においては、高周波強調フィルターのあらゆる方向からの逆投影によって正負が打ち消されて零になる。例えばある方向からは正の値が逆投影され、また別の方向からは負の値が逆投影されることになる。
【００４２】
しかしながら、不等角度間隔の投影データに対する逆投影においては、逆投影時に加算される正負のバランスが崩れてしまい、その結果、アーチファクトを生じてしまう。すなわち逆投影フィルターが掛かるバランスが崩れている。
【００４３】
そこで本実施形態においては、逆投影フィルターが掛かるバランスの適切化を図るために、不等角度間隔で収集された投影データから等角度間隔の投影データを補間演算により作成する。
【００４４】
特に本実施形態においては、補間処理部２０が１次補間を実施する場合について説明する。ここで、θ＝θ0の時点においてＮフレーム目のデータが収集され、θ＝θ1の時点において（Ｎ＋１）フレーム目のデータが収集されているとする。そして現在、θ＝θ2（θ0≦θ2≦θ1）の時点における投影データを補間によって作成する場合、θ＝θ2の時点における投影データｐ（ｘ，ｙ，θ2）は次式（３）に従って算出される。
【００４５】
ｐ（ｘ，ｙ，θ2）＝ｍ×ｐ（ｘ，ｙ，θ0）＋ｎ×ｐ（ｘ，ｙ，θ1）…（３）
ここで、（ｘ，ｙ）は画像上の任意のピクセルを示す座標系を示し、ｍ，ｎは次式（４）および（５）によって与えられる。
【００４６】
ｍ＝（θ1−θ2）／（θ1−θ0）  …（４）
ｎ＝（θ2−θ0）／（θ1−θ0）  …（５）
このような補間により新たに作成された等角度間隔の投影データは、再構成部２１に送られる。この実施形態では、補間関数として一次補間関数を用いた例を紹介したが、本発明は一次補間のみに限定されるものではなく、例えば三次のB-Spline関数のような高次関数を用いても良い。
【００４７】
再構成部２１における再構成処理では、離散化された再構成領域の再構成を行う。ここでの再構成処理は、当該処理が第１実施形態の係数補正処理を行わない点を除いては、第１実施形態の再構成処理と同様である。
【００４８】
以上述べたように第３実施形態は、不等角度間隔で収集された投影データから等角度間隔の投影データを補間演算により作成するものであり、逆投影フィルターが掛かるバランスの適切化を図ることができる。このため投影角度間隔の不均一性を再構成演算前に改善できる。したがって、第１実施形態と同様に短時間で、且つ効果的に再構成画像の劣化を抑制できる。
【００４９】
（第４実施形態）
次に本発明の第４実施形態を説明する。第４実施形態は上述した第３実施形態の変形に係り、上述したＤＳＡを実施するものである。図７は本実施形態の構成を示すブロック図である。
【００５０】
本実施形態では、被検体に対し造影剤を注入しながらＸ線撮影を行う。例えば血管のみを表す画像の再構成を行う場合は、投影データの収集を造影剤の注入前と注入後の２回行い、それぞれを第１（マスク像）の投影データ用メモリ１１１、第２（ライブ像）の投影データ用メモリ１１２に記憶する。双方のメモリに所要のデータが蓄積されると、サブトラクション処理部１２においてサブトラクション（差分化）処理が行われる。すなわち、マスク像の投影データとライブ像の投影データとをそれぞれｌｏｇ（対数）変換した後、ｌｏｇ変換結果同士を減算処理する。
【００５１】
サブトラクション処理結果は補間処理部２０に送られる。補間処理部２０における補間処理の内容及び、再構成部２１における画像の再構成処理は第３実施形態のものとほぼ同様であるが若干異なっている。すなわち、サブトラクション処理部１２において既にｌｏｇ変換が施されているので、本実施形態の再構成部７は第３実施形態のｌｏｇ変換部７２を備えない。
【００５２】
そして補間処理部２０は、第３実施形態と同様に、逆投影フィルターが掛かるバランスの適切化を図るために、不等角度間隔で収集された投影データから等角度間隔の投影データを補間演算により作成する。このため投影角度間隔の不均一性を再構成演算前に改善できる。したがって、短時間で、且つ効果的に再構成画像の劣化を抑制してサブトラクション像を得ることができる。
【００５３】
（第５実施形態）
上記した第１乃至第４の実施形態では、現在開発段階にある三次元ＣＴ装置において、撮影系として回転アーム（Ｃアームなど）からなる循環器用Ｘ線撮影装置を用いてデータ収集を行った後、画像再構成する場合について説明した。一般に普及している二次元ＣＴ（被検体の断面毎に撮影を行う）や、被検体の周囲でＸ線管球を回転させながら寝台を移動することにより、図８のように被検体を螺旋状にスキャンするヘリカルスキャンＣＴ（スパイラルスキャンＣＴ、螺旋スキャンＣＴとも呼ばれる）等に対しても本発明は実施可能である。
【００５４】
前者、すなわちヘリカルスキャンを行わないＣＴの場合、再構成領域は、全投影方向のＸ線束に内接する円として定義され、第１実施形態と同様に断面内における離散化を行なう。また後者の場合も、寝台の天板の移動方向に垂直な断面内での再構成空間は二次元ＣＴと同様に定義され、離散化も同様に行われる。ヘリカルスキャンＣＴは、基本的には任意のスライス位置における断面像を再構成することができる。
【００５５】
そして何れの場合においても第１乃至第４実施形態と同様に投影角度間隔の不均一性によって生じる再構成画像の劣化を抑制できる。
再構成方法の具体例としては、ヘリカルスキャンＣＴでは参考文献（東木  祐介他著、「ヘリカルスキャンに於ける補間再構成法の検討」、Med.Imag.Tech., Vol.8, No.3, 253-254(1990)）に示されている対向ビーム補間法や３６０度補間法により投影データを補間し、ファンビーム再構成法と同様の方法で再構成を行えば良い。
【００５６】
なお、本発明は上述した実施形態に限定されず、種々変形して実施可能である。例えば、本発明は図１に示した循環器用Ｘ線診断システムのみに限定されるものではなく、被検体周辺を回転しながらＸ線像の収集を行うことができる種々の循環器用Ｘ線診断システムに適用可能である。
【００５７】
【発明の効果】
本発明によれば、投影データ収集時における角度ピッチの不均一性を起因として発生する偽像（アーチファクト）による再構成画像の劣化を、効果的に、且つ処理時間を大幅に増やすことなく抑制し得る画像再構成装置を提供できる。
【００５８】
したがって、実際の再構成装置あるいは実用段階の実験装置等で使用されているフィルター補正逆投影法（Feldkamp法のような三次元再構成法を含む）において発生するアーチファクトを、演算時間をほとんど増やすことなく抑止できる。
【図面の簡単な説明】
【図１】本発明の画像再構成装置の第１実施形態に係る循環器用Ｘ線撮影システムの回転アームの外観を示す図。
【図２】画像の再構成部を含む上記第１実施形態の構成を示すブロック図。
【図３】上記第１実施形態に係り、投影データの収集角度間隔が変化している様子を示す図。
【図４】上記第１実施形態に係り、投影データのパワーの差異を説明するための模式図。
【図５】本発明の画像再構成装置の第２実施形態に係り、ＤＳＡを実施する場合の構成を示すブロック図。
【図６】本発明の画像再構成装置の第３実施形態の構成を示すブロック図。
【図７】本発明の画像再構成装置の第４実施形態に係り、ＤＳＡを実施する場合の構成を示すブロック図。
【図８】本発明の第５実施形態に係るヘリカルスキャン方式を説明するための図。
【符号の説明】
１…Ｘ線管球
２…二次元検出器
３…回転アーム
４…データ収集部
５…Ａ／Ｄ変換器
６…投影データ用メモリ
７…再構成部
８…画像メモリ
９…Ｄ／Ａ変換器
１０…モニタ
７１…歪補正部
７２…対数（ｌｏｇ）変換部
７３…フィルタリング部
７４…係数補正部
７５…逆投影部
【図１】
0004047420.tif 000002
【図２】
0004047420.tif 000003
【図３】
0004047420.tif 000004
【図４】
0004047420.tif 000005
【図５】
0004047420.tif 000006
【図６】
0004047420.tif 000007
【図７】
0004047420.tif 000008
【図８】
0004047420.tif 000009