特定の順番で関数を呼び出すとセグフォが起きる

Question

前提・実現したいこと

cで2値画像のラベリングをし，ラベリングをもとに画像からノイズを消去しようとしています．
ラベリングの関数が完成したので，ノイズを消去する関数を作成し，その関数内からラベリング関数を呼び出そうとすると，segmentation faultが発生しました．
どうやらラベリング関数内でセグフォは起きているようですが，原因がわかりません．

64*64で0,1で構成された画像をcharの二次元配列に読み込み，ラベリングでは画像の1となっているピクセルを2以上の数字でラベリングしています．ノイズの消去ではラベリングしたものをラベルごとに見て，ピクセル数が5以下のラベルのパーツをノイズとみなして消去し，最後にラベリングされたピクセルを1に戻すという作業を行なっています．

該当のソースコード

#include <stdio.h>
#include <stdlib.h>

#define	IMGSIZE	512	
#define NOISESIZE 5

void PrintImg( unsigned char data[512] )
{
	unsigned char mask;
	for(int x=0;x<64;x++){
		for(int y=0;y<8;y++){
			mask=0x80;
			for(int i=0;i<8;i++){
				if(data[8*x+y] & mask)
					printf("*");
				else
					printf("o");

				mask=mask>>1;
			}
		}
		printf("\n");
	}

}

void ExpandImg(unsigned char data[512],char img[64][64]){
    	unsigned char mask;
	for(int y=0;y<64;y++){
		for(int x=0;x<8;x++){
			mask=0x80;
			for(int i=0;i<8;i++){
                if(data[8*y+x] & mask)
				    img[y][8*x+i]=1;
                else
                    img[y][8*x+i]=0;
				mask=mask>>1;
			}
        }
	}
}

void PrintPattern( char img[64][64] )
{
	int x, y, i ;

	for( y = 0 ; y < 64 ; y++ ) {
		for( x = 0 ; x < 64 ; x++ ) {
			putchar( img[y][x] ? '*' : 'o' ) ;
		}
		putchar( '\n' ) ;
	}
}

void NormalizeImg(char img[64][64]){
    int i=0;
    int j=0;
    int origin_x=-1;
    int origin_y=-1;
    int width=-1;
    int height=-1;

    for(i=0;i<64;i++){
        for(j=0;j<64;j++){
            if(img[j][i]==1 && origin_x<0)
                origin_x=i;

            if(img[i][j]==1 && origin_y<0)
                origin_y=i;

            if(origin_x>=0 && origin_y>=0)
                break;
        }
    }

    for(i=(64-1);i>=0;i--){
        for(j=0;j<64;j++){
            if(img[j][i]==1 && width<0)
                width=i-origin_x;

            if(img[i][j]==1 && height<0)
                height=i-origin_y;

            if(width>0 && height>0)
                break;
        }
    }

    double x_ratio=0.0;
    double y_ratio=0.0;
    int point_x=0;
    int point_y=0;
    char img_normalized[64][64];

    x_ratio=(double)width/64.0;
    y_ratio=(double)height/64.0;

    for(i=0;i<64;i++){
        for(j=0;j<64;j++){
            point_x=origin_x+(int)((double)j*x_ratio+0.5);
            point_y=origin_y+(int)((double)i*y_ratio+0.5);
            img_normalized[i][j]=img[point_y][point_x];
        }
    }

    for(i=0;i<64;i++)
        for(j=0;j<64;j++)
            img[i][j]=img_normalized[i][j];
}

