JPH0434867B2 - - Google Patents

Info

Publication number
JPH0434867B2
JPH0434867B2 JP58110285A JP11028583A JPH0434867B2 JP H0434867 B2 JPH0434867 B2 JP H0434867B2 JP 58110285 A JP58110285 A JP 58110285A JP 11028583 A JP11028583 A JP 11028583A JP H0434867 B2 JPH0434867 B2 JP H0434867B2
Authority
JP
Japan
Prior art keywords
pixel value
pixel
value
frame
video signal
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 - Lifetime
Application number
JP58110285A
Other languages
English (en)
Other versions
JPS5940836A (ja
Inventor
Ei Kuruugaa Robaato
Enu Baaron Sutanree
Yuu Ryuu Pii
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.)
TOMUSON SEE ESU EFU BUROOKAASUTO Inc
Original Assignee
TOMUSON SEE ESU EFU BUROOKAASUTO 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 TOMUSON SEE ESU EFU BUROOKAASUTO Inc filed Critical TOMUSON SEE ESU EFU BUROOKAASUTO Inc
Publication of JPS5940836A publication Critical patent/JPS5940836A/ja
Publication of JPH0434867B2 publication Critical patent/JPH0434867B2/ja
Granted legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/30Cameras or camera modules comprising electronic image sensors; Control thereof for generating image signals from X-rays

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Closed-Circuit Television Systems (AREA)
  • Image Processing (AREA)

Description

