Circle detection in CTscan

2014 4 7
2 个回答
回答已采纳
更新时间：2015 4 16
2 次查看（30 天）

0 个投票

Dear specialists, Currently, I have to detect some circles from CT scan images. The sample is as below
The radius range is from 5-110, so I tried to divide this range into several smaller intervals. However, the results is not adequate. Could anybody please help me to detect them ? Thank you very much.

2 个评论
显示 无 隐藏 无

Peter To
Peter To 2014-4-7
Sorry, I forget to mention that I use imfindcircles command. This command advices rmax <3*rmin and rmax-rmin <100.
Spandan Tiwari
Spandan Tiwari 2014-4-11
Did you try changing the Sensitivity parameter in imfindcircles? Can you show the code you used?

请先登录，再进行评论。

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

 采纳的回答

Image Analyst
Image Analyst 2014-4-7

0 个投票

Did you try imfindcircles()?

11 个评论
显示 9更早的评论 隐藏 9更早的评论

Peter To
Peter To 2014-4-8
As I mentioned in my comment, I used imfindcircles, but the results wasn't good. I have tried to play with parameters, but no improvement gained. I was surprised because I thought those circles should be easy to be detected.
Image Analyst
Image Analyst 2014-4-8
Sorry - I must not have seen your comment. If you can't get imfindcircles to work, then you can draw them manually with imellipse - let me know if you need a demo.
Or else you should try to do a background correction to brighten the center and dim the outside. Take the centroid and then take the average intensity profile radially from the centroid. Fit a quadratic to it and use it to create a model background then divide the image by it.
Then use a noise reduction filter, maybe start simple, like with medfilt2(). If that doesn't work, then try more sophisticated ones like mean shift, non-local means, or things like that.
Then try to use an edge detector like imgradient or a texture filter like entropyfilt() or stdfilt(). It will definitely take some experimentation.
Peter To
Peter To 2014-4-8
@Image Analyst: Thank you very much for your reply. Because I have thousands of images like this to analyse, I think I prefer an automatic solution with imfindcircles or something similar. I also had tried to correct the contrast, but the noise became more clear that reduced the accuracy of imfindcircles. Therefore, I tried to reduce the noise first. I haven't used non-local means but I used medfilt2 and wiener2. However, both of them do not increase the accuracy of detection. Maybe I did not use the right parameters, but I don't know how because I have tried several values in vain. As I mentioned above, I tried to divide the radius range into several small interval with different sensitivity. However, the big circles can be detected with only very high sensitivity 0.98, which makes some fake circles. It surprised me because I thought the big circles are very easy to be detected.
Image Analyst
Image Analyst 2014-4-9
Can you average several slices together to reduce noise, then use the same circles for each slice?
Peter To
Peter To 2014-4-9
编辑：Peter To 2014-4-9
@ Image Analyst: Yes, I did it. I tried to read 5 slices at once, and correlated the images. When a point had its grayscale values more than 150 in 4 slices, I set the value to 255. However, it did not improved the results. Do you have any suggestion ?
Image Analyst
Image Analyst 2014-4-9
Have you tried averaging them together? Maybe try texture segmentation. Check out standard ISO 25178 on Surface Texture.
Peter To
Peter To 2014-4-11
@Image Analyst: The radius of circles should change by slices, so the averaging process can make some extra difficulties for the circle detection. I have no idea about the texture segmentation. Could you please tell me more precisely?
Image Analyst
Image Analyst 2014-4-11
Try entropyfilt() and stdfilt() and see what you get.
Peter To
Peter To 2014-4-12
I tried, and both of them give me ugly results like
Image Analyst
Image Analyst 2014-4-12
The standard deviation filter did a pretty good job on the bigger circles. See attached m-file.
Peter To
Peter To 2014-4-14
It is awesome, I will study your code right now. Although this code cannot detect the small circles, it helps me to start over. Thank you very much.

请先登录，再进行评论。

更多回答（1 个）

Julia
Julia 2015-4-15

0 个投票

This code looks really useful. I tried to use it for the detection of pores in the attached picture. As you can see, especially in the center, the pore spaces are widely connected, so a distinction would be needed. Basically i would like to fit in circles into the pore spaces to get a pore size distribution. Using your proposed code as it is does not help in this case. Could you suggest which parameters to play with in order to make it fit better for less obvious cases like mine?

2 个评论
显示 无 隐藏 无

Image Analyst
Image Analyst 2015-4-15
There is no picture attached. I suggest you start a new question with it. You can reference this one if you want.
Julia
Julia 2015-4-16
编辑：Julia 2015-4-16
I attached the picture. I will also start a new question. Thanks. http://www.mathworks.com/matlabcentral/answers/203771-detect-highly-connected-pores-in-a-ct-scan

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Image Filtering 的更多信息

产品

标签

Community Treasure Hunt

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

Start Hunting!

Translated by