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

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

ただいまの
回答率

87.49%

c言語でのテンプレートマッチングについて

受付中

回答 1

投稿 編集

  • 評価
  • クリップ 0
  • VIEW 1,633
退会済みユーザー

退会済みユーザー

学校課題で与えられた画像にテンプレートと一致する部分がいくつあるか数え、また一致している部分のみ抽出してファイル出力せよという課題が
ありました。SSD法での解決を試みたのですが、どうやらssdの計算の部分が上手く進んでいないようです。自分でどこがおかしいか気づけなかったので
ここで質問致します。また、その他おかしい点があればご指摘ください。
以下がコードです。

include <stdio.h>

include <stdlib.h>

typedef struct bmp{
short int r;
short int g;
short int b;
short int t;
} bmp;

int main(void) {
FILE *fp,*fp2;
bmp pic,pic2,**sad;
char data[50];
int i,j,w=255;
int x,y,ofset,pad=0,x2,y2,pad2=0;
int w=255;
short int pb,pb2;
printf("file name:");
scanf("%s",data);
printf("%s\n",data);
if((fp=fopen(data,"rb+"))==NULL){
printf("File not found.\n");
return 0;
}
fseek(fp,0x0012,SEEK_SET);
fread(&x,4,1,fp);
fread(&y,4,1,fp);
fseek(fp,0x001c,SEEK_SET);
fread(&pb,2,1,fp);
printf("Horizontal size=%d  Vertical size=%d,\n",x,y);
printf("Bit length per pixel=%d\n",pb);
fseek(fp,0x000a,SEEK_SET);
fread(&ofset,4,1,fp);
printf("offset:%d\n",ofset);
if(x*pb/8%4!=0){
pad=(4-x*pb/8%4)%4;
}

fseek(fp,ofset,SEEK_SET);

pic=(bmp **)malloc(sizeof(bmp *)*y);
for(i=0;i<y;i++){
pic[i]=(bmp *)malloc(sizeof(bmp)*x);
}
for(i=0;i<y;i++){
for(j=0;j<x;j++){
fread(&pic[i][j].b,1,1,fp);
fread(&pic[i][j].g,1,1,fp);
fread(&pic[i][j].r,1,1,fp);

if((pic[i][j].b==255)&&(pic[i][j].g==255)&&(pic[i][j].r==255)){
pic[i][j].t=1;
}
else{
pic[i][j].t=0;
}
}
fseek(fp,pad,SEEK_CUR);
}

if((fp2=fopen("example256_2.bmp","rb+"))==NULL){
printf("File not found.\n");
return 0;
}
fseek(fp2,0x0012,SEEK_SET);
fread(&x2,4,1,fp2);
fread(&y2,4,1,fp2);
fseek(fp2,0x001c,SEEK_SET);
fread(&pb2,2,1,fp2);
printf("bh265.bmp:Horizontal size=%d  Vertical size=%d,\n",x2,y2);
printf("Bit length per pixel=%d\n",pb2);
fseek(fp2,0x000a,SEEK_SET);

if((x2)*(pb2)/8%4!=0){
pad2=(4-(x2)*(pb2)/8%4)%4;
}

fseek(fp2,ofset,SEEK_SET);

pic2=(bmp **)malloc(sizeof(bmp *)*y2);
for(i=0;i<y2;i++){
pic2[i]=(bmp *)malloc(sizeof(bmp)*x2);
}
for(i=0;i<y2;i++){
for(j=0;j<x2;j++){
fread(&pic2[i][j].b,1,1,fp2);
fread(&pic2[i][j].g,1,1,fp2);
fread(&pic2[i][j].r,1,1,fp2);
pic2[i][j].t=0;
}

fseek(fp2,pad2,SEEK_CUR);
}
fseek(fp2,ofset,SEEK_SET);

int l,m,c=0;
bmp min;
min.b=1000;
min.g=1000;
min.r=1000;

sad=(bmp **)malloc(sizeof(bmp *)*y2);
for(i=0;i<y2;i++){
sad[i]=(bmp *)malloc(sizeof(bmp)*x2);
}

for(l=0;l<y2-y;l++){
for(m=0;m<x2-x;m++){
for(i=0;i<y;i++){
for(j=0;j<x;j++){
sad[l][m].b=abs(pic2[i+l][j+m].b-pic[i][j].b);
sad[l][m].b=abs(pic2[i+l][j+m].g-pic[i][j].g);
sad[l][m].b=abs(pic2[i+l][j+m].r-pic[i][j].r);
printf("test");
}
}
}
}

for(i=0;i<y2;i++){
for(j=0;j<x2;j++){
if(sad[i][j].b<=min.b){
min.b=sad[i][j].b;
}
if(sad[i][j].g<=min.g){
min.g=sad[i][j].g;
}
if(sad[i][j].r<=min.r){
min.r=sad[i][j].r;
}
}
}

for(l=0;l<y2;l++){
for(m=0;m<x2;m++){
if((sad[l][m].b==min.b)&&(sad[l][m].g==min.g)&&(sad[l][m].r==min.r)){
c++;

for(i=0;i<y;i++){
for(j=0;j<x;j++){
if((pic[i][j].b==pic2[i+l][j+m].b)&&(pic[i][j].g==pic2[i+l][j+m].g)&&(pic[i][j].r==pic2[i+l][j+m].r)&&(pic2[i+l][j+m]==0)){
pic2[i+l][j+m].t=1;
}
}
}

}
}
}

fseek(fp2,ofset,SEEK_SET);

for(i=0;i<y2;i++){
for(j=0;j<x2;j++){
if(pic2[i][j].t==1){
fwrite(&pic2[i][j].b,1,1,fp2);
fwrite(&pic2[i][j].g,1,1,fp2);
fwrite(&pic2[i][j].r,1,1,fp2);
}
else{
fwrite(&w,1,1,fp2);
fwrite(&w,1,1,fp2);
fwrite(&w,1,1,fp2);  
}
}
}

printf("number of mushrooms :%d",c);

fclose(fp);
fclose(fp2);

free(pic);
free(pic2);
free(sad);
return 0;
}

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

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

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

    クリップを取り消します

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

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

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

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

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

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

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

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

    質問の評価を下げる

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

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

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

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

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

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

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

    詳細な説明はこちら

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

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

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

質問への追記・修正、ベストアンサー選択の依頼

  • y_waiwai

    2018/08/03 07:55

    コードは、編集ボラン押して<code>ボタン、’’’の枠に中にコピーしてください

    キャンセル

  • rocket

    2018/09/13 00:50 編集

    <code>ボタンの事など、https://teratail.com/help#about-markdown ここのヘルプを見ると分かると思います

    キャンセル

回答 1

0

ざっくり見て気になった点。

  • pic pic2は実体で宣言されているのに、mallocで領域を割り当てようとしている。
  • RGBをfreadで1バイトずつ読み込んでいるが、構造体の定義はshort型。
    例えば0xFFを読み込んだ場合、0xFF??という値になり、??の部分は初期化していないので不定値です。
    たまたま0が入っていてリトルエンディアンなら数値化したときに255になるが、0以外の値だと想定した値にはならない。
    ビッグエンディアンなら0xFF00がそのまま数値化されるので、65280という値になりif文で正しく判定できない。
    unsigned charで宣言するべきかと思われます。

SSD法とかは知らんのであとは先生に聞きましょう。

投稿

  • 回答の評価を上げる

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

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

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

  • 回答の評価を下げる

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

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

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

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

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

関連した質問

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