JPS60251475A - サブトラクシヨン画像の階調処理最適化方法および装置 - Google Patents

サブトラクシヨン画像の階調処理最適化方法および装置

Info

Publication number
JPS60251475A
JPS60251475A JP59107912A JP10791284A JPS60251475A JP S60251475 A JPS60251475 A JP S60251475A JP 59107912 A JP59107912 A JP 59107912A JP 10791284 A JP10791284 A JP 10791284A JP S60251475 A JPS60251475 A JP S60251475A
Authority
JP
Japan
Prior art keywords
signal
image
value
minimum
subtraction
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
JP59107912A
Other languages
English (en)
Inventor
Kazuo Shimura
一男 志村
Nobuyoshi Nakajima
中島 延淑
Masamitsu Ishida
石田 正光
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.)
Fujifilm Holdings Corp
Original Assignee
Fuji Photo Film 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 Fuji Photo Film Co Ltd filed Critical Fuji Photo Film Co Ltd
Priority to JP59107912A priority Critical patent/JPS60251475A/ja
Priority to DE8585104941T priority patent/DE3583065D1/de
Priority to EP85104941A priority patent/EP0163903B1/en
Publication of JPS60251475A publication Critical patent/JPS60251475A/ja
Priority to US07/116,196 priority patent/US5028781A/en
Pending legal-status Critical Current

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/40—Picture signal circuits
    • H04N1/407—Control or modification of tonal gradation or of extreme levels, e.g. background level
    • H04N1/4072—Control or modification of tonal gradation or of extreme levels, e.g. background level dependent on the contents of the original
    • H04N1/4074—Control or modification of tonal gradation or of extreme levels, e.g. background level dependent on the contents of the original using histograms
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/20—Measuring radiation intensity with scintillation detectors
    • G01T1/2012—Measuring radiation intensity with scintillation detectors using stimulable phosphors, e.g. stimulable phosphor sheets
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

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

Description

