JP2004523037A

JP2004523037A - 収集された投影から空間体積の画像再構成をする方法および装置

Info

Publication number: JP2004523037A
Application number: JP2002561778A
Authority: JP
Inventors: シュレーゲル，ヴォルフガング; エーベルト，マティアス; ヘサー，ユルゲン; シェドラー，ボリス
Original assignee: ドイチェスクレブスフォルシュンクスツェントルムスチフトゥングデスエッフェントリヒェンレヒツ
Priority date: 2001-01-29
Filing date: 2002-01-26
Publication date: 2004-07-29
Also published as: WO2002061686A1; DE10290311D2; EP1356432A1; US20050151736A1; US20040114728A1

Abstract

本発明は、収集された投影から空間体積の画像再構成をする方法および装置に関するものであり、再構成ステップの間に、収集された各々の投影ないし収集された投影の各々の区域が、収集された投影のための記憶装置からいったんデータ処理部に送られ、再構成された空間体積画像のボクセルの強度が、再構成ステップの間に、それぞれボクセルに関連する投影ないし区域について更新される。

Description

【技術分野】
【０００１】
本発明は、収集された投影から空間体積の画像再構成をする方法および装置に関するものである。
【背景技術】
【０００２】
投影から三次元データを生成することは、医療分野において、たとえばコンピュータ断層撮影システム（ＣＴ）で具体化されている。その際に記録された投影から画像を構成するには、専用のマルチプロセッサシステムを用いて、再構成される層１つあたり数秒の、是認できる再構成時間を得るようにするるのが通常であり、従来式のＣＴ機器では層ごとに記録が行われるのに対して、新型のＣＴ機器では多重ライン型の検出器によって、最大４つの層を同時に走査することができる。
【０００３】
最近、従来式のＸ線撮影システムを、断層撮影の画像作成に利用することに大きな関心が寄せられている。一方では、円錐ビームジオメトリーの使用、および延長された二次元の検出器の使用により、体積全体を一度に走査することができる。また他方では、特にＣＴをわざわざ調達しなくてもすむので、従来式のＸ線撮影システムに要する費用は安価であり、実現可能な解像度がＣＴよりも高いのが普通である。
【０００４】
その際に収集された投影から画像再構成をするために、通常、逆投影アルゴリズムが用いられる。この場合、各々の投影画像について、投影中に含まれているグレー値情報が、投影ピクセルから放射源までビームに沿って均等に配分され、幾何学的な距離係数で重みづけされる。このような種類の逆投影は、特に、投影データがまず位置依存的に重みづけされてフィルタリングにかけられ、次いで、ゼロで初期設定された体積データセットに逆投影される、フィルタ補正投影との関連で適用されている。さらに逆投影は、各々のステップに投影操作と逆投影操作が含まれる反復方式を含んでいてよい。画像再構成が終了した後、再構成された空間体積画像の記憶装置から体積データを読み出して、用途に応じてすぐに断面図の形で二次元に視覚化したり、あるいは三次元に視覚化することができる。
【０００５】
その際に必要となる時間は著しく長いことが判明しており、たとえば２５６３のサイズのデータボリュームならば、Ｆｅｌｄｋａｍｐアルゴリズムを用いておよそ１５分で再構成することができる。反復式の再構成法では、典型的には数十回のステップのオーダーである所要の反復回数の分だけ、この時間が長くなる。しかしながら、このような計算時間はリアルタイム再構成とは程遠いものであり、特に、相応の装置を位置管理装置や位置監視装置と連動させて利用しようとする場合には、是認できるものでない。
【０００６】
必要な計算時間を短縮しようとする取組が進められている。たとえばＴ．Ｂｏｒｔｆｅｌｄ著”ｏｐｔｉｍｉｚｅｄｐｌａｎｎｉｎｇｕｓｉｎｇｐｈｙｓｉｃａｌｏｂｊｅｃｔｉｖｅｓａｎｄｃｏｎｓｔｒａｉｎｓ”（物理的な目標物と拘束を用いた最適化プランニング），Ｓｅｍ．ＩｎＲａｄ．Ｏｎｃ．，９：２０−３４（１９９９）では、平行ビームジオメトリーと共通の回転軸の場合、フーリエ空間における投影は体積のフーリエ変換の平面になるので、フーリエ空間での新たな走査と逆変換とによって逆変換を得られるというラドン変換の特性を活用することが提案されている（フーリエ断面定理）。しかしこの方法は、前述したジオメトリーだけに限定されている。さらに、アルゴリズムを並行化することによって、計算速度を高める方法がある。ＳａｓａｋｉＴ，ＦｕｋｕｄａＹ著”ｒｅｃｏｎｓｔｒｕｃｔｉｏｎｏｆ３−ＤＸ−ｒａｙｃｏｍｐｕｔｅｒｉｚｅｄｔｏｍｏｇｒａｐｈｙｉｍａｇｅｓｕｓｉｎｇａｄｉｓｔｒｉｂｕｔｅｄｍｅｍｏｒｙｍｕｌｔｉｐｒｏｃｅｓｓｏｒｓｙｓｔｅｍ”（分散メモリマルチプロセッサシステムを用いたコンピュータ断層撮影三次元Ｘ線画像の再構成），ｔｒａｎｓａｃｔｉｏｎｏｆｔｈｅｉｎｆｏｒｍａｔｉｏｎｐｒｏｃｅｓｓｉｎｇｓｏｃｉｅｔｙｏｆＪａｐａｎ，ｖｏｌ．３８，ｎｏ．９，Ｓｅｐｔ．１９９７，ｐ．１６８１−１６９３およびＤｅｈｎｅｒＧ，ＨｅｒｂｅｒｔＭ著”ｖｅｃｔｏｒｃｏｍｐｕｔｉｎｇｉｎＣＴｉｍａｇｅｒｅｃｏｎｓｔｒｕｃｔｉｏｎ−ａｌｇｏｒｉｔｈｍ，ｄａｔａｒａｔｅ，ｃｏｍｐｕｔｅｐｏｗｅｒ，ｐａｒａｌｌｅｌｐｒｏｃｅｓｓｉｎｇ”（ＣＴ画像再構成アルゴリズムにおけるベクトル計算、データ転送速度、計算能力、並行処理），ＳＰＥＥＤＵＰ，ｖｏｌ．９，ｎｏ．２，Ｄｅｃ．１９９５，ｐ．４１−４７（スイス国）によれば、このような種類の並行化が問題なく可能である。しかしながら、このような種類の並行化では、残念ながらパフォーマンスがほぼ直線的に増大することが判明している。この方法は基本的にメモリに制限されるからである。本来の計算オペレーションは、そのために必要なデータの読み取りよりも、高速で進行するのが普通である。プロセッサの計算速度は早くなる一方であるが、メモリ帯域幅（ライン１本ごとに読み出し可能なデータの数）は徐々にしか増えていないので、この比率がいっそう不都合になっている。その意味で、理論的には並行化によって任意のパフォーマンス向上が可能ではあるが、このことは、ハードウェアが極端に高価になるという代償を払わなくては実現不可能である。別案として、再構成のために専用のハードウェアが利用される。たとえばＡｊａｋｕｉｊａｌａＪ，ＪａｓｋｅＵＭ，ＳａｌｌｉｎｅｎＳ，ＨｅｈｍｉｎｅｎＨ，ＬａｉｔｉｎｅｎＪ著”ｒｅｃｏｎｓｔｒｕｃｔｉｏｎｏｆｄｉｇｉｔａｌｒａｄｉｏｇｒａｐｈｓｂｙｔｅｘｔｕｒｅｍａｐｐｉｎｇ，ｒａｙｃａｓｔｉｎｇａｎｄｓｐｌａｔｔｉｎｇ”（テクスチャマッピング、光線キャスティング、およびスプラッティングにおけるデジタル放射線写真の再構成），ｐｒｏｃｅｅｄｉｎｇｓｏｆｔｈｅ１８ｔｈＡｎｎｕａｌＩｎｔｅｒｎａｔｉｏｎａｌＣｏｎｆｅｒｅｎｃｅｏｆｔｈｅＩＥＥＥＥｎｇｉｎｅｅｒｉｎｇｉｎＭｅｄｉｃｉｎｅａｎｄＢｉｏｌｏｇｙＳｏｃｉｅｔｙ，”ＢｒｉｄｇｉｎｇＤｉｓｃｉｐｌｉｎｅｓｆｏｒＢｉｏｍｅｄｉｃｉｎｅ” Ｃａｔ．Ｎｏ．９６ＣＨ３６０３６）ＩＥＥＥ，Ｐａｒｔｖｏｌ．２，１９９７，ｐ６４３−６４５ｖｏｌ．２．ＮｅｗＹｏｒｋ，ＮＹ，ＵＳＡには、３Ｄグラフィックカードへのテクスチャマッピングによって（ＯｐｅｎＧＬインプリメンテーション）、ＣＴ再構成を高速化できることが記載されている。しかしながら、このような種類のシステムは限定的にしか並行処理に適しておらず、スケーリングすることができない。さらに、ＴｒｅｓｐＶ，ＳｎｅｌｌＲ，ＧｍｉｔｒｏＡＦ著”ｖｉｄｅｏｇｒａｐｈｉｃｔｏｍｏｇｒａｐｈｙＩＩｒｅｃｏｎｓｔｒｕｃｔｉｏｎｗｉｔｈｆａｎ−ｂｅａｍｐｒｏｊｅｃｔｉｏｎｄａｔａ”（ビデオグラフィック・トモグラフィＩＩ扇形ビーム投影データによる再構成），ＩＥＥＥｔｒａｎｓａｃｔｉｏｎｏｎｍｅｄｉｃａｌｉｍａｇｉｎｇ，ｖｏｌ．１３，ｎｏ．１，Ｍａｒｃｈ１９９４，ｐ．１３７−１４３，ＵＳＡによれば、高い計算コストがかかる部分を光学的に解決する方法があるが、これは現実には普及していない。
【発明の開示】
【発明が解決しようとする課題】
【０００７】
本発明の課題は、最低限のハードウェアコストでできるだけ高い処理能力が達成されるような、収集された投影から空間体積の画像再構成をする方法ないし装置を提供することにある。
【課題を解決するための手段】
【０００８】
この課題は、一方では、再構成ステップの間に、収集された各々の投影が、収集された投影のための記憶装置からいったんデータ処理部に送られ、再構成された空間体積画像のボクセルの強度が、再構成ステップの間に、それぞれ関連する投影のボクセルについて更新される、収集された投影から空間体積の画像再構成をする方法によって解決される。このようにして、ロードされるべき投影の数を最低限に抑え、それにより、もっとも時間コストのかかる記憶装置のアクションの回数を減らし、すなわち投影データの読み込み回数を減らして、方法を実施するための総時間を最低限に減らすことが可能である。
【０００９】
ただし本方法では、たとえばドイツ特許出願公開４２２４５６８Ａ１号明細書に記載されているように、三次元データセットから陰影のある二次元の画像を展開することが目的なのではなく、個々の投影から三次元データセットを求めることを目指すものである。
【００１０】
投影面の互いに相対的な位置によっては、それぞれの投影面全体ではなく、投影面の各区域だけを相応に取り扱うのが好ましい場合がある。その場合、ボクセルの位置が変わったときには、それぞれの投影から別の区域を選択し、それぞれの再構成ステップについて、これらの区域を利用することが必要になる。この場合、これらの区域の位置によっては、投影面の特定の領域がこのようなやり方で何度も取り扱われることは言うまでもないが、本発明によれば、各々の区域はその特別な形状で一度だけしかロードされない。特に、それによって全体的な所要ハードウェアが減る。なぜなら、投影に含まれる関連情報を記憶するのに、そのために絶対に必要な所要メモリだけを準備すればよいからである。特にこのような手法は、このような情報を準備することができる中間記憶装置ないしキャッシュの所要メモリを低減する。
【００１１】
他方、前述の課題は、収集された投影のための記憶装置と、再構成された空間体積画像のための記憶装置とを備え、これらの記憶装置がデータ処理部によって相互に組み合わされている、収集された投影から空間体積の画像再構成をする装置によって解決され、この装置は、データ処理部が少なくとも２つの処理パイプラインを含んでおり、これらの処理パイプラインは、一方では、投影または投影の区域のための投影キャッシュの少なくとも１つの記憶領域とそれれぞれ接続されるとともに、他方では、ボクセル中間記憶装置の少なくとも１つの記憶領域と接続されており、ボクセル中間記憶装置は、再構成された空間体積画像のための記憶装置と組み合わされており、投影キャッシュは、収集された投影のための記憶装置と組み合わされていることを特徴としている。
【００１２】
このような種類の装置で、たとえば前述した方法を実施可能であり、このときの処理速度は従来のシステムの速度を明らかに上回っている。他方、このような利点を維持しながら、このような種類の設備で別の方法を実施することも考えられる。
【００１３】
このように、複数の投影、ないしこれらの投影の複数の区域が並行して処理される利点があるので、これらの投影または区域を１回しかロードしなくてよく、特に、このような種類の構成では、いったんロードされた投影ないし投影の区域について、関連するすべてのボクセルが再構成の設定事項に準じて処理される。
【００１４】
特に、再構成のサブステップで、ボクセル空間に由来するサブ立方体を処理することができ、効率の理由およびこれに伴う速度の理由から、このような種類のサブ立方体はボクセル中間記憶装置に保存される。サブ立方体を中間記憶装置に保存することで、この強度データを評価のために中間記憶装置で利用することができ、再構成された空間体積画像のための本来の記憶装置から、そのつどのデータを取り出してくる必要がない。このようなデータの取り出しは、そのサイズが大きいので非常に時間コストがかかるのが普通である。このように実行される方法では、中間記憶装置を満たすために別個のアルゴリズムないし別個のハードウェア構造を利用することができる。
【００１５】
特に、複数の処理パイプライン、ないしは事前設定された並行計算機構造との連係で（あるいはこれらと関わりなく）、適当に選択されたサブ立方体の形態での空間体積画像全体の処理が、その他の方法手順とは関わりなく処理速度を高めるという利点がある。このことは、特に相応のボクセル中間記憶装置がシフトレジスタとして構成されている場合に該当する。
【００１６】
すでに上に述べたように、まず再構成サブステップで、投影面の同じ区域に投影が結像されるボクセルの選択を処理するのが好ましい場合がある。それにより、この区域に存在している投影情報を、比較的小型の記憶装置で、特にキャッシュで利用することができ、キャッシュによって、本発明による他の方法手順に関わりなく、ないしは前述した装置の他の構成要件に関わりなく、処理時間を縮めることができる。なぜなら、このような種類のキャッシュは、これと接続された計算ユニットに対して、はるかに少ないアクセス時間しか有していないからである。
【００１７】
このような区域は、たとえば、それぞれのボクセルを投影面にそれぞれ投影し、それぞれ覆われた面を利用することによって求めることができる。
【００１８】
こうしたボクセルの選択は前述したサブ立方体であるのが好ましく、それにより、前述した解決方法のそれぞれの利点が合算される。当然ながら、サブ立方体は必ずしも立方体の形状を有している必要はなく、むしろ、特定の再構成ステップのシーケンスのために順次必要となる投影ないし区域の数を減らすことができるのであれば、あらゆるボクセル量を適用することができ、特に、単純に連関するあらゆるボクセル量を適用することができる。
【００１９】
ボクセル中間記憶装置の使用は、投影キャッシュないし並行計算機構造と関わりがなくても好ましい。なぜなら、再構成された全部の画像情報を記録するのに必要となるような大型の主記憶装置は、はるかに低速のアクセス速度しか有していないので、このような種類の中間記憶装置によって計算ユニットのアクセス時間を短縮できるからである。
【００２０】
処理時間をいっそう短くするために、第１のボクセル中間記憶装置の代わりに処理パイプラインと接続することができ、第１のボクセル中間記憶装置とは関わりなく、再構成された空間体積画像のための記憶装置とデータを交換することができる、第２のボクセル中間記憶装置が設けられていてよい。その意味で、両方のボクセル中間記憶装置の一方は、再構成された空間体積画像のための記憶装置とデータを交換することができるのに対して、他方のボクセル中間記憶装置は計算操作のために利用される。計算過程が終わると、単純な切換器によってこれらの接続部を交換することができる。このようにして、ボクセル中間記憶装置と、再構成された空間体積画像のための記憶装置との間のデータ交換によるむだ時間が回避される。
【００２１】
さらに本発明は、収集された投影のための記憶装置と、再構成された空間体積画像のための記憶装置とを備え、これらの記憶装置がデータ処理部によって相互に組み合わされている、収集された投影から空間体積の画像再構成をする装置であって、メモリ帯域幅がシステム全体の処理能力よりも下回っているものを提案する。このように構成された装置は、記憶装置にとって可能であるよりも迅速に作業を進めることができ、したがって、ハードウェアが最善に利用されるという結果をもたらす。メモリ帯域幅と処理能力はボクセル／秒を単位として比較するのが好ましく、この場合、再構成された空間体積画像のための記憶装置の能力、および収集された投影のための記憶装置の能力と、処理能力との比較を可能にする別の基準も考えられる。
【００２２】
当然ながら、本発明の関連において「空間体積画像」という用語は、投影に含まれる情報が三次元で復元されて記憶される、あらゆる表現を含むものである。特に、これはボクセル空間内の強度分布であり得る。同様のことは「投影」についても当てはまる。
【発明を実施するための最良の形態】
【００２３】
本発明のその他の特性、目的、および利点は、以下において、本発明の方法ステップならびに計算機アーキテクチャを示す添付図面の説明のなかに記載されている。
【００２４】
図１と図２に模式的に図示しているＸ線撮影設備では、人物１が放射源２によってＸ線透視を受ける。それにより、相応の検出器によって投影３を撮影することができ、この投影は最終的に、相応のビーム円錐４と、Ｘ線透視を受ける人物の身体との間の相互作用を反映するものである。放射源２と相応の検出器は、人物を中心として回転可能に配置されているので、さまざまな投影方向を撮影することができる。当然ながら、このような種類の装置に代えて、逆投影が必要となる別の装置も適用することができる。特に、物品もこのようにして検査することができる。
【００２５】
求められた投影５は、収集された投影のための相応の記憶装置６に保存される。これらの投影を基にして、空間体積画像を求めることが意図され、この空間体積画像は、再構成された空間体積画像のための記憶装置７に保存される（図１から図５を参照）。
【００２６】
両方の記憶装置６および７はデータ処理部８を介して相互に組み合わされており、図３と図４に示す実施例では、データ処理部８は投影キャッシュ９を含んでいる。投影キャッシュ９では、必要に応じて、投影記憶装置６からのデータを保存しておくことができる。投影キャッシュ９は個々の記憶セグメント１０を有しており、伝送された投影面のピクセルデータがこれらのセグメントにそれぞれ保存される。記憶装置単位１０ごとに、シフトレジスタとして構成されたボクセル中間記憶装置１３のセル１２に割り当てられたハードウェアパイプラインないし処理パイプライン１１が設けられている。ハードウェアパイプライン１１は、記憶装置１２に保存されている体積要素について、必要な投影ピクセルを投影記憶装置９から読み出して、この体積要素の寄与を計算する。この寄与が、体積要素のそれまでの寄与に加算される。このとき本実施例では逆投影をするために、「光源−投影ピクセル」のビームに沿って、該当するボクセル全体にわたって強度が均等に配分されるが、具体的な実施形態に応じて、幾何学的な減衰を考慮するさらに別の機能が設けられていてよい。
【００２７】
このとき、各々のボクセルについて、ボクセルの中心がどの点に結像されるかを、各々の投影を前提として計算しなければならない。投影面５における点の強度は、通常、隣接するピクセル強度の双一次補間によって算出される。そして、求められた値が、「光源−ボクセル」の距離の逆２乗と乗算されて、ボクセルにおけるそれまでの寄与に加算される。
【００２８】
すべてのハードウェアパイプライン１１が同期して作動するので、これらは同じ時点で完了する。次いで、シフトレジスタ１３の結果が記憶素子１２をさらにシフトさせ、その結果、各々の体積要素について次の投影面の寄与が規定される。
【００２９】
すべての投影面の寄与が算出されて累積されると、それぞれのデータが再び記憶装置７に戻されて書き込まれる。
【００３０】
それぞれのデータをボクセル中間記憶装置１３から記憶装置７へ戻して新たに保存する時間を無駄に経過させないようにするために、図４に示す実施形態では、第１のボクセル中間記憶装置１３の再保存の進行中にハードウェアパイプライン１１と接続される第２のボクセル中間記憶装置１４が設けられている。このようにして、ハードウェアパイプライン１１のパラレルプロセッサが最善に活用される。中間記憶装置１４に保存されているボクセルが計算されている間に、ボクセル中間記憶装置１３が相応に空になって新規にロードされるので、以後の計算のために利用できるようになり、ならびに、ボクセル記憶装置１４のデータについての計算が終結する。次いで、ボクセル記憶装置１３が再びハードウェアパイプライン１１と接続され、それに対してボクセル記憶装置１４については、記憶装置７との間でデータ交換が行われる。
【００３１】
投影面５のすべてのデータを投影キャッシュ９にそのつど読み込むことも考えられるが、ボクセル中間記憶装置１３のボリュームは限られているので、多数の冗長的なデータがロードされることになる。なぜなら１つの投影ごとに、特定のボクセルについて非常に小さな区域しか関与させることができないからである。その意味で、たとえば図５に符号１５を付しているような区域だけをロードするのが好ましい。本発明との関連で「区域」という用語は、比較的少ない量の、特に単純に連関するピクセルを意味しており、この量の大きさは、そのピクセル量を投影キャッシュ９へ容易にロードできるように選択されている。
【００３２】
原則として、投影キャッシュ９およびボクセル中間記憶装置１３ないし１４には、任意のボクセルないし投影５または区域１５をロードすることができる。しかしながら、投影キャッシュ９に読み込まれる区域１５は相互に相関関係にあるのが好ましい。この相関関係は、再構成されるべき画像空間１６からサブ立方体１７が選択されるように選択することができ、このサブ立方体のボクセルの数は、ボクセル中間記憶装置１３ないし１４の記憶素子１２の数に一致しているのが好ましい。一例として図７には、辺の長さが４つのボクセルであるサブ立方体１７が図示されており、このサブ立方体は、６４の記憶素子を備えるボクセル中間記憶装置１３に保存することができる。
【００３３】
そして投影キャッシュ９には、それぞれ関連する画像データをそれぞれの投影方向と対応するように含んでいる区域１５がすべてロードされ、その様子は一例として図５と図９に示されている。相応の投影データがロードされていれば、それぞれの計算を実行することができ、その際にはシフトレジスタによって、投影キャッシュ１０にロードされている区域１５についてのサブ立方体１７の全ボクセルが並行して計算される。投影キャッシュ９にすべての投影方向が保存されるのが理想的ではあるが、必要な投影５の一部だけを投影キャッシュ９に保存し、その間に相応のデータ交換を意図してから、コストと有益性の勘案を行うこともできる。
【００３４】
このようにして実施される方法手順が図６に示されている。まず、投影面５のどの区域１５が、対応するサブ立方体１７に寄与しているかが判定される。これは、プログラムをインストールするとき、ないしはハードウェア構造を構想するときに、すでに行うことができる。次いで、それぞれサブ立方体について内部のループで投影の全寄与が相応のボクセルについて求められ、このときサブ立方体１７はシフトレジスタ１３によって処理される。次いで別のサブ立方体１７が選択され、その際に考慮された区域はできるだけ維持されて、選択的にいくつかの区域１５を交換するだけでよい。ただし場合によっては、区域１５のセットがそっくり交換される。
【００３５】
本実施例では、投影５の数が投影キャッシュ９の記憶装置１０の数を上回っていることを前提としているので、体積１６のすべてのサブ立方体１７が処理されると、そっくり新しいセットの投影方向ないし投影５が選択され、これらの投影の寄与が同様にして規定される。
【００３６】
具体的な１つの実施態様では、再構成方法が実質的に２つのステップに基づいている。すなわち、最初にデータがフィルタリングされてから、逆投影が実行される。このとき、一般性を制限するものではないが、Ｘ線撮影設備が、体積立方体ないし再構成されるべき画像空間１６の１つの側面の法線と平行な軸を中心として回転したと想定することができる（図１０参照）。そして逆投影ステップでは、この側面と平行に、厚さＮのボクセルスライス１７Ａを切り離して観察し、投影５の対応する投影ライン５Ａ（一例を図１０に符号で示す）を再構成することができる。これらのスライスから、図１１に示すように、ＮｘＮｘＮの大きさのボクセル立方体１７を読み取ることができ、これはボクセルライン１７Ｂに従って反復的に行うのが好ましい。この読み取られる立方体の各々のボクセルに対するＮｘＮｘＮ投影の積み重なり（投影ブロック）の寄与を計算し、ボクセルの最新の値にこれを加算する。
【００３７】
そして投影５を、それぞれボクセルスライス１７Ａに関連するライン５Ａで、個々のボクセルを算出するために利用することができる。特に、それぞれボクセル立方体スライス１７Ａに該当する投影５の部分が、すべての投影５についてロードされている場合には、すべてのボクセル立方体スライス１７Ａへの寄与を算出することができる。したがって各々のパイプライン１１は、これらのデータにとって十分なメモリを含んでいるのが好ましい。体積に対するこのような投影の適切な配置に基づき、これは投影画像における複数の隣接するラインであると考えることができる。この方法がいっそう効率的に進むようにするために、個々のラインが計算中にリロード（２ウェイメモリ）されてアクセスのコンフリクトを回避するか、ないしは複数のメモリバンクを利用するのが好ましい。
【００３８】
メモリのサイズは、この実施形態の場合、およそボクセル立方体幅投影画像の辺の長さ（典型的には５１２または１０２４）の値（典型的には１６ビット）であるのが好ましい。ボクセル立方体幅は典型的には４であり（６４パイプラインは４×４×４に相当しているから）、すなわち、各々のパイプラインについて、４×２バイト×１０２４＝８キロバイトのサイズのバッファが必要であり、つまり合計では５１２キロバイトのバッファメモリとなる。６４のパイプライン段の各々について８つの乗算器が用いられるのが好ましく、すなわち乗算器は合計５１２個となる。
【００３９】
そして、計算手順は次のように進めるのが好ましい：

