質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

90.34%

C言語でガウシアンフィルタの実装をしたい

解決済

回答 1

投稿 編集

  • 評価
  • クリップ 1
  • VIEW 2,036

Trebey

score 5

 前提・実現したいこと

c言語で、(24bit)ビットマップを対象に、3*3のガウシアンフィルタを実装したいと考えています。
OpenCVなどの、画像処理用ライブラリは使わずに実装しようとしています。
画像の処理の流れとしては、
ビットマップ画像をグレースケール化→さらにガウシアンフィルタをかける
というものです。
ただし、端の4辺は、グレースケール化したものをそのまま出力するものとしています。

 発生している問題・エラーメッセージ

画像を読み込み、グレースケール化するところまではできている(と思います)。

しかし、ガウシアンフィルタをかけると、「ぼやける(平滑化される)」はずが、
左右に分身するような現象が起きてしまいます(「試したこと」欄に記載)。
コードに間違いがあると思いますので、ご指摘お願いいたします。

 該当のソースコード

/**********************************************
gaussian: ビットマップ画像のグレースケール化、ガウシアンフィルタ
概要:入力画像をグレースケール化し、ガウシアンフィルタを掛ける
**********************************************/
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <string.h>

int main(void){
//unsigned char img[512][512][3]; //読み込んだ画素情報を保存する(512*512まで可能)
unsigned char BitMapFileHeader[14]; //BMPのファイルヘッダーを保存する
unsigned int biSize; //BMPのサイズを保存する
int biWidth; //BMPの幅を保存する
int biHeight; //BMPの高さを保存する
unsigned char BitMapInfoHeader[28]; //上記3つ以外のBMPの情報ヘッダーを保存する
int i,j,c; //rgb入力for文用



char filename[50];
printf("please write a Bitmap filename(without [.bmp])\n");
scanf("%s",filename);

printf("\n");
if(filename==NULL){
    printf("please write a Bitmap filename(without [.bmp])\n");
    scanf("%s",filename);
}


char toGray[50];
char toGaussian[50];
strcpy(toGray,filename);
strcpy(toGaussian,filename);
FILE *fp; //ファイルポインタ
fp = fopen(strcat(filename,".bmp"),"rb");


//読み込み
fread(&BitMapFileHeader,sizeof(char),14,fp); //ファイルヘッダーを読み込む
fread(&biSize,sizeof(int),1,fp); //情報ヘッダーにあるサイズを読み込む
fread(&biWidth,sizeof(int),1,fp); //情報ヘッダーにある幅を保存
fread(&biHeight,sizeof(int),1,fp); //情報ヘッダーにある高さを保存
fread(&BitMapInfoHeader,sizeof(char),28,fp); //残りの情報ヘッダーを保存

//ビットマップは、横の幅が4の倍数になるまで0を加えるので、横幅を調整
while(1){
    if(biWidth%4 == 0){
        break;
    }
biWidth ++;
}

unsigned char img[biWidth][biHeight][3];//画像幅、高さ、rgb分を確保
unsigned char gray[biWidth][biHeight];//画像幅、高さ、rgb分を確保
unsigned char gaussian[biWidth][biHeight];//ガウシアンフィルタ処理後の値を格納 
for(i = 0; i < biWidth; i++){ //0から幅まで
    for(j = 0; j < biHeight; j++){ //0から高さまで
        for(c = 0; c < 3; c++){ //RGBのそれぞれ
            fread(&img[i][j][c],sizeof(char),1,fp); //画素の情報を読み込んで保存する

        }
    }
}

fclose(fp); //ファイルを閉じる

printf("width: %d\nheight: %d\n",biWidth,biHeight);

//グレースケール書き出し
char *grayfile = strcat(toGray,"_gray.bmp");
fp = fopen(grayfile,"wb"); //パスのファイルへ書き込み
fwrite(&BitMapFileHeader,sizeof(char),14,fp); //ファイルヘッダーを書き込む

fwrite(&biSize,sizeof(int),1,fp); //情報ヘッダーへファイルサイズを書き込む
fwrite(&biWidth,sizeof(int),1,fp); //情報ヘッダーへ幅を書き込む
fwrite(&biHeight,sizeof(int),1,fp); //情報ヘッダーへ高さを書き込む
fwrite(&BitMapInfoHeader,sizeof(char),28,fp); //残りの情報ヘッダーを書き込む


for(i = 0; i < biWidth; i++){ //0から幅まで
    for(j = 0; j < biHeight; j++){ //0から高さまで
        //gray = (img[i][j][0] + img[i][j][1] + img[i][j][2]) / 3;
        gray[i][j] = img[i][j][2] * 0.3 + img[i][j][1] * 0.59 + img[i][j][0] * 0.11;//ビットマップは、BGRの順で記述される

        for(c = 0; c < 3; c++){
            fwrite(&gray[i][j],sizeof(char),1,fp); //画素の情報を保存する
        }
    }
}


fclose(fp); //ファイルを閉じる


//ガウシアンフィルタ書き出し

int side;
int vertical;
char *bifile = strcat(toGaussian,"_gaussian.bmp");
fp = fopen(bifile,"wb"); //パスのファイルへ書き込み
fwrite(&BitMapFileHeader,sizeof(char),14,fp); //ファイルヘッダーを書き込む

fwrite(&biSize,sizeof(int),1,fp); //情報ヘッダーへファイルサイズを書き込む
fwrite(&biWidth,sizeof(int),1,fp); //情報ヘッダーへ幅を書き込む
fwrite(&biHeight,sizeof(int),1,fp); //情報ヘッダーへ高さを書き込む
fwrite(&BitMapInfoHeader,sizeof(char),28,fp); //残りの情報ヘッダーを書き込む


for(i = 0; i < biWidth; i++){ //0から幅まで
    for(j = 0; j < biHeight; j++){ //0から高さまで
    gaussian[i][j] = 0;//初期化
    }
}


for(i = 0; i < biWidth; i++){ //0から幅まで
    for(j = 0; j < biHeight; j++){ //0から高さまで
    //ガウシアンフィルタ

    if( (i>0) && (i<biWidth-1) && (j>0) && (j<biHeight-1) ){
    //gaussian[i][j] = gray[i][j]/16*4 + (gray[i-1][j-1]/16)*1 + (gray[i-1][j]/16)*2 + (gray[i-1][j+1]/16)*2 + (gray[i][j-1]/16)*2 + (gray[i][j+1]/16)*2 + (gray[i+1][j-1]/16)*1 + (gray[i+1][j]/16)*2 + (gray[i+1][j+1]/16)*1;

        //周辺のものと加算
        for(side = i-1; side < i+2; side++){
            for(vertical = j-1; vertical < j+2; vertical++){

                if( (side == i) && (vertical == j) ){
                    //自分自身を指す場合
                    gaussian[i][j] += (gray[side][vertical] /16) * 4;
                }
                else if( ( (side == i) && (vertical != j) ) || ( (side != i) && (vertical == j) ) ){
                    //上下左右のどれか
                    gaussian[i][j] += (gray[side][vertical] /16) * 2;

                }
                else if( (side != i) && (vertical != j) ){
                    //四隅のどれか
                    gaussian[i][j] += (gray[side][vertical] /16) * 1;
                }
            }
        }
    }
    else{//端の4辺
        gaussian[i][j] = gray[i][j];
    }

    for(c = 0; c < 3; c++){ //RGBのそれぞれ(順番は、[0]B [1]G [2]R)
        fwrite(&gaussian[i][j],sizeof(char),1,fp); //画素の情報を保存する
    }

    }
}

fclose(fp); //ファイルを閉じる

return 0;
}

試したこと

画像(jpeg)を示します。
元画像
イメージ説明
グレースケール化
イメージ説明
ガウシアンフィルタ(左右に分身がでてしまう)
イメージ説明

 補足情報

コンパイラ:gcc(win)
エディタ:Visual Studio Code

  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

回答 1

checkベストアンサー

+1

こんにちは。

ざっと見ただけなので外しているかもしれませんが、通常ビットマップ・ファイルは、横方向に走査されます。
それに対して、縦方向に走査したような処理になっています。縦横サイズが異なる画像ですので、操作方向の相違でズレているのかも知れません。
グレイスケール変換時は、1ピクセルだけで処理が完結するので走査方向の相違の影響は受けない筈です。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2018/10/26 00:28

    ようはとりあえずループ回数に指定されるheightとwidth入れ替えてみたら?ってことかな(適当)

    キャンセル

  • 2018/11/01 11:52

    読み込みと、フィルタをかけるところの走査方向を横向きに(行ごとに)したところ、それらしくなりました!ありがとうございました!

    キャンセル

15分調べてもわからないことは、teratailで質問しよう!

  • ただいまの回答率 90.34%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

同じタグがついた質問を見る