JPH03172975A - 像形成装置および方法 - Google Patents

像形成装置および方法

Info

Publication number
JPH03172975A
JPH03172975A JP2299680A JP29968090A JPH03172975A JP H03172975 A JPH03172975 A JP H03172975A JP 2299680 A JP2299680 A JP 2299680A JP 29968090 A JP29968090 A JP 29968090A JP H03172975 A JPH03172975 A JP H03172975A
Authority
JP
Japan
Prior art keywords
filter
image
filter function
bandwidth
selection means
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.)
Pending
Application number
JP2299680A
Other languages
English (en)
Inventor
Dominic J Heuscher
ドミニク ジェイ.ホーシャー
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Philips Nuclear Medicine Inc
Original Assignee
Picker International Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Picker International Inc filed Critical Picker International Inc
Publication of JPH03172975A publication Critical patent/JPH03172975A/ja
Pending legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5205Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/005Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 〔産業上の利用分野〕 本発明hm形成装置および方法に関する。本発明は、特
にコンピュータを使用したトモグラフィー(C’l’)
走査における像の再構成に用途があるので、特にこれを
参照して説明する。しかしながら、本発明は、磁気共鳴
データ、デジタルX線データ、パネル電子管(PET 
)データ等からの像を再構成することにも利用できる。
〔従来技術および発明が解決しようとする課題〕CTス
キャナーは、像データ値のラインまたは像を発生する。
6像のデータ値は、患者のスライス(断面)を通る扇形
のサンプルの放射線に沿つ次数射線の減衰を一般を示す
。患者についての異なる角度での異なる像により同一の
扇形のサンプルを示す。像データに、データメモリに収
集され、フィルタ処理されたパックプロジェクションと
して知られているたたみこみおよびバックプロジェクシ
ョンにより最も一般的に像へと再構成される。すべての
像を共通フィルタによりフィルタ処理すると、像中の人
為的に生じたノイズ、ノイズ劣化等を低減または解消で
きる。従来はすべての像を同一のフィルタ関数でフィル
タ処理していた。
米国特許第4,555,145号は、像を再構成する前
にデータの一部のみをフィルタ処理することを開示して
いる。患者の腕を通過し、他の放射線よシも多くの骨を
通過した放射線から収集し次データは、バンド幅が異な
る。米国特許第4.33へ143号に、患者の腕を通過
した放射線データのみを所定のフィルタ関数でフィルタ
処理し、得られる像中のバンド幅に関連した人為的ノイ
ズを低減または解消することを開示している。
他のスキャナーからのデータも、収集、および表示像へ
の再構成の際にフィルタ処理されてい次が、かかるフィ
ルタに、所定像に対し、固足されている。
よって、本発明の目的は、新規でかつ改良され几像形成
装置および方法を提供することにある。
〔課題を解決するための手段および作用〕本発明の第1
の特徴によれば、1つの表示像に再構成するための複数
の像を発生するスキャナーと、フィルタ関数で6像をフ
ィルタ処理するためのフィルタ手段と、フィルタ処理さ
れた像から表示像t−再構成するための再構成手段とか
ら成る像形成装置であって、6像を処理するフィルタ関
数は一連のフィルタ関数のうちの1りであり、本装置t
は上記一連のフィルタ関数から6像をフィルタ処理する
ためのフィルタ関数を選択するフィルタ選択手段を含む
ことを特徴とする像形成装置が提供される。
本発明の第2の特徴によれは、検査領域内のサンプルの
内部を非侵入的に検査し、被検査領域の表示像に再構成
するための複数の像を含むサンプルのデータを発生し、
6像に一つのフィルタ関数による演it施し、フィルタ
処理された像を表示像に再構成する工程から成る、像形
成方法であって、6像のフィルタ処理をするフィルタ関
数は一連のフィルタ関数の一つであり、本方法は上記一
連のフィルタ関連から6像のフィルタ処理をするための
フィルタ関数を選択するための工程を含むことを特徴と
する像形成方法が提供される。
本発明の利点の一つは、再構成された像のうちの多数の
源の各々に対しフィルタ関数を最適化できることにある
以下、添附図面を参照して例によ〕本発明の像形成およ
び方法について説明する。
第1図を参照すると、CTスキャナー、出猟共鳴スキャ
ナー、PETスキャナー、デジタルX線スキャナー等の
スキャナーA[、検査領域内にあるサンプルの当該領域
すなわち、スライスまたは空間を非侵入的に検査する。
このスヤヤナーは、データ像全発生するが、6像は当該
領域の性質および6像の識別を表示する。第5世代OC
Tスキャナーでは、6像は一連のデータ値であって、こ
れらのデータfllIU放射線源と共に走査円のまわり
e回転する放射線検出器の弧に対応している。検出器は
、線源および検出器が所定の回転角距離を回転する度に
再サンプリングされ、更に像データを発生する。磁気共
鳴スキャナーでは、各像ri複数の位相二ンコーディン
グの一つでサンプリングされた周波数スペクトルの6各
に対応するデータ値を含む。データの発生は。
位相二ンコーディングの各々で反復され、複数の像を発
生するようになっている。第4世代のスΦヤナーでは、
走査円のまわシに放射線検出器が固定的に暇付けられ、
各検出器は、X線源が検出器からの走査円のまわシを回
転する際に多数回サンプリングされる。6像は、放射線
源が走査円の後方を回転する際の一つの検出器のサンプ
リングに対応する。多数の検出器の各々は、複数の像の
一つを与えるようになっている。
フィルタ手段BVi各像を選択し、これを対応するフィ
ルタ関数でフィルタ処理し、一つ以上のノイズ、喪失像
、高濃度対象、ノイズテクスチュア等を補正する。像再
構成手段Cは、一義的にフィルタ処理された像を再構成
して、再構成像を得る。
フィルタ手段Biよシ詳細に説明すると、データは、像
の表示に従って像メモリ手段10に格納される。喪失デ
ータ内挿手段12は隣接データを内挿し、喪失像のデー
タまたは他の喪失データと置換する。こうして内挿され
たデータは、像メモリ10へ戻され、喪失データにょシ
空のままになっているメモリセルを満す。この内挿手段
12は、内挿像または他の内挿データをフィルタ制御手
段14にも表示する。
読み出し手段16に、像メモリ10からのへθのサンプ
リングされた像の6像をシリアルに読み出しする。これ
とは異なシ、メモリ10は、パイプライン処理のため像
が通過するバッファ等でもよい。6像が読み出しされる
際、フィルタ制御手段14は、対応するフィルタ関数を
指定して読み出され7’C1m!をフィルタ処理する。
好ましい実施態様をより詳細に説明すると、フィルタ制
御手段にバンド幅スケーリングファクタBWFot′f
i−e’!他のフィルタ指定値を発生し、このフィルタ
指定値は、選択可能なフィルタ関数すなワチ曲線のレン
ジの各々のデジタル記述値を格納するフィルタ関数すな
わち曲線ルックアップテーブル18をアドレス指定する
。フィルタ関数は、ノツチフィルタ、ローパスフィルタ
、バイパスフィルタ、バンドパスフィルタ等ヲ含む。こ
れらフィルタは、コサイン関数に基づくものであること
が好ましい。当然ながら種々のフィルタ関数を便用する
こともできる。
フィルタテーブル18は、フィルタ曲線の限定された数
の点または値NTを記憶する。内挿手段20は、記憶さ
れたフィルタ点の間で内挿をする。すなわち、線形内挿
をし、フィルタ関数をサンプリングされた像N8のデー
タ値の数と同じ数のデータまで広げる。乗算手段22は
、ポイントごとに各フィルタ関数に対応する保全乗算す
る。読み出し手段16が次の1!’kgみ出すと、フィ
ルタ関数制御手段14は、その像に対する適当なバンド
幅スケーリングファクタを決定し、フィルタルックアッ
プテーブル18のアドレス指定をする。このように6像
に対し、フィルタ関数が別々に自動的に選択される。
一般に、6像は、実数成分および虚数成分、サインおよ
びコサインまたは相一致および相ずれ成分を含む。読み
出し手段16は、6像に対するこれらの2つのデータ成
分を別々に読み出しする。第1フイルタ制御手段14#
:t、実数成分に対する適当なフィルタを決定し、虚数
成分に対する適当なフィルタを決定する。第2フィルタ
制御手$14’は、第2バンド幅スケーリングファクタ
を選択し、このスケーリングファクタは、第2フィルタ
曲線テーブル18′のアドレス指定をし、選択したフィ
ルタ曲線の値を検索する。第2内挿手段20′は、第2
乗算手段22′がN8データ値の各々に対応フィルタ値
を乗算できるようN8値にフィルタ曲線を内挿する。
像再構成手段Cは、再構成アルゴリズム30を含む。C
Tスキャナーに対しては、従来のたたみこみおよびバッ
クプロジェクションまたはバックグロジエクションアル
ゴリズムヲ利用テきる。磁気共鳴イメージヤ−に対して
は、一般に2次元逆変換アルゴリズムを利用できる。他
の従来の再構成アルゴリズムは、スキャナーのタイプお
よび診断手続きの性質により選択できる。再構成手段は
、すべての像が処理されて完全な電子像の表示素形成す
るまで像メモリ32内に再構成データを累積する。像メ
モリから各表示像を選択的に読み出し、ビデオモニタ3
4にデイスプレィできる。種々の後処理技術を用いて像
メモリ32内の像の質を改善することもできる。増感さ
れfI:、′または生の表示像をテープまたはディスク
に記憶し、他の処理を行うこともできる。
フィルタ処理されたバックプロジェクションではフィル
タ関数によってデータ像を演算する。
このフィルタ処理演算の一部はフィルタ処理手段B内で
行うことができる。ポイント毎の乗算手段22.22’
が再構成アルゴリズムの一部を実施するように再構成フ
ィルタの全部または一部をフィルタ曲線ルックアップテ
ーブル18内に格納されたフィルタ曲線と組み合わせる
ことができる。また乗算手段36は内挿または外挿され
たフィルタ関数201Cシステムカーネル38を乗算す
ること屯できる。再構成手段Cから他の再構成アルゴリ
ズムの同様なフィルタ用部分をシフトしてもよい。
好ましい実施態様ではフィルタ制御手段14゜14′は
数徨の規準に基づいて冬休に対するフィルタ関数を選択
する。フィルタ制御手段14.14’は実質的に同一の
構造であるので以下、フィルタ制御手段14について記
述し、下記の説明は双方九同じようにあてはまると解す
べきである。
このフィルタ制御手段は、処理された冬休のノイズに従
って適当なフィルタ関数を選択するための適応型ノイズ
等化フィルタ選択手段40を含む。このノイズは先の像
のノイズ、先のスライスまたは操作の対応する像からの
ノイズ等に基づいて評価できる。またこれとは異なシ、
ノイズフィルタ選択手段が操作および冬休の実際のノイ
ズに基づきノイズフィルタを選択するよう別の処理を開
始する前に、像データの全セットを像メモリ10に収集
してもよい。
喪失像バンド幅修正フィルタ手段42は内挿法によって
形成された適当なフィルタ隣接像を選択する。即ち内挿
された像は実際に収集された像よりもバンド幅が狭く、
内挿された像の1つからの別のノイズを有することがあ
る。喪失像バンド幅修正フィルタ手段42は、内挿され
た像に対し、大小のバンド幅フィルタを検索するバンド
スケーリングファクタを形成する。保をつなげるために
選択されたバンド幅スケーリングファクタは、喪失像よ
シも幅の狭いバンド幅のフィルタ関数を検索する。
予め定めたバンド幅の修正手段44は高密度の対象ま九
は当該領域に隣接する他の材料全補償するようバンド幅
スケーリングファクタを選択する。このフィルタqCT
操作像の角度位置に従って選択でき、この角度位置は米
国特許第4.334145号におけるように像メモリ1
0から生じたデータまたはフィルタ処理されていないパ
イロット再構成像からの情報等に基づく。
ノイズテキスチユアフィルタ手段46は操作円のまわり
の角度位置と共に変化するフィルタ関数をアドレス指定
するバンド幅スケーリングファクタを選択し、よシ望ま
しいノイズテキスチユアを与える。よシ詳細には、円対
称の二次元核に基づく再構成像は混乱を生じさせうる円
対象ノイズテキスチユアを有する傾向がある。
ノイズテキスチユアフィルタ選択手段は適当な角度依存
フィルタ関数を選択し、コーナ一部でのノイズテキスチ
ユアをストレッチし、正方形または長方形デイスプレィ
によ)近く近似させる。
他の規準VC基づくフィルタ関数を選択するよう、別の
フィルタ関数選択手段48を設けてもよい。選択手段5
0は、フィルタ関数選択手段40〜48により指定され
るバンド幅スケーリングファクタの各々と所定規準とを
比較し、そのうちの一つ金アドレスとしてフィルタルッ
クアップテーブル18に送る。好ましい実施態様では、
フィルタ曲線比較手段50は、冬体に対する最大フィル
タ関数を選択する。すなわち最小バンド幅を通過する。
第2図に、代表的ローパスフィルタ関数60を示す。こ
の曲線の形状は、この関数が適用される周波数レンジt
−調節するよう62に示すよう調節される。曲線62は
、フィルタ曲線60と異なる傾きすなわち広いバンド幅
を有する。
フィルタ関数のリスケーリングに加えて、フィルタ関数
は、フィルタ曲線64jC示されるオフセット量だけシ
フトしてもよい。フィルタ曲線64は、実質的には曲線
60と同一であるがよす狭い周波数レンジをカットオフ
するようシフトされる。このフィルタ関数は、広いレン
ジ、例えばフィルタ曲線66の極小レンジまたはフィル
タ曲線68および70のほぼすべての包括的周波数レン
ジにわ次ってこれら2つの制御ファクタと共に変動する
ルックアップテーブル18は、再構成前に周波数または
7−リエ領域で使用すべきバンド制限フィルタ関数の形
状を定める。フィルタ関数制御手段14は、フィルタ関
数を適用すべき周波数レンジおよび選択されたフィルタ
関数の中心または他の周波数特性を定めるバンド幅オフ
セットBWOi定める。バンド幅スケールファクタは、
バンド無制限、すなわちアボダイズされた関数がDC値
にマツピングされるすべての値に対応するゼロ値からD
Cまたは平均化されたデータ値のみが通過されるフィル
タ値への極小ローパスフィルタリングに対応する任意の
大ファクタまで変動する。アポダイズされた関数は、数
式により次のように表わされる。
B’WF 6 = 0 (D場合 A、 (f)= T
、 (Q   (11)0≦f < N T/BWF、
の場合人a(f)=Tl(f −BWF、 )(1b) f > N T/ BWF a (D場合 A θ(f
) =T1(NT) (Ic)ここで T、(財)は、NT−t(40)を値を含むバンド幅テ
ーブル1Bから内挿手段20により内挿された値であり
、 BWFθは、表示像0の関数としてのバンド幅スケール
ファクタであシ、 サンプリングされた実関数に対してfはf=0.Δf、
 2Δf、・・・ fN(ここでΔf=fN/N8)を
とシ、 fNt−t、スペクトルが2N8個のサンプル(実〕の
離散フーリエ変換により計算されるサンプリング関数の
ナイキスト周波数である。
バンド幅ファクタBWFθは、コンピュータトモグラフ
ィーでは投射角であるθの関数として定義される。この
バンド幅すなわち解像度は、投影角θの連続関数として
変化できる。CTでは、バンド(22)は、N0個の離
散像上でサンプリングされる。第4世代OCTスキャナ
ーでは、像の像は、検出器の数に対応する。第3世代の
システムでは、像の数は、検出器がサンプリングされる
X線源の位置の数により決まる。
CTスキャナーの実m態様では、再構成手段30により
実施されるたたみこみ一バックプロジェクション貴構成
は、数式により次のように表記される。
ALG e (’) ”= ALGP(’) ・A5 
(f)        (2)ここでf;0.Δf、 
2Δf、・・・・・・f9で、ALGPは、た比みこみ
一バックプロジェクションのため予め定めたシステムの
た几みこみ核である。非線形位相補正をしないとき、こ
れらFi実関数であり、フーリエのたたみこみを行う際
に複素7−リエデータの乗*tするのに使用される。
CDO(f) = DB (f) ・ALGO(f) 
      (3a )=Dθ(f)・[ALGP(f
)・Ao(f))     (3b)ここで、D#(f
)は、角度θでの写影のフーリエ変換値であシ、ALG
θ(f) #:t、角度θで印加されるべきアルゴリズ
ムであシ、CDθ(f)は、角度θでのたたみこみ写影
の7−リエ変換値である。
DoFi複素数であるので、その結果も複素数である。
このようにALGθ(f)は、単一バンド幅のスケール
ファクタBWFθで定義され、この値は選択手段40,
42.44および46の上記像角度依存補正を同時に実
施してCT再構成像を得ることを可能にする。
バンド幅フィルタ発生と再構成核とを組み合わせるため
、ルックアップテーブル18内のアルゴリズムは2つの
クラスに分かれている。第1クラスでは(Δf*BWF
o<1)では、読み出された各新テーブル値に対し、リ
ニア進行により出力値Aθ(f)が発生される。この結
果、内部ループは極めて効率的になる。i2クラス(Δ
f*BWF O) 1の場合)では、テーブルの限度N
Tを越える点に対応する値に対する点T(NT)を外挿
することによりはとんどの値が発生される。
テーブルNT内の点の数は、サンプリングされるデータ
点N8の数よシも大幅に少ないことが好ましいので、外
挿が必要となる。上記のようにテーブルの範囲内の点は
、上記のように線形内挿法により発生される。内挿値お
よび/または外挿値を発生しなから式(2)および(3
b)の乗Xtするとき、像に依存するフィルタ処理に対
して増加した時間普は、無視できること、すなわち再構
成に時間が実質的にかからないことに注目できる。
バンド幅オフセット値は、選択されたフィルタ関数の中
心周波数、カットオフ周波数その他の周波数特性を定め
る。好ましい実施態様では、BWFoスケーリングファ
クタを有する選択されたフィルタ関数は、所定のカット
オフ周波数金有する。フィルタ選択手段は、フィルタ関
数のバンド幅をシフトさせるよう所定周波数から加減さ
れるオフセット周波数BWO=i指定することもできる
次に適応性ノイズ等化バンド幅スケールファクタ選択手
段40について説明すると、適応ノイズの等化により再
構成された対象すなわち実質的に非対称の対象内でのノ
イズの分布はよシ均一となる。多くの用途では、コント
ラストの検出度が良好に低い場合、対象全体にわたって
均一なノイズ統計値を得ることが有利である。
上記のダイナミックバンド幅のアルゴリズムでこれを達
成するには、投影角の関数として生じるノイズ変動を補
正するバンド幅スケールファクタが発生される。eTに
おけるノイズ統計値は、基本的にはX線−フォトン統計
値に基づくものであると解すべきである。
バンド幅ファクタは、次のようにパワースペクトルまた
はノイズ7アクタに関連している。
ノイズファクタBw#1/(BWB’(θ)*N3/N
T+13  (4)CTでは、ノイズパワースペクトル
は予対数減衰量に比例する。サンプリングされたデータ
からは、すべての0に対し、角度θでの最大減衰量の概
算値Am(0)が得られる。A oI’S sすべての
Am(0)の最小値を示す。
すなわち、 ノイズファクターA#、 (A−(′)鳩)(5)減衰
ノイズファクタを調節するためのバンド幅補正7アクタ
に対し、 ノイズファクター・ノイズファクターBW〜1=単位ノ
イズファクタ (6) 従って、ノイズフィルタ選択手段40は、次の方程式に
より各θに対するバンド幅スケールファクターBWFX
θ)t−決定する。
BWF(’)=(e(Am(′)−Ao)−1〕・%/
Ns   (7)方程式(7)は、BWF(θ)が常に
0以上であること、すなわち負の値を暇らないことを保
証する。
喪失gRフィルタ関数選択手段42を検討すると、喪失
像を有する領域は、投影角θに関し、アンダーサンプリ
ングされたものとみなすことができる。第4世代OCT
スキャナーでは、喪失像は、一般に不良検出器または検
出器モジュールに対応している。よって、喪失像の数お
よび位置が4見られれば、バンド幅スケール7アフタ用
の式は容易に得られる。一つの喪失像を例とすると、2
NT/N8のバンド幅スケールファクタは外挿された喪
失像に対する像サンプリング条件を満すが、隣接像は約
NT/N8のバンド幅選択7アクタを有し、(他の喪失
像をとらない)他のすべての像はゼロのファクタを有す
る。
すなわちフィルタリングがない。喪失像の数は、像のす
べての総計数と比較して通常極めて少ないので、喪失像
に関係する人為的ノイズは、解像度に測定できる程の影
響を与えることなく解消できる。
最も簡単な例における高密度対象のバンド幅フィルタ選
択手段44は、患者の体を通る像を嘔ったときのゼロ値
と、患者の腕を通る像に対する所定値との間でバンド幅
スケールファクタを変え得る。これと異なシ、他の規準
、例えば患者および患者の腕の大きさ、院内の骨組織と
首部の背柱との整合匿、他の骨組織と、撮像すべき患者
の背骨または他の領域との整合性、当該領域に隣接する
他の骨組織への代償作用等に従って、ダイナミックに選
択できる。撮像領域外の高密度対象を含む像では、干渉
する対象の高解像度および距離により像サンプリング規
準がやぶられてしまう。バンド幅スケールノアクタは、
干渉対象の所定角度および位置範囲内の適当な像が充分
にフィルタ処理されて、フィルタ処理されたデータがサ
ンプリング規準を満足するよう選択される。これは、実
施中の走査タイプに基づき予め概算できるし、その時収
集されているデータ先の隣接スライスからのデータ等か
ら計算できる。
ノイズテキスチユア手段46は、再構成された像中の非
円形対称点レスポンス関数を得るバンド幅修正ファクタ
を選択する。高次になれは、公称ノイズおよび解像度の
測定値を変えることなくノイズ統計値すなわちテキスチ
ュアを変えることができる。例えば、0°、90° 1
80°および270@で全バンド幅を維持しなから45
’。
135@ 225″′および315°でバンド@を低減
すると、正方形ま次は長方形に近似するようにコーナー
(45”、135°、225°および515°)でノイ
ズテキスチユアを引き、中間点(0090°、180°
、270°)で維持するこ七にょ夛通常円形のノイズテ
キスチユアを歪ませる。このようなノイズテキスチユア
修正を行うために適するバンド幅7アクタは、 BWF(θ)=(NT/N8) ・ (1+sin”(
H心である・・・・ (8) 再構成手段30で円形対称の再構成核を使用するとき、
式(勢のフィルタ関数は非対称点のレスポンスを与える
が、このレスポンスは360’にわたって測定した公称
レスポンスが元の像に匹敵するような一般に矩形をして
いる。このタイプのレスポンスは矩形マトリックス上で
像形成する上で特に有利である。
同じように他の像依存バンド幅スケーリングファクタも
選択できる。バンド幅スケ−リング7アクタ選択十段5
0にバンド幅ファクタ選択手段40〜48の各々により
発生されたバンド幅ファクタの内から選択する。好まし
い実施態様では選択されるバンド幅スケーリ/グアアク
タBWF(θ)は発生されるバンド幅関数のうち最大の
値である。即ち BWF(θ)==MAX(BWF、、(θ)、BWF、
、、(θ)。
BWFfi、(#)、 13WFM@(θ))−(9)
このように像表示関数としての単一バンド幅スケールフ
ァクタ0が得られる。
【図面の簡単な説明】
第1図は本発明の装置の略図、 第2図は本装置の作動を示すグラフである。 A−スキャナ、14.14’−フィルタ関数制御手段、
18.18’−フィルタ曲線ルックアップテーブル、2
0.20’−内挿/外挿手段、50−フィルタ関数選択
手段。

Claims (17)

    【特許請求の範囲】
  1. (1)1つの表示像に再構成するための複数の像を発生
    するスキャナー(A)と、フィルタ関数で各像をフィル
    タ処理するためのフィルタ手段(22)と、フィルタ処
    理された像から表示像を再構成するための再構成手段(
    C)とから成る像形成装置であつて、各像を処理するフ
    ィルタ関数は一連のフィルタ関数のうちの1つであり、
    本装置は上記一連のフィルタ関数から各像をフィルタ処
    理するためのフィルタ関数を選択するフィルタ選択手段
    (14)を含むことを特徴とする像形成装置。
  2. (2)フィルタ選択手段(14)は、各像に対するバン
    ド幅スケーリングファクタ(BWF_θ)およびバンド
    幅オフセット(BWO)の少なくともーつを発生し、本
    装置は、対応するフィルタ関数値を検索するようバンド
    幅スケーリングファクタ (BWF_θ)およびバンド幅オフセット(BWO)の
    発生された一方によりアドレス指定されるフィルタ関数
    ルックアップテーブル(18)を更に含む請求項(1)
    に記載の装置。
  3. (3)各像は第1の複数のデジタルデータ値を有し、フ
    ィルタ関数ルックアップテーブル(18)により記憶さ
    れた各フィルタ関数はより少数の複数のデジタルフィル
    タ値を有し、本装置は、各フィルタ関数が各像と同数の
    デジタル値を有するよう各フィルタ関数のフィルタ値を
    内挿するための内挿手段(20)を更に含む請求項(2
    )に記載の装置。
  4. (4)前記フィルタ手段(22)は、ポイントごとに、
    フィルタ値とデータ値とを乗算するための乗算手段(2
    2)を含む請求項(3)に記載の装置。
  5. (5)各像は、実数成分(Re)と虚数(Im)成分を
    有し、前記フィルタ選択手段(14)と、前記フィルタ
    関数ルックアップテーブル(18)と、前記内挿手段(
    20)と、前記乗算手段(22)は、フィルタ関数の値
    により前記実数成分(Re)に演算を施し、本装置は別
    のフィルタ関数の値により前記虚数成分(Im)に演算
    を施す第2フィルタ選択手段(14′)と、第2フィル
    タ関数ルックアップテーブル(18′)と、第2内挿手
    段(20′)と、第2乗算手段(22′)とを更に含む
    請求項(4)に記載の装置。
  6. (6)フィルタ選択手段(14)は、ノイズを補償する
    ための像依存バンド幅スケーリングファクタを選択する
    ためのノイズフィルタ選択手段(40)と、喪失像を補
    償するためのバンド幅スケーリングファクタを選択する
    ための喪失像フィルタ選択手段(42)と、ノイズテキ
    スチユアを調節するためのバンド幅スケーリングファク
    タを選択するためのノイズテキスチユア選択手段(46
    )のうちの少なくとも一つを含む請求項(1)〜(5)
    のいずれかに記載の装置。
  7. (7)フィルタ選択手段(14)は、各々が複数のモニ
    タ条件の異なる一つに従つてバンド幅スケーリングファ
    クタを発生する複数のバンド幅フィルタ関数選択手段(
    40、42、44、46、48)と、複数のバンド幅ス
    ケーリングファクタと所定条件とを比較し、バンド幅ス
    ケーリングファクタの一つを選択してフィルタ関数ルッ
    クアップテーブル(18)のアドレス指定をするバンド
    幅スケーリングファクタ選択手段(50)とから成る請
    求項(1)から(5)のうちのいずれかに記載の装置。
  8. (8)前記複数のバンド幅フィルタ関数選択手段(40
    、42、44、46、48)は、ノイズを補償するため
    の像依存バンド幅スケーリングファクタを選択するため
    のノイズフィルタ選択手段(40)と、喪失像を補償す
    るためのバンド幅スケーリングファクタを選択するため
    の喪失像フィルタ選択手段時と、ノイズテキスチユアを
    調節するためのバンド幅スケーリングファクタを選択す
    るためのノイズテキスチユア手段(46)と、当該像領
    域に隣接する材料を含むよう選択された像のバンド幅を
    調節する当該材料の領域の外のバンド幅スケーリングフ
    ァクタを選択するための手段(44)のうちの少なくと
    も2つを含む請求項(7)に記載の装置。
  9. (9)前記フィルタ選択手段(14)により選択される
    フィルタ関数は、選択されたフィルタ関数に関連する像
    にそれぞれ対応する請求項(1)に記載の装置。
  10. (10)フィルタ選択手段(14)は、像のノイズに従
    つてフィルタ関数を選択するためのノイズフィルタ選択
    手段(40)を含む請求項(9)に記載の装置。
  11. (11)フィルタ選択手段(14)は、喪失像と置換す
    るより内挿された像および喪失像に隣接する像を選択す
    るための喪失像フィルタ関数選択手段(42)を含む請
    求項(9)または(10)に記載の装置。
  12. (12)フィルタ関数選択手段(14)は、再構成され
    た表示像のノイズテキスチユアが変更されるよう像のた
    めのフィルタ関数を選択するためのノイズテキスチユア
    フイルタ関数選択手段(46)を含む請求項(9)、(
    10)または(11)に記載の装置。
  13. (13)検査領域内のサンプルの内部を非侵入的に検査
    し、被検査領域の表示像に再構成するための複数の像を
    含むサンプルのデータを発生し、各像に一つのフィルタ
    関数による演算を施し、フィルタ処理された像を表示像
    に再構成する工程から成る、像形成方法であつて、 各像のフィルタ処理をするフィルタ関数は一連のフィル
    タ関数の一つであり、本方法は上記一連のフィルタ関連
    から各像のフィルタ処理をするためのフィルタ関数を選
    択するための工程を含むことを特徴とする像形成方法。
  14. (14)フィルタ関数選択工程は、像ごとのバンド幅ス
    ケーリングファクタ(BWF_θ)およびバンド幅オフ
    セット(BWO)を決定し、バンド幅スケーリングファ
    クタ(BWF_θ)およびバンド幅オフセット(BWO
    )のうちの少なくとも一つによりフィルタ関数ルックア
    ップテーブル(18)のアドレス指定をし、このテーブ
    ルから対応するフィルタ関数曲線を検索することを含む
    請求項時に記載の方法。
  15. (15)フィルタ関数選択工程は、異なる規準に基づい
    て決定された、各像のための複数のバンド幅スケーリン
    グファクタを選択し、複数のバンド幅スケーリングファ
    クタと所定の規準とを比較し、バンド幅スケーリングフ
    ァクタのうちの一つを選択してルックアップテーブル(
    18)のアドレス指定をする請求項(14)に記載の方
    法。
  16. (16)複数のバンド幅スケーリングファクタは、ノイ
    ズ、喪失像、当該領域の周辺の材料およびノイズテキス
    チユアを補償する請求項四に記載の方法。
  17. (17)フィルタ関数選択工程は、フィルタ関数を選択
    し、ノイズ、喪失像、当該領域の周辺材料およびノイズ
    テキスチユアのうちの少なくとも一つを補正するフィル
    タ関数を選択することを含む請求項(13)に記載の方
    法。
JP2299680A 1989-11-17 1990-11-05 像形成装置および方法 Pending JPH03172975A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US438,687 1982-11-03
US07/438,687 US5276614A (en) 1989-11-17 1989-11-17 Dynamic bandwidth reconstruction

Publications (1)

Publication Number Publication Date
JPH03172975A true JPH03172975A (ja) 1991-07-26

Family

ID=23741619

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2299680A Pending JPH03172975A (ja) 1989-11-17 1990-11-05 像形成装置および方法

Country Status (3)

Country Link
US (1) US5276614A (ja)
EP (1) EP0429191A3 (ja)
JP (1) JPH03172975A (ja)

Families Citing this family (40)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5170439A (en) * 1991-06-11 1992-12-08 Picker International, Inc. Cone beam reconstruction using combined circle and line orbits
US5416815A (en) * 1993-07-02 1995-05-16 General Electric Company Adaptive filter for reducing streaking artifacts in x-ray tomographic images
US5424946A (en) * 1994-04-28 1995-06-13 General Electric Company Adaptive angular transmission filter for pet scanner
US5840026A (en) * 1994-09-21 1998-11-24 Medrad, Inc. Patient specific dosing contrast delivery systems and methods
US6397098B1 (en) * 1994-09-21 2002-05-28 Medrad, Inc. Data communication and control for medical imaging systems
US5739691A (en) * 1995-11-28 1998-04-14 The Regents Of The University Of California Multi-frequency digital low pass filter for magnetic resonance imaging
US6018600A (en) * 1998-04-15 2000-01-25 Arch Development Corp. Sampling and reconstruction of signals and images including MR images of multiple regions
FR2799028B1 (fr) * 1999-09-27 2002-05-03 Ge Medical Syst Sa Procede de reconstitution d'une image tridimentionnelle d'elements de fort contraste
US6438195B1 (en) * 2001-01-26 2002-08-20 Ge Medical Systems Global Technology Company, Llc Methods and apparatus for compensating for view aliasing artifacts
US6633162B2 (en) * 2001-08-22 2003-10-14 Ge Medical Systems Global Technology Company, Llc System and method for filtering frequency encoded imaging signals
US7050618B2 (en) * 2002-02-08 2006-05-23 Eastman Kodak Company Method for antiscatter stationary grid artifacts detection and attenuation in digital radiographic images
US6587537B1 (en) 2002-04-01 2003-07-01 Ge Medical Systems Global Technology Company, Llc Methods and apparatus for multi-slice image reconstruction
US7813473B2 (en) * 2002-07-23 2010-10-12 General Electric Company Method and apparatus for generating temporally interpolated projections
US6904118B2 (en) * 2002-07-23 2005-06-07 General Electric Company Method and apparatus for generating a density map using dual-energy CT
US6810102B2 (en) 2002-10-04 2004-10-26 Ge Medical Systems Global Technology Company, Llc Methods and apparatus for truncation compensation
US6856666B2 (en) * 2002-10-04 2005-02-15 Ge Medical Systems Global Technology Company, Llc Multi modality imaging methods and apparatus
US7254259B2 (en) * 2002-10-04 2007-08-07 Ge Medical Systems Global Technology Company Llc Methods and apparatus for truncation compensation
US7639774B2 (en) * 2003-12-23 2009-12-29 General Electric Company Method and apparatus for employing multiple axial-sources
US7333587B2 (en) 2004-02-27 2008-02-19 General Electric Company Method and system for imaging using multiple offset X-ray emission points
WO2005121836A1 (en) * 2004-06-09 2005-12-22 Philips Intellectual Property & Standards Gmbh Computerized tomography method with helical relative movement and conical beam
EP1812101A4 (en) 2004-11-16 2014-04-23 Medrad Inc MODELING PHARMACEUTICAL DISSIPATION
JP5232473B2 (ja) 2004-11-24 2013-07-10 メドラッド インコーポレーテッド 注入処理中に生成される圧力をモデル化するためのシステム及び装置
DE102006036327A1 (de) * 2006-08-03 2008-02-14 Siemens Ag Verfahren zum Bereitstellen von 3D-Bilddaten und System zum Aufnehmen von Röntgenbildern
US7835486B2 (en) * 2006-08-30 2010-11-16 General Electric Company Acquisition and reconstruction of projection data using a stationary CT geometry
US7706499B2 (en) * 2006-08-30 2010-04-27 General Electric Company Acquisition and reconstruction of projection data using a stationary CT geometry
US7616731B2 (en) * 2006-08-30 2009-11-10 General Electric Company Acquisition and reconstruction of projection data using a stationary CT geometry
US20080056432A1 (en) * 2006-08-30 2008-03-06 General Electric Company Reconstruction of CT projection data
EP3376504A1 (en) 2006-12-29 2018-09-19 Bayer Healthcare, LLC Patient-based parameter generation systems for medical injection procedures
EP2170165B1 (en) 2007-07-17 2018-12-05 Bayer Healthcare LLC Systems for determination of parameters for a procedure, for estimation of cardiopulmonary function and for fluid delivery
US9421330B2 (en) * 2008-11-03 2016-08-23 Bayer Healthcare Llc Mitigation of contrast-induced nephropathy
WO2011163578A2 (en) 2010-06-24 2011-12-29 Medrad, Inc. Modeling of pharmaceutical propagation and parameter generation for injection protocols
US8538114B2 (en) * 2011-06-06 2013-09-17 Kabushiki Kaisha Toshiba Method and system utilizing parameter-less filter for substantially reducing streak and or noise in computer tomography (CT) images
US9949704B2 (en) 2012-05-14 2018-04-24 Bayer Healthcare Llc Systems and methods for determination of pharmaceutical fluid injection protocols based on x-ray tube voltage
US9555379B2 (en) 2013-03-13 2017-01-31 Bayer Healthcare Llc Fluid path set with turbulent mixing chamber, backflow compensator
WO2017152036A1 (en) 2016-03-03 2017-09-08 Bayer Healthcare Llc System and method for improved fluid delivery in multi-fluid injector systems
US11598664B2 (en) 2017-08-31 2023-03-07 Bayer Healthcare Llc Injector pressure calibration system and method
WO2019046299A1 (en) 2017-08-31 2019-03-07 Bayer Healthcare Llc FLUID PATH IMPEDANCE ASSESSMENT FOR IMPROVING FLUID DISTRIBUTION PERFORMANCE
CN110809482B (zh) 2017-08-31 2023-03-07 拜耳医药保健有限公司 流体注入器系统体积补偿系统和方法
CA3067740A1 (en) 2017-08-31 2019-03-07 Bayer Healthcare Llc Method for drive member position and fluid injection system mechanical calibration
CN110869071B (zh) 2017-08-31 2023-05-02 拜耳医药保健有限公司 流体注入器系统中的动态压力控制的方法

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4288695A (en) * 1979-04-13 1981-09-08 Technicare Corporation Computerized tomographic scanner with shaped radiation filter
US4333145A (en) * 1979-11-29 1982-06-01 Technicare Corporation Method of high resolution partial area scan involving concentrated high density material outside the partial area
US4305127A (en) * 1979-11-29 1981-12-08 Technicare Corporation Projection completion method of partial area scan
NL8401739A (nl) * 1983-06-02 1985-01-02 Elscint Ltd Stelsel en werkwijze voor het opnieuw projecteren van beelden.
US4714997A (en) * 1983-06-02 1987-12-22 Elscint Incorporated Data reduction in reprojection systems
DE3342353A1 (de) * 1983-11-23 1985-05-30 Siemens AG, 1000 Berlin und 8000 München Verfahren zum betrieb eines computertomographen
US4712178A (en) * 1984-11-28 1987-12-08 Picker International Inc. Malfunctioning computed tomography detector correction method
JPS61185256A (ja) * 1985-02-13 1986-08-18 株式会社日立メデイコ X線ct画像処理装置
US4716368A (en) * 1985-08-09 1987-12-29 Picker International, Inc. Magnetic resonance reconstruction and scanning techniques using known information, constraints, and symmetry relations
NL8601004A (nl) * 1986-04-21 1987-11-16 Philips Nv Computertomografie-inrichting met gekombineerd recursief- en convolutiefilter.
US4792900A (en) * 1986-11-26 1988-12-20 Picker International, Inc. Adaptive filter for dual energy radiographic imaging
US4761819A (en) * 1987-02-27 1988-08-02 Picker International, Inc. Adaptive noise reduction filter for reconstructed images
US4887306A (en) * 1987-11-04 1989-12-12 Advanced Technology Laboratories, Inc. Adaptive temporal filter for ultrasound imaging system
US5065444A (en) * 1988-02-08 1991-11-12 Northrop Corporation Streak removal filtering method and apparatus

Also Published As

Publication number Publication date
EP0429191A2 (en) 1991-05-29
EP0429191A3 (en) 1992-06-10
US5276614A (en) 1994-01-04

Similar Documents

Publication Publication Date Title
JPH03172975A (ja) 像形成装置および方法
US6035012A (en) Artifact correction for highly attenuating objects
US9245320B2 (en) Method and system for correcting artifacts in image reconstruction
Lauritsch et al. Theoretical framework for filtered back projection in tomosynthesis
US7023951B2 (en) Method and apparatus for reduction of artifacts in computed tomography images
US7623691B2 (en) Method for helical windmill artifact reduction with noise restoration for helical multislice CT
JPH0445689Y2 (ja)
US7940884B2 (en) Method and image reconstruction device for reconstructing image data
US6108575A (en) Helical weighting algorithms for fast reconstruction
US20050123215A1 (en) Method and apparatus for the reduction of artifacts in computed tomography images
JP2000308637A (ja) Ct装置の像再構成方法
JP3049558B2 (ja) Ctスキャナによるイメージング装置の高速結合逆/順投影器
US6285732B1 (en) Methods and apparatus for adaptive interpolation reduced view CT scan
JPH05258040A (ja) 画像補間方法及び装置
US6332013B1 (en) Methods and apparatus for tilted helical reconstruction multislice CT
JP4301866B2 (ja) 物体の画像を再構成する方法及び装置
US6078639A (en) Real time continuous CT imaging
JPH10243942A (ja) 物体の基本断層画像及び後続の断層画像を形成するシステム並びに物体の画像を再構成する方法
Yu et al. A backprojection-filtration algorithm for nonstandard spiral cone-beam CT with an n-PI-window
JP4558858B2 (ja) 物体の基本断層撮影画像及び後続の断層撮影画像を形成するシステム及び方法
US5546439A (en) Systems, methods and apparatus for incrementally reconstructing overlapped images in a CT system implementing a helical scan
US7379525B2 (en) Method and system for efficient helical cone-beam reconstruction
Hsieh A practical cone beam artifact correction algorithm
US6101235A (en) Methods and apparatus for altering spatial characteristics of a digital image
US6327325B1 (en) Methods and apparatus for adaptive interpolation reduced view CT scan