JP2010523047A - 信号処理システムにおけるノイズを解析する方法及び装置 - Google Patents
信号処理システムにおけるノイズを解析する方法及び装置 Download PDFInfo
- Publication number
- JP2010523047A JP2010523047A JP2010500918A JP2010500918A JP2010523047A JP 2010523047 A JP2010523047 A JP 2010523047A JP 2010500918 A JP2010500918 A JP 2010500918A JP 2010500918 A JP2010500918 A JP 2010500918A JP 2010523047 A JP2010523047 A JP 2010523047A
- Authority
- JP
- Japan
- Prior art keywords
- noise
- subbands
- channel
- white
- subband
- 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 106
- 238000012545 processing Methods 0.000 title claims abstract description 41
- 230000001902 propagating effect Effects 0.000 claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims description 48
- 238000001914 filtration Methods 0.000 claims description 32
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000002059 diagnostic imaging Methods 0.000 claims description 7
- 238000007405 data analysis Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 abstract description 25
- 238000003384 imaging method Methods 0.000 abstract description 12
- 230000000644 propagated effect Effects 0.000 abstract description 10
- 230000001419 dependent effect Effects 0.000 abstract description 5
- 238000007796 conventional method Methods 0.000 abstract description 2
- 238000004458 analytical method Methods 0.000 description 19
- 238000004088 simulation Methods 0.000 description 10
- 238000007792 addition Methods 0.000 description 9
- 238000013459 approach Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 238000000844 transformation Methods 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000013461 design Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000007246 mechanism Effects 0.000 description 5
- 239000006185 dispersion Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000005311 autocorrelation function Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000005192 partition Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 101100289061 Drosophila melanogaster lili gene Proteins 0.000 description 1
- 101000860173 Myxococcus xanthus C-factor Proteins 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- 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
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/30—Noise filtering
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N9/00—Details of colour television systems
- H04N9/64—Circuits for processing colour signals
- H04N9/646—Circuits for processing colour signals for image enhancement, e.g. vertical detail restoration, cross-colour elimination, contour correction, chrominance trapping filters
-
- 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/10024—Color 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/20016—Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
-
- 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/30168—Image quality inspection
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
- G10L25/03—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
- G10L25/18—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters the extracted parameters being spectral information of each sub-band
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N23/00—Cameras or camera modules comprising electronic image sensors; Control thereof
- H04N23/80—Camera processing pipelines; Components thereof
- H04N23/81—Camera processing pipelines; Components thereof for suppressing or minimising disturbance in the image signal generation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Signal Processing (AREA)
- Image Processing (AREA)
- Noise Elimination (AREA)
- Picture Signal Circuits (AREA)
- Processing Of Color Television Signals (AREA)
Abstract
信号処理システムにおけるノイズを解析する技術を提供する。本明細書では、信号処理システムにおけるノイズを解析するための高速で正確な多重チャンネル周波数依存手法を説明する。ノイズは、各チャンネル内で複数の周波数バンドに分解され、サブバンドノイズが伝播される。畳み込みの計算複雑性を回避するために、従来の方法は、信号処理パイプラインにおけるあらゆる点でノイズを白色と仮定するか、又は単純に空間演算を無視するかのいずれかを行う。各周波数バンド内でノイズが白色であると仮定することにより、あらゆる種類のノイズ(白色、有色、ガウス分布、非ガウス分布、及びその他)は、空間変換を通じて非常に高速で正確な方式で伝播させることができる。この技術の有効性を明らかにするために、画像化パイプラインにおける様々な空間演算を通じてノイズ伝播を考える。更に、計算複雑性は、信号処理システムを通じて画像を伝播させる計算コストのごく僅かな部分でしかない。
【選択図】図1
【選択図】図1
Description
本発明は、信号処理の分野に関する。より具体的には、本発明は、信号処理システムにおけるノイズを解析することに関する。
一連の演算の後の多変量ノイズの解析、伝播、及び予測は、音声信号処理、画像ノイズ除去、デジタル画像システム設計、地震波データ解析、及び医療撮像のような多くの用途で発生する基本的な問題である。これらの演算を構成する信号処理段階は、線形、非線形、及び空間的という3つのカテゴリに分類することができる。既存の方法は、線形及び非線形変換を処理することができる。しかし、これらの方法は、空間信号処理段階を無視して不正確なノイズ伝播をもたらすか、又はこれらの方法は、計算上非常に高価な畳み込みを用いる。従って、これらは、この手法を多重パラメータ最適化、音声におけるノイズ相殺、実時間ノイズ除去におけるサブバンドノイズ予測のような用途に対して非実用的にする。
空間変換の入力におけるノイズが白色であると見なされた場合には、フィルタリングの後のノイズ分散は、入力ノイズエネルギとフィルタエネルギの積として単純に計算することができる。残念ながらこの仮定は、殆どのシステムでは妥当ではない。ノイズは、開始時点では白色である可能性があるが、その特性は、ノイズが信号処理パイプラインを通過する時に変化することになる。従って、白色及び有色ノイズ伝播の両方を処理することができる機構が必要である。
本明細書では、信号処理システムにおけるノイズを解析するための高速で正確な多重チャンネル周波数依存手法を説明する。ノイズは、各チャンネル内で複数の周波数バンドに分解され、サブバンドノイズが伝播される。畳み込みの計算複雑性を回避するために、従来の方法は、信号処理パイプラインにおけるあらゆる点でノイズを白色と仮定するか、又は単純に空間演算を無視するかのいずれかを行う。各周波数バンド内でノイズが白色であると仮定することにより、あらゆる種類のノイズ(白色、有色、ガウス分布、非ガウス分布、及びその他)は、空間変換を通じて非常に高速で正確な方式で伝播させることができる。この技術の有効性を明らかにするために、画像化パイプラインにおける様々な空間演算を通じてノイズ伝播を考える。更に、計算複雑性は、信号処理システムを通じて画像を伝播させる計算コストのごく僅かな部分でしかない。
一態様では、信号内のノイズを解析する方法は、ノイズを複数のサブバンドに分離する段階、複数のサブバンド内のサブバンドの各々に対してノイズ平均値を判断する段階、及びノイズを伝播させる段階を含む。本方法は、ノイズを含む信号を受信する段階を更に含む。複数のサブバンドは、多重チャンネルシステム内の各チャンネルに対するものである。ノイズを伝播させる段階は、実時間で行われる。代替的に、ノイズを伝播させる段階は、オフラインで行われる。本方法は、デジタル画像処理、音処理、地震波データ解析、及び医療撮像から成る群から選択される用途に適用される。好ましくは、サブバンドの各々は、等しい幅のものである。サブバンド数は、望ましい精度に依存する。画像ノイズ除去では、好ましくは、複数のサブバンドは、包含的に2及び4の間の個数である。しかし、あらゆるサブバンド数が可能である。ノイズは、白色ノイズ、有色ノイズ、ガウス及び非ガウス分布ノイズといったあらゆる種類のものとすることができるであろう。
別の態様では、ノイズを解析する方法は、ノイズを含む信号を受信する段階、各チャンネルにおけるノイズを複数のサブバンドに分解する段階、複数のサブバンドの各々に対するノイズ平均値を判断する段階、及びノイズを伝播させる段階を含む。複数のサブバンドは、多重チャンネルシステム内の各チャンネルに対するものである。ノイズを伝播させる段階は、実時間で行われる。代替的に、ノイズを伝播させる段階は、オフラインで行われる。本方法は、デジタル画像処理、音処理、地震波データ解析、及び医療撮像から成る群から選択される用途に適用される。好ましくは、サブバンドの各々は、等しい幅のものである。好ましくは、複数のサブバンドは、包含的に2及び4の間の個数である。ノイズは、白色ノイズ、有色ノイズ、ガウス及び非ガウス分布ノイズといったあらゆる種類のものとすることができるであろう。
別の態様では、空間演算の後のノイズを予測する方法は、ノイズを複数のサブバンドに分解する段階、及びサブバンドノイズを予測する段階を含み、サブバンドノイズは、空間演算の前のサブバンドノイズと、空間演算のサブバンドエネルギと、複数のサブバンドの個数との積である。本方法は、ノイズがバンド毎に白色であると仮定する段階を更に含む。ノイズは、実時間で予測される。好ましくは、サブバンドの各々は、等しい幅のものである。好ましくは、複数のサブバンドは、包含的に2及び4の間の個数である。ノイズは、白色ノイズ、有色ノイズ、ガウス及び非ガウス分布ノイズといったあらゆる種類のものとすることができるであろう。
別の態様では、多重チャンネルシステムから多重チャンネル/多重バンド共分散行列を得る方法は、ノイズをフィルタを用いてチャンネル毎に複数のサブバンドに分割する段階、及び各々がM×Mである部分行列から成るN×N行列を発生させる段階を含む。Nは、サブバンド数であり、Mは、チャンネル数である。フィルタは、理想的及びほぼ理想的な帯域通過フィルタから成る群から選択される。好ましくは、サブバンドの各々は、等しい幅のものである。好ましくは、複数のサブバンドは、包含的に2及び4の間の個数である。ノイズは、白色ノイズ、有色ノイズ、ガウス及び非ガウス分布ノイズといったあらゆる種類のものとすることができるであろう。
更に別の態様では、複数のサブバンドのノイズを伝播させるためのシステムは、ノイズを複数のサブバンドに分割するための分割構成要素、及び複数のサブバンドのノイズを受信して伝播させるための1つ又はそれよりも多くの処理構成要素を含む。分割構成要素は、フィルタである。フィルタは、理想的及びほぼ理想的な帯域通過フィルタから成る群から選択される。代替的に、分割構成要素は、ルックアップテーブルである。複数のサブバンドは、各々がチャンネル毎のものである。ノイズは、白色ノイズ、有色ノイズ、ガウス及び非ガウス分布ノイズといったあらゆる種類のものとすることができるであろう。システムには、「人間視覚システム(HVS)」モデル及び「人間聴覚システム(HAS)」モデルのうちの少なくとも一方を組み込むことができる。
別の態様では、信号処理システムは、ノイズを含む信号を受信するための受信機、ノイズを複数のサブバンドに分割するためのフィルタリング構成要素、及び複数のサブバンドにおけるサブバンドの各々に対して閾値を判断し、閾値よりも低いノイズをフィルタリングするための処理構成要素を含む。フィルタリング構成要素は、理想的及びほぼ理想的な帯域通過フィルタから成る群から選択される。複数のサブバンドは、各々がチャンネル毎のものである。ノイズは、白色ノイズ、有色ノイズ、ガウス及び非ガウス分布ノイズといったあらゆる種類のものとすることができるであろう。システムには、「人間視覚システム(HVS)」モデル及び「人間聴覚システム(HAS)」モデルのうちの少なくとも一方を組み込むことができる。
全ての種類の演算(線形、非線形、及び空間的)を通じて有色及び白色の両方のノイズを伝播させるためのノイズ解析フレームワークを本明細書で説明する。入力ノイズは、サブバンドに分解され、サブバンドノイズは、信号処理パイプラインを通じて伝播される。これを用いると、線形又は非線形(ヤコビアンによって線形化される)変換の後のサブバンドノイズを正確に予測することができる。空間演算の後のノイズを予測するために、空間演算のエネルギは、事前にサブバンドに分解される。ノイズをバンド毎に白色であると見なし、理想的(又はほぼ理想的)帯域通過フィルタリングを用いると、空間演算又は変換の後のサブバンドノイズは、単純に、この演算の前のサブバンドノイズエネルギと、空間変換のサブバンドエネルギと、サブバンド数との積である。
各チャンネル内のノイズをサブバンドに分解し、サブバンドノイズが各バンド内で白色であると仮定することにより、全ての信号処理演算(線形、非線形、及び空間的)を通じてノイズを正確に伝播させることができる。
サブバンドノイズを伝播させることにより、「人間視覚システム(HVS)」モデル及び「人間聴覚システム(HAS)」モデルを画像及び音声システムの設計に組み込むことができる。
本明細書に説明する手法は、計算上非常に効率的であり、多くの計算を必要としない。従って、この手法は、多重パラメータ最適化及び実時間用途において用いることができる。また本明細書に説明する手法は、オフライン用途において用いることができる。
サブバンドノイズを伝播させることにより、「人間視覚システム(HVS)」モデル及び「人間聴覚システム(HAS)」モデルを画像及び音声システムの設計に組み込むことができる。
本明細書に説明する手法は、計算上非常に効率的であり、多くの計算を必要としない。従って、この手法は、多重パラメータ最適化及び実時間用途において用いることができる。また本明細書に説明する手法は、オフライン用途において用いることができる。
下記には、ノイズ解析システムの異なる構成要素及び方法を説明する節が続く。「帯域通過フィルタリング」の節では、入力ノイズ及び信号処理段階の両方をサブバンドに分割する処理を説明する。線形演算にわたって多重チャンネル/多重バンドノイズの伝播させる処理は、「線形変換にわたるノイズ伝播」の節で説明する。「分割フィルタ概念」の節では、分割フィルタ概念を呈示し、理想的(又はほぼ理想的)帯域通過フィルタリングを用いることにより、空間演算の後のサブバンドノイズを非常に高速な方法で予測することができることを示す。「空間変換にわたるノイズ伝播」の節では、ノイズがバンド毎に白色であるという仮定を有する分割フィルタ概念により、モザイク除去、FIRフィルタリング、及びエッジ強調のような空間演算の後のノイズの効率的な予測が可能になる。更に、続く節は、本明細書に説明するノイズ解析のシステム及び方法の計算複雑性及び可能な用途を含む。
帯域通過フィルタリング
ノイズ伝播計算を実行する前に、入力ノイズ及び信号処理段階の両方をサブバンドに分割する必要がある。多重チャンネル入力ノイズ共分散行列Σは、多重バンド/多重チャンネル共分散行列に分割する必要がある。例えば、RGB(3チャンネル)画像処理システムでは、全バンド多重チャンネルノイズ共分散行列は、次式のように書くことができる。
ノイズ伝播計算を実行する前に、入力ノイズ及び信号処理段階の両方をサブバンドに分割する必要がある。多重チャンネル入力ノイズ共分散行列Σは、多重バンド/多重チャンネル共分散行列に分割する必要がある。例えば、RGB(3チャンネル)画像処理システムでは、全バンド多重チャンネルノイズ共分散行列は、次式のように書くことができる。
表記を簡易化するために、ノイズをサブバンドに分割した後に得られる多重バンド/多重チャンネル共分散行列を対称ブロック共分散行列に以下のように書き込む。
ここで、Σij=Σjiである。このサブバンド共分散行列(Σband)はN×N個の部分行列を含む。Nがサブバンド数であり、Mがチャンネル数である時に、各部分行列は、M×Mである。対角部分行列は、バンド内共分散(特定の周波数バンド内での)を含み、それに対して非対角部分行列は、バンド間(周波数バンドの間の)共分散を含む。共分散Σを有する多重チャンネルノイズを3チャンネルシステムにおける多重チャンネル/多重バンド共分散行列Σbandへと変換する処理を図1に例示している。
分割フィルタは、あらゆる種類のものとすることができる。しかし、理想的(又はほぼ理想的)帯域通過フィルタ(BPF)が用いられる場合には、サブバンド間の共分散はゼロになる。従って、
は、ブロック対角になり、ノイズ伝播計算は有意に単純化される。更に、理想的(又はほぼ理想的)BPFでは、「有限インパルス応答(FIR)」フィルタリング、モザイク除去、及びエッジ強調のような空間演算を通じてノイズを伝播させることができる。
次に、線形及び空間演算に対する多重チャンネル/多重バンドノイズ伝播を説明する。非線形演算は、ヤコビアン行列によって線形化することができ、従って、解析は、線形式のものと同様である。
次に、線形及び空間演算に対する多重チャンネル/多重バンドノイズ伝播を説明する。非線形演算は、ヤコビアン行列によって線形化することができ、従って、解析は、線形式のものと同様である。
線形変換にわたるノイズ伝播
この小節では、線形演算の多重バンド/多重チャンネルシステムに対する影響が解析される。図2に示しているように、線形行列が、共分散行列の各々に適用される。線形変換Aの後のブロック共分散行列は、次式のように書くことができる。
この小節では、線形演算の多重バンド/多重チャンネルシステムに対する影響が解析される。図2に示しているように、線形行列が、共分散行列の各々に適用される。線形変換Aの後のブロック共分散行列は、次式のように書くことができる。
コンパクトな形式では、この式は次式へと書き直される。
ここで次式が成り立つ。
変換Aは、周波数依存性を持たないので、Aband内の部分行列は、サブバンドからサブバンドへと変化しない。行列ΣAは、ΣA bandから再構成することができ、ここで、ΣA=AΣATは、線形変換Aによって全バンド入力共分散Σを変換することによって得られる全バンド共分散行列である。
分割フィルタ概念
ここまで、多重チャンネル/多重バンドノイズ伝播機構が、線形演算の後のノイズを正確に予測することができることを示した。上述のように、非線形演算(画像化パイプラインにおけるガンマ等)は、ヤコビアンによって線形化される。従って、その解析は、線形演算のものと同様である。入力サブバンドノイズエネルギ及びフィルタサブバンドエネルギによる空間変換の後の出力サブバンドノイズエネルギを計算するための理論表現をここで導出する。
ここまで、多重チャンネル/多重バンドノイズ伝播機構が、線形演算の後のノイズを正確に予測することができることを示した。上述のように、非線形演算(画像化パイプラインにおけるガンマ等)は、ヤコビアンによって線形化される。従って、その解析は、線形演算のものと同様である。入力サブバンドノイズエネルギ及びフィルタサブバンドエネルギによる空間変換の後の出力サブバンドノイズエネルギを計算するための理論表現をここで導出する。
ノイズ及びフィルタは、図3に示しているように2つのバンドに分割される。第1のサブバンドの入力ノイズ自己相関関数は、次式のように書くことができる。
ゼロ遅延における自己相関関数は、エネルギを表し、このエネルギは、周波数バンドの下にある面積でもある。このエネルギは、サブバンド分散である。従って、最初のサブバンド内のノイズ分散は、σx1x1=2Afcであり、フィルタサブバンドエネルギは、σh1h1=2Bfcであり、出力サブバンドエネルギはσy1y1=2ABfcである。全バンドがN個のサブバンドに分割される場合には、各バンドの幅は、2fc=1/Nである。従って、
である。帯域通過フィルタリングが理想的(又はほぼ理想的、例えば、非常に長いタップフィルタ)である限り、出力サブバンドノイズエネルギは、入力サブバンドノイズエネルギと、空間変換サブバンドエネルギと、サブバンド数との積として単純に推定することができる。一般性を損なうことなくこれを次式のように書くことができる。
ここでi=1、…、Nである。推定される全バンドノイズ分散は、単純にサブバンドノイズ分散の総和である。
この結果を用いると、空間変換の後のノイズを非常に高速で正確な方式で予測することができる。
この手法の有効性を評価するために、厳密な白色ノイズ伝播機構及びバンド毎の白色ノイズ伝播機構が、FIR低域通過フィルタ(単一チャンネル空間変換)を通じて有色ノイズを伝播させると考える。厳密な(畳み込み)白色伝播技術及びバンド毎の白色伝播技術を図4A〜図4Bに例示しており、図4Aは、厳密な手法であり、図4Bは、バンド毎の白色の方法である。図4Aでは、σxxは、入力シーケンスxの分散であり、hは、フィルタであり、σyyは、フィルタリング済みシーケンスyの分散である。図4Bでは、σxixi及びσhihiは、それぞれ入力サブバンド分散及びフィルタサブバンド分散である。推定出力サブバンド分散及び推定出力全バンド分散は、それぞれ、σ^yiyi及びσ^yyであり、ここで、i=1、…、Nであり、Nは、サブバンド数である。有色ノイズは、白色ノイズを0.25の正規化遮断周波数を有する低域通過フィルタによって発生させる。理想的(又はほぼ理想的)帯域通過フィルタリングを用いて、入力ノイズ分散σxx及びFIRフィルタ分散σhhの両方をサブバンド分散σxixi及びσhihiに分解し、ここで、i=1、…、Nである。次に、式(6)及び(7)を用いることにより、サブバンドノイズ分散σ^yiyi及び全バンドノイズ分散σ^yyをフィルタリングの後に推定する。フィルタリングの後の厳密なノイズ分散σyyを得るために、分散σxxの入力ノイズシーケンスxをフィルタhで畳み込む。この畳み込みによって出力シーケンスyが与えられ、このシーケンスからσyyが得られる。最後に、ノイズ伝播精度の基準値として正規化誤差を用いる。
空間変換にわたるノイズ伝播
ノイズがバンド毎に白色であると仮定し、上述の分割フィルタ概念を適用することにより、空間演算の後のノイズを予測することができる。式(5)と同様に、ブロック対角空間変換行列が定められる。
ノイズがバンド毎に白色であると仮定し、上述の分割フィルタ概念を適用することにより、空間演算の後のノイズを予測することができる。式(5)と同様に、ブロック対角空間変換行列が定められる。
ここで、Siは、各周波数バンドにおけるM×M変換行列であり、i=1、…、Nであり、Mは、チャンネル数であり、Nは、サブバンド数である。ブロック対角空間変換行列Sbandがある状態で、空間演算の後のサブバンドノイズ共分散ΣS bandは、次式のようにフィルタリングの前のサブバンドノイズ共分散の項で表すことができる。
空間変換がチャンネル間相関を持たない、例えば、FIRフィルタリングの場合には、行列Siは、対角行列になる。これに反して、この演算がチャンネル間で空間相関を有する、例えば、モザイク除去の場合には、行列Siは、非ゼロの非対角項を有することになる。従って、あらゆる空間演算の後のノイズを予測するための鍵は、変換の周波数領域挙動を正確に反映する行列Siを案出することである。
以下では、フィルタリング、モザイク除去、及びエッジ強調(又は鮮明化)のためのサブバンド変換行列Siを導出する。これらの例は、例示目的のみのものである。同様の方式で、あらゆる空間演算に対して行列Siを形成することができる。
以下では、フィルタリング、モザイク除去、及びエッジ強調(又は鮮明化)のためのサブバンド変換行列Siを導出する。これらの例は、例示目的のみのものである。同様の方式で、あらゆる空間演算に対して行列Siを形成することができる。
FIRフィルタリング
この演算は、チャンネル間にいかなる空間相関も持たないので、行列Siは、M×Mの対角行列である。
この演算は、チャンネル間にいかなる空間相関も持たないので、行列Siは、M×Mの対角行列である。
ここで、σj hihiは、j番目のチャンネルのFIRフィルタをサブバンドに分割することによって得られるi番目のサブバンドの分散である。同じFIRフィルタが全てのチャンネルに用いられる場合には、σ1 hihi=σ2 hihi=…=σM hihi=σ1 hihiであり、式(11)は、更に次式へと単純化することができる。
モザイク除去
単一センサカメラにおける「カラーフィルタアレイ(CFA)」は、一般的にBayerパターンの形態にある。RGBカメラでは、画像処理アレイは、25%の赤色ピクセル、50%の緑色ピクセル、及び25%の青色ピクセルを含む。従って、CFAから得られる部分サンプリングされた多重チャンネルデータ内には、緑色チャンネルは、高低両方の周波数成分を有し、それに対して赤色及び青色のチャンネルは、低周波数情報しか持たない。残りのピクセルを推定するために、研究者は、いくつかのモザイク除去アルゴリズムを提案している。これらの方法は有意な差を有するが、これら全ての手法の固有の性質は、緑色チャンネルからの高周波数情報が、赤色及び青色チャンネルへと複写されることである。モザイク除去のこの周波数領域効果を図6に例示している。
単一センサカメラにおける「カラーフィルタアレイ(CFA)」は、一般的にBayerパターンの形態にある。RGBカメラでは、画像処理アレイは、25%の赤色ピクセル、50%の緑色ピクセル、及び25%の青色ピクセルを含む。従って、CFAから得られる部分サンプリングされた多重チャンネルデータ内には、緑色チャンネルは、高低両方の周波数成分を有し、それに対して赤色及び青色のチャンネルは、低周波数情報しか持たない。残りのピクセルを推定するために、研究者は、いくつかのモザイク除去アルゴリズムを提案している。これらの方法は有意な差を有するが、これら全ての手法の固有の性質は、緑色チャンネルからの高周波数情報が、赤色及び青色チャンネルへと複写されることである。モザイク除去のこの周波数領域効果を図6に例示している。
これらの考察を用いると、式(10)は次式のように書き直される。
上記から、RGBカメラにおけるモザイク除去の本質を捉える行列Siに関する表現は次式になる。
モザイク除去アルゴリズムに依存して、高域通過フィルタの特性は、全く異なるものとすることができ、又は付加的なチャンネル間相関を考察する必要がある場合がある。
エッジ強調又は鮮明化
鮮明化演算は、画像内のエッジ又は他の高周波成分を強調するものである。例えば、高周波成分を強調するために、元の画像xから不鮮明又は平滑なバージョンの画像を減算する手順を図7に例示している。自由パラメータkは、鮮明化の程度を制御する。ノイズをバンド毎に白色であると仮定すると、鮮明化済み画像yのi番目のサブバンドは、次式のように書かれる。
鮮明化演算は、画像内のエッジ又は他の高周波成分を強調するものである。例えば、高周波成分を強調するために、元の画像xから不鮮明又は平滑なバージョンの画像を減算する手順を図7に例示している。自由パラメータkは、鮮明化の程度を制御する。ノイズをバンド毎に白色であると仮定すると、鮮明化済み画像yのi番目のサブバンドは、次式のように書かれる。
従って、鮮明化演算の後のノイズを予測するための行列Siは、次式になる(全てのチャンネルに適用された場合)。
図7の鮮明化の例は、低域通過フィルタと高域通過フィルタリングを達成するための減算器との組合せを用いる。この例は、エッジ強調/鮮明化の一実施形態に過ぎない。同じく本発明の技術によって容易に処理することができるこの演算の他の具現化では、増強する必要がある周波数を抽出する帯域通過フィルタを利用することができる。
計算複雑性
ノイズ伝播は、以下の3つの方式のうちの1つで達成することができ、すなわち、a)「Burns及びBerns(B&B)」方法としても公知の単純な多変量誤差伝播であり。この方法は、ノイズに対する線形及び非線形演算の影響を解析するための用いることができる。しかし、B&B法は空間変換を無視するので、ノイズ解析は、空間演算を含むシステムでは誤差を発生させやすく、次に、b)ノイズを含む1組のカラーパッチが信号処理パイプラインを通じて伝播される画像シミュレーションであり、この方法は厳密であり、次に、c)拡張B&B方法と呼ぶ本明細書で説明するノイズ解析構造によるものであり、拡張B&Bの計算複雑性を一般的なRGB(M=3)デジタルカメラ画像化パイプライン(図8に示している)における画像シミュレーションと比較するものである。
ノイズ伝播は、以下の3つの方式のうちの1つで達成することができ、すなわち、a)「Burns及びBerns(B&B)」方法としても公知の単純な多変量誤差伝播であり。この方法は、ノイズに対する線形及び非線形演算の影響を解析するための用いることができる。しかし、B&B法は空間変換を無視するので、ノイズ解析は、空間演算を含むシステムでは誤差を発生させやすく、次に、b)ノイズを含む1組のカラーパッチが信号処理パイプラインを通じて伝播される画像シミュレーションであり、この方法は厳密であり、次に、c)拡張B&B方法と呼ぶ本明細書で説明するノイズ解析構造によるものであり、拡張B&Bの計算複雑性を一般的なRGB(M=3)デジタルカメラ画像化パイプライン(図8に示している)における画像シミュレーションと比較するものである。
画像シミュレーションは、画像化パイプラインを通じて、所定の分散のノイズを含むチャンネル(この場合はRGB)を有するカラーパッチを伝播させる処理を意味する。それとは対照的に、上述のように、拡張B&B法は、各チャンネルに対するノイズ分散を複数のサブバンドに分割し、画像化パイプラインを通じて多重チャンネル/多重バンド分散を伝播させる。式(7)を用いると、画像化パイプライン内のあらゆる点(YCbCr、L*a*b*等)においてサブバンド分散から全バンド分散を推定することができる。
画像シミュレーション及び拡張B&Bの手法では、解析的表現は、浮動小数点演算の回数(フロップ)で計算される。拡張B&B法の計算負荷は、画像シミュレーションの計算コストのごく僅かな部分でしかない。
画像シミュレーション
利得ブロック及び白色均衡ブロックは、3×3の対角行列を表している。これらのブロックは、画像を伝播させるのにピクセル当たり1回の乗算を必要とする。チャンネル間にいかなる相関も持たない単純な4近傍双線形モザイク除去演算を仮定する。この演算は、ピクセル当たり3回の加算及び1回の乗算を必要とする。この行列及びRGB2YCbCr行列は密であり、ピクセル当たり3回の乗算及び2回の加算を必要とする。ガンマ演算は、ルックアップテーブルによって実行され、浮動小数点計算を必要としない。長さPのFIRフィルタによるフィルタリングは、ピクセル当たりP回の乗算及びP−1回の加算を必要とする。同様に、図7に示している鮮明化演算は、ピクセル当たりP+1回の乗算及びP+3回の加算を必要とする。
利得ブロック及び白色均衡ブロックは、3×3の対角行列を表している。これらのブロックは、画像を伝播させるのにピクセル当たり1回の乗算を必要とする。チャンネル間にいかなる相関も持たない単純な4近傍双線形モザイク除去演算を仮定する。この演算は、ピクセル当たり3回の加算及び1回の乗算を必要とする。この行列及びRGB2YCbCr行列は密であり、ピクセル当たり3回の乗算及び2回の加算を必要とする。ガンマ演算は、ルックアップテーブルによって実行され、浮動小数点計算を必要としない。長さPのFIRフィルタによるフィルタリングは、ピクセル当たりP回の乗算及びP−1回の加算を必要とする。同様に、図7に示している鮮明化演算は、ピクセル当たりP+1回の乗算及びP+3回の加算を必要とする。
伝播される画像のノイズ共分散行列の各エントリを計算するには、平均値を得るのにピクセル当たり1回の加算が必要とされ、分散に対してピクセル当たり2回の加算及び1回の乗算が必要である。この行列は対称であるから、M(M+1)/2個の固有エントリしか存在しない。従って、ノイズ共分散行列を計算するためのピクセル当たりのフロップ数は、2M(M+1)である。平均値及び分散を画像内のピクセル数で評価するには、ノイズ分散−共分散行列内の各エントリにおいて2回(平均値及び分散の各々に対して1回)の乗算が必要である。これらの計算回数は、ピクセル当たりのフロップよりも非常に少ないので、無視することができる。
M個のチャンネルのカメラシステムにおける全演算の加算及び乗算を組み合わせると、ピクセル当たり2M(M+1)+16M+2(2P−1)+2P+4=2M2+18M+6P+2フロップが生じる(Cb及びCrの各チャンネルに対して1つづつ、2つのFIRフィルタが存在する)。長さP=25及びM=3のFIRフィルタの場合には、画像シミュレーションは、ピクセル当たり224フロップを必要とする。従って、K×Lの画像を図8に示している単純な画像化パイプラインを通じて伝播させるには、224KLフロップを必要とすることになる。K、Lは、伝播される画像が、統計的に有意なピクセル数を含むように選択すべきである。K=L=100である場合には、画像シミュレーションの計算オーバーヘッドは224万フロップである。この計算コストは、いかなる実用システムにおいても実行不可能である。
拡張B&B法
全てのフィルタ(図6のモザイク除去高域通過フィルタ、図7の鮮明化低域通過フィルタ、及び図8のFIR低域通過フィルタ)の分散は、画像化パイプラインにおいて事前にN個のサブバンドに分解される。センサ出力におけるN個のサブバンドノイズ分散も同様に既知である。M×M対角行列を通じてノイズを伝播させるためには、サブバンド当たり2M回の乗算が必要である。同様に、密なM×M行列では、各行列要素においてサブバンド当たり3回の乗算及び2回の加算が必要である。これは、密な行列を通じてノイズを伝播させるのにサブバンド当たり5M2フロップが必要であることを意味する。モザイク除去の後のノイズを予測するために、上述の「モザイク除去」節に概説した手順が用いられる。式(13)のSiのまばらな性質に起因して、モザイク除去の後のノイズを計算するのに、サブバンド当たり8Mフロップ(4M回の乗算及び4M回の加算)が含まれる。Cb/Crフィルタリング及びYチャンネルの鮮明化を通じてのノイズ伝播を単純化するために、結合サブバンド変換行列が形成される。
全てのフィルタ(図6のモザイク除去高域通過フィルタ、図7の鮮明化低域通過フィルタ、及び図8のFIR低域通過フィルタ)の分散は、画像化パイプラインにおいて事前にN個のサブバンドに分解される。センサ出力におけるN個のサブバンドノイズ分散も同様に既知である。M×M対角行列を通じてノイズを伝播させるためには、サブバンド当たり2M回の乗算が必要である。同様に、密なM×M行列では、各行列要素においてサブバンド当たり3回の乗算及び2回の加算が必要である。これは、密な行列を通じてノイズを伝播させるのにサブバンド当たり5M2フロップが必要であることを意味する。モザイク除去の後のノイズを予測するために、上述の「モザイク除去」節に概説した手順が用いられる。式(13)のSiのまばらな性質に起因して、モザイク除去の後のノイズを計算するのに、サブバンド当たり8Mフロップ(4M回の乗算及び4M回の加算)が含まれる。Cb/Crフィルタリング及びYチャンネルの鮮明化を通じてのノイズ伝播を単純化するために、結合サブバンド変換行列が形成される。
ここで、
は、鮮明化及びCb、Crフィルタリングにおける低域通過フィルタに対するそれぞれのサブバンドフィルタ分散であり、i=1、…、Nである。この結合変換にわたるノイズ伝播は、サブバンド当たり6Mフロップを必要とする。
従って、拡張B&B手法によってノイズを伝播させるサブバンド当たりのオーバーヘッドは、2対角行列(利得及び白色均衡)に対して4Mフロップ、モザイク除去に対して8Mフロップ、2つの密な行列(行列及びRGBtoYCbCr行列)に対して10M2フロップ、ヤコビアンによって線形化されるガンマ演算に対して6Mフロップ、並びに結合鮮明化及びCb/Crフィルタリング変換に対して6Mフロップである。N個のサブバンドでは、このオーバーヘッドは、合計(10M2+24M)Nフロップである。RGBカメラ(M=3)では、このオーバーヘッドは、162Nフロップになる。2個のサブバンドでは、324フロップが必要であり、これに対して4つのサブバンドは、648フロップを必要とする。従って、拡張B&B法の計算コストは、画像シミュレーションの0.03%よりも少なく、実時間用途に対して間違いなく妥当性がある。
サブバンド分散の事前知識が、信号処理パイプラインの開始時点で得られる用途が存在する。そのような状況では、ノイズ分散をサブバンドに分割する計算コストを含める必要がある。各チャンネルのノイズ分散をN個のサブバンドに分解するには、N+1回の「高速フーリエ変換(FFT)」、つまりノイズシーケンスを周波数領域へと持ち込むのに1回のFFT、更に、分割されたバンドを空間領域へと変換し戻すのにN回の逆FFTが必要である。順FFTの長さがQである場合には、逆FFTの長さはQ/N(N個のバンドに分割される)である。M個のチャンネルに対する計算コストは、
フロップである。64点のFFTでは、ノイズをバンドに分割するコストは、(168N+1160)Mになる。ノイズ分散をサブバンドに分割するための累積的フロップ計数としての計算コストは、4812フロップであり、これに対して4個のサブバンドでは、計算オーバーヘッドは、6144フロップである。明らかに、主な複雑性は、ノイズをサブバンドに分割する手順にある。ルックアップテーブルを用いるなどのより効率的な分割処理は、この機構を一層高速にすることになる。いずれにせよ、バンド分割の計算コストがあったとしても、拡張B&Bは、画像シミュレーションの計算複雑性の0.2%よりも少ない。
用途
デジタル画像システム設計
ノイズの解析、予測、及び抑制は、デジタル画像システム設計においては大きな問題である。取り込み画像は、画像センサ内のノイズ、並びにシステムの電子器械によって発生するノイズによって損なわれる。ノイズは、信号独立及び信号依存の両方の成分を有する。画像は、出力デバイスへと書き込まれる前に、色変換、白色均衡、ガンマ補正、及びフィルタリングのような一連の画像処理段階を受ける。これらの演算は、ノイズ特性を変化させる。Burns及びBerns(B&B)手法としても公知である多変量誤差伝播は、ノイズがデジタル画像化パイプラインを通じて伝播する時のノイズに対する線形及び非線形変換の影響を解析するのに用いられている。上述のように、これらの手法は、空間変換を無視しており、従って、これらのノイズ解析は正確ではない。それとは対照的に、本発明の手法は、全ての変換を処理することができ、非常に正確である。
デジタル画像システム設計
ノイズの解析、予測、及び抑制は、デジタル画像システム設計においては大きな問題である。取り込み画像は、画像センサ内のノイズ、並びにシステムの電子器械によって発生するノイズによって損なわれる。ノイズは、信号独立及び信号依存の両方の成分を有する。画像は、出力デバイスへと書き込まれる前に、色変換、白色均衡、ガンマ補正、及びフィルタリングのような一連の画像処理段階を受ける。これらの演算は、ノイズ特性を変化させる。Burns及びBerns(B&B)手法としても公知である多変量誤差伝播は、ノイズがデジタル画像化パイプラインを通じて伝播する時のノイズに対する線形及び非線形変換の影響を解析するのに用いられている。上述のように、これらの手法は、空間変換を無視しており、従って、これらのノイズ解析は正確ではない。それとは対照的に、本発明の手法は、全ての変換を処理することができ、非常に正確である。
デジタル画像システム設計において、本発明の手法の長所をB&B法と比較するために、有色ノイズを発生させ、一般的なデジタルカメラ画像化パイプラインを代表する線形、非線形、及び空間(FIR低域通過フィルタ)演算を通過させる。有意な比較を得るために、出力を知覚的に均一なL*a*b*色空間へと変換する。この構成のブロック図を図10に提供している。両方の方法(拡張B&B及びB&B)を評価するための基準は、図10の演算の全てを通じて画像を伝播させることによって得られるノイズである(σLL、σaa、及びσbb)。計算複雑性の節で示したように、画像を伝播させることは計算上非常に高負担であり、評価目的のためにのみ行われる。B&B法では、全バンドノイズ分散が伝播し、Cb/CrFIR低域通過フィルタは無視される。一方、拡張B&B技術では、ノイズは、サブバンドに分割され、サブバンドノイズ分散(σ^LiLi、σ^aiai、及びσ^bibi、ここでi=1、…、Nである)が、L*a*b*色空間へと伝播される。サブバンド分散を加算すると(式(7)のように)、L*a*b*における推定全バンドノイズ(σ^LL、σ^aa、及びσ^bb)になる。
図8は、カメラパイプライン800のブロック図を示している。画像センサ802は、ガウス分布、白色、かつ無相関のノイズを送信するが、信号レベル依存性を有する。またこのノイズは、欠損ピクセル(モザイク化された)を有する。利得構成要素804は、低強度光に対して高利得を生成し、高強度光に対して低い利得を生成する。信号ノイズ挙動は、それに従って変化する。白色均衡構成要素806は、照明に依存してR、G、及びBに対する利得を変更する。更に、白色均衡構成要素806の後には、チャンネル依存性が存在する。モザイク除去構成要素808を用いて完全なRGB平面をモザイク除去した後には、周波数依存のチャンネル間相関が存在する。具体的に、Gチャンネルの高周波数ノイズは、B及びRチャンネルへと複写され、低周波数におけるものよりも高いチャンネル間相関が維持される。行列構成要素810の後では、チャンネル間相関は、より複雑である。ガンマ構成要素812の後には、強いレベルの依存性が加算され、ノイズは、もはやガウス分布ではない。RGBからYCbCrへの行列814は、色信号バンド制限に起因して付加的な周波数依存性を加算する。鮮明化/フィルタリング構成要素816は、Y信号の高周波数を増強し、C信号は、バンド制限を受ける。通常、最終的なカメラ出力は、チャンネル毎(YCbCr)、信号レベル毎、及びバンド毎のノイズを有する。このノイズ解析及び方法を用いることにより、得られるノイズを迅速かつ正確に予測することができる。
図8を利用して、B&B及び拡張B&Bの両方の手法に対してa*チャンネルに対する正規化ノイズ伝播誤差εσaaを計算する。サブバンド数のノイズ伝播精度に対する影響を解析するために、周波数バンド数も同様に拡張B&B技術において1個から4個へと変化させる。
図9からは、本発明の手法のノイズ伝播精度が、B&B法のものよりも大幅に優れていることが分る。本発明の手法は、理想的(又はほぼ理想的)帯域通過フィルタリングを通じてノイズ及びフィルタをサブバンドに分解し、ノイズをバンド毎に白色であると仮定することによってCb/Cr低域通過フィルタを効率的に組み込んでいる。それとは対照的に、B&B法は、低域通過FIRフィルタを無視する。単純な設定では、2個のサブバンドで十分である。
図9からは、本発明の手法のノイズ伝播精度が、B&B法のものよりも大幅に優れていることが分る。本発明の手法は、理想的(又はほぼ理想的)帯域通過フィルタリングを通じてノイズ及びフィルタをサブバンドに分解し、ノイズをバンド毎に白色であると仮定することによってCb/Cr低域通過フィルタを効率的に組み込んでいる。それとは対照的に、B&B法は、低域通過FIRフィルタを無視する。単純な設定では、2個のサブバンドで十分である。
会話及び音声処理
ノイズの存在は、会話の明瞭性を阻害する。ノイズ低減アルゴリズムが機能するには、これらのアルゴリズムは、音声信号内に存在するノイズの正確な推定を必要とする。一般的に、干渉ノイズは白色ではない。多重バンド/多重チャンネル手法は、ノイズ除去アルゴリズムが最大潜在力を実時間で発揮することができることを保証する高速方式でサブバンドノイズを正確に予測することができる。
ノイズの存在は、会話の明瞭性を阻害する。ノイズ低減アルゴリズムが機能するには、これらのアルゴリズムは、音声信号内に存在するノイズの正確な推定を必要とする。一般的に、干渉ノイズは白色ではない。多重バンド/多重チャンネル手法は、ノイズ除去アルゴリズムが最大潜在力を実時間で発揮することができることを保証する高速方式でサブバンドノイズを正確に予測することができる。
地震波データ解析
多重チャンネル地震データ解釈は、地球内部を画像化するのに用いられ、地震予測、石油及びガス探査、又は海洋学において用途を有する。そのようなデータの解析における重要な段階はノイズ除去である。ノイズは、コヒーレント、非コヒーレント、及びインパルス性とすることができる。このノイズは、爆風、地面のロール、うねりノイズ、船舶/掘削装置の多重回折ノイズ、表面/回折ノイズ等々のような異なる物理過程に起因する。ノイズの多い観測物から画像を分離するのに、多重チャンネルデータ取得及び処理の方法が用いられる。信号処理演算は、とりわけ、逆畳み込み、フィルタリング、及びスタッキングを含む。様々なサブバンドにおける多重チャンネルノイズを効率的かつ正確に予測して正確なノイズ除去を保証するために本発明の手法を用いることができる。
多重チャンネル地震データ解釈は、地球内部を画像化するのに用いられ、地震予測、石油及びガス探査、又は海洋学において用途を有する。そのようなデータの解析における重要な段階はノイズ除去である。ノイズは、コヒーレント、非コヒーレント、及びインパルス性とすることができる。このノイズは、爆風、地面のロール、うねりノイズ、船舶/掘削装置の多重回折ノイズ、表面/回折ノイズ等々のような異なる物理過程に起因する。ノイズの多い観測物から画像を分離するのに、多重チャンネルデータ取得及び処理の方法が用いられる。信号処理演算は、とりわけ、逆畳み込み、フィルタリング、及びスタッキングを含む。様々なサブバンドにおける多重チャンネルノイズを効率的かつ正確に予測して正確なノイズ除去を保証するために本発明の手法を用いることができる。
医療撮像
医療撮像用途では、画像品質は、極めて重要である。重大な考察事項は、ノイズの解析及び除去である。画像ノイズは、異なる特性を有する可能性がある。例えば、磁気共鳴画像処理におけるノイズは、信号依存及び信号独立性のものとすることができる。一般的に、ノイズは、定常的又は白色ではない。本発明のノイズ解析手法は、ノイズを除去することができるように高速方式でノイズを正確に予測することができる。
図11は、ノイズを判断する処理の流れ図を示している。段階1100では、ノイズを含む信号が受信又は取得される。信号は、デジタルカメラ、カムコーダー、又はいずれか他のコンピュータデバイスのようないずれかのデバイスにおいて受信又は取得される。段階1102では、ノイズは、複数のサブバンドに分離される。段階1104では、バンド毎のノイズ平均値が判断される。
医療撮像用途では、画像品質は、極めて重要である。重大な考察事項は、ノイズの解析及び除去である。画像ノイズは、異なる特性を有する可能性がある。例えば、磁気共鳴画像処理におけるノイズは、信号依存及び信号独立性のものとすることができる。一般的に、ノイズは、定常的又は白色ではない。本発明のノイズ解析手法は、ノイズを除去することができるように高速方式でノイズを正確に予測することができる。
図11は、ノイズを判断する処理の流れ図を示している。段階1100では、ノイズを含む信号が受信又は取得される。信号は、デジタルカメラ、カムコーダー、又はいずれか他のコンピュータデバイスのようないずれかのデバイスにおいて受信又は取得される。段階1102では、ノイズは、複数のサブバンドに分離される。段階1104では、バンド毎のノイズ平均値が判断される。
ノイズ解析のシステム及び方法を利用するために、デバイスは、ノイズをより正確かつ効率的な方式で自動的に解析する。このデバイスは、内部又は外部的のいずれかで本明細書に説明する信号処理法へのアクセスを有する。このノイズ解析法は、ハードウエア、ソフトウエア、又はこれら両方の組合せで実施することができ、ノイズに直面すると、フィルタリングのようなあらゆる望ましい演算がこのノイズに対して実施される。本発明の方法がハードウエア内で実施される時に利用されるハードウエア構成要素は、ノイズを含む信号を受信するための受信機、ノイズをサブバンドに分割するための分割/フィルタリング構成要素、並びにノイズを受信して伝播させるための1つ又はそれよりも多くの処理構成要素を含む。一部の実施形態では、処理構成要素は、複数のサブバンドにおけるサブバンドの各々に対して閾値を判断し、この閾値よりも低いノイズをフィルタリングすることができる。
作動時には、ノイズ解析のシステム及び方法は、ノイズを含む信号を受信又は取得する。ノイズは、各チャンネルにおいて2個又はそれよりも多くのサブバンドに分解される。次に、サブバンドに対する解析を行うことができる。ノイズ解析のシステム及び方法は、白色、有色、ガウス及び非ガウス分布のようなあらゆる種類のノイズの高速で正確な伝播を可能にする。
作動時には、ノイズ解析のシステム及び方法は、ノイズを含む信号を受信又は取得する。ノイズは、各チャンネルにおいて2個又はそれよりも多くのサブバンドに分解される。次に、サブバンドに対する解析を行うことができる。ノイズ解析のシステム及び方法は、白色、有色、ガウス及び非ガウス分布のようなあらゆる種類のノイズの高速で正確な伝播を可能にする。
本発明の構成及び作動の原理の理解を容易にするために、本発明を詳細内容を含む特定的な実施形態に関して説明した。特定的な実施形態及びその詳細内容への本明細書におけるそのような参照は、本明細書に添付する特許請求の範囲を限定するように意図していない。当業者には、特許請求の範囲によって定められる本発明の精神及び範囲から逸脱することなく、例示のために選択された実施形態に他の様々な修正を行うことができることが容易に明らかであろう。
Σ 共分散
ΣNN 多重チャンネル/多重バンド共分散行列
ΣNN 多重チャンネル/多重バンド共分散行列
Claims (41)
- 信号のノイズを解析する方法であって、
a.ノイズを複数のサブバンドに分離する段階、
b.前記複数のサブバンド内のサブバンドの各々に対してノイズ平均値を判断する段階、及び
c.前記ノイズを伝播させる段階、
を含むことを特徴とする方法。 - 前記ノイズを含む前記信号を受信する段階を更に含むことを特徴とする請求項1に記載の方法。
- 前記複数のサブバンドは、多重チャンネルシステム内の各チャンネルに対するものであることを特徴とする請求項1に記載の方法。
- 前記ノイズを伝播させる段階は、実時間で行われることを特徴とする請求項1に記載の方法。
- 前記ノイズを伝播させる段階は、オフラインで行われることを特徴とする請求項1に記載の方法。
- デジタル画像処理、音処理、地震波データ解析、及び医療撮像から成る群から選択される用途に適用されることを特徴とする請求項1に記載の方法。
- 前記サブバンドの各々が、等しい幅のものであることを特徴とする請求項1に記載の方法。
- 前記複数のサブバンドは、包含的に2及び4の間の個数であることを特徴とする請求項1に記載の方法。
- 前記ノイズは、白色ノイズ、有色ノイズ、ガウス分布ノイズ、及び非ガウス分布ノイズから成る群から選択されることを特徴とする請求項1に記載の方法。
- ノイズを解析する方法であって、
a.ノイズを含む信号を受信する段階、
b.各チャンネルの前記ノイズを複数のサブバンドに分解する段階、
c.前記複数のサブバンドの各々に対してノイズ平均値を判断する段階、及び
d.前記ノイズを伝播させる段階、
を含むことを特徴とする方法。 - 前記複数のサブバンドは、多重チャンネルシステム内の各チャンネルに対するものであることを特徴とする請求項10に記載の方法。
- 前記ノイズを伝播させる段階は、実時間で行われることを特徴とする請求項10に記載の方法。
- 前記ノイズを伝播させる段階は、オフラインで行われることを特徴とする請求項10に記載の方法。
- デジタル画像処理、音処理、地震波データ解析、及び医療撮像から成る群から選択される用途に適用されることを特徴とする請求項10に記載の方法。
- 前記サブバンドの各々が、等しい幅のものであることを特徴とする請求項10に記載の方法。
- 前記複数のサブバンドは、包含的に2及び4の間の個数であることを特徴とする請求項10に記載の方法。
- 前記ノイズは、白色ノイズ、有色ノイズ、ガウス分布ノイズ、及び非ガウス分布ノイズから成る群から選択されることを特徴とする請求項10に記載の方法。
- 空間演算の後のノイズを予測する方法であって、
a.ノイズを複数のサブバンドに分解する段階、及び
b.サブバンドノイズを予測する段階、
を含み、
前記サブバンドノイズは、空間演算の前の該サブバンドノイズと、空間演算のサブバンドエネルギと、前記複数のサブバンドの個数との積である、
ことを特徴とする方法。 - 前記ノイズがバンド毎に白色であると仮定する段階を更に含むことを特徴とする請求項18に記載の方法。
- 前記ノイズは、実時間で予測されることを特徴とする請求項18に記載の方法。
- 前記サブバンドの各々が、等しい幅のものであることを特徴とする請求項18に記載の方法。
- 前記複数のサブバンドは、包含的に2及び4の間の個数であることを特徴とする請求項18に記載の方法。
- 前記ノイズは、白色ノイズ、有色ノイズ、ガウス分布ノイズ、及び非ガウス分布ノイズから成る群から選択されることを特徴とする請求項18に記載の方法。
- 多重チャンネルシステムから多重チャンネル/多重バンド共分散行列を取得する方法であって、
a.フィルタを用いてノイズをチャンネル毎に複数のサブバンドに分割する段階、及び
b.各々がM×Mである部分行列から成るN×N行列を発生させる段階、
を含むことを特徴とする方法。 - Nは、サブバンドの数であり、Mは、チャンネルの数であることを特徴とする請求項24に記載の方法。
- 前記フィルタは、理想的及びほぼ理想的な帯域通過フィルタから成る群から選択されることを特徴とする請求項24に記載の方法。
- 前記サブバンドの各々が、等しい幅のものであることを特徴とする請求項24に記載の方法。
- 前記複数のサブバンドは、包含的に2及び4の間の個数であることを特徴とする請求項24に記載の方法。
- 前記ノイズは、白色ノイズ、有色ノイズ、ガウス分布ノイズ、及び非ガウス分布ノイズから成る群から選択されることを特徴とする請求項24に記載の方法。
- 複数のサブバンドのノイズを伝播させるためのシステムであって、
a.ノイズを複数のサブバンドに分割するための分割構成要素、及び
b.前記複数のサブバンドの前記ノイズを受信して伝播させるための1つ又はそれよりも多くの処理構成要素、
を含むことを特徴とするシステム。 - 前記分割構成要素は、フィルタであることを特徴とする請求項30に記載のシステム。
- 前記フィルタは、理想的及びほぼ理想的な帯域通過フィルタから成る群から選択されることを特徴とする請求項31に記載のシステム。
- 前記分割構成要素は、ルックアップテーブルであることを特徴とする請求項30に記載のシステム。
- 前記複数のサブバンドは、各々がチャンネル毎のものであることを特徴とする請求項30に記載のシステム。
- 前記ノイズは、白色ノイズ、有色ノイズ、ガウス分布ノイズ、及び非ガウス分布ノイズから成る群から選択されることを特徴とする請求項30に記載のシステム。
- システムに「人間視覚システム(HVS)」モデル及び「人間聴覚システム(HAS)」モデルのうちの少なくとも一方が組み込まれることを特徴とする請求項30に記載のシステム。
- a.ノイズを含む信号を受信するための受信機、
b.前記ノイズを複数のサブバンドに分割するためのフィルタリング構成要素、及び
c.前記複数のサブバンドにおけるサブバンドの各々に対して閾値を判断し、該閾値よりも低い前記ノイズをフィルタリングするための処理構成要素、
を含むことを特徴とする信号処理システム。 - 前記フィルタリング構成要素は、理想的及びほぼ理想的な帯域通過フィルタから成る群から選択されることを特徴とする請求項37に記載のシステム。
- 前記複数のサブバンドは、各々がチャンネル毎のものであることを特徴とする請求項37に記載のシステム。
- 前記ノイズは、白色ノイズ、有色ノイズ、ガウス分布ノイズ、及び非ガウス分布ノイズから成る群から選択されることを特徴とする請求項37に記載のシステム。
- システムに「人間視覚システム(HVS)」モデル及び「人間聴覚システム(HAS)」モデルのうちの少なくとも一方が組み込まれることを特徴とする請求項37に記載のシステム。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11/731,261 US8108211B2 (en) | 2007-03-29 | 2007-03-29 | Method of and apparatus for analyzing noise in a signal processing system |
PCT/US2008/003468 WO2008121225A1 (en) | 2007-03-29 | 2008-03-13 | Method of and apparatus for analyzing noise in a signal processing system |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2010523047A true JP2010523047A (ja) | 2010-07-08 |
Family
ID=39794262
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010500918A Pending JP2010523047A (ja) | 2007-03-29 | 2008-03-13 | 信号処理システムにおけるノイズを解析する方法及び装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US8108211B2 (ja) |
EP (1) | EP2132887A4 (ja) |
JP (1) | JP2010523047A (ja) |
CN (1) | CN101647215B (ja) |
WO (1) | WO2008121225A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012063394A (ja) * | 2010-09-14 | 2012-03-29 | Casio Comput Co Ltd | 雑音抑制装置、雑音抑制方法およびプログラム |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8436909B2 (en) * | 2008-10-21 | 2013-05-07 | Stmicroelectronics S.R.L. | Compound camera sensor and related method of processing digital images |
JP5267794B2 (ja) * | 2008-12-26 | 2013-08-21 | 株式会社リコー | 画像処理装置及び車載カメラ装置 |
US20110286305A1 (en) * | 2009-01-20 | 2011-11-24 | Spectraseis Ag | Time reverse imaging operators for source location |
US20110075935A1 (en) * | 2009-09-25 | 2011-03-31 | Sony Corporation | Method to measure local image similarity based on the l1 distance measure |
JP5640370B2 (ja) * | 2009-12-18 | 2014-12-17 | ソニー株式会社 | 画像処理装置,画像処理方法及び撮像装置 |
JP5574792B2 (ja) * | 2010-04-13 | 2014-08-20 | キヤノン株式会社 | 撮影装置 |
JP5679791B2 (ja) * | 2010-12-14 | 2015-03-04 | キヤノン株式会社 | 画像投影装置及びその制御方法、プログラム |
US8918197B2 (en) | 2012-06-13 | 2014-12-23 | Avraham Suhami | Audio communication networks |
SE536510C2 (sv) * | 2012-02-21 | 2014-01-14 | Flir Systems Ab | Bildbehandlingsmetod för detaljförstärkning och brusreduktion |
CN103076028A (zh) * | 2013-01-21 | 2013-05-01 | 新疆美特智能安全工程股份有限公司 | 光相位振动的小波消噪方法 |
US20160314800A1 (en) * | 2013-12-23 | 2016-10-27 | Analog Devices, Inc. | Computationally efficient method for filtering noise |
US9667842B2 (en) | 2014-08-30 | 2017-05-30 | Apple Inc. | Multi-band YCbCr locally-adaptive noise modeling and noise reduction based on scene metadata |
US9525804B2 (en) | 2014-08-30 | 2016-12-20 | Apple Inc. | Multi-band YCbCr noise modeling and noise reduction based on scene metadata |
JP6080190B2 (ja) * | 2014-09-15 | 2017-02-15 | エスゼット ディージェイアイ テクノロジー カンパニー リミテッドSz Dji Technology Co.,Ltd | デモザイク化の方法 |
US9641820B2 (en) | 2015-09-04 | 2017-05-02 | Apple Inc. | Advanced multi-band noise reduction |
US11143774B2 (en) * | 2015-12-22 | 2021-10-12 | Shell Oil Company | Method and system for separating blended seismic data |
CN109842379B (zh) * | 2017-11-29 | 2023-09-19 | 北京振兴计量测试研究所 | 一种宽带噪声产生方法 |
DE102018107291A1 (de) * | 2018-03-27 | 2019-10-02 | Voestalpine Automotive Components Linz Gmbh | Verfahren zum Schweißen beschichteter Stahlbleche |
US11481677B2 (en) * | 2018-09-30 | 2022-10-25 | Shearwater Geoservices Software Inc. | Machine learning based signal recovery |
CN109614744B (zh) * | 2018-12-28 | 2022-11-22 | 浙江理工大学 | 一种基于大数据的降水量检测方法及系统 |
CN113253336B (zh) * | 2021-07-02 | 2021-10-01 | 深圳市翩翩科技有限公司 | 一种基于深度学习的地震预测方法和系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10254499A (ja) * | 1997-03-14 | 1998-09-25 | Nippon Telegr & Teleph Corp <Ntt> | 帯域分割型雑音低減方法及び装置 |
JP2004502977A (ja) * | 2000-07-12 | 2004-01-29 | アンドレア エレクトロニクス コーポレイション | サブバンド指数平滑雑音消去システム |
Family Cites Families (43)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3971065A (en) * | 1975-03-05 | 1976-07-20 | Eastman Kodak Company | Color imaging array |
FR2389277A1 (fr) * | 1977-04-29 | 1978-11-24 | Ibm France | Procede de quantification a allocation dynamique du taux de bits disponible, et dispositif de mise en oeuvre dudit procede |
JPS60140399A (ja) * | 1983-12-28 | 1985-07-25 | 松下電器産業株式会社 | 雑音除去装置 |
US5012333A (en) * | 1989-01-05 | 1991-04-30 | Eastman Kodak Company | Interactive dynamic range adjustment system for printing digital images |
US6289308B1 (en) * | 1990-06-01 | 2001-09-11 | U.S. Philips Corporation | Encoded wideband digital transmission signal and record carrier recorded with such a signal |
JPH07249101A (ja) | 1994-03-08 | 1995-09-26 | Hitachi Eng Co Ltd | ニューロ利用の文字認識装置及び輝度正規化方法 |
US6098038A (en) * | 1996-09-27 | 2000-08-01 | Oregon Graduate Institute Of Science & Technology | Method and system for adaptive speech enhancement using frequency specific signal-to-noise ratio estimates |
FI104521B (fi) * | 1997-01-30 | 2000-02-15 | Nokia Multimedia Network Termi | Menetelmä kuvasekvenssien kontrastin parantamiseksi |
EP0878790A1 (en) * | 1997-05-15 | 1998-11-18 | Hewlett-Packard Company | Voice coding system and method |
US6847737B1 (en) * | 1998-03-13 | 2005-01-25 | University Of Houston System | Methods for performing DAF data filtering and padding |
US6226321B1 (en) * | 1998-05-08 | 2001-05-01 | The United States Of America As Represented By The Secretary Of The Air Force | Multichannel parametric adaptive matched filter receiver |
US7209567B1 (en) * | 1998-07-09 | 2007-04-24 | Purdue Research Foundation | Communication system with adaptive noise suppression |
US6643398B2 (en) * | 1998-08-05 | 2003-11-04 | Minolta Co., Ltd. | Image correction device, image correction method and computer program product in memory for image correction |
US6600833B1 (en) * | 1999-07-23 | 2003-07-29 | Intel Corporation | Methodology for color correction with noise regulation |
US6480610B1 (en) * | 1999-09-21 | 2002-11-12 | Sonic Innovations, Inc. | Subband acoustic feedback cancellation in hearing aids |
US6963544B1 (en) * | 1999-12-10 | 2005-11-08 | Lucent Technologies Inc. | System for statistically multiplexing real-time and non-real-time voice and data traffic in a wireless system |
US6813389B1 (en) * | 1999-12-15 | 2004-11-02 | Eastman Kodak Company | Digital image processing method and system including noise reduction and tone scale adjustments |
EP1124384A1 (en) * | 2000-02-10 | 2001-08-16 | Sony Corporation | Information processing apparatus for compressed image data using class-classification adaptive processing |
JP3758452B2 (ja) * | 2000-02-28 | 2006-03-22 | コニカミノルタビジネステクノロジーズ株式会社 | 記録媒体、並びに、画像処理装置および画像処理方法 |
US7173657B2 (en) | 2000-03-24 | 2007-02-06 | Rochester Institute Of Technology | System and method for analyzing, optimizing and/or designing spectral sensitivity curves for an imaging device |
US7062108B2 (en) * | 2001-06-12 | 2006-06-13 | Eastman Kodak Company | Method for estimating the appearance of noise in images |
TW533406B (en) * | 2001-09-28 | 2003-05-21 | Ind Tech Res Inst | Speech noise elimination method |
JP4406195B2 (ja) * | 2001-10-09 | 2010-01-27 | セイコーエプソン株式会社 | 画像データの出力画像調整 |
US7236195B2 (en) * | 2001-11-06 | 2007-06-26 | Rochester Institute Of Technology | Method and system for optimizing a selection of spectral sensitivities |
EP1345172A1 (en) * | 2002-02-26 | 2003-09-17 | Sony International (Europe) GmbH | Contrast enhancement for digital images |
JP2003284084A (ja) | 2002-03-20 | 2003-10-03 | Sony Corp | 画像処理装置および方法、並びに画像処理装置の製造方法 |
US20030190090A1 (en) * | 2002-04-09 | 2003-10-09 | Beeman Edward S. | System and method for digital-image enhancement |
US6987881B2 (en) * | 2002-04-26 | 2006-01-17 | Keio University | Automatic visibility improvement method for digital image |
JP4133029B2 (ja) * | 2002-06-25 | 2008-08-13 | 富士フイルム株式会社 | 画像処理方法および装置 |
JP4224996B2 (ja) * | 2002-07-25 | 2009-02-18 | パナソニック株式会社 | 撮像装置 |
JP3762725B2 (ja) | 2002-08-22 | 2006-04-05 | オリンパス株式会社 | 撮像システムおよび画像処理プログラム |
JP3861808B2 (ja) | 2002-12-25 | 2006-12-27 | ソニー株式会社 | 撮像装置 |
US7636436B2 (en) * | 2003-05-15 | 2009-12-22 | At&T Intellectual Property Ii, L.P. | Method for the performance testing of echo cancellers using an artificial segmented test signal |
KR100510532B1 (ko) * | 2003-06-10 | 2005-08-26 | 삼성전자주식회사 | 베이어 패턴 컬러 신호에 대한 적응형 노이즈 제거 필터,이를 구비한 디지털 영상 신호 처리 장치, 및 그 방법 |
JP2007506348A (ja) | 2003-09-23 | 2007-03-15 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | インバンド動き補償時間フィルタリングを利用したビデオノイズ除去アルゴリズム |
JP2005303803A (ja) | 2004-04-14 | 2005-10-27 | Olympus Corp | 撮像装置と画像記録媒体、および画像処理装置ならびに画像処理プログラムとその記録媒体 |
WO2007002406A2 (en) | 2005-06-20 | 2007-01-04 | The Trustees Of Columbia University In The City Of New York | Interactive diagnostic display system |
US7791656B2 (en) * | 2005-08-16 | 2010-09-07 | Konica Minolta Holdings, Inc. | Image sensing apparatus and image processing method |
KR101128454B1 (ko) * | 2005-11-10 | 2012-03-23 | 삼성전자주식회사 | 콘트라스트 향상 방법 및 장치 |
KR100756318B1 (ko) * | 2006-01-16 | 2007-09-06 | 삼성전자주식회사 | 이득 조절 가능한 절단 빈도함수 평활화를 이용한 영상처리 방법 및 시스템 |
US8711249B2 (en) * | 2007-03-29 | 2014-04-29 | Sony Corporation | Method of and apparatus for image denoising |
KR20090008745A (ko) * | 2007-07-18 | 2009-01-22 | 삼성전자주식회사 | 히스토그램 등화 장치 및 방법 |
US20090161982A1 (en) * | 2007-12-19 | 2009-06-25 | Nokia Corporation | Restoring images |
-
2007
- 2007-03-29 US US11/731,261 patent/US8108211B2/en not_active Expired - Fee Related
-
2008
- 2008-03-13 CN CN2008800104323A patent/CN101647215B/zh not_active Expired - Fee Related
- 2008-03-13 JP JP2010500918A patent/JP2010523047A/ja active Pending
- 2008-03-13 WO PCT/US2008/003468 patent/WO2008121225A1/en active Application Filing
- 2008-03-13 EP EP08726878.5A patent/EP2132887A4/en not_active Ceased
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10254499A (ja) * | 1997-03-14 | 1998-09-25 | Nippon Telegr & Teleph Corp <Ntt> | 帯域分割型雑音低減方法及び装置 |
JP2004502977A (ja) * | 2000-07-12 | 2004-01-29 | アンドレア エレクトロニクス コーポレイション | サブバンド指数平滑雑音消去システム |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012063394A (ja) * | 2010-09-14 | 2012-03-29 | Casio Comput Co Ltd | 雑音抑制装置、雑音抑制方法およびプログラム |
Also Published As
Publication number | Publication date |
---|---|
EP2132887A4 (en) | 2015-04-01 |
US20080240203A1 (en) | 2008-10-02 |
EP2132887A1 (en) | 2009-12-16 |
WO2008121225A1 (en) | 2008-10-09 |
CN101647215B (zh) | 2013-12-25 |
US8108211B2 (en) | 2012-01-31 |
CN101647215A (zh) | 2010-02-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2010523047A (ja) | 信号処理システムにおけるノイズを解析する方法及び装置 | |
Aissa-El-Bey et al. | Underdetermined blind separation of nondisjoint sources in the time-frequency domain | |
Argenti et al. | Multiresolution MAP despeckling of SAR images based on locally adaptive generalized Gaussian pdf modeling | |
US7603401B2 (en) | Method and system for on-line blind source separation | |
CN106233382B (zh) | 一种对若干个输入音频信号进行去混响的信号处理装置 | |
US20110002473A1 (en) | Dereverberation apparatus, dereverberation method, dereverberation program, and recording medium | |
Chen et al. | Enhancing seismic reflections using empirical mode decomposition in the flattened domain | |
Dörfler et al. | Nonstationary Gabor frames-approximately dual frames and reconstruction errors | |
Rao et al. | A denoising approach to multisensor signal estimation | |
Gosse et al. | Perfect reconstruction versus MMSE filter banks in source coding | |
Neo et al. | Polynomial eigenvalue decomposition for multichannel broadband signal processing: a mathematical technique offering new insights and solutions | |
CN117111000A (zh) | 一种基于双通道注意力残差网络的sar梳状谱干扰抑制方法 | |
Belaid et al. | A new multi-scale framework for convolutive blind source separation | |
KR100454886B1 (ko) | 독립 성분 분석을 이용한 여파기 적응 알고리즘의 필터뱅크 접근 방법 | |
RU2282209C1 (ru) | Способ и устройство обнаружения сложных широкополосных частотно-модулированных сигналов с фильтрацией в масштабно-временной области | |
RU2439601C1 (ru) | Устройство обнаружения сложных широкополосных частотно-модулированных сигналов с фильтрацией в масштабно-временной области на основе дискретного вейвлет-преобразования | |
Sakiyama et al. | Edge-aware image graph expansion methods for oversampled graph Laplacian matrix | |
CN110491408B (zh) | 一种基于稀疏元分析的音乐信号欠定混叠盲分离方法 | |
Rashobh et al. | Multichannel equalization in the KLT and frequency domains with application to speech dereverberation | |
CN108074580B (zh) | 一种噪声消除方法及装置 | |
Mirarab et al. | A novel wavelet based ICA technique using Kurtosis | |
US20160148361A1 (en) | Method and apparatus for filtering an array of pixels | |
Tang et al. | Multiscale blind image restoration using a wavelet decomposition | |
Koh et al. | Broadband GSC beamformer with spatial and temporal decorrelation | |
Shu et al. | Improvement of DOA estimation using wavelet denoising |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20101019 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20120913 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20120918 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20130304 |