developing clusters in a 3D image

16 次查看（过去 30 天）
C.
C. 2014-10-9
评论： Matt J 2014-10-13
I am very new to MATLAB (and programming) and I'm probably out of my skill level in attempting this, but...
I need to write a script that can take a 3D image and given an x,y,z coordinate within that image test all adjacent coordinates to the given coordinate to see if any of them are greater than 0. If an adjacent coordinate is greater than zero I need to repeat the aforementioned process for this coordinate as well until all coordinates > 0 and branching from the original coordinate are exhausted. Could anyone perhaps give me any hints as to a basic method (not necessarily efficient) to accomplish this? I suspect this process has to be recursive?
Thanks
  1 个评论
Image Analyst
Image Analyst 2014-10-9
What do you want as an output after you're done doing the region growing (which is what it's called, not clustering). And it can be recursive but can be done in a single line of code by calling bwconncomp() or bwlabel().

请先登录，再进行评论。

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

回答（2 个）

Matt J
Matt J 2014-10-9
编辑：Matt J 2014-10-9
You could do this pretty easily with bwconncomp()
BW=Image>0;
BW(x,y,z)=true;
CC=bwconncomp(BW);
lidx=sub2ind(size(BW),x,y,z);
Now just loop through all CC(i).PixelIdxList and see which region contains lidx.
  2 个评论
C.
C. 2014-10-13
编辑：C. 2014-10-13
Basically I am trying to find ROIs in an fMRI image using previous data about where these ROIs are (The coordinate is the max intensity voxel within an activation cluster). The script will start at this voxel and find all active neighbors until they are exhausted and produce a binary image with 1s in active voxels and 0s elsewhere. I'm a little confused by what the last line is doing, however.
lidx=sub2ind(size(BW),x,y,z);
Matt J
Matt J 2014-10-13
It's converting the x,y,z coordinates to linear indeces, which is the format that the pixel coordinates are held in in CC(i).PixelIdxList.
See,
http://www.mathworks.com/help/matlab/ref/sub2ind.html

请先登录，再进行评论。

Image Analyst
Image Analyst 2014-10-9
编辑：Image Analyst 2014-10-9
In addition to the way Matt said, you can use bwlabel and ismember
labeledImage = bwlabel(yourArray);
regionIWant = labeledImage(row, col, slice);
outputImage = ismember(labeledImage, regionIWant);
Or you could use imreconstruct to do a morphological reconstruction:
seedImage = false(size(yourArray));
seedImage(row, col, slice) = true;
outputImage = imreconstruct(seedImage, yourArray);

类别

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

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by