ISOSURFACE: identifying the connected faces

AP
2012 10 19
2 个回答
更新时间：2022 10 31
41 次查看（30 天）

0 个投票

Dear All,
I have extracted the isosurfaces of variable "A" in a 3D volume data. For a specific value of "A", say A=1, I get three connected series of faces which are separated from each other. In other words, there are three different and isolated regions in the domain where A=1.
I am trying to identify these three regions and save them into three different variables. I would be grateful if someone could give me a hint on how I can achieve this goal.
Thanks,
Ahmad

2 个评论
显示 无 隐藏 无

Sean de Wolski
Sean de Wolski 2012-10-19
What do you plan to do with these variables after?
AP
AP 2012-10-19
编辑：AP 2012-10-19
The domain that I have has noise in it. So, when I plot an isosurface of some value, I get a few little islands in the plot beside the main one. I want to have the biggest isosurface. This may be solved by smoothing but I want to avoid using it as it changes my results.

请先登录，再进行评论。

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

回答（2 个）

Matt J
Matt J 2012-10-19

0 个投票

Use BWCONNCOMP or REGIONPROPS

10 个评论
显示 8更早的评论 隐藏 8更早的评论

AP
AP 2012-10-19
Matt,
These are for images. I have 3D objects.
Matt J
Matt J 2012-10-19
编辑：Matt J 2012-10-19
Not sure why you think they won't handle n-dimensional data (unless possibly you have a really old version of MATLAB). Here's a simple example on a 3D volume
>> v=false(3,3,3); v(1)=1; v(end)=1
v(:,:,1) =
1 0 0
0 0 0
0 0 0
v(:,:,2) =
0 0 0
0 0 0
0 0 0
v(:,:,3) =
0 0 0
0 0 0
0 0 1
>> S=regionprops(v,'PixelList'); S.PixelList
ans =
1 1 1
ans =
3 3 3
AP
AP 2012-10-19
编辑：AP 2012-10-19
Thank you for clarification. I would be grateful if you could let me know how I should interpret the results. For my case, I get an array of size 1250x3. How this result can be used to differentiate between isolated objects? I can see in the plot that there are three separated objects but I don't see its connection with the result of S.PixelList. One of the isolated objects is much bigger than the other two but the results of S.PixelList is an array of size 1250x3.
Matt J
Matt J 2012-10-19
编辑：Matt J 2012-10-19
The output should not be a 1250x3 numeric array. It should be a structure array. Make sure your input to regionprops is a logical array, e.g., as in
S=regionprops( (A==1), 'PixelList');
Each S(i) will correspond to a different connected region where A=1.
AP
AP 2012-10-19
I think this is good only for images or logical arrays. It does not work well for real valued arrays. We may have different regions very close to each other which prevents regionprops in resulting an accurate result.
Matt J
Matt J 2012-10-19
编辑：Matt J 2012-10-19
But you will inevitably have a logical array. You said you were looking for regions for example where A equals 1. If so, you will have to compute the logical array A==1
More generally, you will always have to apply some logical test that identifies the regions you're looking for.
AP
AP 2012-10-22
编辑：AP 2012-10-22
I think the way isosurface function works is different from setting A==1.
Can I send you the MAT file? For A==1 I get nothing as the result of S.PixelList although I can see two separate surfaces in the plot.
Matt J
Matt J 2012-10-22
Well how do you find the isosurface? If it isn't A==1, I assume it's something like
A<=1.5 & A>=0.5
or some other logical test.
AP
AP 2012-10-22
I had tested this before by defining a tolerance but unfortunately I don't get the results of S.PixelList similar to that of isosurface function.
Can I send you the MAT file? You may find the source of discrepancy.
naomy_yf
naomy_yf 2022-10-31
regionprops is the precondition for using BWCONNCOMP

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Scalar Volume Data 的更多信息

标签

提问:

AP
AP
2012-10-19

评论：

naomy_yf
naomy_yf
2022-10-31

Community Treasure Hunt

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

Start Hunting!

Translated by