JPH06274614A - 画像処理方法 - Google Patents
画像処理方法Info
- Publication number
- JPH06274614A JPH06274614A JP5057590A JP5759093A JPH06274614A JP H06274614 A JPH06274614 A JP H06274614A JP 5057590 A JP5057590 A JP 5057590A JP 5759093 A JP5759093 A JP 5759093A JP H06274614 A JPH06274614 A JP H06274614A
- Authority
- JP
- Japan
- Prior art keywords
- image
- signal
- wavelet transform
- signals
- image 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.)
- Pending
Links
- 238000003672 processing method Methods 0.000 title claims description 16
- 230000005855 radiation Effects 0.000 claims description 31
- 230000001131 transforming effect Effects 0.000 claims 1
- 230000000007 visual effect Effects 0.000 abstract description 6
- 230000009466 transformation Effects 0.000 abstract 4
- 229910052824 wavellite Inorganic materials 0.000 abstract 1
- 230000006870 function Effects 0.000 description 56
- 238000000034 method Methods 0.000 description 20
- 238000001914 filtration Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 13
- 230000008569 process Effects 0.000 description 11
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 5
- 239000000463 material Substances 0.000 description 5
- 230000005284 excitation Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- NIXOWILDQLNWCW-UHFFFAOYSA-N acrylic acid group Chemical group C(C=C)(=O)O NIXOWILDQLNWCW-UHFFFAOYSA-N 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013144 data compression Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000004020 luminiscence type Methods 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- 238000001454 recorded image Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 229910052709 silver Inorganic materials 0.000 description 1
- 239000004332 silver Substances 0.000 description 1
- -1 silver halide Chemical class 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Complex Calculations (AREA)
- Image Processing (AREA)
Abstract
(57)【要約】
【目的】 放射線画像を表す画像信号に画像処理を施す
画像処理方法において、観察読影適性に優れ、視覚的な
印象にあった自然な強調等が行えるようにする。 【構成】 放射線画像を表す画像信号１に対して、ウェ
ーブレット変換２を施して、この画像信号１を複数の周
波数帯域の信号３に分解する。次いで、この信号３のう
ち、少なくとも１つの周波数帯域の信号３に対して所定
の画像処理４を施し、この画像処理４が施された信号３
と他の信号に逆ウェーブレット変換５を施す。
画像処理方法において、観察読影適性に優れ、視覚的な
印象にあった自然な強調等が行えるようにする。 【構成】 放射線画像を表す画像信号１に対して、ウェ
ーブレット変換２を施して、この画像信号１を複数の周
波数帯域の信号３に分解する。次いで、この信号３のう
ち、少なくとも１つの周波数帯域の信号３に対して所定
の画像処理４を施し、この画像処理４が施された信号３
と他の信号に逆ウェーブレット変換５を施す。
Description
【０００１】
【産業上の利用分野】本発明は、放射線画像を表す画像
信号に基づいて可視像を再生する際の画像処理方法に関
するものである。
信号に基づいて可視像を再生する際の画像処理方法に関
するものである。
【０００２】
【従来の技術】記録された画像を読み取って画像信号を
得、この画像信号に適切な画像処理を施した後、画像を
再生表示することは種々の分野で行われている。たとえ
ば、後の画像処理に適合するように設計されたガンマ値
の低いフイルムを用いてＸ線画像を記録し、このＸ線画
像が記録されたフイルムからＸ線画像を読み取って電気
信号（画像信号）に変換し、この画像信号に画像処理を
施した後コピー写真等に可視像として再生することによ
り、コントラスト，シャープネス，粒状性等の画質性能
の良好な再生画像を得ることが行われている（特公昭61
-5193 号公報参照）。
得、この画像信号に適切な画像処理を施した後、画像を
再生表示することは種々の分野で行われている。たとえ
ば、後の画像処理に適合するように設計されたガンマ値
の低いフイルムを用いてＸ線画像を記録し、このＸ線画
像が記録されたフイルムからＸ線画像を読み取って電気
信号（画像信号）に変換し、この画像信号に画像処理を
施した後コピー写真等に可視像として再生することによ
り、コントラスト，シャープネス，粒状性等の画質性能
の良好な再生画像を得ることが行われている（特公昭61
-5193 号公報参照）。
【０００３】また本願出願人により、放射線（Ｘ線，α
線，β線，γ線，電子線，紫外線等）を照射するとこの
放射線エネルギーの一部が蓄積され、その後可視光等の
励起光を照射すると蓄積されたエネルギーに応じて輝尽
発光を示す蓄積性蛍光体（輝尽性蛍光体）を利用して、
人体等の被写体の放射線画像情報をシート状の蓄積性蛍
光体に一旦記録し、この蓄積性蛍光体シートをレーザー
光等の励起光で走査して輝尽発光光を生ぜしめ、得られ
た輝尽発光光を光電的に読み取って画像信号を得、この
画像信号に基づき被写体の放射線画像を写真感光材料等
の記録材料、ＣＲＴデイスプレイ装置等に可視画像とし
て出力させる放射線画像記録再生システムがすでに提案
されている（特開昭55-12429号，同56-11395号，同55-1
63472 号，同56-104645 号，同55- 116340号等）。
線，β線，γ線，電子線，紫外線等）を照射するとこの
放射線エネルギーの一部が蓄積され、その後可視光等の
励起光を照射すると蓄積されたエネルギーに応じて輝尽
発光を示す蓄積性蛍光体（輝尽性蛍光体）を利用して、
人体等の被写体の放射線画像情報をシート状の蓄積性蛍
光体に一旦記録し、この蓄積性蛍光体シートをレーザー
光等の励起光で走査して輝尽発光光を生ぜしめ、得られ
た輝尽発光光を光電的に読み取って画像信号を得、この
画像信号に基づき被写体の放射線画像を写真感光材料等
の記録材料、ＣＲＴデイスプレイ装置等に可視画像とし
て出力させる放射線画像記録再生システムがすでに提案
されている（特開昭55-12429号，同56-11395号，同55-1
63472 号，同56-104645 号，同55- 116340号等）。
【０００４】このシステムは、従来の銀塩写真を用いる
放射線写真システムと比較して極めて広い放射線露出域
にわたって画像を記録しうるという実用的な利点を有し
ている。すなわち、蓄積性蛍光体においては、放射線露
光量に対して蓄積後に励起によって輝尽発光する発光光
の光量が極めて広い範囲にわたって比例することが認め
られており、従って種々の撮影条件により放射線露光量
がかなり大幅に変動しても、蓄積性蛍光体シートより放
射される輝尽発光光の光量を読取ゲインを適当な値に設
定して光電変換手段により読み取って電気信号（画像信
号）に変換し、この画像信号を用いて写真感光材料等の
記録材料、ＣＲＴデイスプレイ装置等の表示装置に放射
線画像を可視画像として出力させることによって、放射
線露光量の変動に影響されない放射線画像を得ることが
できる。
放射線写真システムと比較して極めて広い放射線露出域
にわたって画像を記録しうるという実用的な利点を有し
ている。すなわち、蓄積性蛍光体においては、放射線露
光量に対して蓄積後に励起によって輝尽発光する発光光
の光量が極めて広い範囲にわたって比例することが認め
られており、従って種々の撮影条件により放射線露光量
がかなり大幅に変動しても、蓄積性蛍光体シートより放
射される輝尽発光光の光量を読取ゲインを適当な値に設
定して光電変換手段により読み取って電気信号（画像信
号）に変換し、この画像信号を用いて写真感光材料等の
記録材料、ＣＲＴデイスプレイ装置等の表示装置に放射
線画像を可視画像として出力させることによって、放射
線露光量の変動に影響されない放射線画像を得ることが
できる。
【０００５】上記各種システムにおいて、前述したよう
に観察・診断に適した可視画像を得るために得られた画
像信号に各種の画像処理を施すことが一般的に行われて
いる。この画像処理としては、画像信号のレベルにより
再生される放射線画像の濃度を調整する階調処理が知ら
れている。
に観察・診断に適した可視画像を得るために得られた画
像信号に各種の画像処理を施すことが一般的に行われて
いる。この画像処理としては、画像信号のレベルにより
再生される放射線画像の濃度を調整する階調処理が知ら
れている。
【０００６】また、放射線画像の診断性能を向上させる
ために、画像信号に対してボケマスク処理等の周波数強
調処理を施す方法が本出願人により提案されている（特
開昭55-163472 号公報、特開昭55-87953号公報等）。こ
の周波数処理は、読取画像信号Ｓorg に、この読取画像
信号Ｓorg からボケマスク信号Ｓusを引いたものに強調
度βを掛けたものを加える処理を施すもので、これによ
り、画像において所定の空間周波数成分を強調するよう
にしたものである。これを式で表すと下記の式のように
なる。
ために、画像信号に対してボケマスク処理等の周波数強
調処理を施す方法が本出願人により提案されている（特
開昭55-163472 号公報、特開昭55-87953号公報等）。こ
の周波数処理は、読取画像信号Ｓorg に、この読取画像
信号Ｓorg からボケマスク信号Ｓusを引いたものに強調
度βを掛けたものを加える処理を施すもので、これによ
り、画像において所定の空間周波数成分を強調するよう
にしたものである。これを式で表すと下記の式のように
なる。
【０００７】 Ｓ＝Ｓorg ＋β（Ｓorg −Ｓus） …(1) （Ｓ：周波数処理された信号、Ｓorg ：読取画像信号、
Ｓus：ボケマスク信号、β：強調度） また、画像信号に対して周波数処理を施す別の方法とし
て、フーリエ変換による周波数解析を用いる方法が知ら
れている。この方法は、画像信号をフーリエ変換して、
画像信号を複数の周波数帯域の信号に分解し、この分解
された信号のうち、所望とする周波数帯域の信号に対し
て強調等の所定の画像処理を施すものである。このフー
リエ変換による周波数解析方法は、画像信号に対して、
Ｓus：ボケマスク信号、β：強調度） また、画像信号に対して周波数処理を施す別の方法とし
て、フーリエ変換による周波数解析を用いる方法が知ら
れている。この方法は、画像信号をフーリエ変換して、
画像信号を複数の周波数帯域の信号に分解し、この分解
された信号のうち、所望とする周波数帯域の信号に対し
て強調等の所定の画像処理を施すものである。このフー
リエ変換による周波数解析方法は、画像信号に対して、
【０００８】
【数１】
【０００９】なる式により、この画像信号中にどの周波
数帯域の成分が含まれているかを解析する方法であり、
振動解析、スペクトル解析等様々の分野における解析に
用いられる方法である。
数帯域の成分が含まれているかを解析する方法であり、
振動解析、スペクトル解析等様々の分野における解析に
用いられる方法である。
【００１０】
【発明が解決しようとする課題】しかしながら、フーリ
エ変換は、画像信号を複数の周波数帯域毎の信号に分解
する際の基底関数が三角関数であるため、信号の振幅等
の性質が信号の部分ごとに変化すると、画像信号のう
ち、本来振動していない部分に振動が発生してしまうと
いう偽振動の問題が発生する。例えば、図７に示すよう
な信号Ｓorg をフーリエ変換し、各周波数帯域毎に逆フ
ーリエ変換した周波数信号で見てみると、原信号Ｓorg
の振動が発生していない部分Ａにおいて、周波数帯域７
の対応する部分Ａ′に偽振動が発生しているのが分か
る。これは、フーリエ変換が原信号Ｓorg を無理に正弦
波に合わせようとしているために、振動していない部分
が他の振動している部分に影響されてしまうことにより
発生するものである。
エ変換は、画像信号を複数の周波数帯域毎の信号に分解
する際の基底関数が三角関数であるため、信号の振幅等
の性質が信号の部分ごとに変化すると、画像信号のう
ち、本来振動していない部分に振動が発生してしまうと
いう偽振動の問題が発生する。例えば、図７に示すよう
な信号Ｓorg をフーリエ変換し、各周波数帯域毎に逆フ
ーリエ変換した周波数信号で見てみると、原信号Ｓorg
の振動が発生していない部分Ａにおいて、周波数帯域７
の対応する部分Ａ′に偽振動が発生しているのが分か
る。これは、フーリエ変換が原信号Ｓorg を無理に正弦
波に合わせようとしているために、振動していない部分
が他の振動している部分に影響されてしまうことにより
発生するものである。
【００１１】このようにフーリエ変換ににおいては、偽
振動が発生するため、フーリエ変換を用いて周波数処理
を行うと、この周波数処理を行った画像信号を再生した
際に、再生画像が視覚的な印象にそぐわず、不自然な強
調等がなされた画像となってしまうという問題がある。
振動が発生するため、フーリエ変換を用いて周波数処理
を行うと、この周波数処理を行った画像信号を再生した
際に、再生画像が視覚的な印象にそぐわず、不自然な強
調等がなされた画像となってしまうという問題がある。
【００１２】また、上述したボケマスク処理は所定の周
波数より高周波の成分をすべて強調してしまう処理であ
るため、このボケマスク処理を用いて周波数強調処理を
行う方法においては、所望とする周波数成分のみならず
画像信号に含まれるノイズをも強調してしまうことがあ
り、観察読影適性に優れた再生画像を得ることができな
かった。
波数より高周波の成分をすべて強調してしまう処理であ
るため、このボケマスク処理を用いて周波数強調処理を
行う方法においては、所望とする周波数成分のみならず
画像信号に含まれるノイズをも強調してしまうことがあ
り、観察読影適性に優れた再生画像を得ることができな
かった。
【００１３】本発明は上記事情に鑑み、観察読影適性に
優れ、視覚的な印象にあった自然な強調等を行うことが
できる放射線画像を表す画像信号に対する画像処理方法
を提供することを目的とするものである。
優れ、視覚的な印象にあった自然な強調等を行うことが
できる放射線画像を表す画像信号に対する画像処理方法
を提供することを目的とするものである。
【００１４】
【課題を解決するための手段】本発明による第１の画像
処理方法は、放射線画像を表す画像信号から該放射線画
像の可視像を再生するに当り、該画像信号に画像処理を
施す画像処理方法において、該画像信号にウェーブレッ
ト変換を施すことにより、該画像信号を複数の周波数帯
域の信号に分解し、該複数の周波数帯域の少なくとも１
つの周波数帯域の信号に対して所定の画像処理を施し、
該画像処理が施された信号および他の信号に逆ウェーブ
レット変換を施すことを特徴とするものである。
処理方法は、放射線画像を表す画像信号から該放射線画
像の可視像を再生するに当り、該画像信号に画像処理を
施す画像処理方法において、該画像信号にウェーブレッ
ト変換を施すことにより、該画像信号を複数の周波数帯
域の信号に分解し、該複数の周波数帯域の少なくとも１
つの周波数帯域の信号に対して所定の画像処理を施し、
該画像処理が施された信号および他の信号に逆ウェーブ
レット変換を施すことを特徴とするものである。
【００１５】また、本発明による第２の画像処理方法
は、上述した本発明による第１の画像処理方法におい
て、前記所定の画像処理が、前記少なくとも１つの周波
数帯域の信号に所定数を乗算する処理であることを特徴
とするものである。
は、上述した本発明による第１の画像処理方法におい
て、前記所定の画像処理が、前記少なくとも１つの周波
数帯域の信号に所定数を乗算する処理であることを特徴
とするものである。
【００１６】ここで、ウェーブレット変換について説明
する。
する。
【００１７】ウェーブレット変換は、周波数解析の方法
として近年開発されたものであり、ステレオのパターン
マッチング、データ圧縮等に応用がなされているもので
ある（OLIVIER RIOUL and MARTIN VETTERLI;Wavelets a
nd Signal Processing,IEEESP MAGAZINE,P.14-38,OCTOB
ER 1991、Stephane Mallat;Zero-Crossings of a Wavel
et Transform,IEEE TRANSACTIONS ON INFORMATION THEO
RY,VOL.37,NO.4,P.1019-1033,JULY 1991 ）。
として近年開発されたものであり、ステレオのパターン
マッチング、データ圧縮等に応用がなされているもので
ある（OLIVIER RIOUL and MARTIN VETTERLI;Wavelets a
nd Signal Processing,IEEESP MAGAZINE,P.14-38,OCTOB
ER 1991、Stephane Mallat;Zero-Crossings of a Wavel
et Transform,IEEE TRANSACTIONS ON INFORMATION THEO
RY,VOL.37,NO.4,P.1019-1033,JULY 1991 ）。
【００１８】このウェーブレット変換は、図８に示すよ
うな関数ｈを基底関数として、
うな関数ｈを基底関数として、
【００１９】
【数２】
【００２０】なる式において信号を複数の周波数帯域の
信号に分解するため、フーリエ変換のような偽振動の問
題が発生しない。すなわち、関数ｈの周期および縮率を
変化させ、原信号を移動させることによりフィルタリン
グ処理を行えば、細かな周波数から粗い周波数までの所
望とする周波数に適合した信号を作成することができ
る。例えば、図９に示すように、前述した図７に示した
信号と同一の信号Sorgをウェーブレット変換し、各周波
数帯域毎に逆ウェーブレット変換した信号で見てみる
と、原信号Ｓorg の振動と対応した周波数帯域の信号を
得ることができる。すなわち、原信号Ｓorg の部分Ｂと
対応する周波数帯域Ｗ７の部分Ｂ′には、原信号と同様
に振動が発生していないものとなる。
信号に分解するため、フーリエ変換のような偽振動の問
題が発生しない。すなわち、関数ｈの周期および縮率を
変化させ、原信号を移動させることによりフィルタリン
グ処理を行えば、細かな周波数から粗い周波数までの所
望とする周波数に適合した信号を作成することができ
る。例えば、図９に示すように、前述した図７に示した
信号と同一の信号Sorgをウェーブレット変換し、各周波
数帯域毎に逆ウェーブレット変換した信号で見てみる
と、原信号Ｓorg の振動と対応した周波数帯域の信号を
得ることができる。すなわち、原信号Ｓorg の部分Ｂと
対応する周波数帯域Ｗ７の部分Ｂ′には、原信号と同様
に振動が発生していないものとなる。
【００２１】
【作用】本発明による画像処理方法は、放射線画像を表
す画像信号に対してウェーブレット変換を施して、画像
信号を複数の周波数帯域の信号に分解するようにしたた
め、この信号には、前述した図７に示すフーリエ変換の
ように偽振動が発生していないものとなる。したがっ
て、この複数の周波数帯域の信号のうち、少なくとも１
つの周波数帯域の信号に対してボケマスク処理、階調処
理、強調処理等の所定の画像処理を施し、この画像処理
が施された信号および他の信号に逆ウェーブレット変換
を施せば、ノイズを強調することなく、所望とする周波
数帯域の画像に対して視覚的な印象にあった自然な強調
等を行うことができる。
す画像信号に対してウェーブレット変換を施して、画像
信号を複数の周波数帯域の信号に分解するようにしたた
め、この信号には、前述した図７に示すフーリエ変換の
ように偽振動が発生していないものとなる。したがっ
て、この複数の周波数帯域の信号のうち、少なくとも１
つの周波数帯域の信号に対してボケマスク処理、階調処
理、強調処理等の所定の画像処理を施し、この画像処理
が施された信号および他の信号に逆ウェーブレット変換
を施せば、ノイズを強調することなく、所望とする周波
数帯域の画像に対して視覚的な印象にあった自然な強調
等を行うことができる。
【００２２】なお、ウェーブレット変換した画像信号に
対して画像処理を施さずに逆ウェーブレット変換をした
場合、原画像が完全に復元できることが数学的に保証さ
れているため、本発明の画像処理方法により、元の放射
線画像の情報が欠落することはない。
対して画像処理を施さずに逆ウェーブレット変換をした
場合、原画像が完全に復元できることが数学的に保証さ
れているため、本発明の画像処理方法により、元の放射
線画像の情報が欠落することはない。
【００２３】
【実施例】以下、図面を参照して本発明の実施例につい
て説明する。
て説明する。
【００２４】図１は、本発明による画像処理方法の基本
的概念を表す図である。図１に示すように、本発明によ
る画像処理方法は、放射線画像を表す画像信号１に対し
て、ウェーブレット変換２を施して、この画像信号１を
複数の周波数帯域の信号３に分解する。次いで、この信
号３のうち、少なくとも１つの周波数帯域の信号３に対
して所定の画像処理４を施し、この画像処理４が施され
た信号と他の信号に逆ウェーブレット変換５を施すよう
にしたものである。
的概念を表す図である。図１に示すように、本発明によ
る画像処理方法は、放射線画像を表す画像信号１に対し
て、ウェーブレット変換２を施して、この画像信号１を
複数の周波数帯域の信号３に分解する。次いで、この信
号３のうち、少なくとも１つの周波数帯域の信号３に対
して所定の画像処理４を施し、この画像処理４が施され
た信号と他の信号に逆ウェーブレット変換５を施すよう
にしたものである。
【００２５】以下、本発明による画像処理方法の実施例
について詳細に説明する。図２は、放射線画像読取装置
の一例を表す図である。
について詳細に説明する。図２は、放射線画像読取装置
の一例を表す図である。
【００２６】図示しない放射線撮影装置で撮影が行わ
れ、放射線画像が記録された蓄積性蛍光体シート14が読
取装置20の所定位置にセットされる。
れ、放射線画像が記録された蓄積性蛍光体シート14が読
取装置20の所定位置にセットされる。
【００２７】蓄積性蛍光体シート14が読取装置20の所定
位置にセットされると、このシート14はモータ21により
駆動されるエンドレスベルト22により、矢印Ｙ方向に搬
送（副走査）される。一方、レーザー光源23から発せら
れた光ビーム24はモータ25により駆動され矢印方向に高
速回転する回転多面鏡26によって反射偏向され、ｆθレ
ンズ等の集束レンズ27を通過した後、ミラー28により光
路を変えて前記シート14に入射し副走査の方向（矢印Ｙ
方向）と略垂直な矢印Ｘ方向に主走査する。シート14の
励起光24が照射された箇所からは、蓄積記録されている
放射線画像情報に応じた光量の輝尽発光光29が発散さ
れ、この輝尽発光光29は光ガイド30によって導かれ、フ
ォトマルチプライヤ（光電子増倍管）31によって光電的
に検出される。上記光ガイド30はアクリル板等の導光性
材料を成形して作られたものであり、直線状をなす入射
端面30a が蓄積性蛍光体シート14上の主走査線に沿って
延びるように配され、円環状に形成された射出端面30b
にフォトマルチプライヤ31の受光面が結合されている。
入射端面30a から光ガイド30内に入射した輝尽発光光29
は、該光ガイド30の内部を全反射を繰り返して進み、射
出端面30b から射出してフォトマルチプライヤ31に受光
され、放射線画像を表わす輝尽発光光29がフォトマルチ
プライヤ31によって電気信号に変換される。
位置にセットされると、このシート14はモータ21により
駆動されるエンドレスベルト22により、矢印Ｙ方向に搬
送（副走査）される。一方、レーザー光源23から発せら
れた光ビーム24はモータ25により駆動され矢印方向に高
速回転する回転多面鏡26によって反射偏向され、ｆθレ
ンズ等の集束レンズ27を通過した後、ミラー28により光
路を変えて前記シート14に入射し副走査の方向（矢印Ｙ
方向）と略垂直な矢印Ｘ方向に主走査する。シート14の
励起光24が照射された箇所からは、蓄積記録されている
放射線画像情報に応じた光量の輝尽発光光29が発散さ
れ、この輝尽発光光29は光ガイド30によって導かれ、フ
ォトマルチプライヤ（光電子増倍管）31によって光電的
に検出される。上記光ガイド30はアクリル板等の導光性
材料を成形して作られたものであり、直線状をなす入射
端面30a が蓄積性蛍光体シート14上の主走査線に沿って
延びるように配され、円環状に形成された射出端面30b
にフォトマルチプライヤ31の受光面が結合されている。
入射端面30a から光ガイド30内に入射した輝尽発光光29
は、該光ガイド30の内部を全反射を繰り返して進み、射
出端面30b から射出してフォトマルチプライヤ31に受光
され、放射線画像を表わす輝尽発光光29がフォトマルチ
プライヤ31によって電気信号に変換される。
【００２８】フォトマルチプライヤ31から出力されたア
ナログ出力信号Ｓ０は対数増幅器32で対数的に増幅さ
れ、Ａ／Ｄ変換器33でディジタル化され、これにより画
像信号Ｓorg が得られ、画像処理装置40に入力される。
この画像処理装置40は、可視画像を再生表示するための
ＣＲＴディスプレイ41、ＣＰＵ，内部メモリ，インター
フェイス等が内蔵された本体部42、フロッピィディスク
が装填され駆動されるフロッピィディスクドライブ部4
3、およびこのＸ線画像読取装置に必要な情報を入力す
るためのキーボード44から構成されている。
ナログ出力信号Ｓ０は対数増幅器32で対数的に増幅さ
れ、Ａ／Ｄ変換器33でディジタル化され、これにより画
像信号Ｓorg が得られ、画像処理装置40に入力される。
この画像処理装置40は、可視画像を再生表示するための
ＣＲＴディスプレイ41、ＣＰＵ，内部メモリ，インター
フェイス等が内蔵された本体部42、フロッピィディスク
が装填され駆動されるフロッピィディスクドライブ部4
3、およびこのＸ線画像読取装置に必要な情報を入力す
るためのキーボード44から構成されている。
【００２９】この画像処理装置40において上述した画像
処理が行われる。
処理が行われる。
【００３０】図３は、画像信号Ｓorg に対するウェーブ
レット変換の詳細を表す図である。
レット変換の詳細を表す図である。
【００３１】なお、本実施例においては、ウェーブレッ
ト変換の各係数が直交する直交ウェーブレット変換を行
うものである。
ト変換の各係数が直交する直交ウェーブレット変換を行
うものである。
【００３２】図３に示すように、画像信号Ｓorg の主走
査方向に基本ウェーブレット関数より求められる関数ｇ
0 と関数ｈ0 とによりフィルタリング処理を行う。すな
わち、このような関数ｇ0 ，ｈ0 による主走査方向に並
ぶ画素の一列毎のフィルタリング処理を副走査方向に一
画素ずつズラしながら行い、画像信号Ｓorg の主走査方
向のウェーブレット変換係数信号Ｗg0，Ｗh0を求めるも
のである。
査方向に基本ウェーブレット関数より求められる関数ｇ
0 と関数ｈ0 とによりフィルタリング処理を行う。すな
わち、このような関数ｇ0 ，ｈ0 による主走査方向に並
ぶ画素の一列毎のフィルタリング処理を副走査方向に一
画素ずつズラしながら行い、画像信号Ｓorg の主走査方
向のウェーブレット変換係数信号Ｗg0，Ｗh0を求めるも
のである。
【００３３】ここで、関数ｇ0 ，ｈ0 は基本ウェーブレ
ット関数より一意に求められるものであり、例えば、図
８に示す基本ウェーブレット関数の場合、以下の表１に
示すものとなる。また、関数ｇ0 ，ｈ0 を図に示すと、
関数ｇ0 は図４(a) に、関数ｈ0 は図５(a) にそれぞれ
示すものとなる。
ット関数より一意に求められるものであり、例えば、図
８に示す基本ウェーブレット関数の場合、以下の表１に
示すものとなる。また、関数ｇ0 ，ｈ0 を図に示すと、
関数ｇ0 は図４(a) に、関数ｈ0 は図５(a) にそれぞれ
示すものとなる。
【００３４】
【表１】
【００３５】このようにして、ウェーブレット変換係数
信号Ｗg0、Ｗh0が求められると、ウェーブレット変換係
数信号Ｗg0、Ｗh0それぞれの副走査方向に関数ｇ0 ，ｈ
0 によりフィルタリング処理を行い、ウェーブレット変
換係数信号ＷＷ0 ，ＷＶ0 ，ＶＷ0 およびＶＶ0 を得
る。次いで、ウェーブレット変換係数信号ＶＶ0 の主走
査方向に関数ｇ1 ，ｈ1 によりフィルタリング処理を行
う。
信号Ｗg0、Ｗh0が求められると、ウェーブレット変換係
数信号Ｗg0、Ｗh0それぞれの副走査方向に関数ｇ0 ，ｈ
0 によりフィルタリング処理を行い、ウェーブレット変
換係数信号ＷＷ0 ，ＷＶ0 ，ＶＷ0 およびＶＶ0 を得
る。次いで、ウェーブレット変換係数信号ＶＶ0 の主走
査方向に関数ｇ1 ，ｈ1 によりフィルタリング処理を行
う。
【００３６】すなわち、前述したｇ0 ，ｈ0 の場合と同
様に、関数ｇ1 ，ｈ1 により主走査方向に並ぶ画素の一
列毎のフィルタリング処理を副走査方向に一画素づつズ
ラながら行い、ウェーブレット変換係数信号ＶＶ0 の主
走査方向のウェーブレット変換係数信号Ｗg1およびＷh1
を求めるものである。
様に、関数ｇ1 ，ｈ1 により主走査方向に並ぶ画素の一
列毎のフィルタリング処理を副走査方向に一画素づつズ
ラながら行い、ウェーブレット変換係数信号ＶＶ0 の主
走査方向のウェーブレット変換係数信号Ｗg1およびＷh1
を求めるものである。
【００３７】ここで、関数ｇ1 ，ｈ1 は、関数ｇ0 ，ｈ
0 の各サンプリング点の間に１個の０を挿入することに
より作成したものである。また、関数ｇ1 ，ｈ1 を図に
示すと、関数ｇ1 は図４(b) に、関数ｈ1 は図５(b) に
それぞれ示すものとなる。
0 の各サンプリング点の間に１個の０を挿入することに
より作成したものである。また、関数ｇ1 ，ｈ1 を図に
示すと、関数ｇ1 は図４(b) に、関数ｈ1 は図５(b) に
それぞれ示すものとなる。
【００３８】このようにして、ウェーブレット変換係数
信号Ｗg1、Ｗh1が求められると、ウェーブレット変換係
数信号Ｗg1、Ｗh1それぞれの副走査方向に関数ｇ1 ，ｈ
1 によりフィルタリング処理を行い、ウェーブレット変
換係数信号ＷＷ1 ，ＷＶ1 ，ＶＷ1 およびＶＶ1 を得
る。
信号Ｗg1、Ｗh1が求められると、ウェーブレット変換係
数信号Ｗg1、Ｗh1それぞれの副走査方向に関数ｇ1 ，ｈ
1 によりフィルタリング処理を行い、ウェーブレット変
換係数信号ＷＷ1 ，ＷＶ1 ，ＶＷ1 およびＶＶ1 を得
る。
【００３９】以下、上述したのと同様にして、ウェーブ
レット変換係数信号ＶＶ1 の主走査方向に関数ｇ2 ，ｈ
2 によりフィルタリング処理を行い、さらに得られたウ
ェーブレット変換係数信号に対し、副走査方向に関数ｇ
2 ，ｈ2 によりフィルタリング処理を行い、ウェーブレ
ット変換係数信号ＷＷ2 ，ＷＶ2 ，ＶＷ2 ，ＶＶ2 を得
る。
レット変換係数信号ＶＶ1 の主走査方向に関数ｇ2 ，ｈ
2 によりフィルタリング処理を行い、さらに得られたウ
ェーブレット変換係数信号に対し、副走査方向に関数ｇ
2 ，ｈ2 によりフィルタリング処理を行い、ウェーブレ
ット変換係数信号ＷＷ2 ，ＷＶ2 ，ＶＷ2 ，ＶＶ2 を得
る。
【００４０】ここで、関数ｇ2 ，ｈ2 は関数ｇ0 ，ｈ0
の各サンプリング点の間に２個の０を入れることにより
作成したものであり、関数ｇ2 は図４(c) に、関数ｈ2
は図５(c) にそれぞれ示すものとなる。
の各サンプリング点の間に２個の０を入れることにより
作成したものであり、関数ｇ2 は図４(c) に、関数ｈ2
は図５(c) にそれぞれ示すものとなる。
【００４１】このようなウェーブレット変換をＮ回繰り
返すことによりウェーブレット変換係数信号ＷＷ0 〜Ｗ
ＷN ，ＷＶ0 〜ＷＶN ，ＶＷ0 〜ＶＷN ，およびＶＶN
を得る。ここで、Ｎ回目のウェーブレット変換の際に用
いられる関数ｇN ，ｈN は、関数ｇ0 ，ｈ0 の各サンプ
リング点の間に２N-1 個の０を入れることにより作成さ
れる。したがって、関数ｇN ，ｈN はＮが大きいほど周
期の大きな関数となる。
返すことによりウェーブレット変換係数信号ＷＷ0 〜Ｗ
ＷN ，ＷＶ0 〜ＷＶN ，ＶＷ0 〜ＶＷN ，およびＶＶN
を得る。ここで、Ｎ回目のウェーブレット変換の際に用
いられる関数ｇN ，ｈN は、関数ｇ0 ，ｈ0 の各サンプ
リング点の間に２N-1 個の０を入れることにより作成さ
れる。したがって、関数ｇN ，ｈN はＮが大きいほど周
期の大きな関数となる。
【００４２】ここで、ウェーブレット変換係数信号ＷＷ
i （ｉ＝０〜Ｎ、以下同様）は、画像信号Ｓorg の主副
両方向の周波数の変化を表すものであり、ｉが大きいほ
ど低周波信号となる。またウェーブレット変換係数信号
ＷＶi は画像信号Ｓorg の主走査方向の周波数の変化を
表すものであり、ｉが大きいほど低周波信号となる。さ
らにウェーブレット変換係数信号ＶＷi は画像信号Ｓor
g の副走査方向の周波数の変化を表すものであり、ｉが
大きいほど低周波信号となる。
i （ｉ＝０〜Ｎ、以下同様）は、画像信号Ｓorg の主副
両方向の周波数の変化を表すものであり、ｉが大きいほ
ど低周波信号となる。またウェーブレット変換係数信号
ＷＶi は画像信号Ｓorg の主走査方向の周波数の変化を
表すものであり、ｉが大きいほど低周波信号となる。さ
らにウェーブレット変換係数信号ＶＷi は画像信号Ｓor
g の副走査方向の周波数の変化を表すものであり、ｉが
大きいほど低周波信号となる。
【００４３】次いで、このようにして得られたウェーブ
レット変換係数信号に対して予め設定されたテーブルβ
WW(i) ，βWV(i) ，βVW(i) を乗算することにより所望
とする周波数帯域の強調を行う。すなわち、ｉ番目のウ
ェーブレット変換係数信号ＷＷi ，ＷＶi ，ＶＷi につ
いて ＷＷi ′＝βWW(i) ×ＷＷi …(4) ＷＶi ′＝βWV(i) ×ＷＶi …(5) ＶＷi ′＝βVW(i) ×ＶＷi …(6) なる演算を行うものである。
レット変換係数信号に対して予め設定されたテーブルβ
WW(i) ，βWV(i) ，βVW(i) を乗算することにより所望
とする周波数帯域の強調を行う。すなわち、ｉ番目のウ
ェーブレット変換係数信号ＷＷi ，ＷＶi ，ＶＷi につ
いて ＷＷi ′＝βWW(i) ×ＷＷi …(4) ＷＶi ′＝βWV(i) ×ＷＶi …(5) ＶＷi ′＝βVW(i) ×ＶＷi …(6) なる演算を行うものである。
【００４４】ここで、放射線画像の高周波成分を強調し
たいときにはβWW(i) ，βWV(i) ，βVW(i) について、
ｉが小さいほど大きい値を設定すればよく、また、低周
波成分を強調したいときには、ｉが大きいほど大きい値
を設定すればよい。さらに、放射線画像の主走査方向成
分のみ強調したいときには、βWV(i) についてβWW(i)
，βVW(i) よりも大きい値を設定すればよい。
たいときにはβWW(i) ，βWV(i) ，βVW(i) について、
ｉが小さいほど大きい値を設定すればよく、また、低周
波成分を強調したいときには、ｉが大きいほど大きい値
を設定すればよい。さらに、放射線画像の主走査方向成
分のみ強調したいときには、βWV(i) についてβWW(i)
，βVW(i) よりも大きい値を設定すればよい。
【００４５】以上のようにしてβWW(i) ，βWV(i) およ
びβVW(i) が乗算されたウェーブレット変換係数信号Ｗ
Ｗi ′，ＷＶi ′およびＶＷi ′が求められると、これ
らの係数信号ＷＷi ′，ＷＶi ′，ＶＷi ′およびＶＶ
i ′を用いて、逆ウェーブレット変換を行う。
びβVW(i) が乗算されたウェーブレット変換係数信号Ｗ
Ｗi ′，ＷＶi ′およびＶＷi ′が求められると、これ
らの係数信号ＷＷi ′，ＷＶi ′，ＶＷi ′およびＶＶ
i ′を用いて、逆ウェーブレット変換を行う。
【００４６】図６は、逆ウェーブレット変換の詳細を表
す図である。
す図である。
【００４７】図６に示すように、まずウェーブレット変
換係数信号ＶＶN を副走査方向に前述した関数ｈN によ
り、ウェーブレット変換係数信号ＶＷN を副走査方向に
前述した関数ｇN によりフィルタリング処理を行う。す
なわち、関数ｇN ，ｈN によるウェーブレット変換係数
信号ＶＶN ，ＶＷN ′の副走査方向に並ぶ一列の画素毎
のフィルタリング処理を主走査方向に一画素ずつズラし
ながら行い、ウェーブレット変換係数信号ＶＶN ，ＶＷ
N ′の逆ウェーブレット変換係数信号を得、これを加算
することにより逆ウェーブレット変換係数信号ＷhN′を
得る。
換係数信号ＶＶN を副走査方向に前述した関数ｈN によ
り、ウェーブレット変換係数信号ＶＷN を副走査方向に
前述した関数ｇN によりフィルタリング処理を行う。す
なわち、関数ｇN ，ｈN によるウェーブレット変換係数
信号ＶＶN ，ＶＷN ′の副走査方向に並ぶ一列の画素毎
のフィルタリング処理を主走査方向に一画素ずつズラし
ながら行い、ウェーブレット変換係数信号ＶＶN ，ＶＷ
N ′の逆ウェーブレット変換係数信号を得、これを加算
することにより逆ウェーブレット変換係数信号ＷhN′を
得る。
【００４８】一方、これと並列して、ウェーブレット変
換係数信号ＷＶN ′を副走査方向に関数ｈN により、ウ
ェーブレット変換係数信号ＷＷN ′を副走査方向に関数
ｇNによりフィルタリング処理を行い、ウェーブレット
変換係数信号ＷＶN ′，ＷＷN ′の逆ウェーブレット変
換係数信号を得、これを加算することにより逆ウェーブ
レット変換係数信号ＷgN′を得る。
換係数信号ＷＶN ′を副走査方向に関数ｈN により、ウ
ェーブレット変換係数信号ＷＷN ′を副走査方向に関数
ｇNによりフィルタリング処理を行い、ウェーブレット
変換係数信号ＷＶN ′，ＷＷN ′の逆ウェーブレット変
換係数信号を得、これを加算することにより逆ウェーブ
レット変換係数信号ＷgN′を得る。
【００４９】次いで、逆ウェーブレット変換係数信号Ｗ
hN′を主走査方向に関数ｈN ′により、逆ウェーブレッ
ト変換係数信号ＷgN′を主走査方向に関数ｇN によりフ
ィルタリング処理し、ウェーブレット変換係数信号Ｗh
N′，ＷgN′の逆ウェーブレット変換係数信号を得、こ
れを加算することにより逆ウェーブレット変換係数信号
ＶＶN-1 ′を得る。
hN′を主走査方向に関数ｈN ′により、逆ウェーブレッ
ト変換係数信号ＷgN′を主走査方向に関数ｇN によりフ
ィルタリング処理し、ウェーブレット変換係数信号Ｗh
N′，ＷgN′の逆ウェーブレット変換係数信号を得、こ
れを加算することにより逆ウェーブレット変換係数信号
ＶＶN-1 ′を得る。
【００５０】次いでこの逆ウェーブレット変換係数信号
ＶＶN-1 ′を副走査方向に前述した関数ｈN-1 により、
ウェーブレット変換係数信号ＶＷN-1 を副走査方向に前
述した関数ｇN-1 によりフィルタリング処理を行う。す
なわち、関数ｇN-1 ，ｈN-1によるウェーブレット変換
係数信号ＶＶN-1 ，ＶＷN-1 ′の副走査方向に並ぶ一列
の画素毎のフィルタリング処理を主走査方向に一画素ず
つズラしながら行い、ウェーブレット変換係数信号ＶＶ
N-1 ，ＶＷN-1 ′の逆ウェーブレット変換係数信号を
得、これを加算することにより逆ウェーブレット変換係
数信号ＷhN-1′を得る。
ＶＶN-1 ′を副走査方向に前述した関数ｈN-1 により、
ウェーブレット変換係数信号ＶＷN-1 を副走査方向に前
述した関数ｇN-1 によりフィルタリング処理を行う。す
なわち、関数ｇN-1 ，ｈN-1によるウェーブレット変換
係数信号ＶＶN-1 ，ＶＷN-1 ′の副走査方向に並ぶ一列
の画素毎のフィルタリング処理を主走査方向に一画素ず
つズラしながら行い、ウェーブレット変換係数信号ＶＶ
N-1 ，ＶＷN-1 ′の逆ウェーブレット変換係数信号を
得、これを加算することにより逆ウェーブレット変換係
数信号ＷhN-1′を得る。
【００５１】一方、これと並列して、ウェーブレット変
換係数信号ＷＶN-1 ′を副走査方向に関数ｈN-1 によ
り、ウェーブレット変換係数信号ＷＷN-1 ′を副走査方
向に関数ｇN-1 によりフィルタリング処理を行い、ウェ
ーブレット変換係数信号ＷＶN-1 ′，ＷＷN-1 ′の逆ウ
ェーブレット変換係数信号を得、これを加算することに
より逆ウェーブレット変換係数信号ＷgN-1′を得る。
換係数信号ＷＶN-1 ′を副走査方向に関数ｈN-1 によ
り、ウェーブレット変換係数信号ＷＷN-1 ′を副走査方
向に関数ｇN-1 によりフィルタリング処理を行い、ウェ
ーブレット変換係数信号ＷＶN-1 ′，ＷＷN-1 ′の逆ウ
ェーブレット変換係数信号を得、これを加算することに
より逆ウェーブレット変換係数信号ＷgN-1′を得る。
【００５２】次いで、逆ウェーブレット変換係数信号Ｗ
hN-1′を主走査方向に関数ｈN-1 ′により、逆ウェーブ
レット変換係数信号ＷgN-1′を主走査方向に関数ｇN-1
によりフィルタリング処理し、ウェーブレット変換係数
信号ＷhN-1′，ＷgN-1′の逆ウェーブレット変換係数信
号を得、これを加算することにより逆ウェーブレット変
換係数信号ＶＶN-2 ′を得る。
hN-1′を主走査方向に関数ｈN-1 ′により、逆ウェーブ
レット変換係数信号ＷgN-1′を主走査方向に関数ｇN-1
によりフィルタリング処理し、ウェーブレット変換係数
信号ＷhN-1′，ＷgN-1′の逆ウェーブレット変換係数信
号を得、これを加算することにより逆ウェーブレット変
換係数信号ＶＶN-2 ′を得る。
【００５３】以下、順次逆ウェーブレット変換係数信号
ＶＶi ′（ｉ＝−１〜Ｎ）を作成し、最終的に逆ウェー
ブレット変換係数信号ＶＶ-1′を得る。この最終的な逆
ウェーブレット変換係数信号ＶＶ-1′が画像信号Ｓorg
の所望とする周波数帯域が強調された画像信号となる。
ＶＶi ′（ｉ＝−１〜Ｎ）を作成し、最終的に逆ウェー
ブレット変換係数信号ＶＶ-1′を得る。この最終的な逆
ウェーブレット変換係数信号ＶＶ-1′が画像信号Ｓorg
の所望とする周波数帯域が強調された画像信号となる。
【００５４】このようにして得られたウェーブレット変
換係数信号ＶＶ-1′は図示しない画像再生装置に送られ
て、放射線画像の再生に供せられる。
換係数信号ＶＶ-1′は図示しない画像再生装置に送られ
て、放射線画像の再生に供せられる。
【００５５】この再生装置は、ＣＲＴ等のディスプレイ
手段でもよいし、感光フイルムに光走査記録を行う記録
装置であってもよいし、あるいはそのために画像信号を
一旦光ディスク、磁気ディスク等の画像ファイルに記憶
させる装置に置き換えられてもよい。
手段でもよいし、感光フイルムに光走査記録を行う記録
装置であってもよいし、あるいはそのために画像信号を
一旦光ディスク、磁気ディスク等の画像ファイルに記憶
させる装置に置き換えられてもよい。
【００５６】このようにして、ウェーブレット変換を行
うことにより得られた信号に強調処理を行うことによ
り、フーリエ変換により得られた信号に強調処理を行っ
た場合よりも、より視覚的な印象に合った画像を得るこ
とができる。例えば、放射線画像に急に濃度が変化する
エッジ（心臓の影、骨等）が存在する場合に、この放射
線画像を表す画像信号をフーリエ変換して、エッジ部分
を強調するとエッジの所で偽輪郭が目立ってしまう。こ
れは、無理にエッジ部分の信号を三角関数に合わせよう
としているため、画像信号のうち振動が発生していない
部分にまで振動が発生してしまうためである。これに対
し、放射線画像を表す画像信号をウェーブレット変換す
ると、偽振動のない信号が得られるため、エッジ部分は
比較的自然に強調され、偽輪郭のない自然な画像を得る
ことができる。
うことにより得られた信号に強調処理を行うことによ
り、フーリエ変換により得られた信号に強調処理を行っ
た場合よりも、より視覚的な印象に合った画像を得るこ
とができる。例えば、放射線画像に急に濃度が変化する
エッジ（心臓の影、骨等）が存在する場合に、この放射
線画像を表す画像信号をフーリエ変換して、エッジ部分
を強調するとエッジの所で偽輪郭が目立ってしまう。こ
れは、無理にエッジ部分の信号を三角関数に合わせよう
としているため、画像信号のうち振動が発生していない
部分にまで振動が発生してしまうためである。これに対
し、放射線画像を表す画像信号をウェーブレット変換す
ると、偽振動のない信号が得られるため、エッジ部分は
比較的自然に強調され、偽輪郭のない自然な画像を得る
ことができる。
【００５７】また、所望とする周波数帯域の信号のみを
強調することができるため、前述した式(1) に示すボケ
マスク処理のようにノイズ成分となる高周波成分を強調
することなく、画像の強調処理を行うことができる。
強調することができるため、前述した式(1) に示すボケ
マスク処理のようにノイズ成分となる高周波成分を強調
することなく、画像の強調処理を行うことができる。
【００５８】なお、上述した実施例においては、ウェー
ブレット変換係数信号にテーブルβWWi ，βWVi ，β
VWi を乗算して、画像Ｓorg のうち所望とする周波数帯
域の信号に対して強調処理を行っているが、これに限定
されるものではなく、ウェーブレット変換係数信号ＷＷ
i ，ＷＶi およびＶＷi に対して、異なるボケマスク処
理を施すようにしてもよい。例えば、放射線画像のう
ち、主要被写体を表現している周波数帯域の信号にのみ
前述した式(1) に示すボケマスク処理を施すことによ
り、主要被写体部分の画像強調を行うことができる。ま
た、ウェーブレット変換係数信号ＷＷi ，ＷＶi および
ＶＷi に対して、階調処理を施すようにしてもよい。
ブレット変換係数信号にテーブルβWWi ，βWVi ，β
VWi を乗算して、画像Ｓorg のうち所望とする周波数帯
域の信号に対して強調処理を行っているが、これに限定
されるものではなく、ウェーブレット変換係数信号ＷＷ
i ，ＷＶi およびＶＷi に対して、異なるボケマスク処
理を施すようにしてもよい。例えば、放射線画像のう
ち、主要被写体を表現している周波数帯域の信号にのみ
前述した式(1) に示すボケマスク処理を施すことによ
り、主要被写体部分の画像強調を行うことができる。ま
た、ウェーブレット変換係数信号ＷＷi ，ＷＶi および
ＶＷi に対して、階調処理を施すようにしてもよい。
【００５９】また、上述した実施例においては、基本ウ
ェーブレット関数として図８に示した関数を用いている
が、これに限定されるものではなく、ウェーブレット変
換に適した関数であれば、いかなる関数を用いてもよ
い。
ェーブレット関数として図８に示した関数を用いている
が、これに限定されるものではなく、ウェーブレット変
換に適した関数であれば、いかなる関数を用いてもよ
い。
【００６０】さらに、上述した実施例において、ウェー
ブレット変換として直交ウェーブレット変換を用いた例
を示したが、これに限定されるものではなく、非直交の
ウェーブレット変換によって周波数帯域に分割された係
数信号を用いてもよく、２次元基本ウェーブレットを拡
大、縮小、移動または回転させてコンボリューションを
行うウェーブレット変換を用いてもよい。
ブレット変換として直交ウェーブレット変換を用いた例
を示したが、これに限定されるものではなく、非直交の
ウェーブレット変換によって周波数帯域に分割された係
数信号を用いてもよく、２次元基本ウェーブレットを拡
大、縮小、移動または回転させてコンボリューションを
行うウェーブレット変換を用いてもよい。
【００６１】
【発明の効果】以上、詳細に説明したように、本発明に
よる画像処理方法は、画像信号にウェーブレット変換を
施すことにより得られた複数の周波数帯域の信号のうち
の、少なくとも１つの周波数帯域の信号に強調等の画像
処理を施すようにしたため、この画像処理が施された画
像信号を逆ウェーブレット変換することにより得られる
画像信号を再生した際に、ノイズを強調することなく、
視覚的な印象にあった自然な再生画像を得ることができ
る。
よる画像処理方法は、画像信号にウェーブレット変換を
施すことにより得られた複数の周波数帯域の信号のうち
の、少なくとも１つの周波数帯域の信号に強調等の画像
処理を施すようにしたため、この画像処理が施された画
像信号を逆ウェーブレット変換することにより得られる
画像信号を再生した際に、ノイズを強調することなく、
視覚的な印象にあった自然な再生画像を得ることができ
る。
【図１】本発明による画像処理方法の基本的概念を表す
図
図
【図２】放射線画像読取装置を表す図
【図３】ウェーブレット変換の詳細を表す図
【図４】関数ｇ0 を表す図
【図５】関数ｈ0 を表す図
【図６】逆ウェーブレット変換の詳細を表す図
【図７】フーリエ変換を説明するための図
【図８】ウェーブレット変換に用いられる基本ウェーブ
レット関数を表す図
レット関数を表す図
【図９】ウェーブレット変換を説明するための図
14 蓄積性蛍光体シート 23 レーザ光源 24 レーザ光 29 輝尽発光光 40 画像処理装置 44 キーボード
Claims (2)
- 【請求項１】 放射線画像を表す画像信号から該放射線
画像の可視像を再生するに当り、該画像信号に画像処理
を施す画像処理方法において、 該画像信号にウェーブレット変換を施すことにより、該
画像信号を複数の周波数帯域の信号に分解し、 該複数の周波数帯域の少なくとも１つの周波数帯域の信
号に対して所定の画像処理を施し、 該画像処理が施された信号および他の信号に逆ウェーブ
レット変換を施すことを特徴とする画像処理方法。 - 【請求項２】 前記所定の画像処理が、前記少なくとも
１つの周波数帯域の信号に所定数を乗算する処理である
ことを特徴とする請求項１記載の画像処理方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP5057590A JPH06274614A (ja) | 1993-03-18 | 1993-03-18 | 画像処理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP5057590A JPH06274614A (ja) | 1993-03-18 | 1993-03-18 | 画像処理方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH06274614A true JPH06274614A (ja) | 1994-09-30 |
Family
ID=13060067
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP5057590A Pending JPH06274614A (ja) | 1993-03-18 | 1993-03-18 | 画像処理方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH06274614A (ja) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001292308A (ja) * | 2000-02-02 | 2001-10-19 | Fuji Photo Film Co Ltd | 画像保存転送方法および装置、並びに画像処理方法、画像処理システム、および画像処理装置 |
| WO2002045020A1 (fr) * | 2000-11-30 | 2002-06-06 | Canon Kabushiki Kaisha | Dispositif et procede de traitement d'images, support de stockage et programme associe |
| WO2002045019A1 (fr) * | 2000-11-30 | 2002-06-06 | Canon Kabushiki Kaisha | Dispositif de traitement d'image, procede de traitement d'image, support d'enregistrement et programme correspondant |
| JP2005005845A (ja) * | 2003-06-10 | 2005-01-06 | Canon Inc | 画像処理装置 |
| US6868187B2 (en) | 2000-05-12 | 2005-03-15 | Fuji Photo Film Co., Ltd. | Image transformation method and apparatus, and storage medium |
| JP2005109867A (ja) * | 2003-09-30 | 2005-04-21 | Konica Minolta Medical & Graphic Inc | 画像処理方法および画像処理装置ならびに画像処理プログラム |
| US7123773B2 (en) | 2000-02-09 | 2006-10-17 | Fuji Photo Film Co., Ltd. | Image coding-decoding method, image coding-decoding system, image coder, image decoder, and storage medium |
| US7423781B2 (en) | 2002-03-20 | 2008-09-09 | Ricoh Company, Ltd. | Image processor and image processing method for image enhancement using edge detection |
| WO2009110260A1 (ja) * | 2008-03-06 | 2009-09-11 | コニカミノルタエムジー株式会社 | 画像処理装置 |
-
1993
- 1993-03-18 JP JP5057590A patent/JPH06274614A/ja active Pending
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001292308A (ja) * | 2000-02-02 | 2001-10-19 | Fuji Photo Film Co Ltd | 画像保存転送方法および装置、並びに画像処理方法、画像処理システム、および画像処理装置 |
| US7123773B2 (en) | 2000-02-09 | 2006-10-17 | Fuji Photo Film Co., Ltd. | Image coding-decoding method, image coding-decoding system, image coder, image decoder, and storage medium |
| US6868187B2 (en) | 2000-05-12 | 2005-03-15 | Fuji Photo Film Co., Ltd. | Image transformation method and apparatus, and storage medium |
| US7289678B2 (en) | 2000-05-12 | 2007-10-30 | Fujifilm Corporation | Image transformation method and apparatus, and storage medium |
| US7248748B2 (en) | 2000-11-30 | 2007-07-24 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method, storage medium, and program |
| US7076111B2 (en) | 2000-11-30 | 2006-07-11 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method, storage medium, and program |
| US7079700B2 (en) | 2000-11-30 | 2006-07-18 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method, storage medium, and program |
| WO2002045020A1 (fr) * | 2000-11-30 | 2002-06-06 | Canon Kabushiki Kaisha | Dispositif et procede de traitement d'images, support de stockage et programme associe |
| WO2002045019A1 (fr) * | 2000-11-30 | 2002-06-06 | Canon Kabushiki Kaisha | Dispositif de traitement d'image, procede de traitement d'image, support d'enregistrement et programme correspondant |
| US7447376B2 (en) | 2000-11-30 | 2008-11-04 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method, storage medium, and program |
| US7558434B2 (en) | 2000-11-30 | 2009-07-07 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method, storage medium, and program |
| US7561750B2 (en) | 2000-11-30 | 2009-07-14 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method, storage medium, and program |
| US7423781B2 (en) | 2002-03-20 | 2008-09-09 | Ricoh Company, Ltd. | Image processor and image processing method for image enhancement using edge detection |
| JP2005005845A (ja) * | 2003-06-10 | 2005-01-06 | Canon Inc | 画像処理装置 |
| JP2005109867A (ja) * | 2003-09-30 | 2005-04-21 | Konica Minolta Medical & Graphic Inc | 画像処理方法および画像処理装置ならびに画像処理プログラム |
| WO2009110260A1 (ja) * | 2008-03-06 | 2009-09-11 | コニカミノルタエムジー株式会社 | 画像処理装置 |
| JP5353876B2 (ja) * | 2008-03-06 | 2013-11-27 | コニカミノルタ株式会社 | 画像処理装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2663189B2 (ja) | 画像のダイナミックレンジ圧縮処理方法 | |
| US4747052A (en) | Radiation image processing | |
| JP3196864B2 (ja) | 画像のダイナミックレンジ圧縮処理方法 | |
| JP2849964B2 (ja) | 画像処理方法および装置 | |
| EP0442468A1 (en) | Method for forming energy substraction images | |
| JP3205417B2 (ja) | 画像処理方法 | |
| JPH0951427A (ja) | 画像データ補間演算方法および装置 | |
| JPH06274614A (ja) | 画像処理方法 | |
| JP2640582B2 (ja) | エネルギーサブトラクション画像生成方法 | |
| JPH0312785A (ja) | 画像処理方法 | |
| JP2754068B2 (ja) | 放射線画像信号生成方法および放射線画像読取装置 | |
| JPH08294006A (ja) | 画像のダイナミックレンジ圧縮処理方法および装置 | |
| JP2952418B2 (ja) | 放射線画像信号生成方法 | |
| US6614044B2 (en) | Image signal generating method, apparatus and program | |
| JPH0944659A (ja) | エネルギーサブトラクション処理方法および装置 | |
| JP2981694B2 (ja) | エネルギーサブトラクション画像生成方法 | |
| JPH11205682A (ja) | エネルギーサブトラクション画像生成方法 | |
| JP3399684B2 (ja) | 画像処理方法および装置 | |
| JPH01237887A (ja) | 放射線画像の粒状雑音低減方法 | |
| US7366358B2 (en) | Method of and system for generating image signal | |
| JPH08214164A (ja) | 階調補正方法および装置 | |
| JPH09160149A (ja) | 画像処理方法および装置 | |
| JP2739385B2 (ja) | 放射線画像読取条件及び／又は画像処理条件決定方法および装置 | |
| JP3241488B2 (ja) | 画像処理方法 | |
| JP3165461B2 (ja) | 画像処理方法、画像読取方法および画像再生方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20010814 |