How to write own standard deviation function.

2 次查看（过去 30 天）
anu
anu 2016-7-20
回答： Image Analyst 2016-7-20
I have read the color image. Then separated RGB values into three different arrays. After that I have written my own function to calculate standard deviation function for each color component. But when I execute my own written function and in built function I got different values? What is wrong int it?
without inbuilt function
im = imread('D:\im112.jpg');
R=im(:,:,1)
[r,c]=size(R);
totmean=sum(R(:))/(r*c);
totdiff=(R-totmean).^2;
totsum=sum(totdiff(:));
nele=(r*c)-1;
totvar=totsum/nele;
totstd=sqrt(totvar);
display(totstd);
Using inbuilt functio
stdr=std(double(R(:)))
  1 个评论
Adam
Adam 2016-7-20
编辑：Adam 2016-7-20
Is your image in integer format? You convert to double to call the builtin, but seemingly not for your own code. I don't know off the top of my head what the result of all those operations is on integers.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Guillaume
Guillaume 2016-7-20
Adam hit the nail on the head in his comment. You're reading a jpg image so most likely the im array and thus R is of type uint8 (8 bit integer). The range for uint8 value is 0 to 255. If any operation involving uint8 overflows the result is clamped to 255 and similarly underflows are clamped to 0.
Therefore, it's highly likely that the result of sum(R(:)) is clamped to 255. Then you do a division. Because the result is still going to be uint8 the result is going to be rounded down to the nearest integer (most likely 0, since 255 divided by the image area is probably less than 1).
To solve this, you indeed need to do the same as you've done for std. First convert your R to double, then all operations will be performed correctly.

更多回答（2 个）

Andrei Bobrov
Andrei Bobrov 2016-7-20
n = numel(R);
yourstd = sqrt(sum((R(:) - sum(R(:))/n).^2)/(n - 1));
Image Analyst
Image Analyst 2016-7-20
Why not simply use std2() - the built in function meant for this????
img = imread('moon.tif');
s = std2(img) % No casting to double needed.

类别

在 Help Center 和 File Exchange 中查找有关 Convert Image Type 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by