我给100分 急...........那位大哥大姐帮我用matlab翻译下
下面的VC 是进行直方图均衡化和规定化操作
那位大哥大姐帮我用matlab翻译下 谢谢~~~
翻译好 能发给我么?
E_mail: ctfysj@126.com
QQ:308463776
void main()
{
BMPFILE bmpfile;
bmpfile.LoadBMPFILE ("实验图象1.bmp");
if(!bmpfile.bIsOk)
{
printf("Error loading image.\n");
return;
}
// do other processing with the imagedata
int H[256]; //原始图象直方图
memset(H,0,sizeof(H));
int S[256]; //变换映射-均衡直方图S
double psnr; //峰值信噪比
int i;
int f;
for(i=0;i<bmpfile.imageh*bmpfile.imagew;i++)
H[*(bmpfile.imagedata+i)]++;
S[0]=H[0];
for(i=1;i<256;i++)
S=S[i-1]+H;
for(i=0;i<256;i++)
S=S*255/(bmpfile.imageh*bmpfile.imagew);
for(i=0;i<bmpfile.imageh*bmpfile.imagew;i++)
*(bmpfile.imagedata+i)=S[*(bmpfile.imagedata+i)];
//原图象的直方图H
printf("H: ");
for(i=1;i<256;i++)
printf("%d ",H);
printf("\n");
//原图象(直方图H)的均衡直方图
printf("S: ");
for(i=1;i<256;i++)
printf("%d ",S);
printf("\n");
//均衡图象
bmpfile.SaveBMPFILE ("实验图象1-均衡图象.bmp");
//PSNR值
psnr=0;
for(i=0;i<255;i++)
psnr+=(S-i)*(S-i)*H;
psnr=10*log10(double(bmpfile.imageh)*bmpfile.imagew*255*255/psnr);
printf("PSNR: %f\n",psnr);
//目标直方图G的均衡直方图
int G[256]; //规定直方图
int GS[256]; //规定直方图的均衡直方图
for(i=0;i<=127;i++)
G=i;
for(i=128;i<=255;i++)
G=(255-i);
GS[0]=G[0];
for(i=1;i<256;i++)
GS=GS[i-1]+G;
f=GS[255]/255+1;
for(i=0;i<256;i++)
GS/=f;
printf("GS: ");
for(i=1;i<256;i++)
printf("%d ",GS);
printf("\n");
//最终增强图象
int i2;
int minNumber;
int SG[256]; //从原图象的均衡图象到最终增强图象的映射
memset(SG,0,sizeof(SG));
for(i=0;i<256;i++)
{
minNumber=255;
for(i2=0;i2<256;i2++)
if(abs(S-GS[i2])<minNumber)
{
minNumber=abs(S-GS[i2]);
SG[S]=i2;
}
}
for(i=0;i<bmpfile.imageh*bmpfile.imagew;i++)
*(bmpfile.imagedata+i)=SG[*(bmpfile.imagedata+i)];
bmpfile.SaveBMPFILE ("实验图象1-最终增强图象.bmp");
//最终增强图象的直方图GR R=real
int GR[256];
memset(GR,0,sizeof(GR));
for(i=1;i<256;i++)
GR[SG[S]]+=H;
printf("GR: ");
for(i=1;i<256;i++)
printf("%d ",GR);
printf("\n");
//最终增强图象的PSNR值
psnr=0;
for(i=0;i<255;i++)
psnr+=(SG[S]-i)*(SG[S]-i)*GR;
psnr=10*log10(double(bmpfile.imageh)*bmpfile.imagew*255*255/psnr);
printf("PSNR of the destination image: %f\n",psnr);
[解决办法]
http://ctfysj.blog.hexun.com/8726749_d.html
直方图均衡化和规定化都有,你自己看吧
