//?FCM_S.cpp:?主項目文件。
//iamzxf@126.com
//QQ?250129131

#include?“stdafx.h“
#include?“cv.h“
#include?“highgui.h“
#include?“cxcore.h“
#include?“time.h“
#include?“WTypes.h“


using?namespace?System;

#define?height?184
#define?width?180
#define?cluster_n?2
#define?LEVEL?256
#define?iters?100
#define?mm?1.75
#define?thres?0.00001
#define?alpha?2
#define?p?0
#define?q?1

double?membership[height][width][cluster_n];
double?center[cluster_n];
double?distance[iters];
double?temp_sumw;



int?main（array?^args）
{
	IplImage?*pImg*newImg;

	pImg=cvLoadImage（“tu2.jpg“0）;

	if（!pImg）{
		printf（“cannot?load?the?file.\n“）;
		return?-1;
	}

	uchar?*data*newdata;
	int?channels?stepdepth;	


	depth=pImg->depth;
	step=pImg->widthStep;
	channels=pImg->nChannels;
	data=（uchar?*）?pImg->imageData;

	int?ijkmnkk;
	int?member_sum;
	double?temp_sum;
	int?it;
	SYSTEMTIME?sys;?

	GetLocalTime（?&sys）;?	
	printf（“\nthe?time?is?%d%d.%d\n“sys.wMinutesys.wSecondsys.wMilliseconds）;

	//initialize?membership?matrix	
	srand（（int）time（0））;

	FILE?*pF*fF;
	pF=fopen（（char?*）“data1.txt““w“）;
	//fF=fopen（（char?*）“data2.txt““w“）;
	if?（NULL==pF）
	{
		printf（“open?write?file?error!!“）;
		return?1;
	}
	
	

	for（i=0;i	{
		for（j=0;j		{
			member_sum=0;
			for（k=0;k			{
				membership[i][j][k]=rand（）+1;
				member_sum+=membership[i][j][k];
			}
			
			for（k=0;k			{
				membership[i][j][k]/=member_sum;
			}
		}
	}
/*	for（i=178;i<=180;i++）
	{
		
		for（j=140;j<=142;j++）
		{
	
			fprintf（fF?“?%03u?“data[i*step+j]?）;
			
		}
		fprintf（fF?“\n“）;
		
	?}
	fclose（fF）;
	*/
	for（i=178;i<=180;i++）
	{
		for（j=140;j<=142;j++）
		{
	
			fprintf（pF?“%lf\n“?membership[i][j][0]）;
			
		}
		
	?}
	fprintf（pF?“\n“）;

	//compute?the?initialized?cluster?centers

	double?sum_downsum_upsum_gray;
	int?num_neigh;
	int?updownleftright;

	for（k=0;k	{
		sum_down=0;
		sum_up=0;
		for（i=0;i		{
			for（j=0;j			{
				sum_down+=（1+alpha）*pow（membership[i][j][k]mm）;

				num_neigh=0;
				sum_gray=0;

				up=i-1;
				down=i+1;
				left=j-1;
				right=j+1;

				if（up<0）?up=0;
				if（down>height）?down=height;
				if（left<0）?left=0;
				if（right>width）?right=width;

				for（m=up;m<=down;m++）
				{
					for（n=left;n<=right;n++）
					{
						sum_gray+=data[m*step+n];
						num_neigh++;
					}
				}

				sum_up+=pow（membership[i][j][k]mm）*（data[i*step+j]+（sum_gray-data[i*step+j]）*alpha/（num_neigh-1））;	

			}
		}

		center[k]=sum_up/sum_down;	
	}
	
	for（i=0;i	{
		printf（“center[%d]=%f?\n“i+1center[i]）;
	}
	//compute?the?distance[0]
	
	for（it=0;it	{
		distance[it]=0;
	}


	for（i=0;i	{
		for（j=0;j		{
			for（k=0;k			{
				distance[0]+=pow（membership[