使用matlab对图象计算均方误差

function rms=rms(base,result) % % 计算均方根误差:含噪图像每个象素点与原始图像间误差的平方和的均值的平方根 % %计算结果，rms等于去噪后的噪声标准差 % % Calculate RMS error between result and original % % Idiff = imsubtract(base,result); % get error values % % Idbl = double(Idiff); % convert to double for arithmetic % d_idbl=block_im(Idbl,8); % for i=1:8 % for j=1:8 % Isqrd{i,j} = d_idbl{i,j}.^2; % square each pixel % end % end % % u = mean(mean(imout)); % get mean of squared errors % for i=1:8 % for j=1:8 % u(i,j) = mean(mean(Isqrd{i,j})); % get mean of squared errors % end % end % mu=mean(u(:)) % rms = mu^0.5; % root of mean of squared errors b_base=block_im(base,8); b_result=block_im(result,8); for i=1:8 for j=1:8 Isqrd{i,j} = (im2double(b_base{i,j})-im2double(b_result{i,j})).^2; % square each pixel end end for i=1:8 for j=1:8 u(i,j) = mean(mean(Isqrd{i,j})); % get mean of squared errors end end mu=mean(u(:)); rms = mu^0.5; % root of mean of squared errors