【００４０】
以上のように、特に、複数の投影の立方体とラインを同時にロードされた状態に保つことができ、シフトレジスタ機構を介して、立方体の各々のボクセルを、ロードされている各々の投影と統合することができる。このことは、算出されたスライスがすべての投影の通過後に完成しており、投影の使用済みラインをほとんどすべて、もはやロードしなくてよいという利点がある。
【００４１】
一例として、図１２と図１３には従来技術に基づく実施例が図示されており、これらの実施例はそれぞれ並列の構造１０８ないし２０８を有しているが、投影面１０５が複数の記憶装置１０６に保存されるか、あるいはボクセル空間２０７が並列メモリに保存されるかのいずれかであり、そのためにコストが著しく増大する。
【００４２】
本発明は、たとえばＣ−Ａｒｍ型血管造影設備や、電子携帯型撮像デバイス（ＥＰＩＤ）と関連する直線加速器で実用化することができる。再構成は、たとえばフィルタ補正逆投影によって、ないしはインペラティブ方式によって行うことができる。フィルタ補正逆投影のためには、投影データがまず位置依存的に重みづけされて、フィルタリングにかけられる。うまく最適化されたフーリエ変換のソフトウェアを用いて周波数空間でフィルタリングを実行すれば、このステップは比較的、時間的にクリティカルでないとみなすことができる。フィルタリングされたプロファイルを逆投影するために、このプロファイルと、ゼロで初期化された体積データセットとが、まず相応の設備カードの記憶装置にロードされる。逆投影は、上に説明したようにして行われる。画像再構成をする反復的な方法では、各々のステップに、投影操作と逆投影操作が含まれる。効率的な投影操作（レイトレーシング）を実現するために、特殊な方法とアーキテクチャがすでに公知である。ボクセルをベースとする逆投影のために、上に説明した方法、ならびに上に説明した装置を用いることができる。
【００４３】
上述したアーキテクチャ、ならびに上述した方法を利用することで、大量のデータを早期に処理して再構成することができるので、上述したシステムを使えば、高エネルギーＣＴを利用して、たとえば直線加速器での患者の位置決めをオンラインで検証することができる。
【図面の簡単な説明】
【００４４】
【図１】Ｘ線撮影設備を示す模式図である。
【図２】図１のＸ線撮影設備を示す断面図である。
【図３】逆投影のための模式的な計算機アーキテクチャである。
【図４】逆投影のための別の計算機アーキテクチャである。
【図５】体積と投影の組み合わせを含む、図３の計算機構造を示す模式図である。
【図６】図３から図５の計算機アーキテクチャによる方法手順である。
【図７】サブ立方体の選択である。
【図８】複数のサブ立方体ないしボクセルについての適当な区域の選択である。
【図９】サブ立方体の複数の投影面である。
【図１０】投影と、再構成されるべき画像空間との考えられる配置である。
【図１１】ボクセルスライスの選択、ボクセルスライスで選択されるボクセル行の選択、ないしボクセル行で選択されるボクセル立方体の選択である。
【図１２】複数の並列の投影記憶装置を備える、従来技術に基づく計算機構造である。
【図１３】強度情報のための複数の並列の記憶装置を備える、従来技術に基づく計算機構造である。

