图像信噪比和峰值信噪比计算

matlab实现计算图像信噪比和峰值信噪比

A = imread（‘3.bmp‘）;%读入图片
A = double（A）;  %把矩阵A定义为double型
B = imread（‘333.bmp‘）;%同上
B = double（B）;

%% 计算MSE

[mn] = size（A）;%计算A的像素（长和宽）
[m1n1] = size（B）;
if m~=m1||n~=n1     %如果A B长宽不同就显示“图像大小不一致”
    error（‘图像大小不一致‘）;
end

msevalue = 0;  %定义初始值
for i = 1:m
    for j = 1:n    %循环用于求和
        msevalue = msevalue+（A（ij）-B（ij））^2;    %每循环一次就加一次
    end
end
msevalue = msevalue/（m*n）;  %用求的和除以m*n
if msevalue == 0      %表示A（ij）=B（ij）
    error（‘图像完全相同‘）;
end

%% 计算峰值信噪比
psnrvalue = 255^2/msevalue;  %这个你懂得     第二个程序类似的
psnrvalue = 10*log10（psnrvalue）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        693  2012-05-27 17:13  信噪比和峰值信噪比计算\test1.m

     文件        735  2012-05-30 16:13  信噪比和峰值信噪比计算\test2.m

     目录          0  2012-06-06 13:13  信噪比和峰值信噪比计算

----------- ---------  ---------- -----  ----

                 1428                    3