JP2952314B2 - サンプル3次元物体の例えばx線による錐形射影での画像の計算及び処理方法と、該計算方法を利用した検査物体の3次元再構成方法 - Google Patents
サンプル3次元物体の例えばx線による錐形射影での画像の計算及び処理方法と、該計算方法を利用した検査物体の3次元再構成方法Info
- Publication number
- JP2952314B2 JP2952314B2 JP2011670A JP1167090A JP2952314B2 JP 2952314 B2 JP2952314 B2 JP 2952314B2 JP 2011670 A JP2011670 A JP 2011670A JP 1167090 A JP1167090 A JP 1167090A JP 2952314 B2 JP2952314 B2 JP 2952314B2
- Authority
- JP
- Japan
- Prior art keywords
- function
- sample
- projection
- gaussian
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims description 42
- 238000004364 calculation method Methods 0.000 title claims description 23
- 238000007689 inspection Methods 0.000 title claims description 8
- 230000006870 function Effects 0.000 claims description 144
- 238000005286 illumination Methods 0.000 claims description 6
- 238000012886 linear function Methods 0.000 claims 1
- 239000013598 vector Substances 0.000 description 17
- 238000001514 detection method Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 6
- 238000005259 measurement Methods 0.000 description 6
- 230000005855 radiation Effects 0.000 description 6
- 229910052704 radon Inorganic materials 0.000 description 5
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 5
- 238000003325 tomography Methods 0.000 description 4
- 238000010521 absorption reaction Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 2
- 238000002583 angiography Methods 0.000 description 2
- 210000004204 blood vessel Anatomy 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 239000002872 contrast media Substances 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 229910052710 silicon Inorganic materials 0.000 description 2
- 239000010703 silicon Substances 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012885 constant function Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000005180 public health Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 230000005428 wave function Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/424—Iterative
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Algebra (AREA)
- Image Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Generation (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Description
【発明の詳細な説明】 産業上の利用分野 本発明は、体積要素でサンプルされた3次元(3D)物
体の例えばX線による錐体射影での画像の計算および処
理方法、ならびにこの計算方法を利用した検査物体の3
次元再構成方法に関する。本発明は、特に、検査すべき
物体が、放射線検査を受ける患者の身体である医療分野
に応用される。このため、本発明は主に3次元再構成法
に関する。しかし、本発明はまた2次元再構成法にも応
用可能である。同様に、本発明はサンプルされた物体の
視覚表示に使用することができる。本発明の目的は、検
査中の3D物体を表すより鮮明な画像の迅速な生成に寄与
することである。
体の例えばX線による錐体射影での画像の計算および処
理方法、ならびにこの計算方法を利用した検査物体の3
次元再構成方法に関する。本発明は、特に、検査すべき
物体が、放射線検査を受ける患者の身体である医療分野
に応用される。このため、本発明は主に3次元再構成法
に関する。しかし、本発明はまた2次元再構成法にも応
用可能である。同様に、本発明はサンプルされた物体の
視覚表示に使用することができる。本発明の目的は、検
査中の3D物体を表すより鮮明な画像の迅速な生成に寄与
することである。
従来の技術 撮像すべき断面で実施された放射線断層密度計測で
(1D)獲得後、物体の断面の画像の2次元(2D)再構成
を実現することはすでに公知である。断層密度計の異な
る世代により、第3世代において扇形ビーム断層密度計
を使用するに到った。この扇形ビーム断層密度計では、
X線点源が、扇形ビームの平面中に並べられた複数のセ
ルを備えたいわゆるマルチディテクタにX線照射する。
検査すべき物体は、X線源とマルチディテクタの間に配
置される。獲得には一続きの照射が行われる。1つの照
射から次の照射へと、X線源とマルチディテクタから構
成される集合体が検査すべき身体の周りを回転する。s
が上記マルチディテクタ上の1セルの縦座標を示し、θ
が照射時点におけるX線源−マルチディテクタ集合体の
位置の角度を示すとすると、P(θ,s)で表される一連
の放射線吸収測定値が得られる。該当する断面におい
て、検査すべき物体の体積要素の座標をx、yとして示
し、物体を通過するX線の線形減衰関数にf(x,y)を
表すとすると、 P(θ,s)=I0exp−∫xf(x,y)dt ・・(1) と表すことができる。この式では、Xは、X線管−マル
チディテクタ集合体が、方向θにあるときセルsに到達
するX線を示す。I0は入射密度である。この式は、線源
が単一エネルギー的である、すなわちtがX線に沿った
一方向の曲線横座標である場合にだけ有効である。
(1D)獲得後、物体の断面の画像の2次元(2D)再構成
を実現することはすでに公知である。断層密度計の異な
る世代により、第3世代において扇形ビーム断層密度計
を使用するに到った。この扇形ビーム断層密度計では、
X線点源が、扇形ビームの平面中に並べられた複数のセ
ルを備えたいわゆるマルチディテクタにX線照射する。
検査すべき物体は、X線源とマルチディテクタの間に配
置される。獲得には一続きの照射が行われる。1つの照
射から次の照射へと、X線源とマルチディテクタから構
成される集合体が検査すべき身体の周りを回転する。s
が上記マルチディテクタ上の1セルの縦座標を示し、θ
が照射時点におけるX線源−マルチディテクタ集合体の
位置の角度を示すとすると、P(θ,s)で表される一連
の放射線吸収測定値が得られる。該当する断面におい
て、検査すべき物体の体積要素の座標をx、yとして示
し、物体を通過するX線の線形減衰関数にf(x,y)を
表すとすると、 P(θ,s)=I0exp−∫xf(x,y)dt ・・(1) と表すことができる。この式では、Xは、X線管−マル
チディテクタ集合体が、方向θにあるときセルsに到達
するX線を示す。I0は入射密度である。この式は、線源
が単一エネルギー的である、すなわちtがX線に沿った
一方向の曲線横座標である場合にだけ有効である。
ラドン理論の適用には、 で表されるP(θ,s)の1次元フーリエ変換を見け出す
必要があった。
必要があった。
で表されるf(x,y)の2次元フーリエ変換を見出すこ
とにより、角度θの直線上で求められる を特定することができる。このことから、 に対する、全て変換 (変数と内挿を変えた後)から得られた2次元フーリエ
逆変換により、測定値P(θ,s)から検査中の断面にお
けるf(x,y)の分布を計算することが可能であること
が導き出された。
とにより、角度θの直線上で求められる を特定することができる。このことから、 に対する、全て変換 (変数と内挿を変えた後)から得られた2次元フーリエ
逆変換により、測定値P(θ,s)から検査中の断面にお
けるf(x,y)の分布を計算することが可能であること
が導き出された。
実際には、過度に長い計算を必要とするフーリエ変換
を避け、後方射影が後に続くいわゆるフィルタリング技
術が利用されている。フィルタリングは、フィルタ関数
q(s)を用いて、測定値を表す関数P(θ,s)のたた
みこみ積分を計算するものである。この積分は次の通り
である。
を避け、後方射影が後に続くいわゆるフィルタリング技
術が利用されている。フィルタリングは、フィルタ関数
q(s)を用いて、測定値を表す関数P(θ,s)のたた
みこみ積分を計算するものである。この積分は次の通り
である。
この式では、tは積分の見掛け変数である。フィルタ
の選択(実質的には、すでに述べた変数および内装の変
更時に生じるのと同じ問題を解決することである)は、
得られた画像の品質を支配する。後方射影操作におい
て、たたみこみ値は該当するX線の通路(θ,s)上に位
置する検査中の物体の断面の体積要素〔ボクセル(voxe
l)〕の全てに割り当てられる。そこから画像を導き出
すために、割り当てられた様々なたたみこみ値が各ボク
セルについて累加される。
の選択(実質的には、すでに述べた変数および内装の変
更時に生じるのと同じ問題を解決することである)は、
得られた画像の品質を支配する。後方射影操作におい
て、たたみこみ値は該当するX線の通路(θ,s)上に位
置する検査中の物体の断面の体積要素〔ボクセル(voxe
l)〕の全てに割り当てられる。そこから画像を導き出
すために、割り当てられた様々なたたみこみ値が各ボク
セルについて累加される。
現在まで、データの獲得および検査された構造の3D再
構成は、検査される身体に沿って、X線源−マルチディ
テクタ集合体を並進運動で移動させ、この身体中の多数
の隣接した2D断面画像を獲得することにより実施されて
きた。しかし、この方法は、影を得るために血管中に造
影剤を注入する血管造影法には全く勧められない。この
ような造影剤の注入は、特に何度も繰り返すと患者に外
傷を負わせる危険性がある。例えば、身体中で256個の
隣接した断面の画像を獲得したいとき、造影剤の注入を
256回行うことは禁止されている。血管造影法における
このような禁止の他にも、2D獲得および対応する2D再構
成の技術は、3D物体を再構成するために使用するにはあ
まりにも時間がかかることは認識すべきである。実際に
は、断面画像中で256点当たり256点の平均解像度では、
現在の断層密度計による獲得時間は約4秒である。そこ
で、256個の断面に必要な獲得は、約30分の検査時間と
なる。この時間の長さは、患者が耐えるのにはあまりに
長すぎ、全体として公衆保健機関(全体コスト)にとっ
ても長すぎる。
構成は、検査される身体に沿って、X線源−マルチディ
テクタ集合体を並進運動で移動させ、この身体中の多数
の隣接した2D断面画像を獲得することにより実施されて
きた。しかし、この方法は、影を得るために血管中に造
影剤を注入する血管造影法には全く勧められない。この
ような造影剤の注入は、特に何度も繰り返すと患者に外
傷を負わせる危険性がある。例えば、身体中で256個の
隣接した断面の画像を獲得したいとき、造影剤の注入を
256回行うことは禁止されている。血管造影法における
このような禁止の他にも、2D獲得および対応する2D再構
成の技術は、3D物体を再構成するために使用するにはあ
まりにも時間がかかることは認識すべきである。実際に
は、断面画像中で256点当たり256点の平均解像度では、
現在の断層密度計による獲得時間は約4秒である。そこ
で、256個の断面に必要な獲得は、約30分の検査時間と
なる。この時間の長さは、患者が耐えるのにはあまりに
長すぎ、全体として公衆保健機関(全体コスト)にとっ
ても長すぎる。
理論では、3D獲得を実施し、検査すべき物体の3D再構
成を直接行うことによりラドン法を一般化することが可
能である。3D再構成とは、検査中の体積体のボクセルを
表すアドレスに位置したメモリセルが、物体中で調べら
れている(放射線学的)現象の分布に対応する情報の項
目を含む数値体積の計算を意味する。しかし、ラドン理
論では、再構成すべき検査中の空間の次元より1次元小
さい次元を有する二次空間アセンブリ、すなわちいわゆ
る超平面(ハイパープレーン)中で撮像すべき物理的特
性の積分に対応する測定値を獲得しなければならない。
言い換えれば、3D空間の場合には、2D超平面について積
分された測定結果を得る必要がある。実際には、点ディ
テクタを用いた放射線獲得は、(1D)直線:X線に関する
特性の積分でしかないのである。実際に線源の軌跡が理
論により定められた条件を満たさない事象においては、
獲得に役立った全ての直線(X線)に沿った射影の認識
に単純に基づいても、再構成に必要な全ての超平面に沿
った射影を計算することはできなくなる。言い換えれ
ば、測定値のラドン空間では点が見失われる。ラドン空
間は、均一には埋められない。その結果、得られた3D再
構成にはアーティファクトが現れることになる。2次元
マルチディテクタが(多数のシリコンフォトディテクタ
を配置し、放射線画像増強スクリーンを使用)が現在考
えられている、この方法による3D再構成は、血管に応用
する場合のように、多数のデータを供給できないときこ
の方法が導く結果の不完全を考慮して、まだ望ましい目
標として留めておくべきである。
成を直接行うことによりラドン法を一般化することが可
能である。3D再構成とは、検査中の体積体のボクセルを
表すアドレスに位置したメモリセルが、物体中で調べら
れている(放射線学的)現象の分布に対応する情報の項
目を含む数値体積の計算を意味する。しかし、ラドン理
論では、再構成すべき検査中の空間の次元より1次元小
さい次元を有する二次空間アセンブリ、すなわちいわゆ
る超平面(ハイパープレーン)中で撮像すべき物理的特
性の積分に対応する測定値を獲得しなければならない。
言い換えれば、3D空間の場合には、2D超平面について積
分された測定結果を得る必要がある。実際には、点ディ
テクタを用いた放射線獲得は、(1D)直線:X線に関する
特性の積分でしかないのである。実際に線源の軌跡が理
論により定められた条件を満たさない事象においては、
獲得に役立った全ての直線(X線)に沿った射影の認識
に単純に基づいても、再構成に必要な全ての超平面に沿
った射影を計算することはできなくなる。言い換えれ
ば、測定値のラドン空間では点が見失われる。ラドン空
間は、均一には埋められない。その結果、得られた3D再
構成にはアーティファクトが現れることになる。2次元
マルチディテクタが(多数のシリコンフォトディテクタ
を配置し、放射線画像増強スクリーンを使用)が現在考
えられている、この方法による3D再構成は、血管に応用
する場合のように、多数のデータを供給できないときこ
の方法が導く結果の不完全を考慮して、まだ望ましい目
標として留めておくべきである。
次のような原理に基づいた根本的に異なる代数的な再
構成が検討されてきた。平面もしくは凹面の2次元マル
チディテクタにより獲得された1組の測定値P(θ,
τ,s)はすでに公知である。また、表現しようとする放
射線吸収現象を表す連続関数f(x,y,z)が存在するこ
とも知られている。再構成によりfを決定する必要があ
る。実際には、全ての計算がデータ処理形式の処理によ
り行われることを考慮に入れて、出力におけるfの認識
はサンプル化された認識である。新しい方法は、で示
され、先験的に設定された離散関数を用いてfを求める
ものである。例えば、は、最初全ボクセルが1(また
はゼロ)に設定された数値体積から成る。次に、この離
散関数は、検査中の体積がこの離散関数に厳密に対応す
るかのように、マルチディテクタ上で射影される。この
ようにして、i(f)として表される推定が得られ
る。これは、 i(f)=Hi ・・・(3) とも表すことができる。
構成が検討されてきた。平面もしくは凹面の2次元マル
チディテクタにより獲得された1組の測定値P(θ,
τ,s)はすでに公知である。また、表現しようとする放
射線吸収現象を表す連続関数f(x,y,z)が存在するこ
とも知られている。再構成によりfを決定する必要があ
る。実際には、全ての計算がデータ処理形式の処理によ
り行われることを考慮に入れて、出力におけるfの認識
はサンプル化された認識である。新しい方法は、で示
され、先験的に設定された離散関数を用いてfを求める
ものである。例えば、は、最初全ボクセルが1(また
はゼロ)に設定された数値体積から成る。次に、この離
散関数は、検査中の体積がこの離散関数に厳密に対応す
るかのように、マルチディテクタ上で射影される。この
ようにして、i(f)として表される推定が得られ
る。これは、 i(f)=Hi ・・・(3) とも表すことができる。
これらの式で、iはマルチディテクタのセル番号iで
あり、Hiは、セルiに対応する射影マトリクスHの行i
を表す。射影マトリクスHは、検査中の身体から独立し
ている。このマトリクスは、射影幾何学にのみ依存して
おり、射影幾何学のモデル化を構成する。次に、
i(f)(推定値)とPi(f)(測定値)を比較する。
同一性がなく、最初から同一性がないことが明らかな場
合には、が訂正される。訂正するためには、再構成す
べき物体上に測定射影と計算射影の間の差の値により後
方射影を行う。これは、k =k-1+λk(Hi *εi k)/|Hi|2 ・・(4) (但し、HiはHi *の随伴演算子であり、εi kは、反復k
−1での測定射影と計算射影の差である)で表される。
あり、Hiは、セルiに対応する射影マトリクスHの行i
を表す。射影マトリクスHは、検査中の身体から独立し
ている。このマトリクスは、射影幾何学にのみ依存して
おり、射影幾何学のモデル化を構成する。次に、
i(f)(推定値)とPi(f)(測定値)を比較する。
同一性がなく、最初から同一性がないことが明らかな場
合には、が訂正される。訂正するためには、再構成す
べき物体上に測定射影と計算射影の間の差の値により後
方射影を行う。これは、k =k-1+λk(Hi *εi k)/|Hi|2 ・・(4) (但し、HiはHi *の随伴演算子であり、εi kは、反復k
−1での測定射影と計算射影の差である)で表される。
このような射影の組は、i(f)とPi(f)の同一
性が充分になるまで繰り返される。この方法は、アカデ
ミックプレス(Academic Press)1980年刊行ハーマン
(G.T.Herman)著「射影からの画像再構成(Images rec
onstruction from projection)」に説明されている。
性が充分になるまで繰り返される。この方法は、アカデ
ミックプレス(Academic Press)1980年刊行ハーマン
(G.T.Herman)著「射影からの画像再構成(Images rec
onstruction from projection)」に説明されている。
さらに、射影マトリクスHが、1および0を持つ2進
マトリクスである場合には、結果は不十分なものであ
る。このように、第1の解決法は、線源Sとマルチディ
テクタ上の各セルの中心とを結ぶ直線Dに沿って取られ
たfの線積分によりを得るというものであった。つま
り、このようにセルの実表面も、該セル上の射影の円錐
特性も考慮にいれられていないので、各セルは無限に寸
法が小さくなると考えられる。この線積分の計算はいく
つかの方法で実施することができる。最も単純な方法
は、X線が横断したボクセルの値の加重和により線積分
を達成することである、これは、 のように表すことができる。
マトリクスである場合には、結果は不十分なものであ
る。このように、第1の解決法は、線源Sとマルチディ
テクタ上の各セルの中心とを結ぶ直線Dに沿って取られ
たfの線積分によりを得るというものであった。つま
り、このようにセルの実表面も、該セル上の射影の円錐
特性も考慮にいれられていないので、各セルは無限に寸
法が小さくなると考えられる。この線積分の計算はいく
つかの方法で実施することができる。最も単純な方法
は、X線が横断したボクセルの値の加重和により線積分
を達成することである、これは、 のように表すことができる。
この式で、fjは、(X線が)透過したボクセル(サン
プルとなった体積要素)の各々における吸収関数を表
し、hijは、このボクセルがセルiに到達するX線によ
り横断された場合に値fj=f(xj,yj,zj)に割り当てら
れる加重を表している。この荷重は、X線とX線が透過
したボクセルとの間の交差の長さから導くことができ
る。この近似法、すなわち引用した最初の文献に記載さ
れいるようないわゆる正方形画素法で得られた結果は、
残念なことに意図していた応用には不十分な品質のもの
であった。
プルとなった体積要素)の各々における吸収関数を表
し、hijは、このボクセルがセルiに到達するX線によ
り横断された場合に値fj=f(xj,yj,zj)に割り当てら
れる加重を表している。この荷重は、X線とX線が透過
したボクセルとの間の交差の長さから導くことができ
る。この近似法、すなわち引用した最初の文献に記載さ
れいるようないわゆる正方形画素法で得られた結果は、
残念なことに意図していた応用には不十分な品質のもの
であった。
アイイーイーイー−エムアイ(IEEE−MI)、第1巻、
第3号の192〜196頁にあるジョセフ(P.M.Joseph)によ
る論文「画素画像を通したX線再照射のための改善され
たアルゴリズム(An improved algorithm for reprojec
ting rays through pixel images)」には、別の加重計
算が考えられ、説明されている。この計算法の基本的な
考えは、X線が続く直線と、マルチディテクタセルの中
心に到達するX線が通ったあるいはこのX線が触れたボ
クセルの中心を通過する水平(または鉛直)直線との交
差点で得られるfの内挿値に関する和によりfの線積分
を達成するというものである。これら交差点におけるf
の値は、X線の直線の勾配に沿った水平方向(もしくは
鉛直方向)の最も接近して位置している2つのサンプル
の間の線形内挿により求められる。
第3号の192〜196頁にあるジョセフ(P.M.Joseph)によ
る論文「画素画像を通したX線再照射のための改善され
たアルゴリズム(An improved algorithm for reprojec
ting rays through pixel images)」には、別の加重計
算が考えられ、説明されている。この計算法の基本的な
考えは、X線が続く直線と、マルチディテクタセルの中
心に到達するX線が通ったあるいはこのX線が触れたボ
クセルの中心を通過する水平(または鉛直)直線との交
差点で得られるfの内挿値に関する和によりfの線積分
を達成するというものである。これら交差点におけるf
の値は、X線の直線の勾配に沿った水平方向(もしくは
鉛直方向)の最も接近して位置している2つのサンプル
の間の線形内挿により求められる。
正方形画素法の一般化が、応用光学(Applied Optic
s)第24巻第23号(1985)4028〜4039頁の「計算断層撮
影のための局所基本−関数(Local basis−functions a
pproach to computed tomography)」において、ハンソ
ン(H.M.Hanson)およびウェクサング(G.W.Wecksung)
により提案されている。この一般化では、連続関数f
(x,y,z)がサンプルf(xj,yj,zj)から構成され、こ
のマルチディテクタのセルi上の連続関数の射影が計算
される。この一般化において、関数fは、一組の基本関
数bj(x,y,z)の線形組合わせとして定義され、 で表すことができる。
s)第24巻第23号(1985)4028〜4039頁の「計算断層撮
影のための局所基本−関数(Local basis−functions a
pproach to computed tomography)」において、ハンソ
ン(H.M.Hanson)およびウェクサング(G.W.Wecksung)
により提案されている。この一般化では、連続関数f
(x,y,z)がサンプルf(xj,yj,zj)から構成され、こ
のマルチディテクタのセルi上の連続関数の射影が計算
される。この一般化において、関数fは、一組の基本関
数bj(x,y,z)の線形組合わせとして定義され、 で表すことができる。
この式では、ajは、最終的な分析で求められる関数f
の既述を表しており、関数bj(x,y,z)は公知である。
局所的では、反復的かつ可分の基本関数が通常利用され
る。「局所的である」とは、各関数bjの支持線が、fの
支持線と比較して小さいことを意味する。「反復的であ
る」とは、各関数bjが単一関数b(x,y,z)からの変換
により演繹されることを意味する。この演繹法は、 bj(x,y,z) =b(x−xj,y−yj,z−zj) ・・(7) で表される。「可分である」とは、 b(x,y,z) =bx(x)・by(y)・bz(z) ・・(8) となるような3つの関数bx(x)、by(y)、bz(Z)
が存在することを意味する。
の既述を表しており、関数bj(x,y,z)は公知である。
局所的では、反復的かつ可分の基本関数が通常利用され
る。「局所的である」とは、各関数bjの支持線が、fの
支持線と比較して小さいことを意味する。「反復的であ
る」とは、各関数bjが単一関数b(x,y,z)からの変換
により演繹されることを意味する。この演繹法は、 bj(x,y,z) =b(x−xj,y−yj,z−zj) ・・(7) で表される。「可分である」とは、 b(x,y,z) =bx(x)・by(y)・bz(z) ・・(8) となるような3つの関数bx(x)、by(y)、bz(Z)
が存在することを意味する。
bxが、byおよびbzに等しく、これら関数が、方形波関
数に等しい事象においては、いわゆる正方形画素方法が
再び出現することに留意されたい。基本関数の方法によ
り得られる一般化は、すでに述べた他の方法への復帰を
可能にするが、これらの方法にさらに明瞭な数学的基礎
を与えると、実施された操作はさらに容易に理解される
と言うことができる。実際には、最良の結果は、B−ス
プライン、主要正弦、ハミング関数、あるいは望ましく
はガウス関数のような基本関数を用いて得られる。この
ような基本関数により得られた結果は、勿論計算時間は
長いが、ジョセフ方法により得られる結果より品質が優
れている。
数に等しい事象においては、いわゆる正方形画素方法が
再び出現することに留意されたい。基本関数の方法によ
り得られる一般化は、すでに述べた他の方法への復帰を
可能にするが、これらの方法にさらに明瞭な数学的基礎
を与えると、実施された操作はさらに容易に理解される
と言うことができる。実際には、最良の結果は、B−ス
プライン、主要正弦、ハミング関数、あるいは望ましく
はガウス関数のような基本関数を用いて得られる。この
ような基本関数により得られた結果は、勿論計算時間は
長いが、ジョセフ方法により得られる結果より品質が優
れている。
fの推定値を射影する操作は、実際には、 を計算することである。
この式では、giは線源sとセルiにより形成される扇
形の外側でゼロ値(または非常に小さい値)を有する加
重関数を示す。関数giは、セルiの照射用の錐形支持線
を表す。これは、 と表すことができるので、の計算は全重みhijの計算
となる。各重みhijはセルi上の物体の射影への基本関
数bjの寄与を表す。
形の外側でゼロ値(または非常に小さい値)を有する加
重関数を示す。関数giは、セルiの照射用の錐形支持線
を表す。これは、 と表すことができるので、の計算は全重みhijの計算
となる。各重みhijはセルi上の物体の射影への基本関
数bjの寄与を表す。
発明が解決しようとする課題 3つの変数の関数fのサンプルfjの認識から出発し
て、本発明の目的は、検出セルiに関するfの射影の数
値計算を可能にすると共に、得られた結果の優れた品質
と、非常に短時間の計算の両方を確実にすることであ
る。これら2つの基準(品質と速度)は、代数的再構成
法による一組の2D錐形射影からの関数f(x,y,z)の再
構成に必須である。従来技術の射影方法では、望ましい
品質と速度を確実にすることができなかったが、本発明
はこれを可能にする。
て、本発明の目的は、検出セルiに関するfの射影の数
値計算を可能にすると共に、得られた結果の優れた品質
と、非常に短時間の計算の両方を確実にすることであ
る。これら2つの基準(品質と速度)は、代数的再構成
法による一組の2D錐形射影からの関数f(x,y,z)の再
構成に必須である。従来技術の射影方法では、望ましい
品質と速度を確実にすることができなかったが、本発明
はこれを可能にする。
課題を解決するための手段 そこで、本発明は、体積要素でサンプルされた3次元
物体の画像の計算および処理方法であって、この画像は
X線による錐形射影で例えば照射され、 ガウス基本関数(bj)の空間についてサンプルとなっ
た3次元物体を分解し、 各射影画像要素ごとに、該当する画像要素に射影され
る体積要素に対応する基本関数の寄与(aj)を計算し、
各寄与が、体積要素に関する基本関数と、該当する画像
要素の照射用の支持線を表す関数(gj)との積の積分に
等しく、 射影された画像が処理される方法において、 支持線を表す関数の代わりにガウス加重関数を用いる
ことを特徴とする方法に関する。
物体の画像の計算および処理方法であって、この画像は
X線による錐形射影で例えば照射され、 ガウス基本関数(bj)の空間についてサンプルとなっ
た3次元物体を分解し、 各射影画像要素ごとに、該当する画像要素に射影され
る体積要素に対応する基本関数の寄与(aj)を計算し、
各寄与が、体積要素に関する基本関数と、該当する画像
要素の照射用の支持線を表す関数(gj)との積の積分に
等しく、 射影された画像が処理される方法において、 支持線を表す関数の代わりにガウス加重関数を用いる
ことを特徴とする方法に関する。
本発明はまた、検査中の物体の3次元再構成方法であ
って、再構成すべき関数fを反復的に見出すために、再
構成すべき物体を先験的に表すサンプル化離散関数から
出発して、反復的に上記計算方法を使用する方法に関す
る。この場合、画像要素は2Dマルチディテクタの1セル
で識別される。
って、再構成すべき関数fを反復的に見出すために、再
構成すべき物体を先験的に表すサンプル化離散関数から
出発して、反復的に上記計算方法を使用する方法に関す
る。この場合、画像要素は2Dマルチディテクタの1セル
で識別される。
実施例 第1図は、マルチディテクタDの検出セルiが、身体
Cを通過後の放射線の減衰を確実に測定できるするよう
な方法で身体Cを照射するX線源Sを示す。以下に行う
説明では、使用されるマルチディテクタは平面(放射線
輝度増強スクリーンの外表面が平面とみなして)である
と仮定する。しかし、本発明は、光検出セルを備えたシ
リコンマルチディテクタと同様に凹状のマルチディテク
タの使用に容易に適用することができる。このタイプの
凹状マルチディテクタは、照射する錐形X線に対してあ
らゆる点において表面が垂直である。
Cを通過後の放射線の減衰を確実に測定できるするよう
な方法で身体Cを照射するX線源Sを示す。以下に行う
説明では、使用されるマルチディテクタは平面(放射線
輝度増強スクリーンの外表面が平面とみなして)である
と仮定する。しかし、本発明は、光検出セルを備えたシ
リコンマルチディテクタと同様に凹状のマルチディテク
タの使用に容易に適用することができる。このタイプの
凹状マルチディテクタは、照射する錐形X線に対してあ
らゆる点において表面が垂直である。
第2図は、3つの直接正規直交基準フレームを示す。
第1正規直交フレームx,y,zは、身体Cと一体である。
線源SとマルチディテクタDから構成される集合体の1
つの位置について、これは線源Sと一体であるとみなす
こともできる。平面と仮定されるディテクタDは、上記
マルチディテクタの平面上の線源Sの正規直交射影であ
るPとして示される中心点により特徴付けられる。それ
に応じて、ディテクタDに属する基準フレームx′,
y′,z′は、z′がSPと同一線上で、かつ同方向にあ
り、x′及びy′(単純化合のため、ディテクタDの直
交する縁に沿って方向付けされている)はz′と共に直
接正規直交基準フレームを形成している。セルiは、デ
ィテクタDの平面上の正方形と仮定され、ziがSPi(Pi
はセルiの中心である)と同一線上にあるような直接正
規直交基準フレームxi,yi,ziに属する。正規直交ベクト
ルxiとyiはziと直交し、ベクトルxiはディテクタDの平
面に含まれている。X線SPiは、ベクトルxiを含む平面
πiに垂直である。
第1正規直交フレームx,y,zは、身体Cと一体である。
線源SとマルチディテクタDから構成される集合体の1
つの位置について、これは線源Sと一体であるとみなす
こともできる。平面と仮定されるディテクタDは、上記
マルチディテクタの平面上の線源Sの正規直交射影であ
るPとして示される中心点により特徴付けられる。それ
に応じて、ディテクタDに属する基準フレームx′,
y′,z′は、z′がSPと同一線上で、かつ同方向にあ
り、x′及びy′(単純化合のため、ディテクタDの直
交する縁に沿って方向付けされている)はz′と共に直
接正規直交基準フレームを形成している。セルiは、デ
ィテクタDの平面上の正方形と仮定され、ziがSPi(Pi
はセルiの中心である)と同一線上にあるような直接正
規直交基準フレームxi,yi,ziに属する。正規直交ベクト
ルxiとyiはziと直交し、ベクトルxiはディテクタDの平
面に含まれている。X線SPiは、ベクトルxiを含む平面
πiに垂直である。
上述した説明に従って、本発明は、基本関数へのfの
分解の一般原理を利用している。選択された基本関数は
局所的、可分、反復的で、ガウス基本関数である。その
結果、 bj(x,y,z)=bx(x−xj)・by(y−yj)・bz(z−z
j) のとき、 と表すことができる。
分解の一般原理を利用している。選択された基本関数は
局所的、可分、反復的で、ガウス基本関数である。その
結果、 bj(x,y,z)=bx(x−xj)・by(y−yj)・bz(z−z
j) のとき、 と表すことができる。
ここで、空間の3つの方向において解像度が同じであ
る格子上にサンプルボクセルが分布されていると仮定す
る。δx、δyおよびδzが、3つの軸x、yおよびz
に沿ったサンプリングピッチを示すとすれば、 δx=δy=δz ・・(12a) と表すことができる。
る格子上にサンプルボクセルが分布されていると仮定す
る。δx、δyおよびδzが、3つの軸x、yおよびz
に沿ったサンプリングピッチを示すとすれば、 δx=δy=δz ・・(12a) と表すことができる。
このとき、一般式: bx(u)=δx・Gσb(u) ・・(12b) を有するガウス関数が選択される。
この式では、Gσb(u)は、 |u|≦Rbのとき、 また、|u|>Rbのとき、 Gσb(u)=0 ・・(13) のような標準偏差σbを有する切頭ガウス関数を示す。
上記式で、切頭ガウス関数の切断の半径の値Rbならび
に標準偏差σbは、ユーザが選択することができる。ガ
ウス関数は、計算がデジタルコンピュータ中で必然的に
制限されるので、切断される。しかし、この切断は、計
算の妥当性に認め得るほどの影響を及ぼさない。実用的
に優れた結果は、中央高さでの幅がσxに等しいガウス
関数、すなわち、標準偏差が、 により与えられるガウス関数を採用することにより得ら
れた。これは、中央高さでの幅とガウス関数の標準偏差
の間に比 が存在するためである。上記式において、nは自然対
数を表す。
に標準偏差σbは、ユーザが選択することができる。ガ
ウス関数は、計算がデジタルコンピュータ中で必然的に
制限されるので、切断される。しかし、この切断は、計
算の妥当性に認め得るほどの影響を及ぼさない。実用的
に優れた結果は、中央高さでの幅がσxに等しいガウス
関数、すなわち、標準偏差が、 により与えられるガウス関数を採用することにより得ら
れた。これは、中央高さでの幅とガウス関数の標準偏差
の間に比 が存在するためである。上記式において、nは自然対
数を表す。
同様に、半径Rbについて満足のゆく実験値はRb=3δ
x/2であった。採用した項δx、δy、δzは、fが一
定関数〔すなわち、f(x,y,z)=a〕である場合に、
係数ajが一定で、aとほとんど同じであることを確実に
する標準化信号である。この標準化により、必要であれ
ば、サンプルfjとガウス関数の係数ajの間の移行を単純
化することが可能になる(近似的にaj=fjを選択するこ
とにより)。しかし、この標準化は本発明において必須
ではない。
x/2であった。採用した項δx、δy、δzは、fが一
定関数〔すなわち、f(x,y,z)=a〕である場合に、
係数ajが一定で、aとほとんど同じであることを確実に
する標準化信号である。この標準化により、必要であれ
ば、サンプルfjとガウス関数の係数ajの間の移行を単純
化することが可能になる(近似的にaj=fjを選択するこ
とにより)。しかし、この標準化は本発明において必須
ではない。
公知のように、加重関数giは、頂点Sと、底面として
のセルiを有する多面体内では負ではなく、この多面体
の外側ではほぼゼロでなければならない。さらに、この
多面体内では、平面zi(軸FPiに対して直角にある平
面)の各々に割り当てられた重みは、これらの面の各々
が、線源から発する同数のX線と交差するので、同じで
なければならない。この制約の結果、セルiに属する基
準フレームxi,yi,zi中で求められる積分: ∬gi(xk,yk,zk)dxkdyk ・・(15) は、zkから独立していなければならない。
のセルiを有する多面体内では負ではなく、この多面体
の外側ではほぼゼロでなければならない。さらに、この
多面体内では、平面zi(軸FPiに対して直角にある平
面)の各々に割り当てられた重みは、これらの面の各々
が、線源から発する同数のX線と交差するので、同じで
なければならない。この制約の結果、セルiに属する基
準フレームxi,yi,zi中で求められる積分: ∬gi(xk,yk,zk)dxkdyk ・・(15) は、zkから独立していなければならない。
これら2つの条件を満たすために、本発明では特殊な
加重関数jiを採用した。これら関数はガウス関数であ
る。これらは切頭ガウス関数であることが望ましい。実
際には、各平面zi=一定において、これら関数は2Dガウ
ス関数である。しかし、3D空間では、これら関数はガウ
ス関数ではない。それでもその数学的公式化にはガウス
関数の使用が必要である。
加重関数jiを採用した。これら関数はガウス関数であ
る。これらは切頭ガウス関数であることが望ましい。実
際には、各平面zi=一定において、これら関数は2Dガウ
ス関数である。しかし、3D空間では、これら関数はガウ
ス関数ではない。それでもその数学的公式化にはガウス
関数の使用が必要である。
第3図は、ディテクタDの平面中での加重関数の1つ
の限定を概略的に示したものである。極く僅かな近似値
以内で、この限定は切頭2Dガウス関数である。このガウ
ス関数の中心はPiであり、この標準偏差はσg、その切
断半径はRgである。標準偏差σgと半径Rgは、ユーザが
選択することができる。ディテクタの平面中の位置を採
用する代わりに、giの数学的表現を与えるためには、セ
ル(xi,yi,zi)の基準フレーム中の位置を採用する方が
好ましいであろう。Oが、中心Pi、半径σgを有するデ
ィテクタ平面中の円を示すとすると、セルiに属し、す
でに定義された平面πi上の該円の方向ziでの射影は楕
円であると仮定することができる。この楕円はxiとyiに
沿った方向の軸を有する。xiがディテクタDの平面およ
び平面πiと共通であることを考慮に入れ、半径SPiが
角度τ(第2図)で半径SPに対して傾斜していることを
測定により確実にすることにより、 と表すことができる。
の限定を概略的に示したものである。極く僅かな近似値
以内で、この限定は切頭2Dガウス関数である。このガウ
ス関数の中心はPiであり、この標準偏差はσg、その切
断半径はRgである。標準偏差σgと半径Rgは、ユーザが
選択することができる。ディテクタの平面中の位置を採
用する代わりに、giの数学的表現を与えるためには、セ
ル(xi,yi,zi)の基準フレーム中の位置を採用する方が
好ましいであろう。Oが、中心Pi、半径σgを有するデ
ィテクタ平面中の円を示すとすると、セルiに属し、す
でに定義された平面πi上の該円の方向ziでの射影は楕
円であると仮定することができる。この楕円はxiとyiに
沿った方向の軸を有する。xiがディテクタDの平面およ
び平面πiと共通であることを考慮に入れ、半径SPiが
角度τ(第2図)で半径SPに対して傾斜していることを
測定により確実にすることにより、 と表すことができる。
次に、giに関して、次の形: gi(xk,yk,zk) =Gσgx(zk)(xk)・Gσgy(ZK)(yk) ・・(17) の可分式を選択することができる。
この式では、Gσ(u)は、 |xk|≦Rgx(zk)のとき、 また、|xk|>Rgx(zk)のとき、 Gσgx(zk)(xk)=0 ・・(18) となるような切頭1Dガウス関数である。
Gσgy(zk)も同様に定義する。これらの式で、Rgx
は、軸xi上の切断の半径である(同様に、Rgyは軸yi上
の切断の半径である)。実際は、ガウス関数はzkにより
パラメータ表示されている。というのは、射影の錐形特
性により、標準偏差が明らかに線源Sに対して加重され
た基本関数の高さに依るものであるからである。このた
め、標準偏差σgx(zk)およびσgy(zk)は、 により定義される。
は、軸xi上の切断の半径である(同様に、Rgyは軸yi上
の切断の半径である)。実際は、ガウス関数はzkにより
パラメータ表示されている。というのは、射影の錐形特
性により、標準偏差が明らかに線源Sに対して加重され
た基本関数の高さに依るものであるからである。このた
め、標準偏差σgx(zk)およびσgy(zk)は、 により定義される。
SPi/zkは、ziに平行な平面πi上の射影の拡大比とし
て示される。
て示される。
FMAG(zk)=SPi/zk ・・・(20) と表すことができる。
同様に、切断半径をFMAGの関数として、 のように定義することができるであろう。
このように定義されたgiは次のような特性を持ってい
る。平面πi中のgiの限定は、定義により、xiに沿った
標準偏差σgxと、yiに沿った標準偏差σgyを有する切頭
2Dガウス関数である。giの支持線は、頂点S、軸SPi、
ならびに楕円形の底面(平面πi中で)を有する錐形で
ある。ディテクタ平面でのgiの限定は、標準偏差σ
g(x′またはy′に沿った)を有する切頭2Dガウス関
数と同じである。この限定の支持線は、中心Piと半径Rg
を持った円とほとんど相違ない。giの制限のための支持
線は、線源Sが実際には無限ではないので、完全な円で
はない。線源Sが無限で、平面πi中での射影が実際に
楕円形であれば、ディテクタ平面での射影は厳密に円と
なる。射影が円錐状であるので、射影された円は若干変
形している。この結果生じたエラーは、二次的なもので
あり、無視することができる。さらに、このように選択
された関数giは、必要条件を満たし、これによって、 ∬gi(xk,yk,zk)dxkdyk ・・(22) が、ziから独立する。実際に、この二重の積は、 ∬Gσgx(zk)(xk)・Gσgy(ZK)(yk)dxkdyk≒
1 ・・(23) のように表すことができ、一定である。使用されたガウ
ス関数は純粋なガウス関数ではなく、切頭関数であるの
で、相等性は厳密ではない。最終的に、基本関数の場合
のように、実際には、ガウス関数の中央高さでの幅が、
ボクセルの幅(属する基準フレーム中の検出セルの位置
で)と等しいような加重関数の標準偏差が選択される。
同様に、この幅の3/2に等しい切断Rgの半径が選択され
る。その結果、射影i(f)の推定のために得られた
式は、次の通りである。
る。平面πi中のgiの限定は、定義により、xiに沿った
標準偏差σgxと、yiに沿った標準偏差σgyを有する切頭
2Dガウス関数である。giの支持線は、頂点S、軸SPi、
ならびに楕円形の底面(平面πi中で)を有する錐形で
ある。ディテクタ平面でのgiの限定は、標準偏差σ
g(x′またはy′に沿った)を有する切頭2Dガウス関
数と同じである。この限定の支持線は、中心Piと半径Rg
を持った円とほとんど相違ない。giの制限のための支持
線は、線源Sが実際には無限ではないので、完全な円で
はない。線源Sが無限で、平面πi中での射影が実際に
楕円形であれば、ディテクタ平面での射影は厳密に円と
なる。射影が円錐状であるので、射影された円は若干変
形している。この結果生じたエラーは、二次的なもので
あり、無視することができる。さらに、このように選択
された関数giは、必要条件を満たし、これによって、 ∬gi(xk,yk,zk)dxkdyk ・・(22) が、ziから独立する。実際に、この二重の積は、 ∬Gσgx(zk)(xk)・Gσgy(ZK)(yk)dxkdyk≒
1 ・・(23) のように表すことができ、一定である。使用されたガウ
ス関数は純粋なガウス関数ではなく、切頭関数であるの
で、相等性は厳密ではない。最終的に、基本関数の場合
のように、実際には、ガウス関数の中央高さでの幅が、
ボクセルの幅(属する基準フレーム中の検出セルの位置
で)と等しいような加重関数の標準偏差が選択される。
同様に、この幅の3/2に等しい切断Rgの半径が選択され
る。その結果、射影i(f)の推定のために得られた
式は、次の通りである。
ただし、 A=Gσgx(zk)(xk)・Gσgy(ZK)(yk) 前記の式で、変数zkを排除できることは明らかであろ
う。これは、近似を犠牲にしても可能である。前述の応
用において十分に確認されたように、検出セルのサイド
δX′は、距離SPまたはSPiと比較して非常に小さいと
仮定する。本発明が今使用されている条件において、例
えば、δX′は0.8mmにほぼ等しく、SPは1200mmにほぼ
等しい。この仮定に関して、giの支持線に対応する円錐
の頂点における角度αは非常に小さい。上記の数値をあ
てはめると、α≒0.04゜となる。
う。これは、近似を犠牲にしても可能である。前述の応
用において十分に確認されたように、検出セルのサイド
δX′は、距離SPまたはSPiと比較して非常に小さいと
仮定する。本発明が今使用されている条件において、例
えば、δX′は0.8mmにほぼ等しく、SPは1200mmにほぼ
等しい。この仮定に関して、giの支持線に対応する円錐
の頂点における角度αは非常に小さい。上記の数値をあ
てはめると、α≒0.04゜となる。
基本関数3Dbjを考慮すると、その支持線がgiの支持円
錐と空交差するか否かがわかる。この交差が空となる場
合には、その射影i(f)への寄与はゼロである。そ
の反対の場合には、交差は、 のような2つの平行な半平面により結合される。
錐と空交差するか否かがわかる。この交差が空となる場
合には、その射影i(f)への寄与はゼロである。そ
の反対の場合には、交差は、 のような2つの平行な半平面により結合される。
この式において、zjは、基本関数bjの中心が置かれる
点Vjの参照セル(xi,yi,zi)中のレベルを示す。第4図
は、点Vjについて限定される関数bjの定義の局所領域を
示す。この領域は、半径Rbによって限定される円である
(Rbは、bjを定義するガウス関数の切断半径である)。
第4図は、式11により与えられる積分を示している。局
所基本関数を考慮に入れる限り、該当する検出セルの加
重関数の支持線の一部分との交差を求める。支持線のこ
の部分は通常楕円錐台であり、その最大の底面は、高さ
zj−Rbにより設定され、またその最小の底面は高さzj−
Rbにより設定される。本発明において、一方でαが非常
に小さく、他方でRbもまた小さい(「局所」底面の仮定
に基づく)ことを考慮に入れると、この円錐台は、軸SP
iを有する円筒要素により、gjの近傍に接近し得ること
は明らかになった。言い換えれば、zkが、間隔(zj -,zj
+)に属するところでは、 σgx(zk)≒σgx(zj) 故に、σgx(zk)≒σg/FMAG(zj) また、σgy(zk)≒σgy(zj) 故に、σgy(zk)≒σgcos(τ)/FMAG(zj) ・・・(26) と表すことができる。
点Vjの参照セル(xi,yi,zi)中のレベルを示す。第4図
は、点Vjについて限定される関数bjの定義の局所領域を
示す。この領域は、半径Rbによって限定される円である
(Rbは、bjを定義するガウス関数の切断半径である)。
第4図は、式11により与えられる積分を示している。局
所基本関数を考慮に入れる限り、該当する検出セルの加
重関数の支持線の一部分との交差を求める。支持線のこ
の部分は通常楕円錐台であり、その最大の底面は、高さ
zj−Rbにより設定され、またその最小の底面は高さzj−
Rbにより設定される。本発明において、一方でαが非常
に小さく、他方でRbもまた小さい(「局所」底面の仮定
に基づく)ことを考慮に入れると、この円錐台は、軸SP
iを有する円筒要素により、gjの近傍に接近し得ること
は明らかになった。言い換えれば、zkが、間隔(zj -,zj
+)に属するところでは、 σgx(zk)≒σgx(zj) 故に、σgx(zk)≒σg/FMAG(zj) また、σgy(zk)≒σgy(zj) 故に、σgy(zk)≒σgcos(τ)/FMAG(zj) ・・・(26) と表すことができる。
この式では、kは加重関数に属するのに対し、jは基
本関数に属する。最後の式によってもたらされる利点
は、σgy(zk)が、近似により、zjだけに依存するの
で、σgy(zk)はzkから独立していることである。この
ような条件下で、式24の近似(f)の計算は、単純化
され、 になる。というのは、関数は可分であるので、従って積
分は可分であり、また、 となるからである。
本関数に属する。最後の式によってもたらされる利点
は、σgy(zk)が、近似により、zjだけに依存するの
で、σgy(zk)はzkから独立していることである。この
ような条件下で、式24の近似(f)の計算は、単純化
され、 になる。というのは、関数は可分であるので、従って積
分は可分であり、また、 となるからである。
xkの積分とykの積分とを分離することにより、 と表すことができる。
実際には、所与の値iの場合には、この最後の式の2
つの積分の各々を2つの有心ガウス関数のたたみこみと
して解釈することができる。2つの有心ガウス関数のた
たみこみは、同じく有心であるガウス関数として公知で
ある。このガウス関数の分散は、初期ガウス関数の2つ
の分散の和に等しい。これによって、計算後、非常に単
純な式: 従って、i(f)の場合には、 が導かれる。
つの積分の各々を2つの有心ガウス関数のたたみこみと
して解釈することができる。2つの有心ガウス関数のた
たみこみは、同じく有心であるガウス関数として公知で
ある。このガウス関数の分散は、初期ガウス関数の2つ
の分散の和に等しい。これによって、計算後、非常に単
純な式: 従って、i(f)の場合には、 が導かれる。
分散σx 2及びσy 2の式は、次の通りである(分散の
和)。
和)。
しかし、用いた関数が切頭ガウス関数である限り、こ
れらの関係が近似にすぎないことに留意しなればならな
い。最後の式は、一切の積分計算が取り除かれているの
で、かなり有利である。その結果、本発明においてこれ
まで証明してきた近似を用いて、非常に単純な計算によ
りサンプル化された関数の射影を計算することができよ
う。
れらの関係が近似にすぎないことに留意しなればならな
い。最後の式は、一切の積分計算が取り除かれているの
で、かなり有利である。その結果、本発明においてこれ
まで証明してきた近似を用いて、非常に単純な計算によ
りサンプル化された関数の射影を計算することができよ
う。
これらの式では、σgxおよびσgyは、FMAGの関数とし
てのそれらの値に置き換えることができる。最後の式
は、 σx 2=σb 2+σg 2/FMGA(zj)2 σy 2=σb 2+σg 2cos2(τ)/FMGA(zj)2 ・・・(32) と表される。
てのそれらの値に置き換えることができる。最後の式
は、 σx 2=σb 2+σg 2/FMGA(zj)2 σy 2=σb 2+σg 2cos2(τ)/FMGA(zj)2 ・・・(32) と表される。
式(30)にある値xjおよびyjは、各セルiについて、
セル基準フレームにiに関するボクセルの座標である。
セル基準フレームにiに関するボクセルの座標である。
射影アルゴリズムfを単純化するためには、ディテク
タ基準フレームに関する座標だけを使用すればさらに有
利である。このため、xi、yi′およびzi′により示され
るディテクタ平面の第2直接正規直交基準フレームを挿
入することができる。この基準フレームの第1ベクトル
は、セル基準フレームiの第1ベクトルと同じである。
これは、すでに述べたように、ベクトルxiもまたディテ
クタ平面に属しているので、可能である。ベクトルzi′
は、ディテクタ平面に属する基準フレームのベクトル
z′に平行である。ベクトルyi′は、直接基準フレーム
の第3ベクトルとして上記ベクトルxiおよびzi′から導
き出される。ベクトルxiおよびyi′が、角度θ(第2
図)(ここでθはベクトルy′と軸PPi′間の角度)を
介した平面回転によりベクトルx′およびy′から導き
出されることがわかる。しかし、実際には、説明した射
影アルゴリズム構造に関する理由から、用いる角度は角
度θではなく、近似値θ′である。角度θ′は軸y′と
軸PVj′間の角度として定義される。ここで、Vj′は
ディテクタ上のVjの射影である。この点は、第2図なら
びに第4図に示されている。Vj′の射影が、加重関数
の局所特性により、Piの近傍に来るときだけVj′が考
慮に入れられる限り、上記近似は説明される。従って、
この近似の結果、過度の再構成エラーを伴うことなく、
上記式においてθをθ′に置き換えることができる。
タ基準フレームに関する座標だけを使用すればさらに有
利である。このため、xi、yi′およびzi′により示され
るディテクタ平面の第2直接正規直交基準フレームを挿
入することができる。この基準フレームの第1ベクトル
は、セル基準フレームiの第1ベクトルと同じである。
これは、すでに述べたように、ベクトルxiもまたディテ
クタ平面に属しているので、可能である。ベクトルzi′
は、ディテクタ平面に属する基準フレームのベクトル
z′に平行である。ベクトルyi′は、直接基準フレーム
の第3ベクトルとして上記ベクトルxiおよびzi′から導
き出される。ベクトルxiおよびyi′が、角度θ(第2
図)(ここでθはベクトルy′と軸PPi′間の角度)を
介した平面回転によりベクトルx′およびy′から導き
出されることがわかる。しかし、実際には、説明した射
影アルゴリズム構造に関する理由から、用いる角度は角
度θではなく、近似値θ′である。角度θ′は軸y′と
軸PVj′間の角度として定義される。ここで、Vj′は
ディテクタ上のVjの射影である。この点は、第2図なら
びに第4図に示されている。Vj′の射影が、加重関数
の局所特性により、Piの近傍に来るときだけVj′が考
慮に入れられる限り、上記近似は説明される。従って、
この近似の結果、過度の再構成エラーを伴うことなく、
上記式においてθをθ′に置き換えることができる。
式30にある値xjおよびyjは、セル基準フレームi内で
は、基本関数bjの中心Vjの横座標および縦座標である。
座標xjおよびyjをベクトルPiVj′の平面基準フレーム
(xi,yi)中の座標の関数として表す。PiVj′の座標は
(x′j,y′j)として示す。ベクトルPiVj′は、ディ
テクタの平面上で方向SPiの座標(xj,yj)のベクトルの
射影に等しい。従って、 が導き出される。これは、 σx′ 2=σb 2・FMAG(zj)2+σg 2 σy′ 2=σb 2・FMAG(zj)2/cos(τ)+σg 2 を設定することにより、前記式を変形する最終式: を導くことができる。
は、基本関数bjの中心Vjの横座標および縦座標である。
座標xjおよびyjをベクトルPiVj′の平面基準フレーム
(xi,yi)中の座標の関数として表す。PiVj′の座標は
(x′j,y′j)として示す。ベクトルPiVj′は、ディ
テクタの平面上で方向SPiの座標(xj,yj)のベクトルの
射影に等しい。従って、 が導き出される。これは、 σx′ 2=σb 2・FMAG(zj)2+σg 2 σy′ 2=σb 2・FMAG(zj)2/cos(τ)+σg 2 を設定することにより、前記式を変形する最終式: を導くことができる。
指数関数の利点は、積分計算の代わりに、予め計算さ
れ、図表作成が可能な関数を使用することである。これ
により、時間が節約される。
れ、図表作成が可能な関数を使用することである。これ
により、時間が節約される。
得られた結果は、射影の計算用に提案されるアルゴ
リズムの全体的説明を行うことによりまとめることがで
きる。サンプルfjは、通常関数の各値jに使用可能で
ある。これらは、射影マトリクスを用いて射影され、射
影i(f)の推定がiの各値について得られることに
なる。この推定値は実際に測定された値と比較し、そこ
から差異(あるいは残り)の後方射影によりfの新しい
サンプリング値を導き出し、推定が次第に所望する関数
に出来る限り近づくまで、この計算全体を繰り返す。従
って、このアルゴリズムは次のような段階を含む。
リズムの全体的説明を行うことによりまとめることがで
きる。サンプルfjは、通常関数の各値jに使用可能で
ある。これらは、射影マトリクスを用いて射影され、射
影i(f)の推定がiの各値について得られることに
なる。この推定値は実際に測定された値と比較し、そこ
から差異(あるいは残り)の後方射影によりfの新しい
サンプリング値を導き出し、推定が次第に所望する関数
に出来る限り近づくまで、この計算全体を繰り返す。従
って、このアルゴリズムは次のような段階を含む。
−サンプルの基礎からガウス基本関数の基数へ移る。こ
れは、fjがわかっているとき、係数aj全てを計算するこ
とである。最初に、基本関数bjが結合された支持線を有
していることを考慮に入れて、この操作を単純化すると
共に、用いたガウス関数が、該方法の第1段階でaj0≒f
j0を選択することによりボクセルの寸法と同じくらいの
中央高さでの幅をもつことを指摘することができる。こ
の関係は近似的であるが、実用面では充分である。
れは、fjがわかっているとき、係数aj全てを計算するこ
とである。最初に、基本関数bjが結合された支持線を有
していることを考慮に入れて、この操作を単純化すると
共に、用いたガウス関数が、該方法の第1段階でaj0≒f
j0を選択することによりボクセルの寸法と同じくらいの
中央高さでの幅をもつことを指摘することができる。こ
の関係は近似的であるが、実用面では充分である。
−全てのi(f)をゼロに初期化し、 −jの各値について、 ・ディテクタ上の基本関数の中心Vjの射影である点V′
jのディテクタ基準フレーム内での座標を計算し、 ・そこからθ′=角度(y′,PV′j)を導き出し、 ・高さzjに位置する上記基本関数の場合に関数FMAG
(zj)を計算し、 ・SPとSPi間で形成される角度τを計算し、 ・σx′およびσy′を計算し、 ・ディテクタ上の基本関数bjの支持線の射影に属するデ
ィテクタ平面の画素(xi′,yi′)の各々の場合に、 ・xj′およびyj′を計算し、 ・基本関数bjの画素xi′、yi′への寄与であるhijを計
算し、 ・インクリメントi(f)=i(f)+hijfjを行
う。
jのディテクタ基準フレーム内での座標を計算し、 ・そこからθ′=角度(y′,PV′j)を導き出し、 ・高さzjに位置する上記基本関数の場合に関数FMAG
(zj)を計算し、 ・SPとSPi間で形成される角度τを計算し、 ・σx′およびσy′を計算し、 ・ディテクタ上の基本関数bjの支持線の射影に属するデ
ィテクタ平面の画素(xi′,yi′)の各々の場合に、 ・xj′およびyj′を計算し、 ・基本関数bjの画素xi′、yi′への寄与であるhijを計
算し、 ・インクリメントi(f)=i(f)+hijfjを行
う。
本発明によるアルゴリズムは、所望の関数fを再構成
するためだけに使用する必要はないことが理解されるで
あろう。適切にサンプル化された関数fがわかっている
とき、このアルゴリズムにより、画定すべき主方向SPで
の錐形射影における上記関数の画像を提供することが可
能となる。従って、本発明の利点は、人為的構造なしに
画像を提供する点であり、この画像は従来の技術より速
く得ることができる。
するためだけに使用する必要はないことが理解されるで
あろう。適切にサンプル化された関数fがわかっている
とき、このアルゴリズムにより、画定すべき主方向SPで
の錐形射影における上記関数の画像を提供することが可
能となる。従って、本発明の利点は、人為的構造なしに
画像を提供する点であり、この画像は従来の技術より速
く得ることができる。
基本関数の高さの関数でもある標準偏差を有するガウ
ス加重関数(上に示した全ての単純化を伴う)を選択す
ることにより、通常照射の円錐性に自然に適合した加重
が実現される。このように、アーティファクトを減少
し、再構成を加速することができる。
ス加重関数(上に示した全ての単純化を伴う)を選択す
ることにより、通常照射の円錐性に自然に適合した加重
が実現される。このように、アーティファクトを減少
し、再構成を加速することができる。
本発明が、あらゆる種類の画像もしくは計算された再
構成の処理に関することは明らかである。このような処
理操作には、特に再構成された構造あるいは計算された
画像の視覚表示が含まれる。
構成の処理に関することは明らかである。このような処
理操作には、特に再構成された構造あるいは計算された
画像の視覚表示が含まれる。
第1図は、X線照射源とマルチディテクタの概略を示
し、 第2図は、検出セルの座標を表すために用いることがで
き、単純化された計算を可能にする異なる基準フレーム
を表示する幾何学的図面であり、 第3図は、本発明に従う加重関数(3次元)を概略的に
示し(単純にするため2次元とした)、 第4図は、射影に関して、採用した単純化された計算を
証明する幾何学的特性を示す。 (主な参照符号) S……X線源、C……身体、 D……マルチディテクタ、 i……検出セル
し、 第2図は、検出セルの座標を表すために用いることがで
き、単純化された計算を可能にする異なる基準フレーム
を表示する幾何学的図面であり、 第3図は、本発明に従う加重関数(3次元)を概略的に
示し(単純にするため2次元とした)、 第4図は、射影に関して、採用した単純化された計算を
証明する幾何学的特性を示す。 (主な参照符号) S……X線源、C……身体、 D……マルチディテクタ、 i……検出セル
───────────────────────────────────────────────────── フロントページの続き (72)発明者 アンヌ ルージュ フランス国 92260 フォントゥネオ ローズ リュ アントワーヌ プティ 4ベー (72)発明者 ケネス ハンソン アメリカ合衆国 87544 ニュー メキ シコ ロス アラモス ラ センダ ロ ード 122 (56)参考文献 特開 昭63−253481(JP,A) (58)調査した分野(Int.Cl.6,DB名) G06T 1/00 A61B 6/03 G01N 23/04
Claims (10)
- 【請求項1】体積要素でサンプルされた3次元物体の画
像の計算および処理方法であって、この画像はX線によ
る錐形射影で例えば照射され、 ガウス基本関数(bj)の空間についてサンプルとなった
3次元物体を分解し、 各射影画像要素ごとに、該当する画像要素に射影される
体積要素に対応する基本関数の寄与(aj)を計算し、各
寄与が、体積要素に関する基本関数と、該当する画像要
素の照射用の支持線を表す関数(gj)との積の積分に等
しく、 射影された画像が処理される方法において 支持線を表す関数の代わりにガウス加重関数を用いるこ
とを特徴とする方法。 - 【請求項2】上記ガウス加重関数が、照射に対して垂直
の方向に標準偏差を有し、該標準偏差の値が照射半径で
の体積要素の位置の関数であることを特徴とする請求項
1記載の方法。 - 【請求項3】上記標準偏差の値が、該当する体積要素の
位置の一次関数であることを特徴とする請求項2記載の
方法。 - 【請求項4】上記標準偏差の値が、該当する体積要素の
近傍で局所的に一定であるとみなすことを特徴とする請
求項2記載の方法。 - 【請求項5】切頭ガウス加重関数が選択されることを特
徴とする請求項1記載の方法。 - 【請求項6】サンプル体積内でサンプルの幾何学的寸法
の3/2にほぼ等しい切断半径(Rg)を有する切頭ガウス
加重関数が選択されることを特徴とする請求項5記載の
方法。 - 【請求項7】反復的、可分でかつ局所的な基本関数が選
択され、 上記寄与の計算が、上記基本関数の支持線の中心(Vj)
の幾何学的射影の座標(xj′,yj′)の関数の指数関数
を計算することにより実施されることを特徴とする請求
項1記載の方法。 - 【請求項8】上記基本関数の寄与を計算するために、 σx′ 2=σb 2・FMAG(zj)2+σg 2 σy′ 2=σb 2・FMAG(zj)2+/cos(τ)+σg 2 を設定することにより、次の式: ただし、 FMAGは、座標zjにより、射影の錐形性を考慮にいれる拡
大係数を示し、 σbとσgは、基本関数と加重関数の標準偏差をそれぞ
れ示し、 δx、δy及びδzはサンプルの幾何学的寸法であり、 τは、直角射影に対して、該当する画像要素への射影の
傾斜を測定する) が与えられることを特徴とする請求項1記載の方法。 - 【請求項9】−中央高さでの幅が、サンプル体積内のサ
ンプルの幾何学的寸法にほぼ等しいガウス加重関数を選
択することを特徴とする請求項1記載の方法。 - 【請求項10】サンプル関数により物体を推定し、 推定されたサンプル関数を射影し、 射影された推定サンプル関数の射影を物体の射影による
検査から得られた測定値と比較し、 そこから物体の新しい推定値を演繹し、 比較結果が満足できるまで上記操作を繰り返す検査中の
物体の3次元再構成の方法であって、 サンプル関数を射影するために、 ガウス基本関数の空間についてサンプルされた3次元物
体を分解し、 推定されたサンプル関数の各射影画像要素について、該
当する射影画像要素に射影される物体の体積要素に対応
する基本関数の寄与を計算し、各寄与が、体積要素に関
する基本関数と、該当する画像要素の照射用の支持線を
表す関数(gj)との積の積分に等しく、 再構成された画像が処理され、 基本支持線を表す関数の代わりにガウス加重関数を用い
る ことを特徴とする方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR8900676 | 1989-01-20 | ||
FR8900676A FR2642198B1 (fr) | 1989-01-20 | 1989-01-20 | Procede de calcul et d'exploitation de l'image en projection conique, par exemple au sens des rayons x, d'un objet tridimensionnel echantillonne, et procede de reconstruction tridimensionnelle d'un objet etudie utilisant ce procede de calcul |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH02230479A JPH02230479A (ja) | 1990-09-12 |
JP2952314B2 true JP2952314B2 (ja) | 1999-09-27 |
Family
ID=9377916
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2011670A Expired - Fee Related JP2952314B2 (ja) | 1989-01-20 | 1990-01-20 | サンプル3次元物体の例えばx線による錐形射影での画像の計算及び処理方法と、該計算方法を利用した検査物体の3次元再構成方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US5218534A (ja) |
EP (1) | EP0379399B1 (ja) |
JP (1) | JP2952314B2 (ja) |
DE (1) | DE69018104T2 (ja) |
FR (1) | FR2642198B1 (ja) |
Families Citing this family (41)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2670039B1 (fr) * | 1990-11-29 | 1993-12-24 | Commissariat A Energie Atomique | Procede et dispositif de reconstruction d'images tridimentionnelles d'un objet en utilisant deux trajectoires circulaires d'acquisition. |
FR2670038B1 (fr) * | 1990-11-29 | 1993-12-24 | Commissariat A Energie Atomique | Procede et dispositif de reconstruction d'images tridimentionnelles d'un objet en utilisant deux trajectoires circulaires d'axe commun. |
US5365560A (en) * | 1991-07-29 | 1994-11-15 | General Electric Company | Method and apparatus for acquiring a uniform distribution of radon data sufficiently dense to constitute a complete set for exact image reconstruction of an object irradiated by a cone beam source |
US5333164A (en) * | 1991-12-11 | 1994-07-26 | General Electric Company | Method and apparatus for acquiring and processing only a necessary volume of radon data consistent with the overall shape of the object for efficient three dimensional image reconstruction |
US5611026A (en) * | 1992-12-21 | 1997-03-11 | General Electric Company | Combining a priori data with partial scan data to project three dimensional imaging of arbitrary objects with computerized tomography |
FR2700909B1 (fr) * | 1993-01-27 | 1995-03-17 | Gen Electric Cgr | Dispositif et procédé automatique de calibration géométrique d'un système d'imagerie par rayons X. |
FR2705223A1 (fr) * | 1993-05-13 | 1994-11-25 | Ge Medical Syst Sa | Procédé d'acquisition d'images d'un corps par placement en rotation. |
JP3490505B2 (ja) * | 1994-09-05 | 2004-01-26 | 株式会社東芝 | X線診断装置 |
US5804881A (en) * | 1995-11-27 | 1998-09-08 | Motorola, Inc. | Method and assembly for providing improved underchip encapsulation |
US6055449A (en) * | 1997-09-22 | 2000-04-25 | Siemens Corporate Research, Inc. | Method for localization of a biopsy needle or similar surgical tool in a radiographic image |
US6373487B1 (en) * | 1999-09-17 | 2002-04-16 | Hewlett-Packard Company | Methods and apparatus for constructing a 3D model of a scene from calibrated images of the scene |
US7129942B2 (en) * | 2002-12-10 | 2006-10-31 | International Business Machines Corporation | System and method for performing domain decomposition for multiresolution surface analysis |
JP4148041B2 (ja) * | 2003-06-27 | 2008-09-10 | ソニー株式会社 | 信号処理装置および信号処理方法、並びにプログラムおよび記録媒体 |
US7693318B1 (en) | 2004-01-12 | 2010-04-06 | Pme Ip Australia Pty Ltd | Method and apparatus for reconstruction of 3D image volumes from projection images |
US20050270298A1 (en) * | 2004-05-14 | 2005-12-08 | Mercury Computer Systems, Inc. | Daughter card approach to employing multiple graphics cards within a system |
US8189002B1 (en) | 2004-10-29 | 2012-05-29 | PME IP Australia Pty, Ltd. | Method and apparatus for visualizing three-dimensional and higher-dimensional image data sets |
US7778392B1 (en) | 2004-11-02 | 2010-08-17 | Pme Ip Australia Pty Ltd | Method of reconstructing computed tomography (CT) volumes suitable for execution on commodity central processing units (CPUs) and graphics processors, and apparatus operating in accord with those methods (rotational X-ray on GPUs) |
US7272205B2 (en) * | 2004-11-17 | 2007-09-18 | Purdue Research Foundation | Methods, apparatus, and software to facilitate computing the elements of a forward projection matrix |
US7609884B1 (en) | 2004-12-23 | 2009-10-27 | Pme Ip Australia Pty Ltd | Mutual information based registration of 3D-image volumes on GPU using novel accelerated methods of histogram computation |
US7623732B1 (en) | 2005-04-26 | 2009-11-24 | Mercury Computer Systems, Inc. | Method and apparatus for digital image filtering with discrete filter kernels using graphics hardware |
US7853061B2 (en) * | 2007-04-26 | 2010-12-14 | General Electric Company | System and method to improve visibility of an object in an imaged subject |
US8019151B2 (en) * | 2007-06-11 | 2011-09-13 | Visualization Sciences Group, Inc. | Methods and apparatus for image compression and decompression using graphics processing unit (GPU) |
CN101371786B (zh) * | 2007-08-24 | 2011-01-12 | 北京师范大学珠海分校 | 一种x射线图像三维重构的方法及*** |
US8392529B2 (en) | 2007-08-27 | 2013-03-05 | Pme Ip Australia Pty Ltd | Fast file server methods and systems |
WO2009067680A1 (en) | 2007-11-23 | 2009-05-28 | Mercury Computer Systems, Inc. | Automatic image segmentation methods and apparartus |
US10311541B2 (en) | 2007-11-23 | 2019-06-04 | PME IP Pty Ltd | Multi-user multi-GPU render server apparatus and methods |
WO2009067675A1 (en) | 2007-11-23 | 2009-05-28 | Mercury Computer Systems, Inc. | Client-server visualization system with hybrid data processing |
US8319781B2 (en) | 2007-11-23 | 2012-11-27 | Pme Ip Australia Pty Ltd | Multi-user multi-GPU render server apparatus and methods |
US9904969B1 (en) | 2007-11-23 | 2018-02-27 | PME IP Pty Ltd | Multi-user multi-GPU render server apparatus and methods |
JP5670050B2 (ja) * | 2009-02-05 | 2015-02-18 | 株式会社東芝 | 画像再構成装置および画像再構成方法 |
US8681902B2 (en) * | 2011-09-16 | 2014-03-25 | Alcatel Lucent | Method and apparatus for low complexity robust reconstruction of noisy signals |
US10540803B2 (en) | 2013-03-15 | 2020-01-21 | PME IP Pty Ltd | Method and system for rule-based display of sets of images |
US11244495B2 (en) | 2013-03-15 | 2022-02-08 | PME IP Pty Ltd | Method and system for rule based display of sets of images using image content derived parameters |
US8976190B1 (en) | 2013-03-15 | 2015-03-10 | Pme Ip Australia Pty Ltd | Method and system for rule based display of sets of images |
US10070839B2 (en) | 2013-03-15 | 2018-09-11 | PME IP Pty Ltd | Apparatus and system for rule based visualization of digital breast tomosynthesis and other volumetric images |
US9509802B1 (en) | 2013-03-15 | 2016-11-29 | PME IP Pty Ltd | Method and system FPOR transferring data to improve responsiveness when sending large data sets |
US11183292B2 (en) | 2013-03-15 | 2021-11-23 | PME IP Pty Ltd | Method and system for rule-based anonymized display and data export |
WO2014207080A1 (en) * | 2013-06-28 | 2014-12-31 | Koninklijke Philips N.V. | Methods for generation of edge-preserving synthetic mammograms from tomosynthesis data |
US9984478B2 (en) | 2015-07-28 | 2018-05-29 | PME IP Pty Ltd | Apparatus and method for visualizing digital breast tomosynthesis and other volumetric images |
US11599672B2 (en) | 2015-07-31 | 2023-03-07 | PME IP Pty Ltd | Method and apparatus for anonymized display and data export |
US10909679B2 (en) | 2017-09-24 | 2021-02-02 | PME IP Pty Ltd | Method and system for rule based display of sets of images using image content derived parameters |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4128877A (en) * | 1977-03-31 | 1978-12-05 | Katz Myron B | Method of obtaining a uniquely determined cross-sectional reconstruction from X-ray projections |
US4630203A (en) * | 1983-12-27 | 1986-12-16 | Thomas Szirtes | Contour radiography: a system for determining 3-dimensional contours of an object from its 2-dimensional images |
US4952922A (en) * | 1985-07-18 | 1990-08-28 | Hughes Aircraft Company | Predictive look ahead memory management for computer image generation in simulators |
US4833698A (en) * | 1986-02-24 | 1989-05-23 | Exxon Research And Engineering Company | Apparatus for three dimensional tomography utilizing an electro-optic x-ray detector |
US4835712A (en) * | 1986-04-14 | 1989-05-30 | Pixar | Methods and apparatus for imaging volume data with shading |
US4751643A (en) * | 1986-08-04 | 1988-06-14 | General Electric Company | Method and apparatus for determining connected substructures within a body |
US4791567A (en) * | 1986-09-15 | 1988-12-13 | General Electric Company | Three dimensional connectivity system employing an equivalence schema for determining connected substructures within a body |
US4821213A (en) * | 1986-12-19 | 1989-04-11 | General Electric Co. | System for the simultaneous display of two or more internal surfaces within a solid object |
US4989142A (en) * | 1988-10-24 | 1991-01-29 | General Electric Company | Three-dimensional images obtained from tomographic slices with gantry tilt |
-
1989
- 1989-01-20 FR FR8900676A patent/FR2642198B1/fr not_active Expired - Fee Related
-
1990
- 1990-01-12 DE DE69018104T patent/DE69018104T2/de not_active Expired - Fee Related
- 1990-01-12 EP EP90400083A patent/EP0379399B1/fr not_active Expired - Lifetime
- 1990-01-17 US US07/466,590 patent/US5218534A/en not_active Expired - Fee Related
- 1990-01-20 JP JP2011670A patent/JP2952314B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
US5218534A (en) | 1993-06-08 |
DE69018104T2 (de) | 1995-10-26 |
FR2642198B1 (fr) | 1991-04-05 |
FR2642198A1 (fr) | 1990-07-27 |
JPH02230479A (ja) | 1990-09-12 |
EP0379399B1 (fr) | 1995-03-29 |
DE69018104D1 (de) | 1995-05-04 |
EP0379399A1 (fr) | 1990-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2952314B2 (ja) | サンプル3次元物体の例えばx線による錐形射影での画像の計算及び処理方法と、該計算方法を利用した検査物体の3次元再構成方法 | |
US5909476A (en) | Iterative process for reconstructing cone-beam tomographic images | |
US9613442B2 (en) | Image reconstruction from limited or incomplete data | |
US5559335A (en) | Rotating and warping projector/backprojector for converging-beam geometries | |
Jia et al. | GPU-based fast low-dose cone beam CT reconstruction via total variation | |
US5930384A (en) | Process for the reconstruction of a 3D image with contrast and resolution improvements and application of said process to the production of an attentuation cartography of an object | |
US6907102B1 (en) | Iterative reconstruction methods for multi-slice computed tomography | |
EP1800264B1 (en) | Image reconstruction with voxel dependent interpolation | |
US20120051626A1 (en) | Three-dimensional forward and back projection methods | |
Wieczorek | The image quality of FBP and MLEM reconstruction | |
JP2011503570A (ja) | 減衰マップを形成するための装置及び方法 | |
Defrise et al. | Image reconstruction from truncated, two-dimensional, parallel projections | |
US8335358B2 (en) | Method and system for reconstructing a medical image of an object | |
Natterer et al. | Past and future directions in x‐ray computed tomography (CT) | |
Brasse et al. | Fast fully 3-D image reconstruction in PET using planograms | |
Liang | Compensation for attenuation, scatter, and detector response in SPECT reconstruction via iterative FBP methods | |
Gindi et al. | Imaging with rotating slit apertures and rotating collimators | |
Saint-Félix et al. | 3D reconstruction of high contrast objects using a multi-scale detection/estimation scheme | |
Stoddart et al. | New multidimensional reconstructions for the 12-detector, scanned focal point, single-photon tomograph | |
Grangeat et al. | Indirect cone-beam three-dimensional image reconstruction | |
Grangeat et al. | Evaluation of the 3-D Radon transform algorithm for cone beam reconstruction | |
Gopi et al. | Iterative computed tomography reconstruction from sparse-view data | |
US6973156B2 (en) | Method and apparatus for reconstruction of the attenuation density of an object from X-ray projection image data | |
Hanson | Computed tomographic (CT) reconstruction from limited projection angles | |
Rajgopal et al. | Image reconstruction from incomplete projection data: A linear prediction approach |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
LAPS | Cancellation because of no payment of annual fees |