JPH0227873B2 - - Google Patents

Info

Publication number
JPH0227873B2
JPH0227873B2 JP58025230A JP2523083A JPH0227873B2 JP H0227873 B2 JPH0227873 B2 JP H0227873B2 JP 58025230 A JP58025230 A JP 58025230A JP 2523083 A JP2523083 A JP 2523083A JP H0227873 B2 JPH0227873 B2 JP H0227873B2
Authority
JP
Japan
Prior art keywords
image
position correction
contrast
sum
opposite polarity
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
JP58025230A
Other languages
English (en)
Other versions
JPS59151941A (ja
Inventor
Kyoshi Okazaki
Katsuji Takasu
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.)
Toshiba Corp
Original Assignee
Tokyo Shibaura Electric Co Ltd
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 Tokyo Shibaura Electric Co Ltd filed Critical Tokyo Shibaura Electric Co Ltd
Priority to JP58025230A priority Critical patent/JPS59151941A/ja
Publication of JPS59151941A publication Critical patent/JPS59151941A/ja
Publication of JPH0227873B2 publication Critical patent/JPH0227873B2/ja
Granted legal-status Critical Current

Links

Landscapes

  • Analysing Materials By The Use Of Radiation (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Description

【発明の詳細な説明】 〔発明の技術分野〕 本発明は、被検体を透過した放射線を検出・処
理することにより、医学的診断に有効な画像を表
示することのできる放射線診断装置に関するもの
である。
〔発明の技術的背景とその問題点〕
従来より、放射線診断装置において、被検体の
撮影目的部位に造影剤が注入される前に撮影した
少なくとも1枚以上積分した画像と、前記撮影目
的部位における関心部位(以下診断目的部位とも
いう)に造影剤が注入された後に撮影した少なく
とも1枚以上積分した画像とを引算器に入力し、
この引算器において前記各画像における各画素毎
の濃度値を引算することにより、造影剤が存在す
る部位のみの画像を抽出することが行なわれてい
る。ところが、この処理に供される2枚の画像を
収集する間に被検体が複雑に動くことが多く、こ
のような場合には造影血管又は造影臓器の影像の
他に背景の位置ずれによる偽像が重なり、診断を
不可能にするという欠点を有していた。
上記欠点を克服するために、従来においてマニ
ユアル位置補正装置を備えた放射線診断装置と最
小2乗誤差法による位置補正装置を有する放射線
診断装置とが考え出された。しかし、前者の場合
には、人間が常に介在するので位置補正係数を求
めるために処理時間が長くなるし、技巧を有する
という問題点があり、一方後者の場合には、造影
前後の画像間では求めた位置補正係数が正しい値
を示さず、位置ずれがある画像を最適なものとし
てしまうという問題点を有していた。
尚、第1図は位置ずれが起きた場合と位置ずれ
が起きない場合とにおいて、各画像の濃度値を位
置の関数として模式的に示したものである。ここ
において、第1図aは造影前の画像であり、第1
図bは造影中の画像で第1図aの位置とは1画素
分ずれている。図中の斜線は診断目的部位内に現
われる造影信号分を示す。第1図cは第1図aの
画像から第1図bの画像間の各画素毎に差分を求
めたものであり、散点部は偽像を示す。第1図d
は位置ずれが起きていない時の造影中の画像であ
る。そして第1図eは第1図aの画像から第1図
dの画像間の各画素毎に差分を求めたものを示
す。
上記のような場合において、後者の最小2乗誤
差法によると、第1図cとeから明らかに第1図
cが最適であると誤認する。
〔発明の目的〕
本発明は、前記事情に鑑みてなされたものであ
り、被検体の動きによりサブトラクト像に出現す
る位置ずれ偽像を、最小負成分差法に基づくアル
ゴリズムを備えた自動位置補正装置により消去
し、造影血管又は造影臓器影像のみをサブトラク
ト画像として表示できる放射線診断装置を提供す
ることを目的とするものである。
〔発明の概要〕
上記目的を達成するための本発明の概要は、造
影剤が被検体の診断目的部位に混入する以前のマ
スク像よりなる第1の放射線透過画像群と、造影
剤が前記診断目的部位に到達した後の造影像より
なる第2の放射線透過画像群とを収集し記憶する
画像収集記憶装置と、第1と第2の画像群間の演
算により特定臓器のサブトラクト像を作成する造
影像作成装置と、該作成したサブトラクト像を表
示する画像表示装置とを備えた放射線診断装置に
おいて、前記画像収集記憶装置における第1と第
2の画像群より各々に選んだ2枚の画像につき、
各画素毎の負成分濃度差信号の和の最小値を求め
ることにより、位置ずれのない一組の画像を検出
し、該画像を前記造影像作成装置に出力する自動
画像位置補正装置を設け、位置ずれによる偽像を
除去したサブトラクト像を表示することを特徴と
するものである。
〔発明の実施例〕
以下、本発明の一実施例を図面を参照しながら
説明する。
第2図は本発明の一実施例である放射線診断装
置を示すブロツク図である。第2図において、1
で示すのは、発生するX線の強度を調節可能なX
線発生器であり、2は該X線発生器1から発生さ
れたX線を照射される被検体である。3で示すの
は、前記被検体2を透過して得られるX線透過像
を電気信号に変換する2次元検出器(例えばイメ
ージインテンシフアイアと光学系と撮像管)であ
る。4は前記2次元検出器3より出力されるアナ
ログ信号をデジタル信号に変換するA/D変換器
である。5で示すのは、造影剤が被検者の診断目
的部位に混入する以前のマスク像よりなる第1の
放射線透過画像群と、造影剤が前記診断目的部位
に到達した後の造影像よりなる第2の放射線透過
画像群とを収集し記憶する画像収集記憶装置であ
る。6は前記画像収集記憶装置5の機器番号とア
ドレス等を制御するメモリコントローラである。
7は画像処理装置であり、当該画像処理装置
は、前記画像収集記憶装置5に記憶されている画
像又は画像間の演算を行ない、表示用画像記憶装
置8若しくは画像収集記憶装置5又はシステムコ
ントローラ12にコントロールデータを送ること
ができるものである。9は表示用画像記憶装置8
のアドレス等の制御を行なう表示用メモリコント
ローラである。10は表示用画像記憶装置8から
出力されるデジタル信号をアナログ信号に変換す
るD/A変換器である。11で示すのは、前記
D/A変換器10から出力されるアナログ信号に
基づき被検体2の画像を表示する画像表示装置で
ある。
12で示すのは、X線発生器1、造影剤注入装
置14、2次元検出器3、メモリコントローラ
6、画像処理装置7、表示用メモリコントローラ
9等の制御を行なう中央演算処理装置を含むシス
テムコントローラである。13で示すのは、前記
システムコントローラ12内で予めプログラムさ
れている被検体2内の診断目的部位に応じた撮影
条件と画像処理モードとシーケンス等を選択する
選択キイを有する操作パネルである。14で示す
のは、システムコントローラ12の制御を受けて
所定時期に所定量の血管造影剤を被検体2内に自
動注入する造影剤注入装置である。
次に、前記画像処理装置7について詳述する。
第3図は該画像処理装置の一例を示すブロツク
図である。第3図に示すように、画像処理装置7
は、位置補正装置20と、造影像作成装置21
と、画像強調装置22と、セレクタ23とより成
る。位置補正装置20は、画像データバス及び
を介して伝送される画像データとシステムコン
トローラ12から出力される制御信号C1とを入
力し、位置補正係数を求め、それを制御データバ
スC2を介してシステムコントローラ12に出力
するように構成されている。造影像作成装置21
は、画像データバス及びから画像データと制
御信号C1とを入力し、対数変換・差分・その他
補正等を実施して造影像をデータバスに出力す
る。画像強調装置22は、造影像作成装置21と
同じように構成されているが、演算がウインドウ
イング・フイルタリング等の画像強調である。セ
レクタ23は、データバス上の画像データを入
力し、画像収集記憶装置5につながるデータバス
又は表示用メモリ8につながるデータバスに
データを出力するかを選択する。
第4図a〜cは上記実施例を用いた場合の動作
例をフローチヤートで示したものである。
先ず、ステツプ30では、造影剤が被検体の診断
目的部位に混入する以前の放射線透過画像(画像
1)における位置補正領域(以下ROI1という)
を設定する。ここではxを水平方向としyを垂直
方向として、ROI1はx方向にXからΔX画素y
方向にYからΔY画素の範囲とする。ステツプ32
に至る。ステツプ32では、造影剤が混入している
放射線透過画像(画像2)における位置補正ROI
(以下ROI2という)に関してROI1に対する初
期値設定用の値(I,J)を定める。ここでは、
先ずx方向にI=−1画素、y方向にJ=−1画
素とする。ステツプ33に到る。ステツプ33では、
画像2におけるROI2を設定する。ここではROI
2の範囲をx方向へX+IからΔX画素y方向へ
Y+JからΔY画素と設定する。ステツプ34に到
る。
ステツプ34では、画像1におけるROI1のアド
レスC1と画像2におけるROI2のアドレスC2
を設定する。ステツプ35に至る。ステツプ35で
は、画像入力行の初期値を設定する。ここでは入
力行の初期値をy=1と設定する。ステツプ36に
至る。
ステツプ36では、画像のパラレル入力を示す。
C1はアドレス信号によつて読み出されたROI1
内の画像データを示し、(t,u)はROI1内の
x,y座標を示す。
ここでは、t=X、X+1、…,X+ΔX−1 u=Y+y−1 となる。
一方、C2はアドレス信号によつて読み出され
たROI2内の画像データを示し、(v,w)は
ROI2内のx,y座標を示す。
ここでは、v=X+I、X+I+1、…、X+
I+ΔX−1 w=Y+J+y−1 となる。
このようにして、ROI1及びROI2におけるx
方向について所定行の一画素毎の画像データが入
力される。ステツプ37に至る。
ステツプ37では、ステツプ36からのROI1及び
ROI2における対応する一画素毎の画像データ間
の濃度差を示す。即ち、 濃度差S(x,y)=C1(t,u)−C2(v,w) となる。ステツプ38に至る。
ステツプ38では、ステツプ37からの濃度差信号
のうち負成分(すなわちサブトラクト信号の逆極
性成分)を判別する。そしてS(x,y)≦0に該
当するものについてはステツプ39に至る。
ステツプ39では、ステツプ38で選別された負成
分濃度差信号の和が求められる。この場合におけ
る負成分濃度差信号の和をε−(I,J)で表わ
すと、 ε-(I,J)=〓Xx=1 S(x,y) となる。
ステツプ40に至る。
ステツプ40では、画像入力行の繰り上げを行な
う。そしてステツプ41で画像入力行がROI内であ
るか否かを判定する。もし画像入力行yが、y≧
ΔY+1でない<NO>ならば、未だROI内であ
るので、ステツプ36に移りその行yの画像のパラ
レル入力が行なわれ、以下同様のステツプを経て
ステツプ41に至る。一方、画像入力行yが、y≧
ΔY+1である<YES>ならば、すでにROI外で
あるので、ステツプ42に移る。ステツプ42では、
以上のような過程を経てROIの全範囲に渡つて検
出した負成分濃度差信号の和ε-(I,J)が記憶
される。ステツプ43に至る。
ステツプ43では、ROI2に関してROI1に対す
る初期値設定用の値(I,J)について、x方向
のI値の繰り上げを行なう。そしてステツプ44で
は、前記I値の判定を行ない、I≧2でない<
NO>ならば位置補正領域内であり、ステツプ33
に移り、ROI2を設定する。その後は同様のステ
ツプを経てステツプ44に至る。一方、I≧2であ
る<YES>ならば位置補正領域外となり、ステ
ツプ45に移る。ステツプ45では、y方向のJ値の
繰り上げを行なう。そしてステツプ46では、前記
J値の判定を行ない、J≧2でない<NO>なら
ば位置補正領域内であり、ステツプ33に移り、
ROI2を設定する。その後は同様のステツプを経
てステツプ46に至る。一方、J≧2である<
YES>ならば位置補正領域外となり、ステツプ
47に移る。尚、この段階でIを初期値に戻す。
ステツプ47では、以上のようにI=−1,0,
1及びJ=−1,0,1の位置補正範囲内におい
て、ROIの位置を平行移動して、繰り返し(9通
り)求めたε-(I,J)の最小値を示す。この際
の最小値(I,J)=(IM,JM)を求める。ステツ
プ48に至る。
ステツプ48では、(IM,JM)の組が2個以上あ
るか否かを判定する。2個以上ある<YES>と
きには正しい位置補正値が判別不能であるので、
ステツプ55に移り位置補正ROIの設定値を変更す
る。例えば、X=X+A、Y=Y+Bとする。そ
して、ステツプ31に移り、ROI1を設定し、以後
は同様のステツプを経て負成分濃度差信号の和を
求め、ステツプ48に至る。一方、最小値(IM
JM)の組が2個以上ない<NO>ときにはステツ
プ49に至り、ステツプ49ではその(IM,JM)値を
位置補正係数(P,Q)としてそれぞれ加算記憶
される。ステツプ50に至る。ステツプ50では記憶
された前記(IM,JM)値がIM及びJMがいずれも零
であるか否かを判別する。もし、IM=JM=0でな
い<NO>ならばステツプ54に移り、ROIの初期
値を(X,Y)を(X+IM、Y+JM)に変更し、
ステツプ31に移つてROI1を設定する。以後は同
様のステツプを経て(IM,JM)のまわりで更に負
成分濃度差信号の和を求め、ステツプ50に至る。
一方、IM=JM=0である<YES>ならばステツプ
51に移り、ステツプ51では、ステツプ49で求めら
れた位置補正係数(P,Q)から位置補正値を算
出しアドレス信号として送り、画像1のアドレス
座標にしての画像2のアドレス座標を指定し、造
影中の画像位置を補正する。これは、IM=JM=0
を示したときが位置ずれのない状態を表わすもの
だからである。ステツプ52に至る。
ステツプ52では、画像1と画像2の対数変換し
てその差をとり、サブトラクト像を生成する。ス
テツプ53に至る。
ステツプ53では、その生成したサブトラクト像
を表示用メモリに出力する。
以上、本発明の一実施例について説明したが、
本発明は前記実施例に限定されるものではなく、
本発明の要旨を変更しない範囲内で適宜に変形し
て実施することができることはいうまでもない。
〔発明の効果〕
以上詳述したように本発明の放射線診断装置
は、造影剤が被検者の診断目的部位に混入する以
前の第1の画像群と、造影剤が前記診断目的部位
に到達した後の第2の画像群より各々に選んだ2
枚の画像につき、各画素毎の負成分濃度差信号の
和の最小値を求めることにより、位置ずれのない
一組の画像を検出し、該画像を前記造影剤作成装
置に出力する自動画像位置補正装置を設けること
によつて、位置ずれによる偽像を除去したサブト
ラクト像を表示できるので、正確且つ効率のよい
診断に貢献でき得るものである。
【図面の簡単な説明】
第1図a,b,c,d,eは本発明の前提とな
る原理を示す模式図、第2図は本発明の一実施例
を示すブロツク図、第3図は本発明の要部となる
画像処理装置の一例を示すブロツク図、第4図
a,b,cは本発明の一実施例における動作例を
フローチヤートで示したものである。 5…画像収集記憶装置、7…画像処理装置、1
1…画像表示装置、20…位置補正装置、21…
造影像作成装置、22…画像強調装置。

Claims (1)

  1. 【特許請求の範囲】 1 造影剤が被検体の診断目的部位に混入する以
    前のマスク像よりなる第1の放射線透過画像群
    と、造影剤が前記診断目的部位に到達した後の造
    影像よりなる第2の放射線透過画像群とを収集し
    記憶する画像収集記憶装置と、第1と第2との画
    像群間の演算により特定臓器のサブトラクト像を
    作成する造影像作成装置と、該作成したサブトラ
    クト像を表示する画像表示装置とを備えた放射線
    診断装置において、前記第1の画像中に設定され
    た第1の位置補正領域ROI1、第2の画像中に設
    定された第2の位置補正領域ROI2における各画
    素間の差分を求め、この差分値のうち逆極性成分
    のみの和を求めると共に、前記ROI2を平行移動
    して各位置における画素間の差分の逆極性成分の
    和を求め、各逆極性成分の和のうち、最小となる
    ROI2の前記ROI1からのずれ成分を位置補正値
    とし、これに基づき画像2の位置を補正する自動
    画像位置補正装置を設けたことを特徴とする放射
    線診断装置。 2 前記逆極性成分の和の最小値が複数存在する
    ときには前記第1の位置補正領域ROI1を平行移
    動して前記演算を繰り返す請求項1記載の放射線
    診断装置。
JP58025230A 1983-02-16 1983-02-16 放射線診断装置 Granted JPS59151941A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP58025230A JPS59151941A (ja) 1983-02-16 1983-02-16 放射線診断装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP58025230A JPS59151941A (ja) 1983-02-16 1983-02-16 放射線診断装置

Publications (2)

Publication Number Publication Date
JPS59151941A JPS59151941A (ja) 1984-08-30
JPH0227873B2 true JPH0227873B2 (ja) 1990-06-20

Family

ID=12160166

Family Applications (1)

Application Number Title Priority Date Filing Date
JP58025230A Granted JPS59151941A (ja) 1983-02-16 1983-02-16 放射線診断装置

Country Status (1)

Country Link
JP (1) JPS59151941A (ja)

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4335427A (en) * 1980-04-21 1982-06-15 Technicare Corporation Method of selecting a preferred difference image
DE3043703C2 (de) * 1980-11-19 1983-01-20 Siemens AG, 1000 Berlin und 8000 München Röntgendiagnostikeinrichtung zur Erstellung von Subtraktionsbildern

Also Published As

Publication number Publication date
JPS59151941A (ja) 1984-08-30

Similar Documents

Publication Publication Date Title
US7761136B2 (en) Medical image processing apparatus for scanning based on a set three-dimensional region of interest
US4674046A (en) Method and apparatus for obtaining three dimensional tomographic images by interpolation of a plurality of projection slice data bind for obtaining projection data for a chosen slice
JP4340334B2 (ja) 対象と対象の3次元表現との間の変換を決定する方法及びその装置
CN108703764B (zh) 血管造影方法、装置、系统、设备及存储介质
JP5379960B2 (ja) 3次元画像処理装置及び再構成領域指定方法
EP0087470A4 (en) DIGITAL FLUOROGRAPHIC SYSTEM AND METHOD.
Van Tran et al. Flexible mask subtraction for digital angiography
US4802093A (en) X-ray image-processing apparatus utilizing grayscale transformation
CN111184523A (zh) 基于dr设备的三维图像重建方法及系统
CN112120722A (zh) X射线断层合成装置、图像处理装置以及计算机可读记录介质
US7377691B2 (en) Radiographic apparatus that removes time lag by recursive computation
JP2004236910A (ja) 医用画像の貼り合わせ装置
JP3540914B2 (ja) X線撮影装置
JPH0572795B2 (ja)
JPH031870B2 (ja)
KR20060061249A (ko) 선량 평가 방법 및 x선 ct 장치
JPH07148155A (ja) コンピュータ断層撮影装置
JP3641499B2 (ja) ディジタルx線装置
US7760856B2 (en) Radiographic apparatus and radiation detection signal processing method
JPH0227873B2 (ja)
JP2005064706A (ja) 放射線撮像装置および放射線検出信号処理方法
US20040228444A1 (en) X-ray diagnostic apparatus with image computer for direction filtering
JPH0879628A (ja) X線診断装置
Tommasini et al. A deterministic approach to automated stenosis quantification
JPH08131428A (ja) ディジタルx線撮影装置