JP2001056856A - 画像処理方法および装置並びに記録媒体 - Google Patents
画像処理方法および装置並びに記録媒体Info
- Publication number
- JP2001056856A JP2001056856A JP2000143166A JP2000143166A JP2001056856A JP 2001056856 A JP2001056856 A JP 2001056856A JP 2000143166 A JP2000143166 A JP 2000143166A JP 2000143166 A JP2000143166 A JP 2000143166A JP 2001056856 A JP2001056856 A JP 2001056856A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- band
- image signal
- image
- frequency
- 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
- 238000000034 method Methods 0.000 title claims abstract description 80
- 238000012545 processing Methods 0.000 title claims description 175
- 238000006243 chemical reaction Methods 0.000 claims abstract description 104
- 230000008569 process Effects 0.000 claims abstract description 53
- 238000003672 processing method Methods 0.000 claims description 26
- 238000007906 compression Methods 0.000 claims description 25
- 230000006835 compression Effects 0.000 claims description 21
- 238000012886 linear function Methods 0.000 claims description 9
- 238000000354 decomposition reaction Methods 0.000 claims description 7
- 230000004044 response Effects 0.000 abstract description 11
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 29
- 238000001914 filtration Methods 0.000 description 23
- 238000004364 calculation method Methods 0.000 description 22
- 101000920618 Homo sapiens Transcription and mRNA export factor ENY2 Proteins 0.000 description 15
- 102100031954 Transcription and mRNA export factor ENY2 Human genes 0.000 description 15
- 230000008859 change Effects 0.000 description 5
- 230000005855 radiation Effects 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 2
- 238000004904 shortening Methods 0.000 description 2
- 102100033040 Carbonic anhydrase 12 Human genes 0.000 description 1
- UNPLRYRWJLTVAE-UHFFFAOYSA-N Cloperastine hydrochloride Chemical compound Cl.C1=CC(Cl)=CC=C1C(C=1C=CC=CC=1)OCCN1CCCCC1 UNPLRYRWJLTVAE-UHFFFAOYSA-N 0.000 description 1
- 101000867855 Homo sapiens Carbonic anhydrase 12 Proteins 0.000 description 1
- 101100339395 Hydra viridissima Hln-1 gene Proteins 0.000 description 1
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000012958 reprocessing Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Editing Of Facsimile Originals (AREA)
- Facsimile Image Signal Circuits (AREA)
- Compression Of Band Width Or Redundancy In Fax (AREA)
- Compression, Expansion, Code Conversion, And Decoders (AREA)
- Image Processing (AREA)
- Picture Signal Circuits (AREA)
Abstract
(57)【要約】
【課題】 原画像信号から帯域制限画像信号を求めて、
周波数処理等の変換処理を行う際に、演算時間を短縮す
る。 【解決手段】 帯域制限画像信号作成手段２において、
原画像信号Ｓorgの各周波数帯域毎の周波数応答特性を
表すとともに、周波数帯域に応じた画素数を有する帯域
制限画像信号Ｂk（ｋ＝１〜ｎ）を作成する。各帯域制
限画像信号Ｂkを変換器２２において変換関数ｆkにより
変換し、変換済み信号ｆkＢkを１段階高周波数帯域の帯
域制限画像信号Ｂk-1と同一の画素数を有するように補
間拡大して拡大高周波信号Ｓｋ′を得る。拡大高周波信
号Ｓｋ′と変換済み信号ｆk-1Ｂk -1とを加算して高周波
信号Ｓｋ−１を得る。この高周波信号Ｓｋ−１の補間拡
大による拡大高周波信号Ｓｋ−１′の取得、および拡大
高周波信号Ｓｋ−１′と変換済み信号ｆk-1Ｂk-1との加
算による高周波信号Ｓｋ−２の取得を繰り返し行って、
最高解像度の高周波信号Ｓ１を得る。この高周波信号Ｓ
１を用いて周波数強調処理を行って処理済み画像信号Ｓ
procを得る。
周波数処理等の変換処理を行う際に、演算時間を短縮す
る。 【解決手段】 帯域制限画像信号作成手段２において、
原画像信号Ｓorgの各周波数帯域毎の周波数応答特性を
表すとともに、周波数帯域に応じた画素数を有する帯域
制限画像信号Ｂk（ｋ＝１〜ｎ）を作成する。各帯域制
限画像信号Ｂkを変換器２２において変換関数ｆkにより
変換し、変換済み信号ｆkＢkを１段階高周波数帯域の帯
域制限画像信号Ｂk-1と同一の画素数を有するように補
間拡大して拡大高周波信号Ｓｋ′を得る。拡大高周波信
号Ｓｋ′と変換済み信号ｆk-1Ｂk -1とを加算して高周波
信号Ｓｋ−１を得る。この高周波信号Ｓｋ−１の補間拡
大による拡大高周波信号Ｓｋ−１′の取得、および拡大
高周波信号Ｓｋ−１′と変換済み信号ｆk-1Ｂk-1との加
算による高周波信号Ｓｋ−２の取得を繰り返し行って、
最高解像度の高周波信号Ｓ１を得る。この高周波信号Ｓ
１を用いて周波数強調処理を行って処理済み画像信号Ｓ
procを得る。
Description
【０００１】
【発明の属する技術分野】本発明は、画像信号に対して
所定の周波数成分を強調する等の画像処理を行うための
画像処理方法および装置並びに画像処理方法をコンピュ
ータに実行させるためのプログラムを記録したコンピュ
ータ読取り可能な記録媒体に関するものである。
所定の周波数成分を強調する等の画像処理を行うための
画像処理方法および装置並びに画像処理方法をコンピュ
ータに実行させるためのプログラムを記録したコンピュ
ータ読取り可能な記録媒体に関するものである。
【０００２】
【従来の技術】画像処理の分野においては、画像信号に
対して各周波数帯域毎に異なる画像処理を施すような場
合に、画像信号を周波数帯域毎に分類する手段としてウ
ェーブレット変換やラプラシアンピラミッドの手法が用
いられている。ここで画像処理としては、例えばノイズ
除去のための高周波の分離、さらにはノイズの多い周波
数帯域のデータを削減することによる圧縮処理などが挙
げられる。本出願人も、ウェーブレット変換を用いて画
像中のエッジ成分のみを強調する等の画像処理を行うよ
うにした画像処理方法を種々提案している（例えば特開
平6-274615号、同6-350989号等）。
対して各周波数帯域毎に異なる画像処理を施すような場
合に、画像信号を周波数帯域毎に分類する手段としてウ
ェーブレット変換やラプラシアンピラミッドの手法が用
いられている。ここで画像処理としては、例えばノイズ
除去のための高周波の分離、さらにはノイズの多い周波
数帯域のデータを削減することによる圧縮処理などが挙
げられる。本出願人も、ウェーブレット変換を用いて画
像中のエッジ成分のみを強調する等の画像処理を行うよ
うにした画像処理方法を種々提案している（例えば特開
平6-274615号、同6-350989号等）。
【０００３】一方、ラプラシアンピラミッドなる方法は
例えば特開平5-244508号、特開平6-96200 、特開平6-30
1766号等に記載されており、このラプラシアンピラミッ
ドは、原画像に対してガウス関数で近似されたようなマ
スクによりマスク処理を施した後、画像をサブサンプリ
ングして画素数を間引いて半分にすることにより、原画
像の１／４のサイズのボケ画像を得、このボケ画像のサ
ンプリングされた画素に値が０の画素を補間して元の大
きさの画像に戻し、この画像に対してさらに上述したマ
スクによりマスク処理を施してボケ画像を得、このボケ
画像を原画像から減算して、原画像信号のある限られた
周波数帯域の周波数成分を表す、すなわち原画像の複数
の周波数帯域毎の周波数応答特性を表す帯域制限画像信
号を得るものである。この処理を得られたボケ画像に対
して繰り返すことにより原画像の１／２2Nのサイズを有
する帯域制限画像信号をＮ個作成するものである。な
お、最低周波数帯域のボケ画像は原画像の低周波成分を
表すものとなる。
例えば特開平5-244508号、特開平6-96200 、特開平6-30
1766号等に記載されており、このラプラシアンピラミッ
ドは、原画像に対してガウス関数で近似されたようなマ
スクによりマスク処理を施した後、画像をサブサンプリ
ングして画素数を間引いて半分にすることにより、原画
像の１／４のサイズのボケ画像を得、このボケ画像のサ
ンプリングされた画素に値が０の画素を補間して元の大
きさの画像に戻し、この画像に対してさらに上述したマ
スクによりマスク処理を施してボケ画像を得、このボケ
画像を原画像から減算して、原画像信号のある限られた
周波数帯域の周波数成分を表す、すなわち原画像の複数
の周波数帯域毎の周波数応答特性を表す帯域制限画像信
号を得るものである。この処理を得られたボケ画像に対
して繰り返すことにより原画像の１／２2Nのサイズを有
する帯域制限画像信号をＮ個作成するものである。な
お、最低周波数帯域のボケ画像は原画像の低周波成分を
表すものとなる。
【０００４】ここで、上記特開平5-244508号において
は、放射線画像をラプラシアンピラミッドにより複数の
周波数帯域の画像に分解して各周波数帯域毎の帯域制限
画像信号を得、各帯域制限画像信号を非線形関数により
変換し、変換後の帯域制限画像信号を最低周波数帯域の
ボケ画像を表すボケ画像信号とともに再構成することに
より、画像中の周波数帯域に応じたコントラストが強調
された処理済み画像信号を得るようにした方法が提案さ
れている。
は、放射線画像をラプラシアンピラミッドにより複数の
周波数帯域の画像に分解して各周波数帯域毎の帯域制限
画像信号を得、各帯域制限画像信号を非線形関数により
変換し、変換後の帯域制限画像信号を最低周波数帯域の
ボケ画像を表すボケ画像信号とともに再構成することに
より、画像中の周波数帯域に応じたコントラストが強調
された処理済み画像信号を得るようにした方法が提案さ
れている。
【０００５】一方、本出願人により、非鮮鋭マスク画像
信号（以下、ボケ画像信号という）を用いて周波数強調
処理あるいはダイナミックレンジ圧縮処理等を行って放
射線画像の診断性能を向上させる数々の画像処理方法お
よび装置が提案されている（特開昭55-163472号、同55-
87953号、特開平3-222577号、同10-75395号、同10-7536
4号、同10-171983号等）。例えば周波数強調処理は、原
画像信号Ｓorgからボケ画像信号Ｓusを引いたものに強
調係数βを乗じたものを、原画像信号Ｓorgに加算する
ことにより、原画像信号の所定の空間周波数成分を強調
するものである。これを式で表すと下記の式（１）のよ
うになる。
信号（以下、ボケ画像信号という）を用いて周波数強調
処理あるいはダイナミックレンジ圧縮処理等を行って放
射線画像の診断性能を向上させる数々の画像処理方法お
よび装置が提案されている（特開昭55-163472号、同55-
87953号、特開平3-222577号、同10-75395号、同10-7536
4号、同10-171983号等）。例えば周波数強調処理は、原
画像信号Ｓorgからボケ画像信号Ｓusを引いたものに強
調係数βを乗じたものを、原画像信号Ｓorgに加算する
ことにより、原画像信号の所定の空間周波数成分を強調
するものである。これを式で表すと下記の式（１）のよ
うになる。
【０００６】 Ｓproc＝Ｓorg＋β×（Ｓorg−Ｓus） （１） （Ｓproc：周波数強調処理された信号、Ｓorg：原画像
信号、Ｓus：ボケ画像信号、β：強調係数） また、特開平10-75395号には、原画像信号に加算する加
算信号の周波数応答特性を調整することにより、周波数
強調処理された信号に対してアーチファクトが発生する
ことを防止する方法が提案されている。この方法とは、
まず鮮鋭度の異なる、すなわち周波数応答特性の異なる
複数のボケ画像信号を作成し、そのボケ画像信号および
原画像信号の中の２つの信号の差分をとることにより、
原画像信号の、ある限られた周波数帯域の周波数成分を
表す複数の帯域制限画像信号（すなわち帯域制限画像信
号）を作成し、さらにその帯域制限画像信号をそれぞれ
異なる変換関数によって所望の大きさとなるように変換
してから、その複数の抑制された帯域制限画像信号を積
算することにより上記加算信号を作成するものである。
この処理は例えば下記の式（２）により表すことができ
る。
信号、Ｓus：ボケ画像信号、β：強調係数） また、特開平10-75395号には、原画像信号に加算する加
算信号の周波数応答特性を調整することにより、周波数
強調処理された信号に対してアーチファクトが発生する
ことを防止する方法が提案されている。この方法とは、
まず鮮鋭度の異なる、すなわち周波数応答特性の異なる
複数のボケ画像信号を作成し、そのボケ画像信号および
原画像信号の中の２つの信号の差分をとることにより、
原画像信号の、ある限られた周波数帯域の周波数成分を
表す複数の帯域制限画像信号（すなわち帯域制限画像信
号）を作成し、さらにその帯域制限画像信号をそれぞれ
異なる変換関数によって所望の大きさとなるように変換
してから、その複数の抑制された帯域制限画像信号を積
算することにより上記加算信号を作成するものである。
この処理は例えば下記の式（２）により表すことができ
る。
【０００７】 Ｓproc＝Ｓorg＋β（Ｓorg）×Ｆusm（Ｓorg,Ｓus1,Ｓus2,…Ｓusn） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…Ｓusn） ＝ｆ1(Ｓorg−Ｓus1)＋ｆ2(Ｓus1−Ｓus2)＋… ＋ｆk(Ｓusk-1−Ｓusk)＋…＋ｆn(Ｓusn-1−Ｓusn) （２） （但し、Ｓproc：処理済み画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜n）：ボケ画像信号 ｆk(k=1〜n）：各帯域制限画像信号を変換する変換関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） さらに、特開平10-75364号には、ダイナミックレンジ圧
縮処理を施す場合において、処理が施された信号に対し
てアーチファクトが発生することを防止する方法が提案
されている。この方法とは、上記特開平10-75395号に記
載されたように、複数の帯域制限画像信号を作成し、こ
の帯域制限画像信号に基づいて原画像信号の低周波成分
に関する信号（低周波成分信号）を得、原画像信号にこ
の低周波成分に関する信号を加算することによりダイナ
ミックレンジ圧縮処理を施すようにしたものである。こ
の処理は例えば下記の式（３）により表すことができ
る。
数） さらに、特開平10-75364号には、ダイナミックレンジ圧
縮処理を施す場合において、処理が施された信号に対し
てアーチファクトが発生することを防止する方法が提案
されている。この方法とは、上記特開平10-75395号に記
載されたように、複数の帯域制限画像信号を作成し、こ
の帯域制限画像信号に基づいて原画像信号の低周波成分
に関する信号（低周波成分信号）を得、原画像信号にこ
の低周波成分に関する信号を加算することによりダイナ
ミックレンジ圧縮処理を施すようにしたものである。こ
の処理は例えば下記の式（３）により表すことができ
る。
【０００８】 Ｓproc＝Ｓorg ＋Ｄ（Ｓorg−Ｆdrc（Ｓorg,Ｓus1,Ｓus2,…Ｓusn）） （３ ） Ｆdrc（Ｓorg,Ｓus1,Ｓus2,…Ｓusn） ＝｛ｆd1(Ｓorg−Ｓus1)＋ｆd2(Ｓus1 −Ｓus2)＋… ＋ｆdk(Ｓusk-1−Ｓusk )＋…＋ｆdn(Ｓusn-1−Ｓusn)｝ （但し、Ｓproc：処理済み画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜n）：ボケ画像信号 ｆdk(k=1〜n）：低周波成分信号を得るために使用する
変換関数 Ｄ（Ｓorg−Ｆdrc）：低周波成分信号に基づいて定めら
れるダイナミックレンジ圧縮係数（ＤはＳorg−Ｆdrcを
変換する関数）） さらにまた、特開平10-171983号には、周波数強調処理
とダイナミックレンジ圧縮処理とを同時に施す場合にお
いて、処理が施された信号に対してアーチファクトが発
生することを防止する方法が提案されている。この方法
とは、上記特開平10-75395号に記載されたように、複数
の帯域制限画像信号を作成し、この帯域制限画像信号に
基づいて原画像信号の高周波成分に関する信号（高周波
成分信号）および低周波成分信号を得、原画像信号にこ
れらの高周波成分に関する信号および低周波成分に関す
る信号を加算することにより周波数強調処理を施すとと
もにダイナミックレンジ圧縮処理を施すようにしたもの
である。この処理は例えば下記の式（４）により表すこ
とができる。
変換関数 Ｄ（Ｓorg−Ｆdrc）：低周波成分信号に基づいて定めら
れるダイナミックレンジ圧縮係数（ＤはＳorg−Ｆdrcを
変換する関数）） さらにまた、特開平10-171983号には、周波数強調処理
とダイナミックレンジ圧縮処理とを同時に施す場合にお
いて、処理が施された信号に対してアーチファクトが発
生することを防止する方法が提案されている。この方法
とは、上記特開平10-75395号に記載されたように、複数
の帯域制限画像信号を作成し、この帯域制限画像信号に
基づいて原画像信号の高周波成分に関する信号（高周波
成分信号）および低周波成分信号を得、原画像信号にこ
れらの高周波成分に関する信号および低周波成分に関す
る信号を加算することにより周波数強調処理を施すとと
もにダイナミックレンジ圧縮処理を施すようにしたもの
である。この処理は例えば下記の式（４）により表すこ
とができる。
【０００９】 Ｓproc＝Ｓorg ＋β（Ｓorg）・Ｆusm（Ｓorg,Ｓus1,Ｓus2,…Ｓusn） ＋Ｄ（Ｓorg−Ｆdrc（Ｓorg,Ｓus1,Ｓus2,…Ｓusn）） （４ ） Ｆusm（Ｓorg,Ｓus1,Ｓus2,…Ｓusn） ＝｛ｆu1(Ｓorg−Ｓus1)＋ｆu2(Ｓus1 −Ｓus2)＋… ＋ｆuk(Ｓusk-1−Ｓusk )＋…＋ｆun(Ｓusn-1−Ｓusn)｝ Ｆdrc（Ｓorg,Ｓus1,Ｓus2,…Ｓusn） ＝｛ｆd1(Ｓorg−Ｓus1)＋ｆd2(Ｓus1 −Ｓus2)＋… ＋ｆdk(Ｓusk-1−Ｓusk )＋…＋ｆdn(Ｓusn-1−Ｓusn)｝ （但し、Ｓproc：処理済み画像信号 Ｓorg：原画像信号 Ｓusk(k=1〜n）：ボケ画像信号 ｆuk(k=1〜n）：高周波成分信号を得るために使用する
変換関数 ｆdk(k=1〜n）：低周波成分信号を得るために使用する
変換関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数 Ｄ（Ｓorg−Ｆdrc）：低周波成分信号に基づいて定めら
れるダイナミックレンジ圧縮係数（ＤはＳorg−Ｆdrcを
変換する関数）） これらの周波数強調処理やダイナミックレンジ圧縮処理
（以下変換処理とする）においては、帯域制限画像信号
を変換する変換関数等の定義を変更することによって原
画像信号に加算する加算信号の周波数応答特性を調整す
ることができる。このため、各変換関数の定義次第で、
アーチファクトの発生防止等所望の周波数応答特性を有
する処理済み画像信号を得ることができる。
変換関数 ｆdk(k=1〜n）：低周波成分信号を得るために使用する
変換関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数 Ｄ（Ｓorg−Ｆdrc）：低周波成分信号に基づいて定めら
れるダイナミックレンジ圧縮係数（ＤはＳorg−Ｆdrcを
変換する関数）） これらの周波数強調処理やダイナミックレンジ圧縮処理
（以下変換処理とする）においては、帯域制限画像信号
を変換する変換関数等の定義を変更することによって原
画像信号に加算する加算信号の周波数応答特性を調整す
ることができる。このため、各変換関数の定義次第で、
アーチファクトの発生防止等所望の周波数応答特性を有
する処理済み画像信号を得ることができる。
【００１０】一方、上述した変換処理に用いられるボケ
画像信号は、まず原画像信号の画素に対して所定間隔毎
に所定のフィルタリング処理を施すことによって画素を
間引きし、そのようにして得た画像信号に対して同様の
フィルタリング処理を繰り返してさらに画素数を少なく
した低解像度の画像信号を複数作成し、低解像度画像信
号のそれぞれに対して、所定の補間方法により原画像と
画素数が同一となるように補間処理を施すことにより作
成される。したがって、ボケ画像信号は画素数は原画像
信号と同一であるが、原画像信号よりも鮮鋭度が低い画
像を表す画像信号となる。
画像信号は、まず原画像信号の画素に対して所定間隔毎
に所定のフィルタリング処理を施すことによって画素を
間引きし、そのようにして得た画像信号に対して同様の
フィルタリング処理を繰り返してさらに画素数を少なく
した低解像度の画像信号を複数作成し、低解像度画像信
号のそれぞれに対して、所定の補間方法により原画像と
画素数が同一となるように補間処理を施すことにより作
成される。したがって、ボケ画像信号は画素数は原画像
信号と同一であるが、原画像信号よりも鮮鋭度が低い画
像を表す画像信号となる。
【００１１】また、帯域制限画像信号は、例えば隣接す
る周波数帯域のボケ画像信号同士で差分をとって、ある
いは原画像信号と各ボケ画像信号の差分をとって作成さ
れる。したがって、帯域制限画像信号は画素数は原画像
信号と同じであり、原画像信号の周波数帯域毎の周波数
応答特性を表す信号となる。
る周波数帯域のボケ画像信号同士で差分をとって、ある
いは原画像信号と各ボケ画像信号の差分をとって作成さ
れる。したがって、帯域制限画像信号は画素数は原画像
信号と同じであり、原画像信号の周波数帯域毎の周波数
応答特性を表す信号となる。
【００１２】
【発明が解決しようとする課題】上記特開平5-244508号
に記載された方法においては、各周波数帯域の帯域制限
画像信号に対して非線形関数により変換処理を施し、変
換された帯域制限画像信号を再構成して処理済み画像信
号を得ているため、変換の程度を修正したい場合には、
再度帯域制限画像信号に対して変換処理を施した後に画
像を再構成する必要があり、この結果、再処理には長時
間を要するものとなる。このため、ＣＲＴを観察しなが
ら画像処理の内容を種々変更して最適な画像処理内容を
設定するような場合には、処理が施された画像が再生さ
れるまでに長時間を要するものとなり、オペレータのス
トレスが大きい。
に記載された方法においては、各周波数帯域の帯域制限
画像信号に対して非線形関数により変換処理を施し、変
換された帯域制限画像信号を再構成して処理済み画像信
号を得ているため、変換の程度を修正したい場合には、
再度帯域制限画像信号に対して変換処理を施した後に画
像を再構成する必要があり、この結果、再処理には長時
間を要するものとなる。このため、ＣＲＴを観察しなが
ら画像処理の内容を種々変更して最適な画像処理内容を
設定するような場合には、処理が施された画像が再生さ
れるまでに長時間を要するものとなり、オペレータのス
トレスが大きい。
【００１３】本発明は上記事情に鑑みなされたものであ
り、帯域制限画像信号を用いた周波数強調処理等の画像
処理を高速に行うことができる画像処理方法および装置
並びに画像処理方法をコンピュータに実行させるための
プログラムを記録したコンピュータ読取り可能な記録媒
体を提供することを目的とするものである。
り、帯域制限画像信号を用いた周波数強調処理等の画像
処理を高速に行うことができる画像処理方法および装置
並びに画像処理方法をコンピュータに実行させるための
プログラムを記録したコンピュータ読取り可能な記録媒
体を提供することを目的とするものである。
【００１４】
【課題を解決するための手段】本発明による画像処理方
法は、原画像を表す原画像信号に対して、該原画像の高
周波成分に関する信号に基づく画像処理を施して処理済
み画像信号を得る画像処理方法において、前記原画像信
号から帯域制限画像信号を作成し、所定の変換関数に基
づいて、前記帯域制限画像信号から前記原画像の高周波
成分に関する信号を得、該高周波成分に関する信号に基
づいて、前記原画像信号に対して前記画像処理を施すこ
とを特徴とするものである。
法は、原画像を表す原画像信号に対して、該原画像の高
周波成分に関する信号に基づく画像処理を施して処理済
み画像信号を得る画像処理方法において、前記原画像信
号から帯域制限画像信号を作成し、所定の変換関数に基
づいて、前記帯域制限画像信号から前記原画像の高周波
成分に関する信号を得、該高周波成分に関する信号に基
づいて、前記原画像信号に対して前記画像処理を施すこ
とを特徴とするものである。
【００１５】ここで、「画像処理」とは、具体的には、
原画像信号に含まれる特定の周波数成分を強調する周波
数強調処理や、原画像の最高濃度と最低濃度との差すな
わちダイナミックレンジを狭めるように高濃度域もしく
は低濃度域あるいは高濃度域および低濃度域の双方のコ
ントラストを下げるダイナミックレンジ圧縮処理等が挙
げられる。さらには、周波数強調処理とダイナミックレ
ンジ圧縮処理とを同時に行うものが挙げられる。なお、
「所定の変換関数」は、画像処理の内容に応じた関数を
用いることが好ましい。
原画像信号に含まれる特定の周波数成分を強調する周波
数強調処理や、原画像の最高濃度と最低濃度との差すな
わちダイナミックレンジを狭めるように高濃度域もしく
は低濃度域あるいは高濃度域および低濃度域の双方のコ
ントラストを下げるダイナミックレンジ圧縮処理等が挙
げられる。さらには、周波数強調処理とダイナミックレ
ンジ圧縮処理とを同時に行うものが挙げられる。なお、
「所定の変換関数」は、画像処理の内容に応じた関数を
用いることが好ましい。
【００１６】なお、本発明による画像処理方法において
は、前記帯域制限画像信号により表される各画像が、そ
の周波数帯域に応じた画素数を有することが好ましい。
は、前記帯域制限画像信号により表される各画像が、そ
の周波数帯域に応じた画素数を有することが好ましい。
【００１７】また、本発明による画像処理方法において
は、前記高周波成分に関する信号が前記原画像信号と同
一画素数であることが好ましい。
は、前記高周波成分に関する信号が前記原画像信号と同
一画素数であることが好ましい。
【００１８】ここで、「原画像信号と同一画素数」と
は、原画像と高周波成分に関する信号により表される画
像との画像サイズが同一であることをいう。
は、原画像と高周波成分に関する信号により表される画
像との画像サイズが同一であることをいう。
【００１９】さらに、本発明による画像処理方法におい
ては、前記原画像信号を多重解像度変換することによ
り、前記帯域制限画像信号を作成し、該帯域制限画像信
号に対して所定の変換関数に基づいて変換処理を施して
変換帯域制限画像信号を得、該変換帯域制限画像信号を
逆多重解像度変換することにより、前記高周波成分に関
する信号を得るものとすることができる。なお、逆多重
解像度変換は、前記多重解像度変換に対応するものであ
って、この逆多重解像度変換を施すことにより、元の信
号を復元（可逆／非可逆のいずれでもよい）することが
できるものであることはいうまでもない。
ては、前記原画像信号を多重解像度変換することによ
り、前記帯域制限画像信号を作成し、該帯域制限画像信
号に対して所定の変換関数に基づいて変換処理を施して
変換帯域制限画像信号を得、該変換帯域制限画像信号を
逆多重解像度変換することにより、前記高周波成分に関
する信号を得るものとすることができる。なお、逆多重
解像度変換は、前記多重解像度変換に対応するものであ
って、この逆多重解像度変換を施すことにより、元の信
号を復元（可逆／非可逆のいずれでもよい）することが
できるものであることはいうまでもない。
【００２０】ここで、「原画像信号を多重解像度変換す
ることにより帯域制限画像信号を作成」するに際して
は、ラプラシアンピラミッドの手法によるラプラシアン
ピラミッド分解により、あるいはウェーブレット変換に
より原画像信号を複数の周波数帯域毎の周波数応答特性
を表す信号に変換する方法などを用いることができる。
この場合、「逆多重解像度変換」としては、ラプラシア
ンピラミッド分解により帯域制限画像信号を得た場合は
ラプラシアンピラミッド再構成の方法が用いられ、ウェ
ーブレット変換により帯域制限画像信号を得た場合は逆
ウェーブレット変換が用いられるのはいうまでもない。
ることにより帯域制限画像信号を作成」するに際して
は、ラプラシアンピラミッドの手法によるラプラシアン
ピラミッド分解により、あるいはウェーブレット変換に
より原画像信号を複数の周波数帯域毎の周波数応答特性
を表す信号に変換する方法などを用いることができる。
この場合、「逆多重解像度変換」としては、ラプラシア
ンピラミッド分解により帯域制限画像信号を得た場合は
ラプラシアンピラミッド再構成の方法が用いられ、ウェ
ーブレット変換により帯域制限画像信号を得た場合は逆
ウェーブレット変換が用いられるのはいうまでもない。
【００２１】なお、ラプラシアンピラミッド分解やウェ
ーブレット変換などにより原画像信号を多重解像度変換
した場合、最低周波数帯域の信号は原画像を縮小した低
周波の情報を表すものであり、原画像信号の周波数帯域
毎の周波数応答特性を表す帯域制限画像信号ではないた
め、本発明においては処理には使用しないか、または値
を０として使用することが好ましい。
ーブレット変換などにより原画像信号を多重解像度変換
した場合、最低周波数帯域の信号は原画像を縮小した低
周波の情報を表すものであり、原画像信号の周波数帯域
毎の周波数応答特性を表す帯域制限画像信号ではないた
め、本発明においては処理には使用しないか、または値
を０として使用することが好ましい。
【００２２】なお、本発明による画像処理方法において
は、前記所定の変換関数は、非線形関数であることが好
ましい。
は、前記所定の変換関数は、非線形関数であることが好
ましい。
【００２３】本発明による画像処理装置は、原画像を表
す原画像信号に対して、該原画像の高周波成分に関する
信号に基づく画像処理を施して処理済み画像信号を得る
画像処理装置において、前記原画像信号から帯域制限画
像信号を作成する帯域制限画像信号作成手段と、所定の
変換関数に基づいて、前記帯域制限画像信号から前記原
画像の高周波成分に関する信号を得る高周波成分取得手
段と、該高周波成分に関する信号に基づいて、前記原画
像信号に対して前記画像処理を施す画像処理手段とを備
えたことを特徴とするものである。
す原画像信号に対して、該原画像の高周波成分に関する
信号に基づく画像処理を施して処理済み画像信号を得る
画像処理装置において、前記原画像信号から帯域制限画
像信号を作成する帯域制限画像信号作成手段と、所定の
変換関数に基づいて、前記帯域制限画像信号から前記原
画像の高周波成分に関する信号を得る高周波成分取得手
段と、該高周波成分に関する信号に基づいて、前記原画
像信号に対して前記画像処理を施す画像処理手段とを備
えたことを特徴とするものである。
【００２４】なお、本発明による画像処理装置において
は、前記帯域制限画像信号作成手段は、前記帯域制限画
像信号により表される各周波数帯域の画像が、その周波
数帯域に応じた画素数を有するよう前記帯域制限画像信
号を作成する手段であることが好ましい。
は、前記帯域制限画像信号作成手段は、前記帯域制限画
像信号により表される各周波数帯域の画像が、その周波
数帯域に応じた画素数を有するよう前記帯域制限画像信
号を作成する手段であることが好ましい。
【００２５】また、本発明による画像処理装置において
は、前記高周波成分取得手段は、前記高周波成分に関す
る信号が前記原画像信号と同一画素数となるように前記
高周波成分に関する信号を取得する手段であることが好
ましい。
は、前記高周波成分取得手段は、前記高周波成分に関す
る信号が前記原画像信号と同一画素数となるように前記
高周波成分に関する信号を取得する手段であることが好
ましい。
【００２６】さらに、本発明による画像処理装置におい
ては、前記帯域制限画像信号作成手段は、前記原画像信
号を多重解像度変換することにより、前記帯域制限画像
信号を作成する手段であり、前記高周波成分取得手段
は、前記帯域制限画像信号に対して前記所定の変換関数
に基づいて変換処理を施して変換帯域制限画像信号を
得、該変換帯域制限画像信号を逆多重解像度変換するこ
とにより、前記高周波成分に関する信号を得る手段であ
ることが好ましい。
ては、前記帯域制限画像信号作成手段は、前記原画像信
号を多重解像度変換することにより、前記帯域制限画像
信号を作成する手段であり、前記高周波成分取得手段
は、前記帯域制限画像信号に対して前記所定の変換関数
に基づいて変換処理を施して変換帯域制限画像信号を
得、該変換帯域制限画像信号を逆多重解像度変換するこ
とにより、前記高周波成分に関する信号を得る手段であ
ることが好ましい。
【００２７】この場合、前記多重解像度変換は、ラプラ
シアンピラミッド分解による変換、またはウェーブレッ
ト変換とすることができる。
シアンピラミッド分解による変換、またはウェーブレッ
ト変換とすることができる。
【００２８】また、本発明による画像処理装置において
は、前記高周波成分取得手段は、最低周波数帯域の帯域
制限画像信号以外の帯域制限画像信号から前記高周波成
分に関する信号を得る手段であることが好ましい。
は、前記高周波成分取得手段は、最低周波数帯域の帯域
制限画像信号以外の帯域制限画像信号から前記高周波成
分に関する信号を得る手段であることが好ましい。
【００２９】さらに、前記所定の変換関数は、非線形関
数であることが好ましい。
数であることが好ましい。
【００３０】さらにまた、前記画像処理は、周波数強調
処理および／またはダイナミックレンジ圧縮処理である
ことが好ましい。
処理および／またはダイナミックレンジ圧縮処理である
ことが好ましい。
【００３１】なお、本発明による画像処理方法をコンピ
ュータに実行させるためのプログラムとして、コンピュ
ータ読取り可能な記録媒体に記録して提供してもよい。
ュータに実行させるためのプログラムとして、コンピュ
ータ読取り可能な記録媒体に記録して提供してもよい。
【００３２】
【発明の効果】本発明によれば、原画像から帯域制限画
像信号を作成し、所定の変換関数に基づいて、帯域制限
画像信号から原画像の高周波成分に関する信号を得、こ
の信号に基づいて原画像信号に対して画像処理を施すよ
うにしたものである。したがって、高周波成分に関する
信号のレベルを変更するのみで原画像信号に対する画像
処理の程度を任意に変更することができるため、上記特
開平5-244508号に記載されたように非線形関数を修正す
ることにより画像処理の程度を変更する方法と比較し
て、画像処理の程度を簡易に変更することができ、これ
により、処理済み画像信号を得るための演算時間を短縮
して、オペレータのストレスを低減することができる。
像信号を作成し、所定の変換関数に基づいて、帯域制限
画像信号から原画像の高周波成分に関する信号を得、こ
の信号に基づいて原画像信号に対して画像処理を施すよ
うにしたものである。したがって、高周波成分に関する
信号のレベルを変更するのみで原画像信号に対する画像
処理の程度を任意に変更することができるため、上記特
開平5-244508号に記載されたように非線形関数を修正す
ることにより画像処理の程度を変更する方法と比較し
て、画像処理の程度を簡易に変更することができ、これ
により、処理済み画像信号を得るための演算時間を短縮
して、オペレータのストレスを低減することができる。
【００３３】また、帯域制限画像信号の各周波数帯域の
画像を、その周波数帯域に応じた画素数を有するものと
することにより、変換処理を行う際の演算量を低減する
ことができ、これにより、処理済み画像信号を得るため
の演算時間を一層短縮することができる。
画像を、その周波数帯域に応じた画素数を有するものと
することにより、変換処理を行う際の演算量を低減する
ことができ、これにより、処理済み画像信号を得るため
の演算時間を一層短縮することができる。
【００３４】さらに、高周波成分に関する信号の画素数
を原画像信号と同一の画素数とすることにより、高周波
成分に関する信号に対して、サイズを変更する処理を施
すことなく、直ちに原画像信号に対して画像処理を施す
ことができるため、画像処理をより効率よく行うことが
できる。
を原画像信号と同一の画素数とすることにより、高周波
成分に関する信号に対して、サイズを変更する処理を施
すことなく、直ちに原画像信号に対して画像処理を施す
ことができるため、画像処理をより効率よく行うことが
できる。
【００３５】
【発明の実施の形態】以下、本発明による画像処理方法
および装置の一実施形態について、図面を参照して詳細
に説明する。以下に示す画像処理装置は、蓄積性蛍光体
シートに記録された人体の放射線画像を読み取って得た
画像信号に対して、その画像が診断に適した画像となる
ように、ボケ画像信号を使用して周波数強調処理を施す
ものであり、処理された画像信号は主としてフィルムに
記録され、診断に用いられる。
および装置の一実施形態について、図面を参照して詳細
に説明する。以下に示す画像処理装置は、蓄積性蛍光体
シートに記録された人体の放射線画像を読み取って得た
画像信号に対して、その画像が診断に適した画像となる
ように、ボケ画像信号を使用して周波数強調処理を施す
ものであり、処理された画像信号は主としてフィルムに
記録され、診断に用いられる。
【００３６】図１は本発明の第１の実施形態による画像
処理装置の構成を示す概略ブロック図である。画像処理
装置１は、読取装置等において得られた所定の解像度を
有する原画像信号Ｓorgから原画像の複数の周波数帯域
毎の周波数応答特性を表す帯域制限画像信号を作成する
帯域制限画像信号作成手段２と、帯域制限画像信号に基
づいて原画像信号Ｓorgに対して特定の周波数を強調す
るための周波数強調処理を行って処理済み画像信号Ｓpr
ocを得る変換処理手段３とを有する。
処理装置の構成を示す概略ブロック図である。画像処理
装置１は、読取装置等において得られた所定の解像度を
有する原画像信号Ｓorgから原画像の複数の周波数帯域
毎の周波数応答特性を表す帯域制限画像信号を作成する
帯域制限画像信号作成手段２と、帯域制限画像信号に基
づいて原画像信号Ｓorgに対して特定の周波数を強調す
るための周波数強調処理を行って処理済み画像信号Ｓpr
ocを得る変換処理手段３とを有する。
【００３７】まず帯域制限画像信号の作成処理について
詳細に説明する。図２は帯域制限画像信号作成処理の概
要を示すブロック図、図３は帯域制限画像信号作成処理
を模式的に示す図である。なお、本実施形態において
は、例えば特開平5-244508号に記載されたラプラシアン
ピラミッドの手法により帯域制限画像信号を作成するも
のとする。まず図２に示すように、図１の帯域制限画像
信号作成手段２は、フィルタリング処理手段１０におい
て原画像信号Ｓorgに対し、原画像のｘ方向およびｙ方
向（図２２参照）に対してフィルタリング処理を施して
原画像信号Ｓorgよりも解像度が低い画像信号Ｌ1（以
下、低解像度画像信号という）を作成し、次にこの低解
像度画像信号Ｌ1に対して同様のフィルタリング処理を
施してこの低解像度画像信号Ｌ1よりもさらに解像度が
低い低解像度画像信号Ｌ2を作成し、以降順次同様のフ
ィルタリング処理を繰り返して各解像度毎の低解像度画
像信号Ｌk（ｋ＝１〜ｎ）を得るものである。そして、
補間処理手段１１において、このフィルタリング処理の
各段において得られる低解像度画像信号Ｌkに対して、
それぞれ２倍の画素数となるように補間処理を施して、
鮮鋭度が異なる複数のボケ画像信号Ｓus1〜Ｓusｎ（以
下Ｓusk（ｋ＝１〜ｎ）で代表させる）を得る。この
後、減算器１２により互いに対応する画素数を有する低
解像度画像信号Ｌk-1とボケ画像信号Ｓuskおよび原画像
信号Ｓorgとボケ画像信号Ｓus１との差分を求め、これ
を帯域制限画像信号Ｂkとする。
詳細に説明する。図２は帯域制限画像信号作成処理の概
要を示すブロック図、図３は帯域制限画像信号作成処理
を模式的に示す図である。なお、本実施形態において
は、例えば特開平5-244508号に記載されたラプラシアン
ピラミッドの手法により帯域制限画像信号を作成するも
のとする。まず図２に示すように、図１の帯域制限画像
信号作成手段２は、フィルタリング処理手段１０におい
て原画像信号Ｓorgに対し、原画像のｘ方向およびｙ方
向（図２２参照）に対してフィルタリング処理を施して
原画像信号Ｓorgよりも解像度が低い画像信号Ｌ1（以
下、低解像度画像信号という）を作成し、次にこの低解
像度画像信号Ｌ1に対して同様のフィルタリング処理を
施してこの低解像度画像信号Ｌ1よりもさらに解像度が
低い低解像度画像信号Ｌ2を作成し、以降順次同様のフ
ィルタリング処理を繰り返して各解像度毎の低解像度画
像信号Ｌk（ｋ＝１〜ｎ）を得るものである。そして、
補間処理手段１１において、このフィルタリング処理の
各段において得られる低解像度画像信号Ｌkに対して、
それぞれ２倍の画素数となるように補間処理を施して、
鮮鋭度が異なる複数のボケ画像信号Ｓus1〜Ｓusｎ（以
下Ｓusk（ｋ＝１〜ｎ）で代表させる）を得る。この
後、減算器１２により互いに対応する画素数を有する低
解像度画像信号Ｌk-1とボケ画像信号Ｓuskおよび原画像
信号Ｓorgとボケ画像信号Ｓus１との差分を求め、これ
を帯域制限画像信号Ｂkとする。
【００３８】本実施形態においては、上記フィルタリン
グ処理のフィルタとして、１次元ガウス分布に略対応し
たフィルタを使用する。すなわちフィルタのフィルタ係
数を、ガウス信号に関する下記の式（５）にしたがって
定める。
グ処理のフィルタとして、１次元ガウス分布に略対応し
たフィルタを使用する。すなわちフィルタのフィルタ係
数を、ガウス信号に関する下記の式（５）にしたがって
定める。
【００３９】
【数１】 これは、ガウス信号は周波数空間および実空間の双方に
おいて局在性がよいためであり、例えば上記（５）式に
おいてσ＝１とした場合の５×１の１次元フィルタは図
４に示すようなものとなる。
おいて局在性がよいためであり、例えば上記（５）式に
おいてσ＝１とした場合の５×１の１次元フィルタは図
４に示すようなものとなる。
【００４０】フィルタリング処理は、図５に示すよう
に、原画像信号Ｓorgに対して、あるいは低解像度画像
信号に対して１画素おきに行う。このような１画素おき
のフィルタリング処理をｘ方向、ｙ方向に行うことによ
り、低解像度画像信号Ｌ1の画素数は原画像の１／４と
なり、フィルタリング処理により得られる低解像度画像
信号に対して繰り返しこのフィルタリング処理を施すこ
とにより、得られるｎ個の低解像度画像信号Ｌk（ｋ＝
１〜ｎ）は、それぞれ画素数が原画像信号Ｓorgの１／
２2kの画像信号となる。
に、原画像信号Ｓorgに対して、あるいは低解像度画像
信号に対して１画素おきに行う。このような１画素おき
のフィルタリング処理をｘ方向、ｙ方向に行うことによ
り、低解像度画像信号Ｌ1の画素数は原画像の１／４と
なり、フィルタリング処理により得られる低解像度画像
信号に対して繰り返しこのフィルタリング処理を施すこ
とにより、得られるｎ個の低解像度画像信号Ｌk（ｋ＝
１〜ｎ）は、それぞれ画素数が原画像信号Ｓorgの１／
２2kの画像信号となる。
【００４１】次に、このようにして得られた低解像度画
像信号Ｌkに対して施される補間処理について説明す
る。補間処理を行うための補間演算の方法としては、Ｂ
スプラインによる方法等種々の方法が挙げられるが、本
実施形態においては、上記フィルタリング処理において
ガウス信号に基づくローパスフィルタを用いているた
め、補間演算についてもガウス信号を用いるものとす
る。具体的には、下記の式（６）において、σ＝２k-1
と近似したものを用いる。
像信号Ｌkに対して施される補間処理について説明す
る。補間処理を行うための補間演算の方法としては、Ｂ
スプラインによる方法等種々の方法が挙げられるが、本
実施形態においては、上記フィルタリング処理において
ガウス信号に基づくローパスフィルタを用いているた
め、補間演算についてもガウス信号を用いるものとす
る。具体的には、下記の式（６）において、σ＝２k-1
と近似したものを用いる。
【００４２】
【数２】 例えば低解像度画像信号Ｌ1を補間する際には、ｋ＝１
であるためσ＝１となる。この場合、補間処理を行うた
めのフィルタは、図６に示すように５×１の１次元フィ
ルタとなる。この補間処理は、まず低解像度画像信号Ｌ
1に対して１画素おきに値が０の画素を１つずつ補間す
ることにより低解像度画像信号Ｌ1を原画像と同一の画
素数となるように拡大し、次に、この補間された低解像
度画像信号Ｌ1に対して上述した図６に示す１次元フィ
ルタによりフィルタリング処理を施すことにより行われ
る。
であるためσ＝１となる。この場合、補間処理を行うた
めのフィルタは、図６に示すように５×１の１次元フィ
ルタとなる。この補間処理は、まず低解像度画像信号Ｌ
1に対して１画素おきに値が０の画素を１つずつ補間す
ることにより低解像度画像信号Ｌ1を原画像と同一の画
素数となるように拡大し、次に、この補間された低解像
度画像信号Ｌ1に対して上述した図６に示す１次元フィ
ルタによりフィルタリング処理を施すことにより行われ
る。
【００４３】同様に、この補間拡大処理を全ての低解像
度画像信号Ｌkに対して行う。低解像度画像信号Ｌkを補
間する際には、上記式（６）に基づいて、３×２k−１
の長さのフィルタを作成し、低解像度画像信号Ｌkの各
画素の間に値が０の画素を１個ずつ補間することによ
り、１段階高解像度の低解像度画像信号Ｌk-1と同一画
素数となるように拡大し、この値が０の画素が補間され
た低解像度画像信号Ｌkに対して３×２k−１の長さのフ
ィルタにより、フィルタリング処理を施すことにより補
間拡大してボケ画像信号Ｓuskを得る。
度画像信号Ｌkに対して行う。低解像度画像信号Ｌkを補
間する際には、上記式（６）に基づいて、３×２k−１
の長さのフィルタを作成し、低解像度画像信号Ｌkの各
画素の間に値が０の画素を１個ずつ補間することによ
り、１段階高解像度の低解像度画像信号Ｌk-1と同一画
素数となるように拡大し、この値が０の画素が補間され
た低解像度画像信号Ｌkに対して３×２k−１の長さのフ
ィルタにより、フィルタリング処理を施すことにより補
間拡大してボケ画像信号Ｓuskを得る。
【００４４】次に、上記のようにして作成されたボケ画
像信号Ｓuskが、対応する画素数を有する低解像度画像
信号Ｌk-1から減算されて、帯域制限画像信号Ｂk（ｋ＝
１〜ｎ）が得られる。なお、帯域制限画像信号Ｂkは下
記の式（７）に示すものとなる。
像信号Ｓuskが、対応する画素数を有する低解像度画像
信号Ｌk-1から減算されて、帯域制限画像信号Ｂk（ｋ＝
１〜ｎ）が得られる。なお、帯域制限画像信号Ｂkは下
記の式（７）に示すものとなる。
【００４５】 Ｂ1＝Ｓorg−Ｓus１ Ｂ2＝Ｌ1−Ｓus２ Ｂ3＝Ｌ2−Ｓus３ ・ （７） ・ Ｂk＝Ｌk−Ｓusｋ 具体的には、図３に示すように５段階の低解像度画像信
号Ｌ1〜Ｌ5が得られた場合、まず最低解像度の低解像度
画像信号Ｌ5に対して補間処理を施して、低解像度画像
信号Ｌ4と同一画素数を有するボケ画像信号Ｓus５を作
成する。そして、低解像度画像信号Ｌ4からボケ画像信
号Ｓus５を減算して帯域制限画像信号Ｂ5を得る。以下
順次Ｌ3−Ｓus４、Ｌ2−Ｓus３、Ｌ1−Ｓus２、Ｓorg−
Ｓus１の演算を行って、帯域制限画像信号Ｂ1〜Ｂ5を得
る。ここで、最低解像度の低解像度画像信号Ｌk（Ｌ5）
は、原画像を縮小した低周波の情報を表すものである
が、これ以降の演算において使用することはない。
号Ｌ1〜Ｌ5が得られた場合、まず最低解像度の低解像度
画像信号Ｌ5に対して補間処理を施して、低解像度画像
信号Ｌ4と同一画素数を有するボケ画像信号Ｓus５を作
成する。そして、低解像度画像信号Ｌ4からボケ画像信
号Ｓus５を減算して帯域制限画像信号Ｂ5を得る。以下
順次Ｌ3−Ｓus４、Ｌ2−Ｓus３、Ｌ1−Ｓus２、Ｓorg−
Ｓus１の演算を行って、帯域制限画像信号Ｂ1〜Ｂ5を得
る。ここで、最低解像度の低解像度画像信号Ｌk（Ｌ5）
は、原画像を縮小した低周波の情報を表すものである
が、これ以降の演算において使用することはない。
【００４６】次に、上述したように算出された帯域制限
画像信号Ｂkを用いて行われる変換処理について説明す
る。本実施形態においては変換処理として周波数強調処
理を行うものとする。図７は変換処理手段３の構成を帯
域制限画像信号作成手段２とともに示す概略ブロック
図、図８は変換処理を模式的に示す図である。図７に示
すように、帯域制限画像信号作成手段２において作成さ
れた帯域制限画像信号Ｂ kが、変換器２２において変換
関数ｆ1〜ｆnにより所望の大きさとなるように抑制され
て、変換済み信号ｆkＢk（ｋ＝１〜ｎ）が得られる。関
数ｆ1の例を図９に示す。この関数ｆ1は、帯域制限画像
信号の絶対値が閾値Ｔｈ１よりも小さい場合は傾きが１
であり、閾値Ｔｈ１よりも大きい場合は傾きが１よりも
小さくなるような非線形関数である。この関数は、各帯
域制限画像信号において同一のものであってもよいが、
各信号毎に異なるものであってもよい。
画像信号Ｂkを用いて行われる変換処理について説明す
る。本実施形態においては変換処理として周波数強調処
理を行うものとする。図７は変換処理手段３の構成を帯
域制限画像信号作成手段２とともに示す概略ブロック
図、図８は変換処理を模式的に示す図である。図７に示
すように、帯域制限画像信号作成手段２において作成さ
れた帯域制限画像信号Ｂ kが、変換器２２において変換
関数ｆ1〜ｆnにより所望の大きさとなるように抑制され
て、変換済み信号ｆkＢk（ｋ＝１〜ｎ）が得られる。関
数ｆ1の例を図９に示す。この関数ｆ1は、帯域制限画像
信号の絶対値が閾値Ｔｈ１よりも小さい場合は傾きが１
であり、閾値Ｔｈ１よりも大きい場合は傾きが１よりも
小さくなるような非線形関数である。この関数は、各帯
域制限画像信号において同一のものであってもよいが、
各信号毎に異なるものであってもよい。
【００４７】そして、この変換済み信号ｆkＢkのうち、
最低解像度の変換済み信号ｆnＢnは高周波信号Ｓｎとさ
れるとともに、高周波信号Ｓｎが１段階高解像度の変換
済み信号ｆn-1Ｂn-1と同一画素数となるように、補間処
理手段２３において上記補間処理手段１１と同様に補間
処理がなされて、拡大高周波信号Ｓｎ′が得られる。こ
の後、変換済み信号ｆn-1Ｂn-1と拡大高周波信号Ｓｎ′
とが加算器２４において加算されて、高周波信号Ｓｎ−
１が得られる。そして、高周波信号Ｓｋ−１の補間拡大
による拡大高周波信号Ｓｋ−１′の取得、および拡大高
周波信号Ｓｋ−１′と変換済み信号ｆk-1Ｂk-1との加算
による高周波信号Ｓｋ−２の取得を繰り返し行って、最
高解像度の高周波信号Ｓ１を得る。
最低解像度の変換済み信号ｆnＢnは高周波信号Ｓｎとさ
れるとともに、高周波信号Ｓｎが１段階高解像度の変換
済み信号ｆn-1Ｂn-1と同一画素数となるように、補間処
理手段２３において上記補間処理手段１１と同様に補間
処理がなされて、拡大高周波信号Ｓｎ′が得られる。こ
の後、変換済み信号ｆn-1Ｂn-1と拡大高周波信号Ｓｎ′
とが加算器２４において加算されて、高周波信号Ｓｎ−
１が得られる。そして、高周波信号Ｓｋ−１の補間拡大
による拡大高周波信号Ｓｋ−１′の取得、および拡大高
周波信号Ｓｋ−１′と変換済み信号ｆk-1Ｂk-1との加算
による高周波信号Ｓｋ−２の取得を繰り返し行って、最
高解像度の高周波信号Ｓ１を得る。
【００４８】具体的には図８に示すように、５段階の帯
域制限画像信号Ｂ1〜Ｂ5が得られた場合には、まず変換
済み信号ｆ1Ｂ1〜ｆ5Ｂ5が得られ、最低解像度の変換済
み信号ｆ5Ｂ5が高周波信号Ｓ５とされる。そして、高周
波信号Ｓ５に対して１段階高解像度の変換済み信号ｆ4
Ｂ4と同一画素数となるように補間処理が施されて拡大
高周波信号Ｓ５′が得られる。そして、変換済み信号ｆ
4Ｂ4と拡大高周波信号Ｓ５′とが加算されて、高周波信
号Ｓ４が得られる。以下同様にして高周波信号Ｓ３，Ｓ
２が得られ、最終的に最高解像度の高周波信号Ｓ１が得
られる。
域制限画像信号Ｂ1〜Ｂ5が得られた場合には、まず変換
済み信号ｆ1Ｂ1〜ｆ5Ｂ5が得られ、最低解像度の変換済
み信号ｆ5Ｂ5が高周波信号Ｓ５とされる。そして、高周
波信号Ｓ５に対して１段階高解像度の変換済み信号ｆ4
Ｂ4と同一画素数となるように補間処理が施されて拡大
高周波信号Ｓ５′が得られる。そして、変換済み信号ｆ
4Ｂ4と拡大高周波信号Ｓ５′とが加算されて、高周波信
号Ｓ４が得られる。以下同様にして高周波信号Ｓ３，Ｓ
２が得られ、最終的に最高解像度の高周波信号Ｓ１が得
られる。
【００４９】なお、変換器２２、補間処理手段２３およ
び加算器２４が本発明による高周波成分取得手段と観念
される。
び加算器２４が本発明による高周波成分取得手段と観念
される。
【００５０】このようにして最高解像度の高周波信号Ｓ
１が得られると、本発明の画像処理手段としての演算器
２５において、下記の式（８）に示すように、高周波信
号Ｓ１に対して原画像信号Ｓorg の値に応じた、周波数
強調の度合いを表すパラメータとしての強調度βが乗じ
られ、さらにこの強調度βが乗じられた高周波信号Ｓ１
が原画像信号Ｓorg と加算されて処理済み画像信号Ｓpr
ocが得られる。
１が得られると、本発明の画像処理手段としての演算器
２５において、下記の式（８）に示すように、高周波信
号Ｓ１に対して原画像信号Ｓorg の値に応じた、周波数
強調の度合いを表すパラメータとしての強調度βが乗じ
られ、さらにこの強調度βが乗じられた高周波信号Ｓ１
が原画像信号Ｓorg と加算されて処理済み画像信号Ｓpr
ocが得られる。
【００５１】 Ｓproc＝Ｓorg ＋β（Ｓorg ）・Ｓ１ （８） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） なお、強調度β（Ｓorg）の値を変更することにより、
処理済み画像信号Ｓprocの画像処理の程度を任意に変更
することができる。
数） なお、強調度β（Ｓorg）の値を変更することにより、
処理済み画像信号Ｓprocの画像処理の程度を任意に変更
することができる。
【００５２】次いで、第１の実施形態の動作について説
明する。図１０は第１の実施形態の動作を示すフローチ
ャートである。まず、読取装置等から原画像信号Ｓorg
が画像処理装置１に入力される（ステップＳ１）。原画
像信号Ｓorgは帯域制限画像信号作成手段２に入力され
てここで原画像信号Ｓorgの周波数帯域毎の周波数応答
特性を表す帯域制限画像信号Ｂkが作成される（ステッ
プＳ２）。帯域制限画像信号Ｂkは、上記図９に示すよ
うな変換関数により変換されて変換済み信号ｆkＢkが得
られ（ステップＳ３）、さらに変換済み信号ｆkＢkの１
段階高周波数帯域への補間処理による高周波信号Ｓｋの
取得、および高周波信号Ｓｋと対応する周波数帯域の変
換済み信号ｆkＢkとの加算による高周波信号Ｓｋ−１の
取得を最高周波数帯域の変換済み信号ｆ1Ｂ1まで繰り返
し行って、高周波信号Ｓ１を得る（ステップＳ４）。そ
して高周波信号Ｓ１を用いて上記式（８）に示す演算を
行って処理済み画像信号Ｓprocを得て（ステップＳ
５）、処理を終了する。
明する。図１０は第１の実施形態の動作を示すフローチ
ャートである。まず、読取装置等から原画像信号Ｓorg
が画像処理装置１に入力される（ステップＳ１）。原画
像信号Ｓorgは帯域制限画像信号作成手段２に入力され
てここで原画像信号Ｓorgの周波数帯域毎の周波数応答
特性を表す帯域制限画像信号Ｂkが作成される（ステッ
プＳ２）。帯域制限画像信号Ｂkは、上記図９に示すよ
うな変換関数により変換されて変換済み信号ｆkＢkが得
られ（ステップＳ３）、さらに変換済み信号ｆkＢkの１
段階高周波数帯域への補間処理による高周波信号Ｓｋの
取得、および高周波信号Ｓｋと対応する周波数帯域の変
換済み信号ｆkＢkとの加算による高周波信号Ｓｋ−１の
取得を最高周波数帯域の変換済み信号ｆ1Ｂ1まで繰り返
し行って、高周波信号Ｓ１を得る（ステップＳ４）。そ
して高周波信号Ｓ１を用いて上記式（８）に示す演算を
行って処理済み画像信号Ｓprocを得て（ステップＳ
５）、処理を終了する。
【００５３】このように、第１の実施形態においては、
強調度β（Ｓorg）の値を変更するのみで、処理済み画
像信号Ｓprocの画像処理の程度を変更することができる
ため、上記特開平5-244508号に記載されたように非線形
関数を修正することにより画像処理の程度を変更する方
法と比較して、画像処理の程度を簡易に変更することが
でき、これにより、処理済み画像信号Ｓprocを得るため
の演算時間を短縮することができる。したがって、例え
ば画像処理の内容を切り替えて連続的に処理済み画像信
号ＳprocをＣＲＴ等に表示する場合にも、処理時間を低
減してオペレータのストレスを低減することができる。
強調度β（Ｓorg）の値を変更するのみで、処理済み画
像信号Ｓprocの画像処理の程度を変更することができる
ため、上記特開平5-244508号に記載されたように非線形
関数を修正することにより画像処理の程度を変更する方
法と比較して、画像処理の程度を簡易に変更することが
でき、これにより、処理済み画像信号Ｓprocを得るため
の演算時間を短縮することができる。したがって、例え
ば画像処理の内容を切り替えて連続的に処理済み画像信
号ＳprocをＣＲＴ等に表示する場合にも、処理時間を低
減してオペレータのストレスを低減することができる。
【００５４】ここで、第１の実施形態による高周波信号
Ｓ１を得る処理の演算量と従来の演算量とを比較する。
図１１は例えば上記特開平10-75395号に記載された従来
の処理を模式的に示す図である。図１１に示すように、
従来の処理は、本実施形態と同様に低解像度画像信号Ｌ
kを求め、この低解像度画像信号Ｌkに対して補間処理を
施すことにより、原画像信号Ｓorgと同一の画素数を有
するボケ画像信号Ｓuskを得、このボケ画像信号Ｓusk間
および原画像信号Ｓorgとボケ画像信号Ｓus１との減算
を行って、原画像と同一画素数の帯域制限画像信号Ｂk
を得、各帯域制限画像信号Ｂkに対して変換関数による
変換処理を行った後、変換後の帯域制限画像信号に対し
て強調係数βを乗じたものを原画像信号Ｓorgに加算す
ることにより、処理済み画像信号Ｓprocを得るものであ
る。この処理を下記の式（９）に示す。
Ｓ１を得る処理の演算量と従来の演算量とを比較する。
図１１は例えば上記特開平10-75395号に記載された従来
の処理を模式的に示す図である。図１１に示すように、
従来の処理は、本実施形態と同様に低解像度画像信号Ｌ
kを求め、この低解像度画像信号Ｌkに対して補間処理を
施すことにより、原画像信号Ｓorgと同一の画素数を有
するボケ画像信号Ｓuskを得、このボケ画像信号Ｓusk間
および原画像信号Ｓorgとボケ画像信号Ｓus１との減算
を行って、原画像と同一画素数の帯域制限画像信号Ｂk
を得、各帯域制限画像信号Ｂkに対して変換関数による
変換処理を行った後、変換後の帯域制限画像信号に対し
て強調係数βを乗じたものを原画像信号Ｓorgに加算す
ることにより、処理済み画像信号Ｓprocを得るものであ
る。この処理を下記の式（９）に示す。
【００５５】 Ｓproc＝Ｓorg ＋β（Ｓorg ）×Ｆusm （Ｓorg,Ｓus1,Ｓus2,…Ｓusn） Ｆusm （Ｓorg,Ｓus1,Ｓus2,…Ｓusn） ＝ｆ1(Ｓorg −Ｓus1)＋ｆ2(Ｓus1 −Ｓus2)＋… ＋ｆk(Ｓusk-1−Ｓusk )＋…＋ｆn(Ｓusn-1−Ｓusn) （９） （但し、Ｓproc：高周波成分が強調された画像信号 Ｓorg ：原画像信号 Ｓusk(k=1〜n）：ボケ画像信号（原画像と同一画素数） ｆk(k=1〜n）：各帯域制限画像信号を変換する変換関数 β（Ｓorg）：原画像信号に基づいて定められる強調係
数） このような従来の処理においては、低解像度画像信号Ｌ
kに対して原画像信号Ｓorgと同一画素数となるような補
間処理を施しているため、例えば原画像の画素数が１０
２４×１０２４であり、この原画像を表す原画像信号Ｓ
orgから６段階の低解像度画像信号Ｌ1〜Ｌ6を得た場
合、補間処理が４×４の範囲における１６画素を用いて
１つの画素を求める演算を行うものであるとすると、補
間処理の演算量としては、(1024×1024)×16×６＝1006
63296回の演算が必要となる。これに対して本実施形態
においては、ボケ画像信号Ｓuskはその周波数帯域に応
じた画素数を有しており、変換済み信号ｆkＢkから得ら
れる高周波信号Ｓｋを拡大するときに補間処理を行うも
のであるため、その演算量としては、((1024×1024)+(5
12×512)+(256×256)+(128×128)+(64×64)+(32×32))
×16=22364160となる。なお、実際にはボケ画像信号Ｓu
skを求める際にも補間処理を行っているため、補間処理
のトータル演算量は44728320回となる。この演算量は上
記従来の演算量の約２．２５倍少なくなっているため、
演算時間を約２．２５倍短縮することができる。
数） このような従来の処理においては、低解像度画像信号Ｌ
kに対して原画像信号Ｓorgと同一画素数となるような補
間処理を施しているため、例えば原画像の画素数が１０
２４×１０２４であり、この原画像を表す原画像信号Ｓ
orgから６段階の低解像度画像信号Ｌ1〜Ｌ6を得た場
合、補間処理が４×４の範囲における１６画素を用いて
１つの画素を求める演算を行うものであるとすると、補
間処理の演算量としては、(1024×1024)×16×６＝1006
63296回の演算が必要となる。これに対して本実施形態
においては、ボケ画像信号Ｓuskはその周波数帯域に応
じた画素数を有しており、変換済み信号ｆkＢkから得ら
れる高周波信号Ｓｋを拡大するときに補間処理を行うも
のであるため、その演算量としては、((1024×1024)+(5
12×512)+(256×256)+(128×128)+(64×64)+(32×32))
×16=22364160となる。なお、実際にはボケ画像信号Ｓu
skを求める際にも補間処理を行っているため、補間処理
のトータル演算量は44728320回となる。この演算量は上
記従来の演算量の約２．２５倍少なくなっているため、
演算時間を約２．２５倍短縮することができる。
【００５６】なお、本発明においては、帯域制限画像信
号の作成方法としては、上記第１の実施形態において説
明した方法に限定されるものではなく、上記特開平10-7
5395号に記載された方法により帯域制限画像信号を作成
してもよいものである。
号の作成方法としては、上記第１の実施形態において説
明した方法に限定されるものではなく、上記特開平10-7
5395号に記載された方法により帯域制限画像信号を作成
してもよいものである。
【００５７】また、上記第１の実施形態においては、ラ
プラシアンピラミッドの手法により原画像信号Ｓorgか
ら各周波数帯域毎の特性を表す帯域制限画像信号を得て
いるが、例えば特開平6-274615号に示すように、ウェー
ブレット変換により帯域制限画像信号を得るようにして
もよい。以下ウェーブレット変換を用いた画像処理の実
施形態を第２の実施形態として説明する。
プラシアンピラミッドの手法により原画像信号Ｓorgか
ら各周波数帯域毎の特性を表す帯域制限画像信号を得て
いるが、例えば特開平6-274615号に示すように、ウェー
ブレット変換により帯域制限画像信号を得るようにして
もよい。以下ウェーブレット変換を用いた画像処理の実
施形態を第２の実施形態として説明する。
【００５８】図１２は本発明の第２の実施形態による画
像処理装置の構成を示す概略ブロック図である。図１２
に示すように、本発明の第２の実施形態による画像処理
装置３１は、読取装置等において得られた所定の解像度
を有する原画像信号Ｓorgをウェーブレット変換するウ
ェーブレット変換手段３２と、ウェーブレット変換によ
り得られた信号に基づいて原画像信号Ｓorgに対して特
定の周波数を強調するための周波数強調処理を行って処
理済み画像信号Ｓprocを得る変換処理手段３３とを有す
る。なお、第２の実施形態における主走査方向および副
走査方向とは、原画像に対して図２２に示す方向とす
る。
像処理装置の構成を示す概略ブロック図である。図１２
に示すように、本発明の第２の実施形態による画像処理
装置３１は、読取装置等において得られた所定の解像度
を有する原画像信号Ｓorgをウェーブレット変換するウ
ェーブレット変換手段３２と、ウェーブレット変換によ
り得られた信号に基づいて原画像信号Ｓorgに対して特
定の周波数を強調するための周波数強調処理を行って処
理済み画像信号Ｓprocを得る変換処理手段３３とを有す
る。なお、第２の実施形態における主走査方向および副
走査方向とは、原画像に対して図２２に示す方向とす
る。
【００５９】図１３はウェーブレット変換手段３２の構
成を示す概略ブロック図である。なお、本実施形態にお
いては、ウェーブレット変換の各係数が直交する直交ウ
ェーブレット変換を行うものである。
成を示す概略ブロック図である。なお、本実施形態にお
いては、ウェーブレット変換の各係数が直交する直交ウ
ェーブレット変換を行うものである。
【００６０】まず、図１３に示すように原画像信号Ｓor
gに対してウェーブレット変換部４１においてウェーブ
レット変換が施される。図１４はウェーブレット変換部
４１において行われる処理を示すブロック図である。図
１４に示すように、原画像信号Ｓorg（信号ＬＬｋ）の
主走査方向に基本ウェーブレット関数Ｈ，Ｇによりフィ
ルタリング処理を行うとともに、主走査方向の画素を１
画素おきに間引き（図中↓２で表す）、主走査方向の画
素数を１／２にする。ここで、関数Ｈはハイパスフィル
タであり、関数Ｇはローパスフィルタである。さらに、
この画素が間引かれた信号のそれぞれに対して副走査方
向に関数Ｈ，Ｇによりフィルタリング処理を行うととも
に、副走査方向の画素を１画素おきに間引き、副走査方
向の画素数を１／２にして、ウェーブレット変換係数信
号（以下単に信号とすることもある）ＨＨ１，ＨＬ１，
ＬＨ１，ＬＬ１（ＨＨｋ＋１，ＨＬｋ＋１，ＬＨｋ＋
１，ＬＬｋ＋１）を得る。ここで、信号ＬＬ１は原画像
の縦横を１／２に縮小した画像を表し、信号ＨＬ１、Ｌ
Ｈ１およびＨＨ１はそれぞれ原画像の１／２縮小画像に
おいて縦エッジ、横エッジおよび斜めエッジ成分の画像
を表すものとなる。
gに対してウェーブレット変換部４１においてウェーブ
レット変換が施される。図１４はウェーブレット変換部
４１において行われる処理を示すブロック図である。図
１４に示すように、原画像信号Ｓorg（信号ＬＬｋ）の
主走査方向に基本ウェーブレット関数Ｈ，Ｇによりフィ
ルタリング処理を行うとともに、主走査方向の画素を１
画素おきに間引き（図中↓２で表す）、主走査方向の画
素数を１／２にする。ここで、関数Ｈはハイパスフィル
タであり、関数Ｇはローパスフィルタである。さらに、
この画素が間引かれた信号のそれぞれに対して副走査方
向に関数Ｈ，Ｇによりフィルタリング処理を行うととも
に、副走査方向の画素を１画素おきに間引き、副走査方
向の画素数を１／２にして、ウェーブレット変換係数信
号（以下単に信号とすることもある）ＨＨ１，ＨＬ１，
ＬＨ１，ＬＬ１（ＨＨｋ＋１，ＨＬｋ＋１，ＬＨｋ＋
１，ＬＬｋ＋１）を得る。ここで、信号ＬＬ１は原画像
の縦横を１／２に縮小した画像を表し、信号ＨＬ１、Ｌ
Ｈ１およびＨＨ１はそれぞれ原画像の１／２縮小画像に
おいて縦エッジ、横エッジおよび斜めエッジ成分の画像
を表すものとなる。
【００６１】次に、信号ＬＬ１に対してさらにウェーブ
レット変換部４１においてウェーブレット変換が施され
て、信号ＨＨ２，ＨＬ２，ＬＨ２，ＬＬ２が得られる。
ここで、信号ＬＬ２は原画像の縦横を１／４に縮小した
画像を表し、信号ＨＬ２、ＬＨ２およびＨＨ２はそれぞ
れ原画像の１／４縮小画像において縦エッジ、横エッジ
および斜めエッジ成分の画像を表すものとなる。
レット変換部４１においてウェーブレット変換が施され
て、信号ＨＨ２，ＨＬ２，ＬＨ２，ＬＬ２が得られる。
ここで、信号ＬＬ２は原画像の縦横を１／４に縮小した
画像を表し、信号ＨＬ２、ＬＨ２およびＨＨ２はそれぞ
れ原画像の１／４縮小画像において縦エッジ、横エッジ
および斜めエッジ成分の画像を表すものとなる。
【００６２】以下、上記と同様にして、各周波数帯域に
おいて得られるウェーブレット変換係数信号ＬＬｋに対
するウェーブレット変換をｎ回繰り返すことによりウェ
ーブレット変換係数信号ＨＨ１〜ＨＨｎ，ＨＬ１〜ＨＬ
ｎ，ＬＨ１〜ＬＨｎ，ＬＬ１〜ＬＬｎを得る。ここで、
ｎ回目のウェーブレット変換により得られるウェーブレ
ット変換係数信号ＨＨｎ，ＨＬｎ，ＬＨｎ，ＬＬｎは、
原画像信号Ｓorgと比較して主副各方向の画素数が(1/
2）n となっているため、各ウェーブレット変換係数信
号はｎが大きいほど周波数帯域が低く、原画像データの
周波数成分のうち低周波数成分を表すデータとなる。し
たがって、ウェーブレット変換係数信号ＨＨｋ（ｋ＝０
〜ｎ、以下同様）は、原画像信号Ｓorgの主副両方向の
周波数の変化を表すものであり、ｋが大きいほど低周波
信号となる。またウェーブレット変換係数信号ＨＬｋは
原画像信号Ｓorgの主走査方向の周波数の変化を表すも
のであり、ｋが大きいほど低周波信号となる。さらにウ
ェーブレット変換係数信号ＬＨｋは原画像信号Ｓorg の
副走査方向の周波数の変化を表すものであり、ｋが大き
いほど低周波信号となる。
おいて得られるウェーブレット変換係数信号ＬＬｋに対
するウェーブレット変換をｎ回繰り返すことによりウェ
ーブレット変換係数信号ＨＨ１〜ＨＨｎ，ＨＬ１〜ＨＬ
ｎ，ＬＨ１〜ＬＨｎ，ＬＬ１〜ＬＬｎを得る。ここで、
ｎ回目のウェーブレット変換により得られるウェーブレ
ット変換係数信号ＨＨｎ，ＨＬｎ，ＬＨｎ，ＬＬｎは、
原画像信号Ｓorgと比較して主副各方向の画素数が(1/
2）n となっているため、各ウェーブレット変換係数信
号はｎが大きいほど周波数帯域が低く、原画像データの
周波数成分のうち低周波数成分を表すデータとなる。し
たがって、ウェーブレット変換係数信号ＨＨｋ（ｋ＝０
〜ｎ、以下同様）は、原画像信号Ｓorgの主副両方向の
周波数の変化を表すものであり、ｋが大きいほど低周波
信号となる。またウェーブレット変換係数信号ＨＬｋは
原画像信号Ｓorgの主走査方向の周波数の変化を表すも
のであり、ｋが大きいほど低周波信号となる。さらにウ
ェーブレット変換係数信号ＬＨｋは原画像信号Ｓorg の
副走査方向の周波数の変化を表すものであり、ｋが大き
いほど低周波信号となる。
【００６３】ここで、図１５にウェーブレット変換係数
信号を複数の周波数帯域毎に示す。なお、図１５におい
ては便宜上２回目のウェーブレット変換を行った状態ま
でを表すものとする。なお、図１５において信号ＬＬ２
は原画像を主副各方向が１／４に縮小した画像を表すも
のとなっている。
信号を複数の周波数帯域毎に示す。なお、図１５におい
ては便宜上２回目のウェーブレット変換を行った状態ま
でを表すものとする。なお、図１５において信号ＬＬ２
は原画像を主副各方向が１／４に縮小した画像を表すも
のとなっている。
【００６４】なお、ウェーブレット変換係数信号ＨＨ
ｋ，ＨＬｋ，ＬＨｋ，ＬＬｋ（ｋ＝１〜ｎ）のうち、信
号ＨＨｋ，ＨＬｋ，ＬＨｋはその周波数帯域におけるエ
ッジ成分を表すものであり、換言すれば原画像における
特定の周波数帯域（帯域制限画像特性）を有する画像を
表すもの、すなわち主にその周波数帯域における画像の
コントラストを表すものとなっている。また、ウェーブ
レット変換係数信号ＬＬｋは上述したように原画像を縮
小した画像を表すものとなっている。なお、本実施形態
においては、ウェーブレット変換係数信号ＨＨｋ，ＨＬ
ｋ，ＬＨｋを帯域制限画像信号と称し、ウェーブレット
変換係数信号ＬＬｋを解像度信号と称し、帯域制限画像
信号および解像度信号を総称してウェーブレット変換係
数信号と称するものとする。ここで、最低解像度の信号
ＬＬｎは帯域制限画像信号を得る目的のためには必要な
いため、値を０とする。
ｋ，ＨＬｋ，ＬＨｋ，ＬＬｋ（ｋ＝１〜ｎ）のうち、信
号ＨＨｋ，ＨＬｋ，ＬＨｋはその周波数帯域におけるエ
ッジ成分を表すものであり、換言すれば原画像における
特定の周波数帯域（帯域制限画像特性）を有する画像を
表すもの、すなわち主にその周波数帯域における画像の
コントラストを表すものとなっている。また、ウェーブ
レット変換係数信号ＬＬｋは上述したように原画像を縮
小した画像を表すものとなっている。なお、本実施形態
においては、ウェーブレット変換係数信号ＨＨｋ，ＨＬ
ｋ，ＬＨｋを帯域制限画像信号と称し、ウェーブレット
変換係数信号ＬＬｋを解像度信号と称し、帯域制限画像
信号および解像度信号を総称してウェーブレット変換係
数信号と称するものとする。ここで、最低解像度の信号
ＬＬｎは帯域制限画像信号を得る目的のためには必要な
いため、値を０とする。
【００６５】変換処理手段３３は、上記第１の実施形態
における変換処理手段３と同様に周波数強調処理を行う
ものである。図１６は変換処理手段３３の構成をウェー
ブレット変換手段３２とともに示す概略ブロック図であ
る。図１６に示すように、ウェーブレット変換手段３２
において得られた帯域制限画像信号ＨＨｋ，ＨＬｋ，Ｌ
Ｈｋ（以下Ｂkで代表させる）が、変換器５１において
周波数帯域毎に異なる変換関数ｆ1〜ｆnにより所望の大
きさとなるように抑制されて、変換済み信号Ｂ k′（Ｈ
Ｈｋ′，ＨＬｋ′，ＬＨｋ′、ｋ＝１〜ｎ）が得られ
る。そして、変換済み信号ＨＨｋ′，ＨＬｋ′，ＬＨ
ｋ′に対して逆ウェーブレット変換手段５２において逆
ウェーブレット変換が施される。図１７は、逆ウェーブ
レット変換手段５２において行われる逆ウェーブレット
変換を説明するための図である。図１７に示すように、
最低周波数帯域の変換済み信号ＨＨｎ′，ＨＬｎ′，Ｌ
Ｈｎ′，ＬＬｎ（＝０）に対して逆ウェーブレット変換
手段５２において逆ウェーブレット変換を施して処理済
み信号ＬＬｎ−１′を得る。
における変換処理手段３と同様に周波数強調処理を行う
ものである。図１６は変換処理手段３３の構成をウェー
ブレット変換手段３２とともに示す概略ブロック図であ
る。図１６に示すように、ウェーブレット変換手段３２
において得られた帯域制限画像信号ＨＨｋ，ＨＬｋ，Ｌ
Ｈｋ（以下Ｂkで代表させる）が、変換器５１において
周波数帯域毎に異なる変換関数ｆ1〜ｆnにより所望の大
きさとなるように抑制されて、変換済み信号Ｂ k′（Ｈ
Ｈｋ′，ＨＬｋ′，ＬＨｋ′、ｋ＝１〜ｎ）が得られ
る。そして、変換済み信号ＨＨｋ′，ＨＬｋ′，ＬＨ
ｋ′に対して逆ウェーブレット変換手段５２において逆
ウェーブレット変換が施される。図１７は、逆ウェーブ
レット変換手段５２において行われる逆ウェーブレット
変換を説明するための図である。図１７に示すように、
最低周波数帯域の変換済み信号ＨＨｎ′，ＨＬｎ′，Ｌ
Ｈｎ′，ＬＬｎ（＝０）に対して逆ウェーブレット変換
手段５２において逆ウェーブレット変換を施して処理済
み信号ＬＬｎ−１′を得る。
【００６６】図１８は逆ウェーブレット変換手段５２に
おいて行われる処理を示すブロック図である。図１８に
示すように変換済み信号ＬＬｎ′（ＬＬｋ′，ｋ＝ｎの
場合ＬＬｎ＝０）および変換済み信号ＬＨｎ′（ＬＨ
ｋ′）の副走査方向に対して画素間に１画素分の間隔を
あける処理を行うとともに（図中↑２で表す）、関数
Ｇ，Ｈに対応する逆ウェーブレット変換関数Ｇ′，Ｈ′
によりフィルタリング処理を副走査方向に施してこれら
を加算し、さらに加算により得られた信号（第１の加算
信号とする）の主走査方向に対して画素間に１画素分の
間隔をあける処理を行うとともに、関数Ｇ′によりフィ
ルタリング処理を主走査方向に施して第１の信号を得
る。一方、信号ＨＬｎ′（ＨＬｋ′）および信号ＨＨ
ｎ′（ＨＨｋ′）の副走査方向に対して画素間に１画素
分の間隔をあける処理を行うとともに、関数Ｇ′，Ｈ′
によりフィルタリング処理を副走査方向に施してこれら
を加算し、さらに加算により得られた信号（第２の加算
信号とする）の主走査方向に対して画素間に１画素分の
間隔をあける処理を行うとともに、関数Ｈ′によりフィ
ルタリング処理を主走査方向に施して第２の信号を得
る。そして第１および第２の信号を加算して変換済み信
号ＬＬｎ−１′（ＬＬｋ−１′）を得る。なお、最低解
像度のウェーブレット変換係数信号ＬＬｎは０とされて
いるため、変換済み信号ＬＬｎ−１′は原画像信号Ｓor
gのバンドパス特性を表すものとなる。
おいて行われる処理を示すブロック図である。図１８に
示すように変換済み信号ＬＬｎ′（ＬＬｋ′，ｋ＝ｎの
場合ＬＬｎ＝０）および変換済み信号ＬＨｎ′（ＬＨ
ｋ′）の副走査方向に対して画素間に１画素分の間隔を
あける処理を行うとともに（図中↑２で表す）、関数
Ｇ，Ｈに対応する逆ウェーブレット変換関数Ｇ′，Ｈ′
によりフィルタリング処理を副走査方向に施してこれら
を加算し、さらに加算により得られた信号（第１の加算
信号とする）の主走査方向に対して画素間に１画素分の
間隔をあける処理を行うとともに、関数Ｇ′によりフィ
ルタリング処理を主走査方向に施して第１の信号を得
る。一方、信号ＨＬｎ′（ＨＬｋ′）および信号ＨＨ
ｎ′（ＨＨｋ′）の副走査方向に対して画素間に１画素
分の間隔をあける処理を行うとともに、関数Ｇ′，Ｈ′
によりフィルタリング処理を副走査方向に施してこれら
を加算し、さらに加算により得られた信号（第２の加算
信号とする）の主走査方向に対して画素間に１画素分の
間隔をあける処理を行うとともに、関数Ｈ′によりフィ
ルタリング処理を主走査方向に施して第２の信号を得
る。そして第１および第２の信号を加算して変換済み信
号ＬＬｎ−１′（ＬＬｋ−１′）を得る。なお、最低解
像度のウェーブレット変換係数信号ＬＬｎは０とされて
いるため、変換済み信号ＬＬｎ−１′は原画像信号Ｓor
gのバンドパス特性を表すものとなる。
【００６７】次に、変換済み信号ＨＨｎ−１′，ＨＬｎ
−１′，ＬＨｎ−１′，ＬＬｎ−１′に対して上記と同
様に逆ウェーブレット変換手段５２において逆ウェーブ
レット変換を行って、変換済み信号ＬＬｎ−２′を得
る。そして、以下上記と同様にして逆ウェーブレット変
換を最高周波数帯域まで繰り返し、さらに変換済み信号
ＨＨ１′，ＨＬ１′，ＬＬ１′を逆ウェーブレット変換
することにより高周波信号Ｓ１が得られる。
−１′，ＬＨｎ−１′，ＬＬｎ−１′に対して上記と同
様に逆ウェーブレット変換手段５２において逆ウェーブ
レット変換を行って、変換済み信号ＬＬｎ−２′を得
る。そして、以下上記と同様にして逆ウェーブレット変
換を最高周波数帯域まで繰り返し、さらに変換済み信号
ＨＨ１′，ＨＬ１′，ＬＬ１′を逆ウェーブレット変換
することにより高周波信号Ｓ１が得られる。
【００６８】得られた高周波信号Ｓ１は、演算器５３に
おいて上記第１の実施形態と同様に式（８）に示すよう
な演算が行われて処理済み画像信号Ｓprocが得られる。
おいて上記第１の実施形態と同様に式（８）に示すよう
な演算が行われて処理済み画像信号Ｓprocが得られる。
【００６９】次いで、第２の実施形態の動作について説
明する。図１９は第２の実施形態の動作を示すフローチ
ャートである。まず、読取装置等から原画像信号Ｓorg
が画像処理装置１に入力される（ステップＳ１１）。原
画像信号Ｓorgはウェーブレット変換手段３２において
ウェーブレット変換が施されて各周波数帯域毎のウェー
ブレット変換係数信号が得られる（ステップＳ１２）。
次に、各ウェーブレット変換係数信号Ｂkが上述した変
換関数により変換されて変換済み信号Ｂk′が得られる
（ステップＳ１３）。この後、変換済み信号Ｂｋ′が、
逆ウェーブレット変換手段５２において逆ウェーブレッ
ト変換されて高周波信号Ｓ１が得られる（ステップＳ１
４）。そして高周波信号Ｓ１を用いて上記式（８）に示
す演算を行って処理済み画像信号Ｓprocを得て（ステッ
プＳ１５）、処理を終了する。
明する。図１９は第２の実施形態の動作を示すフローチ
ャートである。まず、読取装置等から原画像信号Ｓorg
が画像処理装置１に入力される（ステップＳ１１）。原
画像信号Ｓorgはウェーブレット変換手段３２において
ウェーブレット変換が施されて各周波数帯域毎のウェー
ブレット変換係数信号が得られる（ステップＳ１２）。
次に、各ウェーブレット変換係数信号Ｂkが上述した変
換関数により変換されて変換済み信号Ｂk′が得られる
（ステップＳ１３）。この後、変換済み信号Ｂｋ′が、
逆ウェーブレット変換手段５２において逆ウェーブレッ
ト変換されて高周波信号Ｓ１が得られる（ステップＳ１
４）。そして高周波信号Ｓ１を用いて上記式（８）に示
す演算を行って処理済み画像信号Ｓprocを得て（ステッ
プＳ１５）、処理を終了する。
【００７０】このように、第２の実施形態においても、
従来の周波数強調処理方法のように、原画像信号Ｓorg
と同一の画素数を有するボケ画像信号を用いて原画像信
号Ｓorgと同一の画素数を有する帯域制限画像信号Ｂkを
作成するものと比較して、周波数帯域毎のボケ画像信号
を原画像信号と同一の画素数にする補間演算を行う必要
がなくなる。したがって、演算量を低減することがで
き、これにより画像処理を高速に行って、オペレータの
ストレスを低減することができる。
従来の周波数強調処理方法のように、原画像信号Ｓorg
と同一の画素数を有するボケ画像信号を用いて原画像信
号Ｓorgと同一の画素数を有する帯域制限画像信号Ｂkを
作成するものと比較して、周波数帯域毎のボケ画像信号
を原画像信号と同一の画素数にする補間演算を行う必要
がなくなる。したがって、演算量を低減することがで
き、これにより画像処理を高速に行って、オペレータの
ストレスを低減することができる。
【００７１】また、強調度β（Ｓorg）の値を変更する
のみで、処理済み画像信号Ｓprocの画像処理の程度を変
更することができるため、上記特開平5-244508号に記載
されたように、非線形関数を修正し、画像を再構成する
ことにより画像処理の程度を変更する方法と比較して、
画像処理の程度を簡易に変更することができ、これによ
り、処理済み画像信号Ｓprocを得るための演算時間を短
縮して、オペレータのストレスを低減することができ
る。
のみで、処理済み画像信号Ｓprocの画像処理の程度を変
更することができるため、上記特開平5-244508号に記載
されたように、非線形関数を修正し、画像を再構成する
ことにより画像処理の程度を変更する方法と比較して、
画像処理の程度を簡易に変更することができ、これによ
り、処理済み画像信号Ｓprocを得るための演算時間を短
縮して、オペレータのストレスを低減することができ
る。
【００７２】なお、上記第１および第２の実施形態にお
いては、逆多重解像度変換を最高周波数帯域まで行っ
て、原画像信号Ｓorgにより表されると同一解像度の画
像を表す高周波信号Ｓ１を得ているが、逆多重解像度変
換を最高周波数帯域まで行うことなく、途中の周波数帯
域まで行うようにしてもよい。この場合、途中の周波数
帯域まで逆多重解像度変換を行うことにより得られた高
周波信号に基づいて、原画像信号Ｓorgに対して画像処
理を行うことにより、原画像信号Ｓorgに含まれるその
周波数帯域までの成分を強調することができることとな
る。
いては、逆多重解像度変換を最高周波数帯域まで行っ
て、原画像信号Ｓorgにより表されると同一解像度の画
像を表す高周波信号Ｓ１を得ているが、逆多重解像度変
換を最高周波数帯域まで行うことなく、途中の周波数帯
域まで行うようにしてもよい。この場合、途中の周波数
帯域まで逆多重解像度変換を行うことにより得られた高
周波信号に基づいて、原画像信号Ｓorgに対して画像処
理を行うことにより、原画像信号Ｓorgに含まれるその
周波数帯域までの成分を強調することができることとな
る。
【００７３】なお、上記第１および第２の実施形態にお
いては、演算器２５，５３において行われる画像処理と
して周波数強調処理を行っているが、ダイナミックレン
ジ圧縮処理を行ってもよい。この場合、上述したように
得られた高周波信号Ｓ１に対して、下記の式（１０）に
示すように演算が施されて、原画像信号Ｓorgのダイナ
ミックレンジが圧縮される。また、式（１０）において
は、Ｄ（Ｓorg−Ｓ１）の値を変更することにより、処
理済み画像信号Ｓprocのダイナミックレンジ圧縮の程度
を変更することができる。
いては、演算器２５，５３において行われる画像処理と
して周波数強調処理を行っているが、ダイナミックレン
ジ圧縮処理を行ってもよい。この場合、上述したように
得られた高周波信号Ｓ１に対して、下記の式（１０）に
示すように演算が施されて、原画像信号Ｓorgのダイナ
ミックレンジが圧縮される。また、式（１０）において
は、Ｄ（Ｓorg−Ｓ１）の値を変更することにより、処
理済み画像信号Ｓprocのダイナミックレンジ圧縮の程度
を変更することができる。
【００７４】 Ｓproc＝Ｓorg ＋Ｄ（Ｓorg−Ｓ１） （１０） （但し、Ｓproc：処理済み画像信号 Ｓorg：原画像信号 Ｄ（Ｓorg−Ｓ１）：低周波成分信号に基づいて定めら
れるダイナミックレンジ圧縮係数（ＤはＳorg−Ｓ１を
変換する関数）） なお、ダイナミックレンジ圧縮処理を行う場合には、帯
域制限画像信号に変換処理を行う関数ｆkとしては、例
えば図２０、図２１に示される変換関数、あるいはこれ
らを組み合わせた関数を使用することが好ましい。
れるダイナミックレンジ圧縮係数（ＤはＳorg−Ｓ１を
変換する関数）） なお、ダイナミックレンジ圧縮処理を行う場合には、帯
域制限画像信号に変換処理を行う関数ｆkとしては、例
えば図２０、図２１に示される変換関数、あるいはこれ
らを組み合わせた関数を使用することが好ましい。
【００７５】ここで、図２０に示す変換関数は、振幅の
大きな帯域制限画像信号を抑制するような変換を行うも
のであり、周波数帯域の高い帯域制限画像信号の抑制の
度合いを、周波数帯域の低い帯域制限画像信号よりも強
くするものであるが、これは実際の放射線画像のエッジ
に含まれている高周波成分が、低周波成分に比べてその
振幅が小さいということを考慮したものである。実際の
放射線画像においては、かなり急峻なエッジでさえも正
確な階段状にはなっておらず、高周波成分になるほどそ
の振幅が小さくなっていることが多い。このため、各周
波数成分の振幅に合わせて、周波数の高い帯域制限画像
信号ほど小さい振幅から抑制を行うことが望ましく、本
関数によりそれを実現することができる。
大きな帯域制限画像信号を抑制するような変換を行うも
のであり、周波数帯域の高い帯域制限画像信号の抑制の
度合いを、周波数帯域の低い帯域制限画像信号よりも強
くするものであるが、これは実際の放射線画像のエッジ
に含まれている高周波成分が、低周波成分に比べてその
振幅が小さいということを考慮したものである。実際の
放射線画像においては、かなり急峻なエッジでさえも正
確な階段状にはなっておらず、高周波成分になるほどそ
の振幅が小さくなっていることが多い。このため、各周
波数成分の振幅に合わせて、周波数の高い帯域制限画像
信号ほど小さい振幅から抑制を行うことが望ましく、本
関数によりそれを実現することができる。
【００７６】また図２１の関数は、帯域制限画像信号
を、帯域制限画像信号の絶対値に基づいて決まる、その
絶対値以下の値となるように変換を行うもので、この関
数が低周波帯域を処理する関数であるほど、帯域制限画
像信号の絶対値が０近傍の所定の範囲内の値である帯域
制限画像信号を変換した際に得られる変換画像信号の絶
対値が小さい値であることを特徴とするものである。言
い換えれば、これらの関数はそれぞれ、原点を通り、関
数の傾きがその関数により処理される値に拘わらず１以
下であり、その関数の０近傍における傾きが、低周波帯
域を処理する関数であるほど小さいことを特徴とするも
のである。これらの関数は、変換画像信号を積算して得
られる信号を、原画像信号Ｓorgに加えた場合、原画像
信号Ｓorgと加算された信号とのつなぎ目、すなわち信
号の立ち上がりをより自然なものとするという効果があ
る。
を、帯域制限画像信号の絶対値に基づいて決まる、その
絶対値以下の値となるように変換を行うもので、この関
数が低周波帯域を処理する関数であるほど、帯域制限画
像信号の絶対値が０近傍の所定の範囲内の値である帯域
制限画像信号を変換した際に得られる変換画像信号の絶
対値が小さい値であることを特徴とするものである。言
い換えれば、これらの関数はそれぞれ、原点を通り、関
数の傾きがその関数により処理される値に拘わらず１以
下であり、その関数の０近傍における傾きが、低周波帯
域を処理する関数であるほど小さいことを特徴とするも
のである。これらの関数は、変換画像信号を積算して得
られる信号を、原画像信号Ｓorgに加えた場合、原画像
信号Ｓorgと加算された信号とのつなぎ目、すなわち信
号の立ち上がりをより自然なものとするという効果があ
る。
【００７７】さらに、演算器２５，５３において、下記
の式（１１）に示すように、周波数強調処理とダイナミ
ックレンジ圧縮処理とを同時に行うようにしてもよい。
なお、式（１１）においては、β（Ｓorg）およびＤ
（Ｓorg−Ｓ１）の値を変更することにより、処理済み
画像信号Ｓprocの周波数強調処理およびダイナミックレ
ンジ圧縮の程度を変更することができる。
の式（１１）に示すように、周波数強調処理とダイナミ
ックレンジ圧縮処理とを同時に行うようにしてもよい。
なお、式（１１）においては、β（Ｓorg）およびＤ
（Ｓorg−Ｓ１）の値を変更することにより、処理済み
画像信号Ｓprocの周波数強調処理およびダイナミックレ
ンジ圧縮の程度を変更することができる。
【００７８】 Ｓproc＝Ｓorg ＋β（Ｓorg）・Ｓ１＋Ｄ（Ｓorg−Ｓ１′） （１１） （但し、Ｓproc：処理済み画像信号 Ｓorg：原画像信号 β（Ｓorg）：原画像信号に基づいて定められる強調係
数 Ｄ（Ｓorg−Ｓ１）：低周波成分信号に基づいて定めら
れるダイナミックレンジ圧縮係数 この場合、周波数強調処理を行うための高周波信号Ｓ１
は、帯域制限画像信号に対して上記図９に示す変換関数
により変換処理を行うことにより得、ダイナミックレン
ジ圧縮処理を行うための高周波信号Ｓ１′は、帯域制限
画像信号に対して上記図２０，２１に示す変換関数によ
り変換処理を行うことにより得ればよい。
数 Ｄ（Ｓorg−Ｓ１）：低周波成分信号に基づいて定めら
れるダイナミックレンジ圧縮係数 この場合、周波数強調処理を行うための高周波信号Ｓ１
は、帯域制限画像信号に対して上記図９に示す変換関数
により変換処理を行うことにより得、ダイナミックレン
ジ圧縮処理を行うための高周波信号Ｓ１′は、帯域制限
画像信号に対して上記図２０，２１に示す変換関数によ
り変換処理を行うことにより得ればよい。
【００７９】このように、周波数強調処理とダイナミッ
クレンジ圧縮処理とを行う場合には、補間処理を行うた
めの演算回数は22364160×３＝67092480回となる。この
演算量は従来の演算量の約１．５倍少なくなっているた
め、演算時間を約１．５倍高速化することができる。
クレンジ圧縮処理とを行う場合には、補間処理を行うた
めの演算回数は22364160×３＝67092480回となる。この
演算量は従来の演算量の約１．５倍少なくなっているた
め、演算時間を約１．５倍高速化することができる。
【００８０】さらに、上記各実施形態においては、変換
関数を非線形関数として帯域制限画像信号に対して非線
形処理を施しているが、非線形処理に限定されるもので
はなく、変換関数を線形関数や定数としてもよい。
関数を非線形関数として帯域制限画像信号に対して非線
形処理を施しているが、非線形処理に限定されるもので
はなく、変換関数を線形関数や定数としてもよい。
【図１】本発明の第１の実施形態による画像処理装置の
構成を示す概略ブロック図
構成を示す概略ブロック図
【図２】帯域制限画像信号作成処理の概要を示すブロッ
ク図
ク図
【図３】帯域制限画像信号作成処理を模式的に示す図
【図４】フィルタリング処理に使用されるフィルタの一
例を示す図
例を示す図
【図５】低解像度画像信号作成処理の詳細を示す図
【図６】補間処理に使用されるフィルタの一例を示す図
【図７】変換処理手段の構成を帯域制限画像信号作成手
段とともに示す概略ブロック図
段とともに示す概略ブロック図
【図８】変換処理を模式的に示す図
【図９】非線形関数の一例を示す図
【図１０】第１の実施形態の動作を示すフローチャート
【図１１】従来の変換処理を模式的に示す図
【図１２】本発明の第２の実施形態による画像処理装置
の構成を示す概略ブロック図
の構成を示す概略ブロック図
【図１３】ウェーブレット変換手段の構成を示す概略ブ
ロック図
ロック図
【図１４】ウェーブレット変換部において行われる処理
を示すブロック図
を示すブロック図
【図１５】ウェーブレット変換係数信号を複数の周波数
帯域毎に示す図
帯域毎に示す図
【図１６】変換処理手段の構成を帯域制限画像信号作成
手段とともに示す概略ブロック図
手段とともに示す概略ブロック図
【図１７】逆ウェーブレット変換を説明するための図
【図１８】逆ウェーブレット変換手段において行われる
処理を示すブロック図
処理を示すブロック図
【図１９】第２の実施形態の動作を示すフローチャート
【図２０】ダイナミックレンジ圧縮処理を行う際の非線
形関数の例を示す図（その１）
形関数の例を示す図（その１）
【図２１】ダイナミックレンジ圧縮処理を行う際の非線
形関数の例を示す図（その２）
形関数の例を示す図（その２）
【図２２】原画像の主走査方向および副走査方向を示す
図
図
１，３１ 画像処理装置 ２ 帯域制限画像信号作成手段 ３，３３ 変換処理手段 １０ フィルタリング処理手段 １１，２３ 補間処理手段 １２ 減算器 ２２，５１ 変換器 ２４ 加算器 ２５，５３ 演算器 ３２ ウェーブレット変換手段 ４１ ウェーブレット変換部 ５２ 逆ウェーブレット変換手段
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.7 識別記号 ＦＩ テーマコート゛(参考） Ｈ０４Ｎ 1/41 Ｈ０４Ｎ 5/202 ５Ｊ０６４ 5/202 5/205 5/205 1/40 １０１Ｄ Ｆターム(参考） 5B057 AA07 BA03 CA08 CA12 CA16 CB08 CB12 CB16 CD05 CD11 CE03 CE06 CG02 5C021 PA38 PA66 PA72 PA99 SA08 XA31 XB01 XB07 XB16 5C076 AA21 AA22 AA36 BA06 BA09 BB24 BB40 5C077 LL18 NN03 PP03 PP04 PP10 PP20 PP49 PQ12 RR21 5C078 BA01 CA31 DB13 5J064 AA02 BA16 BB04 BC18 BC29 BD03
Claims (27)
- 【請求項１】 原画像を表す原画像信号に対して、該
原画像の高周波成分に関する信号に基づく画像処理を施
して処理済み画像信号を得る画像処理方法において、 前記原画像信号から帯域制限画像信号を作成し、 所定の変換関数に基づいて、前記帯域制限画像信号から
前記原画像の高周波成分に関する信号を得、 該高周波成分に関する信号に基づいて、前記原画像信号
に対して前記画像処理を施すことを特徴とする画像処理
方法。 - 【請求項２】 前記帯域制限画像信号により表される
各周波数帯域の画像が、その周波数帯域に応じた画素数
を有することを特徴とする請求項１記載の画像処理方
法。 - 【請求項３】 前記高周波成分に関する信号が前記原
画像信号と同一画素数であることを特徴とする請求項１
または２記載の画像処理方法。 - 【請求項４】 前記原画像信号を多重解像度変換する
ことにより、前記帯域制限画像信号を作成し、 該帯域制限画像信号に対して前記所定の変換関数に基づ
いて変換処理を施して変換帯域制限画像信号を得、 該変換帯域制限画像信号を逆多重解像度変換することに
より、前記高周波成分に関する信号を得ることを特徴と
する請求項１から３のいずれか１項記載の画像処理方
法。 - 【請求項５】 前記多重解像度変換は、ラプラシアン
ピラミッド分解による変換、またはウェーブレット変換
であることを特徴とする請求項４記載の画像処理方法。 - 【請求項６】 最低周波数帯域の帯域制限画像信号以
外の帯域制限画像信号から前記高周波成分に関する信号
を得ることを特徴とする請求項１から５のいずれか１項
記載の画像処理方法。 - 【請求項７】 前記所定の変換関数は、非線形関数で
あることを特徴とする請求項１から６のいずれか１項記
載の画像処理方法。 - 【請求項８】 前記画像処理は、周波数強調処理であ
ることを特徴とする請求項１から７のいずれか１項記載
の画像処理方法。 - 【請求項９】 前記画像処理は、ダイナミックレンジ
圧縮処理であることを特徴とする請求項１から８のいず
れか１項記載の画像処理方法。 - 【請求項１０】 原画像を表す原画像信号に対して、
該原画像の高周波成分に関する信号に基づく画像処理を
施して処理済み画像信号を得る画像処理装置において、 前記原画像信号から帯域制限画像信号を作成する帯域制
限画像信号作成手段と、 所定の変換関数に基づいて、前記帯域制限画像信号から
前記原画像の高周波成分に関する信号を得る高周波成分
取得手段と、 該高周波成分に関する信号に基づいて、前記原画像信号
に対して前記画像処理を施す画像処理手段とを備えたこ
とを特徴とする画像処理装置。 - 【請求項１１】 前記帯域制限画像信号作成手段は、
前記帯域制限画像信号により表される各周波数帯域の画
像が、その周波数帯域に応じた画素数を有するよう前記
帯域制限画像信号を作成する手段であることを特徴とす
る請求項１０記載の画像処理装置。 - 【請求項１２】 前記高周波成分取得手段は、前記高
周波成分に関する信号が前記原画像信号と同一画素数と
なるように前記高周波成分に関する信号を取得する手段
であることを特徴とする請求項１０または１１記載の画
像処理装置。 - 【請求項１３】 前記帯域制限画像信号作成手段は、
前記原画像信号を多重解像度変換することにより、前記
帯域制限画像信号を作成する手段であり、 前記高周波成分取得手段は、前記帯域制限画像信号に対
して前記所定の変換関数に基づいて変換処理を施して変
換帯域制限画像信号を得、該変換帯域制限画像信号を逆
多重解像度変換することにより、前記高周波成分に関す
る信号を得る手段であることを特徴とする請求項１０か
ら１２のいずれか１項記載の画像処理装置。 - 【請求項１４】 前記多重解像度変換は、ラプラシア
ンピラミッド分解による変換、またはウェーブレット変
換であることを特徴とする請求項１３記載の画像処理装
置。 - 【請求項１５】 前記高周波成分取得手段は、最低周
波数帯域の帯域制限画像信号以外の帯域制限画像信号か
ら前記高周波成分に関する信号を得る手段であることを
特徴とする請求項１０から１４のいずれか１項記載の画
像処理装置。 - 【請求項１６】 前記所定の変換関数は、非線形関数
であることを特徴とする請求項１０から１５のいずれか
１項記載の画像処理装置。 - 【請求項１７】 前記画像処理は、周波数強調処理で
あることを特徴とする請求項１０から１６のいずれか１
項記載の画像処理装置。 - 【請求項１８】 前記画像処理は、ダイナミックレン
ジ圧縮処理であることを特徴とする請求項１０から１７
のいずれか１項記載の画像処理装置。 - 【請求項１９】 原画像を表す原画像信号に対して、
該原画像の高周波成分に関する信号に基づく画像処理を
施して処理済み画像信号を得る画像処理方法をコンピュ
ータに実行させるためのプログラムを記録したコンピュ
ータ読取り可能な記録媒体において、 前記プログラムは、前記原画像信号から帯域制限画像信
号を作成する手順と、 所定の変換関数に基づいて、前記帯域制限画像信号から
前記原画像の高周波成分に関する信号を得る手順と、 該高周波成分に関する信号に基づいて、前記原画像信号
に対して前記画像処理を施す手順とを有することを特徴
とするコンピュータ読取り可能な記録媒体。 - 【請求項２０】 前記帯域制限画像信号を作成する手
順は、前記帯域制限画像信号により表される各周波数帯
域の画像が、その周波数帯域に応じた画素数を有するよ
う前記帯域制限画像信号を作成する手順であることを特
徴とする請求項１９記載のコンピュータ読取り可能な記
録媒体。 - 【請求項２１】 前記高周波成分に関する信号を得る
手順は、前記高周波成分に関する信号が前記原画像信号
と同一画素数となるように前記高周波成分に関する信号
を得る手順であることを特徴とする請求項１９または２
０記載のコンピュータ読取り可能な記録媒体。 - 【請求項２２】 前記帯域制限画像信号を作成する手
順は、前記原画像信号を多重解像度変換することによ
り、前記帯域制限画像信号を作成する手順であり、 前記高周波成分に関する信号を得る手順は、該帯域制限
画像信号に対して前記所定の変換関数に基づいて変換処
理を施して変換帯域制限画像信号を得る手順と、該変換
帯域制限画像信号を逆多重解像度変換することにより、
前記高周波成分に関する信号を得る手順とを有すること
を特徴とする請求項１９から２１のいずれか１項記載の
コンピュータ読取り可能な記録媒体。 - 【請求項２３】 前記多重解像度変換は、ラプラシア
ンピラミッド分解による変換、またはウェーブレット変
換であることを特徴とする請求項２２記載のコンピュー
タ読取り可能な記録媒体。 - 【請求項２４】 前記高周波成分に関する信号を得る
手順は、最低周波数帯域の帯域制限画像信号以外の帯域
制限画像信号から前記高周波成分に関する信号を得る手
順であることを特徴とする請求項１９から２３のいずれ
か１項記載のコンピュータ読取り可能な記録媒体。 - 【請求項２５】 前記所定の変換関数は、非線形関数
であることを特徴とする請求項１９から２４のいずれか
１項記載のコンピュータ読取り可能な記録媒体。 - 【請求項２６】 前記画像処理は、周波数強調処理で
あることを特徴とする請求項１９から２５のいずれか１
項記載のコンピュータ読取り可能な記録媒体。 - 【請求項２７】 前記画像処理は、ダイナミックレン
ジ圧縮処理であることを特徴とする請求項１９から２６
のいずれか１項記載のコンピュータ読取り可能な記録媒
体。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000143166A JP2001056856A (ja) | 1999-06-10 | 2000-05-16 | 画像処理方法および装置並びに記録媒体 |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP11-163670 | 1999-06-10 | ||
| JP16367099 | 1999-06-10 | ||
| JP2000143166A JP2001056856A (ja) | 1999-06-10 | 2000-05-16 | 画像処理方法および装置並びに記録媒体 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2001056856A true JP2001056856A (ja) | 2001-02-27 |
Family
ID=26489052
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2000143166A Pending JP2001056856A (ja) | 1999-06-10 | 2000-05-16 | 画像処理方法および装置並びに記録媒体 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2001056856A (ja) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006506021A (ja) * | 2002-11-06 | 2006-02-16 | ディジビジョン・インコーポレイテッド | 多次元画像の画質向上のためのシステムおよび方法 |
| JP2009118337A (ja) * | 2007-11-08 | 2009-05-28 | Olympus Corp | 映像信号処理装置、映像信号処理プログラム |
| JP2014524084A (ja) * | 2011-06-27 | 2014-09-18 | イスム リサーチ ディべロップメント カンパニー オブ ザ ヒーブリュー ユニバーシティ オブ エルサレム エルティーディー | 画像処理目的のためのフィルタによるコンボルーションの高速数値近似の適用 |
| JP2021184249A (ja) * | 2020-05-22 | 2021-12-02 | ペキン シャオミ モバイル ソフトウェア カンパニー, リミテッドBeijing Xiaomi Mobile Software Co., Ltd. | 画像処理方法、端末デバイスにおけるイメージシグナルプロセッサ |
| WO2024057768A1 (ja) * | 2022-09-14 | 2024-03-21 | 富士フイルム株式会社 | 画像生成装置、学習装置、画像処理装置、画像生成方法、学習方法及び画像処理方法 |
-
2000
- 2000-05-16 JP JP2000143166A patent/JP2001056856A/ja active Pending
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006506021A (ja) * | 2002-11-06 | 2006-02-16 | ディジビジョン・インコーポレイテッド | 多次元画像の画質向上のためのシステムおよび方法 |
| US7668390B2 (en) | 2002-11-06 | 2010-02-23 | Digivision, Inc. | Systems and methods for image enhancement in multiple dimensions |
| JP2010213284A (ja) * | 2002-11-06 | 2010-09-24 | Digivision Inc | 多次元画像の画質向上のためのシステムおよび方法 |
| JP2009118337A (ja) * | 2007-11-08 | 2009-05-28 | Olympus Corp | 映像信号処理装置、映像信号処理プログラム |
| JP2014524084A (ja) * | 2011-06-27 | 2014-09-18 | イスム リサーチ ディべロップメント カンパニー オブ ザ ヒーブリュー ユニバーシティ オブ エルサレム エルティーディー | 画像処理目的のためのフィルタによるコンボルーションの高速数値近似の適用 |
| JP2021184249A (ja) * | 2020-05-22 | 2021-12-02 | ペキン シャオミ モバイル ソフトウェア カンパニー, リミテッドBeijing Xiaomi Mobile Software Co., Ltd. | 画像処理方法、端末デバイスにおけるイメージシグナルプロセッサ |
| US11625817B2 (en) | 2020-05-22 | 2023-04-11 | Beijing Xiaomi Mobile Software Co., Ltd. | Pyramid-based tone mapping |
| JP7297010B2 (ja) | 2020-05-22 | 2023-06-23 | ペキン シャオミ モバイル ソフトウェア カンパニー， リミテッド | 画像処理方法、端末デバイスにおけるイメージシグナルプロセッサ |
| WO2024057768A1 (ja) * | 2022-09-14 | 2024-03-21 | 富士フイルム株式会社 | 画像生成装置、学習装置、画像処理装置、画像生成方法、学習方法及び画像処理方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3995854B2 (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP3816151B2 (ja) | 画像処理方法および装置 | |
| US20100142790A1 (en) | Image processing method capable of enhancing contrast and reducing noise of digital image and image processing device using same | |
| JP4004562B2 (ja) | 画像処理方法および装置 | |
| US20070188785A1 (en) | Image processing apparatus, image processing method, storage medium, and program | |
| US6771793B1 (en) | Image processing method and apparatus | |
| JP3700804B2 (ja) | 画像処理方法および装置 | |
| JP3738788B2 (ja) | 画像のダイナミックレンジ圧縮処理方法および装置 | |
| JP2007018379A (ja) | 画像処理方法及び画像処理装置 | |
| JP4244094B2 (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP2001167264A (ja) | 画像処理方法および装置並びに記録媒体 | |
| JPH09212623A (ja) | 画像処理方法 | |
| JPH0944656A (ja) | 画像処理方法および装置 | |
| JP3700798B2 (ja) | 画像処理方法および装置 | |
| JP3738791B2 (ja) | 画像処理方法および装置 | |
| JP3999432B2 (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP4004561B2 (ja) | 画像処理方法および装置 | |
| JP2002074356A (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP2001056856A (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP3783116B2 (ja) | 放射線画像強調処理方法および装置 | |
| JP4014671B2 (ja) | 多重解像度変換方法および装置 | |
| JP2001167261A (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP3669452B2 (ja) | 画像処理方法および装置 | |
| JPH0944651A (ja) | 画像処理方法および装置 | |
| JP3870173B2 (ja) | 画像処理方法、画像処理装置、プログラム及びコンピュータ記録媒体 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050912 |
|
| A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20061201 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20080730 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080805 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20081118 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20090414 |