MATLABにおける信号処理のコード

Question

###前提・実現したいこと
MATLABで書かれた画像からより高い周波数の成分を抑制するFIRフィルタに関するコードのフィルタの係数(タップ数)を変更しようとしています。

###発生している問題・エラーメッセージ
デジタル信号処理とコードの理解がそれぞれ十分ではないため、フィルタの係数を減らすために以下のようにコードを書き換えると、画像の色が濃くなるという状況が起きています。

%  ローパスフィルターを画像に適用する
picture = imread('sample.bmp');
[a,b,c] = size(picture);

% フィルター係数4つ（真ん中の2つを削除)
fir = [];
dv = double(zeros(a,b,c));

picture = uint8((dv+dh)/2.0);
image(picture);

###該当のソースコード

%  ローパスフィルターを画像に適用する
picture = imread('sample.bmp');
[a,b,c] = size(picture);

% フィルター係数4つ
fir = [];
dv = double(zeros(a,b,c));
dh = double(zeros(a,b,c));
for i=4:a-3
    for j=4:b-3
        for k=1:c
 
        end
    end
end
% 符号なしのint型をもつと考えられる
picture = uint8((dv+dh)/2.0);
image(picture);

Comments

[ozwk]

元のコード（係数6）がrのループでr使ってないなど、ちゃんと動いてなさそうなんですが、大丈夫ですか？

[退会済みユーザー]

ご指摘いただきましてありがとうございます。元のコードはMATLABで動くことを確認しました。

[ozwk]

6つの方の係数どうやって決めました？

[退会済みユーザー]

詳しくはわからないのですが、MATLABのフィルターデザイナーかと思われます。

Answer 1

元の画像の全部の画素値が1だったとしましょう
係数6つ版だったらフィルターかけると出力画像の各画素の値は大体1.1ちょいになりますが
4つ版だと0.1ぐらいですね。
暗くなるのはそういう理由です。

均一な画像にフィルターかけたとき1倍になるように係数を調整しましょう

Comments

[退会済みユーザー]

ご回答いただきましてありがとうございます。
係数6つ版の「1.1」という数字はどのように計算されたのか教えていただけますと幸いです。

[ozwk]

どのようにも何もコードに沿って計算したんですが
そうなりませんか?

[退会済みユーザー]

MATLABでコードを実行したところ、ワークスペースには1.1という数字は見当たらず、「1,1」という数字は何を指し示しているのか教えていただきたいです。

[退会済みユーザー]

仮に「1.1」という数字が係数を足して出された値として、試したことの追記を行いました。どうぞよろしくお願いいたします。

[ozwk]

出力結果に満足してるならそれでいいんじゃないですかね
満足してないなら理想と現状の差が書いてないので答えようがないです

[退会済みユーザー]

デジタル信号処理に関して理解が及んでおらず、このように安易にコードを書き換えただけで正しく「2次元FIRフィルタを適用して画像の高周波成分を抑圧」できているのか分からない状態です。間違いであれば、ご指摘いただきたいです。

[ozwk]

とりあえず1次元から勉強を始めたらいいと思います。

[退会済みユーザー]

ご回答とご説明ありがとうございました。最後に、「1.1」という値が何を指し示しているのかお聞きできましたら幸いです。長期的に勉強をしますが、このコードに関しては早急に理解する必要があるため、お手数ですが、教えていただけますと大変助かります。

[ozwk]

全画素値がxである画像をこのフィルター通したら各画素値が何倍になるかです。

[退会済みユーザー]

ありがとうございました。