RIGHT:寄稿：東條遼平 * 2次元フーリエ変換を使う [#a466e092] 縦でも横でもいいのですが、画像の1行もしくは1列を時間で変化する波と考えると画像に対しフーリエ変換をすることができそうです。画像の縦もしくは横がnピクセルであれば、その行や列は #ref(img18.png) の関数の離散値であると考えられそうです。とすれば画像のRGBデータに対し 2次元フーリエ変換を掛けてやることで周波数データにすることができます。フーリエ変換については1次元フーリエ変換、 2次元フーリエ変換を参照してください。 の関数の離散値であると考えられそうです。とすれば画像のRGBデータに対し 2次元フーリエ変換を掛けてやることで周波数データにすることができます。フーリエ変換については [[1次元フーリエ変換]] 、 [[2次元フーリエ変換]] を参照してください。 周波数データにしてなにが嬉しいかというと、この変換をすることで 2次元配列は次のような意味を持ちます。 #ref(bmpfourier.png) 左上が配列の(0,0)右下を配列の(n,n)と考えてください。ざっくりと考えて上の画像で示す中央程、低周波がどれほどの重みを持っているのか、外側程、高周波がどれほどの重みを持っているのかという事になります。そのことを利用すれば様々な事ができるのですが、加工は次に回すとして、今回は変換してそのまま戻すという作業をしたいと思います。ちなみに理論上は可逆変化なので全く元通りになるハズですが、計算途中で誤差が出るため実際には不可逆になります。 void RgbFrequencyTrans(Image *img, double *re, double *im, int flag) { int i, j; if(flag==DFT){ for(i=0; i<img->height*img->width; i++){ re[i] = img->data[i].r; im[i] = 0; } } dft_swap2(re, im, img->width, img->height); dft_idft2(re, im, img->width, img->height, flag); dft_swap2(re, im, img->width, img->height); if(flag==IDFT){ for(i=0; i<img->height*img->width; i++){ if(re[i]<0) re[i] = 0; if(re[i]>255) re[i] = 255; img->data[i].r = img->data[i].g = img->data[i].b = re[i]; } } } 今回作った関数があまりに短いので全部載せています。既にフーリエ変換部分は完成しているため簡単に実装できました。 dft_swap2などについては2次元フーリエ変換 に詳しく書いています。特に言う事はないのですが、強いて言えば、 double型からunsigned char型に変換する際に0から255の間に入るようにしてやる必要があることです。また、今回この関数はヘッダに宣言せずに、ソースファイルの中で宣言しています。これはこの関数がソースファイル内で使う為に作られた関数であることをアピールするためにしています。今回は実際に動かす為にmainで使っていますが、以降は外部からは直接使えないようにstaticを指定することにします。 今回作った関数があまりに短いので全部載せています。既にフーリエ変換部分は完成しているため簡単に実装できました。 dft_swap2などについては [[2次元フーリエ変換]] に詳しく書いています。特に言う事はないのですが、強いて言えば、 double型からunsigned char型に変換する際に0から255の間に入るようにしてやる必要があることです。また、今回この関数はヘッダに宣言せずに、ソースファイルの中で宣言しています。これはこの関数がソースファイル内で使う為に作られた関数であることをアピールするためにしています。今回は実際に動かす為にmainで使っていますが、以降は外部からは直接使えないようにstaticを指定することにします。 - &ref(main.c); - &ref(bitmap.c); - &ref(bitmap.h); - &ref(calculation.c); - &ref(calculation.h);
