JPH11113892A - 画像再構成装置 - Google Patents
画像再構成装置Info
- Publication number
- JPH11113892A JPH11113892A JP9275864A JP27586497A JPH11113892A JP H11113892 A JPH11113892 A JP H11113892A JP 9275864 A JP9275864 A JP 9275864A JP 27586497 A JP27586497 A JP 27586497A JP H11113892 A JPH11113892 A JP H11113892A
- Authority
- JP
- Japan
- Prior art keywords
- projection
- image
- reconstruction
- projection data
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000012937 correction Methods 0.000 claims abstract description 36
- 238000012545 processing Methods 0.000 claims description 32
- 230000006866 deterioration Effects 0.000 claims description 13
- 230000001629 suppression Effects 0.000 claims description 7
- 238000010521 absorption reaction Methods 0.000 claims description 4
- 238000000034 method Methods 0.000 abstract description 39
- 230000015654 memory Effects 0.000 abstract description 13
- 230000015556 catabolic process Effects 0.000 abstract description 4
- 238000006731 degradation reaction Methods 0.000 abstract description 4
- 230000001678 irradiating effect Effects 0.000 abstract description 4
- 230000000452 restraining effect Effects 0.000 abstract 2
- 238000002591 computed tomography Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 12
- 238000006243 chemical reaction Methods 0.000 description 9
- 210000004204 blood vessel Anatomy 0.000 description 7
- 238000013480 data collection Methods 0.000 description 7
- 230000007246 mechanism Effects 0.000 description 7
- 239000011295 pitch Substances 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 210000000056 organ Anatomy 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 239000002872 contrast media Substances 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 238000002059 diagnostic imaging Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000002583 angiography Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- HEFNNWSXXWATRW-UHFFFAOYSA-N Ibuprofen Chemical compound CC(C)CC1=CC=C(C(C)C(O)=O)C=C1 HEFNNWSXXWATRW-UHFFFAOYSA-N 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
(57)【要約】
【課題】本発明は投影データ収集時における角度ピッチ
の不均一性を起因として発生する偽像による再構成画像
の劣化を、処理時間を大幅に増やすことなく効果的に抑
制し得る画像再構成装置を提供することを目的とする。 【解決手段】再構成部７は、離散化された再構成領域の
再構成演算を実施する際に、Ｘ線の照射角度の不均一性
により生じる再構成画像の劣化を抑制する抑制手段を備
えている。この抑制手段は、投影データに対しフィルタ
ー補正逆投影法に係るコンボリューションフィルタを掛
けた後、投影角度間隔の不均一性に伴うパワーの差異を
補正するための補正係数を掛けて逆投影する。これによ
りパワーの差異に起因するアーチファクトを低減するこ
とができ、再構成画像の劣化を防止できる。
の不均一性を起因として発生する偽像による再構成画像
の劣化を、処理時間を大幅に増やすことなく効果的に抑
制し得る画像再構成装置を提供することを目的とする。 【解決手段】再構成部７は、離散化された再構成領域の
再構成演算を実施する際に、Ｘ線の照射角度の不均一性
により生じる再構成画像の劣化を抑制する抑制手段を備
えている。この抑制手段は、投影データに対しフィルタ
ー補正逆投影法に係るコンボリューションフィルタを掛
けた後、投影角度間隔の不均一性に伴うパワーの差異を
補正するための補正係数を掛けて逆投影する。これによ
りパワーの差異に起因するアーチファクトを低減するこ
とができ、再構成画像の劣化を防止できる。
Description
【０００１】
【発明の属する技術分野】本発明は、被検体内部のＸ線
の吸収係数に応じた分布画像を、異なる角度からの複数
の投影データに基づいて再構成する画像再構成装置に関
する。
の吸収係数に応じた分布画像を、異なる角度からの複数
の投影データに基づいて再構成する画像再構成装置に関
する。
【０００２】
【従来の技術】被検体を透過した放射線の線量分布を複
数方向で計測した投影データから、被検体のＸ線吸収係
数の分布画像を再構成する方法としては、参考文献
（「医用画像工学ハンドブック」日本医用画像工学会監
修、日本医用画像工学会・医用画像工学ハンドブック編
集委員会編集、篠原出版、Ｉ．基礎論理、第５章、画像
再構成）に幾つか紹介されている。その中でも直接的
（解析的）再構成法の１つであるフィルター補正逆投影
法（Filtered Backprojection ）がよく知られている。
数方向で計測した投影データから、被検体のＸ線吸収係
数の分布画像を再構成する方法としては、参考文献
（「医用画像工学ハンドブック」日本医用画像工学会監
修、日本医用画像工学会・医用画像工学ハンドブック編
集委員会編集、篠原出版、Ｉ．基礎論理、第５章、画像
再構成）に幾つか紹介されている。その中でも直接的
（解析的）再構成法の１つであるフィルター補正逆投影
法（Filtered Backprojection ）がよく知られている。
【０００３】現在のＣＴ（コンピュータ断層撮影）装置
において適用される画像再構成法は、ハードウェア化を
容易にするために多少の近似が行われていたり、空間フ
ィルターを補正フィルターに加味している等の相違点を
有しているが、基本的にはこののような方法を採用して
いる。
において適用される画像再構成法は、ハードウェア化を
容易にするために多少の近似が行われていたり、空間フ
ィルターを補正フィルターに加味している等の相違点を
有しているが、基本的にはこののような方法を採用して
いる。
【０００４】この画像再構成法は、投影データが等角度
間隔、すなわち各々が等しいピッチで収集されることを
前提としている。しかしながら、Ｘ線を用いて画像再構
成を行う装置として最もメジャーなＣＴ装置において
も、Ｘ線管球の回転速度あるいはＸ線の照射角度にはＸ
線管球の重み等に起因する微小なブレが生じるため、現
実には不等角度間隔で投影データが収集される場合が起
こり得る。さらに、近年話題になっている三次元ＣＴ
を、回転アーム機構を有する循環器Ｘ線撮影装置によっ
て実現する場合においては、このような影響は顕著とな
る。具体的に言えば、回転開始から加速している期間及
び減速から回転終了までの期間においては、回転速度が
動的に（ダイナミックに）変化する。ところがＸ線曝射
間隔は常に等時間間隔である。このため投影データ収集
の角度ピッチがこの間にダイナミックに変化してしま
う。
間隔、すなわち各々が等しいピッチで収集されることを
前提としている。しかしながら、Ｘ線を用いて画像再構
成を行う装置として最もメジャーなＣＴ装置において
も、Ｘ線管球の回転速度あるいはＸ線の照射角度にはＸ
線管球の重み等に起因する微小なブレが生じるため、現
実には不等角度間隔で投影データが収集される場合が起
こり得る。さらに、近年話題になっている三次元ＣＴ
を、回転アーム機構を有する循環器Ｘ線撮影装置によっ
て実現する場合においては、このような影響は顕著とな
る。具体的に言えば、回転開始から加速している期間及
び減速から回転終了までの期間においては、回転速度が
動的に（ダイナミックに）変化する。ところがＸ線曝射
間隔は常に等時間間隔である。このため投影データ収集
の角度ピッチがこの間にダイナミックに変化してしま
う。
【０００５】また、たとえ回転速度がほぼ一定であって
も、Ｘ線管球や検出器の重さによる回転速度への直接的
な影響や、回転アームのたわみ等による間接的な影響に
よって、投影データ収集の角度ピッチが変化することも
考えられる。
も、Ｘ線管球や検出器の重さによる回転速度への直接的
な影響や、回転アームのたわみ等による間接的な影響に
よって、投影データ収集の角度ピッチが変化することも
考えられる。
【０００６】したがって、このような投影データ収集に
おける角度ピッチの不均一性により、フィルター補正逆
投影法による再構成処理により得られた再構成像が著し
く劣化してしまうという問題点がある。
おける角度ピッチの不均一性により、フィルター補正逆
投影法による再構成処理により得られた再構成像が著し
く劣化してしまうという問題点がある。
【０００７】上記した三次元ＣＴに係る画像再構成法に
参考文献（L.A.Feldkamp他著、"Practical cone-beam a
lgorithm" 、J. Opt. Soc. Am. A, No.6, 612-619(198
4) ）あるいは（大石 悟他著、「コーンビーム投影系
を用いた三次元再構成」、Med.Imag.Tech.,Vol.8, No.
2, 127-130(1990) ）に記載のFeldkampのフィルター補
正逆投影法がある。この方法は、Feldkamp等によって提
案され、三次元再構成法において最も有名で且つ実用的
な方法の一つとされており、投影データに対してコンボ
リューションフィルターを掛けた後、コーンビームの広
がり補正を行いながら三次元的な逆投影演算を行うもの
である。
参考文献（L.A.Feldkamp他著、"Practical cone-beam a
lgorithm" 、J. Opt. Soc. Am. A, No.6, 612-619(198
4) ）あるいは（大石 悟他著、「コーンビーム投影系
を用いた三次元再構成」、Med.Imag.Tech.,Vol.8, No.
2, 127-130(1990) ）に記載のFeldkampのフィルター補
正逆投影法がある。この方法は、Feldkamp等によって提
案され、三次元再構成法において最も有名で且つ実用的
な方法の一つとされており、投影データに対してコンボ
リューションフィルターを掛けた後、コーンビームの広
がり補正を行いながら三次元的な逆投影演算を行うもの
である。
【０００８】しかしながらこのFeldkampのフィルター補
正逆投影法においても上述した問題は解消されることな
く、投影データ収集における不均一な角度ピッチが再構
成像を劣化させるように働いてしまう。
正逆投影法においても上述した問題は解消されることな
く、投影データ収集における不均一な角度ピッチが再構
成像を劣化させるように働いてしまう。
【０００９】
【発明が解決しようとする課題】本発明は上述した事情
を考慮してなされたものであり、その目的とするところ
は、投影データ収集時における角度ピッチの不均一性を
起因として発生する偽像（アーチファクト）による再構
成画像の劣化を、処理時間を大幅に増やすことなく効果
的に抑制し得る画像再構成装置を提供することにある。
を考慮してなされたものであり、その目的とするところ
は、投影データ収集時における角度ピッチの不均一性を
起因として発生する偽像（アーチファクト）による再構
成画像の劣化を、処理時間を大幅に増やすことなく効果
的に抑制し得る画像再構成装置を提供することにある。
【００１０】
【課題を解決するための手段】上記課題を解決し目的を
達成するために本発明は以下の手段を備えている。 （１）本発明の画像再構成装置は、被検体に対し多方向
からＸ線を照射することによって所定の投影角度間隔に
おける複数の投影データを収集する収集手段と、前記被
検体内部のＸ線吸収係数分布を表す画像を、前記投影デ
ータに基づいて再構成する再構成手段と、前記投影角度
間隔の不均一性により生じる前記再構成画像の劣化を抑
制する抑制手段と、を具備する。
達成するために本発明は以下の手段を備えている。 （１）本発明の画像再構成装置は、被検体に対し多方向
からＸ線を照射することによって所定の投影角度間隔に
おける複数の投影データを収集する収集手段と、前記被
検体内部のＸ線吸収係数分布を表す画像を、前記投影デ
ータに基づいて再構成する再構成手段と、前記投影角度
間隔の不均一性により生じる前記再構成画像の劣化を抑
制する抑制手段と、を具備する。
【００１１】このため投影角度間隔の不均一性により生
じる再構成画像の劣化を抑制できる。 （２）本発明の画像再構成装置は、上記（１）に記載の
装置であって、且つ前記抑制手段は、前記再構成手段に
おいて行われる逆投影演算に際し、前記投影角度間隔の
不均一性に応じたパワーの差異を補正するための補正係
数を各々の逆投影データに掛ける手段を具備することを
特徴とする。
じる再構成画像の劣化を抑制できる。 （２）本発明の画像再構成装置は、上記（１）に記載の
装置であって、且つ前記抑制手段は、前記再構成手段に
おいて行われる逆投影演算に際し、前記投影角度間隔の
不均一性に応じたパワーの差異を補正するための補正係
数を各々の逆投影データに掛ける手段を具備することを
特徴とする。
【００１２】このため投影角度間隔の不均一性に応じた
パワーの差異を補正して逆投影演算を行うことができ
る。 （３）本発明の画像再構成装置は、上記（２）に記載の
装置であって、且つ前記補正係数は、基準となる投影角
度間隔に対する所定の投影角度における投影角度間隔の
比に基づいて算出されることを特徴とする。
パワーの差異を補正して逆投影演算を行うことができ
る。 （３）本発明の画像再構成装置は、上記（２）に記載の
装置であって、且つ前記補正係数は、基準となる投影角
度間隔に対する所定の投影角度における投影角度間隔の
比に基づいて算出されることを特徴とする。
【００１３】このため、各々の投影角度毎の補正係数を
算出できる。 （４）本発明の画像再構成装置は、上記（１）に記載の
装置であって、且つ前記抑制手段は、前記収集手段によ
り収集された不等角度間隔の投影データを補間処理する
ことにより等角度間隔の投影データを作成する補間手段
を具備し、前記再構成手段は、前記補間手段により作成
された投影データに基づいて前記画像を再構成すること
を特徴とする。このため、投影角度間隔の不均一性を再
構成演算前に改善できる。
算出できる。 （４）本発明の画像再構成装置は、上記（１）に記載の
装置であって、且つ前記抑制手段は、前記収集手段によ
り収集された不等角度間隔の投影データを補間処理する
ことにより等角度間隔の投影データを作成する補間手段
を具備し、前記再構成手段は、前記補間手段により作成
された投影データに基づいて前記画像を再構成すること
を特徴とする。このため、投影角度間隔の不均一性を再
構成演算前に改善できる。
【００１４】
【発明の実施の形態】以下、図面を参照しながら本発明
の画像再構成装置の実施形態を説明する。 （第１実施形態）図１及び図２は、本発明による画像再
構成装置の第１実施形態に係る循環器用Ｘ線診断システ
ムの構成を示す外観図及びブロック図である。図１に示
すＸ線管球１からコーンビーム状に照射されたＸ線は図
示しない被検体を透過し、二次元検出器２に入射する。
この二次元検出器２は被検体を透過した透過Ｘ線の強度
分布を検出する。Ｘ線管球１と二次元検出器２とは、Ｃ
型の回転アーム３に接続されており、この回転アーム３
を回転させながらＸ線を複数回曝射（又は連続的に曝
射）することにより、被検体の周囲多方向からの投影デ
ータを収集するように構成されている。
の画像再構成装置の実施形態を説明する。 （第１実施形態）図１及び図２は、本発明による画像再
構成装置の第１実施形態に係る循環器用Ｘ線診断システ
ムの構成を示す外観図及びブロック図である。図１に示
すＸ線管球１からコーンビーム状に照射されたＸ線は図
示しない被検体を透過し、二次元検出器２に入射する。
この二次元検出器２は被検体を透過した透過Ｘ線の強度
分布を検出する。Ｘ線管球１と二次元検出器２とは、Ｃ
型の回転アーム３に接続されており、この回転アーム３
を回転させながらＸ線を複数回曝射（又は連続的に曝
射）することにより、被検体の周囲多方向からの投影デ
ータを収集するように構成されている。
【００１５】回転アーム３をスライド回転可能に支持す
るホルダ３ａは、駆動機構３ｂ（不図示）を備えてお
り、これにより回転アーム３は矢印Ｒ１の方向にスライ
ド回転される。また、ホルダ３ａは天井に対し回転可能
に支持されており、駆動機構３ｃ（不図示）によって矢
印Ｒ２の方向に回転駆動される。なお、これら駆動機構
３ｂ、３ｃは別途設けられた制御装置（不図示）によっ
て制御されるものとなっている。
るホルダ３ａは、駆動機構３ｂ（不図示）を備えてお
り、これにより回転アーム３は矢印Ｒ１の方向にスライ
ド回転される。また、ホルダ３ａは天井に対し回転可能
に支持されており、駆動機構３ｃ（不図示）によって矢
印Ｒ２の方向に回転駆動される。なお、これら駆動機構
３ｂ、３ｃは別途設けられた制御装置（不図示）によっ
て制御されるものとなっている。
【００１６】図２の（ａ）に示すように、本実施形態の
Ｘ線診断システムはデータ収集部４、Ａ／Ｄ変換器５、
投影データ用メモリ６、再構成部７、画像メモリ８、Ｄ
／Ａ変換器９、そしてモニタ１０を具備している。なお
データ収集部４は回転アーム３の二次元検出器２を備え
ている。また、図２の（ｂ）は再構成部７の内部構成を
示すブロック図であって、この再構成部７は歪補正部７
１、対数（ｌｏｇ）変換部７２、フィルタリング部７
３、係数補正部７４、そして逆投影部７５から構成され
ている。
Ｘ線診断システムはデータ収集部４、Ａ／Ｄ変換器５、
投影データ用メモリ６、再構成部７、画像メモリ８、Ｄ
／Ａ変換器９、そしてモニタ１０を具備している。なお
データ収集部４は回転アーム３の二次元検出器２を備え
ている。また、図２の（ｂ）は再構成部７の内部構成を
示すブロック図であって、この再構成部７は歪補正部７
１、対数（ｌｏｇ）変換部７２、フィルタリング部７
３、係数補正部７４、そして逆投影部７５から構成され
ている。
【００１７】図３は、アーム３の回転速度に伴い投影デ
ータの収集角度間隔が変化している様子を示す図であ
る。アーム３の回転速度は厳密には一定ではなく、この
ためＸ線の照射角度間隔も一定ではない。したがって、
投影データの収集角度間隔はアーム３のスライド回転の
始まりから終わりまでの間でダイナミックに変化する。
例えば０度〜１５度においてはアーム３の回転速度は比
較的遅く、このときの収集角度間隔はΔθ1であり、１
４５度〜１６０度においてはアーム３の回転速度は比較
的速く、このときの収集角度間隔はΔθ2 である。な
お、図３は理解のためにΔθ1 とΔθ2 の差を誇張
して示したものであって、実際にはΔθ1 とΔθ2
の差はわずかである。
ータの収集角度間隔が変化している様子を示す図であ
る。アーム３の回転速度は厳密には一定ではなく、この
ためＸ線の照射角度間隔も一定ではない。したがって、
投影データの収集角度間隔はアーム３のスライド回転の
始まりから終わりまでの間でダイナミックに変化する。
例えば０度〜１５度においてはアーム３の回転速度は比
較的遅く、このときの収集角度間隔はΔθ1であり、１
４５度〜１６０度においてはアーム３の回転速度は比較
的速く、このときの収集角度間隔はΔθ2 である。な
お、図３は理解のためにΔθ1 とΔθ2 の差を誇張
して示したものであって、実際にはΔθ1 とΔθ2
の差はわずかである。
【００１８】投影データを収集する場合、制御装置は駆
動機構３ｂに対し、回転アーム３をスライド回転させる
ように制御信号を出力する。二次元検出器２は、このス
ライド回転中に一定の時間間隔で周期的にＸ線像を撮影
してこれを出力する。
動機構３ｂに対し、回転アーム３をスライド回転させる
ように制御信号を出力する。二次元検出器２は、このス
ライド回転中に一定の時間間隔で周期的にＸ線像を撮影
してこれを出力する。
【００１９】さて、以上述べたような回転アーム機構に
よるＸ線曝射により、例えば二次元検出器２に対する１
２０パターン分の投影に応じた１８０度分の投影データ
が収集された場合を考える。実際には１８０度分の投影
データだけでは完全な再構成像を得ることは出来ない。
しかし回転アームを回転機構として有するこのような装
置では、１８０度以上回転することは難しいことが多
い。一方、このような再構成処理の目標は血管等の形状
を観察したり計測したりすることであり、その内部のＣ
Ｔ値分布までを正確に求めることではない。従ってこの
ような不十分な投影データの影響によるアーチファクト
（偽像）の発生は、例えば閾値処理等によって無視する
ことができる。
よるＸ線曝射により、例えば二次元検出器２に対する１
２０パターン分の投影に応じた１８０度分の投影データ
が収集された場合を考える。実際には１８０度分の投影
データだけでは完全な再構成像を得ることは出来ない。
しかし回転アームを回転機構として有するこのような装
置では、１８０度以上回転することは難しいことが多
い。一方、このような再構成処理の目標は血管等の形状
を観察したり計測したりすることであり、その内部のＣ
Ｔ値分布までを正確に求めることではない。従ってこの
ような不十分な投影データの影響によるアーチファクト
（偽像）の発生は、例えば閾値処理等によって無視する
ことができる。
【００２０】二次元検出器２により検出された１２０パ
ターン分の投影データは、イメージインテンシファイア
を含むＴＶ系でアナログ信号に変換され、データ収集部
４に送られた後、Ａ／Ｄ変換器５によってディジタル信
号に変換され、投影データ用メモリ６に格納される。
ターン分の投影データは、イメージインテンシファイア
を含むＴＶ系でアナログ信号に変換され、データ収集部
４に送られた後、Ａ／Ｄ変換器５によってディジタル信
号に変換され、投影データ用メモリ６に格納される。
【００２１】再構成部７は、投影データ用メモリ６に格
納された投影データを読み出して再構成処理を行う。こ
こでは、フィルター補正逆投影法による画像再構成を行
うこととする。なお、再構成領域は全投影方向のＸ線束
に内接する円筒若しくは球として定義される。
納された投影データを読み出して再構成処理を行う。こ
こでは、フィルター補正逆投影法による画像再構成を行
うこととする。なお、再構成領域は全投影方向のＸ線束
に内接する円筒若しくは球として定義される。
【００２２】再構成部７は、Ｘ線管球１の回転軸と垂直
をなす平面内において、例えば二次元検出器２の１つの
検出素子の幅に投影される再構成領域中心部における長
さｄによる離散化を行い、これにより得られる離散点の
データの集合によって表される再構成像を得る。なお、
離散間隔は特定の長さｄに限定されるものではなく、装
置やメーカーによって相違する。したがって、原則とし
て装置毎で定義された離散間隔を用いる。さらに、二次
元検出器の１素子の（回転軸方向の）幅を考慮し、回転
軸方向について断層像内と同様の離散化を行う。かくし
て再構成領域は三次元の格子状に離散化される。この格
子状に離散化されたデータの各素子は、一般にボクセル
と称される。
をなす平面内において、例えば二次元検出器２の１つの
検出素子の幅に投影される再構成領域中心部における長
さｄによる離散化を行い、これにより得られる離散点の
データの集合によって表される再構成像を得る。なお、
離散間隔は特定の長さｄに限定されるものではなく、装
置やメーカーによって相違する。したがって、原則とし
て装置毎で定義された離散間隔を用いる。さらに、二次
元検出器の１素子の（回転軸方向の）幅を考慮し、回転
軸方向について断層像内と同様の離散化を行う。かくし
て再構成領域は三次元の格子状に離散化される。この格
子状に離散化されたデータの各素子は、一般にボクセル
と称される。
【００２３】ここで、再構成部７における再構成処理を
より詳細に説明する。再構成部７は、二次元検出器２か
ら得られた１２０フレームの投影データに対しフィルタ
ー補正逆投影法のための例えばShepp & Logan やRamach
andranのような適当なフィルターを掛けてコンボリュー
ションし、かかるコンボリューション結果に対し逆投影
演算を行うことにより再構成データを得る。より具体的
には、歪補正部７１によるイメージインテンシファイア
の歪補正および対数（ｌｏｇ）変換部７２による投影デ
ータの対数変換を行い、フィルタリング部７３により、
投影データに対しフィルター補正逆投影法に係る補正フ
ィルタを用いてコンボリューションし、係数補正部７４
により、コンボリューション結果に対し後述するパワー
補正のための補正係数を掛ける。逆投影部７５は、補正
係数が掛け合わされた投影データを逆投影して画像を再
構成する。なお、ＴＶ系としてイメージインテンシファ
イアを具備しない場合、歪補正部７１は不要である。
より詳細に説明する。再構成部７は、二次元検出器２か
ら得られた１２０フレームの投影データに対しフィルタ
ー補正逆投影法のための例えばShepp & Logan やRamach
andranのような適当なフィルターを掛けてコンボリュー
ションし、かかるコンボリューション結果に対し逆投影
演算を行うことにより再構成データを得る。より具体的
には、歪補正部７１によるイメージインテンシファイア
の歪補正および対数（ｌｏｇ）変換部７２による投影デ
ータの対数変換を行い、フィルタリング部７３により、
投影データに対しフィルター補正逆投影法に係る補正フ
ィルタを用いてコンボリューションし、係数補正部７４
により、コンボリューション結果に対し後述するパワー
補正のための補正係数を掛ける。逆投影部７５は、補正
係数が掛け合わされた投影データを逆投影して画像を再
構成する。なお、ＴＶ系としてイメージインテンシファ
イアを具備しない場合、歪補正部７１は不要である。
【００２４】特に本実施形態の再構成部７は、離散化さ
れた再構成領域の再構成演算を実施する際に、投影デー
タの収集角度間隔の不均一性により生じる再構成画像の
劣化を抑制する抑制手段を備えている。
れた再構成領域の再構成演算を実施する際に、投影デー
タの収集角度間隔の不均一性により生じる再構成画像の
劣化を抑制する抑制手段を備えている。
【００２５】この抑制手段は、投影データに対しフィル
ター補正逆投影法に係るコンボリューションフィルタを
掛けた後、投影角度間隔の不均一性に伴うパワーの差異
を補正するための補正係数を掛けて逆投影を行う。これ
によりパワーの差異に起因するアーチファクトを低減す
ることができ、再構成画像の劣化を防止できる。
ター補正逆投影法に係るコンボリューションフィルタを
掛けた後、投影角度間隔の不均一性に伴うパワーの差異
を補正するための補正係数を掛けて逆投影を行う。これ
によりパワーの差異に起因するアーチファクトを低減す
ることができ、再構成画像の劣化を防止できる。
【００２６】投影角度によっては、逆投影されるデータ
のパワーは異なる。例えば説明を簡単にするために、図
４に示すようにΔθ1 ＝Δθ（Δθ ０）間隔の投影
データとΔθ2 （＝Δθ／２）間隔の投影データとが
混在する場合を考える。３６０度すべてがΔθ間隔で収
集されたデータであればアーチファクトは発生しない。
しかしながら、Δθ／２間隔で収集された任意の方向Ｎ
・Δθ（Ｎは整数）近傍においては、本来Δθ間隔で収
集された場合の次の方向のデータ（Ｎ＋１）・Δθ以外
に、（Ｎ＋１／２）・Δθが存在する。ここでΔθ ０
であるから、（Ｎ＋１／２）・ΔθのデータはＮ・Δθ
のデータとほぼ等しいと考えられる。このことからＮ・
Δθと（Ｎ＋１／２）・Δθのデータを逆投影すること
は、Ｎ・Δθ投影データを２倍して逆投影することに等
しく、言い替えれば投影角度Ｎ・Δθにおいて本来の投
影データの２倍のデータが収集されたことになる。この
ように不等角度間隔で収集された投影方向では、その投
影角度間隔によってパワーの異なる投影データが収集さ
れていることと等しくなっているので、アーチファクト
が生じる。
のパワーは異なる。例えば説明を簡単にするために、図
４に示すようにΔθ1 ＝Δθ（Δθ ０）間隔の投影
データとΔθ2 （＝Δθ／２）間隔の投影データとが
混在する場合を考える。３６０度すべてがΔθ間隔で収
集されたデータであればアーチファクトは発生しない。
しかしながら、Δθ／２間隔で収集された任意の方向Ｎ
・Δθ（Ｎは整数）近傍においては、本来Δθ間隔で収
集された場合の次の方向のデータ（Ｎ＋１）・Δθ以外
に、（Ｎ＋１／２）・Δθが存在する。ここでΔθ ０
であるから、（Ｎ＋１／２）・ΔθのデータはＮ・Δθ
のデータとほぼ等しいと考えられる。このことからＮ・
Δθと（Ｎ＋１／２）・Δθのデータを逆投影すること
は、Ｎ・Δθ投影データを２倍して逆投影することに等
しく、言い替えれば投影角度Ｎ・Δθにおいて本来の投
影データの２倍のデータが収集されたことになる。この
ように不等角度間隔で収集された投影方向では、その投
影角度間隔によってパワーの異なる投影データが収集さ
れていることと等しくなっているので、アーチファクト
が生じる。
【００２７】そこで本実施形態は、逆投影演算に際し、
投影角度間隔ｄθの差異によるパワーの変化を補正する
ために、以下のような補正係数を逆投影データに掛けた
後に逆投影演算を行うものとなっている。
投影角度間隔ｄθの差異によるパワーの変化を補正する
ために、以下のような補正係数を逆投影データに掛けた
後に逆投影演算を行うものとなっている。
【００２８】 Ｗ（ｄθ）＝ｄθ／ｄθ0（ｄθ0は基準となる投影角度間隔）…（１） 上記（１）式を適用することにより、投影角度間隔がΔ
θ0／２で収集された投影データには、逆投影演算に際
し基準投影角度との比に応じた補正係数すなわち１／２
が掛けられることになる。したがって、不均一な投影収
集角度に起因して２倍となったパワーを適切な値に補正
できる。
θ0／２で収集された投影データには、逆投影演算に際
し基準投影角度との比に応じた補正係数すなわち１／２
が掛けられることになる。したがって、不均一な投影収
集角度に起因して２倍となったパワーを適切な値に補正
できる。
【００２９】なお、投影角度間隔ｄθは、投影データ収
集のたびに回転アーム３のスライド回転量から収集角度
を検出するなどして計測により求めてもよいし、事前に
計測された収集角度に基づいて各投影データの間の値を
固定値としてテーブル等に記憶しておき既知の値として
用いても良い。
集のたびに回転アーム３のスライド回転量から収集角度
を検出するなどして計測により求めてもよいし、事前に
計測された収集角度に基づいて各投影データの間の値を
固定値としてテーブル等に記憶しておき既知の値として
用いても良い。
【００３０】次に、逆投影演算により得られた再構成デ
ータをコントラストスケールで正規化する。一般には逆
投影により得られた再構成データを投影方向数Ｎ（一例
として挙げると１２０）で除算するが、特に本実施形態
においては、次式（２）で示されるようなＮdで除算す
る。
ータをコントラストスケールで正規化する。一般には逆
投影により得られた再構成データを投影方向数Ｎ（一例
として挙げると１２０）で除算するが、特に本実施形態
においては、次式（２）で示されるようなＮdで除算す
る。
【００３１】Ｎd＝θa／ｄθ0 …（２） ここで、θaは投影データを収集した角度であって、本
実施形態では例えば１８０度である。
実施形態では例えば１８０度である。
【００３２】再構成部７における画像再構成処理によっ
て得られた画像、すなわち被検体の断層像は画像メモリ
８に送られて記憶保持される。この断層像はモニタ１０
に表示させることもできるし、また図示しない他の画像
処理装置において閾値処理等により血管領域を抽出し、
血管の表面画像に対する陰影付けのためのシェーディン
グ処理等を行うことにより、任意の方向から観察した血
管の表面画像を表示できる。
て得られた画像、すなわち被検体の断層像は画像メモリ
８に送られて記憶保持される。この断層像はモニタ１０
に表示させることもできるし、また図示しない他の画像
処理装置において閾値処理等により血管領域を抽出し、
血管の表面画像に対する陰影付けのためのシェーディン
グ処理等を行うことにより、任意の方向から観察した血
管の表面画像を表示できる。
【００３３】以上説明したように、本実施形態によれ
ば、投影角度間隔の不均一性に応じたパワーの差異を補
正しながらフィルター補正逆投影法に係る再構成処理を
行うので、例えば代数的再構成法（ＡＲＴ：algebraic
reconstruction technique）のような逐次近似法に比べ
て処理時間を大幅に増加させることなく短時間で、且つ
効果的に再構成画像の劣化を抑制できる。
ば、投影角度間隔の不均一性に応じたパワーの差異を補
正しながらフィルター補正逆投影法に係る再構成処理を
行うので、例えば代数的再構成法（ＡＲＴ：algebraic
reconstruction technique）のような逐次近似法に比べ
て処理時間を大幅に増加させることなく短時間で、且つ
効果的に再構成画像の劣化を抑制できる。
【００３４】（第２実施形態）次に本発明の第２実施形
態を説明する。第２実施形態は上述した第１実施形態の
変形に係り、循環器用Ｘ線診断システムにおいてディジ
タル・サブトラクション・アンギオグラフィ（ＤＳＡ：
Digital Subtraction Angiography ）と称される撮影手
法を実施するものである。図４は本実施形態の構成を示
すブロック図である。
態を説明する。第２実施形態は上述した第１実施形態の
変形に係り、循環器用Ｘ線診断システムにおいてディジ
タル・サブトラクション・アンギオグラフィ（ＤＳＡ：
Digital Subtraction Angiography ）と称される撮影手
法を実施するものである。図４は本実施形態の構成を示
すブロック図である。
【００３５】本実施形態では、被検体に対し造影剤を注
入しながらＸ線撮影を行う。例えば血管のみを表す画像
の再構成を行う場合は、投影データの収集を造影剤の注
入前と注入後の２回行い、それぞれを第１（マスク像）
の投影データ用メモリ１１１、第２（ライブ像）の投影
データ用メモリ１１２に記憶する。双方のメモリに所要
のデータが蓄積されると、サブトラクション処理部１２
においてサブトラクション（差分化）処理が行われる。
すなわち、マスク像の投影データとライブ像の投影デー
タとをそれぞれｌｏｇ（対数）変換した後、ｌｏｇ変換
結果同士を減算処理する。
入しながらＸ線撮影を行う。例えば血管のみを表す画像
の再構成を行う場合は、投影データの収集を造影剤の注
入前と注入後の２回行い、それぞれを第１（マスク像）
の投影データ用メモリ１１１、第２（ライブ像）の投影
データ用メモリ１１２に記憶する。双方のメモリに所要
のデータが蓄積されると、サブトラクション処理部１２
においてサブトラクション（差分化）処理が行われる。
すなわち、マスク像の投影データとライブ像の投影デー
タとをそれぞれｌｏｇ（対数）変換した後、ｌｏｇ変換
結果同士を減算処理する。
【００３６】サブトラクション処理結果は再構成部７に
送られる。再構成部７における画像の再構成処理は第１
実施形態のものとほぼ同様であるが若干異なっている。
すなわち、サブトラクション処理部１２において既にｌ
ｏｇ変換が施されているので、本実施形態の再構成部７
はｌｏｇ変換部７２を備えない。
送られる。再構成部７における画像の再構成処理は第１
実施形態のものとほぼ同様であるが若干異なっている。
すなわち、サブトラクション処理部１２において既にｌ
ｏｇ変換が施されているので、本実施形態の再構成部７
はｌｏｇ変換部７２を備えない。
【００３７】そして本実施形態は、第１実施形態と同様
に、離散化された再構成領域の再構成演算を実施する際
に、投影データの収集角度間隔の不均一性により生じる
再構成画像の劣化を抑制する抑制手段を備えており、フ
ィルタリング部７３により投影データに対しフィルター
補正逆投影法に係る補正フィルタを用いてコンボリュー
ションし、係数補正部７４により、コンボリューション
結果に対し投影データの収集角度間隔の不均一性に伴う
パワーの差異を補正するための補正係数を掛ける。逆投
影部７５は補正係数が掛け合わされた投影データを逆投
影して画像を再構成する。
に、離散化された再構成領域の再構成演算を実施する際
に、投影データの収集角度間隔の不均一性により生じる
再構成画像の劣化を抑制する抑制手段を備えており、フ
ィルタリング部７３により投影データに対しフィルター
補正逆投影法に係る補正フィルタを用いてコンボリュー
ションし、係数補正部７４により、コンボリューション
結果に対し投影データの収集角度間隔の不均一性に伴う
パワーの差異を補正するための補正係数を掛ける。逆投
影部７５は補正係数が掛け合わされた投影データを逆投
影して画像を再構成する。
【００３８】したがって第２実施形態によれば、血管の
みが抽出された画像（サブトラクション像）を得ること
ができるとともに、パワーの差異に起因するアーチファ
クトを低減することができ、再構成画像の劣化を防止で
きる。
みが抽出された画像（サブトラクション像）を得ること
ができるとともに、パワーの差異に起因するアーチファ
クトを低減することができ、再構成画像の劣化を防止で
きる。
【００３９】（第３実施形態）次に本発明の第３実施形
態を説明する。上記第１及び第２実施形態では、投影角
度間隔の不均一性に起因する再構成画像の劣化を抑制す
るために、コンボリューション結果に対し投影角度間隔
の差異に伴うパワーの変化を補正するための補正係数を
掛けて逆投影するという再構成処理の構成について説明
した。これに対し第３実施形態の画像再構成装置は、投
影角度間隔が均等な投影データを補間演算によって求
め、この補間により得られた投影データをコンボリュー
ションして逆投影する。その他の構成については、第１
実施形態のものと同様でありその詳細な説明は省略す
る。
態を説明する。上記第１及び第２実施形態では、投影角
度間隔の不均一性に起因する再構成画像の劣化を抑制す
るために、コンボリューション結果に対し投影角度間隔
の差異に伴うパワーの変化を補正するための補正係数を
掛けて逆投影するという再構成処理の構成について説明
した。これに対し第３実施形態の画像再構成装置は、投
影角度間隔が均等な投影データを補間演算によって求
め、この補間により得られた投影データをコンボリュー
ションして逆投影する。その他の構成については、第１
実施形態のものと同様でありその詳細な説明は省略す
る。
【００４０】図６は本発明の第３実施形態の概略構成を
示すブロック図である。図６（ａ）に示すように、投影
データ２０と再構成部２１との間に、補間処理部２０が
設けられている点において第１実施形態のものとは構成
が異なっている。補間処理部２０の具体的な構成及び動
作については後述する。なお、本実施形態の再構成部２
１は、図６（ｂ）に示すように、第１実施形態に係る再
構成部７の係数補正部７４を具備していない。
示すブロック図である。図６（ａ）に示すように、投影
データ２０と再構成部２１との間に、補間処理部２０が
設けられている点において第１実施形態のものとは構成
が異なっている。補間処理部２０の具体的な構成及び動
作については後述する。なお、本実施形態の再構成部２
１は、図６（ｂ）に示すように、第１実施形態に係る再
構成部７の係数補正部７４を具備していない。
【００４１】フィルター補正逆投影法において用いられ
る逆投影フィルターは、高周波強調フィルターであっ
て、ディラックのδ関数を再構成する場合を考えると、
当該δ関数が存在する位置以外においては、高周波強調
フィルターのあらゆる方向からの逆投影によって正負が
打ち消されて零になる。例えばある方向からは正の値が
逆投影され、また別の方向からは負の値が逆投影される
ことになる。
る逆投影フィルターは、高周波強調フィルターであっ
て、ディラックのδ関数を再構成する場合を考えると、
当該δ関数が存在する位置以外においては、高周波強調
フィルターのあらゆる方向からの逆投影によって正負が
打ち消されて零になる。例えばある方向からは正の値が
逆投影され、また別の方向からは負の値が逆投影される
ことになる。
【００４２】しかしながら、不等角度間隔の投影データ
に対する逆投影においては、逆投影時に加算される正負
のバランスが崩れてしまい、その結果、アーチファクト
を生じてしまう。すなわち逆投影フィルターが掛かるバ
ランスが崩れている。
に対する逆投影においては、逆投影時に加算される正負
のバランスが崩れてしまい、その結果、アーチファクト
を生じてしまう。すなわち逆投影フィルターが掛かるバ
ランスが崩れている。
【００４３】そこで本実施形態においては、逆投影フィ
ルターが掛かるバランスの適切化を図るために、不等角
度間隔で収集された投影データから等角度間隔の投影デ
ータを補間演算により作成する。
ルターが掛かるバランスの適切化を図るために、不等角
度間隔で収集された投影データから等角度間隔の投影デ
ータを補間演算により作成する。
【００４４】特に本実施形態においては、補間処理部２
０が１次補間を実施する場合について説明する。ここ
で、θ＝θ0の時点においてＮフレーム目のデータが収
集され、θ＝θ1の時点において（Ｎ＋１）フレーム目
のデータが収集されているとする。そして現在、θ＝θ
2（θ0≦θ2≦θ1）の時点における投影データを補間に
よって作成する場合、θ＝θ2の時点における投影デー
タｐ（ｘ，ｙ，θ2）は次式（３）に従って算出され
る。
０が１次補間を実施する場合について説明する。ここ
で、θ＝θ0の時点においてＮフレーム目のデータが収
集され、θ＝θ1の時点において（Ｎ＋１）フレーム目
のデータが収集されているとする。そして現在、θ＝θ
2（θ0≦θ2≦θ1）の時点における投影データを補間に
よって作成する場合、θ＝θ2の時点における投影デー
タｐ（ｘ，ｙ，θ2）は次式（３）に従って算出され
る。
【００４５】 ｐ（ｘ，ｙ，θ2）＝ｍ×ｐ（ｘ，ｙ，θ0）＋ｎ×ｐ（ｘ，ｙ，θ1）… （３） ここで、（ｘ，ｙ）は画像上の任意のピクセルを示す座
標系を示し、ｍ，ｎは次式（４）および（５）によって
与えられる。
標系を示し、ｍ，ｎは次式（４）および（５）によって
与えられる。
【００４６】 ｍ＝（θ1−θ2）／（θ1−θ0） …（４） ｎ＝（θ2−θ0）／（θ1−θ0） …（５） このような補間により新たに作成された等角度間隔の投
影データは、再構成部２１に送られる。この実施形態で
は、補間関数として一次補間関数を用いた例を紹介した
が、本発明は一次補間のみに限定されるものではなく、
例えば三次のB-Spline関数のような高次関数を用いても
良い。
影データは、再構成部２１に送られる。この実施形態で
は、補間関数として一次補間関数を用いた例を紹介した
が、本発明は一次補間のみに限定されるものではなく、
例えば三次のB-Spline関数のような高次関数を用いても
良い。
【００４７】再構成部２１における再構成処理では、離
散化された再構成領域の再構成を行う。ここでの再構成
処理は、当該処理が第１実施形態の係数補正処理を行わ
ない点を除いては、第１実施形態の再構成処理と同様で
ある。
散化された再構成領域の再構成を行う。ここでの再構成
処理は、当該処理が第１実施形態の係数補正処理を行わ
ない点を除いては、第１実施形態の再構成処理と同様で
ある。
【００４８】以上述べたように第３実施形態は、不等角
度間隔で収集された投影データから等角度間隔の投影デ
ータを補間演算により作成するものであり、逆投影フィ
ルターが掛かるバランスの適切化を図ることができる。
このため投影角度間隔の不均一性を再構成演算前に改善
できる。したがって、第１実施形態と同様に短時間で、
且つ効果的に再構成画像の劣化を抑制できる。
度間隔で収集された投影データから等角度間隔の投影デ
ータを補間演算により作成するものであり、逆投影フィ
ルターが掛かるバランスの適切化を図ることができる。
このため投影角度間隔の不均一性を再構成演算前に改善
できる。したがって、第１実施形態と同様に短時間で、
且つ効果的に再構成画像の劣化を抑制できる。
【００４９】（第４実施形態）次に本発明の第４実施形
態を説明する。第４実施形態は上述した第３実施形態の
変形に係り、上述したＤＳＡを実施するものである。図
７は本実施形態の構成を示すブロック図である。
態を説明する。第４実施形態は上述した第３実施形態の
変形に係り、上述したＤＳＡを実施するものである。図
７は本実施形態の構成を示すブロック図である。
【００５０】本実施形態では、被検体に対し造影剤を注
入しながらＸ線撮影を行う。例えば血管のみを表す画像
の再構成を行う場合は、投影データの収集を造影剤の注
入前と注入後の２回行い、それぞれを第１（マスク像）
の投影データ用メモリ１１１、第２（ライブ像）の投影
データ用メモリ１１２に記憶する。双方のメモリに所要
のデータが蓄積されると、サブトラクション処理部１２
においてサブトラクション（差分化）処理が行われる。
すなわち、マスク像の投影データとライブ像の投影デー
タとをそれぞれｌｏｇ（対数）変換した後、ｌｏｇ変換
結果同士を減算処理する。
入しながらＸ線撮影を行う。例えば血管のみを表す画像
の再構成を行う場合は、投影データの収集を造影剤の注
入前と注入後の２回行い、それぞれを第１（マスク像）
の投影データ用メモリ１１１、第２（ライブ像）の投影
データ用メモリ１１２に記憶する。双方のメモリに所要
のデータが蓄積されると、サブトラクション処理部１２
においてサブトラクション（差分化）処理が行われる。
すなわち、マスク像の投影データとライブ像の投影デー
タとをそれぞれｌｏｇ（対数）変換した後、ｌｏｇ変換
結果同士を減算処理する。
【００５１】サブトラクション処理結果は補間処理部２
０に送られる。補間処理部２０における補間処理の内容
及び、再構成部２１における画像の再構成処理は第３実
施形態のものとほぼ同様であるが若干異なっている。す
なわち、サブトラクション処理部１２において既にｌｏ
ｇ変換が施されているので、本実施形態の再構成部７は
第３実施形態のｌｏｇ変換部７２を備えない。
０に送られる。補間処理部２０における補間処理の内容
及び、再構成部２１における画像の再構成処理は第３実
施形態のものとほぼ同様であるが若干異なっている。す
なわち、サブトラクション処理部１２において既にｌｏ
ｇ変換が施されているので、本実施形態の再構成部７は
第３実施形態のｌｏｇ変換部７２を備えない。
【００５２】そして補間処理部２０は、第３実施形態と
同様に、逆投影フィルターが掛かるバランスの適切化を
図るために、不等角度間隔で収集された投影データから
等角度間隔の投影データを補間演算により作成する。こ
のため投影角度間隔の不均一性を再構成演算前に改善で
きる。したがって、短時間で、且つ効果的に再構成画像
の劣化を抑制してサブトラクション像を得ることができ
る。
同様に、逆投影フィルターが掛かるバランスの適切化を
図るために、不等角度間隔で収集された投影データから
等角度間隔の投影データを補間演算により作成する。こ
のため投影角度間隔の不均一性を再構成演算前に改善で
きる。したがって、短時間で、且つ効果的に再構成画像
の劣化を抑制してサブトラクション像を得ることができ
る。
【００５３】（第５実施形態）上記した第１乃至第４の
実施形態では、現在開発段階にある三次元ＣＴ装置にお
いて、撮影系として回転アーム（Ｃアームなど）からな
る循環器用Ｘ線撮影装置を用いてデータ収集を行った
後、画像再構成する場合について説明した。一般に普及
している二次元ＣＴ（被検体の断面毎に撮影を行う）
や、被検体の周囲でＸ線管球を回転させながら寝台を移
動することにより、図８のように被検体を螺旋状にスキ
ャンするヘリカルスキャンＣＴ（スパイラルスキャンＣ
Ｔ、螺旋スキャンＣＴとも呼ばれる）等に対しても本発
明は実施可能である。
実施形態では、現在開発段階にある三次元ＣＴ装置にお
いて、撮影系として回転アーム（Ｃアームなど）からな
る循環器用Ｘ線撮影装置を用いてデータ収集を行った
後、画像再構成する場合について説明した。一般に普及
している二次元ＣＴ（被検体の断面毎に撮影を行う）
や、被検体の周囲でＸ線管球を回転させながら寝台を移
動することにより、図８のように被検体を螺旋状にスキ
ャンするヘリカルスキャンＣＴ（スパイラルスキャンＣ
Ｔ、螺旋スキャンＣＴとも呼ばれる）等に対しても本発
明は実施可能である。
【００５４】前者、すなわちヘリカルスキャンを行わな
いＣＴの場合、再構成領域は、全投影方向のＸ線束に内
接する円として定義され、第１実施形態と同様に断面内
における離散化を行なう。また後者の場合も、寝台の天
板の移動方向に垂直な断面内での再構成空間は二次元Ｃ
Ｔと同様に定義され、離散化も同様に行われる。ヘリカ
ルスキャンＣＴは、基本的には任意のスライス位置にお
ける断面像を再構成することができる。
いＣＴの場合、再構成領域は、全投影方向のＸ線束に内
接する円として定義され、第１実施形態と同様に断面内
における離散化を行なう。また後者の場合も、寝台の天
板の移動方向に垂直な断面内での再構成空間は二次元Ｃ
Ｔと同様に定義され、離散化も同様に行われる。ヘリカ
ルスキャンＣＴは、基本的には任意のスライス位置にお
ける断面像を再構成することができる。
【００５５】そして何れの場合においても第１乃至第４
実施形態と同様に投影角度間隔の不均一性によって生じ
る再構成画像の劣化を抑制できる。再構成方法の具体例
としては、ヘリカルスキャンＣＴでは参考文献（東木
祐介他著、「ヘリカルスキャンに於ける補間再構成法の
検討」、Med.Imag.Tech.,Vol.8, No.3, 253-254(199
0)）に示されている対向ビーム補間法や３６０度補間法
により投影データを補間し、ファンビーム再構成法と同
様の方法で再構成を行えば良い。
実施形態と同様に投影角度間隔の不均一性によって生じ
る再構成画像の劣化を抑制できる。再構成方法の具体例
としては、ヘリカルスキャンＣＴでは参考文献（東木
祐介他著、「ヘリカルスキャンに於ける補間再構成法の
検討」、Med.Imag.Tech.,Vol.8, No.3, 253-254(199
0)）に示されている対向ビーム補間法や３６０度補間法
により投影データを補間し、ファンビーム再構成法と同
様の方法で再構成を行えば良い。
【００５６】なお、本発明は上述した実施形態に限定さ
れず、種々変形して実施可能である。例えば、本発明は
図１に示した循環器用Ｘ線診断システムのみに限定され
るものではなく、被検体周辺を回転しながらＸ線像の収
集を行うことができる種々の循環器用Ｘ線診断システム
に適用可能である。
れず、種々変形して実施可能である。例えば、本発明は
図１に示した循環器用Ｘ線診断システムのみに限定され
るものではなく、被検体周辺を回転しながらＸ線像の収
集を行うことができる種々の循環器用Ｘ線診断システム
に適用可能である。
【００５７】
【発明の効果】本発明によれば、投影データ収集時にお
ける角度ピッチの不均一性を起因として発生する偽像
（アーチファクト）による再構成画像の劣化を、効果的
に、且つ処理時間を大幅に増やすことなく抑制し得る画
像再構成装置を提供できる。
ける角度ピッチの不均一性を起因として発生する偽像
（アーチファクト）による再構成画像の劣化を、効果的
に、且つ処理時間を大幅に増やすことなく抑制し得る画
像再構成装置を提供できる。
【００５８】したがって、実際の再構成装置あるいは実
用段階の実験装置等で使用されているフィルター補正逆
投影法（Feldkamp法のような三次元再構成法を含む）に
おいて発生するアーチファクトを、演算時間をほとんど
増やすことなく抑止できる。
用段階の実験装置等で使用されているフィルター補正逆
投影法（Feldkamp法のような三次元再構成法を含む）に
おいて発生するアーチファクトを、演算時間をほとんど
増やすことなく抑止できる。
【図１】本発明の画像再構成装置の第１実施形態に係る
循環器用Ｘ線撮影システムの回転アームの外観を示す
図。
循環器用Ｘ線撮影システムの回転アームの外観を示す
図。
【図２】画像の再構成部を含む上記第１実施形態の構成
を示すブロック図。
を示すブロック図。
【図３】上記第１実施形態に係り、投影データの収集角
度間隔が変化している様子を示す図。
度間隔が変化している様子を示す図。
【図４】上記第１実施形態に係り、投影データのパワー
の差異を説明するための模式図。
の差異を説明するための模式図。
【図５】本発明の画像再構成装置の第２実施形態に係
り、ＤＳＡを実施する場合の構成を示すブロック図。
り、ＤＳＡを実施する場合の構成を示すブロック図。
【図６】本発明の画像再構成装置の第３実施形態の構成
を示すブロック図。
を示すブロック図。
【図７】本発明の画像再構成装置の第４実施形態に係
り、ＤＳＡを実施する場合の構成を示すブロック図。
り、ＤＳＡを実施する場合の構成を示すブロック図。
【図８】本発明の第５実施形態に係るヘリカルスキャン
方式を説明するための図。
方式を説明するための図。
【符号の説明】 １…Ｘ線管球 ２…二次元検出器 ３…回転アーム ４…データ収集部 ５…Ａ／Ｄ変換器 ６…投影データ用メモリ ７…再構成部 ８…画像メモリ ９…Ｄ／Ａ変換器 １０…モニタ ７１…歪補正部 ７２…対数（ｌｏｇ）変換部 ７３…フィルタリング部 ７４…係数補正部 ７５…逆投影部
Claims (4)
- 【請求項１】 被検体に対し多方向からＸ線を照射する
ことによって所定の投影角度間隔における複数の投影デ
ータを収集する収集手段と、 前記被検体内部のＸ線吸収係数分布を表す画像を、前記
投影データに基づいて再構成する再構成手段と、 前記投影角度間隔の不均一性により生じる前記再構成画
像の劣化を抑制する抑制手段と、 を具備することを特徴とする画像再構成装置。 - 【請求項２】 前記抑制手段は、前記再構成手段におい
て行われる逆投影演算に際し、前記投影角度間隔の不均
一性に応じたパワーの差異を補正するための補正係数を
各々の逆投影データに掛ける手段を具備することを特徴
とする請求項１に記載の画像処理装置。 - 【請求項３】 前記補正係数は、基準となる投影角度間
隔に対する所定の投影角度における投影角度間隔の比に
基づいて算出されることを特徴とする請求項２に記載の
画像再構成装置。 - 【請求項４】 前記抑制手段は、前記収集手段により収
集された不等角度間隔の投影データを補間処理すること
により等角度間隔の投影データを作成する補間手段を具
備し、 前記再構成手段は、前記補間手段により作成された投影
データに基づいて前記画像を再構成することを特徴とす
る請求項１に記載の画像再構成装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP27586497A JP4047420B2 (ja) | 1997-10-08 | 1997-10-08 | 画像再構成装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP27586497A JP4047420B2 (ja) | 1997-10-08 | 1997-10-08 | 画像再構成装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH11113892A true JPH11113892A (ja) | 1999-04-27 |
| JP4047420B2 JP4047420B2 (ja) | 2008-02-13 |
Family
ID=17561500
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP27586497A Expired - Lifetime JP4047420B2 (ja) | 1997-10-08 | 1997-10-08 | 画像再構成装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4047420B2 (ja) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000322557A (ja) * | 1999-05-12 | 2000-11-24 | Hitachi Medical Corp | 画像処理装置 |
| JP2007175194A (ja) * | 2005-12-27 | 2007-07-12 | Natl Inst Of Radiological Sciences | Ｃｔ投影データの差分画像再構成方法及び装置 |
| JP2009254428A (ja) * | 2008-04-14 | 2009-11-05 | Canon Inc | Ｘ線動画撮影装置 |
| JP2012170736A (ja) * | 2011-02-23 | 2012-09-10 | Toshiba Corp | Ｘ線コンピュータ断層撮影装置 |
| WO2013051402A1 (ja) * | 2011-10-06 | 2013-04-11 | 株式会社 東芝 | Ｘ線診断装置 |
| WO2016002034A1 (ja) * | 2014-07-03 | 2016-01-07 | 株式会社ニコン | Ｘ線装置、画像形成方法、構造物の製造方法、及び構造物製造システム |
| CN113229834A (zh) * | 2021-05-19 | 2021-08-10 | 有方（合肥）医疗科技有限公司 | 锥形束ct系统的重建图像获取方法、系统及存储介质 |
-
1997
- 1997-10-08 JP JP27586497A patent/JP4047420B2/ja not_active Expired - Lifetime
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000322557A (ja) * | 1999-05-12 | 2000-11-24 | Hitachi Medical Corp | 画像処理装置 |
| JP2007175194A (ja) * | 2005-12-27 | 2007-07-12 | Natl Inst Of Radiological Sciences | Ｃｔ投影データの差分画像再構成方法及び装置 |
| JP2009254428A (ja) * | 2008-04-14 | 2009-11-05 | Canon Inc | Ｘ線動画撮影装置 |
| US8842807B2 (en) | 2008-04-14 | 2014-09-23 | Canon Kabushiki Kaisha | X-ray moving image radiographing apparatus |
| JP2012170736A (ja) * | 2011-02-23 | 2012-09-10 | Toshiba Corp | Ｘ線コンピュータ断層撮影装置 |
| WO2013051402A1 (ja) * | 2011-10-06 | 2013-04-11 | 株式会社 東芝 | Ｘ線診断装置 |
| JP2013081527A (ja) * | 2011-10-06 | 2013-05-09 | Toshiba Corp | Ｘ線診断装置 |
| US9084543B2 (en) | 2011-10-06 | 2015-07-21 | Kabushiki Kaisha Toshiba | X-ray diagnostic apparatus |
| WO2016002034A1 (ja) * | 2014-07-03 | 2016-01-07 | 株式会社ニコン | Ｘ線装置、画像形成方法、構造物の製造方法、及び構造物製造システム |
| JPWO2016002034A1 (ja) * | 2014-07-03 | 2017-05-25 | 株式会社ニコン | Ｘ線装置、画像形成方法、構造物の製造方法、及び構造物製造システム |
| CN113229834A (zh) * | 2021-05-19 | 2021-08-10 | 有方（合肥）医疗科技有限公司 | 锥形束ct系统的重建图像获取方法、系统及存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP4047420B2 (ja) | 2008-02-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4384749B2 (ja) | 高減衰性物体のためのアーティファクト補正 | |
| US6292530B1 (en) | Method and apparatus for reconstructing image data acquired by a tomosynthesis x-ray imaging system | |
| JP4579971B2 (ja) | 動き補償方法および装置 | |
| EP1588664B1 (en) | X-Ray imaging apparatus and its control method | |
| JP3894993B2 (ja) | Ｘ線源に供給されるｘ線管電流を変調させるシステム | |
| US5663995A (en) | Systems and methods for reconstructing an image in a CT system performing a cone beam helical scan | |
| JP4441095B2 (ja) | コーン・ビーム・マルチスライス式ｃｔ補正方法及び装置 | |
| US6285732B1 (en) | Methods and apparatus for adaptive interpolation reduced view CT scan | |
| JP2009529977A (ja) | コンピュータ断層撮影データ収集装置及び方法 | |
| JP2002345808A (ja) | Ｃｔスカウト画像処理のための方法及び装置 | |
| US20130251229A1 (en) | System and method for partial scan artifact reduction in myocardial ct perfusion | |
| JP5324883B2 (ja) | Ｃｔ装置および金属形状抽出方法 | |
| US6332013B1 (en) | Methods and apparatus for tilted helical reconstruction multislice CT | |
| EP1372115A2 (en) | Methods and apparatus for reconstructing an image of an object | |
| IL154862A (en) | Methods and apparatus for multi-slice image reconstruction | |
| JPH11226003A (ja) | 断層画像形成システム | |
| JP4509255B2 (ja) | 透視画像作成方法及び装置 | |
| EP0969414B1 (en) | Computerized tomographic multi-frame image reconstruction method and apparatus for helical scanning | |
| US6980681B1 (en) | Methods and apparatus for helical reconstruction for multislice CT scan | |
| JP4047420B2 (ja) | 画像再構成装置 | |
| JP2007159878A (ja) | Ｘ線ｃｔ装置およびそのｘ線ｃｔ画像再構成方法 | |
| US6647084B1 (en) | Method and apparatus for filtering projection data of a helical scan | |
| JPH119582A (ja) | Ｘ線コンピュータ断層撮影装置 | |
| US7734079B2 (en) | Methods and apparatus for image reconstruction | |
| US6999550B2 (en) | Method and apparatus for obtaining data for reconstructing images of an object |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20040928 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070607 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070619 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070820 |
|
| 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: 20071120 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20071122 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101130 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101130 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111130 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111130 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121130 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131130 Year of fee payment: 6 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313114 Free format text: JAPANESE INTERMEDIATE CODE: R313111 Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
| R371 | Transfer withdrawn |
Free format text: JAPANESE INTERMEDIATE CODE: R371 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313117 Free format text: JAPANESE INTERMEDIATE CODE: R313111 Free format text: JAPANESE INTERMEDIATE CODE: R313114 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| EXPY | Cancellation because of completion of term |