【発明の詳細な説明】 (技術分野) この発明は、放射線写真像を得る方法と装置に
関し、より具体的には、ビデオ(映像)プロセツ
サと関連して蛍光透視法形の装置を用いて、不透
明化された組織の改良された像を得る方法と装置
に関する。
(背景技術) 典型的なX線の蛍光透視法形装置は、X線源と
X線を検出する検出器(image intensifier:像
増倍管)を含む。検出器の出力はテレビカメラで
観察され、その結果得られるデレビ信号はモニタ
に表示され、かつ/または記録される。患者等の
人体がX線源と検出器の間に配置されると、人体
の種々の部分の厚さや組成に応じてX線の吸収さ
れる程度が変わる。この結果得られる2次元の像
により、例えば人体内部の構造的な異常を診断す
ることができる。
この技術を用いて人体の構造を観察できるかど
うかは、当該する構造のX線吸収特性とその構造
の隣接部分の物質のX線吸収特性の関係によつて
決まる。両特性の差が大きい程、テレビ像に現わ
れる当該構造の「コントラスト」が大きくなる。
コントラストが大きい程、像における構造の明確
さが大きくなる。従つて、撮像処理において高い
コントラストを得ることが品質として要求され
る。
X線吸収作用が元来全くないかあるいはわずか
しかない場合には、放射線写真のコントラスト剤
を用いて大きな差を創り出す。例えば(胸部以外
の)血管は、血液、筋肉、脂肪および柔軟組織が
全て同じようなX線吸収作用を呈するため、実際
に蛍光透視法的な像では見ることができない。放
射線写真のコントラスト剤は、X線吸収特性が血
液、筋肉、脂肪および柔軟組織とは異質である物
質(例えば、空気、バリウム、ヨウ素)を含んで
いる。例えば、ヨウ素化液によるコントラスト物
質の大丸薬を動脈または静脈に注入した場合は、
このコントラスト物質が特定の導管体節内に存在
する間は、その導管構造は人工的に高いコストラ
ストを呈する。
従来、放射線写真像を改良するために、デイジ
タル映像処理が採用されている。像コントラスト
の増強のためのデイジタル処理技術の従来の成功
例は、放射線を通さないコントラスト剤の経時作
用のプリオリ(priori)知識を利用したものであ
る。そのようなコストラスト剤が人体内の予め選
定された動脈や器官を流れる時に、これを他と区
別して撮像することにより、器官の機能に関する
解剖学的かつ生理学的な情報を、特定の下剤を用
いることなくかつそれに付随する危険を伴なうこ
となく得ることができる。
いわゆるマスクモードの撮像は、デイジタル処
理により実行される経時的な減算撮像の直接的な
形式である。患者がX線テーブル上に乗せられ、
例えば頚動脈や心室などの当該部分が選択され
る。細い針が腕の静脈に刺され、30〜50mlのヨウ
素化コントラスト剤が注射される。注射の前に、
1つのデイジタル像がいくつか(典型的には1な
いし4)の映像フレーム上に形成され、かつデイ
ジタルメモリに記憶される。次いでコントラスト
剤が急速(例えば3ないし5秒間)に注入され、
このコントラスト剤は右心室、肺、左心室と流
れ、左心室から人体の動脈系全体へ圧送される。
コントラスト物質が当該部分を通過する時、一連
の付加的な像が第2のデイジタルメモリに累積さ
れる。これらの不透明化像のそれぞれは、注入前
の像から連続的に減算される。このようにして形
成された減算像は、比較的静止している動脈構造
に対して1秒間当たり約1つの像の割合で合成さ
れる。心臓のような急激に動く構造に対しては、
その割合は1秒間当たり15ないし16の像に増大す
る。
このようにして形成された経時的な減算像は、
理想的には不透明化された心臓導管の組織のみを
表示する。患者が全く動かない場合には、例えば
骨などの不透明化されない組織に基づく像コント
ラストは除去される。無関係な像の情報を取り除
くことにより、不透明化された組織のコントラス
トを増強することができる。増強された減算像
は、映像フオーマツトに再び変換され、例えばビ
デオテープやビデオデイスク等へ記憶される。全
ての処理と外部記憶は、テーブル上の患者と実時
間で行なうことができる。
本出願人に譲渡された同時係属の米国特許出願
第342376号明細書には、循環形のフイルタリング
を用い、撮像している人体の部分を通過するコン
トラスト物質の造影剤の移動について予測される
一時的な周波数に実質的に対応する一時的な周波
数を持つフイルタ機能により、放射線写真で生成
した映像信号をフイルタリングする技術が開示さ
れている。このようにして処理された映像は、静
止組織および通過帯域ピークより実質的に高い一
時的周波数で動く組織からの影響を減衰させなが
ら、コントラスト物質が通過する部分の像を増強
させる傾向を示す。
上述したように、減算技術と循環形フイルタリ
ング技術の双方によつて、放射線写真で生成され
た像の改良が助長されるが、この発明の目的は、
さらに改良を達成する映像増強技術を提供するこ
とにある。
(発明の開示) 上述した減算技術は、ほとんどのデータ(すな
わち減算を行なう映像フレーム以外の映像フレー
ム)を棄てる遡及的なフイルタリングと考えられ
る。これに対して、上上述した循環形フイルタリ
ング技術は、予測的なフイルタ特性を選択し、フ
イルタの選択はコントラスト希釈曲線の形状に関
するプリオリ仮定に基づいている。こらの観点か
ら、従来技術は妥協を示している。
この発明は、収集したデータを実質的に全て使
用し、コントラスト物質が通過する際の人体の当
該部分の実際の撮像特性の予測的な知識を用いて
いる。
この発明の実施例によれば、人体の像を表わす
処理された映像信号の発生方法およびその装置が
提供される。放射線源が人体に向けられ、人体か
らの放射線を検出する検出手段と、検出された放
射線を電子的な映像信号に変換する変換手段が設
けられる。映像信号の一連のフレームをデイジタ
ルなピクセル値(pixel value)として記憶する
記憶手段が設けられる。このデイジタルなピクセ
ル値は一般的に、フレームの要素的なピクセル位
置に対応する人体部分の放射線吸収特性を表わす
が、必要であれば、この発明の原理に従つて、透
過特性などの他の特性を用いて処理してもよいこ
とは理解されよう。一連のフレームiのそれぞれ
に対して平均ピクセル値Siを決定する平均ピクセ
ル値決定手段が設けられる。また、一連のフレー
ムiのそれぞれに対して、Siの関数として増強値
(またはフイルタ係数)kiを決定する増強値決定
手段が設けられる。好ましくは、増強値kiはSiと
Sの差として決定され、ここには全レームに亘
るSiの平均値である。増強値kiは正または負のど
ちらの値をもとり得る。さらに、記憶された各フ
レームiのピクセルのピクセル値を増強値Kiと
結合して、増強ピクセル値を得るピクセル値結合
手段が設けられる。さらに、一連のフレームの対
応するピクセルの増強ピクセル値を結合して、各
ピクセル位置における処理されたピクセル値を得
る増強ピクセル値結合手段が設けられる。処理さ
れたピクセル値は表示のためにアナログ値に変換
してもよい。
本実施例では、各フレームiに対して増強値ki
を決定するために使用される吸収密度は、窓領域
で規定されたフレームの予め選定された部分の範
囲内で得られる。もしも窓領域が、フレームの処
理操作中にコントラスト物質の大丸薬が通過する
人体の当該武分のほとんどまたは全てを含むよう
に選択された場合は、kiを窓領域内のピクセル値
によつて決定するように制限することにより、ki
の決定に際してフレーム領域を考慮する場合より
も、有効なフイルタ特性が大丸薬の一時的応答に
対してより良好な「調和」を呈する。
本出願人は、この発明の技術によれば信号対雑
音比が改良され、そのため従来技術に比較して、
所定の品質の像を得るために必要なコントラスト
物質の投与量が著しく減少することを確認した。
この発明のその他の特徴と効果は、図面を参照
した以下の詳細な説明から明らかになろう。
(発明を実施するための最良の形態) 第1図のブロツク図は、人体20の内部構造の
表示された像を得るための装置100を示す。装
置100は従来と同様に、放射線源110、典型
的にはX線源と、検出器120と映像発生器13
0を含む。検出器と映像発生器の組合せは、例え
ばX線像増倍管とテレビカメラでよい。映像発生
器130の出力端子はビデオレコーダ150と
A/D(アナログ/デイジタル)変換器151に
接続され、A/D変換器151はテレビ信号をデ
イジタル信号に変換する。ビデオレコーダ150
の出力端子もA/D変換器151に接続される。
デイジタル化されたテレビ信号を得る装置は周知
かつ市販されており、例えばトムソン−セーエス
エフブロードカスト社製のモデルAD−964310が
ある。このテレビ信号のデイジタル化装置は、映
像フレームの各ピクセル(pixel)について、256
段階の輝度レベル(本実施例では1つの単色信号
を考慮している)のうちの1段階を表わす8ビツ
トのデイジタル信号を発生する。ビデオレコーダ
150は、例えばビデオテープレコーダやビデオ
デイスクレコーダなどの適当な記録装置でよい。
映像発生器130とビデオレコーダ150とA/
D変換器151は従来と同様に同期信号を受信
し、A/D変換器151はピクセル速度でクロツ
ク信号を受信する。
A/D変換器151の出力端子は処理回路(ま
たはプロセツサ)200に接続され、この処理回
路200は汎用デイジタル計算機、汎用マイクロ
プロセツサまたは特殊デイジタル電子装置で構成
される。A/D変換器151の出力端子は四則累
積器160は処理回路200に接続されてその制
御を受ける。映像発生器130とビデオレコーダ
150も処理回路200により制御される。処理
回路200は例えばモトローラ社製のモデル
68000型のマイクロプロセツサ、またはいわゆる
その「ビツトスライス(bit slice)」実施で構成
されるものでよい。
A/D変換器151の出力端子はまた、処理回
路200により制御される乗算回路170に接続
される。乗算回路170の出力端子は、フレーム
メモリ305と加算回路(加算器)306を含む
もう1つの累積器300に接続される。累積器3
00の出力端子はD/A(デイジタル/アナログ)
変換器195に接続され、D/A変換器195の
出力端子は表示器198およびビデオレコーダに
接続され、ビデオレコーダはビデオレコーダ15
0でよい。さらに、垂直同期信号と水平同期信号
およびクロツク信号が、従来と同様に同期および
クロツク発生回路(図示しない)から得られる。
システムの動作と制御を説明する前に、いくつ
かの理論的考察を述べる。
一時的フイルタリング法は全て、周波数領域に
おいてこれに関連する一時的周波数応答(振幅お
よび位相)によるか、あるいは時間領域において
そのインパルス応答により表わすことができる。
フイルタの入力信号をs(t)で表わし、そのフ
イルタのインパルス応答をf(t)で表わすと、
フイルタの出力s′(t)は、 s′(t)=s(t)*f(t) ……(1) =∫ -∞s(t′)f(t−t′)dt′ ……(2) で与えられる。但し*はたたみこみを示す。信号
が間隔T秒のs(t)の1組のN個の離散的なサ
ンプルである場合、フイルタの出力は、 s′〔i〕=j=-∞ s〔j〕f〔i−j〕 ……(3) で与えられる離散的な和であり、ここにs〔i〕
はs(t)のi番目のサンプルであり、f〔i〕は
f(t)のi番目のサンプルである。原因となる
入力信号と原因となるフイルタに対して、(3)式は s′〔i〕=ij=0 s〔j〕f〔i−j〕 ……(3′) となる。n番目の出力は、 s′〔N〕=Nj=0 s〔j〕f〔N−j〕 ……(4) または、 s′〔N〕=Nj=0 kjs〔j〕 ……(4′) となり、但し、 kj=f〔N−j〕 ……(5) である。(4′)式は入力サンプルを全て重み付け
してN番目のサンプルまで加算したものである。
1組の離散的なデータサンプルに対して施され
るフイルタリングの方法は全て、データサンプル
の線形結合と考えることができる。この点を留意
すると、本明細書の背景技術において述べたマス
クモードの減算撮像は、kjが kj=1,J=l −1,J=m 0,他の全てのj によつて選択された遡及的なフイルタリングの特
別な場合である。l番目のサンプルは不透明化が
最大の時刻の近くに選ばれ、m番目のサンプルは
不透明化が低い時刻に選ばれる。他のN−2個の
サンプルは棄てられる。
前述した同時係属の米国特許出願第342376号明
細書における循環形フイルタリングに対しては、
ほとんどのkjはゼロではない。従つて各データサ
ンプルは使用されるが、しかし、フイルタインパ
ルス応答に関連するkjは予測的に(例えばコント
ラスト物質の一時的応答に調和するように)選択
され、kjは1組のサンプルs〔j〕のサンプル値
が実際に何であるかということに関する知識には
影響を受けない。
従つて減算技術は、ほとんどのデータを棄てる
遡及的フイルタリング過程と考えることができ
る。ほとんどのデータを使用する循環形フイルタ
リングは、コントラスト希釈曲線の形状に関する
プリオリ仮定に基づいて選択を行なう、固定され
た予測的フイルタを採用する。各方法は妥協を示
す。
この発明のフイルタリング法は、収集されたデ
ータを実質的に全て使用し、かつs〔j〕の既に
収集したデータに基づく知識を組み込むもので、
両方の撮像方法の長所を併せ持つ。1組のサンプ
ルs〔j〕が与えられると、選択された規準に従
つてNj=0 kis〔j〕を最適にするような方法でkjを選
択することができる。N番目の出力の信号対雑音
比が最大になると共に、同時に静止した背景組織
が消去されるようにkjが選択される場合は、フイ
ルタリング手順に対する1つのピクセルに関連す
る信号Sと雑音nは、 S≡s′〔N〕=Nj=0 kis〔j〕 ……(6) n≡〔Nj=0 kl2σ21/2 ……(7) で与えられる。ここで、σは各サンプルに対する
ピクセル当たりの雑音の標準偏差である。(7)式
は、サンプル雑音σはサンプルからサンプルへは
相関せず、かつσはjとは独立していると仮定し
ている。N番目の出力における静止した背景組織
が削除されることを保証するため、Nj=0 kj=0 という条件を課す。この条件により、(6)式は、 S=Nj=0 kjs〔j〕=Nj=1 kj(s〔j〕−c) ……(8) と書き直せる。但しcは任意の定数である。S/
nは、 ∂(S/n)/∂ki i=0,1…N ……(9) と設定することにより最大となる。(7)式と(8)式を
(9)式に代入してkiについて解くと、 が得られる。(10)式は、 ki=a({ki})(s〔i〕−c) i=0,1…N
……(11) と書き直すことができ、ここでa({ki})はkjに
よつて決まる比例定数である。
a({ki})の大きさがS/nにどの程度影響す
るかを調べるため、(12)式を(7)式と(8)式に代入して
S/nを求めると、 となり、これはa({kj})には無関係である。簡
単のため、a({kj})=1と選択する。従つて、
S/nを最大にするkiの1つの選択は、 ki=s〔i〕−c ……(15) となる。Ni=0 ki=0であることを思い起こすと、C
は、Ni=0 ki=Ni=0 (s〔i〕−c)=0 〓c=1/N+1Ni=0 s〔i〕≡ ……(16) 但し、はs〔i〕の平均値である。従つて最
終的に、S/nを最大にするkiの選択は、 ki=s〔i〕− i=0,1,2…N
……(17) となる。(5)式から、フイルタのインパルス応答
は、 f〔i〕=kN-i=S〔N−i〕−……(17) (17)式を用いて、(7)式と(8)式からS/nを、 あるいは、 S/n=(N+1)1/2(s〔i〕−)rns/σ……
(18) で表わすことができる。但し、 (s〔i〕−)rns≡{1/N+1Ni=0 (s〔i〕−)21/2 ……(19) であり、s〔i〕−)rnsはフイルタリングに使用
されるサンプルに対してのみ表現される。N+1
は処理のために保持された入力サンプルの数であ
る。
第2図は基本的なフローチヤートで、同図中に
示される第3図ないし第5図のフローチヤートを
参照して、本実施例による第1図の装置100の
回路の動作を制御するために処理回路200がど
のようにプログラムされるかを示す。ブロツク2
01は、各フレームiに対する平均ピクセル値Si
の決定を示す。この操作の“パス”と指定され
るルーチンは第3図に示される。一連のフレーム
のそれぞれについて求められた平均ピクセル値Si
は、ブロツク202に示すように、平均をとられ
て全体平均ピクセル値を求める。この操作のル
ーチンは第4図に示される。次いでブロツク20
3に示すように、前述した(17)式に従つて、各フ
レームに対する増強値kiが計算される。次にブロ
ツク204に示すように、パスルーチンの実行
により、一連のフレームiのそれぞれについて各
ピクセルのピクセル値にkiが乗算される。この操
作のルーチンは第5図に示される。一連のフレー
ムのそれぞれからの対応する増強ピクセル値は、
前述の(4′)式に従つて累積器300(第1図)
において累積される。各増強ピクセル値が乗算回
路170から出力されると、その増強ピクセル値
は加算回路305により、フレームメモリ306
から送り返された特定のピクセルに対する累積値
に加算される。この新しい累積値にフレームメモ
リ306の特定のピクセル位置に再び記憶され
る。全てのフレームの処理が終ると、フレームメ
モリ306は所望の累積ピクセル値を各ピクセル
位置に保持する。
第3図は、一連のフレームiのそれぞれに対す
る平均ピクセル値Siを求めるためのパスルーチ
ンを示す。ブロツク211は操作者による窓領域
の設定を示す。前述したように、窓は、本実施例
の全体平均ピクセル値を得る領域を規定する。
操作者は、例えば公知のカーソルを窓領域の向か
い合う角部を対角線状に移動して、その角部の座
標を入力することにより、窓領域を設定すること
ができる。あるいは、周知のごとくキーボードに
より、窓領域の角部の座標を入力することもでき
る。また操作者は、処理しようとする一連のフレ
ームのフレーム数を入力することができる(ブロ
ツク212)。四則累積器160(第1図)がク
リヤされ、フレームカウントiが1に初期設定さ
れる(ブロツク213)。ブロツク214は、映
像発生器130からA/D変換器151を介して
次の映像フレームを入力することを示す。次いで
ループ220の実行が開始されるが、このループ
220は、フレームiについて平均ピクセル値Si
を求める手順の実行を示す。具体的には、個々の
ピクセル値がA/D変換器151から読み込まれ
(ブロツク215)、状況に応じて、ピクセル要素
またはラインカウントがインクリメントされる
(ブロツク216)。現在のピクセルが予め設定さ
れた窓領域の範囲内にあるものか否かが判定され
る(ダイヤモンド217)。yesであればブロツク
218に入り、累積器のカウントを現在のピクセ
ル値によりインクリメントする。ピクセルが窓領
域の範囲外であれば、ブロツク218は飛ばされ
る。次いで、窓領域の範囲外であれば、ブロツク
218は飛ばされる。次いで、窓領域の最後のピ
クセルが処理され終つたか否かが判定される(ダ
イヤモンド221)。最後のピクセルでなければ、
フレームiが完了するまでループ220が繰り返
される。(本実施例では、窓領域内のピクセルの
みがSiの決定に用いられるので、窓領域の最右下
端のピクセルに到達した後は、ピクセルの判定を
継続する必要はない。)次いでSiの値がフレーム
番号iと共にメモリに記憶される(ブロツク22
2)。次に累積器がクリヤされ(ブロツク22
3)、一連のフレームの最後のフレームの処理が
終つたか否かが判定される(ダイヤモンド22
4)。最後のフレームでなければ、フレームカウ
ントをインクリメントし(ブロツク235)かつ
ブロツク214に再び入ることにより、ループ2
30が継続して実行される。ループ230は、一
連のフレームについて全てのSiが決定されかつ記
憶され終るまで、継続して実行される。
第4図のフローチヤートは、第2図のブロツク
202の詳細すなわち個々の平均ピクセル値Siの
平均がとられて全体平均ピクセル値を求めるル
ーチンを示す。フレームカウントiが初期設定さ
れ(ブロツク411)、Siの値を累積するための
変数SBARが0に初期設定される(ブロツク41
2)。記憶されたSiの値とiが呼び出される(ブ
ロツク413)。次いでブロツク414に示すよ
うに、変数SBARが呼び出されたSiだけ増加され
る。一連のフレームの最後のフレームに到達した
か否かが判定される(ダイヤモンド415)。最
後のフレームでなければ、フレームカウントiを
インクリメントし(ブロツク416)、全てのSi
の値を変数SBARに加え終るまで、ループ410
を継続する。ダイヤモンド415で最後のフレー
ムに到達したと判断したら、ブロツク421に入
り、和SBARを一連のフレームの数で割つて全体
平均ピクセル値を求め、次いでこの値を記憶
する(ブロツク422)。
第5図のフローチヤートは第2図のブロツク2
04のパスルーチンを示し、このパスルーチ
ンでは、一連のフレームiのそれぞれについて、
各ピクセルのピクセル値に増強値kiが乗算され、
増強されたピクセル値が累積器300(第1図)
によつて一連のフレームの中の他のフレームの対
応するピクセルのピクセル値に加算される。フレ
ーム番号iが初期設定され(ブロツク511)、
増強値kiが乗算回路170に出力される(ブロツ
ク512)。i番目のフレームの各ピクセルが、
A/D変換器151から1時に1個ずつ(ブロツ
ク514)乗算回路170に入力される(ブロツ
ク513)。そのフレームの最後のピクセルが処
理される(ダイヤモンド515)まで、ループ5
20が続く。次いで、一連のフレームの最後のフ
レームの処理が終つたか否かが判定される(ダイ
ヤモンド525)。最後のフレームでなければ、
フレーム番号iがインクリメントされ(ブロツク
526)、一連のフレームが全て処理され終るま
で、ループ530が繰り返される。
以上、この発明を特定の実施例を参照して説明
したが、この発明の精神と特許請求の範囲の範囲
内の変形例は、当業者には明らかであろう。例え
ば、必要であれば、開示した技術の手順のうちの
特定のいくつかは、手動制御により、あるいは特
殊制御回路を用いることにより、実行することが
できる。
【図面の簡単な説明】
第1図はこの発明の方法および装置を説明する
ブロツク図、第2図は第1図の処理回路を制御す
るための基本的なフローチヤート、第3図は第2
図のパスルーチンの詳細を示すフローチヤー
ト、第4図は第2図のを決定するルーチンの詳
細を示すフローチヤート、第5図は第2図のパス
ルーチンの詳細を示すフローチヤートである。 20…人体、110…放射線源、120…検出
器、130…映像発生器、150…ビデオレコー
ダ、151…A/D変換器、160…四則累積
器、170…乗算回路、195…A/D変換器、
200…処理回路、300…累積器、305…加
算回路、306…イメージメモリ、i…フレー
ム、Si…平均ピクセル値、…全体平均ピクセル
値、ki…増強値。

