[論文レビュー] Optimum window length of Savitzky-Golay filters with arbitrary order

Mohammad Sadeghi, Fereidoon Behnia|arXiv (Cornell University)|Aug 30, 2018

Image and Signal Denoising Methods参考文献 17被引用数 23

ひとこと要約

本稿では、バイアスとばらつきのバランスをとることで平均二乗誤差（MSE）を最小化する任意の次数のSavitzky-Golayフィルタの最適なウィンドウ長を決定する最適化アルゴリズムを提案する。チェビシェフ直交多項式を用いてフィルタをモデル化することで、多様な信号タイプにおいて従来手法を上回る精度とロバスト性を実現する。

ABSTRACT

One of the widely used denoising methods in different domains is the Savitzky-Golay (SG) filter. The SG filter has two design parameters: window length and the filter order. As the length of the window increases, the estimation variance decreases, but the bias error increases at the same time. Mean square error (MSE) measure includes both bias and variance criteria. In this paper, we obtain the optimal window length of an SG filter with arbitrary order which minimizes the MSE. To achieve the optimal window length, we propose an algorithm whose performance is better than the existing methods. In this paper, we follow the viewpoint proposed by Persson and Strang and design the filter on the basis of Chebyshev orthogonal polynomials

研究の動機と目的

• 任意の次数のSavitzky-Golayフィルタに対して、平均二乗誤差（MSE）を最小化する最適なウィンドウ長を特定すること。
• SGフィルタ設計における推定ばらつきとバイアス誤差のトレードオフを扱うこと。
• 従来手法よりも精度とロバスト性に優れた、最適なウィンドウ長を決定するアルゴリズムの開発。
• 数値的安定性と性能向上を図るため、Savitzky-Golayフィルタ設計をチェビシェフ直交多項式に基づいて再定式化すること。

提案手法

• 著者らはPerssonとStrangが提唱したフレームワークを採用し、チェビシェフ重み関数に基づく直交多項式を用いてSavitzky-Golayフィルタを再定式化する。
• 信号を対称的なウィンドウ上での直交多項式基底に射影することにより、フィルタ係数の解析的表現を導出する。
• 最適なウィンドウ長は、信号の曲率とノイズパワーの関数として表されるMSEを最小化することで決定される。
• MSEをさまざまなウィンドウ長について計算し、最小の誤差を生じるものを選択するアルゴリズムが、直交多項式展開の性質を活用して実装される。
• 標準的な単項式ベースのアプローチで一般的に見られる悪条件な多項式基底集合を避けることで、数値的安定性を確保する。
• 提案されたアルゴリズムは、合成信号および実世界の信号を用いて評価され、MSEを最小化する点で従来手法を上回ることを検証する。

実験結果

リサーチクエスチョン

• RQ1 任意の次数のSavitzky-Golayフィルタに対して、平均二乗誤差（MSE）を最小化する最適なウィンドウ長は何か？
• RQ2 チェビシェフ直交多項式の使用は、標準的な多項式基底と比較して、SGフィルタ設計の精度と安定性をどのように向上させるか？
• RQ3 提案されたアルゴリズムは、多様な信号特性にわたって、従来手法を上回って最適なウィンドウ長を決定できるか？
• RQ4 ウィンドウ長、フィルタ次数、およびそれによるバイアス-ばらつきトレードオフの関係は何か？

主な発見

• 提案されたアルゴリズムは、さまざまな信号タイプにおいて、従来手法を上回る低い平均二乗誤差（MSE）を達成する。
• チェビシェフ直交多項式を用いることで、標準的な単項式ベースの多項式フィッティングと比較して、数値的安定性が著しく向上し、条件数が低減される。
• 最適なウィンドウ長は、フィルタ次数と信号の曲率の両方に依存しており、高次フィルタでは低バイアスを維持するためにはより長いウィンドウが必要であることが示された。
• MSEをウィンドウ長とフィルタ次数の関数として閉形式で表現でき、正確な最適化が可能になる。
• ノイズレベルや曲率が異なる合成信号において、本手法は従来のヒューリスティックなウィンドウ選択ルールを上回るロバストな性能を示した。
• 結果から、本手法によるMSEの最小化が、固定値または経験的に選択されたウィンドウ長と比較して、優れたノイズ除去性能をもたらすことが確認された。