フィルタ処理の高速化アルゴリズム(縦横に処理を分ける)

シェアする

  • このエントリーをはてなブックマークに追加
最近の記事
  • 7/6 【参考書籍】画像処理・機械学習プログラミング OpenCV 3対応
  • 6/20 【Python,matplotlib】動くグラフをAnimationGifに保存する方法
  • 6/17 シグモイド関数の微分
  • 6/15 シグモイド関数
  • 6/13 合成関数の微分
  • 6/12 WordPressで数式エディタ風に数式を入力したい
  • 6/11 PythonをVisual Studioでインストールする方法
  • 6/9 【Python】OpenCVをAnacondaでインストール(Windows編)
  • 6/6 【Python】Anacondaで複数バージョンの環境切り替え
  • 6/6 画像センシング展2017に出展します。
  • 6/1 【Office365】Web版Outlookのフォントサイズ変更
  • 6/1 【Anaconda】モジュールのアップデートでエラー発生
  • 6/1 【Anaconda】コマンドリストの表示
  • 5/29 Windows10パソコン購入
  • 5/24 Anacondaのアンインストール
  • 5/24 【Jupyter Notebook】新規プログラムの作成
  • 5/23 【Python】開発環境の構築
  • 5/23 Pythonはじめました
  • 4/6 【Office365】Web版Outlookのスレッド表示を解除する方法
  • 4/5 【Excel】フーリエ解析(FFT)
  • 3/20 Canny edge detection
  • 3/20 【Excel2016】分析ツールの表示
  • 3/5 【Visual Studio】黒い背景色を白に変更する方法
  • 2/8 【Windows10】拡張モニタに表示されたウィンドウを元に戻す
  • 2/7 複素数の計算
  • 1/18 【Excel】棒グラフの横軸の目盛を0始まりにする
  • 1/16 【Excel】フーリエ変換
  • 1/6 【OpenCV】疑似カラー(カラーマップ)
  • 11/8 【Visual Studio】検索結果のウィンドウ表示
  • 11/3 ニコン 一眼レフカメラ D5500レビュー
  • 10/26 カラーカメラはモノクロカメラを兼ねない
  • 9/6 (Free Soft)Animation GIF Builder
  • 8/30 【C#】タブの無いTabControlっぽいものを作る
  • 8/29 【OpenCvSharp】サンプルプログラムの公開
  • 8/28 【PowerPoint】部分的にカラーにする(セレクトカラー処理)
  • 8/27 【C#】引数の値渡し、参照渡し(ref, out)
  • 8/26 【Word/Excel】図形内に文字を挿入する
  • 8/25 【C#】NumericUpDownコントロール
  • 8/24 ニコン D3400 VS D5500 仕様比較
  • 8/22 【C#】MDIフォームにリサイズ可能なPanelを設置する
  • 8/20 【Visual Studio 2015】文字の色が変、かすむ
  • 8/20 【Visual Studio】行のコピー、切り取り、貼り付け
  • 8/20
  • 8/16
  • 7/7 標準偏差のよくある誤解
  • 6/17 ホーム
  • 6/15
  • 6/15
  • 6/15
  • 6/15

  • 前回、フィルタ処理の高速化アルゴリズム(重複した計算を行わない)で紹介した方法ではカーネルの値が全て同じでないと使えないので、今回はフィルタ処理を縦方向と横方向に分けて行う事でフィルタ処理の高速化を行う方法をガウシアンフィルタを例にとって紹介します。

    ガウシアンフィルタのカーネルには、

    ガウシアンフィルタのカーネル

    が良く用いられますが、この処理を注目画素の周辺の輝度値をI0~I8とした場合、
    ガウシアンフィルタの処理を行列で

    ガウシアンフィルタ

    と、表すこともでき、この事は縦方向に3×1のガウシアンフィルタ処理をおこなってから、
    横方向に1×3のガウシアンフィルタ処理を行うことを意味しています。
    (横方向に処理をしてから縦方向に処理をしても同じです。)

    このように処理を縦と横に分けることで、カーネルのサイズm×nの場合、通常の処理では
    m×n回の掛け算を行うところ、m+n回の掛け算で済む事になります。
    (ただし、縦横に処理を分ける事で全画素を2回参照することになるので、カーネルのサイズが
    小さいと効果はあまりありません。)

    他にも、移動平均フィルタの場合

    移動平均フィルタ

    ソーベルフィルタの場合

    ソーベルフィルタ

    となります。
    ソーベルフィルタの行列を見ると、縦方向にガウシアンフィルタ処理をしてから、横方向に微分処理している事が分かりやすくなっているかと思います。

    また、比較的処理の重いメディアンフィルタにおいても、処理を縦と横に分けることによって、
    ほぼ、同様な効果を得ることができます。
    厳密には同じ結果にはならないのですが、スパイクノイズを除去するという意味では
    十分な結果を得る事が出来ます。

    試しに何回か、メディアンフィルタ処理を縦方向に1列分の処理を行ってから、横方向に1列分の
    処理を行ってみましたが、ほぼ、良好な結果を得る事ができていると思います。
    (画像にするともう少し分かりやすいかと思いますが、プログラムが無いもので...)

    メディアンフィルタの例

    画像処理アルゴリズムへ戻る