JPH05232608A - 放射線画像の重ね合せ処理方法および装置 - Google Patents
放射線画像の重ね合せ処理方法および装置Info
- Publication number
- JPH05232608A JPH05232608A JP4045529A JP4552992A JPH05232608A JP H05232608 A JPH05232608 A JP H05232608A JP 4045529 A JP4045529 A JP 4045529A JP 4552992 A JP4552992 A JP 4552992A JP H05232608 A JPH05232608 A JP H05232608A
- Authority
- JP
- Japan
- Prior art keywords
- image
- signal
- addition
- image signal
- radiation
- 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.)
- Granted
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 89
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000035945 sensitivity Effects 0.000 claims abstract description 68
- 238000003672 processing method Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 4
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 34
- 230000005284 excitation Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- NIXOWILDQLNWCW-UHFFFAOYSA-N acrylic acid group Chemical group C(C=C)(=O)O NIXOWILDQLNWCW-UHFFFAOYSA-N 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 238000001454 recorded image Methods 0.000 description 1
- GGCZERPQGJTIQP-UHFFFAOYSA-N sodium;9,10-dioxoanthracene-2-sulfonic acid Chemical compound [Na+].C1=CC=C2C(=O)C3=CC(S(=O)(=O)O)=CC=C3C(=O)C2=C1 GGCZERPQGJTIQP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Radiography Using Non-Light Waves (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
- Image Analysis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
(57)【要約】
【目的】 放射線画像の重ね合せ処理方法および装置に
おいて、画像信号のノイズを最小とし、画質の良好な放
射線画像を得る。 【構成】 同一被写体の放射線画像情報を担持する２つ
の画像信号Ｓ1 、Ｓ2それぞれのノイズをＮ1 、Ｎ2 と
した際に、 Ｓ＝（Ｎ2 2 ／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ1＋（Ｎ1 2
／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ2 なる重み付け加算を行なって、加算信号Ｓを得る。また
２つの画像信号Ｓ1 、Ｓ2 の感度Ｓk1、Ｓk2を求め、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算を行なうようにしてもよい。また、画
素毎に感度Ｓk1，Ｓk2を求めるようにしてもよく、画素
毎の感度Ｓk1，Ｓk2を画像信号Ｓ1 ，Ｓ2 のボケ信号に
より求めてもよい。
おいて、画像信号のノイズを最小とし、画質の良好な放
射線画像を得る。 【構成】 同一被写体の放射線画像情報を担持する２つ
の画像信号Ｓ1 、Ｓ2それぞれのノイズをＮ1 、Ｎ2 と
した際に、 Ｓ＝（Ｎ2 2 ／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ1＋（Ｎ1 2
／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ2 なる重み付け加算を行なって、加算信号Ｓを得る。また
２つの画像信号Ｓ1 、Ｓ2 の感度Ｓk1、Ｓk2を求め、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算を行なうようにしてもよい。また、画
素毎に感度Ｓk1，Ｓk2を求めるようにしてもよく、画素
毎の感度Ｓk1，Ｓk2を画像信号Ｓ1 ，Ｓ2 のボケ信号に
より求めてもよい。
Description
【０００１】
【産業上の利用分野】本発明は、同一被写体の放射線画
像情報を担持する２つの画像信号の加算処理を行う放射
線画像の重ね合せ処理方法および装置に関するものであ
る。
像情報を担持する２つの画像信号の加算処理を行う放射
線画像の重ね合せ処理方法および装置に関するものであ
る。
【０００２】
【従来の技術】記録された放射線画像を読み取って画像
信号を得、この画像信号に適切な画像処理を施した後、
画像を再生記録することが種々の分野で行われている。
たとえば、後の画像処理に適合するように設計されたガ
ンマ値の低いＸ線フイルムを用いてＸ線画像を記録し、
このＸ線画像が記録されたフイルムからＸ線画像を読み
取って電気信号に変換し、この電気信号（画像信号）に
画像処理を施した後コピー写真等に可視像として再生す
ることにより、コントラスト，シャープネス，粒状性等
の画質性能の良好な再生画像を得ることの出来るシステ
ムが開発されている（特公昭61−5193号公報参照）。
信号を得、この画像信号に適切な画像処理を施した後、
画像を再生記録することが種々の分野で行われている。
たとえば、後の画像処理に適合するように設計されたガ
ンマ値の低いＸ線フイルムを用いてＸ線画像を記録し、
このＸ線画像が記録されたフイルムからＸ線画像を読み
取って電気信号に変換し、この電気信号（画像信号）に
画像処理を施した後コピー写真等に可視像として再生す
ることにより、コントラスト，シャープネス，粒状性等
の画質性能の良好な再生画像を得ることの出来るシステ
ムが開発されている（特公昭61−5193号公報参照）。
【０００３】また本出願人により、放射線（Ｘ線，α
線，β線，γ線，電子線，紫外線等）を照射するとこの
放射線エネルギーの一部が蓄積され、その後可視光等の
励起光を照射すると蓄積されたエネルギーに応じた光量
の輝尽発光光を放射する蓄積性蛍光体（輝尽性蛍光体）
を利用して、人体等の被写体の放射線画像を一旦シート
状の蓄積性蛍光体に撮影記録し、蓄積性蛍光体シートを
レーザ光等の励起光で走査して輝尽発光光を生ぜしめ、
得られた輝尽発光光を光電的に読み取って画像信号を
得、この画像信号に基づいて被写体の放射線画像を写真
感光材料等の記録材料、ＣＲＴ等に可視像として出力さ
せる放射線記録再生システムがすでに提案されている
（特開昭55-12429号，同56-11395号，同55-0163472号，
同56-164645 号，同55-116340 号等）。
線，β線，γ線，電子線，紫外線等）を照射するとこの
放射線エネルギーの一部が蓄積され、その後可視光等の
励起光を照射すると蓄積されたエネルギーに応じた光量
の輝尽発光光を放射する蓄積性蛍光体（輝尽性蛍光体）
を利用して、人体等の被写体の放射線画像を一旦シート
状の蓄積性蛍光体に撮影記録し、蓄積性蛍光体シートを
レーザ光等の励起光で走査して輝尽発光光を生ぜしめ、
得られた輝尽発光光を光電的に読み取って画像信号を
得、この画像信号に基づいて被写体の放射線画像を写真
感光材料等の記録材料、ＣＲＴ等に可視像として出力さ
せる放射線記録再生システムがすでに提案されている
（特開昭55-12429号，同56-11395号，同55-0163472号，
同56-164645 号，同55-116340 号等）。
【０００４】このシステムは、従来の銀塩写真を用いる
放射線写真システムと比較して極めて広い放射線露光域
にわたって画像を記録し得るという実用的な利点を有し
ている。すなわち、放射線露光量に対する、蓄積後に励
起によって発光する輝尽発光光の光量が極めて広い範囲
に渡って比例することが認められており、従って種々の
撮影条件により放射線露光量がかなり大幅に変動して
も、蓄積性蛍光体シートより放射される輝尽発光光を読
取りゲインを適当な値に設定して光電変換手段により読
み取って電気信号（画像信号）に変換し、この画像信号
を用いて写真感光材料、ＣＲＴ等の表示装置に放射線画
像を可視像として出力することによって、放射線露光量
の変動に影響されない放射線画像を得ることができる。
放射線写真システムと比較して極めて広い放射線露光域
にわたって画像を記録し得るという実用的な利点を有し
ている。すなわち、放射線露光量に対する、蓄積後に励
起によって発光する輝尽発光光の光量が極めて広い範囲
に渡って比例することが認められており、従って種々の
撮影条件により放射線露光量がかなり大幅に変動して
も、蓄積性蛍光体シートより放射される輝尽発光光を読
取りゲインを適当な値に設定して光電変換手段により読
み取って電気信号（画像信号）に変換し、この画像信号
を用いて写真感光材料、ＣＲＴ等の表示装置に放射線画
像を可視像として出力することによって、放射線露光量
の変動に影響されない放射線画像を得ることができる。
【０００５】一方、従来より放射線画像の重ね合せ処理
が公知となっている（例えば特開昭56-11399号参照）。
一般に、放射線画像は診断用その他の目的に使われる
が、その使用に当たっては被写体の微小な放射線吸収差
を良好に検出することが要求される。放射線画像におけ
るこの検出の程度をコントラスト検出能または単に検出
能と呼ぶが、この検出能の高いもの程診断性能も高く、
実用的価値が高い放射線画像であると言うことができ
る。したがって診断性能を高めるため、この検出能を高
くすることが望まれるが、その最も大きな障害要因は各
種ノイズである。
が公知となっている（例えば特開昭56-11399号参照）。
一般に、放射線画像は診断用その他の目的に使われる
が、その使用に当たっては被写体の微小な放射線吸収差
を良好に検出することが要求される。放射線画像におけ
るこの検出の程度をコントラスト検出能または単に検出
能と呼ぶが、この検出能の高いもの程診断性能も高く、
実用的価値が高い放射線画像であると言うことができ
る。したがって診断性能を高めるため、この検出能を高
くすることが望まれるが、その最も大きな障害要因は各
種ノイズである。
【０００６】例えば、前述した蓄積性蛍光体シートを使
用する放射線画像記録方式においては、放射線画像を蓄
積性蛍光体シートに蓄積記録し、読み出すステップにお
いて次のようなノイズの存在が認められている。
用する放射線画像記録方式においては、放射線画像を蓄
積性蛍光体シートに蓄積記録し、読み出すステップにお
いて次のようなノイズの存在が認められている。
【０００７】(1) 放射線の量子ノイズ (2) 蓄積性蛍光体シートの蛍光体塗布分布もしくは蛍光
体粒子分布の不均一によるノイズ (3) 蓄積性蛍光体シートに蓄積記録された画像を輝尽発
光させる励起光のノイズ (4) 蓄積性蛍光体シートから発せられ、集光、検出され
る輝尽発光光のノイズ (5) 電気信号を増幅、処理する系における電気的ノイズ 重ね合せ処理は、これらのノイズを大幅に減少させ、被
写体の僅かな放射線吸収差も最終画像において明確に観
察可能にする、すなわち検出能を大幅に向上させる方法
である。重ね合せ処理の一般的な手法および作用は、次
の通りである。
体粒子分布の不均一によるノイズ (3) 蓄積性蛍光体シートに蓄積記録された画像を輝尽発
光させる励起光のノイズ (4) 蓄積性蛍光体シートから発せられ、集光、検出され
る輝尽発光光のノイズ (5) 電気信号を増幅、処理する系における電気的ノイズ 重ね合せ処理は、これらのノイズを大幅に減少させ、被
写体の僅かな放射線吸収差も最終画像において明確に観
察可能にする、すなわち検出能を大幅に向上させる方法
である。重ね合せ処理の一般的な手法および作用は、次
の通りである。
【０００８】複数枚重ねた記録媒体に放射線画像を撮影
（記録）し、この複数枚の記録媒体を読取処理にかけて
得た複数の画像信号を重ね合わせる。このことにより、
前述の各種ノイズを減少させることができる。すなわ
ち、前述の蓄積性蛍光体シートのノイズ(1) 〜(5) は各
シートの画像毎に異なった分布を示す場合が多いので、
これらのシートの画像を重ね合わせることにより各ノイ
ズは平均化され、重ね合せ処理をした画像ではノイズが
目立たなくなる。つまり、Ｓ／Ｎの良い画像信号が得ら
れる。Ｘ線フイルムに記録された放射線画像を読み取っ
た場合にも、これと同様のことがいえる。さらに詳しく
は、ノイズ(1) 〜(5) には、ポアソン統計で近似できる
ノイズが多く、特に放射線画像のノイズの中で支配的な
要因の１つである(1) 放射線のノイズはその一例であ
る。ここで、ノイズがポアソン統計で近似できるとし、
２枚の放射線画像がそれぞれ同等の大きさの信号Ｓ1 、
Ｓ2 およびノイズＮ1 、Ｎ2 を持つと考えた場合、２枚
の画像を重ね合せた場合の信号とノイズの大きさは、信
号がＳ1 ＋Ｓ2 、ノイズが
（記録）し、この複数枚の記録媒体を読取処理にかけて
得た複数の画像信号を重ね合わせる。このことにより、
前述の各種ノイズを減少させることができる。すなわ
ち、前述の蓄積性蛍光体シートのノイズ(1) 〜(5) は各
シートの画像毎に異なった分布を示す場合が多いので、
これらのシートの画像を重ね合わせることにより各ノイ
ズは平均化され、重ね合せ処理をした画像ではノイズが
目立たなくなる。つまり、Ｓ／Ｎの良い画像信号が得ら
れる。Ｘ線フイルムに記録された放射線画像を読み取っ
た場合にも、これと同様のことがいえる。さらに詳しく
は、ノイズ(1) 〜(5) には、ポアソン統計で近似できる
ノイズが多く、特に放射線画像のノイズの中で支配的な
要因の１つである(1) 放射線のノイズはその一例であ
る。ここで、ノイズがポアソン統計で近似できるとし、
２枚の放射線画像がそれぞれ同等の大きさの信号Ｓ1 、
Ｓ2 およびノイズＮ1 、Ｎ2 を持つと考えた場合、２枚
の画像を重ね合せた場合の信号とノイズの大きさは、信
号がＳ1 ＋Ｓ2 、ノイズが
【０００９】
【数１】
【００１０】となる。一方、放射線画像の検出能を表わ
す一つの指標であるＳ／Ｎを考えた場合、重ね合せる前
の各画像のＳ／Ｎはそれぞれ、Ｓ1 ／Ｎ1 、Ｓ2 ／Ｎ2
であるが、重ね合せ処理を行なうことによりＳ／Ｎは、
す一つの指標であるＳ／Ｎを考えた場合、重ね合せる前
の各画像のＳ／Ｎはそれぞれ、Ｓ1 ／Ｎ1 、Ｓ2 ／Ｎ2
であるが、重ね合せ処理を行なうことによりＳ／Ｎは、
【００１１】
【数２】
【００１２】となり、Ｓ／Ｎが向上する。また、重ね合
せ処理を行なう際に、それぞれの信号に重み付けを行な
うことにより、Ｓ／Ｎ向上の最適化が可能である。
せ処理を行なう際に、それぞれの信号に重み付けを行な
うことにより、Ｓ／Ｎ向上の最適化が可能である。
【００１３】従来、実際にこの重ね合せ処理を行なうた
めには、例えば、蓄積性蛍光体シートを用いた場合に
は、カセッテに蓄積性蛍光体シートを２枚重ねて入れて
被写体の撮影を行ない、２枚の蓄積性蛍光体シートに対
して通常の読取処理と同様の読取処理を逐次行なって、
２組の画像信号を得る、という方法が用いられている。
めには、例えば、蓄積性蛍光体シートを用いた場合に
は、カセッテに蓄積性蛍光体シートを２枚重ねて入れて
被写体の撮影を行ない、２枚の蓄積性蛍光体シートに対
して通常の読取処理と同様の読取処理を逐次行なって、
２組の画像信号を得る、という方法が用いられている。
【００１４】
【発明が解決しようとする課題】上述した放射線画像の
重ね合せ処理は、同一被写体の放射線画像情報を担持す
る２組の画像信号のそれぞれを適当な値で重み付けした
上で、対応する画素間で加算して、加算信号を求めてい
る。
重ね合せ処理は、同一被写体の放射線画像情報を担持す
る２組の画像信号のそれぞれを適当な値で重み付けした
上で、対応する画素間で加算して、加算信号を求めてい
る。
【００１５】ところが、この画像信号を重み付けする適
当な値は、画像信号毎にマニュアルで適宜設定した値ま
たは予め決められたある一定値であったため、重ね合せ
られる放射線画像に適した重み付けを行なうことができ
ず、得られた重ね合せ画像の画質を向上させることは困
難であった。
当な値は、画像信号毎にマニュアルで適宜設定した値ま
たは予め決められたある一定値であったため、重ね合せ
られる放射線画像に適した重み付けを行なうことができ
ず、得られた重ね合せ画像の画質を向上させることは困
難であった。
【００１６】本発明は上記事情に鑑み、２組の画像信号
の最適な重み付けの係数を求め、重ね合せ画像の画質を
向上させることができる放射線画像の重ね合せ処理方法
および装置を提供することを目的とするものである。
の最適な重み付けの係数を求め、重ね合せ画像の画質を
向上させることができる放射線画像の重ね合せ処理方法
および装置を提供することを目的とするものである。
【００１７】
【課題を解決するための手段】本発明による第１の放射
線画像の重ね合せ処理方法は、同一被写体の放射線画像
情報を担持する第１および第２の画像信号に対して相対
応する画素についての信号間で加算を行って加算信号を
得る放射線画像の重ね合せ処理方法において、前記第１
の画像信号をＳ1 、該第１の画像信号のノイズをＮ1 、
前記第２の画像信号をＳ2 、該第２の画像信号のノイズ
をＮ2 および前記加算信号をＳとしたときに、 Ｓ＝（Ｎ2 2 ／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ1＋（Ｎ1 2
／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ2 なる重み付け加算により、前記加算信号Ｓを求めること
を特徴とするものである。
線画像の重ね合せ処理方法は、同一被写体の放射線画像
情報を担持する第１および第２の画像信号に対して相対
応する画素についての信号間で加算を行って加算信号を
得る放射線画像の重ね合せ処理方法において、前記第１
の画像信号をＳ1 、該第１の画像信号のノイズをＮ1 、
前記第２の画像信号をＳ2 、該第２の画像信号のノイズ
をＮ2 および前記加算信号をＳとしたときに、 Ｓ＝（Ｎ2 2 ／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ1＋（Ｎ1 2
／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ2 なる重み付け加算により、前記加算信号Ｓを求めること
を特徴とするものである。
【００１８】本発明による第１の放射線画像の重ね合せ
処理装置は、上述した本発明による第１の放射線画像の
重ね合せ処理方法を実施するための装置であって、同一
被写体の放射線画像情報を担持する第１および第２の画
像信号に対して相対応する画素についての信号間で加算
を行って加算信号を得る加算手段からなる放射線画像の
重ね合せ処理装置において、前記加算手段が前記第１の
画像信号をＳ1 、該第１の画像信号のノイズをＮ1、前
記第２の画像信号をＳ2 、該第２の画像信号のノイズを
Ｎ2 および前記加算信号をＳとしたときに、 Ｓ＝（Ｎ2 2 ／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ1＋（Ｎ1 2
／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ2 なる重み付け加算を行う加算手段であることを特徴とす
るものである。
処理装置は、上述した本発明による第１の放射線画像の
重ね合せ処理方法を実施するための装置であって、同一
被写体の放射線画像情報を担持する第１および第２の画
像信号に対して相対応する画素についての信号間で加算
を行って加算信号を得る加算手段からなる放射線画像の
重ね合せ処理装置において、前記加算手段が前記第１の
画像信号をＳ1 、該第１の画像信号のノイズをＮ1、前
記第２の画像信号をＳ2 、該第２の画像信号のノイズを
Ｎ2 および前記加算信号をＳとしたときに、 Ｓ＝（Ｎ2 2 ／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ1＋（Ｎ1 2
／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ2 なる重み付け加算を行う加算手段であることを特徴とす
るものである。
【００１９】本発明による第２の放射線画像の重ね合せ
処理方法は、同一被写体の放射線画像情報を担持する第
１および第２の画像信号に対して相対応する画素につい
ての信号間で加算を行って加算信号を得る放射線画像の
重ね合せ処理方法において、前記第１の画像信号を
Ｓ1 、該第１の画像信号の感度をＳk1、前記第２の画像
信号をＳ２、該第２の画像信号の感度をＳk2、前記加算
信号をＳとしたときに、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算により、前記加算信号Ｓを求めること
を特徴とするものである。
処理方法は、同一被写体の放射線画像情報を担持する第
１および第２の画像信号に対して相対応する画素につい
ての信号間で加算を行って加算信号を得る放射線画像の
重ね合せ処理方法において、前記第１の画像信号を
Ｓ1 、該第１の画像信号の感度をＳk1、前記第２の画像
信号をＳ２、該第２の画像信号の感度をＳk2、前記加算
信号をＳとしたときに、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算により、前記加算信号Ｓを求めること
を特徴とするものである。
【００２０】ここで、画像信号の感度とあるが、この感
度は、例えば蓄積性蛍光体シートから読み取った画像信
号の場合は、その画像信号の感度を意味し、またＸ線画
像が記録されたＸ線フイルム等から読み取った画像信号
の場合は、そのＸ線フイルムの感度を意味するものとす
る。
度は、例えば蓄積性蛍光体シートから読み取った画像信
号の場合は、その画像信号の感度を意味し、またＸ線画
像が記録されたＸ線フイルム等から読み取った画像信号
の場合は、そのＸ線フイルムの感度を意味するものとす
る。
【００２１】また、本発明による第２の放射線画像の重
ね合せ処理装置は、上述した本発明による第２の放射線
画像の重ね合せ処理方法を実施するための装置であっ
て、同一被写体の放射線画像情報を担持する第１および
第２の画像信号に対して相対応する画素についての信号
間で加算を行って加算信号を得る加算手段からなる放射
線画像の重ね合せ処理装置において、前記加算手段が前
記第１の画像信号をＳ1 、該第１の画像信号の感度をＳ
k1、前記第２の画像信号をＳ2 、該第２の画像信号の感
度をＳk2および前記加算信号をＳとしたときに、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算を行う加算手段であることを特徴とす
るものである。
ね合せ処理装置は、上述した本発明による第２の放射線
画像の重ね合せ処理方法を実施するための装置であっ
て、同一被写体の放射線画像情報を担持する第１および
第２の画像信号に対して相対応する画素についての信号
間で加算を行って加算信号を得る加算手段からなる放射
線画像の重ね合せ処理装置において、前記加算手段が前
記第１の画像信号をＳ1 、該第１の画像信号の感度をＳ
k1、前記第２の画像信号をＳ2 、該第２の画像信号の感
度をＳk2および前記加算信号をＳとしたときに、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算を行う加算手段であることを特徴とす
るものである。
【００２２】ここで、本発明による第１の重ね合せ処理
方法および装置において、 Ｓ＝（Ｎ2 2 ／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ1 ＋（Ｎ1 2 ／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ2 …(1) なる重み付け加算により、加算信号Ｓを求めているが、
この式(1) における重み付け係数Ｎ2 2 ／（Ｎ1 2 ＋Ｎ
2 2 ）およびＮ1 2 ／（Ｎ1 2 ＋Ｎ2 2 ）は、以下のよ
うにして求めるものである。
方法および装置において、 Ｓ＝（Ｎ2 2 ／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ1 ＋（Ｎ1 2 ／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ2 …(1) なる重み付け加算により、加算信号Ｓを求めているが、
この式(1) における重み付け係数Ｎ2 2 ／（Ｎ1 2 ＋Ｎ
2 2 ）およびＮ1 2 ／（Ｎ1 2 ＋Ｎ2 2 ）は、以下のよ
うにして求めるものである。
【００２３】すなわち、 Ｓ＝ａ×Ｓ1 ＋（１−ａ）×Ｓ2 …(2) 但しａは任意の定数 なる重み付け加算により、加算信号Ｓを求めるとし、画
像信号Ｓ1 のノイズをＮ 1 、画像信号Ｓ2 のノイズをＮ
2 とすると、加算信号Ｓのノイズは、
像信号Ｓ1 のノイズをＮ 1 、画像信号Ｓ2 のノイズをＮ
2 とすると、加算信号Ｓのノイズは、
【００２４】
【数３】
【００２５】で表わされる。
【００２６】したがって、ノイズＮを最小にするような
定数ａを求め、この定数ａを重み付け係数として式(2)
に代入すれば、加算信号Ｓにより得られる放射線画像は
最適なものとなる。
定数ａを求め、この定数ａを重み付け係数として式(2)
に代入すれば、加算信号Ｓにより得られる放射線画像は
最適なものとなる。
【００２７】ここで、ノイズＮを最小にするにはＮ2 を
最小にすればよいから、式(3) を変形して Ｎ2 ＝（ａ×Ｎ1 ）2 ＋（（１−ａ）×Ｎ2 ）2 …(4) とし、式(4) を定数ａで偏微分すると、
最小にすればよいから、式(3) を変形して Ｎ2 ＝（ａ×Ｎ1 ）2 ＋（（１−ａ）×Ｎ2 ）2 …(4) とし、式(4) を定数ａで偏微分すると、
【００２８】
【数４】
【００２９】ノイズＮの２乗を定数ａで偏微分した値が
０となる点においてノイズＮ2 は最小となるから、式
(5) より、 ２ａ×Ｎ1 2 −２（１−ａ）×Ｎ2 2 ＝０ …(6) となる定数ａの値を求めると、 ａ＝Ｎ2 2 ／（Ｎ1 2 ＋Ｎ2 2 ）， （１−ａ）＝Ｎ1 2 ／（Ｎ1 2 ＋Ｎ2 2 ） …
（７） となる。
０となる点においてノイズＮ2 は最小となるから、式
(5) より、 ２ａ×Ｎ1 2 −２（１−ａ）×Ｎ2 2 ＝０ …(6) となる定数ａの値を求めると、 ａ＝Ｎ2 2 ／（Ｎ1 2 ＋Ｎ2 2 ）， （１−ａ）＝Ｎ1 2 ／（Ｎ1 2 ＋Ｎ2 2 ） …
（７） となる。
【００３０】この定数ａの値を式（２） に代入する
と、式(1) が得られる。したがって、式(1) に示す重み
付け加算により、加算信号Ｓを求めれば、加算信号Ｓの
ノイズは最小となるのである。
と、式(1) が得られる。したがって、式(1) に示す重み
付け加算により、加算信号Ｓを求めれば、加算信号Ｓの
ノイズは最小となるのである。
【００３１】また、本発明による第２の重ね合せ処理方
法および装置において、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1 ＋（Ｓk1／（Ｓk1＋Ｓk2））×Ｓ2 …(8) なる重み付け加算により、加算信号Ｓを求めているが、
式(8) において、画像信号に画像処理等を施すためのパ
ラメータである感度Ｓk1およびＳk2を用いているのは以
下の理由によるものである。
法および装置において、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1 ＋（Ｓk1／（Ｓk1＋Ｓk2））×Ｓ2 …(8) なる重み付け加算により、加算信号Ｓを求めているが、
式(8) において、画像信号に画像処理等を施すためのパ
ラメータである感度Ｓk1およびＳk2を用いているのは以
下の理由によるものである。
【００３２】すなわち、画像信号のレベルが最も支配的
な領域においては、ノイズの２乗値が記録媒体に照射さ
れる放射線量とほぼ線形の関係にあり、さらに、感度
は、記録媒体に照射される放射線量との相関関係が強い
ものである。したがって、本発明による第１の重ね合せ
処理方法および装置におけるＮ1 2 およびＮ2 2 を感度
Ｓk1およびＳk2で置き換えることにより、加算信号Ｓの
ノイズを最小となるようにしているのである。
な領域においては、ノイズの２乗値が記録媒体に照射さ
れる放射線量とほぼ線形の関係にあり、さらに、感度
は、記録媒体に照射される放射線量との相関関係が強い
ものである。したがって、本発明による第１の重ね合せ
処理方法および装置におけるＮ1 2 およびＮ2 2 を感度
Ｓk1およびＳk2で置き換えることにより、加算信号Ｓの
ノイズを最小となるようにしているのである。
【００３３】また、ここにいうノイズの２乗値と線形の
関係にある感度とは、図４のグラフ30に示すように、画
像信号から求められる感度の対数をとったものであるた
め、上述した式(8) における感度とは、この対数値をと
った感度を表わすものとなる。したがって、対数値をと
る前の感度を使用して重み付け加算を行う場合、対数値
をとる前の感度をＳk1′, Ｓk2′とすれば、式(8) は、 Ｓ＝（log Ｓk2′／（log Ｓk1′＋log Ｓk2′））×Ｓ1 ＋（log Ｓk1′／（log Ｓk1′＋log Ｓk2′））×Ｓ2 …(8′) のように変形される。
関係にある感度とは、図４のグラフ30に示すように、画
像信号から求められる感度の対数をとったものであるた
め、上述した式(8) における感度とは、この対数値をと
った感度を表わすものとなる。したがって、対数値をと
る前の感度を使用して重み付け加算を行う場合、対数値
をとる前の感度をＳk1′, Ｓk2′とすれば、式(8) は、 Ｓ＝（log Ｓk2′／（log Ｓk1′＋log Ｓk2′））×Ｓ1 ＋（log Ｓk1′／（log Ｓk1′＋log Ｓk2′））×Ｓ2 …(8′) のように変形される。
【００３４】また、本発明による第３の放射線画像の重
ね合せ処理方法は同一被写体の放射線画像情報を担持す
る第１および第２の画像信号に対して相対応する画素に
ついての信号間で加算を行って加算信号Ｓを得る放射線
画像の重ね合せ処理方法において、前記第１の画像信号
Ｓ1 および前記第２の画像信号Ｓ2 に基づいて、前記各
画素毎の前記第１の画像信号の感度Ｓk1および前記第２
の画像信号の感度Ｓk2をそれぞれ求め、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算を前記相対応する画素毎に行うことに
より前記加算信号Ｓを求めることを特徴とするものであ
る。
ね合せ処理方法は同一被写体の放射線画像情報を担持す
る第１および第２の画像信号に対して相対応する画素に
ついての信号間で加算を行って加算信号Ｓを得る放射線
画像の重ね合せ処理方法において、前記第１の画像信号
Ｓ1 および前記第２の画像信号Ｓ2 に基づいて、前記各
画素毎の前記第１の画像信号の感度Ｓk1および前記第２
の画像信号の感度Ｓk2をそれぞれ求め、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算を前記相対応する画素毎に行うことに
より前記加算信号Ｓを求めることを特徴とするものであ
る。
【００３５】さらに、本発明による第３の放射線画像の
重ね合せ処理装置は、同一被写体の放射線画像情報を担
持する第１および第２の画像信号に対して相対応する画
素についての信号間で加算を行って加算信号Ｓを得る加
算手段からなる放射線画像の重ね合せ処理装置におい
て、前記加算手段が、前記第１の画像信号Ｓ1 および前
記第２の画像信号Ｓ2 に基づいて、前記画素毎に前記第
１の画像信号の感度Ｓk1および前記第２の画像信号Ｓk2
をそれぞれ求める演算部と、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算を前記相対応する画素毎に行う加算部
とからなることを特徴とするものである。
重ね合せ処理装置は、同一被写体の放射線画像情報を担
持する第１および第２の画像信号に対して相対応する画
素についての信号間で加算を行って加算信号Ｓを得る加
算手段からなる放射線画像の重ね合せ処理装置におい
て、前記加算手段が、前記第１の画像信号Ｓ1 および前
記第２の画像信号Ｓ2 に基づいて、前記画素毎に前記第
１の画像信号の感度Ｓk1および前記第２の画像信号Ｓk2
をそれぞれ求める演算部と、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算を前記相対応する画素毎に行う加算部
とからなることを特徴とするものである。
【００３６】さらに、本発明による第４の放射線画像の
重ね合せ処理方法は、同一被写体の放射線画像情報を担
持する第１および第２の画像信号に対して相対応する画
素についての信号間で加算を行って加算信号Ｓを得る放
射線画像の重ね合せ処理方法において、前記第１の画像
信号Ｓ1 および前記第２の画像信号Ｓ2 のボケ信号をそ
れぞれ求め、前記第１の画像信号Ｓ1 のボケ信号および
前記第２の画像信号Ｓ2 のボケ信号に基づいて、前記各
画素毎の前記第１の画像信号の感度Ｓk1および前記第２
の画像信号の感度Ｓk2をそれぞれ求め、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算を前記相対応する画素毎に行うことに
より前記加算信号Ｓを求めることを特徴とするものであ
る。
重ね合せ処理方法は、同一被写体の放射線画像情報を担
持する第１および第２の画像信号に対して相対応する画
素についての信号間で加算を行って加算信号Ｓを得る放
射線画像の重ね合せ処理方法において、前記第１の画像
信号Ｓ1 および前記第２の画像信号Ｓ2 のボケ信号をそ
れぞれ求め、前記第１の画像信号Ｓ1 のボケ信号および
前記第２の画像信号Ｓ2 のボケ信号に基づいて、前記各
画素毎の前記第１の画像信号の感度Ｓk1および前記第２
の画像信号の感度Ｓk2をそれぞれ求め、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算を前記相対応する画素毎に行うことに
より前記加算信号Ｓを求めることを特徴とするものであ
る。
【００３７】また、本発明による第４の放射線画像の重
ね合せ処理装置は、同一被写体の放射線画像情報を担持
する第１および第２の画像信号に対して相対応する画素
についての信号間で加算を行って加算信号Ｓを得る加算
手段からなる放射線画像の重ね合せ処理装置において、
前記加算手段が、前記第１の画像信号Ｓ1 および前記第
２の画像信号Ｓ2 のボケ信号を求め、該第１の画像信号
Ｓ1 のボケ信号および前記第２の画像信号Ｓ2のボケ信
号に基づいて、前記画素毎に前記第１の画像信号の感度
Ｓk1および前記第２の画像信号Ｓk2をそれぞれ求める演
算部と、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算を前記相対応する画素毎に行う加算部
とからなることを特徴とするものである。
ね合せ処理装置は、同一被写体の放射線画像情報を担持
する第１および第２の画像信号に対して相対応する画素
についての信号間で加算を行って加算信号Ｓを得る加算
手段からなる放射線画像の重ね合せ処理装置において、
前記加算手段が、前記第１の画像信号Ｓ1 および前記第
２の画像信号Ｓ2 のボケ信号を求め、該第１の画像信号
Ｓ1 のボケ信号および前記第２の画像信号Ｓ2のボケ信
号に基づいて、前記画素毎に前記第１の画像信号の感度
Ｓk1および前記第２の画像信号Ｓk2をそれぞれ求める演
算部と、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算を前記相対応する画素毎に行う加算部
とからなることを特徴とするものである。
【００３８】
【作用】本発明による放射線画像の重ね合せ処理方法お
よび装置は、加算信号を求める際の画像信号の重み付け
係数として、加算信号Ｓのノイズが最小となる重み付け
係数を用いるようにした。このため、加算信号を再生す
ることにより得られる画像のノイズが最小となり、画質
の良い画像を得ることができる。
よび装置は、加算信号を求める際の画像信号の重み付け
係数として、加算信号Ｓのノイズが最小となる重み付け
係数を用いるようにした。このため、加算信号を再生す
ることにより得られる画像のノイズが最小となり、画質
の良い画像を得ることができる。
【００３９】また、加算信号Ｓのノイズが最小となる重
み付け係数を求めるための感度を各画素毎に求めるよう
にすれば、被写体を透過し記録シート等に到達した放射
線の線量の場所毎の違いに影響されず、各画素毎にノイ
ズを最小とすることができるため、より画質を向上させ
ることができる。
み付け係数を求めるための感度を各画素毎に求めるよう
にすれば、被写体を透過し記録シート等に到達した放射
線の線量の場所毎の違いに影響されず、各画素毎にノイ
ズを最小とすることができるため、より画質を向上させ
ることができる。
【００４０】さらに、画像信号のボケ信号を求め、この
ボケ信号に基づいて画素毎の感度を求めるようにすれ
ば、ノイズに影響されない大域的な画素毎の感度を求め
て重み付け係数を決定することができる。
ボケ信号に基づいて画素毎の感度を求めるようにすれ
ば、ノイズに影響されない大域的な画素毎の感度を求め
て重み付け係数を決定することができる。
【００４１】
【実施例】以下、図面を参照して本発明の実施例につい
て説明する。
て説明する。
【００４２】なお、本実施例においては記録媒体として
蓄積性蛍光体シートを用いるものとする。
蓄積性蛍光体シートを用いるものとする。
【００４３】図１は２枚の蓄積性蛍光体シート４Ａ，４
Ｂに、同一の被写体１を透過した放射線２を照射する状
態を示す図である。
Ｂに、同一の被写体１を透過した放射線２を照射する状
態を示す図である。
【００４４】すなわち、第１の蓄積性蛍光体シート４Ａ
および第２の蓄積性蛍光体シート４Ｂとを重ねて配置し
て放射線源３を駆動させて、放射線２を発せしめ、被写
体１を透過した放射線２は第１および第２の蓄積性蛍光
体シート４Ａ，４Ｂに照射され、被写体１の放射線画像
情報が蓄積性蛍光体シート４Ａおよび４Ｂに蓄積記録さ
れる。
および第２の蓄積性蛍光体シート４Ｂとを重ねて配置し
て放射線源３を駆動させて、放射線２を発せしめ、被写
体１を透過した放射線２は第１および第２の蓄積性蛍光
体シート４Ａ，４Ｂに照射され、被写体１の放射線画像
情報が蓄積性蛍光体シート４Ａおよび４Ｂに蓄積記録さ
れる。
【００４５】次にこれら２枚の蓄積性蛍光体シート４
Ａ，４Ｂから、図２に示すような画像読取手段によって
放射線画像を読み取り、放射線画像を表わす画像信号を
得る。まず蓄積性蛍光体シート４Ａをエンドレスベルト
等の副走査手段９により矢印Ｙの方向に移動させなが
ら、レーザー光源10からのレーザー光（励起光）11を走
査ミラー12によって偏向させ、シート４Ａ上をＸ方向に
主走査させる。この励起光走査により蓄積性蛍光体シー
ト４Ａからは、蓄積記録されている放射線画像情報に応
じた光量の輝尽発光光13が発散する。輝尽発光光13は、
透明なアクリル板を成形して作られた光ガイド14の一端
面からこの光ガイド14の内部に入射し、その中を全反射
を繰返しながら進行して、フォトマルチプライヤー（光
電子増倍管）15に受光される。このフォトマルチプライ
ヤー15からは、輝尽発光光13の発光量に対応した、つま
り上記画像情報を示す出力信号ＳA が出力される。
Ａ，４Ｂから、図２に示すような画像読取手段によって
放射線画像を読み取り、放射線画像を表わす画像信号を
得る。まず蓄積性蛍光体シート４Ａをエンドレスベルト
等の副走査手段９により矢印Ｙの方向に移動させなが
ら、レーザー光源10からのレーザー光（励起光）11を走
査ミラー12によって偏向させ、シート４Ａ上をＸ方向に
主走査させる。この励起光走査により蓄積性蛍光体シー
ト４Ａからは、蓄積記録されている放射線画像情報に応
じた光量の輝尽発光光13が発散する。輝尽発光光13は、
透明なアクリル板を成形して作られた光ガイド14の一端
面からこの光ガイド14の内部に入射し、その中を全反射
を繰返しながら進行して、フォトマルチプライヤー（光
電子増倍管）15に受光される。このフォトマルチプライ
ヤー15からは、輝尽発光光13の発光量に対応した、つま
り上記画像情報を示す出力信号ＳA が出力される。
【００４６】この出力信号ＳA は、対数増幅器16により
対数増幅され、次いでＡ／Ｄ変換器17に通されて、デジ
タルの画像信号Ｓ1 に変換される。この画像信号Ｓ1 は
例えば磁気ディスク等の記憶媒体18に記憶される。次
に、全く同様にして、もう１枚の蓄積性蛍光体シート４
Ｂの記録画像情報が読み出され、その画像情報を示すデ
ジタルの画像信号Ｓ2 が同様に記憶媒体18に記憶され
る。
対数増幅され、次いでＡ／Ｄ変換器17に通されて、デジ
タルの画像信号Ｓ1 に変換される。この画像信号Ｓ1 は
例えば磁気ディスク等の記憶媒体18に記憶される。次
に、全く同様にして、もう１枚の蓄積性蛍光体シート４
Ｂの記録画像情報が読み出され、その画像情報を示すデ
ジタルの画像信号Ｓ2 が同様に記憶媒体18に記憶され
る。
【００４７】次に、上述のようにして得られた画像信号
Ｓ1 ，Ｓ2 を用いて重ね合せ処理を行なう。図３はこの
重ね合せ処理を行なう装置を示している。まず記憶媒体
18内の画像ファイル18Ａと、画像ファイル18Ｂから、前
記画像信号Ｓ1 ，Ｓ2 が読み出され、重ね合せ演算回路
19に入力される。重ね合せ演算回路19は、入力された２
つの画像信号Ｓ1 ，Ｓ2 それぞれの感度Ｓk1，Ｓk2が求
められる。ここで感度とは、所定の光量の輝尽発光光を
どのレベルの画像信号とするかを定める光電変換率をい
う。この感度は、画像信号を読み取る際の読取条件や、
画像信号に画像処理を施す際の画像処理条件を決定する
ためのパラメータであり、これを求める演算はあらかじ
め多数の放射線画像を統計的に処理した結果からそのア
ルゴリズムが定められている。このアルゴリズムのひと
つとして、画像信号のヒストグラムを求め、このヒスト
グラムに基づいて読取条件および／または画像処理条件
を求める方法が多数知られている（例えば、特開昭60-1
85944 号公報、特開昭61-280163 号公報等）。したがっ
て、ここでは感度を求めるための詳細な説明は省略す
る。また、ここで求められた感度Ｓk1，Ｓk2は、前述し
た図４のグラフ30に示すように、画像信号から求められ
た感度のＳk1′，Ｓk2′の対数をとったものである。す
なわち Ｓk1＝log Ｓk1′，Ｓk2＝log Ｓk2′ となるものである。
Ｓ1 ，Ｓ2 を用いて重ね合せ処理を行なう。図３はこの
重ね合せ処理を行なう装置を示している。まず記憶媒体
18内の画像ファイル18Ａと、画像ファイル18Ｂから、前
記画像信号Ｓ1 ，Ｓ2 が読み出され、重ね合せ演算回路
19に入力される。重ね合せ演算回路19は、入力された２
つの画像信号Ｓ1 ，Ｓ2 それぞれの感度Ｓk1，Ｓk2が求
められる。ここで感度とは、所定の光量の輝尽発光光を
どのレベルの画像信号とするかを定める光電変換率をい
う。この感度は、画像信号を読み取る際の読取条件や、
画像信号に画像処理を施す際の画像処理条件を決定する
ためのパラメータであり、これを求める演算はあらかじ
め多数の放射線画像を統計的に処理した結果からそのア
ルゴリズムが定められている。このアルゴリズムのひと
つとして、画像信号のヒストグラムを求め、このヒスト
グラムに基づいて読取条件および／または画像処理条件
を求める方法が多数知られている（例えば、特開昭60-1
85944 号公報、特開昭61-280163 号公報等）。したがっ
て、ここでは感度を求めるための詳細な説明は省略す
る。また、ここで求められた感度Ｓk1，Ｓk2は、前述し
た図４のグラフ30に示すように、画像信号から求められ
た感度のＳk1′，Ｓk2′の対数をとったものである。す
なわち Ｓk1＝log Ｓk1′，Ｓk2＝log Ｓk2′ となるものである。
【００４８】感度Ｓk1，Ｓk2が求められると２つの画像
信号Ｓ1 ，Ｓ2 に対して感度Ｓk1，Ｓk2で重み付けをし
て、式 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1 ＋（Ｓk1／（Ｓk1＋Ｓk2））×Ｓ2 …(8) により加算信号Ｓを求める。これは前述した式(1) にお
いて、Ｎ1 2 を感度Ｓk1で、Ｎ2 2 を感度Ｓk2で置換し
たものであるため、この加算信号Ｓはノイズが最小のも
のとなっている。
信号Ｓ1 ，Ｓ2 に対して感度Ｓk1，Ｓk2で重み付けをし
て、式 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1 ＋（Ｓk1／（Ｓk1＋Ｓk2））×Ｓ2 …(8) により加算信号Ｓを求める。これは前述した式(1) にお
いて、Ｎ1 2 を感度Ｓk1で、Ｎ2 2 を感度Ｓk2で置換し
たものであるため、この加算信号Ｓはノイズが最小のも
のとなっている。
【００４９】一方、感度Ｓk1，Ｓk2として対数値をとる
前の感度Ｓk1′，Ｓk2′を使用する場合は、 Ｓ＝（log Ｓk2′／（log Ｓk1′＋log Ｓk2′））×Ｓ1 ＋（log Ｓk1′／（log Ｓk1′＋log Ｓk2′））×Ｓ2 …(8′) により加算信号Ｓを求めるようにしてもよい。
前の感度Ｓk1′，Ｓk2′を使用する場合は、 Ｓ＝（log Ｓk2′／（log Ｓk1′＋log Ｓk2′））×Ｓ1 ＋（log Ｓk1′／（log Ｓk1′＋log Ｓk2′））×Ｓ2 …(8′) により加算信号Ｓを求めるようにしてもよい。
【００５０】さらに、図４の感度Ｓk に対して、便宜上
の目盛りを使用する場合（例えば、図４の横軸Ｓk ″の
ように目盛りを使用する場合）は、Ｓk ″とＳk ′との
関係は、
の目盛りを使用する場合（例えば、図４の横軸Ｓk ″の
ように目盛りを使用する場合）は、Ｓk ″とＳk ′との
関係は、
【００５１】
【数５】
【００５２】となるため、 Ｓ＝（４−Ｓk2″）／（（４−Ｓk1″）＋（４−Ｓk2″））×Ｓ1 ＋（４−Ｓk1″）／（（４−Ｓk1″）＋（４−Ｓk2″））×Ｓ2 …(8″) によって加算信号Ｓを求めるようにしてもよい。
【００５３】このように加算信号Ｓが求められると、加
算信号Ｓは、画像処理回路20において階調処理、周波数
処理等の画像処理を受けた後、画像再生装置21に送られ
て、放射線画像の再生に供せられる。この加算信号Ｓに
基づいて放射線画像を再生すれば、ノイズが最も少ない
画質の良好な放射線画像を得ることができる。
算信号Ｓは、画像処理回路20において階調処理、周波数
処理等の画像処理を受けた後、画像再生装置21に送られ
て、放射線画像の再生に供せられる。この加算信号Ｓに
基づいて放射線画像を再生すれば、ノイズが最も少ない
画質の良好な放射線画像を得ることができる。
【００５４】一方、上述した感度Ｓk1，Ｓk2を画像信号
が得られた各画素毎に求めるようにしてもよい。以下、
画像信号Ｓ1 について各画素毎の感度を求める方法につ
いて説明する。
が得られた各画素毎に求めるようにしてもよい。以下、
画像信号Ｓ1 について各画素毎の感度を求める方法につ
いて説明する。
【００５５】まず、画像信号Ｓ1 から画像信号Ｓ1 の感
度Ｓk1およびラチチュードＧp1を前述した特開昭60-185
944 号公報等に開示された方法により求める。ここでラ
チチュードとは、最も微弱な輝尽発光光に対する最も強
大な輝尽発光光の光量比に対応するものである。例え
ば、感度Ｓk1が２、ラチチュードＧp1が２というように
求められたとすると、感度Ｓk1と画像信号Ｓ1 との関係
は図５に示すものとなる。ここで、画素値が10bit の場
合、Ｓmin ＝０，Ｓmax ＝1023となる。蓄積性蛍光体シ
ート４Ａについて図６に示すようなx-y 座標を定め、蓄
積性蛍光体シート４Ａ上の任意の点(x,y) における感度
をＳk1(x,y) とすると、
度Ｓk1およびラチチュードＧp1を前述した特開昭60-185
944 号公報等に開示された方法により求める。ここでラ
チチュードとは、最も微弱な輝尽発光光に対する最も強
大な輝尽発光光の光量比に対応するものである。例え
ば、感度Ｓk1が２、ラチチュードＧp1が２というように
求められたとすると、感度Ｓk1と画像信号Ｓ1 との関係
は図５に示すものとなる。ここで、画素値が10bit の場
合、Ｓmin ＝０，Ｓmax ＝1023となる。蓄積性蛍光体シ
ート４Ａについて図６に示すようなx-y 座標を定め、蓄
積性蛍光体シート４Ａ上の任意の点(x,y) における感度
をＳk1(x,y) とすると、
【００５６】
【数６】
【００５７】ただし Ｓ1 (x,y) ：蓄積性蛍光体シート
４Ａ上の任意の点(x,y) における画像信号Ｓ1 の値 により求めることができる。
４Ａ上の任意の点(x,y) における画像信号Ｓ1 の値 により求めることができる。
【００５８】式(9) により各画素における感度Ｓk1(x,
y) が求められると、同様に画像信号Ｓ1 についても蓄
積性蛍光体シート４Ｂ上の各画素における感度Ｓk2(x,
y) が求められ、式 Ｓ(x,y) ＝Ｓk2(x,y) ／（Ｓk1(x,y) ＋Ｓk2(x,y) ）×Ｓ1 ＋Ｓk1(x,y) ／（Ｓk1(x,y) ＋Ｓk2(x,y) ）×Ｓ2 …(10) により加算信号Ｓが求められる。
y) が求められると、同様に画像信号Ｓ1 についても蓄
積性蛍光体シート４Ｂ上の各画素における感度Ｓk2(x,
y) が求められ、式 Ｓ(x,y) ＝Ｓk2(x,y) ／（Ｓk1(x,y) ＋Ｓk2(x,y) ）×Ｓ1 ＋Ｓk1(x,y) ／（Ｓk1(x,y) ＋Ｓk2(x,y) ）×Ｓ2 …(10) により加算信号Ｓが求められる。
【００５９】このように、重み付け係数を各画素毎に求
めるようにすれば、被写体を透過し、蓄積性蛍光体シー
トに到達した放射線線量の場所毎の違いに影響されず、
画素毎にノイズを最小とすることができる。
めるようにすれば、被写体を透過し、蓄積性蛍光体シー
トに到達した放射線線量の場所毎の違いに影響されず、
画素毎にノイズを最小とすることができる。
【００６０】また、画素毎の感度を、各画像信号のボケ
信号に基づいて求めるようにすれば、画像信号の細かい
変化やノイズに影響されずに、大域的な感度を求めるこ
とができるため、より好ましい。ここで、ボケ信号は、
ボケマスク信号を用いる方法やメジアンフィルタ、ＦＦ
Ｔフィルタ、Ｖ−フィルタを用いる方法等、各種平滑化
処理を用いて求められる。
信号に基づいて求めるようにすれば、画像信号の細かい
変化やノイズに影響されずに、大域的な感度を求めるこ
とができるため、より好ましい。ここで、ボケ信号は、
ボケマスク信号を用いる方法やメジアンフィルタ、ＦＦ
Ｔフィルタ、Ｖ−フィルタを用いる方法等、各種平滑化
処理を用いて求められる。
【００６１】上述した実施例においては、２つの画像信
号Ｓ1 ，Ｓ2 それぞれの感度Ｓk1，Ｓk2に基づいて重み
付けを行なう際の重み付け係数を定めているが、２つの
画像信号Ｓ1 ，Ｓ2 のそれぞれのノイズＮ1 ，Ｎ2 を検
出し、このノイズＮ1 ，Ｎ2に基づいて前述した式(1)
により重み付け加算を行なうようにしてもよい。
号Ｓ1 ，Ｓ2 それぞれの感度Ｓk1，Ｓk2に基づいて重み
付けを行なう際の重み付け係数を定めているが、２つの
画像信号Ｓ1 ，Ｓ2 のそれぞれのノイズＮ1 ，Ｎ2 を検
出し、このノイズＮ1 ，Ｎ2に基づいて前述した式(1)
により重み付け加算を行なうようにしてもよい。
【００６２】また、上述した実施例においては、特開昭
60-185944 号公報、特開昭61-280163 号公報等に開示さ
れた方法により感度を求めるようにしているが、これに
限定されるものではなく、ニューラルネットワークを用
いた方法等いかなる方法で求めるようにしてもよい。
60-185944 号公報、特開昭61-280163 号公報等に開示さ
れた方法により感度を求めるようにしているが、これに
限定されるものではなく、ニューラルネットワークを用
いた方法等いかなる方法で求めるようにしてもよい。
【００６３】また、上記実施例は蓄積性蛍光体シートか
ら読み取られた画像信号を用いた重ね合せ処理について
説明したが、本発明の対象とする同一被写体の２つの画
像信号はこれに限られるものではなく、Ｘ線フイルムか
ら読み取られた画像信号、イメージインテンシファイア
から得られた画像信号等、その他の記録媒体から得られ
た画像信号を用いるようにしてもよい。
ら読み取られた画像信号を用いた重ね合せ処理について
説明したが、本発明の対象とする同一被写体の２つの画
像信号はこれに限られるものではなく、Ｘ線フイルムか
ら読み取られた画像信号、イメージインテンシファイア
から得られた画像信号等、その他の記録媒体から得られ
た画像信号を用いるようにしてもよい。
【００６４】
【発明の効果】以上、詳細に説明したように、本発明に
よる放射線画像の重ね合せ処理方法および装置は、重み
付け加算処理により得られる加算信号のノイズを最小に
することが可能である。したがって、この加算信号を再
生することにより、粒状性の向上した画質が良好な放射
線画像を得ることが可能となる。
よる放射線画像の重ね合せ処理方法および装置は、重み
付け加算処理により得られる加算信号のノイズを最小に
することが可能である。したがって、この加算信号を再
生することにより、粒状性の向上した画質が良好な放射
線画像を得ることが可能となる。
【００６５】また、重み付け加算処理の重み付け係数に
用いる感度を、各画素毎に求めるようにすれば、被写体
の場所毎の線量の変化に影響されず、より画質が良好な
放射線画像を得ることが可能となる。
用いる感度を、各画素毎に求めるようにすれば、被写体
の場所毎の線量の変化に影響されず、より画質が良好な
放射線画像を得ることが可能となる。
【００６６】さらに、画像信号のボケ信号を求め、この
ボケ信号に基づいて各画素毎の感度を求めるようにすれ
ば、画像信号の細かい変化やノイズに影響されない大域
的な感度を求めることができ、より好ましいものであ
る。
ボケ信号に基づいて各画素毎の感度を求めるようにすれ
ば、画像信号の細かい変化やノイズに影響されない大域
的な感度を求めることができ、より好ましいものであ
る。
【図１】本発明の放射線画像の重ね合せ処理方法におけ
る、放射線画像の蓄積記録ステップを表す図
る、放射線画像の蓄積記録ステップを表す図
【図２】蓄積記録がなされた蓄積性蛍光体シートからの
放射線画像情報読取りを説明する概略図
放射線画像情報読取りを説明する概略図
【図３】重ね合せ処理を実施する装置の概要を示すブロ
ック図
ック図
【図４】感度とノイズとの関係を表わすグラフ
【図５】感度と画像信号との関係を表わすグラフ
【図６】蓄積性蛍光体シートをｘ，ｙ座標として表わし
た図
た図
１ 被写体 ２ 放射線 ３ 放射線源 ４Ａ，４Ｂ 蓄積性蛍光体シート 10 レーザー光源 11 レーザー光 12 走査ミラー 13 輝尽発光光 14 光ガイド 15 フォトマルチプライヤー 16 対数増幅器 17 Ａ／Ｄ変換器 19 重ね合せ演算回路 20 画像処理回路 21 画像再生装置 Ｓ1 ，Ｓ2 デジタル画像信号 Ｓ 加算信号
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.5 識別記号 庁内整理番号 ＦＩ 技術表示箇所 Ｇ０６Ｆ 15/68 ３５０ 8420−5Ｌ Ｈ０４Ｎ 1/40 １０１ Ｃ 9068−5Ｃ
Claims (8)
- 【請求項１】 同一被写体の放射線画像情報を担持する
第１および第２の画像信号に対して相対応する画素につ
いての信号間で加算を行って加算信号を得る放射線画像
の重ね合せ処理方法において、 前記第１の画像信号をＳ1 、該第１の画像信号のノイズ
をＮ1 、前記第２の画像信号をＳ2 、該第２の画像信号
のノイズをＮ2 および前記加算信号をＳとしたときに、 Ｓ＝（Ｎ2 2 ／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ1＋（Ｎ1 2
／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ2 なる重み付け加算により、前記加算信号Ｓを求めること
を特徴とする放射線画像の重ね合せ処理方法。 - 【請求項２】 同一被写体の放射線画像情報を担持する
第１および第２の画像信号に対して相対応する画素につ
いての信号間で加算を行って加算信号を得る加算手段か
らなる放射線画像の重ね合せ処理装置において、 前記加算手段が前記第１の画像信号をＳ1 、該第１の画
像信号のノイズをＮ1、前記第２の画像信号をＳ2 、該
第２の画像信号のノイズをＮ2 および前記加算信号をＳ
としたときに、 Ｓ＝（Ｎ2 2 ／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ1＋（Ｎ1 2
／（Ｎ1 2 ＋Ｎ2 2 ））×Ｓ2 なる重み付け加算を行う加算手段であることを特徴とす
る放射線画像の重ね合せ処理装置。 - 【請求項３】 同一被写体の放射線画像情報を担持する
第１および第２の画像信号に対して相対応する画素につ
いての信号間で加算を行って加算信号を得る放射線画像
の重ね合せ処理方法において、 前記第１の画像信号をＳ1 、該第１の画像信号の感度を
Ｓk1、前記第２の画像信号をＳ２、該第２の画像信号の
感度をＳk2、前記加算信号をＳとしたときに、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1+(Sk1／（Ｓk1＋Ｓ
k2））×Ｓ2 なる重み付け加算により、前記加算信号Ｓを求めること
を特徴とする放射線画像の重ね合せ処理方法。 - 【請求項４】 同一被写体の放射線画像情報を担持する
第１および第２の画像信号に対して相対応する画素につ
いての信号間で加算を行って加算信号を得る加算手段か
らなる放射線画像の重ね合せ処理装置において、 前記加算手段が前記第１の画像信号をＳ1 、該第１の画
像信号の感度をＳk1、前記第２の画像信号をＳ2 、該第
２の画像信号の感度をＳk2および前記加算信号をＳとし
たときに、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算を行う加算手段であることを特徴とす
る放射線画像の重ね合せ処理装置。 - 【請求項５】 同一被写体の放射線画像情報を担持する
第１および第２の画像信号に対して相対応する画素につ
いての信号間で加算を行って加算信号Ｓを得る放射線画
像の重ね合せ処理方法において、 前記第１の画像信号Ｓ1 および前記第２の画像信号Ｓ2
に基づいて、前記各画素毎の前記第１の画像信号の感度
Ｓk1および前記第２の画像信号の感度Ｓk2をそれぞれ求
め、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算を前記相対応する画素毎に行うことに
より前記加算信号Ｓを求めることを特徴とする放射線画
像の重ね合せ処理方法。 - 【請求項６】 同一被写体の放射線画像情報を担持する
第１および第２の画像信号に対して相対応する画素につ
いての信号間で加算を行って加算信号Ｓを得る加算手段
からなる放射線画像の重ね合せ処理装置において、 前記加算手段が、前記第１の画像信号Ｓ1 および前記第
２の画像信号Ｓ2 に基づいて、前記画素毎に前記第１の
画像信号の感度Ｓk1および前記第２の画像信号Ｓk2をそ
れぞれ求める演算部と、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算を前記相対応する画素毎に行う加算部
とからなることを特徴とする放射線画像の重ね合せ処理
装置。 - 【請求項７】 同一被写体の放射線画像情報を担持する
第１および第２の画像信号に対して相対応する画素につ
いての信号間で加算を行って加算信号Ｓを得る放射線画
像の重ね合せ処理方法において、 前記第１の画像信号Ｓ1 および前記第２の画像信号Ｓ2
のボケ信号をそれぞれ求め、 前記第１の画像信号Ｓ1 のボケ信号および前記第２の画
像信号Ｓ2 のボケ信号に基づいて、前記各画素毎の前記
第１の画像信号の感度Ｓk1および前記第２の画像信号の
感度Ｓk2をそれぞれ求め、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算を前記相対応する画素毎に行うことに
より前記加算信号Ｓを求めることを特徴とする放射線画
像の重ね合せ処理方法。 - 【請求項８】 同一被写体の放射線画像情報を担持する
第１および第２の画像信号に対して相対応する画素につ
いての信号間で加算を行って加算信号Ｓを得る加算手段
からなる放射線画像の重ね合せ処理装置において、 前記加算手段が、前記第１の画像信号Ｓ1 および前記第
２の画像信号Ｓ2 のボケ信号を求め、該第１の画像信号
Ｓ1 のボケ信号および前記第２の画像信号Ｓ2のボケ信
号に基づいて、前記画素毎に前記第１の画像信号の感度
Ｓk1および前記第２の画像信号Ｓk2をそれぞれ求める演
算部と、 Ｓ＝（Ｓk2／（Ｓk1＋Ｓk2））×Ｓ1＋（Ｓk1／（Ｓk1
＋Ｓk2））×Ｓ2 なる重み付け加算を前記相対応する画素毎に行う加算部
とからなることを特徴とする放射線画像の重ね合せ処理
装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US07/997,533 US5402338A (en) | 1991-12-26 | 1992-12-28 | Method for forming energy subtraction images |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP34704491 | 1991-12-27 | ||
| JP3-347044 | 1991-12-27 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH05232608A true JPH05232608A (ja) | 1993-09-10 |
| JP2715212B2 JP2715212B2 (ja) | 1998-02-18 |
Family
ID=18387545
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP4045529A Expired - Fee Related JP2715212B2 (ja) | 1991-12-26 | 1992-03-03 | 放射線画像の重ね合せ処理方法および装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2715212B2 (ja) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009544001A (ja) * | 2006-05-24 | 2009-12-10 | デュール デンタル アクチェンゲゼルシャフト | 放射線画像を作成するための記憶ユニット及びかかる記憶ユニットを読みとる方法 |
| JP2010243470A (ja) * | 2009-04-09 | 2010-10-28 | Toshiba It & Control Systems Corp | Ｃｔ装置 |
| JP2018202009A (ja) * | 2017-06-07 | 2018-12-27 | キヤノンメディカルシステムズ株式会社 | 医用画像診断装置、医用画像処理装置及び医用画像処理プログラム |
-
1992
- 1992-03-03 JP JP4045529A patent/JP2715212B2/ja not_active Expired - Fee Related
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009544001A (ja) * | 2006-05-24 | 2009-12-10 | デュール デンタル アクチェンゲゼルシャフト | 放射線画像を作成するための記憶ユニット及びかかる記憶ユニットを読みとる方法 |
| JP2010243470A (ja) * | 2009-04-09 | 2010-10-28 | Toshiba It & Control Systems Corp | Ｃｔ装置 |
| JP2018202009A (ja) * | 2017-06-07 | 2018-12-27 | キヤノンメディカルシステムズ株式会社 | 医用画像診断装置、医用画像処理装置及び医用画像処理プログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2715212B2 (ja) | 1998-02-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US4931644A (en) | Method of adjusting radiation image read-out conditions | |
| US5402338A (en) | Method for forming energy subtraction images | |
| JPH07159910A (ja) | 放射線画像の重ね合わせ処理方法および装置 | |
| JP2715212B2 (ja) | 放射線画像の重ね合せ処理方法および装置 | |
| JPH04180056A (ja) | 放射線画像処理方法 | |
| JP3454318B2 (ja) | 画像重ね合せ方法およびエネルギーサブトラクション方法 | |
| JP3432917B2 (ja) | 画像重ね合せ方法およびエネルギーサブトラクション方法 | |
| JPH02266477A (ja) | 偽画像信号の検出方法 | |
| JPH01237887A (ja) | 放射線画像の粒状雑音低減方法 | |
| JP2824878B2 (ja) | 放射線画像のエネルギーサブトラクション撮影処理方法 | |
| JP3454319B2 (ja) | 画像重ね合せ方法およびエネルギーサブトラクション方法 | |
| JPH01238650A (ja) | 放射線画像の被写体認識方法 | |
| JPH0533379B2 (ja) | ||
| CA1237540A (en) | Method of adjusting radiation image read-out conditions | |
| JPH0314372A (ja) | 放射線画像の時間差サブトラクション方法および放射線画像処理装置 | |
| JPH0584504B2 (ja) | ||
| JPH0786925B2 (ja) | 放射線画像処理方法および装置 | |
| JPH04212949A (ja) | 放射線画像読取条件及び／又は画像処理条件決定方法および装置 | |
| JPH01238649A (ja) | 放射線画像の被写体認識方法 | |
| JPH04324438A (ja) | 放射線画像読取方法 | |
| JPH07111733B2 (ja) | 放射線画像のエネルギーサブトラクション方法および装置 | |
| JPH02272681A (ja) | 放射線画像処理装置 | |
| JPH06130518A (ja) | 放射線画像読取条件決定方法および装置 | |
| JPH0584500B2 (ja) | ||
| JPH04310944A (ja) | 放射線照射野認識方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 19970924 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20071107 Year of fee payment: 10 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20071107 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081107 Year of fee payment: 11 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091107 Year of fee payment: 12 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091107 Year of fee payment: 12 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101107 Year of fee payment: 13 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111107 Year of fee payment: 14 |
|
| LAPS | Cancellation because of no payment of annual fees |