均匀量化（lena图）

#include<stdio.h>
#include<math.h>
#define ROW 512
#define COL 512
typedef unsigned char  BYTE;
double cal_psnr(BYTE origin[ROW][COL],BYTE output[ROW][COL])
{
    int i, j;
    double PSNR = , MSE = , MAXI = ;
    for (i = ; i < COL; i++) {
        for (j = ; j < ROW; j++) {
            MSE += (origin[i][j] - output[i][j]) * (origin[i][j] - output[i][j]);
        }
    }
    MSE = MSE/(ROW*COL);
    printf("MSE: %f\n", MSE);
    PSNR = *log10(MAXI) - *log10(MSE);
    return PSNR;
}
int main()
{
    FILE *f = NULL;
    BYTE output[ROW][COL];
    BYTE origin[ROW][COL];
    f = fopen("D:\lena512.raw","rb");
    fread(origin,sizeof(BYTE),ROW*COL,f);
    printf("has already read\n");
    for(int i =;i<ROW;i++)
        for(int j=;j<COL;j++)
        {
            if(origin[i][j]>=&&origin[i][j]<=)
                output[i][j] = ;
            else if (origin[i][j]>&&origin[i][j]<=)
                output[i][j] = ;
            else if (origin[i][j]>&&origin[i][j]<=)
                output[i][j] = ;
            else
                output[i][j] = ;
        }
    printf("PSNR: %lf",cal_psnr(origin,output));
    fclose(f);
 
    return ;
}

lena512.raw 下载地址：https://files.cnblogs.com/files/jzcbest1016/lena512_20171219131444306.rar

.raw文件可以用photoshop打开