智能平台实验报告2.docxVIP

实验 2 基于 FCM 的图像分割黄灏卓远 一、实验目的： 通过 FCM 实现图像的分割。 二、流程图： 开始 开始 在 0,1 间随机初始化隶属度矩阵 U ， 使 得 满 足 ?c i?1 U ? 1, ( j ? 1,2...N ) ij 计算C 个聚类中心： ?n c ? i u m x ij j j ?1 ?n j ?1 u m ij 计 算 目 标 函 数 j (U , c .... c ) ? ?? c n 1 2 u d ij 2 ij i ? 1 j 如果目标函数是否小于阈 值，或对上此目标函数改变是否小于某个阈值。 是 否 重新计算新的 U 重新计算新的 U 矩阵 结束 三．具体代码 /*bmpimage.imgbuf 为图像数据首地址，指针类型*/ /*图像宽度变量为 bmpimage.infohead.biWidth，高度为 bmpimage.infohead.biHeight*/ /*分割结果（不同区域标识对应的灰度值）要求仍然存储在bmpimage.imgbuf 为首地址的区域*/ //模糊 k 均值聚类程序unsigned int d1,d2,i,j,flag=0; unsigned char jun1,jun2,**u; float zsh0=0,zsh1=0; float J0=0,J1=0; float geshu0=0,geshu1=0,temp,temp0,temp1; u=(unsigned char**)malloc(2*sizeof(unsigned char*)); u[0]=(unsigned char*)malloc (bmpimage.infohead.biWidth *bmpimage.infohead.biHeight); u[1]=(unsigned char*)malloc (bmpimage.infohead.biWidth *bmpimage.infohead.biHeight); srand((unsigned)time(NULL));//随机初始化隶属矩阵 for(i = 0; i bmpimage.infohead.biHeight*bmpimage.infohead.biWidth; i++) { u[0][i] = (double)rand()/RAND_MAX; u[1][i] = 1 - u[0][i]; } while(flag==0) { zsh0=0; zsh1=0; geshu0=0;geshu1=0; for(i=0;ibmpimage.infohead.biHeight*bmpimage.infohead.biWidth;i++) { temp0=u[0][i]*u[0][i]; temp1=u[1][i]*u[1][i]; zsh0+=temp0; zsh1+=temp1; geshu0+=temp0*bmpimage.imgbuf[i]; geshu1+=temp1*bmpimage.imgbuf[i]; } jun1=zsh0/geshu0; jun2=zsh1/geshu1; for(i=0;ibmpimage.infohead.biHeight*bmpimage.infohead.biWidth;i++) { if(bmpimage.imgbuf[i]jun1) d1=bmpimage.imgbuf[i]-jun1; else d1=jun1-bmpimage.imgbuf[i]; if(bmpimage.imgbuf[i]jun2) d2=bmpimage.imgbuf[i]-jun2; else d2=jun2-bmpimage.imgbuf[i]; temp=(d1*1.0/d2); u[0][i]=1.0/(1+temp*temp); u[1][i]=1-u[0][i]; } J1=0; for(i=0;ibmpimage.infohead.biHeight*bmpimage.infohead.biWidth;i++) { J1+=u[0][i]*(bmpimage.imgbuf[i]-jun1)*(bmpimage.imgbuf[i]-jun1); J1+=u[1][i]*(bmpimage.imgbuf[i]-jun2)*(bmpimage.imgbuf[i]-jun2); } if(abs(J0-J1)0.01) flag=1; else J0=J1; } for(i=0;ibmpimage.infohead.biHeight*bmpimage.infohead.biWidth;i++) { if (u[0][i]=0.5) bmpimage.imgbuf[i]=255; else bmp