#include<stdio.h>
#include"hensu.h"

#define L_BASE 100

void labelset(LPDWORD image,int xs,int ys,int label)
{
	int x,y,cnt,xm,xp,ym,yp;

	image[xs+ys*iWidth]=RGB(label,label,label);
	for(;;){
		cnt=0;
		for(y=0;y<iHeight;y++) 			for(x=0;x<iWidth;x++) 				if(image[x+y*iWidth]==RGB(label,label,label)){ 					xm=x-1; xp=x+1; ym=y-1; yp=y+1; 					if(xm<0) xm=0; if(xp>iWidth-1) xp=iWidth-1;
					if(ym<0) ym=0; if(yp>iHeight-1) yp=iHeight-1;
					if(image[xm+ym*iWidth]==HIGH){
						image[xm+ym*iWidth]=RGB(label,label,label);
						cnt++;
					}
					if(image[x+ym*iWidth]==HIGH){
						image[x+ym*iWidth]=RGB(label,label,label);
						cnt++;
					}
					if(image[xp+ym*iWidth]==HIGH){
						image[xp+ym*iWidth]=RGB(label,label,label);
						cnt++;
					}
					if(image[xm+y*iWidth]==HIGH){
						image[xm+y*iWidth]=RGB(label,label,label);
						cnt++;
					}
					if(image[xp+y*iWidth]==HIGH){
						image[xp+y*iWidth]=RGB(label,label,label);
						cnt++;
					}
					if(image[xm+yp*iWidth]==HIGH){
						image[xm+yp*iWidth]=RGB(label,label,label);
						cnt++;
					}
					if(image[x+yp*iWidth]==HIGH){
						image[x+yp*iWidth]=RGB(label,label,label);
						cnt++;
					}
					if(image[xp+yp*iWidth]==HIGH){
						image[xp+yp*iWidth]=RGB(label,label,label);
						cnt++;
					}
				}
		if(cnt==0) break;
	}
}

int labeling(LPDWORD image_in,LPDWORD image_label,int *cnt,char *buf)
{
	int x,y,label;

	for(y=0;y<iHeight;y++) 		for(x=0;x<iWidth;x++) 			image_label[x+y*iWidth]=image_in[x+y*iWidth]; 	label=L_BASE; 	for(y=0;y<iHeight;y++) 		for(x=0;x<iWidth;x++) 			if(image_label[x+y*iWidth]==HIGH){ 				if(label>=(BYTE)HIGH){
					sprintf(buf,"--Error! too many labels.\n");
					return -1;
				}
				labelset(image_label,x,y,label);
				label++;
			}
	*cnt=label-L_BASE;
	sprintf(buf,"  number of labels : %d \n",*cnt);
	return 0;
}

#define PI 3.14159265
#define ROOT2 1.41421356

double calc_size(LPDWORD image_label,int label,int *cx,int *cy)
{
	int x,y;
	double tx,ty,total;

	tx=0; ty=0; total=0;
	for(y=0;y<iHeight;y++) 		for(x=0;x<iWidth;x++) 			if(image_label[x+y*iWidth]==RGB(label,label,label)){ 				tx+=x; ty+=y; total++; 			} 	if(total==0) return 0; 	*cx=(int)(tx/total); *cy=(int)(ty/total); 	return total; }  double trace(LPDWORD image_label,int xs,int ys) { 	int x,y,vec; 	double l; 	DWORD no;  	l=0; x=xs; y=ys; no=image_label[x+1+y*iWidth]; vec=5; 	for(;;){ 		if(x==xs && y==ys && l!=0) return l; 		image_label[x+y*iWidth]=HIGH; 		switch(vec){ 			case 3: 				if(image_label[x+1+y*iWidth]!=no && image_label[x+1+(y-1)*iWidth]==no){ 					x=x+1; y=y; l++; vec=0; continue;} 			case 4: 				if(image_label[x+1+(y-1)*iWidth]!=no && image_label[x+(y-1)*iWidth]==no){ 					x=x+1; y=y-1; l+=ROOT2; vec=1; continue;} 			case 5: 				if(image_label[x+(y-1)*iWidth]!=no && image_label[x-1+(y-1)*iWidth]==no){ 					x=x; y=y-1; l++; vec=2; continue;} 			case 6: 				if(image_label[x-1+(y-1)*iWidth]!=no && image_label[x-1+y*iWidth]==no){ 					x=x-1; y=y-1; l+=ROOT2; vec=3; continue;} 			case 7: 				if(image_label[x-1+y*iWidth]!=no && image_label[x-1+(y+1)*iWidth]==no){ 					x=x-1; y=y; l++; vec=4; continue;} 			case 0: 				if(image_label[x-1+(y+1)*iWidth]!=no && image_label[x+(y+1)*iWidth]==no){ 					x=x-1; y=y+1; l+=ROOT2; vec=5; continue;} 			case 1: 				if(image_label[x+(y+1)*iWidth]!=no && image_label[x+1+(y+1)*iWidth]==no){ 					x=x; y=y+1; l++; vec=6; continue;} 			case 2: 				if(image_label[x+1+(y+1)*iWidth]!=no && image_label[x+1+y*iWidth]==no){ 					x=x+1; y=y+1; l+=ROOT2; vec=7; continue;} 				vec=3; 		} 	} }  double calc_length(LPDWORD image_label,int label) { 	int x,y;  	for(y=0;y<iHeight;y++) 		for(x=0;x<iWidth;x++) 			if(image_label[x+y*iWidth]==RGB(label,label,label)) 				return trace(image_label,x-1,y); 	return 0; }  void features(LPDWORD image_label_in,LPDWORD image_label_out, 			  int cnt,double size[],double ratio[],char *buf) { 	int i,x,y,center_x,center_y; 	double l; 	int posi,m;  	posi=0; 	for(y=0;y<iHeight;y++) 		for(x=0;x<iWidth;x++) 			image_label_out[x+y*iWidth]=image_label_in[x+y*iWidth]; 	m=sprintf(buf," 番号  面積  周囲長  円形度  重点(x,y)\n"); 	posi+=m; 	for(i=0;i<cnt;i++){ 		size[i]=calc_size(image_label_out,i+L_BASE,&center_x,&center_y); 		l=calc_length(image_label_out,i+L_BASE); 		ratio[i]=4*PI*size[i]/(l*l); 		image_label_out[center_x+center_y*iWidth]=HIGH; 		m=sprintf(&buf[posi],"  %3d %6d %8.2lf %8.4lf (%3d,%3d)\n", 			i,(int)size[i],l,ratio[i],center_x,center_y); 		posi+=m; 	} }  void ratio_extract(LPDWORD image_label_in,LPDWORD image_label_out, 				   int cnt,double ratio[],double ratio_min,double ratio_max) { 	int i,j,x,y; 	int lno[256];  	for(i=0,j=0;i<cnt;i++) 		if(ratio[i]>=ratio_min && ratio[i]<=ratio_max)
			lno[j++]=L_BASE+i;
	for(y=0;y