Determine the intensity value, "T", in a 2D image for which 99.9% of all intensity values are less than "T"

Question

Ed Principe 2020-10-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/630984-determine-the-intensity-value-t-in-a-2d-image-for-which-99-9-of-all-intensity-values-are-less-t

评论： Ed Principe 2020-10-30

采纳的回答： Akira Agata

% I can read in an image and calculate a threshold for which the pixels are above some value, i.e., 6000.

Threshold=Q>6000;

% I can count the number of those pixels above that value

PixelCount = sum(Threshold(:))

% ratio the PixelCount over the total # of pixels in an image with dimensions X1, Y1 to determine a fraction, or percent of pixels above that value

Fraction=PixelCount/(X1*Y1)

But if I wish to calculate the intensity value that represents a desired fraction, i.e. 0.001, - instead of choosing an arbitrary number, like 6000........... how do I program to extract the intensity value that represents cutoff for which 99.9% of pixels have a lower intensity ?

I should mention these are uint16 images.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Ed Principe 2020-10-30

I think I have something that works for me..... Q is an image from an image stack extracted in a for loop wih index 'p':

Q=l.images1(:,:,p);

SortedQ=sort(Q(:));

figure;plot(SortedQ);

title(['Sorted Pixel Intensities:',F,': image #: ', num2str(p),'/',num2str(l.noimages)],'Interpreter', 'none');

SortedQ(round(0.999*X1*Y1))

请先登录，再进行评论。

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

Answer 1

Akira Agata 2020-10-30

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/630984-determine-the-intensity-value-t-in-a-2d-image-for-which-99-9-of-all-intensity-values-are-less-t#answer_528624

在 MATLAB Online 中打开

I believe prctile function will be helpful to this task, like:

% Read sample gray scale image
I = imread('cameraman.tif');
% Calculate threshold value (99.9% cut-off)
th = prctile(I(:),99.9);

Check the result:

>> nnz(I<th)/numel(I)

ans =

0.9990

3 个评论
显示 1更早的评论隐藏 1更早的评论

Akira Agata 2020-10-30

在 MATLAB Online 中打开

OK, then how about the following workaround?

Isort = sort(I(:));
pt = round(numel(x)*0.999);
th = Isort(pt);

If the numel(I) (= total number of pixel) is sufficiently large, the result will be the same.

Ed Principe 2020-10-30

Yes, that works, that is essentially what I came up with in my comment on the original post above. Thank you very much!!

请先登录，再进行评论。

Determine the intensity value, "T", in a 2D image for which 99.9% of all intensity values are less than "T"

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

采纳的回答

3 个评论
显示 1更早的评论隐藏 1更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

Determine the intensity value, "T", in a 2D image for which 99.9% of all intensity values are less than "T"

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

采纳的回答

3 个评论 显示 1更早的评论隐藏 1更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

3 个评论
显示 1更早的评论隐藏 1更早的评论