JPH09212623A - 画像処理方法 - Google Patents
画像処理方法Info
- Publication number
- JPH09212623A JPH09212623A JP8014512A JP1451296A JPH09212623A JP H09212623 A JPH09212623 A JP H09212623A JP 8014512 A JP8014512 A JP 8014512A JP 1451296 A JP1451296 A JP 1451296A JP H09212623 A JPH09212623 A JP H09212623A
- Authority
- JP
- Japan
- Prior art keywords
- image
- wavelet transform
- signal
- frequency band
- transform coefficient
- 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.)
- Withdrawn
Links
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
(57)【要約】
【目的】 原画像データを多重解像変換してノイズを低
減する画像処理を施す画像処理方法において、画像中必
要な領域を低減することなく、ノイズ成分のみを低減さ
せる。 【構成】 原画像データ１に対してウェーブレット変換
３を施し、複数の周波数帯域毎の画像データ３を得る。
次いで、各画像データ３に対して所定の閾値以下の信号
値を０とする画像処理４を行う。その後画像処理４が施
された画像データ３に対して逆ウェーブレット変換５を
施し、処理済画像データ６を得る。
減する画像処理を施す画像処理方法において、画像中必
要な領域を低減することなく、ノイズ成分のみを低減さ
せる。 【構成】 原画像データ１に対してウェーブレット変換
３を施し、複数の周波数帯域毎の画像データ３を得る。
次いで、各画像データ３に対して所定の閾値以下の信号
値を０とする画像処理４を行う。その後画像処理４が施
された画像データ３に対して逆ウェーブレット変換５を
施し、処理済画像データ６を得る。
Description
【０００１】
【発明の属する技術分野】本発明は放射線画像における
所定の周波数帯域に画像処理を施す画像処理方法に関す
るものである。
所定の周波数帯域に画像処理を施す画像処理方法に関す
るものである。
【０００２】
【従来の技術】画像を表す画像信号を得、この画像信号
に適切な画像処理を施した後、画像を再生表示すること
が種々の分野で行われている。例えば放射線画像の診断
性能を向上させるために、画像信号に対してボケマスク
処理等の周波数強調処理を施す方法が本出願人により提
案されている（特開昭55-163772 等）。この周波数処理
は、原画像を表す画像信号からボケマスク信号を減算し
たものに強調度を乗じたものを加える処理を施すもの
で、これにより画像において所定の空間周波数成分を強
調するようにしたものである。
に適切な画像処理を施した後、画像を再生表示すること
が種々の分野で行われている。例えば放射線画像の診断
性能を向上させるために、画像信号に対してボケマスク
処理等の周波数強調処理を施す方法が本出願人により提
案されている（特開昭55-163772 等）。この周波数処理
は、原画像を表す画像信号からボケマスク信号を減算し
たものに強調度を乗じたものを加える処理を施すもの
で、これにより画像において所定の空間周波数成分を強
調するようにしたものである。
【０００３】一方、上述した画像信号を処理するための
方法として、画像を複数の周波数帯域毎の多重解像度画
像に変換し、各周波数帯域の画像に対して所定の処理を
行って、再度これを逆多重解像度変換することにより、
最終的な処理済画像を得るための多重解像度変換なる方
法が提案されている。この多重解像度変換の方法として
はウェーブレット変換、ラプラシアンピラミッドあるい
はフーリエ変換等の方法が知られている。
方法として、画像を複数の周波数帯域毎の多重解像度画
像に変換し、各周波数帯域の画像に対して所定の処理を
行って、再度これを逆多重解像度変換することにより、
最終的な処理済画像を得るための多重解像度変換なる方
法が提案されている。この多重解像度変換の方法として
はウェーブレット変換、ラプラシアンピラミッドあるい
はフーリエ変換等の方法が知られている。
【０００４】ここで、ウェーブレット変換について説明
する。
する。
【０００５】ウェーブレット変換は、周波数解析の方法
として近年開発されたものであり、ステレオのパターン
マッチング、データ圧縮等に応用がなされているもので
ある（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 ）。
【０００６】このウェーブレット変換は、図９に示すよ
うな関数ｈを基底関数として、
うな関数ｈを基底関数として、
【０００７】
【数１】
【０００８】なる式においてデータを複数の周波数帯域
毎の周波数信号に変換するため、フーリエ変換のような
偽振動の問題が発生しない。すなわち、関数ｈの周期お
よび縮率を変化させ、原信号を移動させることによりフ
ィルタリング処理を行えば、細かな周波数から粗い周波
数までの所望とする周波数に適合した周波数信号を作成
することができる。例えば、図１０に示すように、信号
Ｓorg をウェーブレット変換し、各周波数帯域毎に逆ウ
ェーブレット変換した信号と、図１１に示すように信号
Ｓorg をフーリエ変換し、各周波数帯域毎に逆フーリエ
変換した信号で見てみると、ウェーブレット変換はフー
リエ変換と比べて原信号Ｓorg の振動と対応した周波数
帯域の周波数信号を得ることができる。すなわち、フー
リエ変換において原信号Ｓorg の部分Ｂと対応する周波
数帯域７の部分Ｂ′には振動が発生しているのに対し、
ウェーブレット変換では原信号Ｓorg の部分Ａと対応す
る周波数帯域Ｗ７の部分Ａ′には原信号と同様に振動は
発生していないものとなる。
毎の周波数信号に変換するため、フーリエ変換のような
偽振動の問題が発生しない。すなわち、関数ｈの周期お
よび縮率を変化させ、原信号を移動させることによりフ
ィルタリング処理を行えば、細かな周波数から粗い周波
数までの所望とする周波数に適合した周波数信号を作成
することができる。例えば、図１０に示すように、信号
Ｓorg をウェーブレット変換し、各周波数帯域毎に逆ウ
ェーブレット変換した信号と、図１１に示すように信号
Ｓorg をフーリエ変換し、各周波数帯域毎に逆フーリエ
変換した信号で見てみると、ウェーブレット変換はフー
リエ変換と比べて原信号Ｓorg の振動と対応した周波数
帯域の周波数信号を得ることができる。すなわち、フー
リエ変換において原信号Ｓorg の部分Ｂと対応する周波
数帯域７の部分Ｂ′には振動が発生しているのに対し、
ウェーブレット変換では原信号Ｓorg の部分Ａと対応す
る周波数帯域Ｗ７の部分Ａ′には原信号と同様に振動は
発生していないものとなる。
【０００９】一方、ラプラシアンピラミッドなる方法は
例えば特開平5-244508号、特開平6-96200 、特開平6-30
1766号に記載されており、このラプラシアンピラミッド
は、原画像に対してガウス関数で近似されたようなマス
クによりマスク処理を施した後、画像をサブサンプリン
グして画素数を間引いて半分にすることにより、原画像
の１／４のサイズのボケ画像を得、このボケ画像のサン
プリングされた画素に値が０の画素を補間して元の大き
さの画像に戻し、この画像に対してさらに上述したマス
クによりマスク処理を施してボケ画像を得、このボケ画
像を原画像から減算して原画像の所定の周波数帯域を表
す細部画像を得るものである。この処理を得られたボケ
画像に対して繰り返すことにより原画像の１／２2Nの大
きさのボケ画像をＮ個作成するものである。ここで、ガ
ウス関数で近似されたようなマスクによりマスク処理を
施した画像に対してサンプリングを行っているため、実
際にはガウシアンフィルタを用いているが、ラプラシア
ンフィルタをかけた場合と同様の処理済画像が得られ
る。そしてこのように原画像サイズの画像から順に１／
２2Nの大きさの低周波数帯域の画像が得られるため、こ
の処理の結果得られた画像はラプラシアンピラミッドと
呼ばれる。
例えば特開平5-244508号、特開平6-96200 、特開平6-30
1766号に記載されており、このラプラシアンピラミッド
は、原画像に対してガウス関数で近似されたようなマス
クによりマスク処理を施した後、画像をサブサンプリン
グして画素数を間引いて半分にすることにより、原画像
の１／４のサイズのボケ画像を得、このボケ画像のサン
プリングされた画素に値が０の画素を補間して元の大き
さの画像に戻し、この画像に対してさらに上述したマス
クによりマスク処理を施してボケ画像を得、このボケ画
像を原画像から減算して原画像の所定の周波数帯域を表
す細部画像を得るものである。この処理を得られたボケ
画像に対して繰り返すことにより原画像の１／２2Nの大
きさのボケ画像をＮ個作成するものである。ここで、ガ
ウス関数で近似されたようなマスクによりマスク処理を
施した画像に対してサンプリングを行っているため、実
際にはガウシアンフィルタを用いているが、ラプラシア
ンフィルタをかけた場合と同様の処理済画像が得られ
る。そしてこのように原画像サイズの画像から順に１／
２2Nの大きさの低周波数帯域の画像が得られるため、こ
の処理の結果得られた画像はラプラシアンピラミッドと
呼ばれる。
【００１０】なお、このラプラシアンピラミッドについ
ては、Burt P.J.,“Fast Filter Transforms for Image
Processing”，Computer Graphics and Image Process
ing16巻、20〜51頁、1981年；Crowley J.L.,Stern R.
M.,“Fast Computation of the Difference of Low・Pa
ss Transform”IEEETrans.on Pattern Analysis andMac
hine Intelligence、６巻、２号、1984年３月、Mallat
S.G.,“A Theory forMultiresolution Signal Decompos
ition ；The Wavelet Representation”IEEETrans.on P
attern Analysis and Machine Intelligence 、11巻、
７号、1989年７月；Ebrahimi T.,Kunt M.,“Image comp
ression by Gabor Expansion”，Optical Engineering,
30巻、７号、873 〜880 頁、1991年７月、およびPieter
Vuylsteke,Emile Schoeters，“Multiscale Image Con
trast Amplification ”SPIEVol.2167 Image Processin
g(1994),pp551 〜560 に詳細が記載されている。
ては、Burt P.J.,“Fast Filter Transforms for Image
Processing”，Computer Graphics and Image Process
ing16巻、20〜51頁、1981年；Crowley J.L.,Stern R.
M.,“Fast Computation of the Difference of Low・Pa
ss Transform”IEEETrans.on Pattern Analysis andMac
hine Intelligence、６巻、２号、1984年３月、Mallat
S.G.,“A Theory forMultiresolution Signal Decompos
ition ；The Wavelet Representation”IEEETrans.on P
attern Analysis and Machine Intelligence 、11巻、
７号、1989年７月；Ebrahimi T.,Kunt M.,“Image comp
ression by Gabor Expansion”，Optical Engineering,
30巻、７号、873 〜880 頁、1991年７月、およびPieter
Vuylsteke,Emile Schoeters，“Multiscale Image Con
trast Amplification ”SPIEVol.2167 Image Processin
g(1994),pp551 〜560 に詳細が記載されている。
【００１１】そしてこのようにして得られたラプラシア
ンピラミッドの全ての周波数帯域の画像に対して、画像
の値を強調するような処理を施し、この強調処理が施さ
れた各周波数帯域の画像を逆変換して処理済画像を得る
方法が上記特開平6-301766号に記載されている。この方
法は、各周波数帯域の画像信号に対して、下記の式 ｙ＝−ｍ×（−ｘ／ｍ）p （ｘ＜０） ｙ＝ｍ×（−ｘ／ｍ）p （ｘ≧０） 但し、ｘ：画像の各画素における画素値 ｙ：強調処理が施された画像の各画素における画素値 ｍ：画素のとりうる値の範囲（例えば、画素のとりうる
値の範囲が１０ビットである場合ｍ＝１０２３とな
る。） により画像の強調を行うものである。すなわち、ｐの値
が小さいほど強調度が大きく、ｐの値が大きいほど強調
度が小さくされて画像の強調が行われる。そしてこのよ
うに処理が施された画像は、各周波数帯域において画像
が強調されているため、実質的に上述したボケマスク処
理において複数のサイズのマスクによりボケマスク処理
を施したような画像となっている。
ンピラミッドの全ての周波数帯域の画像に対して、画像
の値を強調するような処理を施し、この強調処理が施さ
れた各周波数帯域の画像を逆変換して処理済画像を得る
方法が上記特開平6-301766号に記載されている。この方
法は、各周波数帯域の画像信号に対して、下記の式 ｙ＝−ｍ×（−ｘ／ｍ）p （ｘ＜０） ｙ＝ｍ×（−ｘ／ｍ）p （ｘ≧０） 但し、ｘ：画像の各画素における画素値 ｙ：強調処理が施された画像の各画素における画素値 ｍ：画素のとりうる値の範囲（例えば、画素のとりうる
値の範囲が１０ビットである場合ｍ＝１０２３とな
る。） により画像の強調を行うものである。すなわち、ｐの値
が小さいほど強調度が大きく、ｐの値が大きいほど強調
度が小さくされて画像の強調が行われる。そしてこのよ
うに処理が施された画像は、各周波数帯域において画像
が強調されているため、実質的に上述したボケマスク処
理において複数のサイズのマスクによりボケマスク処理
を施したような画像となっている。
【００１２】一方、放射線画像においては、放射線量が
少なく濃度が低い部分において、放射線の量子ノイズが
目立ってしまう。このため、上記特開平6-96200 号にお
いては、放射線画像をラプラシアンピラミッドにより複
数の周波数帯域の画像に分解し、各周波数帯域の画像に
対して、各周波数帯域の画像信号の局所分散値を算出
し、この局所分散値の大きさに応じて画像信号に対して
非線形処理を施し、画像中のノイズ成分を低減するよう
にしている。
少なく濃度が低い部分において、放射線の量子ノイズが
目立ってしまう。このため、上記特開平6-96200 号にお
いては、放射線画像をラプラシアンピラミッドにより複
数の周波数帯域の画像に分解し、各周波数帯域の画像に
対して、各周波数帯域の画像信号の局所分散値を算出
し、この局所分散値の大きさに応じて画像信号に対して
非線形処理を施し、画像中のノイズ成分を低減するよう
にしている。
【００１３】
【発明が解決しようとする課題】上記特開平6-96200 号
に記載された方法においては、各周波数帯域の画像信号
の局所分散値に基づいて処理を行うようにしているた
め、放射線画像中にノイズ以外の必要な成分が含まれて
いる場合には、分散値が大きくなり、よりノイズを低減
することができることとなるが、同時に必要な成分をも
除去することとなってしまう。このため、処理を施すこ
とにより得られた画像において、ノイズだけではなく画
像中必要な成分までもが除去されてしまい、処理がなさ
れた放射線画像を観察する上で好ましくない。
に記載された方法においては、各周波数帯域の画像信号
の局所分散値に基づいて処理を行うようにしているた
め、放射線画像中にノイズ以外の必要な成分が含まれて
いる場合には、分散値が大きくなり、よりノイズを低減
することができることとなるが、同時に必要な成分をも
除去することとなってしまう。このため、処理を施すこ
とにより得られた画像において、ノイズだけではなく画
像中必要な成分までもが除去されてしまい、処理がなさ
れた放射線画像を観察する上で好ましくない。
【００１４】本発明は上記事情に鑑み、放射線画像中の
必要な成分を除去することなく、ノイズ成分のみを除去
することができる画像処理方法を提供することを目的と
するものである。
必要な成分を除去することなく、ノイズ成分のみを除去
することができる画像処理方法を提供することを目的と
するものである。
【００１５】
【課題を解決するための手段】本発明による画像処理方
法は、放射線画像を表す画像信号を多重解像度空間に変
換することにより、該画像信号を複数の周波数帯域毎の
画像を表す周波数帯域画像信号に分解し、該複数の周波
数帯域画像信号のうち所定の周波数帯域画像信号に対し
て、該各信号値の絶対値が所定の閾値以下の信号値を０
とする処理を施し、該処理が施された周波数帯域画像信
号および他の周波数帯域画像信号を逆変換することによ
り処理済画像信号を得ることを特徴とするものである。
法は、放射線画像を表す画像信号を多重解像度空間に変
換することにより、該画像信号を複数の周波数帯域毎の
画像を表す周波数帯域画像信号に分解し、該複数の周波
数帯域画像信号のうち所定の周波数帯域画像信号に対し
て、該各信号値の絶対値が所定の閾値以下の信号値を０
とする処理を施し、該処理が施された周波数帯域画像信
号および他の周波数帯域画像信号を逆変換することによ
り処理済画像信号を得ることを特徴とするものである。
【００１６】また、上記画像処理方法においては、前記
周波数帯域画像信号のうち比較的中間周波数帯域の周波
数帯域画像信号について、前記所定の閾値を他の周波数
帯域の閾値よりも大きくすることが好ましい。
周波数帯域画像信号のうち比較的中間周波数帯域の周波
数帯域画像信号について、前記所定の閾値を他の周波数
帯域の閾値よりも大きくすることが好ましい。
【００１７】さらに、前記周波数帯域画像信号の比較的
低濃度部分の前記閾値を前記周波数帯域画像信号の比較
的高濃度部の前記閾値よりも大きくすることが好まし
い。
低濃度部分の前記閾値を前記周波数帯域画像信号の比較
的高濃度部の前記閾値よりも大きくすることが好まし
い。
【００１８】
【発明の効果】放射線画像において、ノイズが目立つの
は比較的信号値が小さい低濃度部分である。本発明によ
る画像処理方法は、多重解像度空間に変換された複数の
周波数帯域の画像信号に対して、この画像信号値の絶対
値が所定の閾値以下の信号を０とする処理を施すように
したため、ノイズが目立つ比較的信号値の低濃度部分は
信号値が０となる。したがって、画像中のノイズと見な
せる低濃度部分を０とすることができ、これにより画像
中のノイズ成分を除去することができる。
は比較的信号値が小さい低濃度部分である。本発明によ
る画像処理方法は、多重解像度空間に変換された複数の
周波数帯域の画像信号に対して、この画像信号値の絶対
値が所定の閾値以下の信号を０とする処理を施すように
したため、ノイズが目立つ比較的信号値の低濃度部分は
信号値が０となる。したがって、画像中のノイズと見な
せる低濃度部分を０とすることができ、これにより画像
中のノイズ成分を除去することができる。
【００１９】また、比較的中間周波数帯域の画像信号に
対しては、他の周波数帯域の画像信号よりも閾値を高く
することにより、ノイズが目立つ中間周波数帯域につい
ては、より濃度の高い部分から信号値が０とされるた
め、ノイズをより一層低減することができる。
対しては、他の周波数帯域の画像信号よりも閾値を高く
することにより、ノイズが目立つ中間周波数帯域につい
ては、より濃度の高い部分から信号値が０とされるた
め、ノイズをより一層低減することができる。
【００２０】さらに、ノイズを低減させる際に、同時に
画像中の必要な信号成分をも低減させてしまうことがあ
るが、低濃度部の閾値を高濃度部の閾値よりも大きくす
ることにより、ノイズを低減する必要のない高濃度域に
おいて、ノイズはほとんど低減されず、ノイズが目立つ
低濃度域においてのみ低減されることとなり、画像中の
ノイズをより効果的に低減することができる。
画像中の必要な信号成分をも低減させてしまうことがあ
るが、低濃度部の閾値を高濃度部の閾値よりも大きくす
ることにより、ノイズを低減する必要のない高濃度域に
おいて、ノイズはほとんど低減されず、ノイズが目立つ
低濃度域においてのみ低減されることとなり、画像中の
ノイズをより効果的に低減することができる。
【００２１】
【発明の実施の形態】以下図面を参照して本発明の実施
の形態について説明する。
の形態について説明する。
【００２２】図１は本発明による画像処理方法の基本的
概念を表す図である。図１に示すように、本発明による
画像処理方法は、放射線画像を表す原画像データ１に対
してウェーブレット変換２を施して複数の周波数帯域毎
の画像データ３を得る。次いで、各画像データ３に対し
て後述するような画像処理４を施す。そしてこの画像処
理が施された画像データ３に対して逆ウェーブレット変
換５を施して最終的な処理済画像データ６を得るもので
ある。
概念を表す図である。図１に示すように、本発明による
画像処理方法は、放射線画像を表す原画像データ１に対
してウェーブレット変換２を施して複数の周波数帯域毎
の画像データ３を得る。次いで、各画像データ３に対し
て後述するような画像処理４を施す。そしてこの画像処
理が施された画像データ３に対して逆ウェーブレット変
換５を施して最終的な処理済画像データ６を得るもので
ある。
【００２３】以下本発明による実施の形態の詳細につい
て説明する。
て説明する。
【００２４】本実施の形態は、例えば特開昭55-12492号
公報や特開昭56-11395号等に記録されている蓄積性蛍光
体シートを利用した放射線画像情報記録再生システムに
おいて、蓄積性蛍光体シートに記録された人体の放射線
画像をレーザビーム走査によりデジタル画像データとし
て読み取ったものを対象としている。なお、放射線画像
の読み取りは、図２に示すように、蓄積性蛍光体シート
10に対して主走査方向（横方向）にレーザビームを走査
させながらシート10を副走査方向（縦方向）に移動させ
てシート10を２次元走査することにより行われたもので
ある。
公報や特開昭56-11395号等に記録されている蓄積性蛍光
体シートを利用した放射線画像情報記録再生システムに
おいて、蓄積性蛍光体シートに記録された人体の放射線
画像をレーザビーム走査によりデジタル画像データとし
て読み取ったものを対象としている。なお、放射線画像
の読み取りは、図２に示すように、蓄積性蛍光体シート
10に対して主走査方向（横方向）にレーザビームを走査
させながらシート10を副走査方向（縦方向）に移動させ
てシート10を２次元走査することにより行われたもので
ある。
【００２５】次いで、原画像データに対してウェーブレ
ット変換がなされる。
ット変換がなされる。
【００２６】図３は、原画像データＳorg に対するウェ
ーブレット変換の詳細を表す図である。
ーブレット変換の詳細を表す図である。
【００２７】なお、本実施の形態においては、ウェーブ
レット変換の各係数が直交する直交ウェーブレット変換
を行うものであり、前述したMarc Antonini らの文献に
記載されているものである。
レット変換の各係数が直交する直交ウェーブレット変換
を行うものであり、前述したMarc Antonini らの文献に
記載されているものである。
【００２８】図３に示すように、原画像データＳorg の
主走査方向に基本ウェーブレット関数より求められる関
数ｇと関数ｈとによりフィルタリング処理を行う。すな
わち、このような関数ｇ，ｈによる主走査方向に並ぶ画
素の一列毎のフィルタリング処理を副走査方向に一画素
ずつズラしながら行い、原画像データＳorg の主走査方
向のウェーブレット変換係数信号Ｗg0，Ｗh0を求めるも
のである。
主走査方向に基本ウェーブレット関数より求められる関
数ｇと関数ｈとによりフィルタリング処理を行う。すな
わち、このような関数ｇ，ｈによる主走査方向に並ぶ画
素の一列毎のフィルタリング処理を副走査方向に一画素
ずつズラしながら行い、原画像データＳorg の主走査方
向のウェーブレット変換係数信号Ｗg0，Ｗh0を求めるも
のである。
【００２９】ここで、関数ｇ，ｈは基本ウェーブレット
関数より一意に求められるものであり、例えば、関数ｈ
は、以下の表１に示すものとなる。なお、表１において
関数ｈ′は、ウェーブレット変換がなされた画像データ
に逆ウェーブレット変換を行う際に用いる関数を表すも
のである。また以下の式(2) に示すように関数ｇは関数
ｈ′から求められ、逆ウェーブレット変換を行うための
関数ｇ′は関数ｈから求められる。
関数より一意に求められるものであり、例えば、関数ｈ
は、以下の表１に示すものとなる。なお、表１において
関数ｈ′は、ウェーブレット変換がなされた画像データ
に逆ウェーブレット変換を行う際に用いる関数を表すも
のである。また以下の式(2) に示すように関数ｇは関数
ｈ′から求められ、逆ウェーブレット変換を行うための
関数ｇ′は関数ｈから求められる。
【００３０】
【表１】
【００３１】 ｇ′＝（−１）n ｈ ｇ ＝（−１）n ｈ′ …(2) このようにして、ウェーブレット変換係数信号Ｗg0、Ｗ
h0が求められると、ウェーブレット変換係数信号Ｗg0、
Ｗh0について、主走査方向の画素を１画素おきに間引
き、主走査方向の画素数を1/2 にする。ついで、この画
素が間引かれたウェーブレット変換係数信号Ｗg0、Ｗh0
それぞれの副走査方向に関数ｇ，ｈによりフィルタリン
グ処理を行い、ウェーブレット変換係数信号ＷＷ0 ，Ｗ
Ｖ0 ，ＶＷ0 およびＶＶ0 を得る。
h0が求められると、ウェーブレット変換係数信号Ｗg0、
Ｗh0について、主走査方向の画素を１画素おきに間引
き、主走査方向の画素数を1/2 にする。ついで、この画
素が間引かれたウェーブレット変換係数信号Ｗg0、Ｗh0
それぞれの副走査方向に関数ｇ，ｈによりフィルタリン
グ処理を行い、ウェーブレット変換係数信号ＷＷ0 ，Ｗ
Ｖ0 ，ＶＷ0 およびＶＶ0 を得る。
【００３２】次いでウェーブレット変換係数信号Ｗ
Ｗ0 ，ＷＶ0 ，ＶＷ0 およびＶＶ0 について、副走査方
向の画素を１画素おきに間引くことを行い、副走査方向
の画素数を1/2 とする処理を行う。これにより、各ウェ
ーブレット変換係数信号ＶＶ0 ，ＷＶ0 ，ＶＷ0 ，ＷＷ
0 の画素数は原画像データＳorg の画素数の1/4 とな
る。次いで、ウェーブレット変換係数信号ＶＶ0 の主走
査方向に関数ｇ，ｈによりフィルタリング処理を行う。
Ｗ0 ，ＷＶ0 ，ＶＷ0 およびＶＶ0 について、副走査方
向の画素を１画素おきに間引くことを行い、副走査方向
の画素数を1/2 とする処理を行う。これにより、各ウェ
ーブレット変換係数信号ＶＶ0 ，ＷＶ0 ，ＶＷ0 ，ＷＷ
0 の画素数は原画像データＳorg の画素数の1/4 とな
る。次いで、ウェーブレット変換係数信号ＶＶ0 の主走
査方向に関数ｇ，ｈによりフィルタリング処理を行う。
【００３３】すなわち、関数ｇ，ｈにより主走査方向に
並ぶ画素の一列毎のフィルタリング処理を副走査方向に
一画素ずつズラしながら行い、ウェーブレット変換係数
信号ＶＶ0 の主走査方向のウェーブレット変換係数信号
Ｗg1およびＷh1を求めるものである。
並ぶ画素の一列毎のフィルタリング処理を副走査方向に
一画素ずつズラしながら行い、ウェーブレット変換係数
信号ＶＶ0 の主走査方向のウェーブレット変換係数信号
Ｗg1およびＷh1を求めるものである。
【００３４】ここでウェーブレット変換係数信号ＶＶ0
は主副両方向について画素数が原画像データの1/2 とな
っているため、画像の解像度は原画像データと比較して
半分となっている。したがって、ウェーブレット変換係
数信号ＶＶ0 を関数ｇ，ｈでフィルタリング処理を施す
ことにより、原画像データの周波数成分のうちウェーブ
レット変換係数信号ＶＶ0 が表す周波数成分よりも低周
波数成分を表すウェーブレット変換係数信号Ｗg1，Ｗh1
が求められる。
は主副両方向について画素数が原画像データの1/2 とな
っているため、画像の解像度は原画像データと比較して
半分となっている。したがって、ウェーブレット変換係
数信号ＶＶ0 を関数ｇ，ｈでフィルタリング処理を施す
ことにより、原画像データの周波数成分のうちウェーブ
レット変換係数信号ＶＶ0 が表す周波数成分よりも低周
波数成分を表すウェーブレット変換係数信号Ｗg1，Ｗh1
が求められる。
【００３５】このようにして、ウェーブレット変換係数
信号Ｗg1，Ｗh1が求められると、ウェーブレット変換係
数信号Ｗg1，Ｗh1について、主走査方向の画素を１画素
おきに間引き、主走査方向の画素数をさらに1/2 とす
る。次いでウェーブレット変換係数信号Ｗg1、Ｗh1それ
ぞれの副走査方向に関数ｇ，ｈによりフィルタリング処
理を行い、ウェーブレット変換係数信号ＷＷ1 ，Ｗ
Ｖ1 ，ＶＷ1 およびＶＶ1 を得る。
信号Ｗg1，Ｗh1が求められると、ウェーブレット変換係
数信号Ｗg1，Ｗh1について、主走査方向の画素を１画素
おきに間引き、主走査方向の画素数をさらに1/2 とす
る。次いでウェーブレット変換係数信号Ｗg1、Ｗh1それ
ぞれの副走査方向に関数ｇ，ｈによりフィルタリング処
理を行い、ウェーブレット変換係数信号ＷＷ1 ，Ｗ
Ｖ1 ，ＶＷ1 およびＶＶ1 を得る。
【００３６】次いでウェーブレット変換係数信号Ｗ
Ｗ1 ，ＷＶ1 ，ＶＷ1 ，ＶＶ1 について、副走査方向の
画素を１画素おきに間引き、副走査方向の画素数を1/2
とする処理を行う。これにより、各ウェーブレット変換
係数信号ＶＶ1 ，ＷＶ1 ，ＶＷ1，ＷＷ1 の画素数は原
画像データＳorg の画素数の1/16となる。
Ｗ1 ，ＷＶ1 ，ＶＷ1 ，ＶＶ1 について、副走査方向の
画素を１画素おきに間引き、副走査方向の画素数を1/2
とする処理を行う。これにより、各ウェーブレット変換
係数信号ＶＶ1 ，ＷＶ1 ，ＶＷ1，ＷＷ1 の画素数は原
画像データＳorg の画素数の1/16となる。
【００３７】以下、上述したのと同様にして、画素が間
引かれたウェーブレット変換係数信号ＶＶ1 の主走査方
向に関数ｇ，ｈによりフィルタリング処理を行い、さら
に得られたウェーブレット変換係数信号の主走査方向の
画素を間引き、この画素を間引いたウェーブレット変換
係数信号について、副走査方向に関数ｇ，ｈによりフィ
ルタリング処理を行い、ウェーブレット変換係数信号Ｗ
Ｗ2 ，ＷＶ2 ，ＶＷ2，ＶＶ2 を得る。
引かれたウェーブレット変換係数信号ＶＶ1 の主走査方
向に関数ｇ，ｈによりフィルタリング処理を行い、さら
に得られたウェーブレット変換係数信号の主走査方向の
画素を間引き、この画素を間引いたウェーブレット変換
係数信号について、副走査方向に関数ｇ，ｈによりフィ
ルタリング処理を行い、ウェーブレット変換係数信号Ｗ
Ｗ2 ，ＷＶ2 ，ＶＷ2，ＶＶ2 を得る。
【００３８】このようなウェーブレット変換をＮ回繰り
返すことによりウェーブレット変換係数信号ＷＷ0 〜Ｗ
ＷN ，ＷＶ0 〜ＷＶN ，ＶＷ0 〜ＶＷN ，およびＶＶN
を得る。ここで、Ｎ回目のウェーブレット変換により得
られるウェーブレット変換係数信号ＷＷN ，ＷＶN ，Ｖ
ＷN ，ＶＶN は、原画像データと比較して主副各方向の
画素数が(1/2）N となっているため、各ウェーブレット
変換係数信号はＮが大きいほど周波数帯域が低く、原画
像データの周波数成分のうち低周波数成分を表すデータ
となる。
返すことによりウェーブレット変換係数信号ＷＷ0 〜Ｗ
ＷN ，ＷＶ0 〜ＷＶN ，ＶＷ0 〜ＶＷN ，およびＶＶN
を得る。ここで、Ｎ回目のウェーブレット変換により得
られるウェーブレット変換係数信号ＷＷN ，ＷＶN ，Ｖ
ＷN ，ＶＶN は、原画像データと比較して主副各方向の
画素数が(1/2）N となっているため、各ウェーブレット
変換係数信号はＮが大きいほど周波数帯域が低く、原画
像データの周波数成分のうち低周波数成分を表すデータ
となる。
【００３９】したがって、ウェーブレット変換係数信号
ＷＷi （ｉ＝０〜Ｎ、以下同様）は、原画像データＳor
g の主副両方向の周波数の変化を表すものであり、ｉが
大きいほど低周波信号となる。またウェーブレット変換
係数信号ＷＶi は画像信号Ｓorg の主走査方向の周波数
の変化を表すものであり、ｉが大きいほど低周波信号と
なる。さらにウェーブレット変換係数信号ＶＷi は画像
信号Ｓorg の副走査方向の周波数の変化を表すものであ
り、ｉが大きいほど低周波信号となる。
ＷＷi （ｉ＝０〜Ｎ、以下同様）は、原画像データＳor
g の主副両方向の周波数の変化を表すものであり、ｉが
大きいほど低周波信号となる。またウェーブレット変換
係数信号ＷＶi は画像信号Ｓorg の主走査方向の周波数
の変化を表すものであり、ｉが大きいほど低周波信号と
なる。さらにウェーブレット変換係数信号ＶＷi は画像
信号Ｓorg の副走査方向の周波数の変化を表すものであ
り、ｉが大きいほど低周波信号となる。
【００４０】ここで、図４にウェーブレット変換係数信
号を複数の周波数帯域毎に表す図を示す。なお、図４に
おいては便宜上３回目のウェーブレット変換を行った状
態までを表すものとする。なお、図４においてウェーブ
レット変換係数信号ＷＷ3 は原画像を主副各方向が(1/
2）3 に縮小したものとなっている。
号を複数の周波数帯域毎に表す図を示す。なお、図４に
おいては便宜上３回目のウェーブレット変換を行った状
態までを表すものとする。なお、図４においてウェーブ
レット変換係数信号ＷＷ3 は原画像を主副各方向が(1/
2）3 に縮小したものとなっている。
【００４１】次いで、ウェーブレット変換係数信号ＷＶ
i ，ＶＷi ，ＷＷi に対して画像処理がなされる。以下
この画像処理について説明する。
i ，ＶＷi ，ＷＷi に対して画像処理がなされる。以下
この画像処理について説明する。
【００４２】本実施の形態によりウェーブレット変換が
施された放射線画像において、ノイズが目立つのは比較
的信号値が小さい低濃度部分である。したがって、本発
明による画像処理方法においては、上述したウェーブレ
ット変換により得られた各ウェーブレット変換係数信号
ＷＶi ，ＶＷi ，ＷＷi に対して、信号値の絶対値が所
定の閾値以下の信号を０とする処理が施される。この処
理は、各ウェーブレット変換係数信号ＷＶi ，ＶＷi ，
ＷＷi に対して、図５に示すテーブルにより信号値を処
理することにより行われる。図５に示すテーブルは、各
ウェーブレット変換係数信号ＷＶi ，ＶＷi ，ＷＷi の
信号値の絶対値が閾値Ｔｈ１より小さい信号値を０とす
るように処理を行うものである。なお、放射線画像の中
間周波数帯域においては他の周波数帯域よりもノイズが
目立つため、この中間周波数帯域においては他の周波数
帯域よりも閾値を大きく設定することが好ましい。例え
ば図５の閾値Ｔｈ２のように閾値の絶対値がＴｈ１より
も大きい値を設定することが好ましい。
施された放射線画像において、ノイズが目立つのは比較
的信号値が小さい低濃度部分である。したがって、本発
明による画像処理方法においては、上述したウェーブレ
ット変換により得られた各ウェーブレット変換係数信号
ＷＶi ，ＶＷi ，ＷＷi に対して、信号値の絶対値が所
定の閾値以下の信号を０とする処理が施される。この処
理は、各ウェーブレット変換係数信号ＷＶi ，ＶＷi ，
ＷＷi に対して、図５に示すテーブルにより信号値を処
理することにより行われる。図５に示すテーブルは、各
ウェーブレット変換係数信号ＷＶi ，ＶＷ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′を得る。
【００４６】このようにウェーブレット変換を行う関数
と逆ウェーブレット変換を行う関数とを異なるものとし
ているのは、以下のような理由からである。ウェーブレ
ット変換と逆ウェーブレット変換で同一の関数となる、
すなわち、直交する関数を設計することは難しく、直交
性、連続性、関数の短さ、対称性のいずれかの条件を緩
める必要がある。そこで、直交性の条件を緩めることに
より他の条件を満たす関数を選択したものである。
と逆ウェーブレット変換を行う関数とを異なるものとし
ているのは、以下のような理由からである。ウェーブレ
ット変換と逆ウェーブレット変換で同一の関数となる、
すなわち、直交する関数を設計することは難しく、直交
性、連続性、関数の短さ、対称性のいずれかの条件を緩
める必要がある。そこで、直交性の条件を緩めることに
より他の条件を満たす関数を選択したものである。
【００４７】以上より、本実施の形態ではウェーブレッ
ト変換を行う関数ｈ，ｇと逆ウェーブレット変換を行う
関数ｈ′，ｇ′とを双直交の異なるものとしている。し
たがって、ウェーブレット変換係数信号ＶＶi ，Ｖ
Ｗi ，ＷＶi ，ＷＷi を関数ｈ′，ｇ′で逆ウェーブレ
ット変換することにより、原画像データを完全に復元で
きることとなる。
ト変換を行う関数ｈ，ｇと逆ウェーブレット変換を行う
関数ｈ′，ｇ′とを双直交の異なるものとしている。し
たがって、ウェーブレット変換係数信号ＶＶi ，Ｖ
Ｗi ，ＷＶi ，ＷＷi を関数ｈ′，ｇ′で逆ウェーブレ
ット変換することにより、原画像データを完全に復元で
きることとなる。
【００４８】一方、これと並列して、ウェーブレット変
換係数信号ＷＶN を副走査方向に関数ｈ′により、ウェ
ーブレット変換係数信号ＷＷN を副走査方向に関数ｇ′
によりフィルタリング処理を行い、ウェーブレット変換
係数信号ＷＶN ，ＷＷN の逆ウェーブレット変換係数信
号を得、これを２倍して加算することにより逆ウェーブ
レット変換係数信号ＷgN′を得る。
換係数信号ＷＶN を副走査方向に関数ｈ′により、ウェ
ーブレット変換係数信号ＷＷN を副走査方向に関数ｇ′
によりフィルタリング処理を行い、ウェーブレット変換
係数信号ＷＶN ，ＷＷN の逆ウェーブレット変換係数信
号を得、これを２倍して加算することにより逆ウェーブ
レット変換係数信号ＷgN′を得る。
【００４９】次いで、逆ウェーブレット変換係数信号Ｗ
hN′，ＷgN′について主走査方向に並ぶ画素間に１画素
分の間隔をあける処理を行う。その後逆ウェーブレット
変換係数信号ＷhN′を主走査方向に関数ｈ′により、逆
ウェーブレット変換係数信号ＷgN′を主走査方向に関数
ｇ′によりフィルタリング処理し、ウェーブレット変換
係数信号ＷhN′，ＷgN′の逆ウェーブレット変換係数信
号を得、これを２倍して加算することにより逆ウェーブ
レット変換係数信号ＶＶN-1 ′を得る。
hN′，ＷgN′について主走査方向に並ぶ画素間に１画素
分の間隔をあける処理を行う。その後逆ウェーブレット
変換係数信号ＷhN′を主走査方向に関数ｈ′により、逆
ウェーブレット変換係数信号ＷgN′を主走査方向に関数
ｇ′によりフィルタリング処理し、ウェーブレット変換
係数信号ＷhN′，Ｗ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 の逆ウェーブレット変
換係数信号を得、これを２倍して加算することにより逆
ウェーブレット変換係数信号ＷgN-1′を得る。
換係数信号ＷＶN-1 を副走査方向に関数ｈ′により、ウ
ェーブレット変換係数信号ＷＷN-1 を副走査方向に関数
ｇ′によりフィルタリング処理を行い、ウェーブレット
変換係数信号ＷＶN-1 ，ＷＷN-1 の逆ウェーブレット変
換係数信号を得、これを２倍して加算することにより逆
ウェーブレット変換係数信号ＷgN-1′を得る。
【００５２】次いで、逆ウェーブレット変換係数信号Ｗ
hN-1′，ＷgN-1′について主走査方向に並ぶ画素間に１
画素分の間隔をあける処理を行う。その後逆ウェーブレ
ット変換係数信号ＷhN-1′を主走査方向に関数ｈ′によ
り、逆ウェーブレット変換係数信号ＷgN-1′を主走査方
向に関数ｇ′によりフィルタリング処理し、ウェーブレ
ット変換係数信号ＷhN-1′，ＷgN-1′の逆ウェーブレッ
ト変換係数信号を得、これを２倍して加算することによ
り逆ウェーブレット変換係数信号ＶＶN-2 ′を得る。
hN-1′，ＷgN-1′について主走査方向に並ぶ画素間に１
画素分の間隔をあける処理を行う。その後逆ウェーブレ
ット変換係数信号ＷhN-1′を主走査方向に関数ｈ′によ
り、逆ウェーブレット変換係数信号ＷgN-1′を主走査方
向に関数ｇ′によりフィルタリング処理し、ウェーブレ
ット変換係数信号ＷhN-1′，ＷgN-1′の逆ウェーブレッ
ト変換係数信号を得、これを２倍して加算することによ
り逆ウェーブレット変換係数信号ＶＶN-2 ′を得る。
【００５３】以下、順次逆ウェーブレット変換係数信号
ＶＶi ′（ｉ＝−１〜Ｎ）を作成し、最終的に逆ウェー
ブレット変換係数信号ＶＶ-1′を得る。この最終的な逆
ウェーブレット変換係数信号ＶＶ-1′が処理済画像デー
タ６を表すものとなる。
ＶＶi ′（ｉ＝−１〜Ｎ）を作成し、最終的に逆ウェー
ブレット変換係数信号ＶＶ-1′を得る。この最終的な逆
ウェーブレット変換係数信号ＶＶ-1′が処理済画像デー
タ６を表すものとなる。
【００５４】このようにして得られたウェーブレット変
換係数信号ＶＶ-1′は図示しない画像再生装置に送られ
て、放射線画像の再生に供せられる。
換係数信号ＶＶ-1′は図示しない画像再生装置に送られ
て、放射線画像の再生に供せられる。
【００５５】この再生装置は、ＣＲＴ等のディスプレイ
手段でもよいし、感光フイルムに光走査記録を行う記録
装置であってもよい。
手段でもよいし、感光フイルムに光走査記録を行う記録
装置であってもよい。
【００５６】ここで、再生された画像における比較的低
濃度の部分は、上述したように信号値を０とする処理が
施されているため、再生された画像におけるノイズ成分
は低減されているものである。したがって、ノイズのな
い高画質の画像を再生されることとなる。
濃度の部分は、上述したように信号値を０とする処理が
施されているため、再生された画像におけるノイズ成分
は低減されているものである。したがって、ノイズのな
い高画質の画像を再生されることとなる。
【００５７】なお、上述した処理において原画像データ
に対してウェーブレット変換を施す前に、図７に示すよ
うなテーブルにより原画像データを変換することが考え
られる。このようなテーブルにより原画像データを変換
することにより、低濃度部の信号コントラストを高濃度
部と比較して相対的に低下させることができる。そして
変換された原画像データに対してウェーブレット変換を
施した後、上述した図５に示すテーブルにより各ウェー
ブレット変換係数信号ＷＶi ，ＶＷi ，ＷＷiに対して
閾値処理を行うと、低濃度部は高濃度部と比較して相対
的に高い閾値処理を施したこととなる。そして、閾値を
調整することにより高濃度部に対してはほとんど閾値処
理が施されないようにすることもできる。そして処理が
施されたウェーブレット変換係数信号ＷＶi ，ＶＷi ，
ＷＷi を逆ウェーブレット変換して逆ウェーブレット変
換信号を得、この逆ウェーブレット変換信号に対して図
７により低下されたコントラストを元のコントラストに
戻す処理を行う。この処理は図８に示すテーブルにより
逆ウェーブレット変換信号を変換することにより行う。
そしてこの図８に示すテーブルにより処理された信号が
最終的な処理済画像データ６となる。この処理済画像デ
ータ６は上述した実施の形態と同様に再生に供される。
に対してウェーブレット変換を施す前に、図７に示すよ
うなテーブルにより原画像データを変換することが考え
られる。このようなテーブルにより原画像データを変換
することにより、低濃度部の信号コントラストを高濃度
部と比較して相対的に低下させることができる。そして
変換された原画像データに対してウェーブレット変換を
施した後、上述した図５に示すテーブルにより各ウェー
ブレット変換係数信号ＷＶi ，ＶＷi ，ＷＷiに対して
閾値処理を行うと、低濃度部は高濃度部と比較して相対
的に高い閾値処理を施したこととなる。そして、閾値を
調整することにより高濃度部に対してはほとんど閾値処
理が施されないようにすることもできる。そして処理が
施されたウェーブレット変換係数信号ＷＶi ，ＶＷi ，
ＷＷi を逆ウェーブレット変換して逆ウェーブレット変
換信号を得、この逆ウェーブレット変換信号に対して図
７により低下されたコントラストを元のコントラストに
戻す処理を行う。この処理は図８に示すテーブルにより
逆ウェーブレット変換信号を変換することにより行う。
そしてこの図８に示すテーブルにより処理された信号が
最終的な処理済画像データ６となる。この処理済画像デ
ータ６は上述した実施の形態と同様に再生に供される。
【００５８】このような処理を行うことにより、低濃度
部のノイズのみがより一層低減され、高濃度部に対して
はほとんどノイズが低減されないようにすることがで
き、これによりノイズが目立つ低濃度領域のみを効果的
に低減することができる。
部のノイズのみがより一層低減され、高濃度部に対して
はほとんどノイズが低減されないようにすることがで
き、これによりノイズが目立つ低濃度領域のみを効果的
に低減することができる。
【００５９】また、上述した実施の形態においては、ウ
ェーブレット変換を行うための関数ｈ，ｈ′として表１
に示すものを用いたが、これに限定されるものではなく
以下に示す表２、表３に示すものを用いてもよい。
ェーブレット変換を行うための関数ｈ，ｈ′として表１
に示すものを用いたが、これに限定されるものではなく
以下に示す表２、表３に示すものを用いてもよい。
【００６０】
【表２】
【００６１】
【表３】
【００６２】また、これ以外にもウェーブレット変換を
行うことのできる関数であれば、いかなる関数を用いて
もよく、例えば双直交ではなく対称ではないが直交する
ものを用いてもよい。
行うことのできる関数であれば、いかなる関数を用いて
もよく、例えば双直交ではなく対称ではないが直交する
ものを用いてもよい。
【００６３】さらに、表１，２および３に示すようにｎ
＝０の軸に関して左右対称な関数のみではなく、ｎ＝０
の軸に関して左右非対称な関数を用いてウェーブレット
変換を行うようにしてもよいものである。このように左
右非対称な関数を用いてウェーブレット変換を行った場
合は、ウェーブレット変換を行った関数をｎ＝０の軸に
関して左右を反転させた関数を用いて逆ウェーブレット
変換を行うものである。すなわち、左右非対称な関数
ｇ，ｈについて、逆ウェーブレット変換を行う関数
ｇ′，ｈ′は、 ｇ［ｎ］＝ｇ′［−ｎ］ ｈ［ｎ］＝ｈ′［−ｎ］ …(3) 但し、［−ｎ］は左右反転を表す。
＝０の軸に関して左右対称な関数のみではなく、ｎ＝０
の軸に関して左右非対称な関数を用いてウェーブレット
変換を行うようにしてもよいものである。このように左
右非対称な関数を用いてウェーブレット変換を行った場
合は、ウェーブレット変換を行った関数をｎ＝０の軸に
関して左右を反転させた関数を用いて逆ウェーブレット
変換を行うものである。すなわち、左右非対称な関数
ｇ，ｈについて、逆ウェーブレット変換を行う関数
ｇ′，ｈ′は、 ｇ［ｎ］＝ｇ′［−ｎ］ ｈ［ｎ］＝ｈ′［−ｎ］ …(3) 但し、［−ｎ］は左右反転を表す。
【００６４】となる。
【００６５】さらに、上述した実施の形態においては、
ウェーブレット変換により、画像を多重解像度画像に変
換するようにしているが、これに限定されるものではな
く、上述したラプラシアンピラミッドの手法やフーリエ
変換により画像を多重解像度に変換するようにしてもよ
いものである。
ウェーブレット変換により、画像を多重解像度画像に変
換するようにしているが、これに限定されるものではな
く、上述したラプラシアンピラミッドの手法やフーリエ
変換により画像を多重解像度に変換するようにしてもよ
いものである。
【図１】本発明による画像処理方法の基本的概念を表す
図
図
【図２】本発明に用いられる画像信号の読み取り方式を
表す図
表す図
【図３】ウェーブレット変換の詳細を表す図
【図４】ウェーブレット変換係数信号を表す図
【図５】ウェーブレット変換係数信号に対して閾値処理
を施すテーブルを表す図
を施すテーブルを表す図
【図６】逆ウェーブレット変換の詳細を表す図
【図７】原画像データのコントラストを低減させるテー
ブルを表す図
ブルを表す図
【図８】原画像データのコントラストを元に戻すテーブ
ルを表す図
ルを表す図
【図９】ウェーブレット変換に用いられる基本ウェーブ
レット関数を表す図
レット関数を表す図
【図１０】ウェーブレット変換を説明するための図
【図１１】フーリエ変換を説明するための図
10 蓄積性蛍光体シート ｈ，ｈ′，ｇ，ｇ′ ウェーブレット変換を行うため
の関数 ＶＶi ，ＶＷi ，ＷＶi ，ＷＷi （ｉ＝１〜ｎ） ウェーブレット変換係数信号
の関数 ＶＶi ，ＶＷi ，ＷＶi ，ＷＷi （ｉ＝１〜ｎ） ウェーブレット変換係数信号
Claims (3)
- 【請求項１】 放射線画像を表す画像信号を多重解像度
空間に変換することにより、該画像信号を複数の周波数
帯域毎の画像を表す周波数帯域画像信号に分解し、 該複数の周波数帯域画像信号のうち所定の周波数帯域画
像信号に対して、該各信号値の絶対値が所定の閾値以下
の信号値を０とする処理を施し、 該処理が施された周波数帯域画像信号および他の周波数
帯域画像信号を逆変換することにより処理済画像信号を
得ることを特徴とする画像処理方法。 - 【請求項２】 前記周波数帯域画像信号のうち比較的中
間周波数帯域の周波数帯域画像信号について、前記所定
の閾値を他の周波数帯域の閾値よりも大きくすることを
特徴とする請求項１記載の画像処理方法。 - 【請求項３】 前記周波数帯域画像信号の比較的低濃度
部分の前記閾値を、前記周波数帯域画像信号の比較的高
濃度部の前記閾値よりも大きくすることを特徴とする請
求項１または２記載の画像処理方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP8014512A JPH09212623A (ja) | 1996-01-30 | 1996-01-30 | 画像処理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP8014512A JPH09212623A (ja) | 1996-01-30 | 1996-01-30 | 画像処理方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH09212623A true JPH09212623A (ja) | 1997-08-15 |
Family
ID=11863145
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP8014512A Withdrawn JPH09212623A (ja) | 1996-01-30 | 1996-01-30 | 画像処理方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH09212623A (ja) |
Cited By (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100335753B1 (ko) * | 1999-09-22 | 2002-05-08 | 이수봉 | 웨이블렛 변환을 이용한 차량 검지방법 |
| JP2002133399A (ja) * | 2000-10-23 | 2002-05-10 | Hitachi Medical Corp | 画像処理装置及びそれを用いたｘ線ｃｔ装置 |
| WO2002045019A1 (en) * | 2000-11-30 | 2002-06-06 | Canon Kabushiki Kaisha | Image processing device, image processing method, recording medium and program |
| WO2002045020A1 (en) * | 2000-11-30 | 2002-06-06 | Canon Kabushiki Kaisha | Image processing device, image processing method, storage medium, and program |
| JP2004127064A (ja) * | 2002-10-04 | 2004-04-22 | Konica Minolta Holdings Inc | 画像処理方法、画像処理装置、画像処理プログラム及び画像記録装置 |
| US6771793B1 (en) | 1999-02-17 | 2004-08-03 | Fuji Photo Film Co., Ltd. | Image processing method and apparatus |
| US6985632B2 (en) | 2000-04-17 | 2006-01-10 | Canon Kabushiki Kaisha | Image processing system, image processing apparatus, and image processing method |
| JP2007141176A (ja) * | 2005-11-22 | 2007-06-07 | Iix Inc | ノイズ低減回路及びノイズ低減方法 |
| US7308152B2 (en) | 2001-10-26 | 2007-12-11 | Konica Minolta Holdings, Inc. | Image processing method, apparatus, and program |
| JP2008541672A (ja) * | 2005-05-20 | 2008-11-20 | 株式会社エヌ・ティ・ティ・ドコモ | ビデオ符号化におけるノイズフィルタリング方法及び装置 |
| JP2010020523A (ja) * | 2008-07-10 | 2010-01-28 | Olympus Imaging Corp | 画像処理装置及び撮像システム |
| KR100965880B1 (ko) * | 2003-09-29 | 2010-06-24 | 삼성전자주식회사 | 영상의 잡음 제거 방법 |
| US8270758B2 (en) | 2009-03-31 | 2012-09-18 | Olympus Corporation | Image processing device and method, program recording medium, and imaging system |
| US8335376B2 (en) | 2008-03-04 | 2012-12-18 | Olympus Corporation | Image processing apparatus, image processing method, and computer-readable storage medium storing image processing program |
| US8391635B2 (en) | 2007-08-31 | 2013-03-05 | Olympus Corporation | Noise removal device, noise removal method, and computer readable recording medium |
| JP2013074306A (ja) * | 2011-09-26 | 2013-04-22 | Canon Inc | 画像処理装置および画像処理方法 |
| JP2021527977A (ja) * | 2018-06-28 | 2021-10-14 | アップル インコーポレイテッドＡｐｐｌｅ Ｉｎｃ． | 動画符号化システム |
| US11451788B2 (en) | 2018-06-28 | 2022-09-20 | Apple Inc. | Rate control for low latency video encoding and transmission |
| US11496758B2 (en) | 2018-06-28 | 2022-11-08 | Apple Inc. | Priority-based video encoding and transmission |
-
1996
- 1996-01-30 JP JP8014512A patent/JPH09212623A/ja not_active Withdrawn
Cited By (27)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6771793B1 (en) | 1999-02-17 | 2004-08-03 | Fuji Photo Film Co., Ltd. | Image processing method and apparatus |
| KR100335753B1 (ko) * | 1999-09-22 | 2002-05-08 | 이수봉 | 웨이블렛 변환을 이용한 차량 검지방법 |
| US6985632B2 (en) | 2000-04-17 | 2006-01-10 | Canon Kabushiki Kaisha | Image processing system, image processing apparatus, and image processing method |
| JP2002133399A (ja) * | 2000-10-23 | 2002-05-10 | Hitachi Medical Corp | 画像処理装置及びそれを用いたｘ線ｃｔ装置 |
| WO2002045019A1 (en) * | 2000-11-30 | 2002-06-06 | Canon Kabushiki Kaisha | Image processing device, image processing method, recording medium and program |
| WO2002045020A1 (en) * | 2000-11-30 | 2002-06-06 | Canon Kabushiki Kaisha | Image processing device, 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 |
| 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 |
| US7558434B2 (en) | 2000-11-30 | 2009-07-07 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method, storage medium, and program |
| US7248748B2 (en) | 2000-11-30 | 2007-07-24 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method, storage medium, and program |
| US7447376B2 (en) | 2000-11-30 | 2008-11-04 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method, storage medium, and program |
| US7308152B2 (en) | 2001-10-26 | 2007-12-11 | Konica Minolta Holdings, Inc. | Image processing method, apparatus, and program |
| JP2004127064A (ja) * | 2002-10-04 | 2004-04-22 | Konica Minolta Holdings Inc | 画像処理方法、画像処理装置、画像処理プログラム及び画像記録装置 |
| US7260266B2 (en) | 2002-10-04 | 2007-08-21 | Konica Corporation | Image processing method, image processing apparatus, image processing program and image recording apparatus |
| KR100965880B1 (ko) * | 2003-09-29 | 2010-06-24 | 삼성전자주식회사 | 영상의 잡음 제거 방법 |
| JP2008541672A (ja) * | 2005-05-20 | 2008-11-20 | 株式会社エヌ・ティ・ティ・ドコモ | ビデオ符号化におけるノイズフィルタリング方法及び装置 |
| JP2007141176A (ja) * | 2005-11-22 | 2007-06-07 | Iix Inc | ノイズ低減回路及びノイズ低減方法 |
| US8391635B2 (en) | 2007-08-31 | 2013-03-05 | Olympus Corporation | Noise removal device, noise removal method, and computer readable recording medium |
| US8335376B2 (en) | 2008-03-04 | 2012-12-18 | Olympus Corporation | Image processing apparatus, image processing method, and computer-readable storage medium storing image processing program |
| JP2010020523A (ja) * | 2008-07-10 | 2010-01-28 | Olympus Imaging Corp | 画像処理装置及び撮像システム |
| US8270758B2 (en) | 2009-03-31 | 2012-09-18 | Olympus Corporation | Image processing device and method, program recording medium, and imaging system |
| JP2013074306A (ja) * | 2011-09-26 | 2013-04-22 | Canon Inc | 画像処理装置および画像処理方法 |
| JP2021527977A (ja) * | 2018-06-28 | 2021-10-14 | アップル インコーポレイテッドＡｐｐｌｅ Ｉｎｃ． | 動画符号化システム |
| US11451788B2 (en) | 2018-06-28 | 2022-09-20 | Apple Inc. | Rate control for low latency video encoding and transmission |
| US11496758B2 (en) | 2018-06-28 | 2022-11-08 | Apple Inc. | Priority-based video encoding and transmission |
| US11653026B2 (en) | 2018-06-28 | 2023-05-16 | Apple Inc. | Video encoding system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH09212623A (ja) | 画像処理方法 | |
| JP4004562B2 (ja) | 画像処理方法および装置 | |
| US6771793B1 (en) | Image processing method and apparatus | |
| US6389176B1 (en) | System, method and medium for increasing compression of an image while minimizing image degradation | |
| JP2001057677A (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP2002522831A (ja) | Ｌｃｄパネルにおける画像表示に適したｄｗｔに基づくアップサンプリング・アルゴリズム | |
| JP4244094B2 (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP5027757B2 (ja) | 動画像雑音除去装置、その方法およびそのプログラム | |
| JP3675896B2 (ja) | 画像処理方法および装置 | |
| JP2000224421A (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP2001167264A (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP3700798B2 (ja) | 画像処理方法および装置 | |
| JP3738791B2 (ja) | 画像処理方法および装置 | |
| JP4004561B2 (ja) | 画像処理方法および装置 | |
| JPH1031737A (ja) | 画像処理方法 | |
| JP2002074356A (ja) | 画像処理方法および装置並びに記録媒体 | |
| JPH0944651A (ja) | 画像処理方法および装置 | |
| US8180169B2 (en) | System and method for multi-scale sigma filtering using quadrature mirror filters | |
| JPH09214967A (ja) | 画像データ圧縮処理方法 | |
| JP2000307868A (ja) | 画像処理方法および装置 | |
| JP2003337942A (ja) | 画像処理装置 | |
| JPH09307897A (ja) | 画像データ圧縮処理方法 | |
| JPH09161061A (ja) | 画像処理方法および装置 | |
| JP2001056856A (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP5022400B2 (ja) | 動画像雑音除去装置、動領域画像雑音除去装置、動画像雑音除去プログラム、及び動領域画像雑音除去プログラム |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A300 | Application deemed to be withdrawn because no request for examination was validly filed |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20030401 |