Claims (1)

  1. 【特許請求の範囲】 1 人体に向けられる放射線源と、該人体からの
    放射線を検出する検出手段と、検出された放射線
    を電子的な映像信号に変換する変換手段とを含む
    X線システムと共に用いられる、前記人体のX線
    像を表わす処理された映像信号の発生装置におい
    て、 前記映像信号による一連のフレームiを各ピク
    セル位置におけるピクセル値として記憶する記憶
    手段と; 前記一連のフレームiのそれぞれについて該フ
    レームの平均ピクセル値Siを決定する平均ピクセ
    ル値決定手段と; 前記一連のフレームiのそれぞれについて前記
    平均ピクセル値Siの関数として該フレームの増強
    値kiを決定する増強値決定手段と; 記憶されたフレームiのそれぞれの各ピクセル
    の前記ピクセル値を該フレームの前記増強値kiと
    結合して、増強ピクセル値を得るピクセル値結合
    手段と; 前記一連のフレームの対応するピクセルの前記
    増強ピクセル値を結合して、各ピクセル位置にお
    ける処理されたピクセル値を得る増強ピクセル値
    結合手段とから構成されるX線像の処理された映
    像信号の発生装置。 2 各増強値kiが差(Si−)(但し、は平均
    ピクセル値(Si)の全体平均値)から得られる特
    許請求の範囲第1項記載のX線像の処理された映
    像信号の発生装置。 3 平均ピクセル値決定手段が、フレームiのピ
    クセル値の和を求めるものである特許請求の範囲
    第1項または第2項記載のX線像の処置された映
    像信号の発生装置。 4 平均ピクセル値決定手段が、フレームiの全
    体より小さな、予め選定された部分の範囲内で該
    フレームiのピクセル値の和を求めるものである
    特許請求の範囲第1項または第2項記載のX線像
    の処理された映像信号の発生装置。 5 各ピクセル値が、該ピクセルの位置における
    人体の要素部分のX線の吸収を表わす特許請求の
    範囲第1項または第2項または第4項記載のX線
    像の処理された映像信号の発生装置。 6 増強ピクセル値結合手段が、対応するピクセ
    ルのピクセル値を加えるものである特許請求の範
    囲第1項または第2項または第4項記載のX線像
    の処理された映像信号の発生装置。 7 ピクセル値結合手段が、フレームiの各ピク
    セル値に増強値kiを乗算するものである特許請求
    の範囲第1項または第2項または第4項記載のX
    線像の処理された映像信号の発生装置。 8 処理されたピクセル値をデイジタル値からア
    ナログ値に変換する変換手段を含む特許請求の範
    囲第1項または第2項または第4項記載のX線像
    の処理された映像信号の発生装置。
