Roipoly for 3D images

12 views (last 30 days)
Jacob
Jacob on 1 Dec 2011
Hello,
I am trying to use the roipoly tool to segment out a region from a stack of medical images. In many cases the ROI is fairly constant between different slices of the total volume. Is there a way to select the ROI in the first slice, copy those points to the adjacent slice, and then edit those points?
Right now I'm trying the following:
[slice_1_roi tmp_xi tmp_yi]=roipoly(slice1);
[slice_2_roi tmp_xi tmp_yi]=roipoly(slice2,tmp_xi,tmp_yi);
[slice_3_roi tmp_xi tmp_yi]=roipoly(slice3,tmp_xi,tmp_yi);...
However, this just copies the ROI from the first slice without allowing any interaction in the next two slices. Help?

Accepted Answer

Amith Kamath
Amith Kamath on 2 Dec 2011
This is interesting too. I had to work on a similar problem of 3D segmentation, where the correlation between slices was significant, although my contours were not manually set for each slice. But you could use the function called impoly. For example,
load mri
[slice_1_roi tmp_xi tmp_yi]=roipoly(D(:,:,:,1));
imshow(D(:,:,:,2)); %To get the figure handle.
h = impoly(gca,[tmp_xi tmp_yi]); %Alter your mask here. DONT close the figure window yet.
BW = createMask(h); %BW contains the mask which you just altered.
pos = getPosition(h); %Stores the points in your new mask.
%Close the figure now, and repeat this for the next slice.
imshow(D(:,:,3));
h = impoly(gca,pos); %Now pos contains the points from your mask.
Let me know if this works!
  2 Comments
Jacob
Jacob on 2 Dec 2011
Great solution, thank you!
Jawahir
Jawahir on 21 Mar 2013
Is this possible in 3D medical image in .obj format? i've already read the file and displayed it in the figure window in mesh form.

Sign in to comment.

More Answers (0)

Categories

Find more on Image Processing Toolbox in Help Center and File Exchange

Products

Community Treasure Hunt

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

Start Hunting!