各種の周期ムラの2次​元画像を作成して、F​FT変換して、対応を​学習したい

周期性のムラのある2次元画像をMatlabで生成して、その画像をFFT変換します。そのFFT変換画像と元の周期性画像との関係を比較して学習したいです。各種、周期性ムラの生成方法と、FFT変換画像の作成の仕方をご教授願いたいです。

[Atsushi Ueno]

> 周期性ムラの生成方法

I = imread('peppers.png');
S = size(I);
func = (sin(10 * pi * (1/S(2):1/S(2):1)) + 3) / 4; % 画像幅と同じ長さの正弦波
strp = repmat(func,[S(1),1,3]); % 画像と同じサイズの縞模様データを作る
I = uint8(double(I) .* strp); % 画像と縞模様の要素を掛け合わせる
subplot(1,2,1), imshow(I);
subplot(1,2,2), imshow(strp);

[Hirofumi]

Uenoさん

ありがとうございます。確認できました。頂いた例をもとに、FFT変換、MTFカーブ算出を検討します。時間が許すなら、下記のレクチャもよろしくお願いします。

X,Y方向のFFTを求めて、画像のムラがX,Y方向のどの周波数のモノかを求めて、その周波数を改善するとMTFカーブどう変わるかを試そうとしています。

[Atsushi Ueno]

> 時間が許すなら、下記のレクチャもよろしくお願いします

時間はいくらでも許されるんですけどね。なんか fft2 関数に画像を入力してみたりしてるんですが、なんか良く分かんないんですよね。わかったら回答に追記します。