JPH0771941A - ３次元像を再構築する装置及び方法 - Google Patents
３次元像を再構築する装置及び方法Info
- Publication number
- JPH0771941A JPH0771941A JP6142197A JP14219794A JPH0771941A JP H0771941 A JPH0771941 A JP H0771941A JP 6142197 A JP6142197 A JP 6142197A JP 14219794 A JP14219794 A JP 14219794A JP H0771941 A JPH0771941 A JP H0771941A
- Authority
- JP
- Japan
- Prior art keywords
- array
- axis
- detectors
- dimensional image
- spiral
- 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.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims description 30
- 238000005259 measurement Methods 0.000 claims description 44
- 230000009466 transformation Effects 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 6
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims description 5
- 229910052704 radon Inorganic materials 0.000 claims description 4
- 230000008521 reorganization Effects 0.000 claims description 3
- 238000012887 quadratic function Methods 0.000 claims 1
- 230000005855 radiation Effects 0.000 abstract description 8
- 238000007689 inspection Methods 0.000 abstract description 6
- 239000000463 material Substances 0.000 abstract description 4
- 239000011159 matrix material Substances 0.000 abstract description 2
- 238000012545 processing Methods 0.000 description 7
- 238000003491 array Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 230000014509 gene expression Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000008707 rearrangement Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 230000001131 transforming effect 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/027—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Pure & Applied Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Stereoscopic And Panoramic Photography (AREA)
Abstract
(57)【要約】 （修正有）
【目的】 焦点が点ではないが、完全に異なる処理動作
を必要とする平行な扇形焦点合わせモードではなく、純
円錐状焦点合わせモードを伴う像再構築アルゴリズムを
ほとんど改造なしに用いることができる装置を提供す
る。 【構成】 検出器２の平面アレイ１によって対象物４の
３次元像を再構築する装置及び方法である。検出器の焦
点は円Ｃｅ上に位置している。単一点に焦点合わせする
場合よりもより広い検査領域を得ることができる。単一
点への円錐焦点合わせによる特性測定結果を反転するア
ルゴリズムを用いる。医学や産業上の非破壊検査に適用
される。
を必要とする平行な扇形焦点合わせモードではなく、純
円錐状焦点合わせモードを伴う像再構築アルゴリズムを
ほとんど改造なしに用いることができる装置を提供す
る。 【構成】 検出器２の平面アレイ１によって対象物４の
３次元像を再構築する装置及び方法である。検出器の焦
点は円Ｃｅ上に位置している。単一点に焦点合わせする
場合よりもより広い検査領域を得ることができる。単一
点への円錐焦点合わせによる特性測定結果を反転するア
ルゴリズムを用いる。医学や産業上の非破壊検査に適用
される。
Description
【０００１】
【産業上の利用分野】本発明は、３次元像を再構築する
装置及び方法に関する。
装置及び方法に関する。
【０００２】
【従来の技術及び発明が解決しようとする課題】特に医
療又は工業上の非破壊検査の分野において、断層法によ
りその像を再構築しようとする対象物は、その対象物の
回りを回転する２次元アレイセンサの前に置かれる。こ
の２次元アレイセンサは、その各検出器によって走査さ
れる空間を単一の線に規定する例えば鉛又はこれに類似
した材料の板からなるコリメータによって覆われてい
る。このアレイセンサは、対象物に関する円又は他の閉
軌道又は螺旋軌道上を通過する。幾つかの異なる入射角
それぞれによる２次元像について得られた検出器の測定
結果を照合し、測定結果の数学的変換の反転を行う解析
手法又は方程式の体系を反転する代数的解法を用いるこ
とによって対象物の３次元像を再構築することが可能と
なる。これら検出器によって測定される対象物の関数、
即ち像を特徴付ける特性は、対象物の各点から放出され
る放射か又は対象物の向こう側の焦線上に存在する放射
源から放出される放射についての対象物の各点による減
衰のどちらかによるものである。
療又は工業上の非破壊検査の分野において、断層法によ
りその像を再構築しようとする対象物は、その対象物の
回りを回転する２次元アレイセンサの前に置かれる。こ
の２次元アレイセンサは、その各検出器によって走査さ
れる空間を単一の線に規定する例えば鉛又はこれに類似
した材料の板からなるコリメータによって覆われてい
る。このアレイセンサは、対象物に関する円又は他の閉
軌道又は螺旋軌道上を通過する。幾つかの異なる入射角
それぞれによる２次元像について得られた検出器の測定
結果を照合し、測定結果の数学的変換の反転を行う解析
手法又は方程式の体系を反転する代数的解法を用いるこ
とによって対象物の３次元像を再構築することが可能と
なる。これら検出器によって測定される対象物の関数、
即ち像を特徴付ける特性は、対象物の各点から放出され
る放射か又は対象物の向こう側の焦線上に存在する放射
源から放出される放射についての対象物の各点による減
衰のどちらかによるものである。
【０００３】検出器の焦点、即ち焦線の交差点、の位置
は、種々の形態に想定できる。このため、各２次元像に
対して単一の焦点が存在する純円錐状のコリメータを用
いることが提案されいるが、これは点放射源のみの使用
を可能としている。本願の発明者の１人及び他の科学者
によってなされた先行する研究、特にヨーロッパ特許公
開公報第２９２４０２号は、像又はより一般的には平面
上の関数の和を再構築するために、関数のラドン（Ｒａ
ｄｏｎ）変換の１次又は最高次微分と称されるものの反
転を用いたときのこの解法の価値を検認している。
は、種々の形態に想定できる。このため、各２次元像に
対して単一の焦点が存在する純円錐状のコリメータを用
いることが提案されいるが、これは点放射源のみの使用
を可能としている。本願の発明者の１人及び他の科学者
によってなされた先行する研究、特にヨーロッパ特許公
開公報第２９２４０２号は、像又はより一般的には平面
上の関数の和を再構築するために、関数のラドン（Ｒａ
ｄｏｎ）変換の１次又は最高次微分と称されるものの反
転を用いたときのこの解法の価値を検認している。
【０００４】検出器アレイの回転軸に平行な直線上に検
出器の焦点を合わせること、即ち焦線を計画された平面
上で扇形にかつ軌道に平行に配置することも提案されて
いる。幾つかの文献に記載されているこの「平行な扇
形」の焦点合わせモードは、回転軸内での拡大ができな
いことに基づく不都合を有している。
出器の焦点を合わせること、即ち焦線を計画された平面
上で扇形にかつ軌道に平行に配置することも提案されて
いる。幾つかの文献に記載されているこの「平行な扇
形」の焦点合わせモードは、回転軸内での拡大ができな
いことに基づく不都合を有している。
【０００５】米国特許第４８２０９２４号及び第４８２
３０１７号には、さらに複雑な焦点位置について提案が
なされている。この場合、検出器の水平及び垂直列（行
及び列）のコリメータ線は、少なくとも幾つかの場合
に、アレイに垂直な線上に収束する。１つ実施態様で
は、焦点位置はアレイに平行でありかつその回転軸に垂
直な直線である。しかしながら、その後の再構築の計算
が非常に複雑となる。
３０１７号には、さらに複雑な焦点位置について提案が
なされている。この場合、検出器の水平及び垂直列（行
及び列）のコリメータ線は、少なくとも幾つかの場合
に、アレイに垂直な線上に収束する。１つ実施態様で
は、焦点位置はアレイに平行でありかつその回転軸に垂
直な直線である。しかしながら、その後の再構築の計算
が非常に複雑となる。
【０００６】１９９３年、カールスルーエ、ＥＣＡＰＴ
報告によって出版されたミューラー、アース及びブラー
クによる「コンピュータ化２次元断層法における打ち切
りアーチファクト及び合成スキャナアレイ」なる文献に
は、３次元像の再構築方法が記載されている。この方法
は、測定における打ち切り（ｔｒｕｎｃａｔｉｏｎ）の
影響を低減するために、投射を虚像投射に変換する段階
を有している。検出器アレイは、その検出器アレイの円
形の移動の結果、円錐幾何学における複数の欠けた投射
を得ることを可能とする。これは、放射源には無関係で
あり、虚像投射は平面の虚像センサアレイ上で生じる。
本発明は、実像投射から虚像投射への変換を含んでいる
が、上述の文献においては焦点合わせが円錐又は平行扇
形であるので、この変換は異なる獲得条件となる。
報告によって出版されたミューラー、アース及びブラー
クによる「コンピュータ化２次元断層法における打ち切
りアーチファクト及び合成スキャナアレイ」なる文献に
は、３次元像の再構築方法が記載されている。この方法
は、測定における打ち切り（ｔｒｕｎｃａｔｉｏｎ）の
影響を低減するために、投射を虚像投射に変換する段階
を有している。検出器アレイは、その検出器アレイの円
形の移動の結果、円錐幾何学における複数の欠けた投射
を得ることを可能とする。これは、放射源には無関係で
あり、虚像投射は平面の虚像センサアレイ上で生じる。
本発明は、実像投射から虚像投射への変換を含んでいる
が、上述の文献においては焦点合わせが円錐又は平行扇
形であるので、この変換は異なる獲得条件となる。
【０００７】本発明は、焦点が点ではないが、完全に異
なる処理動作を必要とする平行な扇形焦点合わせモード
ではなく、純円錐状焦点合わせモードを伴う像再構築ア
ルゴリズムをほとんど改造なしに用いることができる装
置に関するものである。
なる処理動作を必要とする平行な扇形焦点合わせモード
ではなく、純円錐状焦点合わせモードを伴う像再構築ア
ルゴリズムをほとんど改造なしに用いることができる装
置に関するものである。
【０００８】点焦点又は点のような焦点固有の欠点は、
検査領域が焦点に向かって非常に急激に狭くなることで
ある。この焦点が対象物に比較的近い場合、この対象物
の一部が全ての測定の放射領域内に通常は含まれず、欠
けた測定が行われることとなって反転式では訂正できな
い再構築エラーの発生を招く。焦点が対象物から離れて
いれば即ち対象物が検出器アレイに近ければ、その対象
物は通常は全て検査円錐内に含まれが、像の拡大倍率が
小さくなる。これは、人体内の臓器のように対象物の一
部を検査する必要がある場合に、ほとんどの検出器が対
象物のそれ以外の部分の測定に用いられるてしまうので
不都合となる。
検査領域が焦点に向かって非常に急激に狭くなることで
ある。この焦点が対象物に比較的近い場合、この対象物
の一部が全ての測定の放射領域内に通常は含まれず、欠
けた測定が行われることとなって反転式では訂正できな
い再構築エラーの発生を招く。焦点が対象物から離れて
いれば即ち対象物が検出器アレイに近ければ、その対象
物は通常は全て検査円錐内に含まれが、像の拡大倍率が
小さくなる。これは、人体内の臓器のように対象物の一
部を検査する必要がある場合に、ほとんどの検出器が対
象物のそれ以外の部分の測定に用いられるてしまうので
不都合となる。
【０００９】
【課題を解決するための手段】本発明によれば、対象物
に向かって配設された検出器の２次元アレイを備えてお
り、検出器アレイが螺旋軌道に沿って軸の回りを回転し
つつ対象物に対して移動するべく構成されている、対象
物の３次元像を再構築する装置が提供される。特に本発
明によれば、検出器は、軸に中心を有する螺旋上にほぼ
位置している複数の点上に軸方向に焦点合わせされてお
り、検出器によってなされた測定結果を格納しかつ測定
結果を組み合わせるための本質的に公知の手段が像再構
築に用いられている。像を組み合わせるとは、例えば測
定結果の再配列及び処理並びに３次元像再構築処理等の
一連の段階を意味している。検査プリズムは円錐よりか
なり幅広である。螺旋軌道が円形軌道であり螺旋が円で
あることも可能である。これは、螺旋のピッチが零の場
合に対応する。これは非常に重要である。即ち、装置は
通常このような円軌道を実現するように設計され、再構
築処理及び演算についても実際には通常の場合と同じで
あるからである。
に向かって配設された検出器の２次元アレイを備えてお
り、検出器アレイが螺旋軌道に沿って軸の回りを回転し
つつ対象物に対して移動するべく構成されている、対象
物の３次元像を再構築する装置が提供される。特に本発
明によれば、検出器は、軸に中心を有する螺旋上にほぼ
位置している複数の点上に軸方向に焦点合わせされてお
り、検出器によってなされた測定結果を格納しかつ測定
結果を組み合わせるための本質的に公知の手段が像再構
築に用いられている。像を組み合わせるとは、例えば測
定結果の再配列及び処理並びに３次元像再構築処理等の
一連の段階を意味している。検査プリズムは円錐よりか
なり幅広である。螺旋軌道が円形軌道であり螺旋が円で
あることも可能である。これは、螺旋のピッチが零の場
合に対応する。これは非常に重要である。即ち、装置は
通常このような円軌道を実現するように設計され、再構
築処理及び演算についても実際には通常の場合と同じで
あるからである。
【００１０】演算を容易にするために、軸に平行な列上
に位置するアレイの全ての検出器が、複数の焦点のうち
の単一の焦点に焦点合わせされていることが好ましい。
に位置するアレイの全ての検出器が、複数の焦点のうち
の単一の焦点に焦点合わせされていることが好ましい。
【００１１】本発明はさらに、特に上述した装置によっ
てなされた測定結果を、純円錐焦点合わせ（虚像焦点合
わせと以下称される）の特性データに変換することによ
って使用可能とした方法にも関している。この変換は、
実際の測定値に関するものではなく、これらを定義する
幾何学的なパラメータに係るものである。その最も一般
的な形態として、本発明は、螺旋軌道に従って軸の回り
を回転しつつ対象物に対して移動する検出器の少なくと
も１つの２次元アレイによって対象物の２次元像を取る
ことを含んでおり、２次元像がコリメータ線に沿った対
象物の特性に関する各検出器の測定結果から構成されて
おり、これら測定結果を用いて式を解くことによって３
次元像を再構築する方法に関する。本発明によれば、特
に、各コリメータ線が軸を中心とする螺旋上をほぼ通っ
ており、本発明の方法は、上述の式を解く前に行われ
る、コリメータ線をパラメータで表わす変換演算段階と
再編制段階とを含んでいる。各コリメータ線が螺旋の単
一の焦点上に焦点合わせされた検出器の虚像２次元アレ
イの始点と見なされており、虚像アレイが軸についてこ
れから一定距離で回転するものであり、変換演算が虚像
アレイの各コリメータ線について始点となる実像アレイ
内における位置及び実像アレイの軸に対する角度位置を
決定するものであり、次いで、虚像円錐投影を得るため
に再編制段階において変換演算の結果に応じて測定結果
を再編制することからなり、上述の解が再編制された測
定結果に関連している。測定結果について解く方法が、
対象物の像を特徴付ける関数のラドン変換の１次微分の
反転を含むことが有利である。
てなされた測定結果を、純円錐焦点合わせ（虚像焦点合
わせと以下称される）の特性データに変換することによ
って使用可能とした方法にも関している。この変換は、
実際の測定値に関するものではなく、これらを定義する
幾何学的なパラメータに係るものである。その最も一般
的な形態として、本発明は、螺旋軌道に従って軸の回り
を回転しつつ対象物に対して移動する検出器の少なくと
も１つの２次元アレイによって対象物の２次元像を取る
ことを含んでおり、２次元像がコリメータ線に沿った対
象物の特性に関する各検出器の測定結果から構成されて
おり、これら測定結果を用いて式を解くことによって３
次元像を再構築する方法に関する。本発明によれば、特
に、各コリメータ線が軸を中心とする螺旋上をほぼ通っ
ており、本発明の方法は、上述の式を解く前に行われ
る、コリメータ線をパラメータで表わす変換演算段階と
再編制段階とを含んでいる。各コリメータ線が螺旋の単
一の焦点上に焦点合わせされた検出器の虚像２次元アレ
イの始点と見なされており、虚像アレイが軸についてこ
れから一定距離で回転するものであり、変換演算が虚像
アレイの各コリメータ線について始点となる実像アレイ
内における位置及び実像アレイの軸に対する角度位置を
決定するものであり、次いで、虚像円錐投影を得るため
に再編制段階において変換演算の結果に応じて測定結果
を再編制することからなり、上述の解が再編制された測
定結果に関連している。測定結果について解く方法が、
対象物の像を特徴付ける関数のラドン変換の１次微分の
反転を含むことが有利である。
【００１２】鮮鋭な像を得るために、全ての焦線が螺旋
又は場合によっては円を通ることが好ましい。しかしな
がら、焦線が円より僅かにずれても像が僅かに劣化する
のみであるから、この条件は必須のものではない。公知
のコリメータによると、全ての焦線が円に接する直線を
通る場合に満足する結果が得られる。
又は場合によっては円を通ることが好ましい。しかしな
がら、焦線が円より僅かにずれても像が僅かに劣化する
のみであるから、この条件は必須のものではない。公知
のコリメータによると、全ての焦線が円に接する直線を
通る場合に満足する結果が得られる。
【００１３】必ずしも必要ではないが、焦点合わせされ
る螺旋のピッチを軌道のピッチと同一にすることができ
る。
る螺旋のピッチを軌道のピッチと同一にすることができ
る。
【００１４】
【実施例】以下本発明について、これに限定されない実
施例及び添付の図面に関連して詳細に説明する。
施例及び添付の図面に関連して詳細に説明する。
【００１５】図１は、回転軸に垂直な行及び回転軸に平
行な列にマトリクス状に配置された検出器２のアレイ１
を示している。これら検出器２は、アレイ１の中心Ｏか
ら始まるアレイ１の平面内の行座標ｐ及び列座標ｑで識
別される。
行な列にマトリクス状に配置された検出器２のアレイ１
を示している。これら検出器２は、アレイ１の中心Ｏか
ら始まるアレイ１の平面内の行座標ｐ及び列座標ｑで識
別される。
【００１６】コリメータ３がアレイ１を覆っている。図
には詳しく示されていないが、コリメータ３は、一般的
な方法により、検出器２によって測定されるべき放射を
透過しない材料で形成されており複数のウィンドウが設
けられている。材料としては、通常は鉛が用いられる。
ウィンドウの数は検出器２の数以上であり、各ウィンド
ウは各検出器が検査すべき対象物４に及びこれを越えて
配向された焦線Ｌ上の空間を走査できるように各検出器
の前に配置されている。コリメータは、検出器の数より
多いウィンドウを備えていることが有利である。ウィン
ドウは、同一のサイズでありかつコリメータの表面に渡
って均等に分布配置されていることが有利である。
には詳しく示されていないが、コリメータ３は、一般的
な方法により、検出器２によって測定されるべき放射を
透過しない材料で形成されており複数のウィンドウが設
けられている。材料としては、通常は鉛が用いられる。
ウィンドウの数は検出器２の数以上であり、各ウィンド
ウは各検出器が検査すべき対象物４に及びこれを越えて
配向された焦線Ｌ上の空間を走査できるように各検出器
の前に配置されている。コリメータは、検出器の数より
多いウィンドウを備えていることが有利である。ウィン
ドウは、同一のサイズでありかつコリメータの表面に渡
って均等に分布配置されていることが有利である。
【００１７】検出器の前に置かれたウィンドウは、検出
器が単一の焦線Ｌ上の空間を走査できるように隣接する
穴の複数の方向を有している。線Ｌの方向は、検出器の
前のウィンドウの複数の方向の平均方向である。
器が単一の焦線Ｌ上の空間を走査できるように隣接する
穴の複数の方向を有している。線Ｌの方向は、検出器の
前のウィンドウの複数の方向の平均方向である。
【００１８】アレイ１の回転軸はＸで表わされる。本発
明の必須の特徴点は、対象物４の回りに位置しており軸
Ｘに垂直な平面に属する円Ｃｅに全ての焦線Ｌが導かれ
ることにある。Ｘ軸が円Ｃｅの中心を通過することが好
ましい。螺旋軌道の場合は、全ての焦線Ｌが螺旋状の弧
に導かれる。この平面は焦平面と呼ばれる。さらに本発
明では、各列に位置する全ての検出器２の焦線Ｌは、そ
れぞれの単一点Ｐに収束する。焦点Ｐのこの構成によれ
ば、走査される領域が大幅に広がり、対象物４をかなり
容易にその領域に含ませ得ることは明らかである。放射
の放出及び減衰に対するこれら焦点の貢献も、より正当
に測定されている。焦線についてのこのような構成を可
能とするコリメータは、多焦点幾何学又は多焦点コリメ
ータと呼ばれる。この種のコリメータは、現存する断層
撮影装置の平行又は円錐コリメータに置換可能である。
明の必須の特徴点は、対象物４の回りに位置しており軸
Ｘに垂直な平面に属する円Ｃｅに全ての焦線Ｌが導かれ
ることにある。Ｘ軸が円Ｃｅの中心を通過することが好
ましい。螺旋軌道の場合は、全ての焦線Ｌが螺旋状の弧
に導かれる。この平面は焦平面と呼ばれる。さらに本発
明では、各列に位置する全ての検出器２の焦線Ｌは、そ
れぞれの単一点Ｐに収束する。焦点Ｐのこの構成によれ
ば、走査される領域が大幅に広がり、対象物４をかなり
容易にその領域に含ませ得ることは明らかである。放射
の放出及び減衰に対するこれら焦点の貢献も、より正当
に測定されている。焦線についてのこのような構成を可
能とするコリメータは、多焦点幾何学又は多焦点コリメ
ータと呼ばれる。この種のコリメータは、現存する断層
撮影装置の平行又は円錐コリメータに置換可能である。
【００１９】図２は、アレイ１及び焦平面上の焦線Ｌの
垂直投影図である。Ｆが焦軸、即ちアレイ１の中心Ｏ及
びＸ軸を通過する円Ｃｅの直径を通る軸であり、そして
ｆが焦線Ｆ上で円Ｃｅ及びアレイ１間の測定した距離で
あるとすると、焦平面上の焦線Ｌｍの投影は複数の横焦
点Ｆｔｍにおいて焦軸Ｆと交差する。焦線Ｌｍについて
考察すると、前述した検出器２ｍの行座標がｐｍである
とすると、横焦点Ｆｔｍ及びアレイ１間の横焦点距離Ｄ
（ｐｍ）は、ｐｍの２次式又はｆを越えて少なくとも加
算する式、即ち次の（１）式で定義される。ただし、ｋ
は一定の係数である。 Ｄ（ｐｍ）＝ｆ＋ｋ・ｐｍ2 （１）
垂直投影図である。Ｆが焦軸、即ちアレイ１の中心Ｏ及
びＸ軸を通過する円Ｃｅの直径を通る軸であり、そして
ｆが焦線Ｆ上で円Ｃｅ及びアレイ１間の測定した距離で
あるとすると、焦平面上の焦線Ｌｍの投影は複数の横焦
点Ｆｔｍにおいて焦軸Ｆと交差する。焦線Ｌｍについて
考察すると、前述した検出器２ｍの行座標がｐｍである
とすると、横焦点Ｆｔｍ及びアレイ１間の横焦点距離Ｄ
（ｐｍ）は、ｐｍの２次式又はｆを越えて少なくとも加
算する式、即ち次の（１）式で定義される。ただし、ｋ
は一定の係数である。 Ｄ（ｐｍ）＝ｆ＋ｋ・ｐｍ2 （１）
【００２０】これにより、対象物４の全幅（この幅は僅
かに収束するのみである外側の焦線Ｌによって定義され
る）をカバーするという要求と対象物４の中心（より大
幅に収束する中央焦線Ｌによって定義される）に関して
より良い解像度及び感度を得ようとする常なる重要な要
求との間の良好な折衷案を得ることができる。
かに収束するのみである外側の焦線Ｌによって定義され
る）をカバーするという要求と対象物４の中心（より大
幅に収束する中央焦線Ｌによって定義される）に関して
より良い解像度及び感度を得ようとする常なる重要な要
求との間の良好な折衷案を得ることができる。
【００２１】再び図１に戻ると、Ｘ軸及び対象物４の回
りのアレイ１による一連の測定が検出器２による多数の
測定結果を提供することは明らかであり、円錐焦点合わ
せについて既に存在する方法と同じ方法によって、対象
物４の３次元像を再構築するためにこれら測定結果を配
列することができる。
りのアレイ１による一連の測定が検出器２による多数の
測定結果を提供することは明らかであり、円錐焦点合わ
せについて既に存在する方法と同じ方法によって、対象
物４の３次元像を再構築するためにこれら測定結果を配
列することができる。
【００２２】この像再構築は、全測定を純円錐の虚像焦
点合わせ法で行われる測定に減らすことによって効果的
に実施される。図３は、アレイ１のある位置における、
焦点Ｐに向かって配置されている検出器２ｍの焦線Ｌｍ
を示している。アレイ１から半径Ｒの円Ｃｅまでの距離
が、再び、ｆで示されている。円錐焦点合わせによる虚
像検出器アレイ１０１を想定することにより全てのこれ
ら検出器は点Ｐに向かって収束し、さらに虚像アレイ１
０１に垂直でありＰ及びＣを通る線上での点Ｐからの距
離ｆを選ぶことによりこの虚像アレイ１０１は単一形態
となる。焦線Ｌｍは、点Ｐの虚像アレイ１０１に対する
垂直投影で定義される中心Ｏ′から行軸に沿って距離ｐ
ｃかつ列軸に沿って距離ｑｃの位置におけるこの虚像ア
レイ１０１まで伸長される。２つのアレイ１及び１０１
間の角度はα（ｐｍ）である。
点合わせ法で行われる測定に減らすことによって効果的
に実施される。図３は、アレイ１のある位置における、
焦点Ｐに向かって配置されている検出器２ｍの焦線Ｌｍ
を示している。アレイ１から半径Ｒの円Ｃｅまでの距離
が、再び、ｆで示されている。円錐焦点合わせによる虚
像検出器アレイ１０１を想定することにより全てのこれ
ら検出器は点Ｐに向かって収束し、さらに虚像アレイ１
０１に垂直でありＰ及びＣを通る線上での点Ｐからの距
離ｆを選ぶことによりこの虚像アレイ１０１は単一形態
となる。焦線Ｌｍは、点Ｐの虚像アレイ１０１に対する
垂直投影で定義される中心Ｏ′から行軸に沿って距離ｐ
ｃかつ列軸に沿って距離ｑｃの位置におけるこの虚像ア
レイ１０１まで伸長される。２つのアレイ１及び１０１
間の角度はα（ｐｍ）である。
【００２３】これにより、ｐｃ及びｐｍ間、並びにｑｃ
及びｑｍ間の関係を見い出すことができラドン変換の１
次微分を用いたアルゴリズムのように公知の反転アルゴ
リズムを適用することができる。なお、アレイ１によっ
てなされた測定結果を虚像アレイ１０１によってなされ
る測定結果に変換するための角度α（ｐｍ）と同様に、
ｑｍはアレイ１及び１０１（図４）の他方の座標に沿っ
た検出器２ｍの列座標である。略三角形状の存在により
次式が得られる。ここで、ｔは中心Ｃ及び線Ｌｍ間の距
離である。 ｔ＝ｐｃ｛Ｒ／√（ｆ2 ＋ｐｃ2 ）｝ （２）
及びｑｍ間の関係を見い出すことができラドン変換の１
次微分を用いたアルゴリズムのように公知の反転アルゴ
リズムを適用することができる。なお、アレイ１によっ
てなされた測定結果を虚像アレイ１０１によってなされ
る測定結果に変換するための角度α（ｐｍ）と同様に、
ｑｍはアレイ１及び１０１（図４）の他方の座標に沿っ
た検出器２ｍの列座標である。略三角形状の存在により
次式が得られる。ここで、ｔは中心Ｃ及び線Ｌｍ間の距
離である。 ｔ＝ｐｃ｛Ｒ／√（ｆ2 ＋ｐｃ2 ）｝ （２）
【００２４】ψｍ及びψｃが焦平面の固定された基準軸
に関してアレイ１及び１０１の向きを指定するのに用い
られるとすると、下式が定義できる。ここで、θは基準
軸及び焦線Ｌｍ上の円Ｃｅの中心の投影間の角度であ
る。 θ＝ψｃ−ａｒｃｔｇ（ｆ／ｐｃ） （３）
に関してアレイ１及び１０１の向きを指定するのに用い
られるとすると、下式が定義できる。ここで、θは基準
軸及び焦線Ｌｍ上の円Ｃｅの中心の投影間の角度であ
る。 θ＝ψｃ−ａｒｃｔｇ（ｆ／ｐｃ） （３）
【００２５】距離Ｄ（ｐｍ）を導くのに（２）式及び
（３）式は下記の（４）式及び（５）式と等価である。 ｔ＝ｐｍ［｛Ｄ（ｐｍ）−（ｆ−Ｒ）｝／√｛Ｄ（ｐｍ）2 ＋ｐｍ2 ｝］ （４） θ＝ψｍ−ａｒｃｔｇ｛Ｄ（ｐｍ）／ｐｍ｝ （５）
（３）式は下記の（４）式及び（５）式と等価である。 ｔ＝ｐｍ［｛Ｄ（ｐｍ）−（ｆ−Ｒ）｝／√｛Ｄ（ｐｍ）2 ＋ｐｍ2 ｝］ （４） θ＝ψｍ−ａｒｃｔｇ｛Ｄ（ｐｍ）／ｐｍ｝ （５）
【００２６】（２）式及び（４）式に基づいてｐｍ及び
ｐｃ間の関係が下式のごとくに求まる。 ｐｍ［｛Ｄ（ｐｍ）−（ｆ−Ｒ）｝／√｛Ｄ（ｐｍ）2
＋ｐｍ2 ｝］＝ｐｃ｛Ｒ／√（ｆ2 ＋ｐｃ2 ）｝ 次いでｐｍがｐｃの関数として計算される。
ｐｃ間の関係が下式のごとくに求まる。 ｐｍ［｛Ｄ（ｐｍ）−（ｆ−Ｒ）｝／√｛Ｄ（ｐｍ）2
＋ｐｍ2 ｝］＝ｐｃ｛Ｒ／√（ｆ2 ＋ｐｃ2 ）｝ 次いでｐｍがｐｃの関数として計算される。
【００２７】あらかじめ作成した対応テーブルを定義す
ることも可能である。
ることも可能である。
【００２８】（３）式及び（５）式から下式が得られ
る。 ψｍ＝ψｃ−ａｒｃｔｇ（ｆ／ｐｃ）＋ａｒｃｔｇ｛Ｄ
（ｐｍ）／ｐｍ｝ａｒｃｔｇ｛Ｄ（ｐｍ）／ｐｍ｝−ａ
ｒｃｔｇ（ｆ／ｐｃ）＝δψ（ｐｃ）
る。 ψｍ＝ψｃ−ａｒｃｔｇ（ｆ／ｐｃ）＋ａｒｃｔｇ｛Ｄ
（ｐｍ）／ｐｍ｝ａｒｃｔｇ｛Ｄ（ｐｍ）／ｐｍ｝−ａ
ｒｃｔｇ（ｆ／ｐｃ）＝δψ（ｐｃ）
【００２９】対応テーブルδψを設けることも可能であ
る。
る。
【００３０】ψｍ及びψｃ間の変数変化は、関数δψ
（ｐｃ）を用いて次のように定義される。 ψｍ＝ψｃ＋δψ（ｐｃ）
（ｐｃ）を用いて次のように定義される。 ψｍ＝ψｃ＋δψ（ｐｃ）
【００３１】図４は、まさに、 ｑｃ／ｑｍ＝ｆｃ／ｆｍ （６） を示している。ここで、ｆｃ及びｆｍは、点Ｐと実像ア
レイ１間及び点Ｐと虚像アレイ１０１間における焦線Ｌ
ｍの焦平面上への投影の長さをそれぞれ表わしている。
レイ１間及び点Ｐと虚像アレイ１０１間における焦線Ｌ
ｍの焦平面上への投影の長さをそれぞれ表わしている。
【００３２】ｆｃ及びｆｍは下記の（７）式及び（８）
式で与えられる。 ｆｍ＝｛ｐｍ−Ｒｓｉｎα（ｐｍ）｝√［１＋｛Ｄ（ｐｍ）／ｐｍ｝2 ］ （７） ｆｃ＝√（ｆ2 ＋ｐｃ2 ） （８）
式で与えられる。 ｆｍ＝｛ｐｍ−Ｒｓｉｎα（ｐｍ）｝√［１＋｛Ｄ（ｐｍ）／ｐｍ｝2 ］ （７） ｆｃ＝√（ｆ2 ＋ｐｃ2 ） （８）
【００３３】対応テーブルδｑが設けられる。ｑｍ及び
ｑｃ間の変数の変化は、関数δｑ（ｐｃ）をｑｍ＝ｑｃ
×δｑ（ｐｃ）及びδｑ（ｐｃ）＝ｆｃ／ｆｍのように
用いることによって確立される。ここで、ｑｍはｑｃの
関数として得られる。実際には、虚像測定結果ψｃ、ｐ
ｃ、ｑｃのボリュームに対応しており上述した式から得
られる実像測定結果ψｍ、ｐｍ、ｑｍのボリュームの座
標は、これら虚像測定結果ψｃ、ｐｃ、ｑｃのボリュー
ムの各点と統合していない。即ち、これらは測定点に対
応していない。点ｃ、ｐｃ及びｑｃにおける実像測定結
果のボリュームの値は、三線補間又は１つのディメンジ
ョンにおける３回の連続する直線補間によって決定され
る。これらの値は、全て、Ｄ（ｐｍ）、螺旋ピッチ、焦
点距離ｆ、検出器及び回転軸Ｒ間距離、並びにｐｍ、ψ
ｍ、ｑｍにそれぞれ応じたものとなる。実像測定結果
は、対象物４の回りを移動する虚像アレイ１０１によっ
てなされψｃ、ｐｃ及びｑｃによって特徴付けられる虚
像測定結果に関連付けられる。これら虚像測定結果は、
ヨーロッパ特許公開公報第２９２４０２号と同様の方法
によって組み合わされ、関数変換され、そして反転され
る。この公報には、組み合わせを行い、ラドン関数変換
の１次微分を反転し、結果を表示する基本的なデータ処
理、測定結果格納手段について簡単に記載されている。
ｑｃ間の変数の変化は、関数δｑ（ｐｃ）をｑｍ＝ｑｃ
×δｑ（ｐｃ）及びδｑ（ｐｃ）＝ｆｃ／ｆｍのように
用いることによって確立される。ここで、ｑｍはｑｃの
関数として得られる。実際には、虚像測定結果ψｃ、ｐ
ｃ、ｑｃのボリュームに対応しており上述した式から得
られる実像測定結果ψｍ、ｐｍ、ｑｍのボリュームの座
標は、これら虚像測定結果ψｃ、ｐｃ、ｑｃのボリュー
ムの各点と統合していない。即ち、これらは測定点に対
応していない。点ｃ、ｐｃ及びｑｃにおける実像測定結
果のボリュームの値は、三線補間又は１つのディメンジ
ョンにおける３回の連続する直線補間によって決定され
る。これらの値は、全て、Ｄ（ｐｍ）、螺旋ピッチ、焦
点距離ｆ、検出器及び回転軸Ｒ間距離、並びにｐｍ、ψ
ｍ、ｑｍにそれぞれ応じたものとなる。実像測定結果
は、対象物４の回りを移動する虚像アレイ１０１によっ
てなされψｃ、ｐｃ及びｑｃによって特徴付けられる虚
像測定結果に関連付けられる。これら虚像測定結果は、
ヨーロッパ特許公開公報第２９２４０２号と同様の方法
によって組み合わされ、関数変換され、そして反転され
る。この公報には、組み合わせを行い、ラドン関数変換
の１次微分を反転し、結果を表示する基本的なデータ処
理、測定結果格納手段について簡単に記載されている。
【００３４】実際には、実像多焦点測定結果から虚像円
錐測定結果に至ることを可能とする上述した式又は式の
他の任意の組を使用する。すべての場合にＤ（ｐｍ）を
知る必要があり、図１及び図２に示す焦点合わせに関す
る不変の幾何学的特性に基づいて、ｐｍの関数であるＤ
（ｐｍ）のテーブル作成を前もって行うことが可能であ
る。ｐｃ及び特にｐｍに応じた値となるＤ（ｐｍ）から
ｐｍを求めることは、常に直ちに行えるものではない。
従って、各測定結果について計算を繰り返すよりも、ｐ
ｃの関数としてのｐｍ、δψ及びδｑのテーブルを前も
って作成しておくことが好ましい。実際には、これらテ
ーブルはディメンジョンＮ3 で再構築されるべき対象物
のために作成され、再構築段階が前述のヨーロッパ特許
公開公報第２９２４０２号に記載されているようにほぼ
Ｎの基本的な動作を必要とするのに対して、再配列段階
はほぼＮ3 の三線補間を必要とする。
錐測定結果に至ることを可能とする上述した式又は式の
他の任意の組を使用する。すべての場合にＤ（ｐｍ）を
知る必要があり、図１及び図２に示す焦点合わせに関す
る不変の幾何学的特性に基づいて、ｐｍの関数であるＤ
（ｐｍ）のテーブル作成を前もって行うことが可能であ
る。ｐｃ及び特にｐｍに応じた値となるＤ（ｐｍ）から
ｐｍを求めることは、常に直ちに行えるものではない。
従って、各測定結果について計算を繰り返すよりも、ｐ
ｃの関数としてのｐｍ、δψ及びδｑのテーブルを前も
って作成しておくことが好ましい。実際には、これらテ
ーブルはディメンジョンＮ3 で再構築されるべき対象物
のために作成され、再構築段階が前述のヨーロッパ特許
公開公報第２９２４０２号に記載されているようにほぼ
Ｎの基本的な動作を必要とするのに対して、再配列段階
はほぼＮ3 の三線補間を必要とする。
【００３５】実像多焦点測定結果を虚像円錐測定結果に
変換する段階において、テーブルｐｍ、δψ及びδｑ
は、変数ｐｍ、ｑｍ及びψｍの変化に対してｐｃ、ｑｃ
及びψｃの関数として後に使用できる。
変換する段階において、テーブルｐｍ、δψ及びδｑ
は、変数ｐｍ、ｑｍ及びψｍの変化に対してｐｃ、ｑｃ
及びψｃの関数として後に使用できる。
【００３６】これらの変換に起因して処理が直ちに複雑
となるものでもなく、また、一般に多くの適用例におい
て、処理が複雑となることよりもより良い像を得ること
の方が重要である。
となるものでもなく、また、一般に多くの適用例におい
て、処理が複雑となることよりもより良い像を得ること
の方が重要である。
【００３７】特定の装置においては、幅４５０ｍｍ、高
さ３５０ｍｍ（Ｘ軸方向において）の矩形の検出器アレ
イを覆うコリメータを設けることができる。これは、図
２及びｋ＝０．０１６の場合の（１）式の状況に係るこ
とを可能にしている。Ｒ＝３５０ｍｍ及びｆ＝６００ｍ
ｍであるとき、Ｄ（ｐｍ）は６００ｍｍから１４１０ｍ
ｍへ変化する。検出器アレイは患者の回りを回転してデ
ータ獲得を可能とするガンマカメラに属するものであ
る。例えばこの種の装置は、ソファメディカル（Ｓｏｐ
ｈａ Ｍｅｄｉｃａｌ）ＤＳＸカメラの機械的部分から
構成される。
さ３５０ｍｍ（Ｘ軸方向において）の矩形の検出器アレ
イを覆うコリメータを設けることができる。これは、図
２及びｋ＝０．０１６の場合の（１）式の状況に係るこ
とを可能にしている。Ｒ＝３５０ｍｍ及びｆ＝６００ｍ
ｍであるとき、Ｄ（ｐｍ）は６００ｍｍから１４１０ｍ
ｍへ変化する。検出器アレイは患者の回りを回転してデ
ータ獲得を可能とするガンマカメラに属するものであ
る。例えばこの種の装置は、ソファメディカル（Ｓｏｐ
ｈａ Ｍｅｄｉｃａｌ）ＤＳＸカメラの機械的部分から
構成される。
【００３８】螺旋軌道によってデータ獲得を行うには、
検出器アレイは、回転移動に加えて患者の一部分を回転
軸に沿って並進移動する。本発明に関する説明を簡易化
するために、以下の式は平面検出器によって取得された
測定結果について与えられるものとする。しかしなが
ら、円形の検出器も使用可能であり、この場合、検出器
の２つの端を結ぶ線と検出器の幅とによって規定される
平面上の獲得線の分布によって測定が行われる。
検出器アレイは、回転移動に加えて患者の一部分を回転
軸に沿って並進移動する。本発明に関する説明を簡易化
するために、以下の式は平面検出器によって取得された
測定結果について与えられるものとする。しかしなが
ら、円形の検出器も使用可能であり、この場合、検出器
の２つの端を結ぶ線と検出器の幅とによって規定される
平面上の獲得線の分布によって測定が行われる。
【００３９】像再構築アルゴリズムは、前に述べたよう
に、グランゲート（Ｇｒａｎｇｅａｔ）アルゴリズムと
呼ばれるラドン変換の１次微分の反転であることが有利
である。
に、グランゲート（Ｇｒａｎｇｅａｔ）アルゴリズムと
呼ばれるラドン変換の１次微分の反転であることが有利
である。
【００４０】焦点Ｐが、有利なことには、アレイ１の回
転軸Ｘに属した中心Ｃを有する円Ｃｅ上に位置している
場合、変化が通常は許容できる品質である像の鮮鋭度を
劣化させるので、前述の条件は必須ではなくなる。この
ため、焦点は、装置の残りの部分におけるアレイ１の位
置及び装着不安定性によって、その中心が回転軸Ｘに一
致せずその近傍に存在している円の上にも位置し得る。
図５には、図２の場合に類似しているが焦点Ｐがアレイ
１に平行でありかつＸ軸に垂直な線Ｙ上に位置している
という異なる状況が示されている。この場合、焦点Ｐが
円上、例えば円Ｃｅ上、に位置していると仮定して近似
が行われる。この条件で得られる投影点はＰ′で示され
る。これは、特に、円Ｃｅを含む平面への各焦線ＬのＸ
軸に平行な投影と円Ｃｅとの交差点に相当する。焦点Ｐ
は、像の中心における投影点Ｐ′に虚像的に一致し、こ
れによって再構築された３次元像は、対象物４をカバー
するこの領域においては像の周辺部よりも鮮鋭度が高く
なる。
転軸Ｘに属した中心Ｃを有する円Ｃｅ上に位置している
場合、変化が通常は許容できる品質である像の鮮鋭度を
劣化させるので、前述の条件は必須ではなくなる。この
ため、焦点は、装置の残りの部分におけるアレイ１の位
置及び装着不安定性によって、その中心が回転軸Ｘに一
致せずその近傍に存在している円の上にも位置し得る。
図５には、図２の場合に類似しているが焦点Ｐがアレイ
１に平行でありかつＸ軸に垂直な線Ｙ上に位置している
という異なる状況が示されている。この場合、焦点Ｐが
円上、例えば円Ｃｅ上、に位置していると仮定して近似
が行われる。この条件で得られる投影点はＰ′で示され
る。これは、特に、円Ｃｅを含む平面への各焦線ＬのＸ
軸に平行な投影と円Ｃｅとの交差点に相当する。焦点Ｐ
は、像の中心における投影点Ｐ′に虚像的に一致し、こ
れによって再構築された３次元像は、対象物４をカバー
するこの領域においては像の周辺部よりも鮮鋭度が高く
なる。
【００４１】上述した処理は、減衰補正動作と同様に、
検出器２の均一性及び感度の補正等の前処理によって実
施することができる。この処理は、サウンダーズによっ
て出版されたジェイ・エイ・ソレンソン及びエム・イー
・フェルプスによる「核医学における物理」なる本に記
載されている従来の再構築においても、使用可能であ
る。
検出器２の均一性及び感度の補正等の前処理によって実
施することができる。この処理は、サウンダーズによっ
て出版されたジェイ・エイ・ソレンソン及びエム・イー
・フェルプスによる「核医学における物理」なる本に記
載されている従来の再構築においても、使用可能であ
る。
【００４２】最も単純な場合について説明をしてきた
が、本発明は検出器アレイが対象物４の回りを徐々に上
昇して行く螺旋軌道についても適用することができる。
図６には、軌道Ｔがその各ターン間のピッチｈを有する
螺旋軌道の例が示されている。この螺旋は、Ｘ軸を中心
としている。コリメータ３の焦点Ｐの位置は、図６にそ
の一部を表わした、同一ピッチｈを有する他の螺旋Ｈを
描いて移動する。この螺旋ＨもＸ軸に中心を有してい
る。この螺旋軌道における像再構築方法は、各焦点Ｐが
異なる列の検出器２に係り合い得るので、前述した方法
と同じままであることは明らかである。各照射の瞬間の
関数としての像再構築について、ピッチｈのみが計算の
最終段階で影響を及ぼす。円軌道及び円Ｃｅは、螺旋軌
道においてそのピッチがゼロの場合に対応している。
が、本発明は検出器アレイが対象物４の回りを徐々に上
昇して行く螺旋軌道についても適用することができる。
図６には、軌道Ｔがその各ターン間のピッチｈを有する
螺旋軌道の例が示されている。この螺旋は、Ｘ軸を中心
としている。コリメータ３の焦点Ｐの位置は、図６にそ
の一部を表わした、同一ピッチｈを有する他の螺旋Ｈを
描いて移動する。この螺旋ＨもＸ軸に中心を有してい
る。この螺旋軌道における像再構築方法は、各焦点Ｐが
異なる列の検出器２に係り合い得るので、前述した方法
と同じままであることは明らかである。各照射の瞬間の
関数としての像再構築について、ピッチｈのみが計算の
最終段階で影響を及ぼす。円軌道及び円Ｃｅは、螺旋軌
道においてそのピッチがゼロの場合に対応している。
【００４３】このように本発明によれば、検出器２の平
面アレイ１によって対象物４の３次元像が再構築され
る。ここで、検出器の焦点は円Ｃｅ上に位置しているか
ら、単一点に焦点合わせする場合よりもより広い検査領
域を得ることができる。単一点への円錐焦点合わせによ
る特性測定結果を反転するアルゴリズムが用いられる。
本発明は、医学や産業上の非破壊検査に適用される。
面アレイ１によって対象物４の３次元像が再構築され
る。ここで、検出器の焦点は円Ｃｅ上に位置しているか
ら、単一点に焦点合わせする場合よりもより広い検査領
域を得ることができる。単一点への円錐焦点合わせによ
る特性測定結果を反転するアルゴリズムが用いられる。
本発明は、医学や産業上の非破壊検査に適用される。
【図１】本発明の一般的な構成を示す図である。
【図２】本発明の好ましい実施例における焦線の分布配
置を軸に沿って示した図である。
置を軸に沿って示した図である。
【図３】円錐焦点合わせ法に基づく像再構築に用いられ
る方法を説明する図である。
る方法を説明する図である。
【図４】円錐焦点合わせ法に基づく像再構築に用いられ
る方法を説明する図である。
る方法を説明する図である。
【図５】円上に焦点合わせする方法と僅かに異なる条件
で本発明がどのように用いられるかを説明するための図
である。
で本発明がどのように用いられるかを説明するための図
である。
【図６】本発明の方法を一般化して表わした図である。
１ アレイ ２ 検出器 ３ コリメータ ４ 対象物
フロントページの続き (51)Int.Cl.6 識別記号 庁内整理番号 ＦＩ 技術表示箇所 // Ｇ０２Ｂ 27/22 9120−2Ｋ
Claims (11)
- 【請求項１】 対象物（４）に向かって配設された検出
器（２）の２次元アレイ（１）を備えており、該検出器
アレイが螺旋軌道（Ｔ）に沿って軸（Ｘ）の回りを回転
しつつ前記対象物に対して移動するべく構成されてい
る、前記対象物（４）の３次元像を再構築する装置であ
って、前記検出器（２）は、前記軸（Ｘ）に中心を有す
る螺旋（Ｈ）上にほぼ位置している複数の点（Ｐ）上に
軸方向に焦点合わせされており、前記検出器によってな
された測定結果を格納しかつ該測定結果を組み合わせる
ための本質的に公知の手段が像再構築に用いられている
ことを特徴とする３次元像を再構築する装置。 - 【請求項２】 前記軸に平行な線上に位置するアレイの
前記検出器が全て単一の焦点に焦点合わせされているこ
とを特徴とする請求項１に記載の装置。 - 【請求項３】 前記螺旋軌道が円形軌道であり前記螺旋
が円（Ｃｅ）であることを特徴とする請求項１に記載の
装置。 - 【請求項４】 前記焦点（Ｐ）の位置する円（Ｃｅ）の
平面上において、前記円の中心（Ｃ）を通る前記検出器
アレイへの垂線（Ｆ）上での前記検出器アレイと前記焦
線（Ｌｍ）及び該垂線（Ｆ）の交点との間の距離（Ｄ
（ｐｍ））が前記円（Ｃｅ）の半径より大きいことを特
徴とする請求項３に記載の装置。 - 【請求項５】 前記距離が、前記平面内において前記垂
線とそれから焦線がそれぞれ生じる検出器（２ｍ）との
間で測定される距離（ｐｍ）の増加関数であることを特
徴とする請求項４に記載の装置。 - 【請求項６】 前記増加関数が２次関数であることを特
徴とする請求項５に記載の装置。 - 【請求項７】 螺旋軌道に従って軸の回りを回転しつつ
対象物（４）に対して移動する検出器（２）の少なくと
も１つの２次元アレイ（１）によって前記対象物（４）
の２次元像を得ることを含んでおり、前記２次元像がコ
リメータ線（Ｌ）に沿った前記対象物（４）の特性に関
する前記各検出器（２）の測定結果から構成されてお
り、該測定結果を用いて式を解くことによって３次元像
を再構築する方法であって、前記各コリメータ線が前記
軸を中心とする螺旋上をほぼ通り、前記方法が前記式を
解く前に行われる、前記コリメータ線をパラメータで表
わす変換演算段階と再編制段階とを含んでおり、該各コ
リメータ線が前記螺旋の単一の焦点上に焦点合わせされ
た前記検出器の虚像２次元アレイ（１０１）の始点と見
なされており、該虚像アレイが前記軸（Ｘ）についてこ
れから一定距離で回転するものであり、前記変換演算が
前記虚像アレイの前記各コリメータ線について始点とな
る実像アレイ内における位置及び該実像アレイの前記軸
に対する角度位置を決定するものであり、次いで、虚像
円錐投影を得るために前記再編制段階において該変換演
算の結果に応じて前記測定結果を再編制することからな
り、前記解が該再編制された測定結果に関連しているこ
とを特徴とする３次元像を再構築する方法。 - 【請求項８】 前記測定結果について解く方法がラドン
変換の１次微分の反転を含むことを特徴とする請求項７
に記載の３次元像を再構築する方法。 - 【請求項９】 前記螺旋軌道が円形軌道であり前記螺旋
が円（Ｃｅ）であることを特徴とする請求項７に記載の
３次元像を再構築する方法。 - 【請求項１０】 全ての前記焦線が前記円の接線を通る
ことを特徴とする請求項９に記載の３次元像を再構築す
る方法。 - 【請求項１１】 前記変換演算が３つの座標（ｐ、ｑ、
ψ）による虚像測定結果の３つの連続する変換段階を備
えており、該変換段階の各々が直線補間を含んでいるこ
とを特徴とする請求項７に記載の３次元像を再構築する
方法。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR9306564 | 1993-06-02 | ||
| FR9306564A FR2706043B1 (fr) | 1993-06-02 | 1993-06-02 | Installation et procédé de reconstruction d'images tridimensionnelles. |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH0771941A true JPH0771941A (ja) | 1995-03-17 |
Family
ID=9447654
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP6142197A Withdrawn JPH0771941A (ja) | 1993-06-02 | 1994-06-02 | ３次元像を再構築する装置及び方法 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US5933517A (ja) |
| EP (1) | EP0627702A1 (ja) |
| JP (1) | JPH0771941A (ja) |
| CA (1) | CA2124905A1 (ja) |
| FR (1) | FR2706043B1 (ja) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6160924A (en) * | 1998-08-12 | 2000-12-12 | Northrop Grumman Corporation | Method for forming a map of a three-dimensional object |
| US6944259B2 (en) * | 2001-09-26 | 2005-09-13 | Massachusetts Institute Of Technology | Versatile cone-beam imaging apparatus and method |
| US6865246B2 (en) * | 2001-09-26 | 2005-03-08 | Massachusetts Institute Of Technology | True 3D cone-beam imaging method and apparatus |
| US7847274B2 (en) * | 2004-08-12 | 2010-12-07 | Navotek Medical Ltd. | Localization of a radioactive source within a body of a subject |
| WO2006043276A2 (en) * | 2004-10-19 | 2006-04-27 | Navotek Medical Ltd. | Locating a catheter tip using a tracked guide |
| US8239002B2 (en) * | 2004-08-12 | 2012-08-07 | Novatek Medical Ltd. | Guiding a tool for medical treatment by detecting a source of radioactivity |
| WO2007017846A2 (en) * | 2005-08-11 | 2007-02-15 | Navotek Medical Ltd. | Localization of a radioactive source |
| ATE555737T1 (de) * | 2005-08-11 | 2012-05-15 | Navotek Medical Ltd | Lokalisation einer radioaktiven quelle |
| JP5539719B2 (ja) * | 2006-08-31 | 2014-07-02 | コーニンクレッカ フィリップス エヌ ヴェ | 画像形成システム |
| US7942843B2 (en) * | 2008-08-18 | 2011-05-17 | Navotek Medical Ltd. | Implantation device for soft tissue markers and other implants |
| CN103047946B (zh) * | 2012-12-10 | 2015-11-04 | 北京航空航天大学 | 一种基于虚拟数据的轴对称参数分布图像重建方法 |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2529949B2 (ja) * | 1986-08-12 | 1996-09-04 | 株式会社東芝 | 同期画像再構成装置 |
| IL80364A (en) * | 1986-10-20 | 1990-03-19 | Elscint Ltd | Three dimensional image construction using binary space interpolation |
| US4823017A (en) * | 1986-12-19 | 1989-04-18 | Siemens Gammasonics, Inc. | Scintillation camera and multifocal fan-beam collimator used therein |
| US4820924A (en) * | 1986-12-19 | 1989-04-11 | Siemens Gammasonics, Inc. | Scintillation camera and three dimensional multifocal collimator used therewith |
| US5150427A (en) * | 1989-09-29 | 1992-09-22 | General Electric Company | Three dimensional disarticulation |
| 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. |
| US5210421A (en) * | 1991-06-10 | 1993-05-11 | Picker International, Inc. | Simultaneous transmission and emission converging tomography |
| US5404293A (en) * | 1991-06-11 | 1995-04-04 | The University Of Utah | Cone beam reconstruction using helical data collection paths |
| US5170439A (en) * | 1991-06-11 | 1992-12-08 | Picker International, Inc. | Cone beam reconstruction using combined circle and line orbits |
| JPH0528243A (ja) * | 1991-07-24 | 1993-02-05 | Toshiba Corp | 画像作成装置 |
-
1993
- 1993-06-02 FR FR9306564A patent/FR2706043B1/fr not_active Expired - Fee Related
-
1994
- 1994-05-30 EP EP94401185A patent/EP0627702A1/fr not_active Withdrawn
- 1994-06-01 CA CA002124905A patent/CA2124905A1/en not_active Abandoned
- 1994-06-02 JP JP6142197A patent/JPH0771941A/ja not_active Withdrawn
-
1997
- 1997-02-07 US US08/797,266 patent/US5933517A/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CA2124905A1 (en) | 1994-12-03 |
| US5933517A (en) | 1999-08-03 |
| EP0627702A1 (fr) | 1994-12-07 |
| FR2706043A1 (fr) | 1994-12-09 |
| FR2706043B1 (fr) | 1995-07-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5838756A (en) | Radiation computed tomography apparatus | |
| EP0500857B1 (en) | Cone beam scanning trajectories for three-dimensional computerized tomography data acquisition where object is larger than the field of view | |
| US7106825B2 (en) | Apparatus and method for reconstruction of volumetric images in a divergent scanning computed tomography system | |
| US6256365B1 (en) | Apparatus and method for reconstruction of images in a computed tomography system using oblique slices | |
| JP2914891B2 (ja) | Ｘ線コンピュータ断層撮影装置 | |
| US5319693A (en) | Three dimensional computerized tomography scanning configuration for imaging large objects with smaller area detectors | |
| US5375156A (en) | Method and apparatus for 3-D computer tomography | |
| JP5133690B2 (ja) | ボクセルに依存する補間を用いる画像再構成 | |
| US8094910B2 (en) | Method of reconstructing an image function from Radon data | |
| US8284892B2 (en) | System and method for image reconstruction | |
| CN1688254B (zh) | X线断层摄影装置 | |
| US5396528A (en) | Tomographic image reconstruction using cross-plane rays | |
| US7366277B2 (en) | Tomographic device and method therefor | |
| JPH0442011B2 (ja) | ||
| Karolczak et al. | Implementation of a cone‐beam reconstruction algorithm for the single‐circle source orbit with embedded misalignment correction using homogeneous coordinates | |
| US6574298B2 (en) | Cone beam scanner using oblique surface reconstructions | |
| IL96050A (en) | A method of reconstructing a computed tomographic image for the local scan | |
| JPH06214031A (ja) | イメージ表示装置およびイメージ表示方法 | |
| JP5180584B2 (ja) | 螺旋相対運動及び円錐ビームを用いたコンピュータ断層撮影方法 | |
| JP2007512034A (ja) | 発散ビームスキャナのための画像再構成方法 | |
| KR20040072466A (ko) | 3차원 백 프로젝션 방법 및 ｘ-레이 ｃｔ 장치 | |
| JPH0771941A (ja) | ３次元像を再構築する装置及び方法 | |
| US6885764B2 (en) | High Speed Z-smoothing method and apparatus for CT imaging system | |
| JP3913798B2 (ja) | 画像再構成処理装置 | |
| US7020317B2 (en) | Computed tomography apparatus |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A300 | Application deemed to be withdrawn because no request for examination was validly filed |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20010904 |