www.pudn.com > 医学算法.rar > tof.c

#include	

main()
{
	int i,j,k,a[400];
	float b[400],th,dth;
	float t[256];
	float max,tt[256];

	fread(a,sizeof(*a),400,fopen("udata","r"));
	for(i=max=0;i<400;i++){
		b[i] = (a[i] < 0)? -a[i] : a[i];
		if(b[i] > max) max = b[i];
	}
	dth = max / 257;
	th = 0.0;

	for(i=0;i<256;i++){
		for(j=0;j<400;j++)if(b[j] > th)break;
		t[i] = j;
		tt[i] = th;
		th += dth;
	}
	fwrite(t,sizeof(*t),256,fopen("th","w"));
	fwrite(tt,sizeof(*tt),256,fopen("tt","w"));
	fwrite(b,sizeof(*b),400,fopen("b","w"));
}