Question

在matlab里怎么求图片的信噪比？

知道这个求信噪比的公式
function snr=estsnr(x)
%SNR estimation of blurred noisy images
[n,m]=size(x);
c=ones(5,5)/25;
a=conv2(x,c);
al=a(3:n+2,3:m+2);
v=(double(x)-double(al)).^2;
a=conv2(double(v),double(c));
v=a(3:n+2,3:m+2);
b=10*log10(max(max(v))/min(min(v)));
snr=1.04*b-7;
该怎么求啊？不懂啊，需要图片的名字么，在哪里写代码什么的~
??? Strings passed to EVAL cannot contain function declarations.
直接运行会提示这个错误

Answer 1

1、在当前文件夹中新建一个SNR1.m文件，内容为：“ function snr=SNR1(I,In) Ps=sum(sum(I.^2)); Pn=sum(sum((In-I).^2)); snr=10*log10(Ps/Pn);”。

2、准备好两个矩阵a和b，这里使用的是两个291*240的double类型的矩阵。

3、在命令行窗口中输入：“x=SNR1(a,b)”，调用已经写好的信噪比函数。

4、回车后计算得a和b的信噪比为25.2021。

Answer 2

这是一个函数，你要么在另一个程序中定义X（即读入一幅图像）以后在调用它，要么在这个函数里面定义好X以后直接运行。

追问

额。。具体怎么操作啊啊 我不懂啊

Answer 3

Matlab中计算信噪比方式：
%===========================Happy===================================%
function snr=SNR(I,In)
% 计算信号噪声比函数
% by Qulei
% I :original signal
% In:noisy signal(ie. Original signal + noise signal)
% snr=10*log10(sigma2(I2)/sigma2(I2-I1))
[row,col,nchannel]=size(I);
snr=0;
if nchannel==1%gray image
Ps=sum(sum((I-mean(mean(I))).^2));%signal power
Pn=sum(sum((I-In).^2));%noise power
snr=10*log10(Ps/Pn);
elseif nchannel==3%color image
for i=1:3
Ps=sum(sum((I(:,:,i)-mean(mean(I(:,:,i)))).^2));%signal power
Pn=sum(sum((I(:,:,i)-In(:,:,i)).^2));%noise power
snr=snr+10*log10(Ps/Pn);
end
snr=snr/3;
end