Signal to noise ratio color image (code) for dB

1 次查看(过去 30 天)
I am trying to compute signal to noise ratio (dB) in color images as follows:
noisydata = imnoise(original,'gaussian',0.0,0.1); % mean = 0.0, variance = 0.1
[var_orgR, var_noiseR] = snr(noisydata(:,:,1),original(:,:,1));
[var_orgR, var_noiseR] = snr(noisydata(:,:,2),original(:,:,2));
[var_orgR, var_noiseR] = snr(noisydata(:,:,3),original(:,:,3));
SNR= 10*log10((var_orgR + var_orgG + var_orgB)/(var_noiseR + var_noiseG + var_noiseB));
function [var_original, var_noise] = snr(noisydata, original)
mean_original = mean(original(:));
tmp = original - mean_original;
var_original = sum(tmp(:).^2);
noise = noisydata - original;
mean_noise = mean(noise(:));
tmp = noise - mean_noise;
var_noise = sum(tmp(:).^2);
end
The problem is I get value of SNR = 0.0330, I would like to get SNR in dB like 1 dB, 5 dB, etc.
Could any one please guide/suggest how to obtain SNR in dB?
Thank you. Dbs

回答(2 个)

Image Analyst
Image Analyst 2013-2-3
Where did you get these formulas? I'm not sure I agree with them. In fact I don't, unless you can show me some reference that say otherwise.
  5 个评论
Image Analyst
Image Analyst 2013-2-3
I don't think that would fix it. I would need proof. Anyway, I don't know what a "fix" is. Like I said, you asked how to get the number in decibels and I told you that it's already in decibels.
Dadhichi
Dadhichi 2013-2-4
What proof do you need? If you have access to papers you'll agree. The code is correct, uint8 was programming error, it is not required. Anyways, 'noisydata' and 'original' should have values in same range i.e. between 0 to 1 (was another programming error) :D .

请先登录,再进行评论。


Toqeer Mahmood
Toqeer Mahmood 2015-10-9
Dear Dadhichi
I am also confused in the same problem, I also want to find SNR with the whole number as 5dB, 10db, 15dB etc. Please send me the working code to find it, if possible for you (my email is toqeer.mahmood@yahoo.com).
Thanks in advance, waiting for your kind reply.

类别

Help CenterFile Exchange 中查找有关 Spectral Measurements 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by