JPH11328358A - 医用画像デ―タを処理する方法及び装置 - Google Patents
医用画像デ―タを処理する方法及び装置Info
- Publication number
- JPH11328358A JPH11328358A JP11081993A JP8199399A JPH11328358A JP H11328358 A JPH11328358 A JP H11328358A JP 11081993 A JP11081993 A JP 11081993A JP 8199399 A JP8199399 A JP 8199399A JP H11328358 A JPH11328358 A JP H11328358A
- Authority
- JP
- Japan
- Prior art keywords
- data
- function data
- input function
- projection
- grid
- 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
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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S378/00—X-ray or gamma ray systems or devices
- Y10S378/901—Computer tomography program or processor
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 Processing (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Analysis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
(57)【要約】
【課題】 より減少された時間で所与の関数データから
任意の投影を計算しうる医用画像データを処理する方法
及び装置を提供することを目的とする。 【解決手段】 本発明による方法は、（ａ）上記入力関
数データ（Ｉc （ｒ））を重畳カーネル（Ｃ（ｋ））の
フーリエ変換（ｃ（ｒ））によって割り算することによ
って変更された入力関数データ（Ｉc ’（ｒ））を決定
する段階と、（ｂ）上記変更された入力関数データ（Ｉ
c ’（ｒ））のフーリエ変換（Ｓc ’（ｋ））を決定す
る段階と、（ｃ）上記変更された入力関数データのフー
リエ変換（Ｓc ’（ｋ））を重畳カーネル（Ｃ（ｋ））
と共に重畳することによって上記出力関数データ（Ｓ
（ｋ））を決定する段階とを含み、上記重畳カーネル
（Ｃ（ｋ））は、上記選択された格子の格子点（８）か
らの距離に従って上記デカルト格子（１，６）の格子点
（２，７）を重み付けする重み関数である。
任意の投影を計算しうる医用画像データを処理する方法
及び装置を提供することを目的とする。 【解決手段】 本発明による方法は、（ａ）上記入力関
数データ（Ｉc （ｒ））を重畳カーネル（Ｃ（ｋ））の
フーリエ変換（ｃ（ｒ））によって割り算することによ
って変更された入力関数データ（Ｉc ’（ｒ））を決定
する段階と、（ｂ）上記変更された入力関数データ（Ｉ
c ’（ｒ））のフーリエ変換（Ｓc ’（ｋ））を決定す
る段階と、（ｃ）上記変更された入力関数データのフー
リエ変換（Ｓc ’（ｋ））を重畳カーネル（Ｃ（ｋ））
と共に重畳することによって上記出力関数データ（Ｓ
（ｋ））を決定する段階とを含み、上記重畳カーネル
（Ｃ（ｋ））は、上記選択された格子の格子点（８）か
らの距離に従って上記デカルト格子（１，６）の格子点
（２，７）を重み付けする重み関数である。
Description
【０００１】
【発明の属する技術分野】本発明は、デカルト格子の格
子点に配置された入力関数データから選択された格子の
格子点に配置された出力関数データを決定する医用画像
データ処理方法、及び係る方法を実行しうる装置に関す
る。
子点に配置された入力関数データから選択された格子の
格子点に配置された出力関数データを決定する医用画像
データ処理方法、及び係る方法を実行しうる装置に関す
る。
【０００２】
【従来の技術】上述の種類の方法は、医用画像データが
処理されうる多くの医用適用において必要とされる。例
えば、２次元／３次元記録のためには、３次元データセ
ットから２次元投影を計算する効率的な方法を有するこ
とが重要である。現在、このために、例えば（Ｘ線撮像
の場合）放射源から検出器平面中の像点までの光線に沿
って手術前２次元データセット中のグレーレベルの積分
を計算する光線追跡アルゴリズムが使用されうる。係る
方法の手術間使用を可能にするために、手術前データセ
ットの複数の投影は異なる角度で計算され２次元／３次
元記録が形成されることを可能にするために現在測定さ
れている投影と比較されねばならないため、計算は充分
に迅速に行われねばならない。上述の種類の方法はＸ線
撮像に使用されるだけでなく、例えば磁気共鳴断層撮影
法、Ｘ線コンピュータ断層撮影法（ＣＴ）、特に投影の
計算といった他の医用撮像技術のために使用される。
処理されうる多くの医用適用において必要とされる。例
えば、２次元／３次元記録のためには、３次元データセ
ットから２次元投影を計算する効率的な方法を有するこ
とが重要である。現在、このために、例えば（Ｘ線撮像
の場合）放射源から検出器平面中の像点までの光線に沿
って手術前２次元データセット中のグレーレベルの積分
を計算する光線追跡アルゴリズムが使用されうる。係る
方法の手術間使用を可能にするために、手術前データセ
ットの複数の投影は異なる角度で計算され２次元／３次
元記録が形成されることを可能にするために現在測定さ
れている投影と比較されねばならないため、計算は充分
に迅速に行われねばならない。上述の種類の方法はＸ線
撮像に使用されるだけでなく、例えば磁気共鳴断層撮影
法、Ｘ線コンピュータ断層撮影法（ＣＴ）、特に投影の
計算といった他の医用撮像技術のために使用される。
【０００３】
【発明が解決しようとする課題】本発明は、特に既知の
方法よりも高速な上述の種類の医用画像データを処理す
る方法を提供することを目的とする。
方法よりも高速な上述の種類の医用画像データを処理す
る方法を提供することを目的とする。
【０００４】
【課題を解決するための手段】この目的は、請求項１記
載の方法によって達成される。本発明は、時間のかかる
補間の計算は、より簡単な計算ステップ（いわゆるリサ
ンプリング）によって置き換えられうる。デカルト格子
の格子点に存在する例えば空間領域中の２次元対象関数
データといった入力関数データは、まず重畳カーネルの
フーリエ変換によって割り算され、変更された入力関数
データが決定される。重畳カーネルは一種の重み関数で
あり、それによりデカルト格子の格子点は、選択された
格子の格子点からの距離に従って重み付けされる。重畳
カーネルは選択された格子の格子点における続く重畳の
ために所望とされる。即ちこの段階において変更された
入力関数データのフーリエ変換は、重畳カーネル自体の
中で重畳され、例えば周波数領域の中の極格子といった
選択された格子上に、例えば投影データといった出力関
数データを生成する。
載の方法によって達成される。本発明は、時間のかかる
補間の計算は、より簡単な計算ステップ（いわゆるリサ
ンプリング）によって置き換えられうる。デカルト格子
の格子点に存在する例えば空間領域中の２次元対象関数
データといった入力関数データは、まず重畳カーネルの
フーリエ変換によって割り算され、変更された入力関数
データが決定される。重畳カーネルは一種の重み関数で
あり、それによりデカルト格子の格子点は、選択された
格子の格子点からの距離に従って重み付けされる。重畳
カーネルは選択された格子の格子点における続く重畳の
ために所望とされる。即ちこの段階において変更された
入力関数データのフーリエ変換は、重畳カーネル自体の
中で重畳され、例えば周波数領域の中の極格子といった
選択された格子上に、例えば投影データといった出力関
数データを生成する。
【０００５】本発明の主な利点は、短い計算時間にあ
る。即ち、１次元投影を計算するために必要とされる計
算ステップの数は、ｎｌｏｇｎ（ｎ＝サンプリング点の
数）に比例し、２次元投影ではｎ2 ｌｏｇに比例し、一
方慣習的な方法では計算ステップの数は１次元投影に対
してはｎ2 に比例し、２次元投影に対してはｎ3 に比例
する。更に、本発明による方法の様々な適用において、
画質は向上されうる。
る。即ち、１次元投影を計算するために必要とされる計
算ステップの数は、ｎｌｏｇｎ（ｎ＝サンプリング点の
数）に比例し、２次元投影ではｎ2 ｌｏｇに比例し、一
方慣習的な方法では計算ステップの数は１次元投影に対
してはｎ2 に比例し、２次元投影に対してはｎ3 に比例
する。更に、本発明による方法の様々な適用において、
画質は向上されうる。
【０００６】本発明による方法は極格子上の投影データ
を計算するために使用されることが望ましい。選択され
た格子の格子点は、周波数領域の原点を通る放射状の線
又は平面上に位置する。画像関数データを形成するため
に、空間領域もまた逆フーリエ変換によって投影データ
から簡単に得られる。請求項３乃至６は、本発明の方法
の望ましい適用を開示する。当該の適用に依存して、出
力関数データが配置されるべき適当な格子が選択され
る。
を計算するために使用されることが望ましい。選択され
た格子の格子点は、周波数領域の原点を通る放射状の線
又は平面上に位置する。画像関数データを形成するため
に、空間領域もまた逆フーリエ変換によって投影データ
から簡単に得られる。請求項３乃至６は、本発明の方法
の望ましい適用を開示する。当該の適用に依存して、出
力関数データが配置されるべき適当な格子が選択され
る。
【０００７】望ましくは、重畳カーネルとしてカイザー
ベッセル関数又はガウス関数が使用されるが、他のカー
ネルの使用もまた可能である。本発明による方法の一連
のステップは、本発明の結果に重大な影響を与えること
なく逆にされうる。
ベッセル関数又はガウス関数が使用されるが、他のカー
ネルの使用もまた可能である。本発明による方法の一連
のステップは、本発明の結果に重大な影響を与えること
なく逆にされうる。
【０００８】
【発明の実施の形態】以下、本発明を図面を参照して詳
述する。図１は、デカルト格子１を示す図であり、格子
点２に、例えばＸ線画像の２次元対象関数データといっ
た入力関数データＩc （ｒ）が配置されている。以下、
添え字「ｃ」は、デカルト格子上のデータが考慮されて
いることを示すものとする。平行ビーム投影は入力関数
データＩc （ｒ）から形成されるべきであり、即ち強度
は平行線５に沿った直線４上に投影されるべきである。
係る投影線４に沿った投影値３は投影を形成する。
述する。図１は、デカルト格子１を示す図であり、格子
点２に、例えばＸ線画像の２次元対象関数データといっ
た入力関数データＩc （ｒ）が配置されている。以下、
添え字「ｃ」は、デカルト格子上のデータが考慮されて
いることを示すものとする。平行ビーム投影は入力関数
データＩc （ｒ）から形成されるべきであり、即ち強度
は平行線５に沿った直線４上に投影されるべきである。
係る投影線４に沿った投影値３は投影を形成する。
【０００９】図１から明らかであるように、平行光線５
は概して入力関数データＩc （ｒ）が存在する格子点２
にぶつからない。平行光線５上に存在する関数データを
決定するために、従って、使用可能な入力関数データＩ
c （ｒ）の間で補間が実行されるべきである。周波数領
域中のフーリエ変換は入力関数データを再びデカルト格
子の格子点７へ移動させる（図２参照）。周波数領域
（ｋ空間）の中の原点Ｏを通る放射状の線９に対応する
投影を形成するために、従って、放射状の線９に沿った
サンプリング点８においてデータを獲得するために、格
子点７に存在するデータの補間は再び実行されねばなら
ない。
は概して入力関数データＩc （ｒ）が存在する格子点２
にぶつからない。平行光線５上に存在する関数データを
決定するために、従って、使用可能な入力関数データＩ
c （ｒ）の間で補間が実行されるべきである。周波数領
域中のフーリエ変換は入力関数データを再びデカルト格
子の格子点７へ移動させる（図２参照）。周波数領域
（ｋ空間）の中の原点Ｏを通る放射状の線９に対応する
投影を形成するために、従って、放射状の線９に沿った
サンプリング点８においてデータを獲得するために、格
子点７に存在するデータの補間は再び実行されねばなら
ない。
【００１０】図３は、空間領域中の対象関数データＩc
（ｒ）から周波数領域中の投影データＳ（ｋ）を決定す
るための本発明による方法のフローチャートを示す図で
ある。第１のステップ１０において、対象関数データＩ
c （ｒ）は重畳カーネルＣ（ｋ）のフーリエ変換ｃ
（ｒ）によって割り算され、変更された対象関数データ
Ｉ c ’（ｒ）を生成する。次のステップ１１において、
そのフーリエ変換Ｓc ’（ｋ）が決定される。このフー
リエ変換は続いて、ステップ１２において選択された格
子の格子点において重畳カーネルＣ（ｋ）によって重畳
され、周波数領域中の投影データＳ（ｋ）を生成する。
このデータはもはやデカルト格子ではなく選択された格
子上に配置されている。このデータは逆フーリエ変換に
よって空間領域の中の投影データへ変換されうる。
（ｒ）から周波数領域中の投影データＳ（ｋ）を決定す
るための本発明による方法のフローチャートを示す図で
ある。第１のステップ１０において、対象関数データＩ
c （ｒ）は重畳カーネルＣ（ｋ）のフーリエ変換ｃ
（ｒ）によって割り算され、変更された対象関数データ
Ｉ c ’（ｒ）を生成する。次のステップ１１において、
そのフーリエ変換Ｓc ’（ｋ）が決定される。このフー
リエ変換は続いて、ステップ１２において選択された格
子の格子点において重畳カーネルＣ（ｋ）によって重畳
され、周波数領域中の投影データＳ（ｋ）を生成する。
このデータはもはやデカルト格子ではなく選択された格
子上に配置されている。このデータは逆フーリエ変換に
よって空間領域の中の投影データへ変換されうる。
【００１１】上述のステップは、数学的には以下のよう
に記述されうる。 Ｉc ’（ｒ）＝Ｉc （ｒ）／ｃ（ｒ） (a) Ｓc ’（ｋ）＝ＦＴ｛Ｉc ’（ｒ）｝＝Ｓc （ｋ）＊１／Ｃ（ｋ） (b) Ｓ（ｋ）＝Ｓc ’（ｋ）＊Ｃ（ｋ）＝（Ｓc （ｋ）＊１／Ｃ（ｋ））＊Ｃ（ｋ ） (c) 方法は、概してｎ次元対象関数データから（ｎ−１）次
元投影を決定するために使用されうる。このため、夫々
のｎ次元フーリエ変換はステップ１１において実行され
ねばならない。ステップ１２は、しばしば「重畳補間」
を使用する「リサンプリング」又は「グリッディング」
と称される。
に記述されうる。 Ｉc ’（ｒ）＝Ｉc （ｒ）／ｃ（ｒ） (a) Ｓc ’（ｋ）＝ＦＴ｛Ｉc ’（ｒ）｝＝Ｓc （ｋ）＊１／Ｃ（ｋ） (b) Ｓ（ｋ）＝Ｓc ’（ｋ）＊Ｃ（ｋ）＝（Ｓc （ｋ）＊１／Ｃ（ｋ））＊Ｃ（ｋ ） (c) 方法は、概してｎ次元対象関数データから（ｎ−１）次
元投影を決定するために使用されうる。このため、夫々
のｎ次元フーリエ変換はステップ１１において実行され
ねばならない。ステップ１２は、しばしば「重畳補間」
を使用する「リサンプリング」又は「グリッディング」
と称される。
【００１２】図４は、本発明による方法を示す他のフロ
ーチャートである。図５を参照してこの螺旋コンピュー
タ断層法のフローチャートの使用について説明する。こ
の方法によれば、Ｘ線露光中、患者は測定野を通るｚ方
向に連続的に移動され、Ｘ線管の異なる角位置に対する
一連の螺旋体積データ（投影データ）ｆ（ｚ）を生成す
る。所望の位置ｚ0 にスライス画像を形成するために、
スライスＳに沿った必要な投影値ｆ0 （ｚ）は補間によ
って、例えば同じ角位置に属する２つの隣接する投影値
から決定される。例えば、投影値ｆ0 （ｚ0 ；３／２
π）は補間によって投影値ｆ（ｚ2 ）及びｆ（ｚ4 ）か
ら決定される。
ーチャートである。図５を参照してこの螺旋コンピュー
タ断層法のフローチャートの使用について説明する。こ
の方法によれば、Ｘ線露光中、患者は測定野を通るｚ方
向に連続的に移動され、Ｘ線管の異なる角位置に対する
一連の螺旋体積データ（投影データ）ｆ（ｚ）を生成す
る。所望の位置ｚ0 にスライス画像を形成するために、
スライスＳに沿った必要な投影値ｆ0 （ｚ）は補間によ
って、例えば同じ角位置に属する２つの隣接する投影値
から決定される。例えば、投影値ｆ0 （ｚ0 ；３／２
π）は補間によって投影値ｆ（ｚ2 ）及びｆ（ｚ4 ）か
ら決定される。
【００１３】しかしながら、本発明による方法は、係る
補間を必要としない。第１のステップ１３において、投
影データｆ（ｚ）は重畳カーネルのフーリエ変換ｃ
（ｒ）と共に重畳される。そこから獲得される変更され
た投影データｆ’（ｚ）は、続いてフーリエ変換（ステ
ップ１４）を受ける。フーリエ変換Ｆ’（ｋ）は続いて
重畳カーネルＣ（ｋ）によって割り算され（ステップ１
５）、スライスＳに沿った所望の投影データは続く逆フ
ーリエ変換（ステップ１６）によって獲得される。本発
明による方法の利点は画像中にボケの影響が生じないこ
とであり、これに対して従来の方法では必要な補間によ
って係る影響がもたらされる。
補間を必要としない。第１のステップ１３において、投
影データｆ（ｚ）は重畳カーネルのフーリエ変換ｃ
（ｒ）と共に重畳される。そこから獲得される変更され
た投影データｆ’（ｚ）は、続いてフーリエ変換（ステ
ップ１４）を受ける。フーリエ変換Ｆ’（ｋ）は続いて
重畳カーネルＣ（ｋ）によって割り算され（ステップ１
５）、スライスＳに沿った所望の投影データは続く逆フ
ーリエ変換（ステップ１６）によって獲得される。本発
明による方法の利点は画像中にボケの影響が生じないこ
とであり、これに対して従来の方法では必要な補間によ
って係る影響がもたらされる。
【００１４】図６は、３次元対象関数データから２次元
円錐ビーム投影又は扇形ビーム投影を決定するためのフ
ローチャートを示す図である。まず（ステップ２０）、
３次元対象関数データに対して３次元フーリエ変換が適
用される。続いて（ステップ２１）、３次元リサンプリ
ングが行われ（図３中ステップ１２に対応）、データは
３次元周波数領域の原点を通る放射状の線上で決定され
る。各線は続いて、１次元傾斜フィルタによって掛け算
され（ステップ２２）、これらの線の１次元フーリエ変
換によってフィルタリングされた線に沿ったラドン変換
が獲得される。ステップ２３において２次元リサンプリ
ングが行われ、それによりデータは球状シェル上の２次
元極格子上にマップされる。ステップ２４において、こ
のデータは続いて、本質的には放射状の線に沿った１次
元フーリエ変換によって、続いてこれらの線に沿ったフ
ィルタリング及びデカルト格子上の逆フーリエ変換によ
って、重み付けされた円錐ビーム投影へ変換される。
所望の円錐ビーム投影は、続いてステップ２５におい
て、仮想平面から実平面（例えば、Ｘ線撮像の場合の検
出器平面）へのデータ投影によって獲得される。本発明
による方法は、上記の全てのリサンプリングステップに
おいて使用されうる。
円錐ビーム投影又は扇形ビーム投影を決定するためのフ
ローチャートを示す図である。まず（ステップ２０）、
３次元対象関数データに対して３次元フーリエ変換が適
用される。続いて（ステップ２１）、３次元リサンプリ
ングが行われ（図３中ステップ１２に対応）、データは
３次元周波数領域の原点を通る放射状の線上で決定され
る。各線は続いて、１次元傾斜フィルタによって掛け算
され（ステップ２２）、これらの線の１次元フーリエ変
換によってフィルタリングされた線に沿ったラドン変換
が獲得される。ステップ２３において２次元リサンプリ
ングが行われ、それによりデータは球状シェル上の２次
元極格子上にマップされる。ステップ２４において、こ
のデータは続いて、本質的には放射状の線に沿った１次
元フーリエ変換によって、続いてこれらの線に沿ったフ
ィルタリング及びデカルト格子上の逆フーリエ変換によ
って、重み付けされた円錐ビーム投影へ変換される。
所望の円錐ビーム投影は、続いてステップ２５におい
て、仮想平面から実平面（例えば、Ｘ線撮像の場合の検
出器平面）へのデータ投影によって獲得される。本発明
による方法は、上記の全てのリサンプリングステップに
おいて使用されうる。
【００１５】図７は、本発明による方法が実施されうる
装置を概略的に示す図である。この装置は、例えばＣＴ
装置、Ｘ線装置、又は磁気共鳴装置であり得る。例えば
Ｘ線撮像の場合は強度値である入力関数データＩ
c （ｒ）は、例えば検出器装置であるデータ捕捉ユニッ
ト３０によって捕捉される。例えば関数空間中の投写デ
ータである出力関数データＳ（ｋ）を形成するために、
このデータは所望の適用に従って上述の方法を実施する
処理ユニット３１によって処理される。画像計算ユニッ
ト３２はそこから、例えばモニタ３３上に表示されるべ
き投影画像といった画像の形成のために画像データを決
定する。
装置を概略的に示す図である。この装置は、例えばＣＴ
装置、Ｘ線装置、又は磁気共鳴装置であり得る。例えば
Ｘ線撮像の場合は強度値である入力関数データＩ
c （ｒ）は、例えば検出器装置であるデータ捕捉ユニッ
ト３０によって捕捉される。例えば関数空間中の投写デ
ータである出力関数データＳ（ｋ）を形成するために、
このデータは所望の適用に従って上述の方法を実施する
処理ユニット３１によって処理される。画像計算ユニッ
ト３２はそこから、例えばモニタ３３上に表示されるべ
き投影画像といった画像の形成のために画像データを決
定する。
【図１】空間領域中の入力関数データを示す図である。
【図２】周波数領域中の出力関数データを示す図であ
る。
る。
【図３】本発明による方法を示すフローチャートであ
る。
る。
【図４】本発明による方法を示す他のフローチャートで
ある。
ある。
【図５】螺旋コンピュータ断層撮影法のための投影デー
タを示す図である。
タを示す図である。
【図６】本発明による方法が使用されるための円錐ビー
ム投影を決定する方法を示す図である。
ム投影を決定する方法を示す図である。
【図７】本発明による方法を実施しうる装置を示す図で
ある。
ある。
１ デカルト格子 ２ 格子点 ３ 投影値 ４ 投影線 ５ 平行光線 ６ デカルト格子 ７ 格子点 ８ サンプリング点 ９ 放射状の線 ３０ データ捕捉ユニット ３１ 処理ユニット ３２ 画像計算ユニット ３３ モニタ
───────────────────────────────────────────────────── フロントページの続き (71)出願人 590000248 Ｇｒｏｅｎｅｗｏｕｄｓｅｗｅｇ １, 5621 ＢＡ Ｅｉｎｄｈｏｖｅｎ， Ｔｈ ｅ Ｎｅｔｈｅｒｌａｎｄｓ (72)発明者 ミヒャエル グラス ドイツ連邦共和国，25474 エレールベッ ク，レーペンカンプスヴェーク 18 (72)発明者 ヘンニング ブラオニーシュ ドイツ連邦共和国，30916 イーゼルンハ ーゲン，トラケーネルシュトラーセ 25
Claims (9)
- 【請求項１】 デカルト格子（１，６）の格子点（２，
７）に配置された入力関数データ（Ｉc （ｒ））から選
択された格子の格子点（８）に配置された出力関数デー
タ（Ｓ（ｋ））を決定する医用画像データ処理方法であ
って、（ａ）上記入力関数データ（Ｉc （ｒ））を重畳
カーネル（Ｃ（ｋ））のフーリエ変換（ｃ（ｒ））によ
って割り算することによって変更された入力関数データ
（Ｉc ’（ｒ））を決定する段階と、（ｂ）上記変更さ
れた入力関数データ（Ｉc ’（ｒ））のフーリエ変換
（Ｓc ’（ｋ））を決定する段階と、（ｃ）上記変更さ
れた入力関数データのフーリエ変換（Ｓc ’（ｋ））を
重畳カーネル（Ｃ（ｋ））と共に重畳することによって
上記出力関数データ（Ｓ（ｋ））を決定する段階とを含
み、 上記重畳カーネル（Ｃ（ｋ））は、上記選択された格子
の格子点（８）からの距離に従って上記デカルト格子
（１，６）の格子点（２，７）を重み付けする重み関数
である方法。 - 【請求項２】 上記入力関数データ（Ｉc （ｒ））は空
間領域中の対象関数データであり、上記出力関数データ
（Ｓ（ｋ））は周波数領域中の極格子上の投影データで
あることを特徴とする、請求項１記載の方法。 - 【請求項３】 上記出力関数データ（Ｓ（ｋ））は、特
にＸ線撮像用の、平行ビーム投影、円錐ビーム投影、又
は扇形ビーム投影の投影データであることを特徴とす
る、請求項１記載の方法。 - 【請求項４】 方法は対象の画像を回転させるために使
用されえ、このとき、回転された画像の対象関数データ
はフーリエ変換によって上記出力関数データ（Ｓ
（ｋ））から決定されることを特徴とする、請求項１記
載の方法。 - 【請求項５】 方法は画像又は画像部分の変形、変形−
消去又は拡大／縮小のために使用されえ、このとき、入
力関数データはフーリエ変換によって空間領域中の画像
データから形成されることを特徴とする請求項１記載の
方法。 - 【請求項６】 方法は螺旋コンピュータ断層撮影法にお
いて、螺旋体積データから補間データ（ｆ0 （ｚ））を
決定するために使用されうることを特徴とする、請求項
１記載の方法。 - 【請求項７】 上記重畳カーネル（Ｃ（ｋ））は、カイ
ザーベッセル関数又はガウス関数であることを特徴とす
る、請求項１記載の方法。 - 【請求項８】 上記段階（ａ），（ｂ），（ｃ）は、逆
の順序（ｃ，ｂ，ａ）で実行されることを特徴とする、
請求項１記載の方法。 - 【請求項９】 入力関数データ（Ｉc （ｒ））を捕捉す
るための手段（３０）と、出力関数データ（Ｓ（ｋ））
の処理及び／又は画像の形成のための手段（３２，３
３）と、請求項１記載の方法によって入力関数データ
（Ｉc （ｒ））から出力関数データ（Ｓ（ｋ））を決定
するためのデータ処理ユニット（３１）とを含む、特に
Ｘ線装置、コンピュータ断層撮影装置、又は磁気共鳴装
置である装置。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| DE19813926:8 | 1998-03-28 | ||
| DE19813926A DE19813926A1 (de) | 1998-03-28 | 1998-03-28 | Verfahren und Anordnung der medizinischen Bilddatenverarbeitung |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH11328358A true JPH11328358A (ja) | 1999-11-30 |
Family
ID=7862787
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP11081993A Withdrawn JPH11328358A (ja) | 1998-03-28 | 1999-03-25 | 医用画像デ―タを処理する方法及び装置 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US6151378A (ja) |
| EP (1) | EP0947958B1 (ja) |
| JP (1) | JPH11328358A (ja) |
| DE (2) | DE19813926A1 (ja) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6983182B2 (en) * | 1998-04-10 | 2006-01-03 | Wisconsin Alumni Research Foundation | Time resolved computed tomography angiography |
| JP4473358B2 (ja) * | 1999-01-21 | 2010-06-02 | 株式会社東芝 | 診断装置 |
| KR100308456B1 (ko) * | 1999-07-09 | 2001-11-02 | 오길록 | 주파수 공간상에서의 질감표현방법 및 질감기반 검색방법 |
| DE19961524A1 (de) * | 1999-12-20 | 2001-07-05 | Siemens Ag | Verfahren zum Betrieb eines CT-Gerätes |
| US6766062B1 (en) * | 2000-02-16 | 2004-07-20 | The Board Of Trustees Of The Leland Stanford Junior University - Office Of Technology | Digital ridgelet transform via digital polar coordinate transform |
| WO2003101300A2 (en) * | 2002-06-04 | 2003-12-11 | Koninklijke Philips Electronics N.V. | Rotational angiography based hybrid 3-d reconstruction of coronary arterial structure |
| DE102005046732B4 (de) * | 2005-04-18 | 2010-04-15 | Siemens Ag | Verbessertes Rekonstruktionsverfahren bei der Propellerbildgebung in der Magnetresonanztomographie |
| JP2007101249A (ja) * | 2005-09-30 | 2007-04-19 | Fujifilm Corp | 光断層画像化方法および装置 |
| CN114897033B (zh) * | 2022-07-13 | 2022-09-27 | 中国人民解放军海军工程大学 | 用于多波束窄带历程数据的三维卷积核组计算方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3637074B2 (ja) * | 1992-12-15 | 2005-04-06 | 株式会社東芝 | ヘリカルスキャン方式のコンピュータ断層撮影装置 |
| US5974110A (en) * | 1997-11-26 | 1999-10-26 | General Electric Company | Helical reconstruction algorithm |
-
1998
- 1998-03-28 DE DE19813926A patent/DE19813926A1/de not_active Withdrawn
-
1999
- 1999-03-23 DE DE59913252T patent/DE59913252D1/de not_active Expired - Fee Related
- 1999-03-23 US US09/274,390 patent/US6151378A/en not_active Expired - Fee Related
- 1999-03-23 EP EP99200884A patent/EP0947958B1/de not_active Expired - Lifetime
- 1999-03-25 JP JP11081993A patent/JPH11328358A/ja not_active Withdrawn
Also Published As
| Publication number | Publication date |
|---|---|
| EP0947958B1 (de) | 2006-03-22 |
| EP0947958A2 (de) | 1999-10-06 |
| US6151378A (en) | 2000-11-21 |
| DE59913252D1 (de) | 2006-05-11 |
| EP0947958A3 (de) | 2003-09-24 |
| DE19813926A1 (de) | 1999-09-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5737725B2 (ja) | 局所化および高度に限定された画像再構成法 | |
| US6418243B1 (en) | Apparatus and method for providing high fidelity reconstruction of an observed sample | |
| US5802133A (en) | Method and apparatus of X-ray computerized tomography | |
| NL192189C (nl) | Stelsel voor het reprojekteren van beelden onder toepassing van omzettechnieken. | |
| Kudo et al. | Virtual thin slice: 3D conditional GAN-based super-resolution for CT slice interval | |
| JPH10323347A (ja) | 対象と対象の３次元表現との間の変換を決定する方法及びその装置 | |
| JP2006015136A (ja) | 断層合成造影における直接再生方法及び装置 | |
| JP2002177254A (ja) | 造影剤に覆われた被検体のダイナミック３次元画像を再構成するための方法及びデバイス | |
| JPH11326243A (ja) | Ｃｔイメ―ジング装置を動作する方法及び走査及びデ―タ収集装置 | |
| JP3987024B2 (ja) | 横方向のフィルタリング処理を用いたトモシンセシス画像を強調する方法及びシステム | |
| JP2016152916A (ja) | Ｘ線コンピュータ断層撮像装置及び医用画像処理装置 | |
| Nilchian et al. | Optimized Kaiser–Bessel window functions for computed tomography | |
| JPH11328358A (ja) | 医用画像デ―タを処理する方法及び装置 | |
| JPH0375045A (ja) | 再生像の視野外の物体に対するｃｔデータの補償方式 | |
| US6885764B2 (en) | High Speed Z-smoothing method and apparatus for CT imaging system | |
| KR20190057122A (ko) | 화상 처리 장치, 화상 처리 방법 및 프로그램 | |
| JPH10320544A (ja) | 診断画像処理方法 | |
| JPH06261895A (ja) | Ｘ線断層撮影方法 | |
| JP2007198866A (ja) | 広義サドルコーンビームｃｔ装置および３次元再構成法 | |
| JPS6017568A (ja) | 画像処理方法および装置 | |
| US7209580B2 (en) | Fast computed tomography method | |
| EP0851393B1 (en) | Methods and apparatus for simplified filtering of scan data in an imaging system | |
| Gong et al. | Deep learning enabled ultra‐fast‐pitch acquisition in clinical x‐ray computed tomography | |
| JPH0647037A (ja) | コンピュータ断層撮影装置 | |
| Zeng et al. | Non‐iterative reconstruction with a prior for undersampled radial MRI data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060323 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060323 |
|
| A761 | Written withdrawal of application |
Free format text: JAPANESE INTERMEDIATE CODE: A761 Effective date: 20071015 |