【発明の詳細な説明】 (発明の分野) 本発明は放射線画像のサブトラクション処理、詳細には
蓄積性螢光体シートを用いて行なう放射線画像のデジタ
ルサブトラクション処理において、濃度範囲、コントラ
ストが常に一定となるようにサブトラクション画像を階
調処理する方法およびその装置に関するものである。
(発明の技術的背景および先行技術) 従来より放射線画像のデジタルサブトラクションが公知
となっている。この放射線画像のデジタルサブトラクシ
ョンとは、異なった条件で撮影した2つの放射線画像を
光電的に読み出してデジタル画像信号を得た後、これら
のデジタル画像信号を両画像の各画素を対応させて減算
処理し、放射線画像中の特定の構造物の画像を形成する
ための差信号を得る方法であり、このようにして得た差
信号を用いて特定構造物のみが抽出された放射線画像を
再生することができる。
このサブトラクション処理には、基本的に次の2つの方
法がある。即ち、造影剤注入により特定の構造物が強調
された放射線画像の画像信号から、造影剤が注入されて
いない放射線画像の画像信号を引き算(サブトラクト)
づ−ることによって特定の構造物を抽出するいわゆる時
間サブトラクション処理と、同一の被写体に対して相異
なるエネルギー分布を有する放射線を照射し、それによ
り特定の構造物が特有の放射線エネルギー吸収特性を有
することを利用して特定構造物が異なる画像を2つの放
射線画像間に存在せしめ、その後この2つの放射線画像
の画像信号間で適当な重みづけをした上で引き算くサブ
トラクト)を行ない特定の構造物の画像を抽出するいわ
ゆるエネルギーサブトラクション処理である。
このサブトラクション処理は特に医療用のX線写真の画
像処理において診断上きわめて有効な方法であるため、
近年大いに注目され、電子工学技術を駆使してその研究
、開発が盛んに進められている。この技術は、特にデジ
タルサブトラクション処理(通常Digital Ra
dioqraphy)と呼ばれ、DRと略称されている
。
さらに最近では例えば特開昭58−163340号公報
に示されるように、きわめて広い放射線露出域を有する
蓄積性螢光体シートを使用し、これらの螢光体シートに
前述のように異なった条件で同一の被写体を透過した放
射線を照射して、これらの螢光体シートに特定構造物の
画像情報が異なる放射線画像を蓄積記録し、これらの蓄
積画像を励起光による走査により読み出してデジタル信
号に変換し、これらデジタル信号により前記デジタルサ
ブトラクションを行なうことも提案されている。上記蓄
積性螢光体シートとは、例えば特開昭55−12429
号公報に開示されているように放射線(X線、α線、β
線、γ線、紫外線等)を照射するとその放射線エネルギ
ーの一部を螢光体中に蓄積し、その後可視光等の励起光
を照射すると蓄積された放射線エネルギー量に応じて螢
光体が輝尽発光を示すもので、きわめて広いラチチュー
ド(露出域)を有し、かつ著しく高い解像力を有するも
のである。したがって、この螢光体シ−トに蓄積記録さ
れた放射線画像情報を利用して前記デジタルザブトラク
ションを行なえば、放射線量が変動しても常に十分な画
像情報を得ることができ、診断能の高い放射線画像を得
ることができる。
しかしながら、前述のような蓄積性螢光体シートを用い
てエネルギーザブトラクションを行なう場合、サブトラ
クションによって得られた差信号をそのまま用いて出力
画像を形成すると、各出力画像によって濃度範囲(ある
いは輝度範囲)やコントラストがまちまちになることが
多い。これは前述した蓄積性螢光体シートの感度のバラ
ツキや、画@読取器の感度バラツキ、撮影放射線エネル
ギー量のバラツキ等に起因して、サブトラクションに供
する2つの画像の画像信号が変動することによるもので
ある。サブトラクションによって得られる差信号は微小
なものであるので、差信号を増幅するような階調処理を
行なうことが有効であるが、その時上記のように画像信
号が変動づ゛ると、差信号におけるバラツキ成分は本来
の信号成分に比べてかなり大きなものとなり、濃度範囲
やコントラストの大きな変動となって表われる。上記の
ようにサブトラクション画像の濃度(あるいは輝度)や
コントラストがまちまちであると、複数のサブトラクシ
ョン画像を比較して診断を下す場合等において適正な診
断が妨げられたり、また人間の視覚特性の良好な領域か
ら画像の濃度がはずれてしまうことが指摘されている。
なお、蓄積性螢光体シートを用いるエネルギーサブトラ
クションにおいて、サブトラクションに供される放射線
画像(すなわち特定構造物の画像情報が互いに異なる放
射線画像)を得るには、放射線源の管電圧を変えて放射
線エネルギーを変化させる他に、放射線源と被写体との
間に放射線のエネルギー分布を変化させるフィルタを出
し入れする等、その他の公知の方法が用いられてもよい
し、さらには本出願人による特願昭57−193765
号に開示されているように、蓄積性螢光体シートとフィ
ルタとの積層体等を使用して、1回の放射線撮影によっ
て得るようにしてもよい。
(発明の目的) 本発明は上記のような事情に鑑みてなされたものであり
、濃度範囲、コントラストが一定で診断能の良いエネル
ギーサブトラクション画像を得ることができる、サブト
ラクション画像の階調処理最適化方法および装置を提供
することを目的とするものである。
(発明の構成) 本発明のサブトラクション画像の階調処理最適化方法は
、前述したように蓄積性螢光体シートを用いて行なうエ
ネルギーサブトラクション処理において、サブトラクシ
ョンによって得た差信号の信号最大値、信号最小値をめ
、次に上記差信号を、上記信号最大値、最小値を予め決
められた所定の出力画像濃度最大値、最小値にそれぞれ
対応させて階調処理するようにしたものである。
上記差信号の信号最大値、信号最小値は、差信号を順次
比較してめることもできるし、あるいは該差信号のヒス
トグラムをめてこのヒストグラムからめることもできる
が、特に順次比較による方法はヒストグラムを作ってか
らめる方法に比べ、簡単にハードウェアのみで構成する
ことができ、また演算スピードも速いという利点を有し
ている。
上述のようにしてめた差信号の信号最大値、最小値を所
定の出力画像濃度最大値、最小値に対応させるためには
、差信号を階調処理前に一律に信号変換してもよいし、
また階調変換テーブルを補正してもよい。
上記方法を実施する本発明のサブトラクション画像の階
調処理最適化装置は、差信号を一律に信号変換する場合
には、前述したような励起光走査と輝尽発光光の光電的
読み出しにより、蓄積性螢光体シートの放射線画像のデ
ジタル信号を得る画像読取手段と、この画像読取手段が
読み取った、互いに特定構造物の画像情報が異なる被写
体のデジタル画像信号とを、対応する画素間で減算して
前記特定構造物の画像を形成する差信号を得るサブトラ
クション演算手段と、この差信号に階調変換テーブルに
基づいて階調処理を施す画像処理手段と、前記差信号の
信号最大値と信号最小値をめる最大最小演算手段と、上
記信号最大値と最小値が予め決められた所定の出力画像
濃度最大値、最小値にそれぞれ対応するように、前記階
調処理を受ける前の前記差信号を一律に変換する信号変
換回路とから構成される。
また上記信号変換回路に代えて、前記差信号の信号最大
値、最小値が予め決められた所定の出力画像濃度最大値
、最小値にそれぞれ対応するように、前記階調変換テー
ブルを補正する階調変換テーブル補正回路を用いれば、
差信号はそのまま用いて階調処理を最適化できる。
(実施態様) 以下、図面に示す実施態様に基づいて本発明の詳細な説
明する。
第1図は2枚の蓄積性螢光体シートA、Bに同一の被写
体1を透過したX線2を、それぞれエネルギーを変えて
照射する状態を示す。すなわち第1の蓄積性螢光体シー
トAに被写体1のX線透過像を蓄積記録し、次いで短時
間内で蓄積性螢光体シートA、Bを素早く取り替えると
同時に、X線源3の管電圧を変えて、透過X線のエネル
ギーが異なる被写体1のX線画像を蓄積性螢光体シート
Bに蓄積記録する。このとき蓄積性螢光体シートAとB
とで被写体1の位置関係は同じとする。
このようにして、少なくとも一部の画像情報が異なる2
つの放射線画像を2枚の蓄積性螢光体シートA、Bに蓄
積記録する。次にこれら2枚の蓄積性螢光体シートA、
Bから、第2図に示すような画像読取手段によってX線
画像を読み取り、画像を表わすデジタル画像信号を得る
。先ず、蓄積性螢光体シートAを矢印Yの方向に副走査
のために移動させながら、レーザー光源10からのレー
ザー光11を走査ミラー12によってX方向に主走査さ
せ、螢光体シートAから蓄積X線エネルギーを、蓄積記
録されたX線画像にしたがって輝尽発光光13として発
散させる。輝尽発光光13は透明なアクリル板を成形し
て作られた集光板14の一端面からこの集光板14の内
部に入射し、中を全反射を繰返しつつフォトマル15に
至り、輝尽発光光13の発光量が画像信@Sとして出力
される。この出力された画像信号Sは増幅器とA/D変
換器を含む対数変換器16により対数値(1ops)の
デジタル画像信号10QSAに変換される。このデジタ
ル画像信号10QsAは例えば磁気ディスク等の記憶媒
体17に記憶される。次に、全く同様にして、もう1枚
の蓄積性螢光体シートBの記録画像が読み出され、その
デジタル画像信号logsBが同様に記憶媒体17に記
憶される。
次に、上述のようにして得られたデジタル画像信号+0
QsA、IoqSBを用いてサブトラクション処理を行
なう。第3図は本発明方法の第1実施態様による階調処
理最適化方法を適用して行なわれるサブトラクション処
理の流れを示している。まず前記記憶媒体17内の画像
ファイル17Aと、画像ファイル17Bからそれぞれ、
前記デジタル画像信号1ogsハ、10g5Bが読み出
され、サブトラクション演算回路18に入力される。こ
のサブトラクション演算回路18は、上記2つのデジタ
ル画像信号1OGsAと1oqsBを適当な重みづけを
した上で対応する画素毎に減算し、デジタルの差信号 5sub=a @ l0QsA b ・ 1OQsB+
c(a、bは重みづけ係数、Cは概略一定濃度にするよ
うなバイアス成分である)をめる。この差信号5Sub
は−たん画像ファイル19に記憶されてから、画像処理
回路20に入力され、該画像処理回路20において後述
の信号変換回路20aを介して階調変換テーブル20b
に基づいて階調処理される。
階調処理を受けた差信号3sub’は、例えばCRT等
のディスプレイ装置や、走査記録装置等の再生記録装置
21に入力され、該差信号5SUb′によってサブトラ
クション画像が再生記録される。第4図はサブトラクシ
ョン画像再生記録システムの一例として、画像走査記録
装置を示すものである。感光フィルム30を矢印Yの副
走査方向へ移動させるとともにレーザービーム31をこ
の感光フィルム30上にX方向に主走査させ、し−ザー
ビーム31をA10変調器32により画像信号供給器3
3からの画像信号によって変調することにより、感光フ
ィルム30上に可視像を形成する。この変調用画像信号
として、前記差信号5sub’を使用すれば、デジタル
サブトラクション処理による所望の特定構造物のみの画
像を感光フィルム30上に再生記録することができる。
このようにして感光フィルム30上に記録再生されるエ
ネルギーサブトラクション画像は、2枚の蓄積性螢光体
シートA、Bの感度のバラツキや撮影X線エネルギーの
バラツキなどにより前記差信号3subが変動するため
に、各画像により濃度範囲、コントラストがまちまちに
なってしまうことが多い。そこで第3図に示されるよう
に上記差信号3SUbをヒストグラム演算回路22に入
力し、該差信号5subのヒストグラムをめる。
そして該ヒストグラム演算回路22は第5図に示すよう
に、このヒストグラムから一例として上記差信号3su
bの最大値Smax、最小値Sm1nをめてその情報を
前記信号変換回路20aに入力する。この信号変換回路
20aはこれら最大値Smax 、最小値51nがそれ
ぞれ常に一定の値Cmax 、 Cm1nをとるように
差信号5subを変換する。すなわち例えば 5sub−+p−8sub+q p−(Cmax −Cmin ) / (SIIlaX
 −8lllin )q = (Cmax + Cm1
n ) / 2− (Smax −+−8m1n ) 
/2なる変換を行なう。
変換された信号p−3Sub+qは前述した階調変換テ
ーブル20bに基づいて階調処理されるが、ここで上記
一定値Cmax 、 Cm1nは第5図に示すように、
サブトラクション画像の濃度最大値DIIla×、濃度
最小値Ominを得る信号となっている。したがって上
記変換された信号p−8sub+qによるザブトラクシ
ョン画像は常に一定の濃度範囲、コントラストを有する
ようになる。
第6図は本発明の第2実施態様による階調処理最適化方
法を適用して行なわれるエネルギーサブトラクション処
理の流れを示すものである。本実施態様においては、前
記第1実施態様において用いられたヒストグラム演算回
路220代わりに、最大最小演算手段として順次比較演
算回路122a、122bが用いられ、また信号変換回
路20aに代わるものとして階調変換テーブル補正回路
120が用いられている。
上記順次比較演算回路122aはサブトラクション演算
回路18から送られた最初の前記差信号5subを最大
値として記憶し、次に送られて来る差信号3subと上
記記憶された差信号3subとの大小を比較する。そし
て信号値が大きい方の差信号3subを最大値として記
憶し、以下送られて来る差信号3subに対して同様の
処理を順次施す。この処理が最後の差信号5subまで
順次行なわれると、差信号3subの全5領域における
真の最大値Smax’がめられる。また使方の順次比較
演算回路122bは、同様に差信号5subを比較して
小さい方の差信号5subを最小値として記憶すること
により、差信号5subの全領域における真の最小値3
m1n ’をめる。
なおこのような順次比較法によって差信号3subの最
大!aSIllax゛、最小値3m1n +をめる場合
、上記のように最初に送られて来る差信号5Subをひ
とまず最大値3max ’ 、最小値5m1n’として
設定する他に、想定される最大値Smax’よりも十分
に小さい任意の信号値を予め最大値として記憶させてお
き、また想定される最小値3m1n′よりも十分に大き
い任意の信号値を予め最小値として記憶させておいても
よい。
上記のようにしてめられた真の最大値3max’、最小
値3m1n ’はそれぞれ減算回路121a、加算回路
121bに送られる。減算回路121aにおいては真の
最大値3max、’から所定信号レベルSΔaが差し引
かれ、(3max ’−8Δa)が階調処理最適化用の
最大値3 maxとして出力される。
また加算回路121bにおいては上記真の最小値3m1
n’に所定信号レベルSΔbが加えられ、(3min 
’ +3Δb〉が階調処理最適化のための最小値3 m
inとして出力される。階調変換テーブル補正回路12
0は上記最大値smax、最小値Sminを受けて、第
7図に示されるように、サブトラクション画像の濃度最
大値[Nnax 、 il1度最小値D minに上記
最大値5Illax1最小値S minが対応するよう
に階調変換テーブル20bを補正する(図中破線表示)
。このように階調変換テーブル20bを補正すると、前
記第1実施態様におけるのと同様に、サブトラクション
画像の濃度範囲、コントラストが一定に保たれるように
なる。
第7図に示されるように差信号3subには、そのレベ
ル最大側、最小側において、信号ノイズや撮影時の散乱
放射線によるノイズ成分のみが含まれる。上記具の最大
値Smax’から差し引く信号レベルSΔa1および真
の最小値31n ’に加える信号レベルSΔbは上述の
ノイズ成分に対応する値に設定されている。したがって
階調処理最適化用の最大値3maX、最小値3 min
としてそれぞれ(Smax ’−8Δa)、(31Il
in ’ +SΔb)の値を使用することにより、前記
信号ノイズや散乱放射線の影響を除去して、サブトラク
ション画像の記録媒体のダイナミックレンジをより有効
に利用することが可能になる。
また別の方法として、真の最大値Smax’から数えで
あるデータ数だけ小さいところの値および真の最小値3
m1n’から数えであるデータ数だけ大きいところの値
をそれぞれ階調処理最適化用の最大値Smaxおよび最
小値Sm1nとして使用するようにしてもよい。
以上2枚のシートA、Bを使用する実施態様に゛ついて
説明したが、その他3枚のシートにそれぞれ異なるエネ
ルギーで放射線撮影し、これらシートから得られる画像
信号を減算処理して(例えばaxlogs +bxlo
gS C× l0Qsこ八 B− +dなど。但し、ここでa、bおよびCは重みづけ係数
、dは概略一定濃度にするようなバイアス成分である。
)サブトラクション画像を得ることもできるが、本発明
はこのような場合にも適用可能である。また差信号の最
大値、最小値は、ある特定の画像部分く例えば関心領域
)のみについてめてもよい。
なお先に述べたように、サブトラクションに供される放
射線画像(すなわち特定構造物の画像情報が互いに異な
る放射線画像)を得るには、前記実施態様におけるよう
に放射線源の管電圧を変えて放射線エネルギーを変化さ
せる他に、放射線源と被写体との間に放射線のエネルギ
ー分布を変化させるフィルタを出し入れする等、その他
の公知の方法が用いられてもよいし、さらには本出願人
による特願昭57−193765号に開示されているよ
うに蓄積性螢光体シートとフィルタとの積層体等を使用
して、1回の放射線撮影によって得るようにしてもよい
。
(発明の効果) 以上詳細に説明した通り本発明によれば、エネルギーサ
ブトラクション画像の全体的な濃度、コントラストを各
画像間で共通に保つことができるから、極めて診断能に
優れたサブトラクション画像を得ることが可能となる。
そして濃度、コントラストを所定の範囲に保てるから、
サブトラクション画像を出力する記録媒体のダイナミッ
クレンジを最も有効に活用することが可能になる。また
本発明方法は電気回路により自動的に実施可能であるか
ら、特にサブトラクション演算処理を連続的に行ない、
サブトラクション画像をCRTを用いてリアルタイムで
観察する場合に特に有効である。
【図面の簡単な説明】
第1図は本発明方法における放射線画像の蓄積記録ステ
ップを示す説明図、 第2図は上記蓄積記録がなされた蓄積性螢光体シートか
らの放射線画像情報読取りを説明する概略図、 第3図は本発明方法の第1実施態様方法を適用して階調
処理を行なうエネルギーサブトラクション処理の概要を
説明する概略図、 第4図はサブトラクション画像の再生記録システムの一
例を示す概略図、 第5図は本発明の第1実施態様による階調処理最適化方
法を説明する説明図、 第6図は本発明方法の第2実施態様方法を適用して階調
処理を行なうエネルギーサブトラクション処理の概要を
説明する概略図、 第7図は本発明の第2実施態様による階調処理最適化方
法を説明する説明図である。 1・・・被写体 2・・・X線 3・・・X線源 10・・・レーザー光源11・・・レ
ーザー光 12・・・走査ミラー13・・・輝尽発光光
 15・・・フォトマル16・・・対数変換器 18・・・サブトラクション演算回路 20・・・画像処理回路 20a・・・信号変換回路2
0b・・・階調変換テーブル 22・・・ヒストグラム演算回路 120・・・階調変換テーブル補正回路122a、12
2b・・・順次比較演算回路A、8・・・蓄積性螢光体
シート I o gSA 、I o’asB・・・デジタル画像
信号5sub・・・デジタル画像信号の差信号5sub
’・・・階調処理された差信号g+図 1!3図 1!5図

Claims (14)

    【特許請求の範囲】
  1. (1) 2枚以上の蓄積性螢光体シートに、放射線エネ
    ルギー吸収特性が他とは異なる特定構造物を含む被写体
    を透過した放射線を照射して、これらの螢光体シートに
    前記被写体の特定構造物の画像情報が互いに異なる放射
    線画像を蓄積記録し、これらの螢光体シートに励起光を
    走査して前記放射線画像を輝尽発光光に変換し、この輝
    尽発光光の発光量を光電的に読み出してデジタル画像信
    号に変換し、各画像の対応する画素間でこのデジタル画
    像信号の減算を行なって放射線画像の特定構造物の画像
    を形成する差信号を得、その後この差信号に階調処理を
    施すエネルギーサブトラクションの階調処理において、
    前記差信号の信号最大値、信号最小値をめ、次に前記差
    信号を、前記信号最大値、信号最小値を予め決められた
    出力画像濃度最大値、濃度最小値のそれぞれに対応させ
    て階調処理することを特徴とするサブトラクション画像
    の階調処理最適化方法。
  2. (2) 前記信号最大値を、前記差信号の全領域におけ
    る真の最大値から所定信号レベルを差し引いた値とし、
    前記信号最小値を、前記差信号の全領域における真の最
    小値に所定信号レベルを加えた値とすることを特徴とす
    る特許請求の範囲第1項記載のサブトラクション画像の
    階調処理最適化方法。
  3. (3) 前記信号最大値を、前記差信号の全領域におけ
    る真の最大値から数えであるデータ数だけ小さいところ
    の値とし、前記信号最小値を、前記差信号の全領域にお
    ける真の最小値から数えであるデータ数だ【プ大きいと
    ころの値とすることを特徴とする特許請求の範囲第1項
    記載のサブ1ヘラクション画像の階調処理最適化方法。
  4. (4) 前記差信号を順次比較することによって前記信
    号最大値、信号最小値をめることを特徴とする特許請求
    の範囲第1項から第3項いずれか1項記載のサブトラク
    ション画像の階調処理最適化方法。
  5. (5) 前記差信号を階調処理前に一律に信号変換する
    ことにより、該差信号の前記最大値、最小値を前記出力
    画像濃度最大値、濃度最小値に対応させることを特徴と
    する特許請求の範囲第1項から第4項いずれか1項記載
    のサブトラクション画像の階調処理最適化方法。
  6. (6) 前記階調処理に用いられる階調変換テーブルを
    修正することにより、前記差信号の前記最大値、最小値
    を前記出力画像濃度最大値、濃度最小値に対応させるこ
    とを特徴とする特許請求の範囲第1項から第4項いずれ
    か1項記載のサブトラクション画像の階調処理最適化方
    法。
  7. (7) 放射線画像が蓄積記録された蓄積性螢光体シー
    トに励起光を走査し、それによって前記蓄積性螢光体シ
    ートから発せられた輝尽発光光を光電的に読み出してデ
    ジタル画像信号に変換する画像読取手段と、この画像読
    取手段が読み取った、互いに特定構造物の画像情報が異
    なる被写体のデジタル画像信号とを、対応する画素間で
    減算して前記特定構造物の画像を形成する差信号を得る
    サブトラクション演算手段と、この差信号に階調変換テ
    ーブルに基づいて階調処理を施す画像処理手段と、前記
    差信号の信号最大値と信号最小値をめる最大最小演算手
    段と、前記信号最大値、最小値が予め決められた所定の
    出力画像濃度最大値、最小値にそれぞれ対応するように
    、前記階調処理を受ける前の前記差信号を一律に変換す
    る信号変換回路とからなるサブトラクション画像の階調
    処理最適化装置。
  8. (8) 前記最大最小演算手段が、前記差信号の全領域
    における真の最大値から所定信号レベルを差し引いた値
    を前記信号最大値とし、前記差信号の全領域におりる真
    の最小値に所定信号レベルを加えた値を前記信号最小値
    とするように構成されていることを特徴とする特許請求
    の範囲第7項記載のサブトラクション画像の階調処理最
    適化方法。
  9. (9) 前記最大最小演算手段が、前記差信号の全領域
    における真の最大値から数えであるデータ数だ【プ小さ
    いところの値を前記信号最大値とし、前記差信号の全領
    域における真の最小値から数えであるデータ数だけ大き
    いところの値を前記信号最小値とするように構成されて
    いることを特徴とする特許請求の範囲第7項記載のサブ
    トラクション画像の階調処理最適化方法。
  10. (10) 前記最大最小演算手段が、前記差信号を順次
    比較して該差信号の前記信号最大値、信号最小値をめる
    順次此較演算手段であることを特徴とする特許請求の範
    囲第7項から第9項いずれが1項記載のサブトラクショ
    ン画像の階調処理最適化装置。
  11. (11) 放射線画像が蓄積記録された蓄積性螢光体シ
    ートに励起光を走査し、それによって前記蓄積性螢光体
    シートから発せられた輝尽発光光を光電的に読み出して
    デジタル画像信号に変換する画像読取手段と、この画像
    読取手段が読み取った、互いに特定構造物の画像情報が
    異なる被写体のデジタル画像信号とを、対応する画素間
    で減算して前記特定構造物の画像を形成する差信号を得
    るサブトラクション演算手段と、この差信号に階調変換
    テーブルに基づいて階調処理を施す画像処理手段と、前
    記差信号の信号最大値と信号最小値をめる最大最小演算
    手段と、前記信号最大値、最小値が予め決められた所定
    の出力画像濃度最大値、最小値にそれぞれ対応するよう
    に、前記階調変換テーブルを補正する階調変換テーブル
    補正回路とからなるサブトラクション画像の階調処理最
    適化装置。
  12. (12) 前記最大最小演算手段が、前記差信号の全領
    域における真の最大値から所定信号レベルを差し引いた
    値を前記信号最大値とし、前記差信号の全領域における
    真の最小値に所定信号レベルを加えた値を前記信号最小
    値とするように構成されていることを特徴とする特許請
    求の範囲第11項記載のサブトラクション画像の階調処
    理最適化方法。
  13. (13) 前記最大最小演算手段が、前記差信号の全領
    域における真の最大値から数えであるデータ数だけ小さ
    いところの値を前記信号最大値とし、前記差信号の全領
    域における真の最小値から数えてあるデータ数だけ大き
    いところの値を前記信号最小値とするように構成されて
    いることを特徴とする特許請求の範囲第11項記載のサ
    ブトラクション画像の階調処理最適化方法。
  14. (14) 前記最大最小演算手段が、前記差信号を順次
    比較して該差信号の前記信号最大値、信号最小値をめる
    順次比較演算手段であることを特徴とする特許請求の範
    囲第11項から第13項いずれか1項記載のサブトラク
    ション画像の階調処理最適化装置。
JP59107912A 1984-04-23 1984-05-28 サブトラクシヨン画像の階調処理最適化方法および装置 Pending JPS60251475A (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP59107912A JPS60251475A (ja) 1984-05-28 1984-05-28 サブトラクシヨン画像の階調処理最適化方法および装置
DE8585104941T DE3583065D1 (de) 1984-04-23 1985-04-23 Schwaerzungsgrad-korrekturverfahren sowie vorrichtung fuer subtraktionsbilder.
EP85104941A EP0163903B1 (en) 1984-04-23 1985-04-23 Density correcting method and apparatus for subtraction image
US07/116,196 US5028781A (en) 1984-04-23 1987-11-03 Density correcting method and apparatus for subtraction image

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP59107912A JPS60251475A (ja) 1984-05-28 1984-05-28 サブトラクシヨン画像の階調処理最適化方法および装置

Publications (1)

Publication Number Publication Date
JPS60251475A true JPS60251475A (ja) 1985-12-12

Family

ID=14471205

Family Applications (1)

Application Number Title Priority Date Filing Date
JP59107912A Pending JPS60251475A (ja) 1984-04-23 1984-05-28 サブトラクシヨン画像の階調処理最適化方法および装置

Country Status (1)

Country Link
JP (1) JPS60251475A (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63225879A (ja) * 1987-03-14 1988-09-20 Fujitsu Ltd 階調補正処理方式

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS55105768A (en) * 1979-02-08 1980-08-13 Nec Corp Picture quantizing unit
JPS58201484A (ja) * 1982-05-20 1983-11-24 Hitachi Medical Corp Dfa画像処理装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS55105768A (en) * 1979-02-08 1980-08-13 Nec Corp Picture quantizing unit
JPS58201484A (ja) * 1982-05-20 1983-11-24 Hitachi Medical Corp Dfa画像処理装置

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63225879A (ja) * 1987-03-14 1988-09-20 Fujitsu Ltd 階調補正処理方式

Similar Documents

Publication Publication Date Title
US4747052A (en) Radiation image processing
CA1223980A (en) Method of adjusting radiation image read-out apparatus
US5048110A (en) Method and apparatus for automatically correcting subtraction image density
US5535289A (en) Method for reducing noise in energy subtraction images
JPH0362411B2 (ja)
JP3901770B2 (ja) エネルギーサブトラクション処理方法および装置
US4761739A (en) Density correcting method and apparatus for energy substraction image
JP2849964B2 (ja) 画像処理方法および装置
US5028781A (en) Density correcting method and apparatus for subtraction image
US7483556B2 (en) Energy subtraction processing method and apparatus
US5604780A (en) Method for processing a radiation image
JPS60250792A (ja) サブトラクシヨン画像の階調処理最適化方法および装置
JPS60224386A (ja) サブトラクシヨン画像の濃度補正方法および装置
US5029083A (en) Energy subtraction processing method using weighted average of plural stacked sheets for creating high energy image
US5237176A (en) Method of energy subtraction for radiation image and apparatus for carrying out the method
JPH0262074B2 (ja)
JPH11205682A (ja) エネルギーサブトラクション画像生成方法
JP2739385B2 (ja) 放射線画像読取条件及び/又は画像処理条件決定方法および装置
JP2824878B2 (ja) 放射線画像のエネルギーサブトラクション撮影処理方法
JPH067669B2 (ja) エネルギ−サブトラクシヨン画像の画質改善方法および装置
JPH0262075B2 (ja)
JPH0363877B2 (ja)
JPS60227735A (ja) サブトラクシヨン画像の自動濃度補正方法および装置
JPH0471395B2 (ja)
JP3165530B2 (ja) 放射線画像読取条件決定方法および装置