JPH0363878B2 - - Google Patents
Info
- Publication number
- JPH0363878B2 JPH0363878B2 JP59107911A JP10791184A JPH0363878B2 JP H0363878 B2 JPH0363878 B2 JP H0363878B2 JP 59107911 A JP59107911 A JP 59107911A JP 10791184 A JP10791184 A JP 10791184A JP H0363878 B2 JPH0363878 B2 JP H0363878B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- signal
- difference signal
- subtraction
- gradation processing
- 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
Links
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
- Conversion Of X-Rays Into Visible Images (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Closed-Circuit Television Systems (AREA)
Description
【発明の詳細な説明】
（発明の分野）
本発明は放射線画像のサブトラクシヨン処理、
詳細には蓄積性蛍光体シートを用いて行なう放射
線画像のデジタルサブトラクシヨン処理におい
て、濃度範囲、コントラストが常に一定になるよ
うにサブトラクシヨン画像を階調処理する方法お
よびその装置に関するものである。
詳細には蓄積性蛍光体シートを用いて行なう放射
線画像のデジタルサブトラクシヨン処理におい
て、濃度範囲、コントラストが常に一定になるよ
うにサブトラクシヨン画像を階調処理する方法お
よびその装置に関するものである。
（発明の技術的背景および先行技術）
従来より放射線画像のデジタルサブトラクシヨ
ンが公知となつている。この放射線画像のデジタ
ルサブトラクシヨンとは、異なつた条件で撮影し
た２つの放射線画像を光電的に読み出してデジタ
ル画像信号を得た後、これらのデジタル画像信号
を両画像の各画素を対応させて減算処理し、放射
線画像中の特定の構造物の画像を形成するための
差信号を得る方法であり、このようにして得た差
信号を用いて特定構造物のみが抽出された放射線
画像を再生することができる。
ンが公知となつている。この放射線画像のデジタ
ルサブトラクシヨンとは、異なつた条件で撮影し
た２つの放射線画像を光電的に読み出してデジタ
ル画像信号を得た後、これらのデジタル画像信号
を両画像の各画素を対応させて減算処理し、放射
線画像中の特定の構造物の画像を形成するための
差信号を得る方法であり、このようにして得た差
信号を用いて特定構造物のみが抽出された放射線
画像を再生することができる。
このサブトラクシヨン処理には、基本的に次の
２つの方法がある。即ち、造影剤注入による特定
の構造物が強調された放射線画像の画像信号か
ら、造影剤が注入されていない放射線画像の画像
信号を引き算（サブトラクト）することによつて
特定の構造物を抽出するいわゆる時間サブトラク
シヨン処理と、同一の被写体に対して相異なるエ
ネルギー分布を有する放射線を照射し、それによ
り特定の構造物が特有の放射線エネルギー吸収特
性を有することを利用して特定構造物が異なる画
像を２つの放射線画像間に存在せしめ、その後こ
の２つの放射線画像の画像信号間で適当な重みづ
けをした上で引き算（サブトラクト）を行ない特
定の構造物の画像を抽出するいわゆるエネルギー
サブトラクシヨン処理である。
２つの方法がある。即ち、造影剤注入による特定
の構造物が強調された放射線画像の画像信号か
ら、造影剤が注入されていない放射線画像の画像
信号を引き算（サブトラクト）することによつて
特定の構造物を抽出するいわゆる時間サブトラク
シヨン処理と、同一の被写体に対して相異なるエ
ネルギー分布を有する放射線を照射し、それによ
り特定の構造物が特有の放射線エネルギー吸収特
性を有することを利用して特定構造物が異なる画
像を２つの放射線画像間に存在せしめ、その後こ
の２つの放射線画像の画像信号間で適当な重みづ
けをした上で引き算（サブトラクト）を行ない特
定の構造物の画像を抽出するいわゆるエネルギー
サブトラクシヨン処理である。
このサブトラクシヨン処理は特に医療用のＸ線
写真の画像処理において診断上きわめて有効な方
法であるため、近年大いに注目され、電子工学技
術を駆使してその研究、開発が盛んに進められて
いる。この技術は、特にデジタルサブトラクシヨ
ン処理（通常Digital Radiography）と呼ばれ、
DRと略称されている。
写真の画像処理において診断上きわめて有効な方
法であるため、近年大いに注目され、電子工学技
術を駆使してその研究、開発が盛んに進められて
いる。この技術は、特にデジタルサブトラクシヨ
ン処理（通常Digital Radiography）と呼ばれ、
DRと略称されている。
さらに最近では例えば特開昭58−163340号公報
に示されるように、きわめて広い放射線露出域を
有する蓄積性螢光体シートを使用し、これらの螢
光体シートに前述のように異なつた条件で同一の
被写体を透過した放射線を照射して、これらの螢
光体シートに特定構造物の画像情報が異なる放射
線画像を蓄積記録し、これらの蓄積画像を励起光
による走査により読み出してデジタル信号に変換
し、これらデジタル信号により前記デジタルサブ
トラクシヨンを行なうことも提案されている。上
記蓄積性螢光体シートとは、例えば特開昭55−
12429号公報に開示されているように放射線（Ｘ
線、α線、β線、γ線、紫外線等）を照射すると
その放射線エネルギーの一部を螢光体中に蓄積
し、その後可視光等の励起光を照射すると蓄積さ
れた放射線エネルギー量に応じて螢光体が輝尽発
光を示すもので、きわめて広いラチチユード（露
出域）を有し、かつ著しく高い解像力を有するも
のである。したがつて、この螢光体シートに蓄積
記録された放射線画像情報を利用して前記デジタ
ルサブトラクシヨンを行なえば、放射線量が変動
しても常に十分な画像情報を得ることができ、診
断能の高い放射線画像を得ることができる。
に示されるように、きわめて広い放射線露出域を
有する蓄積性螢光体シートを使用し、これらの螢
光体シートに前述のように異なつた条件で同一の
被写体を透過した放射線を照射して、これらの螢
光体シートに特定構造物の画像情報が異なる放射
線画像を蓄積記録し、これらの蓄積画像を励起光
による走査により読み出してデジタル信号に変換
し、これらデジタル信号により前記デジタルサブ
トラクシヨンを行なうことも提案されている。上
記蓄積性螢光体シートとは、例えば特開昭55−
12429号公報に開示されているように放射線（Ｘ
線、α線、β線、γ線、紫外線等）を照射すると
その放射線エネルギーの一部を螢光体中に蓄積
し、その後可視光等の励起光を照射すると蓄積さ
れた放射線エネルギー量に応じて螢光体が輝尽発
光を示すもので、きわめて広いラチチユード（露
出域）を有し、かつ著しく高い解像力を有するも
のである。したがつて、この螢光体シートに蓄積
記録された放射線画像情報を利用して前記デジタ
ルサブトラクシヨンを行なえば、放射線量が変動
しても常に十分な画像情報を得ることができ、診
断能の高い放射線画像を得ることができる。
しかしながら、前述のような蓄積性螢光体シー
トを用いてエネルギーサブトラクシヨンを行なう
場合、サブトラクシヨンによつて得られた差信号
をそのまま用いて出力画像を形成すると、各出力
画像によつて濃度範囲（あるいは輝度範囲）やコ
ントラストがまちまちになることが多い。これは
前述した蓄積性螢光体シートの感度のバラツキ
や、画像読取器の感度バラツキ、撮影放射線エネ
ルギー量のバラツキ等に起因して、サブトラクシ
ヨンに供する２つの画像の画像信号が変動するこ
とによるものである。サブトラクシヨンによつて
得られる差信号は微小なものであるので、差信号
を増幅するような階調処理を行なうことが有効で
あるが、その時上記のように画像信号が変動する
と、差信号におけるバラツキ成分は本来の信号成
分に比べてかなり大きなものとなり、濃度範囲や
コントラストの大きな変動となつて表われる。上
記のようにサブトラクシヨン画像の濃度（あるい
は輝度）やコントラストがまちまちであると、複
数のサブトラクシヨン画像を比較して診断を下す
場合等において適正な診断が妨げられたり、また
人間の視覚特性の良好な領域から画像の濃度がは
ずれてしまうことが指摘されている。
トを用いてエネルギーサブトラクシヨンを行なう
場合、サブトラクシヨンによつて得られた差信号
をそのまま用いて出力画像を形成すると、各出力
画像によつて濃度範囲（あるいは輝度範囲）やコ
ントラストがまちまちになることが多い。これは
前述した蓄積性螢光体シートの感度のバラツキ
や、画像読取器の感度バラツキ、撮影放射線エネ
ルギー量のバラツキ等に起因して、サブトラクシ
ヨンに供する２つの画像の画像信号が変動するこ
とによるものである。サブトラクシヨンによつて
得られる差信号は微小なものであるので、差信号
を増幅するような階調処理を行なうことが有効で
あるが、その時上記のように画像信号が変動する
と、差信号におけるバラツキ成分は本来の信号成
分に比べてかなり大きなものとなり、濃度範囲や
コントラストの大きな変動となつて表われる。上
記のようにサブトラクシヨン画像の濃度（あるい
は輝度）やコントラストがまちまちであると、複
数のサブトラクシヨン画像を比較して診断を下す
場合等において適正な診断が妨げられたり、また
人間の視覚特性の良好な領域から画像の濃度がは
ずれてしまうことが指摘されている。
なお、蓄積性螢光体シートを用いるエネルギー
サブトラクシヨンにおいて、サブトラクシヨンに
供される放射線画像（すなわち特定構造物の画像
情報が互いに異なる放射線画像）を得るには、放
射線源の管電圧を変えて放射線エネルギーを変化
させる他に、放射線源と被写体との間に放射線の
エネルギー分布を変化させるフイルタを出し入れ
する等、その他の公知の方法が用いられてもよい
し、さらには本出願人による特願昭57−193765号
に開示されているように、蓄積性螢光体シートと
フイルタとの積層体等を使用して、１回の放射線
撮影によつて得るようにしてもよい。
サブトラクシヨンにおいて、サブトラクシヨンに
供される放射線画像（すなわち特定構造物の画像
情報が互いに異なる放射線画像）を得るには、放
射線源の管電圧を変えて放射線エネルギーを変化
させる他に、放射線源と被写体との間に放射線の
エネルギー分布を変化させるフイルタを出し入れ
する等、その他の公知の方法が用いられてもよい
し、さらには本出願人による特願昭57−193765号
に開示されているように、蓄積性螢光体シートと
フイルタとの積層体等を使用して、１回の放射線
撮影によつて得るようにしてもよい。
（発明の目的）
本発明は上記のような事情に鑑みてなされたも
のであり、濃度範囲、コントラストが一定で診断
能の良いエネルギーサブトラクシヨン画像を得る
ことができる、サブトラクシヨン画像の階調処理
最適化方法および装置を提供することを目的とす
るものである。
のであり、濃度範囲、コントラストが一定で診断
能の良いエネルギーサブトラクシヨン画像を得る
ことができる、サブトラクシヨン画像の階調処理
最適化方法および装置を提供することを目的とす
るものである。
（発明の構成）
本発明のサブトラクシヨン画像の階調処理最適
化方法は、前述したように蓄積性螢光体シートを
用いて行なうエネルギーサブトラクシヨン処理に
おいて、サブトラクシヨンによつて得た差信号の
ヒストグラムを求めてそのヒストグラムから互い
に異なる少なくとも２つの信号値を求め、次に上
記差信号を、上記信号値を予め決められた少なく
とも２つの所定出力画像濃度に対応させて階調処
理するようにしたものである。
化方法は、前述したように蓄積性螢光体シートを
用いて行なうエネルギーサブトラクシヨン処理に
おいて、サブトラクシヨンによつて得た差信号の
ヒストグラムを求めてそのヒストグラムから互い
に異なる少なくとも２つの信号値を求め、次に上
記差信号を、上記信号値を予め決められた少なく
とも２つの所定出力画像濃度に対応させて階調処
理するようにしたものである。
上記ヒストグラムの少なくとも２つの信号値と
して、信号最大値、最小値を用い、これらを予め
決められた出力画像濃度最大値、最小値に対応さ
せれば、最も広い範囲で濃度範囲、コントラスト
を一定に設定することができる。
して、信号最大値、最小値を用い、これらを予め
決められた出力画像濃度最大値、最小値に対応さ
せれば、最も広い範囲で濃度範囲、コントラスト
を一定に設定することができる。
上述のように差信号ヒストグラムのある信号値
を所定出力画像濃度に対応させるためには、差信
号を階調処理前に一律に信号変換してもよいし、
また階調変換テーブルを補正してもよい。
を所定出力画像濃度に対応させるためには、差信
号を階調処理前に一律に信号変換してもよいし、
また階調変換テーブルを補正してもよい。
上記方法を実施する本発明のサブトラクシヨン
画像の階調処理最適化装置は、差信号を一律に信
号変換する場合には、前述したような励起光走査
と輝尽発光光の光電的読み出しにより、蓄積性螢
光体シートの放射線画像のデジタル信号を得る画
像読取手段と、この画像読取手段が読み取つた、
互いに特定構造物の画像情報が異なる被写体のデ
ジタル画像信号とを、対応する画素間で減算して
前記特定構造物の画像を形成する差信号を得るサ
ブトラクシヨン演算手段と、この差信号に階調変
換テーブルに基づいて階調処理を施す画像処理手
段と、前記差信号のヒストグラムを求めて該ヒス
トグラムの互いに異なる少なくとも２つの信号値
を求めるヒストグラム演算手段と、前記少なくと
も２つの信号値が予め決められた少なくとも２つ
の所定出力画像濃度に対応するように、前記階調
処理を受ける前の前記差信号を一律に変換する信
号変換回路とから構成される。
画像の階調処理最適化装置は、差信号を一律に信
号変換する場合には、前述したような励起光走査
と輝尽発光光の光電的読み出しにより、蓄積性螢
光体シートの放射線画像のデジタル信号を得る画
像読取手段と、この画像読取手段が読み取つた、
互いに特定構造物の画像情報が異なる被写体のデ
ジタル画像信号とを、対応する画素間で減算して
前記特定構造物の画像を形成する差信号を得るサ
ブトラクシヨン演算手段と、この差信号に階調変
換テーブルに基づいて階調処理を施す画像処理手
段と、前記差信号のヒストグラムを求めて該ヒス
トグラムの互いに異なる少なくとも２つの信号値
を求めるヒストグラム演算手段と、前記少なくと
も２つの信号値が予め決められた少なくとも２つ
の所定出力画像濃度に対応するように、前記階調
処理を受ける前の前記差信号を一律に変換する信
号変換回路とから構成される。
また上記信号変換回路に代えて、前記少なくと
も２つの信号値が予め決められた少なくとも２つ
の所定出力画像濃度に対応するように、前記階調
変換テーブルを補正する階調変換テーブル補正回
路を用いれば、差信号はそのまま用いて階調処理
を最適化できる。
も２つの信号値が予め決められた少なくとも２つ
の所定出力画像濃度に対応するように、前記階調
変換テーブルを補正する階調変換テーブル補正回
路を用いれば、差信号はそのまま用いて階調処理
を最適化できる。
（実施態様）
以下、図面に示す実施態様に基づいて本発明を
詳細に説明する。
詳細に説明する。
第１図は２枚の蓄積性螢光体シートＡ，Ｂに同
一の被写体１を透過したＸ線２を、それぞれエネ
ルギーを変えて照射する状態を示す。すなわち第
１の蓄積性螢光体シートＡに被写体１のＸ線透過
像を蓄積記録し、次いで短時間内で蓄積性螢光体
シートＡ，Ｂを素早く取り替えると同時に、Ｘ線
源３の管電圧を変えて、透過Ｘ線のエネルギーが
異なる被写体１のＸ線画像を蓄積性螢光体シート
Ｂに蓄積記録する。このとき蓄積性螢光体シート
ＡとＢとで被写体１の位置関係は同じとする。
一の被写体１を透過したＸ線２を、それぞれエネ
ルギーを変えて照射する状態を示す。すなわち第
１の蓄積性螢光体シートＡに被写体１のＸ線透過
像を蓄積記録し、次いで短時間内で蓄積性螢光体
シートＡ，Ｂを素早く取り替えると同時に、Ｘ線
源３の管電圧を変えて、透過Ｘ線のエネルギーが
異なる被写体１のＸ線画像を蓄積性螢光体シート
Ｂに蓄積記録する。このとき蓄積性螢光体シート
ＡとＢとで被写体１の位置関係は同じとする。
このようにして、少なくとも一部の画像情報が
異なる２つの放射線画像を２枚の蓄積性螢光体シ
ートＡ，Ｂに蓄積記録する。次にこれら２枚の蓄
積性螢光体シートＡ，Ｂから、第２図に示すよう
な画像読取手段によつてＸ線画像を読み取り、画
像を表わすデジタル画像信号を得る。先ず、蓄積
性螢光体シートＡを矢印Ｙの方向に副走査のため
に移動させながら、レーザー光源１０からのレー
ザー光１１を走査ミラー１２によつてＸ方向に主
走査させ、螢光体シートＡから蓄積Ｘ線エネルギ
ーを、蓄積記録されたＸ線画像にしたがつて輝尽
発光光１３として発散させる。輝尽発光光１３は
透明なアクリル板を成形して作られた集光板１４
の一端面からこの集光板１４の内部に入射し、中
を全反射を繰返しつつフオトマル１５に至り、輝
尽発光光１３の発光量が画像信号Ｓとして出力さ
れる。この出力された画像信号Ｓは増幅器とＡ／
Ｄ変換器を含む対数変換器１６により対数値
（logS）のデジタル画像信号logSAに変換される。
このデジタル画像信号logSAは例えば磁気デイス
ク等の記憶媒体１７に記憶される。次に、全く同
様にして、もう１枚の蓄積性螢光体シートＢの記
録画像が読み出され、そのデジタル画像信号
logSBが同様に記憶媒体１７に記憶される。
異なる２つの放射線画像を２枚の蓄積性螢光体シ
ートＡ，Ｂに蓄積記録する。次にこれら２枚の蓄
積性螢光体シートＡ，Ｂから、第２図に示すよう
な画像読取手段によつてＸ線画像を読み取り、画
像を表わすデジタル画像信号を得る。先ず、蓄積
性螢光体シートＡを矢印Ｙの方向に副走査のため
に移動させながら、レーザー光源１０からのレー
ザー光１１を走査ミラー１２によつてＸ方向に主
走査させ、螢光体シートＡから蓄積Ｘ線エネルギ
ーを、蓄積記録されたＸ線画像にしたがつて輝尽
発光光１３として発散させる。輝尽発光光１３は
透明なアクリル板を成形して作られた集光板１４
の一端面からこの集光板１４の内部に入射し、中
を全反射を繰返しつつフオトマル１５に至り、輝
尽発光光１３の発光量が画像信号Ｓとして出力さ
れる。この出力された画像信号Ｓは増幅器とＡ／
Ｄ変換器を含む対数変換器１６により対数値
（logS）のデジタル画像信号logSAに変換される。
このデジタル画像信号logSAは例えば磁気デイス
ク等の記憶媒体１７に記憶される。次に、全く同
様にして、もう１枚の蓄積性螢光体シートＢの記
録画像が読み出され、そのデジタル画像信号
logSBが同様に記憶媒体１７に記憶される。
次に、上述のようにして得られたデジタル画像
信号logSA，logSBを用いてサブトラクシヨン処理
を行なう。第３図は本発明方法の第１実施態様に
よる階調処理最適化方法を適用して行なわれるサ
ブトラクシヨン処理の流れを示している。まず前
記記憶媒体１７内の画像フアイル１７Ａと、画像
フアイル１７Ｂからそれぞれ、前記デジタル画像
信号logSA，logSBが読み出され、サブトラクシヨ
ン演算回路１８に入力される。このサブトラクシ
ヨン演算回路１８は、上記２つのデジタル画像信
号logSAとlogSBを適当な重みづけをした上で対応
する画素毎に減算し、デジタルの差信号 Ssub＝ａ・logSA−ｂ・logSB＋ｃ（ａ，ｂは重み
づけ係数、ｃは概略一定濃度にするようなバイア
ス成分である） を求める。この差信号Ssubは一たん画像フアイ
ル１９に記憶されてから、画像処理回路２０に入
力され、該画像処理回路２０において後述の信号
変換回路２０ａを介して階調変換テーブル２０ｂ
に基づいて階調処理される。
信号logSA，logSBを用いてサブトラクシヨン処理
を行なう。第３図は本発明方法の第１実施態様に
よる階調処理最適化方法を適用して行なわれるサ
ブトラクシヨン処理の流れを示している。まず前
記記憶媒体１７内の画像フアイル１７Ａと、画像
フアイル１７Ｂからそれぞれ、前記デジタル画像
信号logSA，logSBが読み出され、サブトラクシヨ
ン演算回路１８に入力される。このサブトラクシ
ヨン演算回路１８は、上記２つのデジタル画像信
号logSAとlogSBを適当な重みづけをした上で対応
する画素毎に減算し、デジタルの差信号 Ssub＝ａ・logSA−ｂ・logSB＋ｃ（ａ，ｂは重み
づけ係数、ｃは概略一定濃度にするようなバイア
ス成分である） を求める。この差信号Ssubは一たん画像フアイ
ル１９に記憶されてから、画像処理回路２０に入
力され、該画像処理回路２０において後述の信号
変換回路２０ａを介して階調変換テーブル２０ｂ
に基づいて階調処理される。
階調処理を受けた差信号Ssub′は、例えばCRT
等のデイスプレイ装置や、走査記録装置等の再生
記録装置２１に入力され、該差信号Ssub′によつ
てサブトラクシヨン画像が再生記録される。第４
図はサブトラクシヨン画像再生記録システムの一
例として、画像走査記録装置を示すものである。
感光フイルム３０を矢印Ｙの副走査方向へ移動さ
せるとともにレーザービーム３１をこの感光フイ
ルム３０上にＸ方向に主走査させ、レーザービー
ム３１をＡ／Ｏ変調器３２により画像信号供給器
３３からの画像信号によつて変調することによ
り、感光フイルム３０上に可視像を形成する。こ
の変調用画像信号として、前記差信号Ssub′を使
用すれば、デジタルサブトラクシヨン処理による
所望の特定構造物のみの画像を感光フイルム３０
上に再生記録することができる。
等のデイスプレイ装置や、走査記録装置等の再生
記録装置２１に入力され、該差信号Ssub′によつ
てサブトラクシヨン画像が再生記録される。第４
図はサブトラクシヨン画像再生記録システムの一
例として、画像走査記録装置を示すものである。
感光フイルム３０を矢印Ｙの副走査方向へ移動さ
せるとともにレーザービーム３１をこの感光フイ
ルム３０上にＸ方向に主走査させ、レーザービー
ム３１をＡ／Ｏ変調器３２により画像信号供給器
３３からの画像信号によつて変調することによ
り、感光フイルム３０上に可視像を形成する。こ
の変調用画像信号として、前記差信号Ssub′を使
用すれば、デジタルサブトラクシヨン処理による
所望の特定構造物のみの画像を感光フイルム３０
上に再生記録することができる。
このようにして感光フイルム３０上に記録再生
されるエネルギーサブトラクシヨン画像は、２枚
の蓄積性螢光体シートＡ，Ｂの感度のバラツキや
撮影Ｘ線エネルギーのバラツキなどにより前記差
信号Ssubが変動するために、各画像により濃度
範囲、コントラストがまちまちになつてしまうこ
とが多い。そこで第３図に示されるように上記差
信号Ssubをヒストグラム演算回路２２に入力し、
該差信号Ssubのヒストグラムを求める。そして
該ヒストグラム演算回路２２は第５図に示すよう
に、このヒストグラムから一例として上記差信号
Ssubの最大値Smax、最小値Sminを求めてその
情報を前記信号変換回路２０ａに入力する。この
信号変換回路２０ａはこれら最大値Smax、最小
値Sminがそれぞれ常に一定の値Cmax、Cminを
とるように差信号Ssubを変換する。すなわち例
えば Ssub→ｐ−Ssub＋ｑ ｐ＝(Cmax-Cmin)／(Smax-Smin) ｑ＝（Cmax＋Cmin）／２ −（Smax＋Smin）／２ なる変換を行なう。
されるエネルギーサブトラクシヨン画像は、２枚
の蓄積性螢光体シートＡ，Ｂの感度のバラツキや
撮影Ｘ線エネルギーのバラツキなどにより前記差
信号Ssubが変動するために、各画像により濃度
範囲、コントラストがまちまちになつてしまうこ
とが多い。そこで第３図に示されるように上記差
信号Ssubをヒストグラム演算回路２２に入力し、
該差信号Ssubのヒストグラムを求める。そして
該ヒストグラム演算回路２２は第５図に示すよう
に、このヒストグラムから一例として上記差信号
Ssubの最大値Smax、最小値Sminを求めてその
情報を前記信号変換回路２０ａに入力する。この
信号変換回路２０ａはこれら最大値Smax、最小
値Sminがそれぞれ常に一定の値Cmax、Cminを
とるように差信号Ssubを変換する。すなわち例
えば Ssub→ｐ−Ssub＋ｑ ｐ＝(Cmax-Cmin)／(Smax-Smin) ｑ＝（Cmax＋Cmin）／２ −（Smax＋Smin）／２ なる変換を行なう。
変換された信号ｐ・Ssub＋ｑは前述した階調
変換テーブル２０ｂに基づいて階調処理される
が、ここで上記一定値Cmax、Cminは第５図に
示すように、サブトラクシヨン画像の濃度最大値
Dmax、濃度最小値Dminを得る信号となつてい
る。したがつて上記変換された信号ｐ・Ssub＋
ｑによるサブトラクシヨン画像は常に一定の濃度
範囲、コントラストを有するようになる。
変換テーブル２０ｂに基づいて階調処理される
が、ここで上記一定値Cmax、Cminは第５図に
示すように、サブトラクシヨン画像の濃度最大値
Dmax、濃度最小値Dminを得る信号となつてい
る。したがつて上記変換された信号ｐ・Ssub＋
ｑによるサブトラクシヨン画像は常に一定の濃度
範囲、コントラストを有するようになる。
なお上述のように、信号変換回路２０ａによつ
て差信号Ssubを変換する代わりに、第６図に示
すような階調変換テーブル補正回路１２０により
階調変換テーブル２０ｂを第７図に示すように修
正することによつても同様の効果が得られる。ま
た差信号のヒストグラムは、ある特定の画像部分
（例えば関心領域）のみについて求めてもよい。
て差信号Ssubを変換する代わりに、第６図に示
すような階調変換テーブル補正回路１２０により
階調変換テーブル２０ｂを第７図に示すように修
正することによつても同様の効果が得られる。ま
た差信号のヒストグラムは、ある特定の画像部分
（例えば関心領域）のみについて求めてもよい。
以上２枚のシートＡ，Ｂを使用する実施態様に
ついて説明したが、その他３枚のシートにそれぞ
れ異なるエネルギーで放射線撮影し、これらシー
トから得られる画像信号を減算処理して（例えば
ａ×logSA＋ｂ×logSB−ｃ×logSC＋ｄなど。但
し、ここでａ，ｂおよびｃは重みづけ係数、ｄは
概略一定濃度にするようなバイアス部分である。）
サブトラクシヨン画像を得ることもできるが、本
発明はこのような場合にも適用可能である。
ついて説明したが、その他３枚のシートにそれぞ
れ異なるエネルギーで放射線撮影し、これらシー
トから得られる画像信号を減算処理して（例えば
ａ×logSA＋ｂ×logSB−ｃ×logSC＋ｄなど。但
し、ここでａ，ｂおよびｃは重みづけ係数、ｄは
概略一定濃度にするようなバイアス部分である。）
サブトラクシヨン画像を得ることもできるが、本
発明はこのような場合にも適用可能である。
先に述べたように、サブトラクシヨンに供され
る放射線画像（すなわち特定構造物の画像情報が
互いに異なる放射線画像）を得るには、前記実施
態様におけるように放射線源の管電圧を変えて放
射線エネルギーを変化させる他に、放射線源と被
写体との間に放射線のエネルギー分布を変化させ
るフイルタを出し入れする等、その他の公知の方
法が用いられてもよいし、さらには本出願人によ
る特願昭57−193765号に開示されているように蓄
積性螢光体シートとフイルタとの積層体等を使用
して、１回の放射線撮影によつて得るようにして
もよい。
る放射線画像（すなわち特定構造物の画像情報が
互いに異なる放射線画像）を得るには、前記実施
態様におけるように放射線源の管電圧を変えて放
射線エネルギーを変化させる他に、放射線源と被
写体との間に放射線のエネルギー分布を変化させ
るフイルタを出し入れする等、その他の公知の方
法が用いられてもよいし、さらには本出願人によ
る特願昭57−193765号に開示されているように蓄
積性螢光体シートとフイルタとの積層体等を使用
して、１回の放射線撮影によつて得るようにして
もよい。
（発明の効果）
以上詳細に説明した通り本発明によれば、エネ
ルギーサブトラクシヨン画像の全体的な濃度、コ
ントラストを各画像間で共通に保つことができる
から、極めて診断能に優れたサブトラクシヨン画
像を得ることが可能となる。そして濃度、コント
ラストを所定の範囲に保てるから、サブトラクシ
ヨン画像を出力する記憶媒体のダイナミツクレン
ジを最も有効に活用することが可能になる。また
本発明方法は電気回路により自動的に実施可能で
あるから、特にサブトラクシヨン演算処理を連続
的に行ない、サブトラクシヨン画像をCRTを用
いてリアルタイムで観察する場合に特に有効であ
る。
ルギーサブトラクシヨン画像の全体的な濃度、コ
ントラストを各画像間で共通に保つことができる
から、極めて診断能に優れたサブトラクシヨン画
像を得ることが可能となる。そして濃度、コント
ラストを所定の範囲に保てるから、サブトラクシ
ヨン画像を出力する記憶媒体のダイナミツクレン
ジを最も有効に活用することが可能になる。また
本発明方法は電気回路により自動的に実施可能で
あるから、特にサブトラクシヨン演算処理を連続
的に行ない、サブトラクシヨン画像をCRTを用
いてリアルタイムで観察する場合に特に有効であ
る。
第１図は本発明方法における放射線画像の蓄積
記録ステツプを示す説明図、第２図は上記蓄積記
録がなされた蓄積性螢光体シートからの放射線画
像情報読取りを説明する概略図、第３図は本発明
方法の第１実施態様方法を適用して階調処理を行
なうエネルギーサブトラクシヨン処理の概要を説
明する概略図、第４図はサブトラクシヨン画像の
再生記録システムの一例を示す概略図、第５図は
本発明の第１実施態様による階調処理最適化方法
を説明する説明図、第６図は本発明方法の第２実
施態様方法を適用して階調処理を行なうエネルギ
ーサブトラクシヨン処理の概要を説明する概略
図、第７図は本発明の第２実施態様による階調処
理最適化方法を説明する説明図である。 １……被写体、２……Ｘ線、３……Ｘ線源、１
０……レーザー光源、１１……レーザー光、１２
……走査ミラー、１３……輝尽発光光、１５……
フオトマル、１６……対数変換器、１８……サブ
トラクシヨン演算回路、２０……画像処理回路、
２０ａ……信号変換回路、２０ｂ……階調変換テ
ーブル、２２……ヒストグラム演算回路、１２０
……階調変換テーブル補正回路、Ａ，Ｂ……蓄積
性螢光体シート、logSA，logSB……デジタル画像
信号、Ssub……デジタル画像信号の差信号、
Ssub′……階調処理された差信号。
記録ステツプを示す説明図、第２図は上記蓄積記
録がなされた蓄積性螢光体シートからの放射線画
像情報読取りを説明する概略図、第３図は本発明
方法の第１実施態様方法を適用して階調処理を行
なうエネルギーサブトラクシヨン処理の概要を説
明する概略図、第４図はサブトラクシヨン画像の
再生記録システムの一例を示す概略図、第５図は
本発明の第１実施態様による階調処理最適化方法
を説明する説明図、第６図は本発明方法の第２実
施態様方法を適用して階調処理を行なうエネルギ
ーサブトラクシヨン処理の概要を説明する概略
図、第７図は本発明の第２実施態様による階調処
理最適化方法を説明する説明図である。 １……被写体、２……Ｘ線、３……Ｘ線源、１
０……レーザー光源、１１……レーザー光、１２
……走査ミラー、１３……輝尽発光光、１５……
フオトマル、１６……対数変換器、１８……サブ
トラクシヨン演算回路、２０……画像処理回路、
２０ａ……信号変換回路、２０ｂ……階調変換テ
ーブル、２２……ヒストグラム演算回路、１２０
……階調変換テーブル補正回路、Ａ，Ｂ……蓄積
性螢光体シート、logSA，logSB……デジタル画像
信号、Ssub……デジタル画像信号の差信号、
Ssub′……階調処理された差信号。
Claims (1)
- 【特許請求の範囲】 １ ２枚以上の蓄積性螢光体シートに、放射線エ
ネルギー吸収特性が他とは異なる特定構造物を含
む被写体を透過した放射線を照射して、これらの
螢光体シートに前記被写体の特定構造物の画像情
報が互いに異なる放射線画像を蓄積記録し、これ
らの螢光体シートに励起光を走査して前記放射線
画像を輝尽発光光に変換し、この輝尽発光光の発
光量を光電的に読み出してデジタル画像信号に変
換し、各画像の対応する画素間でこのデジタル画
像信号の減算を行なつて放射線画像の特定構造物
の画像を形成する差信号を得、その後この差信号
に階調処理を施すエネルギーサブトラクシヨンの
階調処理において、前記差信号のヒストグラムを
求めて該ヒストグラムから互いに異なる少なくと
も２つの信号値を求め、次に前記差信号を、前記
少なくとも２つの信号値を予め決められた少なく
とも２つの所定出力画像濃度のそれぞれに対応さ
せて階調処理することを特徴とするサブトラクシ
ヨン画像の階調処理最適化方法。 ２ 前記少なくとも２つの信号値が信号最大値、
信号最小値であり、これらの信号最大値、信号最
小値をそれぞれ予め決められた出力画像濃度最大
値、濃度最小値に対応させることを特徴とする特
許請求の範囲第１項記載のサブトラクシヨン画像
の階調処理最適化方法。 ３ 前記差信号を階調処理前に一律に信号変換す
ることにより、該差信号の前記少なくとも２つの
信号値を前記少なくとも２つの所定出力画像濃度
に対応させることを特徴とする特許請求の範囲第
１項または第２項記載のサブトラクシヨン画像の
階調処理最適化方法。 ４ 前記階調処理に用いられる階調変換テーブル
を修正することにより、前記差信号の前記少なく
とも２つの信号値を前記少なくとも２つの所定出
力画像濃度に対応させることを特徴とする特許請
求の範囲第１項または第２項記載のサブトラクシ
ヨン画像の階調処理最適化方法。 ５ 放射線画像が蓄積記録された蓄積性蛍光体シ
ートに励起光を走査し、それによつて前記蓄積性
螢光体シートから発せられた輝尽発光光を光電的
に読み出してデジタル画像信号に変換する画像読
取手段と、この画像読取手段が読み取つた、互い
に特定構造物の画像情報が異なる被写体のデジタ
ル画像信号とを、対応する画素間で減算して前記
特定構造物の画像を形成する差信号を得るサブト
ラクシヨン演算手段と、この差信号に階調変換テ
ーブルに基づいて階調処理を施す画像処理手段
と、前記差信号のヒストグラムを求めて該ヒスト
グラムの互いに異なる少なくとも２つの信号値を
求めるヒストグラム演算手段と、前記少なくとも
２つの信号値が予め決められた少なくとも２つの
所定出力画像濃度に対応するように、前記階調処
理を受ける前の前記差信号を一律に変換する信号
変換回路とからなるサブトラクシヨン画像の階調
処理最適化装置。 ６ 放射線画像が蓄積記録された蓄積性蛍光体シ
ートに励起光を走査し、それによつて前記蓄積性
螢光体シートから発せられた輝尽発光光を光電的
に読み出してデジタル画像信号に変換する画像読
取手段と、この画像読取手段が読み取つた、互い
に特定構造物の画像情報が異なる被写体のデジタ
ル画像信号とを、対応する画素間で減算して前記
特定構造物の画像を形成する差信号を得るサブト
ラクシヨン演算手段と、この差信号に階調変換テ
ーブルに基づいて階調処理を施す画像処理手段
と、前記差信号のヒストグラムを求めて該ヒスト
グラムの互いに異なる少なくとも２つの信号値を
求めるヒストグラム演算手段と、前記少なくとも
２つの信号値が予め決められた少なくとも２つの
所定出力画像濃度に対応するように、前記階調変
換テーブルを補正する階調変換テーブル補正回路
とからなるサブトラクシヨン画像の階調処理最適
化装置。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59107911A JPS60250792A (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 |
|---|---|---|---|
| JP59107911A JPS60250792A (ja) | 1984-05-28 | 1984-05-28 | サブトラクシヨン画像の階調処理最適化方法および装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS60250792A JPS60250792A (ja) | 1985-12-11 |
| JPH0363878B2 true JPH0363878B2 (ja) | 1991-10-02 |
Family
ID=14471181
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP59107911A Granted JPS60250792A (ja) | 1984-04-23 | 1984-05-28 | サブトラクシヨン画像の階調処理最適化方法および装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS60250792A (ja) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0661325B2 (ja) * | 1986-07-24 | 1994-08-17 | 富士写真フイルム株式会社 | 放射線画像情報処理方法 |
| JPH01262453A (ja) * | 1988-04-13 | 1989-10-19 | Rigaku Corp | 蓄積性蛍光板を利用したｘ線回折装置における回折ｘ線強度の読取り方法 |
| WO2010050031A1 (ja) * | 2008-10-30 | 2010-05-06 | 株式会社島津製作所 | 放射線撮影装置 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS55116340A (en) * | 1979-02-28 | 1980-09-06 | Fuji Photo Film Co Ltd | Method and device for processing gradation of radiation picture |
| JPS5867239A (ja) * | 1981-10-20 | 1983-04-21 | 株式会社東芝 | 画像表示装置 |
| JPS58163340A (ja) * | 1982-03-20 | 1983-09-28 | 富士写真フイルム株式会社 | 放射線画像のサブトラクシヨン処理方法 |
-
1984
- 1984-05-28 JP JP59107911A patent/JPS60250792A/ja active Granted
Also Published As
| Publication number | Publication date |
|---|---|
| JPS60250792A (ja) | 1985-12-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA1223980A (en) | Method of adjusting radiation image read-out apparatus | |
| US4747052A (en) | Radiation image processing | |
| US5048110A (en) | Method and apparatus for automatically correcting subtraction image density | |
| US5535289A (en) | Method for reducing noise in energy subtraction images | |
| US5049746A (en) | Method and apparatus for displaying energy subtraction images | |
| JPH0362411B2 (ja) | ||
| US5485371A (en) | Method for forming energy subtraction radiation images, and method and apparatus for smoothing radiation images | |
| US4761739A (en) | Density correcting method and apparatus for energy substraction image | |
| JP3901770B2 (ja) | エネルギーサブトラクション処理方法および装置 | |
| US5028781A (en) | Density correcting method and apparatus for subtraction image | |
| JPH0363878B2 (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 | |
| JPS60251475A (ja) | サブトラクシヨン画像の階調処理最適化方法および装置 | |
| JPH0262074B2 (ja) | ||
| JPH0262075B2 (ja) | ||
| JPH0363877B2 (ja) | ||
| US5857037A (en) | Image processing method and apparatus | |
| US4806759A (en) | Method of adjusting radiation image read-out conditions | |
| JPH067669B2 (ja) | エネルギ−サブトラクシヨン画像の画質改善方法および装置 | |
| JPH0471395B2 (ja) | ||
| JP2824878B2 (ja) | 放射線画像のエネルギーサブトラクション撮影処理方法 | |
| JPH0262076B2 (ja) | ||
| JPS6175966A (ja) | 放射線画像情報処理方法および装置 |