JP58110285A 1982-06-21 1983-06-21 X線像の処理された映像信号の発生方法およびその装置 Granted JPS5940836A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US390512 1982-06-21
US06/390,512 US4456926A (en) 1982-06-21 1982-06-21 Enhancement of fluroscopically generated images

Publications (2)

Publication Number Publication Date
JPS5940836A JPS5940836A (ja) 1984-03-06
JPH0434867B2 true JPH0434867B2 (ja) 1992-06-09

Family

ID=23542763

Family Applications (1)

Application Number Title Priority Date Filing Date
JP58110285A Granted JPS5940836A (ja) 1982-06-21 1983-06-21 X線像の処理された映像信号の発生方法およびその装置

Country Status (5)

Country Link
US (1) US4456926A (ja)
EP (1) EP0101336B1 (ja)
JP (1) JPS5940836A (ja)
CA (1) CA1205213A (ja)
DE (1) DE3380324D1 (ja)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4792854A (en) * 1982-09-14 1988-12-20 New York Institute Of Technology Apparatus for temporally processing a video signal
US4542459A (en) * 1982-11-26 1985-09-17 General Electric Company Matched filter for x-ray hybrid subtraction
JPS59103648A (ja) * 1982-12-06 1984-06-15 株式会社東芝 X線撮像装置
IL70214A (en) * 1983-11-13 1987-10-20 Elscint Ltd Image contrast enhancement arrangement
NL8801982A (nl) * 1988-08-09 1990-03-01 Philips Nv Werkwijze en inrichting voor het bepalen van grootheden die kenmerkend zijn voor de stroming van een vloeistof in een vatenstelsel.
DE19625727C2 (de) * 1996-06-27 2000-11-09 Bernd Porr Einrichtung zur Bestimmung der Disparität in einem Stereobildpaar
US6226350B1 (en) * 1998-12-31 2001-05-01 General Electric Company Methods and apparatus for cardiac scoring with a multi-beam scanner
EP1090306A1 (en) * 1999-04-20 2001-04-11 Koninklijke Philips Electronics N.V. Image processing method
US6256368B1 (en) * 1999-10-15 2001-07-03 General Electric Company Methods and apparatus for scout-based cardiac calcification scoring
DE10238322A1 (de) * 2002-08-21 2004-03-11 Siemens Ag Retrospektive bzw. fenstergesteuerte Filterung von Bildern zur Adaption von Schärfe und Rauschen in der Computer-Tomographie
DE10345073A1 (de) * 2003-09-26 2005-05-04 Siemens Ag Betriebsverfahren für ein tomographiefähiges bildgebendes Untersuchungsgeräts und Röntgen-Computertomographiegerät
FR2884340B1 (fr) * 2005-04-11 2012-03-02 Gen Electric Procede et dispositif de traitement d'images en angiographie soustraite
EP4139881A4 (en) * 2020-05-18 2023-03-01 Shanghai United Imaging Healthcare Co., Ltd. SYSTEMS AND METHODS FOR IMAGE OPTIMIZATION

