- ベストアンサー
空間周波数特性について
- 空間周波数特性についての質問文章です。
- ローパスフィルタに関する質問です。空間周波数特性AとBについての影響や出力画像についての疑問を持っています。
- 質問内容は、f(x)をAのフィルタに通した後、さらにBのフィルタに通すとどうなるか、逆にf(x)をBのフィルタに通した後、さらにAのフィルタに通すとどうなるか、そしてf(x)をBのフィルタに通した後、さらにBのフィルタに通すとどうなるかというものです。
- みんなの回答 （2）
- 専門家の回答
質問者が選んだベストアンサー
(1) まず「フーリエ解析」のデータ列の点数は2のべき乗（・・・,32,64,128,256,・・・）でなくてはいけないので、ここでは例えばN=256にします。エクセルのワークシートのA列に縦に、0,1,2,・・・,254,255と入力します。 (2) f(x)は例えばL=64として、B列に縦に上から、0を16個、1を32個、0を16個入力し、この64点を縦方向に3回繰り返しコピーします（これが256点の元画像！）。 (3) B列に作った元画像データを「フーリエ解析」。「入力範囲」に"B1からB256まで"を指定（$B$1:$B$256）し、出力先に"C1"を指定（$C$1）して実行。すると、C列に縦に元画像の空間周波数スペクトルができます。 (4) D列に空間周波数フィルタH(u)のデータを入力。u_0は周期Lの繰り返し信号の周波数なので、ここではL=64点であり、256点内で4回繰り返す周波数に相当するので"4"、特性Aでは"u_0*20"なので、D1からD81まで81個1を入力、それ以降は0を入力。そして256点のうち最後の80点、すなわちD177からD256には再び1を入力（何故こんなことをするかはカテゴリー「数学」の質問No.1517764を参照）。 (5) E列に、C列の元画像のスペクトルとD列のフィルタをかけ算した結果を作る。ただし、C列は複素数なのでかけ算には関数"IMPRODUCT(x,y)"を使用。例えば、セルE1には、"=IMPRODUCT(C1,D1)"と入力。 (6) E1からE256のデータ列を(3)同様の手順で「フーリエ解析」して結果をF列に作る。これが"画像装置出力後の画像"になります。 (7) ただし、(6)でできたF列のデータはやっぱり複素数なので、関数"IMABS(x)"を使って実数化したデータ列をG列に作る。例えば、セルG1には、"=IMABS(F1)/256"と入力（何故256で割るかについてはカテゴリー「数学」の質問No.1488438を参照）。 これで、A列とG列を選択して「散布図」グラフにすれば、結果の画像が視覚化できます。 あとは、 (4)のフィルタデータを特性Aにしたときと特性Bにしたときの比較 結果画像を再び入力元画像に与えて、特性A→特性Bの順に処理したときと特性B→特性Aの順に処理したときの比較 などをやればいいでしょう。細かいところで何度も引っかかるかも知れませんが、試行錯誤で自己解決に励みましょう！ また結果の見方についても自分で良く考えてみましょう！
その他の回答 (1)
- paddler
- ベストアンサー率53% (176/330)
元画像のフーリエ変換に画像装置の周波数特性をかけ算して、その結果をもう一度 逆フーリエ変換すると画像装置を通った後の画像になることはお分かりですよね？ こんな一次元画像の範囲であれば、Excelに入っている「フーリエ解析」を使えば 簡単にシミュレーションできますよ！ 「Lやu」を具体的な数値にしてやってごらんなさい。人に聞くより１００倍良く理解 できること請け合いです！ "フーリエ解析"はメニューから、"ツール"－"分析ツール"－"フーリエ解析"に あります。ただし予めExcelには、"アドイン"で"分析ツール"がインストールされて いることが必要です。
お礼
ご回答ありがとうございます。 Excelでフーリエ解析をしようと試みたのですが、 どのように数値を入力していいのかが分かりませんでした。 Excelのヘルプを見てもさっぱり分かりません。 この問題ではどのようにExcelを使用したらよろしいでしょうか？ 重ね重ね申し訳ありませんが、ご教授ください。
お礼
エクセルの使い方の説明ありがとうございます。 ようやく空間周波数の特性を理解することができました。本当に有難うございました。