学校課題で与えられた画像にテンプレートと一致する部分がいくつあるか数え、また一致している部分のみ抽出してファイル出力せよという課題が
ありました。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(xpb/8%4!=0){
pad=(4-xpb/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;
}