Claims

収集された投影から空間体積の画像再構成をする方法において、再構成ステップの間に、収集された各々の投影ないし収集された投影の各々の区域が、収集された投影のための記憶装置（６）からいったんデータ処理部（８）に送られ、再構成された空間体積画像のボクセルの強度が、再構成ステップの間に、それぞれ関連する投影ないし区域のボクセルについて更新される方法。
複数の投影ないし投影の区域が並行して処理される、請求項１に記載の方法。
複数の投影ないし投影の区域が処理のために投影キャッシュ（９）にロードされる、請求項１または２に記載の方法。
画像再構成のサブステップでサブ立方体（１７）が処理され、このサブ立方体の強度情報がボクセル中間記憶装置（１３）に保存されて、処理中に投影のデータを用いて更新される、請求項１から３までのいずれか１項に記載の方法。
画像再構成のサブステップで、投影面の同じ区域（１５）に投影が結像されるボクセルの選択（１７）が処理される、請求項１から３までのいずれか１項に記載の方法。
収集された投影（５）のための記憶装置（６）と、再構成された空間体積画像のための記憶装置（７）とを備え、これらの記憶装置がデータ処理部（８）によって相互に組み合わされている、収集された投影から空間体積（７）の画像再構成をする装置において、データ処理部（８）が少なくとも２つの処理パイプライン（１１）を含んでおり、これらの処理パイプラインは、一方では、投影キャッシュ（９）の少なくとも１つの記憶領域（１０）とそれれぞれ接続されるとともに、他方では、ボクセル中間記憶装置（１３）の少なくとも１つの記憶領域（１２）とそれぞれ接続されており、ボクセル中間記憶装置（１３）は、再構成された空間体積画像のための記憶装置（７）と組み合わされており、投影キャッシュ（９）は、収集された投影のための記憶装置（６）と組み合わされていることを特徴とする装置。
ボクセル中間記憶装置（１３）がシフトレジスタとして構成されている、請求項６に記載の装置。
第１のボクセル中間記憶装置（１３）と交互に処理パイプライン（１１）と接続され、第１のボクセル中間記憶装置（１３）と関わりなく、再構成された空間体積画像のための記憶装置（７）とデータを交換することができる、少なくとも１つの第２のボクセル中間記憶装置（１４）を備えている、請求項６または７に記載の装置。
収集された投影のための記憶装置（６）と、再構成された空間体積画像のための記憶装置（７）とを備え、これらの記憶装置がデータ処理部（８）によって相互に組み合わされている、収集された投影（５）から空間体積（７）の画像再構成をする装置において、メモリ帯域幅が処理能力よりも下回っていることを特徴とする装置。