Also Published As

Publication number Publication date
EP0101336A2 (fr) 1984-02-22
JPS5940836A (ja) 1984-03-06
CA1205213A (en) 1986-05-27
US4456926A (en) 1984-06-26
EP0101336B1 (fr) 1989-08-02
EP0101336A3 (en) 1985-11-06
DE3380324D1 (en) 1989-09-07

Similar Documents

Publication Publication Date Title
JP4683721B2 (ja) 器官の厚みを補償する方法および装置
US4436095A (en) Method and apparatus for imaging a body
US4662379A (en) Coronary artery imaging system using gated tomosynthesis
EP2320801B1 (en) Mask construction for cardiac subtraction
US20050207538A1 (en) Method of determining an image from an image sequence
JPH0332950B2 (ja)
US4456926A (en) Enhancement of fluroscopically generated images
US20050277830A1 (en) X-ray CT apparatus and myocardial perfusion image generating system
US4677478A (en) Apparatus and method for imaging a body
KR20060028682A (ko) 심혈관 음향을 분석하기 위한 방법과 시스템
JP2006289094A (ja) 作成した画像のノイズを処理するための方法及びシステム
JP2670635B2 (ja) エネルギーサブトラクション画像の表示方法および装置
JPH0214051B2 (ja)
JP2009054013A (ja) 画像処理装置
JP3400063B2 (ja) コンピュータ断層撮影装置
EP0097355A1 (en) Method and apparatus for performing digital intravenous subtraction angiography
Kiuru et al. Dynamic imaging of pulmonary ventilation: Description of a novel digital fluoroscopic system
JP2544638Y2 (ja) 画像処理装置
JPH10229519A (ja) X線画像デジタル処理装置
RU2633949C1 (ru) Способ формирования составного изображения из серии ангиографических цифровых субтракционных кадров
Garland et al. Synthetic mask energy subtraction angiography: phantom study
JPS62104388A (ja) X線画像処理装置
RU2644928C1 (ru) Способ определения кадров, соответствующих границам фаз кровообращения, при проведении ангиографического исследования (варианты)
Reinke et al. Pulmonary videodensitometry in the diagnosis of carcinoma of the lung
JPH06125499A (ja) ディジタル・サブトラクション・アンギオグラフィ装置