JP2002216126A - 画像処理装置、及びその方法、並びにプログラム、記憶媒体 - Google Patents
画像処理装置、及びその方法、並びにプログラム、記憶媒体Info
- Publication number
- JP2002216126A JP2002216126A JP2001328064A JP2001328064A JP2002216126A JP 2002216126 A JP2002216126 A JP 2002216126A JP 2001328064 A JP2001328064 A JP 2001328064A JP 2001328064 A JP2001328064 A JP 2001328064A JP 2002216126 A JP2002216126 A JP 2002216126A
- Authority
- JP
- Japan
- Prior art keywords
- threshold
- coefficient
- processing
- sub
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title description 32
- 238000012545 processing Methods 0.000 claims abstract description 107
- 238000006243 chemical reaction Methods 0.000 claims description 51
- 230000007850 degeneration Effects 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims 5
- 238000003672 processing method Methods 0.000 claims 3
- 230000009466 transformation Effects 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 11
- 238000003384 imaging method Methods 0.000 description 7
- 238000007781 pre-processing Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000013480 data collection Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 101150016011 RR11 gene Proteins 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/10—Image enhancement or restoration using non-spatial domain filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20064—Wavelet transform [DWT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Picture Signal Circuits (AREA)
- Closed-Circuit Television Systems (AREA)
- Compression, Expansion, Code Conversion, And Decoders (AREA)
- Image Analysis (AREA)
Abstract
(57)【要約】
【課題】 効率的にノイズ除去を行うこと。
【解決手段】 入力画像に対して離散ウェーブレット変
換を施し、各サブバンドのウェーブレット係数を出力す
る（Ｓ３０１）。高周波成分を示すＨＬ，ＬＨ，ＨＨの
サブバンドに対して、夫々のサブバンド毎に適当な閾値
を設定する(Ｓ３０２ａ，Ｓ３０２ｂ，Ｓ３０２ｃ)。次
に設定した各サブバンドの閾値を用いてＨＬ，ＬＨ，Ｈ
Ｈのウェーブレット係数に対して閾値処理を行う(Ｓ３
０３ａ，Ｓ３０３ｂ，Ｓ３０３ｃ)。次に各サブバンド
の閾値処理結果に基づき、係数変換処理の処理対象とす
る画素を決定する（Ｓ３０４）。次にＳ３０４で決定さ
れた処理対象画素に応じた各サブバンドのウェーブレッ
ト係数に対して係数変換を施し（Ｓ３０５）、この変換
係数に対して逆離散ウェーブレット変換を施し、ノイズ
の除去された画像を復元し、出力する(Ｓ３０６)。
換を施し、各サブバンドのウェーブレット係数を出力す
る（Ｓ３０１）。高周波成分を示すＨＬ，ＬＨ，ＨＨの
サブバンドに対して、夫々のサブバンド毎に適当な閾値
を設定する(Ｓ３０２ａ，Ｓ３０２ｂ，Ｓ３０２ｃ)。次
に設定した各サブバンドの閾値を用いてＨＬ，ＬＨ，Ｈ
Ｈのウェーブレット係数に対して閾値処理を行う(Ｓ３
０３ａ，Ｓ３０３ｂ，Ｓ３０３ｃ)。次に各サブバンド
の閾値処理結果に基づき、係数変換処理の処理対象とす
る画素を決定する（Ｓ３０４）。次にＳ３０４で決定さ
れた処理対象画素に応じた各サブバンドのウェーブレッ
ト係数に対して係数変換を施し（Ｓ３０５）、この変換
係数に対して逆離散ウェーブレット変換を施し、ノイズ
の除去された画像を復元し、出力する(Ｓ３０６)。
Description
【０００１】
【発明の属する技術分野】本発明は、画像を変更し、特
に周波数変換の際の変換係数を変更する画像処理装置、
及びその方法、並びにプログラム、記憶媒体に関するも
のである。
に周波数変換の際の変換係数を変更する画像処理装置、
及びその方法、並びにプログラム、記憶媒体に関するも
のである。
【０００２】
【従来の技術】近年のデジタル技術の進歩により、放射
線画像をデジタル画像信号に変換し、該デジタル画像信
号に対して画像処理を施し、ＣＲＴ等に表示、あるいは
プリント出力することが行われている。ところで、放射
線画像の撮影においては、患者に与える影響などから撮
影時のＸ線量は少ないことが望ましい。しかし、Ｘ線量
を少なくして撮影された画像には量子ノイズが多く含ま
れることが知られており、このような量子ノイズは診断
を行う際の妨害となる可能性がある。
線画像をデジタル画像信号に変換し、該デジタル画像信
号に対して画像処理を施し、ＣＲＴ等に表示、あるいは
プリント出力することが行われている。ところで、放射
線画像の撮影においては、患者に与える影響などから撮
影時のＸ線量は少ないことが望ましい。しかし、Ｘ線量
を少なくして撮影された画像には量子ノイズが多く含ま
れることが知られており、このような量子ノイズは診断
を行う際の妨害となる可能性がある。
【０００３】そのため、従来からこのようなノイズを除
去する処理が検討されており、例えば、単純なメディア
ンフィルタによるノイズ除去処理や、平滑化画像を用い
て高周波成分を抽出し、ノイズを除去する方法（いわゆ
るフィルタ処理）などが行われてきた。また近年、入力
画像を複数の周波数帯域に分割し、各周波数帯域ごとに
独立した処理を施すことによりノイズの除去を行う多重
周波数処理が検討されている。
去する処理が検討されており、例えば、単純なメディア
ンフィルタによるノイズ除去処理や、平滑化画像を用い
て高周波成分を抽出し、ノイズを除去する方法（いわゆ
るフィルタ処理）などが行われてきた。また近年、入力
画像を複数の周波数帯域に分割し、各周波数帯域ごとに
独立した処理を施すことによりノイズの除去を行う多重
周波数処理が検討されている。
【０００４】
【発明が解決しようとする課題】平滑化画像を用いて高
周波成分を抽出し、ノイズの除去を行うフィルタ処理で
は、単一の周波数帯域を扱うため、ノイズ成分が広い周
波数帯域に分布している場合には効果的にノイズ除去を
行えないという問題がある。これを避けるために、大き
さの異なる（つまり周波数の異なる）複数のフィルタを
同時に利用することも行われている。しかし、処理に必
要な計算コストが大幅に増加する問題がある。また、ノ
イズ除去のためにフィルタの周波数特性を最適に保つた
めには、被写体に応じたフィルタサイズの調整が必須と
なり、汎用性に乏しいという問題がある。
周波成分を抽出し、ノイズの除去を行うフィルタ処理で
は、単一の周波数帯域を扱うため、ノイズ成分が広い周
波数帯域に分布している場合には効果的にノイズ除去を
行えないという問題がある。これを避けるために、大き
さの異なる（つまり周波数の異なる）複数のフィルタを
同時に利用することも行われている。しかし、処理に必
要な計算コストが大幅に増加する問題がある。また、ノ
イズ除去のためにフィルタの周波数特性を最適に保つた
めには、被写体に応じたフィルタサイズの調整が必須と
なり、汎用性に乏しいという問題がある。
【０００５】また、ノイズ除去において離散ウェーブレ
ット変換などに代表される多重周波数処理を用いること
により上記の問題は大きく軽減されるが、一般に多重周
波数処理では空間的にすべての方向についての一律の周
波数処理を行うため、例えば方向成分を持つエッジと持
たない孤立点ノイズを区別することができず、エッジを
保存しながらノイズを除去することは困難であった。
ット変換などに代表される多重周波数処理を用いること
により上記の問題は大きく軽減されるが、一般に多重周
波数処理では空間的にすべての方向についての一律の周
波数処理を行うため、例えば方向成分を持つエッジと持
たない孤立点ノイズを区別することができず、エッジを
保存しながらノイズを除去することは困難であった。
【０００６】本発明は以上の問題に鑑みてなされたもの
であり、効率的にノイズ除去を行うことを目的とする。
であり、効率的にノイズ除去を行うことを目的とする。
【０００７】
【課題を解決するための手段】本発明の目的を達成する
ために、例えば本発明の画像処理装置は以下の構成を備
える。すなわち、画像に対して周波数変換を施すことで
複数のサブバンドを得る周波数変換手段と、前記周波数
変換手段による前記複数のサブバンドのうち少なくとも
２つ以上のサブバンドを用いて変更対象のサブバンド内
の変換係数を決定する決定手段と、前記決定手段により
決定した変換係数に対して変換を施す係数変換手段とを
備え、更に前記周波数変換手段は、前記係数変換手段に
より変換された変換係数を含む全てのサブバンドの変換
係数を用いて画像を復元する。
ために、例えば本発明の画像処理装置は以下の構成を備
える。すなわち、画像に対して周波数変換を施すことで
複数のサブバンドを得る周波数変換手段と、前記周波数
変換手段による前記複数のサブバンドのうち少なくとも
２つ以上のサブバンドを用いて変更対象のサブバンド内
の変換係数を決定する決定手段と、前記決定手段により
決定した変換係数に対して変換を施す係数変換手段とを
備え、更に前記周波数変換手段は、前記係数変換手段に
より変換された変換係数を含む全てのサブバンドの変換
係数を用いて画像を復元する。
【０００８】前記決定手段は更に、前記周波数変換手段
による前記複数のサブバンドのうち少なくとも２つ以上
のサブバンドに対して閾値を設定する閾値設定手段と、
前記閾値設定手段により閾値を設定したサブバンド内の
変換係数に対して当該閾値を用いて閾値処理を施す閾値
処理手段とを備え、前記閾値処理手段による閾値処理結
果を用いて変更対象のサブバンド内の変換係数を決定す
る。
による前記複数のサブバンドのうち少なくとも２つ以上
のサブバンドに対して閾値を設定する閾値設定手段と、
前記閾値設定手段により閾値を設定したサブバンド内の
変換係数に対して当該閾値を用いて閾値処理を施す閾値
処理手段とを備え、前記閾値処理手段による閾値処理結
果を用いて変更対象のサブバンド内の変換係数を決定す
る。
【０００９】
【発明の実施の形態】以下添付図面に従って、本発明を
好適な実施形態に従って詳細に説明する。
好適な実施形態に従って詳細に説明する。
【００１０】［第１の実施形態］図１は本実施形態にお
けるＸ線撮影装置１００の基本構成とＸ線により撮影す
る被写体を示す。Ｘ線撮影装置１００はＸ線により被写
体を撮影し、撮影された画像のノイズ除去を行う機能を
有する装置であり、前処理回路１０６、ＣＰＵ１０８、
メインメモリ１０９、操作パネル１１０、画像表示器１
１１、画像処理回路１１２を備える。なお画像処理回路
１１２は離散ウェーブレット（ＤＷＴ）回路１１６、係
数変換回路１１７、特徴抽出回路１１８、階調変換回路
１１９、閾値処理回路１２０、処理判定回路１２１を備
えており、ＣＰＵバス１０７を介して互いにデータ授受
されるようになされている。
けるＸ線撮影装置１００の基本構成とＸ線により撮影す
る被写体を示す。Ｘ線撮影装置１００はＸ線により被写
体を撮影し、撮影された画像のノイズ除去を行う機能を
有する装置であり、前処理回路１０６、ＣＰＵ１０８、
メインメモリ１０９、操作パネル１１０、画像表示器１
１１、画像処理回路１１２を備える。なお画像処理回路
１１２は離散ウェーブレット（ＤＷＴ）回路１１６、係
数変換回路１１７、特徴抽出回路１１８、階調変換回路
１１９、閾値処理回路１２０、処理判定回路１２１を備
えており、ＣＰＵバス１０７を介して互いにデータ授受
されるようになされている。
【００１１】またＸ線撮影装置１００は、前処理回路１
０６に接続されたデータ収集回路１０５と、データ収集
回路１０５に接続された２次元Ｘ線センサ１０４及びＸ
線発生回路１０１とを備えており、これらの各回路はＣ
ＰＵバス１０７にも接続されている。
０６に接続されたデータ収集回路１０５と、データ収集
回路１０５に接続された２次元Ｘ線センサ１０４及びＸ
線発生回路１０１とを備えており、これらの各回路はＣ
ＰＵバス１０７にも接続されている。
【００１２】図２は上述の構成を備える本実施形態にお
けるＸ線撮影装置が行う大まかな処理のフローチャート
である。なおメインメモリ１０９はＣＰＵ１０８での処
理に必要な各種のデータなどが記憶されるものであると
共に、ＣＰＵ１０８の作業用のワークメモリを含む。Ｃ
ＰＵ１０８はメインメモリ１０９を用いて操作パネル１
１０からの操作に従った装置全体の動作制御等を行う。
これによりＸ線撮影装置１００は、以下のように動作す
る。
けるＸ線撮影装置が行う大まかな処理のフローチャート
である。なおメインメモリ１０９はＣＰＵ１０８での処
理に必要な各種のデータなどが記憶されるものであると
共に、ＣＰＵ１０８の作業用のワークメモリを含む。Ｃ
ＰＵ１０８はメインメモリ１０９を用いて操作パネル１
１０からの操作に従った装置全体の動作制御等を行う。
これによりＸ線撮影装置１００は、以下のように動作す
る。
【００１３】先ずＸ線発生回路１０１は、被写体１０３
に対してＸ線ビーム１０２を放射する。Ｘ線発生回路１
０１から放射されたＸ線ビーム１０２は、被写体１０３
を減衰しながら透過して、２次元Ｘ線センサ１０４に到
達し、２次元Ｘ線センサ１０４によりＸ線画像として出
力される（ステップＳ２０１）。ここでは、２次元Ｘ線
センサ１０４から出力されるＸ線画像を、以下の説明で
は人体画像とするが、これに限定されるものではない。
に対してＸ線ビーム１０２を放射する。Ｘ線発生回路１
０１から放射されたＸ線ビーム１０２は、被写体１０３
を減衰しながら透過して、２次元Ｘ線センサ１０４に到
達し、２次元Ｘ線センサ１０４によりＸ線画像として出
力される（ステップＳ２０１）。ここでは、２次元Ｘ線
センサ１０４から出力されるＸ線画像を、以下の説明で
は人体画像とするが、これに限定されるものではない。
【００１４】データ収集回路１０５は２次元Ｘ線センサ
１０４から出力されたＸ線画像を電気信号に変換して前
処理回路１０６に供給する。前処理回路１０６はデータ
収集回路１０５からの信号（Ｘ線画像信号）に対して、
オフセット補正処理やゲイン補正処理等の前処理を行う
（ステップＳ２０２）。この前処理回路１０６で前処理
が行われたＸ線画像信号は原画像として、ＣＰＵ１０８
の制御により、ＣＰＵバス１０７を介してメインメモリ
１０９、画像処理回路１１２に転送され、後述のノイズ
除去処理（ステップＳ２０３）が行われる。
１０４から出力されたＸ線画像を電気信号に変換して前
処理回路１０６に供給する。前処理回路１０６はデータ
収集回路１０５からの信号（Ｘ線画像信号）に対して、
オフセット補正処理やゲイン補正処理等の前処理を行う
（ステップＳ２０２）。この前処理回路１０６で前処理
が行われたＸ線画像信号は原画像として、ＣＰＵ１０８
の制御により、ＣＰＵバス１０７を介してメインメモリ
１０９、画像処理回路１１２に転送され、後述のノイズ
除去処理（ステップＳ２０３）が行われる。
【００１５】画像処理回路１１２において離散ウェーブ
レット変換回路１１６は入力画像に対して離散ウェーブ
レット変換を施し、サブバンドごとにウェーブレット係
数を出力する。閾値処理回路１２０は離散ウェーブレッ
ト係数に対してサブバンドごとに適当な閾値を設定し、
設定した閾値で閾値処理する。処理判定回路１２１は各
サブバンドの閾値処理結果に基づき各画素の離散ウェー
ブレット係数を処理するかどうか判定する。係数変換回
路１１７は離散ウェーブレット係数を後述の規則に従い
変換する。また画像処理回路１１２は、階調変換に必要
な特徴量を抽出する特徴量抽出回路１１８と、得られた
特徴量に応じて階調変換を行う階調変換回路１１９を備
える。
レット変換回路１１６は入力画像に対して離散ウェーブ
レット変換を施し、サブバンドごとにウェーブレット係
数を出力する。閾値処理回路１２０は離散ウェーブレッ
ト係数に対してサブバンドごとに適当な閾値を設定し、
設定した閾値で閾値処理する。処理判定回路１２１は各
サブバンドの閾値処理結果に基づき各画素の離散ウェー
ブレット係数を処理するかどうか判定する。係数変換回
路１１７は離散ウェーブレット係数を後述の規則に従い
変換する。また画像処理回路１１２は、階調変換に必要
な特徴量を抽出する特徴量抽出回路１１８と、得られた
特徴量に応じて階調変換を行う階調変換回路１１９を備
える。
【００１６】図３は画像処理回路１１２においてノイズ
除去に関する部分についての処理の流れを示すフローチ
ャートであり、図４は同処理の様子を示す図である。こ
れらの図を用いてノイズ除去処理の流れを説明する。
除去に関する部分についての処理の流れを示すフローチ
ャートであり、図４は同処理の様子を示す図である。こ
れらの図を用いてノイズ除去処理の流れを説明する。
【００１７】前処理回路１０６で前処理された原画像
（図４（ａ）に示す画像４０１）はＣＰＵバス１０７を
介して画像処理装置１１２に転送される。画像処理装置
１１２では、まず離散ウェーブレット変換回路１１６が
入力画像に対して離散ウェーブレット変換を施し、各サ
ブバンドの離散ウェーブレット係数を出力する（ステッ
プＳ３０１）。上述の通り図４（ａ）において４０１は
入力された画像であり、（ｂ）において４０２は１レベ
ルの離散ウェーブレット変換が施された画像である。ま
た画像４０２において４０３から４０６の領域は、一般
にそれぞれＬＬ，ＨＬ，ＬＨ，ＨＨと呼ばれるサブバン
ドである。離散ウェーブレット変換の詳細については後
述する。
（図４（ａ）に示す画像４０１）はＣＰＵバス１０７を
介して画像処理装置１１２に転送される。画像処理装置
１１２では、まず離散ウェーブレット変換回路１１６が
入力画像に対して離散ウェーブレット変換を施し、各サ
ブバンドの離散ウェーブレット係数を出力する（ステッ
プＳ３０１）。上述の通り図４（ａ）において４０１は
入力された画像であり、（ｂ）において４０２は１レベ
ルの離散ウェーブレット変換が施された画像である。ま
た画像４０２において４０３から４０６の領域は、一般
にそれぞれＬＬ，ＨＬ，ＬＨ，ＨＨと呼ばれるサブバン
ドである。離散ウェーブレット変換の詳細については後
述する。
【００１８】次に、ステップＳ３０１で得られたサブバ
ンドのうち高周波成分を示すＨＬ，ＬＨ，ＨＨの３つの
サブバンドに対して、それぞれのサブバンドごとに適当
な閾値を設定する(ステップＳ３０２ａ，Ｓ３０２ｂ，
Ｓ３０２ｃ)。閾値の設定方法は特に限定しないが、離
散ウェーブレット変換の分解レベルに応じて経験的に求
めたある定数を設定しても良いし、各サブバンドにおけ
る係数値の平均値や分散値などの統計量をもとに決めて
も良い。
ンドのうち高周波成分を示すＨＬ，ＬＨ，ＨＨの３つの
サブバンドに対して、それぞれのサブバンドごとに適当
な閾値を設定する(ステップＳ３０２ａ，Ｓ３０２ｂ，
Ｓ３０２ｃ)。閾値の設定方法は特に限定しないが、離
散ウェーブレット変換の分解レベルに応じて経験的に求
めたある定数を設定しても良いし、各サブバンドにおけ
る係数値の平均値や分散値などの統計量をもとに決めて
も良い。
【００１９】次にステップＳ３０２ａ，Ｓ３０２ｂ，Ｓ
３０２ｃの夫々で設定した各サブバンドの閾値を用い
て、ＨＬ，ＬＨ，ＨＨのウェーブレット係数に対して閾
値処理を行う(ステップＳ３０３ａ，Ｓ３０３ｂ，Ｓ３
０３ｃ)。図４（ｃ）において４０７はステップＳ３０
１で得られた離散ウェーブレット変換後の画像であり、
４０８はステップＳ３０２の閾値処理を施した画像であ
る。このとき画像４０７と画像４０８は別のメモリ空間
に保持するものとする。閾値処理の方法としては、例え
ば符号付２値化処理が考えられ、以下のように実現でき
る。 if (-THHL ＜ HL(x,y) ＜ THHL) then ２値化画像HL(x,y) = 1; else ２値化画像HL(x,y) = 0; if (-THLH ＜ LH(x,y) ＜ THLH) then ２値化画像LH(x,y) = 1; else ２値化画像LH(x,y) = 0; if (-THHH ＜ HH(x,y) ＜ THHH) then ２値化画像HH(x,y) = 1; else ２値化画像HH(x,y) = 0; ここで、THHL,THLH,THHHは夫々ステップＳ３０２ａ，Ｓ
３０２ｂ，Ｓ３０２ｃで設定された各サブバンドの閾値
であり、HL(x,y),LH(x,y),HH(x,y)は、各サブバンド内
の位置（x,y）における画素値、すなわちウェーブレッ
ト係数値を表す。
３０２ｃの夫々で設定した各サブバンドの閾値を用い
て、ＨＬ，ＬＨ，ＨＨのウェーブレット係数に対して閾
値処理を行う(ステップＳ３０３ａ，Ｓ３０３ｂ，Ｓ３
０３ｃ)。図４（ｃ）において４０７はステップＳ３０
１で得られた離散ウェーブレット変換後の画像であり、
４０８はステップＳ３０２の閾値処理を施した画像であ
る。このとき画像４０７と画像４０８は別のメモリ空間
に保持するものとする。閾値処理の方法としては、例え
ば符号付２値化処理が考えられ、以下のように実現でき
る。 if (-THHL ＜ HL(x,y) ＜ THHL) then ２値化画像HL(x,y) = 1; else ２値化画像HL(x,y) = 0; if (-THLH ＜ LH(x,y) ＜ THLH) then ２値化画像LH(x,y) = 1; else ２値化画像LH(x,y) = 0; if (-THHH ＜ HH(x,y) ＜ THHH) then ２値化画像HH(x,y) = 1; else ２値化画像HH(x,y) = 0; ここで、THHL,THLH,THHHは夫々ステップＳ３０２ａ，Ｓ
３０２ｂ，Ｓ３０２ｃで設定された各サブバンドの閾値
であり、HL(x,y),LH(x,y),HH(x,y)は、各サブバンド内
の位置（x,y）における画素値、すなわちウェーブレッ
ト係数値を表す。
【００２０】次に、ステップＳ３０３ａ，Ｓ３０３ｂ，
Ｓ３０３ｃによる各サブバンドの閾値処理結果に基づ
き、後段の係数変換処理の処理対象とする画素を決定す
る（ステップＳ３０４）。処理対象画素の決定方法はい
くつか考えられるが、本実施形態ではステップＳ３０３
ａ，Ｓ３０３ｂ，Ｓ３０３ｃによる各サブバンドの２値
化画像の論理積を採用する。即ち以下のように実現され
る。 if ((２値化画像HL(x,y) = 1) AND (２値化画像LH(x,y) = 1) AND (２値化画像HH(x,y) = 1) ) then ２値化画像HL(x,y) = ２値化画像LH(x,y) = ２値
化画像HH(x,y) = 1;else ２値化画像HL(x,y) = ２値化
画像LH(x,y) = ２値化画像HH(x,y) = 0;図４（ｄ）にお
いて４０９はステップＳ３０１で得られた離散ウェーブ
レット変換後の画像であり、４１０はステップＳ３０４
の処理対象画素を示す画像である。このとき画像４１０
は同図（ｃ）に示す画像４０８が格納されるメモリと同
じメモリ空間に保持しても構わない。また、処理対象画
素の決定方法は前記の他にも、各サブバンドの閾値処理
結果を多数決する方法なども考えられる。
Ｓ３０３ｃによる各サブバンドの閾値処理結果に基づ
き、後段の係数変換処理の処理対象とする画素を決定す
る（ステップＳ３０４）。処理対象画素の決定方法はい
くつか考えられるが、本実施形態ではステップＳ３０３
ａ，Ｓ３０３ｂ，Ｓ３０３ｃによる各サブバンドの２値
化画像の論理積を採用する。即ち以下のように実現され
る。 if ((２値化画像HL(x,y) = 1) AND (２値化画像LH(x,y) = 1) AND (２値化画像HH(x,y) = 1) ) then ２値化画像HL(x,y) = ２値化画像LH(x,y) = ２値
化画像HH(x,y) = 1;else ２値化画像HL(x,y) = ２値化
画像LH(x,y) = ２値化画像HH(x,y) = 0;図４（ｄ）にお
いて４０９はステップＳ３０１で得られた離散ウェーブ
レット変換後の画像であり、４１０はステップＳ３０４
の処理対象画素を示す画像である。このとき画像４１０
は同図（ｃ）に示す画像４０８が格納されるメモリと同
じメモリ空間に保持しても構わない。また、処理対象画
素の決定方法は前記の他にも、各サブバンドの閾値処理
結果を多数決する方法なども考えられる。
【００２１】次に係数変換回路１１７がステップＳ３０
４で決定された処理対象画素に応じた各サブバンドのウ
ェーブレット係数に対して係数変換を施す（ステップＳ
３０５）。ここでステップＳ３０５の処理の一例とし
て、一般によく知られるウェーブレット縮退とよばれる
ノイズ除去処理などが考えられる。ウェーブレット縮退
とは、各サブバンドの係数をある条件に基づいて抑制す
るものであり、最も簡単な例としては、係数値の絶対値
がある範囲に存在する場合、下記に示すようにその係数
を０で置き換える処理が行われる。 if (TH1 ≦ |HH(x,y)| ≦ TH2) then HH(x,y) = 0; else DO NOTHING; if (TH3 ≦ |HL(x,y)| ≦ TH4) then HL(x,y) = 0; else DO NOTHING; if (TH5 ≦ |LH(x,y)| ≦ TH6) then LH(x,y) = 0; else DO NOTHING; ここで、TH1〜6は別に定める閾値であり、 0 ≦ TH1 ≦ TH2 0 ≦ TH3 ≦ TH4 0 ≦ TH5 ≦ TH6 の条件を満たす処理パラメータである。
４で決定された処理対象画素に応じた各サブバンドのウ
ェーブレット係数に対して係数変換を施す（ステップＳ
３０５）。ここでステップＳ３０５の処理の一例とし
て、一般によく知られるウェーブレット縮退とよばれる
ノイズ除去処理などが考えられる。ウェーブレット縮退
とは、各サブバンドの係数をある条件に基づいて抑制す
るものであり、最も簡単な例としては、係数値の絶対値
がある範囲に存在する場合、下記に示すようにその係数
を０で置き換える処理が行われる。 if (TH1 ≦ |HH(x,y)| ≦ TH2) then HH(x,y) = 0; else DO NOTHING; if (TH3 ≦ |HL(x,y)| ≦ TH4) then HL(x,y) = 0; else DO NOTHING; if (TH5 ≦ |LH(x,y)| ≦ TH6) then LH(x,y) = 0; else DO NOTHING; ここで、TH1〜6は別に定める閾値であり、 0 ≦ TH1 ≦ TH2 0 ≦ TH3 ≦ TH4 0 ≦ TH5 ≦ TH6 の条件を満たす処理パラメータである。
【００２２】最後に、ステップＳ３０５で係数変換回路
１１７によって変換された各サブバンドのウェーブレッ
ト係数を用いて、離散ウェーブレット変換回路１１６が
逆離散ウェーブレット変換を施し、ノイズの除去された
画像を復元し、出力する(ステップＳ３０６)。
１１７によって変換された各サブバンドのウェーブレッ
ト係数を用いて、離散ウェーブレット変換回路１１６が
逆離散ウェーブレット変換を施し、ノイズの除去された
画像を復元し、出力する(ステップＳ３０６)。
【００２３】以上図３及び図４を用いて説明した本実施
形態におけるノイズ除去処理では１レベルの離散ウェー
ブレット変換を扱ったが、これに限定されるものではな
く、さらに高レベルのウェーブレット係数に対して同様
の処理を行うことも可能であり、また複数のレベルのウ
ェーブレット係数に対して同様の処理を行うことももち
ろん可能である。
形態におけるノイズ除去処理では１レベルの離散ウェー
ブレット変換を扱ったが、これに限定されるものではな
く、さらに高レベルのウェーブレット係数に対して同様
の処理を行うことも可能であり、また複数のレベルのウ
ェーブレット係数に対して同様の処理を行うことももち
ろん可能である。
【００２４】離散ウェーブレット変換回路１１６による
離散ウェーブレット変換及び逆離散ウェーブレット変換
は一般によく知られる変換処理であり以下のように動作
する。
離散ウェーブレット変換及び逆離散ウェーブレット変換
は一般によく知られる変換処理であり以下のように動作
する。
【００２５】離散ウェーブレット変換回路１１６は入力
した画像信号に対して２次元の離散ウェーブレット変換
処理を行い、変換係数を計算して出力するものである。
メインメモリ１０９に記憶された画像データは、離散ウ
ェーブレット変換回路１１６により順次読み出されて変
換処理が行われ、再びメインメモリ１０９に書きこまれ
る。本実施形態における離散ウェーブレット変換回路１
１６の処理の構成は図５（ａ）に示すものとする。同図
において、入力された画像信号は遅延素子およびダウン
サンプラの組み合わせにより、偶数アドレスおよび奇数
アドレスの信号に分離され、２つのフィルタｐおよびu
によりフィルタ処理が施される。同図sおよびdは、各々
1次元の画像信号に対して１レベルの分解を行った際の
ローパス係数およびハイパス係数を表しており、次式に
より計算されるものとする。 d(n) = x(2n+1) - floor((x(2n) + x(2n+2))/2) （式1） s(n) = x(2n) + floor((d(n-1) + d(n))/4) （式2） ただし、x(n)は変換対象となる画像信号である。以上の
処理により、画像信号に対する１次元の離散ウェーブレ
ット変換処理が行われる。２次元の離散ウェーブレット
変換は、１次元の変換を画像の水平・垂直方向に対して
順次行うものであり、その詳細は公知であるのでここで
は説明を省略する。図５（ｂ）は２次元の変換処理によ
り得られる２レベルの変換係数群の構成例であり、画像
信号は異なる周波数帯域の係数列ＨＨ１，ＨＬ１，ＬＨ
１，．．．，ＨＨ２，ＨＬ２，ＬＨ２，ＬＬに分解され
る。
した画像信号に対して２次元の離散ウェーブレット変換
処理を行い、変換係数を計算して出力するものである。
メインメモリ１０９に記憶された画像データは、離散ウ
ェーブレット変換回路１１６により順次読み出されて変
換処理が行われ、再びメインメモリ１０９に書きこまれ
る。本実施形態における離散ウェーブレット変換回路１
１６の処理の構成は図５（ａ）に示すものとする。同図
において、入力された画像信号は遅延素子およびダウン
サンプラの組み合わせにより、偶数アドレスおよび奇数
アドレスの信号に分離され、２つのフィルタｐおよびu
によりフィルタ処理が施される。同図sおよびdは、各々
1次元の画像信号に対して１レベルの分解を行った際の
ローパス係数およびハイパス係数を表しており、次式に
より計算されるものとする。 d(n) = x(2n+1) - floor((x(2n) + x(2n+2))/2) （式1） s(n) = x(2n) + floor((d(n-1) + d(n))/4) （式2） ただし、x(n)は変換対象となる画像信号である。以上の
処理により、画像信号に対する１次元の離散ウェーブレ
ット変換処理が行われる。２次元の離散ウェーブレット
変換は、１次元の変換を画像の水平・垂直方向に対して
順次行うものであり、その詳細は公知であるのでここで
は説明を省略する。図５（ｂ）は２次元の変換処理によ
り得られる２レベルの変換係数群の構成例であり、画像
信号は異なる周波数帯域の係数列ＨＨ１，ＨＬ１，ＬＨ
１，．．．，ＨＨ２，ＨＬ２，ＬＨ２，ＬＬに分解され
る。
【００２６】同様に逆離散ウェーブレット変換について
は次のように行われる。メインメモリ１０９に記憶され
た変換係数は離散ウェーブレット変換回路１１６により
順次読み出されて変換処理が行われ、再びメインメモリ
１０９に書きこまれる。本実施形態における離散ウェー
ブレット変換回路１１６による逆離散ウェーブレット変
換処理の構成は図５（ｃ）に示すものとする。入力され
た変換係数はuおよびpの２つのフィルタを用いたフィル
タ処理が施され、アップサンプリングされた後に重ね合
わされて画像信号x'が出力される。これらの処理は次式
により行われる。 x'(2n) = s'(n) - floor ((d'(n-1) + d'(n))/4) （式3） x'(2n+1) = d'(n) + floor ((x'(2n) + x'(2n+2))/2) （式4） 以上の処理により、変換係数に対する１次元の逆離散ウ
ェーブレット変換処理が行われる。２次元の逆離散ウェ
ーブレット変換は、１次元の逆変換を画像の水平・垂直
方向に対して順次行うものであり、その詳細は公知であ
るのでここでは説明を省略する。
は次のように行われる。メインメモリ１０９に記憶され
た変換係数は離散ウェーブレット変換回路１１６により
順次読み出されて変換処理が行われ、再びメインメモリ
１０９に書きこまれる。本実施形態における離散ウェー
ブレット変換回路１１６による逆離散ウェーブレット変
換処理の構成は図５（ｃ）に示すものとする。入力され
た変換係数はuおよびpの２つのフィルタを用いたフィル
タ処理が施され、アップサンプリングされた後に重ね合
わされて画像信号x'が出力される。これらの処理は次式
により行われる。 x'(2n) = s'(n) - floor ((d'(n-1) + d'(n))/4) （式3） x'(2n+1) = d'(n) + floor ((x'(2n) + x'(2n+2))/2) （式4） 以上の処理により、変換係数に対する１次元の逆離散ウ
ェーブレット変換処理が行われる。２次元の逆離散ウェ
ーブレット変換は、１次元の逆変換を画像の水平・垂直
方向に対して順次行うものであり、その詳細は公知であ
るのでここでは説明を省略する。
【００２７】以上説明したように本実施形態によれば、
離散ウェーブレット変換による多重周波数処理を行うた
め、単一の周波数帯域を扱ういわゆるフィルタ処理と比
べて、ノイズ成分が広い周波数帯域に分布している場合
にも効果的にノイズ除去を行える効果がある。また、大
きさの異なる（つまり周波数の異なる）複数のフィルタ
を同時に利用する処理に比べて、処理に必要な計算コス
トを大幅に低減する効果が有り、フィルタサイズの調整
など煩雑な作業を行わずに汎用性の高い処理が行える効
果がある。また多重周波数処理において高周波成分の空
間的な方向を考慮した処理を行っているため、例えば方
向成分を持つエッジと持たない孤立点ノイズを区別する
ことが可能となる効果があり、結果として従来手法に比
べて高性能にノイズ除去された画像を得ることが可能と
なる効果がある。
離散ウェーブレット変換による多重周波数処理を行うた
め、単一の周波数帯域を扱ういわゆるフィルタ処理と比
べて、ノイズ成分が広い周波数帯域に分布している場合
にも効果的にノイズ除去を行える効果がある。また、大
きさの異なる（つまり周波数の異なる）複数のフィルタ
を同時に利用する処理に比べて、処理に必要な計算コス
トを大幅に低減する効果が有り、フィルタサイズの調整
など煩雑な作業を行わずに汎用性の高い処理が行える効
果がある。また多重周波数処理において高周波成分の空
間的な方向を考慮した処理を行っているため、例えば方
向成分を持つエッジと持たない孤立点ノイズを区別する
ことが可能となる効果があり、結果として従来手法に比
べて高性能にノイズ除去された画像を得ることが可能と
なる効果がある。
【００２８】［第２の実施形態］なお、第１の実施形態
におけるステップＳ３０２〜Ｓ３０４では、ステップＳ
３０５における処理対象画素を決定するために閾値処理
を行っているものとして説明したが、これに限定される
ものではなく、本実施形態では３つのサブバンドの大小
関係を比較して、その大小のパターンに応じて、ステッ
プＳ３０５における処理対象画素を決定する。即ち例え
ば、ランダムノイズがＨＨサブバンドにより多く現れる
という傾向を利用して、ステップＳ３０２〜Ｓ３０４の
処理を以下のようにしても良い。 if ( (HL(x,y) ＜ HH(x,y)) AND (LH(x,y) ＜ HH(x,y)) ) then ２値化画像HL(x,y) = ２値化画像LH(x,y) = ２値化画像HH(x,y) = 1; else ２値化画像HL(x,y) = ２値化画像LH(x,y) = ２値化画像HH(x,y) = 0; また、閾値処理と大小比較を組み合わせて、ステップＳ
３０２〜Ｓ３０４の処理を以下のようにしても上述した
ものと同様の目的を達成できる。 if (-THHH ＜ HH(x,y) ＜ THHH) then ２値化画像HH(x,y) = 1; else ２値化画像HH(x,y) = 0; if ((HL(x,y) ＜ HH(x,y)) AND (LH(x,y) ＜ HH(x,y)) AND (２値化画像HH(x,y) = 1) ) then ２値化画像HL(x,y) = ２値化画像LH(x,y) = ２値化画像HH(x,y) = 1; else ２値化画像HL(x,y) = ２値化画像LH(x,y) = ２値化画像HH(x,y) = 0; さらにステップＳ３０２〜Ｓ３０４の処理において、閾
値からの距離による重み付け評価を導入したり、各サブ
バンドの平均値を閾値処理する方法なども考えられる
が、各サブバンドの係数を相互に参照する方法について
は特に限定しない。
におけるステップＳ３０２〜Ｓ３０４では、ステップＳ
３０５における処理対象画素を決定するために閾値処理
を行っているものとして説明したが、これに限定される
ものではなく、本実施形態では３つのサブバンドの大小
関係を比較して、その大小のパターンに応じて、ステッ
プＳ３０５における処理対象画素を決定する。即ち例え
ば、ランダムノイズがＨＨサブバンドにより多く現れる
という傾向を利用して、ステップＳ３０２〜Ｓ３０４の
処理を以下のようにしても良い。 if ( (HL(x,y) ＜ HH(x,y)) AND (LH(x,y) ＜ HH(x,y)) ) then ２値化画像HL(x,y) = ２値化画像LH(x,y) = ２値化画像HH(x,y) = 1; else ２値化画像HL(x,y) = ２値化画像LH(x,y) = ２値化画像HH(x,y) = 0; また、閾値処理と大小比較を組み合わせて、ステップＳ
３０２〜Ｓ３０４の処理を以下のようにしても上述した
ものと同様の目的を達成できる。 if (-THHH ＜ HH(x,y) ＜ THHH) then ２値化画像HH(x,y) = 1; else ２値化画像HH(x,y) = 0; if ((HL(x,y) ＜ HH(x,y)) AND (LH(x,y) ＜ HH(x,y)) AND (２値化画像HH(x,y) = 1) ) then ２値化画像HL(x,y) = ２値化画像LH(x,y) = ２値化画像HH(x,y) = 1; else ２値化画像HL(x,y) = ２値化画像LH(x,y) = ２値化画像HH(x,y) = 0; さらにステップＳ３０２〜Ｓ３０４の処理において、閾
値からの距離による重み付け評価を導入したり、各サブ
バンドの平均値を閾値処理する方法なども考えられる
が、各サブバンドの係数を相互に参照する方法について
は特に限定しない。
【００２９】［他の実施形態］なお、本発明は、複数の
機器（例えばホストコンピュータ、インタフェイス機
器、リーダ、プリンタなど）から構成されるシステムに
適用しても、一つの機器からなる装置（例えば、複写
機、ファクシミリ装置など）に適用してもよい。
機器（例えばホストコンピュータ、インタフェイス機
器、リーダ、プリンタなど）から構成されるシステムに
適用しても、一つの機器からなる装置（例えば、複写
機、ファクシミリ装置など）に適用してもよい。
【００３０】また、本発明の目的は、前述した実施形態
の機能を実現するソフトウェアのプログラムコードを記
録した記憶媒体（または記録媒体）を、システムあるい
は装置に供給し、そのシステムあるいは装置のコンピュ
ータ（またはCPUやMPU）が記憶媒体に格納されたプログ
ラムコードを読み出し実行することによっても、達成さ
れることは言うまでもない。この場合、記憶媒体から読
み出されたプログラムコード自体が前述した実施形態の
機能を実現することになり、そのプログラムコードを記
憶した記憶媒体は本発明を構成することになる。また、
コンピュータが読み出したプログラムコードを実行する
ことにより、前述した実施形態の機能が実現されるだけ
でなく、そのプログラムコードの指示に基づき、コンピ
ュータ上で稼働しているオペレーティングシステム(OS)
などが実際の処理の一部または全部を行い、その処理に
よって前述した実施形態の機能が実現される場合も含ま
れることは言うまでもない。
の機能を実現するソフトウェアのプログラムコードを記
録した記憶媒体（または記録媒体）を、システムあるい
は装置に供給し、そのシステムあるいは装置のコンピュ
ータ（またはCPUやMPU）が記憶媒体に格納されたプログ
ラムコードを読み出し実行することによっても、達成さ
れることは言うまでもない。この場合、記憶媒体から読
み出されたプログラムコード自体が前述した実施形態の
機能を実現することになり、そのプログラムコードを記
憶した記憶媒体は本発明を構成することになる。また、
コンピュータが読み出したプログラムコードを実行する
ことにより、前述した実施形態の機能が実現されるだけ
でなく、そのプログラムコードの指示に基づき、コンピ
ュータ上で稼働しているオペレーティングシステム(OS)
などが実際の処理の一部または全部を行い、その処理に
よって前述した実施形態の機能が実現される場合も含ま
れることは言うまでもない。
【００３１】さらに、記憶媒体から読み出されたプログ
ラムコードが、コンピュータに挿入された機能拡張カー
ドやコンピュータに接続された機能拡張ユニットに備わ
るメモリに書込まれた後、そのプログラムコードの指示
に基づき、その機能拡張カードや機能拡張ユニットに備
わるCPUなどが実際の処理の一部または全部を行い、そ
の処理によって前述した実施形態の機能が実現される場
合も含まれることは言うまでもない。
ラムコードが、コンピュータに挿入された機能拡張カー
ドやコンピュータに接続された機能拡張ユニットに備わ
るメモリに書込まれた後、そのプログラムコードの指示
に基づき、その機能拡張カードや機能拡張ユニットに備
わるCPUなどが実際の処理の一部または全部を行い、そ
の処理によって前述した実施形態の機能が実現される場
合も含まれることは言うまでもない。
【００３２】本発明を上記記憶媒体に適用する場合、そ
の記憶媒体には、先に説明した（図２，及び／又は図３
に示す）フローチャートに対応するプログラムコードが
格納されることになる。
の記憶媒体には、先に説明した（図２，及び／又は図３
に示す）フローチャートに対応するプログラムコードが
格納されることになる。
【００３３】
【発明の効果】以上の説明により本発明によれば効率的
にノイズ除去を行うことができる。
にノイズ除去を行うことができる。
【図１】本発明の第１の実施形態におけるＸ線撮影装置
の基本構成と、Ｘ線により撮影する被写体を示す図であ
る。
の基本構成と、Ｘ線により撮影する被写体を示す図であ
る。
【図２】本発明の第１の実施形態におけるＸ線撮影装置
が行う大まかな処理のフローチャートである。
が行う大まかな処理のフローチャートである。
【図３】画像処理回路１１２においてノイズ除去に関す
る部分についての処理の流れを示すフローチャートであ
る。
る部分についての処理の流れを示すフローチャートであ
る。
【図４】図３に示すフローチャートに従った処理の様子
を示す図である。
を示す図である。
【図５】（ａ）は離散ウェーブレット変換回路１１６の
処理の構成を示す図であり、（ｂ）は２次元の変換処理
により得られる２レベルの変換係数群の構成例を示す図
であり、（ｃ）は離散ウェーブレット変換回路１１６に
よる逆離散ウェーブレット変換処理の構成を示す図であ
る。
処理の構成を示す図であり、（ｂ）は２次元の変換処理
により得られる２レベルの変換係数群の構成例を示す図
であり、（ｃ）は離散ウェーブレット変換回路１１６に
よる逆離散ウェーブレット変換処理の構成を示す図であ
る。
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.7 識別記号 ＦＩ テーマコート゛(参考） Ｈ０４Ｎ 7/18 Ａ６１Ｂ 6/00 ３５０Ｎ ５Ｊ０６４ Ｆターム(参考） 4C093 CA06 FD04 FD05 FD09 FF03 FF08 FF18 5B057 AA08 BA03 BA30 CE02 CE06 CG05 CH01 CH11 5C021 YA01 5C054 AA01 AA06 EJ05 HA12 5C077 LL02 MP01 PP01 PP46 PQ12 PQ18 PQ22 RR11 5J064 AA00 BA16 BB07 BC01 BC22 BC29 BD04
Claims (15)
- 【請求項１】 画像に対して周波数変換を施すことで複
数のサブバンドを得る周波数変換手段と、 前記周波数変換手段による前記複数のサブバンドのうち
少なくとも２つ以上のサブバンドを用いて変更対象のサ
ブバンド内の変換係数を決定する決定手段と、前記決定
手段により決定した変換係数に対して変換を施す係数変
換手段とを備え、 更に前記周波数変換手段は、前記係数変換手段により変
換された変換係数を含む全てのサブバンドの変換係数を
用いて画像を復元することを特徴とする画像処理装置。 - 【請求項２】 前記決定手段は更に、 前記周波数変換手段による前記複数のサブバンドのうち
少なくとも２つ以上のサブバンドに対して閾値を設定す
る閾値設定手段と、 前記閾値設定手段により閾値を設定したサブバンド内の
変換係数に対して当該閾値を用いて閾値処理を施す閾値
処理手段とを備え、 前記閾値処理手段による閾値処理結果を用いて変更対象
のサブバンド内の変換係数を決定することを特徴とする
請求項１に記載の画像処理装置。 - 【請求項３】 前記閾値設定手段は前記閾値を、各サブ
バンドにおける変換係数値の平均値や分散値を含む統計
量を元に設定することを特徴とする請求項２に記載の画
像処理装置。 - 【請求項４】 前記閾値設定手段は所定のサブバンド内
の変換係数を閾値として設定し、 前記閾値処理手段は前記閾値を設定したサブバンド内の
変換係数に対して、当該変換係数に空間的に対応する前
記所定のサブバンド内の変換係数を用いて閾値処理を施
すことを特徴とする請求項２に記載の画像処理装置。 - 【請求項５】 前記所定のサブバンド内の変換係数は所
定の閾値により閾値処理されたものであることを特徴と
する請求項４に記載の画像処理装置。 - 【請求項６】 前記閾値処理は２値化処理であることを
特徴とする請求項２乃至５のいずれか１項に記載の画像
処理装置。 - 【請求項７】 前記係数変換手段はウェーブレット縮退
によりノイズ除去を行うことを特徴とする請求項１に記
載の画像処理装置。 - 【請求項８】 前記周波数変換手段は離散ウェーブレッ
ト変換、もしくは逆離散ウェーブレット変換のいずれか
を実行することを特徴とする請求項１乃至７のいずれか
１項に記載の画像処理装置。 - 【請求項９】 画像に対して周波数変換を施すことで複
数のサブバンドを得る周波数変換工程と、 前記周波数変換工程による前記複数のサブバンドのうち
少なくとも２つ以上のサブバンドを用いて変更対象のサ
ブバンド内の変換係数を決定する決定工程と、 前記決定工程で決定した変換係数に対して変換を施す係
数変換工程とを備え、 更に前記周波数変換工程では、前記係数変換工程で変換
された変換係数を含む全てのサブバンドの変換係数を用
いて画像を復元することを特徴とする画像処理方法。 - 【請求項１０】 前記決定工程は更に、 前記周波数変換工程による前記複数のサブバンドのうち
少なくとも２つ以上のサブバンドに対して閾値を設定す
る閾値設定工程と、 前記閾値設定工程により閾値を設定したサブバンド内の
変換係数に対して当該閾値を用いて閾値処理を施す閾値
処理工程とを備え、 前記閾値処理工程による閾値処理結果を用いて変更対象
のサブバンド内の変換係数を決定することを特徴とする
請求項９に記載の画像処理方法。 - 【請求項１１】 前記閾値設定工程では所定のサブバン
ド内の変換係数を閾値として設定し、 前記閾値処理工程では前記閾値を設定したサブバンド内
の変換係数に対して、当該変換係数に空間的に対応する
前記所定のサブバンド内の変換係数を用いて閾値処理を
施すことを特徴とする請求項１０に記載の画像処理方
法。 - 【請求項１２】 画像に対して周波数変換を施すことで
複数のサブバンドを得る周波数変換工程のプログラム
と、 前記周波数変換工程による前記複数のサブバンドのうち
少なくとも２つ以上のサブバンドを用いて変更対象のサ
ブバンド内の変換係数を決定する決定工程のプログラム
と、 前記決定工程で決定した変換係数に対して変換を施す係
数変換工程のプログラムとを備え、 更に前記周波数変換工程では、前記係数変換工程で変換
された変換係数を含む全てのサブバンドの変換係数を用
いて画像を復元することを特徴とするプログラム。 - 【請求項１３】 前記決定工程のプログラムは更に、 前記周波数変換工程による前記複数のサブバンドのうち
少なくとも２つ以上のサブバンドに対して閾値を設定す
る閾値設定工程のプログラムと、 前記閾値設定工程により閾値を設定したサブバンド内の
変換係数に対して当該閾値を用いて閾値処理を施す閾値
処理工程のプログラムとを備え、 前記閾値処理工程による閾値処理結果を用いて変更対象
のサブバンド内の変換係数を決定することを特徴とする
請求項１２に記載のプログラム。 - 【請求項１４】 前記閾値設定工程では所定のサブバン
ド内の変換係数を閾値として設定し、 前記閾値処理工程では前記閾値を設定したサブバンド内
の変換係数に対して、当該変換係数に空間的に対応する
前記所定のサブバンド内の変換係数を用いて閾値処理を
施すことを特徴とする請求項１３に記載のプログラム。 - 【請求項１５】 請求項１２乃至１４のいずれか１項に
記載のプログラムを格納し、コンピュータが読み取り可
能な記憶媒体。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001328064A JP4046969B2 (ja) | 2000-11-09 | 2001-10-25 | 画像処理装置、及びその方法、並びにプログラム、記憶媒体 |
| US09/986,252 US6823090B2 (en) | 2000-11-09 | 2001-11-08 | Image processing apparatus and its method, program, and storage medium |
| US10/941,879 US7245776B2 (en) | 2000-11-09 | 2004-09-16 | Image processing apparatus and its method, program and storage medium |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000342341 | 2000-11-09 | ||
| JP2000-342341 | 2000-11-09 | ||
| JP2001328064A JP4046969B2 (ja) | 2000-11-09 | 2001-10-25 | 画像処理装置、及びその方法、並びにプログラム、記憶媒体 |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2007260362A Division JP4630893B2 (ja) | 2000-11-09 | 2007-10-03 | 画像処理装置、及びその方法、並びにプログラム、記憶媒体 |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2002216126A true JP2002216126A (ja) | 2002-08-02 |
| JP2002216126A5 JP2002216126A5 (ja) | 2005-04-21 |
| JP4046969B2 JP4046969B2 (ja) | 2008-02-13 |
Family
ID=26603682
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001328064A Expired - Fee Related JP4046969B2 (ja) | 2000-11-09 | 2001-10-25 | 画像処理装置、及びその方法、並びにプログラム、記憶媒体 |
Country Status (2)
| Country | Link |
|---|---|
| US (2) | US6823090B2 (ja) |
| JP (1) | JP4046969B2 (ja) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004127064A (ja) * | 2002-10-04 | 2004-04-22 | Konica Minolta Holdings Inc | 画像処理方法、画像処理装置、画像処理プログラム及び画像記録装置 |
| WO2005032122A1 (en) * | 2003-09-29 | 2005-04-07 | Samsung Electronics Co., Ltd. | Denoising method and apparatus |
| JP2007094555A (ja) * | 2005-09-27 | 2007-04-12 | Canon Inc | 画像処理装置、画像処理方法 |
| JP2008245033A (ja) * | 2007-03-28 | 2008-10-09 | Institute Of National Colleges Of Technology Japan | 被写界深度制御方法、画像処理装置及びプログラムと記録媒体 |
| JP2008541672A (ja) * | 2005-05-20 | 2008-11-20 | 株式会社エヌ・ティ・ティ・ドコモ | ビデオ符号化におけるノイズフィルタリング方法及び装置 |
| JP2010068508A (ja) * | 2008-08-12 | 2010-03-25 | Canon Inc | 画像処理装置、撮像装置、画像処理方法、画像処理プログラム |
| JP2011255032A (ja) * | 2010-06-10 | 2011-12-22 | Shimadzu Corp | 画像処理方法およびそれを用いた放射線撮影装置 |
| US8478058B2 (en) | 2009-03-09 | 2013-07-02 | Samsung Electronics Co., Ltd. | Apparatus and method for reducing noise from an image |
| JP2018173857A (ja) * | 2017-03-31 | 2018-11-08 | 日本放送協会 | 帯域制限装置及びプログラム |
| JP2019205137A (ja) * | 2018-05-25 | 2019-11-28 | 日本放送協会 | 帯域制限装置及びプログラム |
Families Citing this family (30)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4046969B2 (ja) * | 2000-11-09 | 2008-02-13 | キヤノン株式会社 | 画像処理装置、及びその方法、並びにプログラム、記憶媒体 |
| WO2004004361A1 (ja) * | 2002-06-28 | 2004-01-08 | Nokia Corporation | 情報端末 |
| JP4045422B2 (ja) * | 2002-07-23 | 2008-02-13 | 株式会社日立メディコ | 画像処理装置 |
| US7548656B2 (en) * | 2003-01-16 | 2009-06-16 | Konica Minolta Holdings, Inc. | Method and apparatus for processing image signals by applying a multi-resolution conversion processing for reducing the image size and applying a dyadic wavelet transform |
| US7260272B2 (en) * | 2003-07-10 | 2007-08-21 | Samsung Electronics Co.. Ltd. | Method and apparatus for noise reduction using discrete wavelet transform |
| JP4072491B2 (ja) * | 2003-10-23 | 2008-04-09 | キヤノン株式会社 | 画像処理装置、画像処理方法、プログラム及びコンピュータ可読媒体 |
| JP4006431B2 (ja) * | 2004-10-28 | 2007-11-14 | キヤノン株式会社 | 画像処理方法及び画像処理装置 |
| CN101090668B (zh) * | 2005-01-11 | 2010-05-19 | 株式会社日立医药 | X射线摄影装置 |
| DE102006005803A1 (de) * | 2006-02-08 | 2007-08-09 | Siemens Ag | Verfahren zur Rauschreduktion in bildgebenden Verfahren |
| JP2010530087A (ja) * | 2006-12-19 | 2010-09-02 | 創太 清水 | 画像処理プロセッサ |
| WO2008077132A1 (en) * | 2006-12-19 | 2008-06-26 | California Institute Of Technology | Imaging model and apparatus |
| DE102007013570A1 (de) * | 2007-03-21 | 2008-09-25 | Siemens Ag | Verfahren zur Rauschverminderung in digitalen Bildern mit lokal unterschiedlichem und gerichtetem Rauschen |
| US8280181B2 (en) * | 2007-12-13 | 2012-10-02 | Seiko Epson Corporation | Efficient filtering of image data |
| US8339474B2 (en) * | 2008-08-20 | 2012-12-25 | Freescale Semiconductor, Inc. | Gain controlled threshold in denoising filter for image signal processing |
| US8004588B2 (en) * | 2008-08-27 | 2011-08-23 | Seiko Epson Corporation | Image processing procedure for receiving mosaic image data and calculating vertical and horizontal-direction color difference components for each pixel |
| US8625006B2 (en) * | 2010-07-26 | 2014-01-07 | Sony Corporation | Active imaging device and method for speckle noise reduction |
| CN101944230B (zh) * | 2010-08-31 | 2012-03-28 | 西安电子科技大学 | 基于多尺度的自然图像非局部均值去噪方法 |
| CN102196155B (zh) * | 2011-04-01 | 2013-03-20 | 西安电子科技大学 | 基于Surfacelet变换的系数自适应收缩视频去噪方法 |
| CN102158637B (zh) * | 2011-04-01 | 2012-12-05 | 西安电子科技大学 | 基于Surfacelet变换域的空间自适应阈值视频去噪方法 |
| KR20130064483A (ko) * | 2011-12-08 | 2013-06-18 | 삼성전자주식회사 | 의료 영상 처리 방법 및 장치 |
| CN102547073B (zh) * | 2012-01-02 | 2014-04-02 | 西安电子科技大学 | 基于表面波变换的自适应阈值视频去噪方法 |
| JP5981768B2 (ja) | 2012-05-15 | 2016-08-31 | キヤノン株式会社 | 色処理装置およびその方法 |
| US10001993B2 (en) | 2013-08-08 | 2018-06-19 | Linear Algebra Technologies Limited | Variable-length instruction buffer management |
| US9727113B2 (en) | 2013-08-08 | 2017-08-08 | Linear Algebra Technologies Limited | Low power computational imaging |
| US9910675B2 (en) | 2013-08-08 | 2018-03-06 | Linear Algebra Technologies Limited | Apparatus, systems, and methods for low power computational imaging |
| US11768689B2 (en) | 2013-08-08 | 2023-09-26 | Movidius Limited | Apparatus, systems, and methods for low power computational imaging |
| CN103557856A (zh) * | 2013-10-25 | 2014-02-05 | 哈尔滨工程大学 | 一种光纤陀螺随机漂移实时滤波方法 |
| US9196017B2 (en) * | 2013-11-15 | 2015-11-24 | Linear Algebra Technologies Limited | Apparatus, systems, and methods for removing noise from an image |
| US10460704B2 (en) | 2016-04-01 | 2019-10-29 | Movidius Limited | Systems and methods for head-mounted display adapted to human visual mechanism |
| US10949947B2 (en) | 2017-12-29 | 2021-03-16 | Intel Corporation | Foveated image rendering for head-mounted display devices |
Family Cites Families (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5870602A (en) * | 1987-11-03 | 1999-02-09 | Compaq Computer Corporation | Multi-processor system with system wide reset and partial system reset capabilities |
| US6400996B1 (en) * | 1999-02-01 | 2002-06-04 | Steven M. Hoffberg | Adaptive pattern recognition based control system and method |
| US6850252B1 (en) * | 1999-10-05 | 2005-02-01 | Steven M. Hoffberg | Intelligent electronic appliance system and method |
| JPH06222414A (ja) * | 1993-01-22 | 1994-08-12 | Canon Inc | 防振システム |
| US5742355A (en) * | 1993-12-02 | 1998-04-21 | U.S. Philips Corporation | Method and apparatus for reducing noise in a video signal |
| US5619998A (en) * | 1994-09-23 | 1997-04-15 | General Electric Company | Enhanced method for reducing ultrasound speckle noise using wavelet transform |
| US5953439A (en) * | 1994-11-04 | 1999-09-14 | Ishihara; Ken | Apparatus for and method of extracting time series image information |
| TW351898B (en) * | 1995-12-05 | 1999-02-01 | Advantest Corp | Portrait processing method |
| EP0811289B1 (en) * | 1995-12-21 | 2001-10-17 | Koninklijke Philips Electronics N.V. | Noise reduction in an image |
| US5870502A (en) * | 1996-04-08 | 1999-02-09 | The Trustees Of Columbia University In The City Of New York | System and method for a multiresolution transform of digital image information |
| SE508972C2 (sv) * | 1996-08-28 | 1998-11-23 | Ralip International Ab | Förfarande för kvalitetssäkring vid scanning/kopiering av bilder/dokument samt anordning för genomförande av förfarandet |
| JP3042459B2 (ja) * | 1997-08-25 | 2000-05-15 | 日本電気株式会社 | 映像表示装置 |
| US6771833B1 (en) * | 1999-08-20 | 2004-08-03 | Eastman Kodak Company | Method and system for enhancing digital images |
| US6741739B1 (en) * | 2000-05-03 | 2004-05-25 | Itt Manufacturing Enterprises | Method and apparatus for improving signal to noise ratio using wavelet decomposition and frequency thresholding |
| US6498863B1 (en) * | 2000-09-20 | 2002-12-24 | Media Cybernetics Inc. | Method, system, and product for analyzing a digitized image of an array to create an image of a grid overlay |
| JP4046969B2 (ja) * | 2000-11-09 | 2008-02-13 | キヤノン株式会社 | 画像処理装置、及びその方法、並びにプログラム、記憶媒体 |
-
2001
- 2001-10-25 JP JP2001328064A patent/JP4046969B2/ja not_active Expired - Fee Related
- 2001-11-08 US US09/986,252 patent/US6823090B2/en not_active Expired - Lifetime
-
2004
- 2004-09-16 US US10/941,879 patent/US7245776B2/en not_active Expired - Fee Related
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004127064A (ja) * | 2002-10-04 | 2004-04-22 | Konica Minolta Holdings Inc | 画像処理方法、画像処理装置、画像処理プログラム及び画像記録装置 |
| KR100965880B1 (ko) | 2003-09-29 | 2010-06-24 | 삼성전자주식회사 | 영상의 잡음 제거 방법 |
| WO2005032122A1 (en) * | 2003-09-29 | 2005-04-07 | Samsung Electronics Co., Ltd. | Denoising method and apparatus |
| JP2008541672A (ja) * | 2005-05-20 | 2008-11-20 | 株式会社エヌ・ティ・ティ・ドコモ | ビデオ符号化におけるノイズフィルタリング方法及び装置 |
| JP2007094555A (ja) * | 2005-09-27 | 2007-04-12 | Canon Inc | 画像処理装置、画像処理方法 |
| JP4671417B2 (ja) * | 2005-09-27 | 2011-04-20 | キヤノン株式会社 | 画像処理装置、画像処理方法 |
| JP2008245033A (ja) * | 2007-03-28 | 2008-10-09 | Institute Of National Colleges Of Technology Japan | 被写界深度制御方法、画像処理装置及びプログラムと記録媒体 |
| JP2010068508A (ja) * | 2008-08-12 | 2010-03-25 | Canon Inc | 画像処理装置、撮像装置、画像処理方法、画像処理プログラム |
| US8478058B2 (en) | 2009-03-09 | 2013-07-02 | Samsung Electronics Co., Ltd. | Apparatus and method for reducing noise from an image |
| JP2011255032A (ja) * | 2010-06-10 | 2011-12-22 | Shimadzu Corp | 画像処理方法およびそれを用いた放射線撮影装置 |
| JP2018173857A (ja) * | 2017-03-31 | 2018-11-08 | 日本放送協会 | 帯域制限装置及びプログラム |
| JP2019205137A (ja) * | 2018-05-25 | 2019-11-28 | 日本放送協会 | 帯域制限装置及びプログラム |
| JP7071216B2 (ja) | 2018-05-25 | 2022-05-18 | 日本放送協会 | 帯域制限装置及びプログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| US7245776B2 (en) | 2007-07-17 |
| US20020054713A1 (en) | 2002-05-09 |
| US6823090B2 (en) | 2004-11-23 |
| US20050036707A1 (en) | 2005-02-17 |
| JP4046969B2 (ja) | 2008-02-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4046969B2 (ja) | 画像処理装置、及びその方法、並びにプログラム、記憶媒体 | |
| JP3833177B2 (ja) | 画像処理装置、画象処理方法、記憶媒体及びプログラム | |
| JP3754933B2 (ja) | 画像処理装置、画像処理システム、画像処理方法、プログラム及び記憶媒体 | |
| US6975753B2 (en) | Image processing apparatus, image processing method, program for implementing said method, and storage medium therefor | |
| US7627169B2 (en) | Image processing apparatus and method | |
| US20060110061A1 (en) | Image processing apparatus, image processing method, storage medium, and program | |
| US8139891B2 (en) | System and method for structure enhancement and noise reduction in medical images | |
| US7630575B2 (en) | Image processing apparatus, image processing method, program, and computer-readable medium | |
| JP2000306089A (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP2000224421A (ja) | 画像処理方法および装置並びに記録媒体 | |
| JP4630893B2 (ja) | 画像処理装置、及びその方法、並びにプログラム、記憶媒体 | |
| JP2003337942A (ja) | 画像処理装置 | |
| JP2002092588A (ja) | 画像処理装置及び画像処理方法並びに記憶媒体 | |
| JP4428868B2 (ja) | 画像処理装置及びその方法並びに記憶媒体 | |
| JP3870173B2 (ja) | 画像処理方法、画像処理装置、プログラム及びコンピュータ記録媒体 | |
| JP4235622B2 (ja) | 画像処理装置、画像処理システム、画像処理方法、プログラム及び記憶媒体 | |
| JP4579507B2 (ja) | 画像処理装置 | |
| JP4810002B2 (ja) | 画像処理装置、画像処理システム、画像処理方法、記憶媒体、及びプログラム | |
| JP2004086456A (ja) | 画像処理装置 | |
| JP2005005846A (ja) | 画像処理装置 | |
| JP2004173078A (ja) | 画像処理装置 | |
| JP2003331286A (ja) | 画像処理装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20040611 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20040611 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070705 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070713 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070726 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070910 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20071003 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20071119 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20071121 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101130 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 4046969 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101130 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111130 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121130 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131130 Year of fee payment: 6 |
|
| LAPS | Cancellation because of no payment of annual fees |