求图像的信噪比PSNR

　　图像的信噪比是衡量一幅图片质量好坏的重要指标。图像的信噪比应该等于信号与噪声的功率谱之比，但功率谱不好计算，可以用信号与噪声的方差之比来估算。首先计算所有像素的局部方差，将局部方差的最大值认为是信号的方差，最小值是噪声的方差，求出比值，再转换成dB数，最后用经验公式修正。

　　均方差定义为：

　　

　　峰值信噪比定义为：

　　

　　其中MAXI是表示图像点颜色的最大数值，如果每个采样点用8位表示，那么就是255。更为通用的表示是，如果每个采样点用B位线性脉冲编码调制表示，那么MAXI就是

　　　　　　　　　　2^B-1

　　对于每点有RGB三个值的彩色图像来说，峰值信噪比的定义类似，只是均方差所有方差之和除以图像尺寸再除以3。

　　图像压缩中典型的信噪比在20～40dB之间，越高越好。

　　

#include <cv.h>
#include <highgui.h>
#include <math.h>
#include <stdio.h>

int main(int argc, char* argv[])
{
    if(argc<2)
    {
        printf("Usage: main <image-file-name>/n/7");
        exit(0);
    }
    IplImage* img = cvLoadImage(argv[1],-1);
    int height=img->height;
    int width=img->width;
    int step=img->widthStep;
    uchar *data=(uchar *)img->imageData;
    int i,j;
    double sum=0;
    for(i=0;i<height;i++) 
    {
        for(j=0;j<width;j++) 
        {
            sum+=data[i*step+j];
        }
    }
    double mean=0;
    mean=sum/(width*height);
    double mse=0;
    for(i=0;i<height;i++) 
    {
        for(j=0;j<width;j++) 
        {
            mse+=(data[i*step+j]-mean)*(data[i*step+j]-mean);
        }
    }
    mse=mse/(width*height);
    double psnr=0;
    psnr=10*log10(255*255/mse);
    printf("sum  is: %lf
",sum);
    printf("mean is: %lf
",mean);
    printf("mse  is: %lf
",mse);
    printf("psnr is: %lf
",psnr);
    cvNamedWindow(argv[1], CV_WINDOW_AUTOSIZE); 
    cvShowImage(argv[1], img );    
    cvWaitKey(0);
    cvDestroyWindow(argv[1]);
    cvReleaseImage(&img );
    return 0;
}