How can I improve the accuracy of identifying local maxima

Question

Jeremy Zhu 2020-7-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/563123-how-can-i-improve-the-accuracy-of-identifying-local-maxima

编辑： Image Analyst 2020-7-11

Hello all,

I used imregionalmax() to identify the position of the nucleus in the embryo image, but the result of the recognition ignores some cells at the edge of the embryo, how can I improve its accuracy?

What confuses me is that after Gaussian noise reduction, using imregionalmax() can not identify the local maximum of the image.

The following image is a comparison.

raw image

raw image after imregionalmax()

Gaussian noise reduction image

Gaussian after imregionalmax()

Thanks in advance!

2 个评论
显示无隐藏无

Matt J 2020-7-11

What confuses me is that after Gaussian noise reduction, using imregionalmax() can not identify the local maximum of the image.

My guess is that you are displaying the image over a gray level range that hides the regional maxima. The output of imregionalmax that I see is far from blank (see below). It is not useful, however, because unless you get rid of absolutely all noise, every single noise spike will show up as a regional max.

Jeremy Zhu 2020-7-11

So is there any methods to get rid of absolutely all noise?

请先登录，再进行评论。

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

Answer 1

Matt J 2020-7-11

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/563123-how-can-i-improve-the-accuracy-of-identifying-local-maxima#answer_464264

编辑：Matt J 2020-7-11

I don't think you have any hope of perfect recognition, but this result looks better to my non-clinical eye.

A=rgb2gray(imread('raw image.png'));
B=medfilt2(A,[5,5]);
NucleiiBW=bwareafilt(B-medfilt2(B,[20,20])>40,[20,inf]);
imshow(NucleiiBW)

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

Matt J 2020-7-11

编辑：Matt J 2020-7-11

Even if I were open to that, I doubt my Inbox could handle an 80 MB attachment. I don't see what there is to gain by processing the original picture. You should resize it.

Jeremy Zhu 2020-7-11

I see,thanks again!

请先登录，再进行评论。

Answer 2

Image Analyst 2020-7-11

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/563123-how-can-i-improve-the-accuracy-of-identifying-local-maxima#answer_464267

编辑：Image Analyst 2020-7-11

testGrayImage.m

imregionalmax() only identifies the local max -- it will not get the whole spot. Within the whole spot, there may be only 1 or 2 pixels that are actually the max.

I think you were actually meaning to run a top hat filter on it using imtophat(), so give that a try with different window sizes. A top hat filter is a morphological closing operation (a local max but not THE local max like imregionalmax) minus the original image. So it basically finds bright spots on a background which may be varying in intensity.

The image looks fairly flat to begin with so you can use some kind of noise reduction, like imgaussfilt() like you did or medfilt2() like Matt did, and then just threshold and so some size filtering with bwareafilt() or bwareaopen() if you still have some noise left.