int label(char img[64][64]){
    unsigned char last_label=1;
    unsigned char pix=0;
    unsigned char op_label=0;
    unsigned char max_label=0;
    unsigned char table[255];
    char unused_labels[4];
    int i,j,k;

    for(i=0;i<255;i++)
        table[i]=i;

    for(i=0;i<64;i++){
        for(int j=0;j<64;j++){
            if(img[i][j]==1){
                for(k=0;k<4;k++) unused_labels[k]=0;
                if((i-1)>=0){
                    k=0;
                    pix=img[i-1][j];
                    if(pix>1){
                        if((op_label!=0 && op_label>pix) || op_label==0){
                            op_label=pix;
                            unused_labels[k]=pix; k++;
                        }
                    }


                    if((j-1)>=0){
                        pix=img[i-1][j-1];
                        if(pix>1){
                            if((op_label!=0 && op_label>pix) || op_label==0){
                                op_label=pix;
                                unused_labels[k]=pix; k++;
                            }
                        }
                    }

                    if((j+1)<64){
                        pix=img[i-1][j+1];
                        if(pix>1){
                            if((op_label!=0 && op_label>pix) || op_label==0){
                                op_label=pix;
                                unused_labels[k]=pix; k++;
                            }
                        }
                    }
                }

                if((j-1)>=0){
                    pix=img[i][j-1];
                    if(pix>1){
                        if((op_label!=0 && op_label>pix) || op_label==0){
                            op_label=pix;
                            unused_labels[k]=pix; k++;
                        }
                    }
                }

                if(op_label==0){
                    op_label=last_label+1;
                    if(op_label>=255)
                        return 1;
                }else{
                    for(k=0;k<4;k++){
                        if(unused_labels[k]==0) break;
                        table[unused_labels[k]]=op_label;
                        if(unused_labels[k]>max_label)
                            max_label=unused_labels[k];
                    }
                }

                img[i][j]=op_label;
                last_label=op_label;
                
                op_label=0;
            }
        }
    }

    unsigned char origin_label;
    for(k=2;k<=max_label;k++){
        last_label=table[k-1];
        if(table[k]!=(last_label+1) || table[k]!=last_label){
            origin_label=table[k];
            i=k;
            while (table[i]==origin_label){
                table[i]=last_label+1;
                i++;
            }
            k=i-1;
        }
    }

    for(i=0;i<64;i++)
        for(j=0;j<64;j++)
            img[i][j]=table[img[i][j]];

    return 0;
}

void NoiseErase(char p[64][64],const int size){
    if(label(p)!=1){
        int label=2;
        int max_label=0;
        int i,j;
        int num=0;
        int *erase_label;
        int count=0;

        for(i=0;i<64;i++){
            for(j=0;j<64;j++){
                if(p[i][j]>max_label)
                     max_label=p[i][j];
            }
        }
        erase_label=(int*)malloc(max_label*sizeof(int));
        for(i=0;i<max_label;i++)
            erase_label[i]=0;
        for(;label<=max_label;label++){
            for(i=0;(i<64 || num>size);i++){
                for(j=0;(j<64 || num>size);j++){
                    if(p[i][j]==label)
                        num++;
                }
            }
            if(num<=size){
                erase_label[count]=label;
                count++;
            }
            num=0;
        }
        for(i=0;i<64;i++){
            for(j=0;j<64;j++){
                if(p[i][j]>=2){
                    for(int k=0;k<count;k++){
                        if(p[i][j]==erase_label[k]){
                            p[i][j]=0;
                            break;
                        }
                    }
                    if(p[i][j]>=2)
                        p[i][j]=1;
                }
            }
        }
        free(erase_label);
    }
}


int main()
{
    FILE		*fp ;
    char		imgfile[0x80] ;
    unsigned char	data[IMGSIZE] ;
    char img[64][64];

    printf( "画像ファイル名を入力してください:" ) ;
    fscanf( stdin,"%s",imgfile ) ;

    if( ( fp = fopen( imgfile, "rb" ) ) == NULL ) {
	fprintf( stderr, "cannot open %s\n", imgfile ) ;
	return 1 ;
    }

    while( fread( data, IMGSIZE, 1, fp ) ) {
        PrintImg( data ) ;
		ExpandImg( data, img ) ;
		NoiseErase( img, NOISESIZE ) ;
		NormalizeImg( img ) ;
		//PrintPattern( img ) ;
		printf( "次のデータを表示するにはリターンキーを押してください" ) ;
		getchar( ) ;    }

    fclose( fp ) ;
    return 0 ;
}

試したこと

label()を呼び出す部分をNoiseErase()からコメントアウトし，label()もmain()から呼び出してみましたが，やはりlabel()内でセグフォが起きました．
ただ，NoiseErase()を先に呼び出すと問題なく動くようです．
NoiseErase()単体でも動きます．
もちろんそれではやらせたいことは実現できませんが．．．

補足情報（FW/ツールのバージョンなど）

gccでmacOS上でbuildしています

Comments

[izmktr]

パット見た感じ、table[img[i][j]] がきになりますね
255を指定しているってことは、imgはunsigned charを想定しているけど、charだから負の値になってる？
(unsigned charの全域をカバーするなら、256を指定するべきですね)

[SHOMI]

このままではデバッグ丸投げです。
まずはご自分でデバッグしてみてください。
「gdb segmentation fault」や「gdb SEGV」で検索すれば方法は出てきます。

Answer 1

各所に printf と fflash 等を組み合わせたログ表示を仕込み, どこで起きているかを確認されては如何でしょうか.

Answer 2

ぱっと目についたところで、

関数 void ExpandImg(unsigned char data[512],char img[64][64])
ここの img[y][8*x+i]=1; 、　xの最大値が 8なので、[8*x+i] > 64 となりませんか?

他にもありそうですが、取り合えず。