Volumetric slice plot of tomographic data

2017 4 20
2 个回答
回答已采纳
更新时间：2017 4 30
5 次查看（30 天）

0 个投票

I am trying to use the slice function to stack the following 79 slices and display a 3D representation of an object.
The problem is that the black regions in each image are obstructing the image of the previous slice. So, I'm only seeing the very last slice on top, rather than a meaningful volumetric presentation of the data. So, how can I include only the non-zero elements of each image in the plot?
Any suggestions would be greatly appreciated.

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

Philip G
Philip G 2017-4-20
编辑：Philip G 2017-4-20
Exactly. As you can see from the documentation of "slice" - it just makes solid planes. So either you reduce the number of planes to display (so not your 1:79 but just a few slices in different orientations).
A second option would be to use one of the many 3D viewer on the FEX. here are some.
If you really would like to get a 3d volume of your data ... that is a little bit more tricky. You could take a look at the MATLAB functions isosurface (that basically relies on a simple threshold that creates a surface) and related helper functions like smooth3 to clean up your data.
Hope this already helps a bit.
Edit: vol3d might be very useful for you. As an alpha map you can basically use a binary mask or even your inverted brightness level.
Sordin
Sordin 2017-4-21
Hi Philip, thank you for your suggestions.
I already have the inverted brightness or complement images for each of the slices (I can also make the binary images). But how can I apply an alpha mask to images based on these masks? Do you know what code I need to use?
Basically, I want the black areas (backgrounds) to be made highly transparent. While the non-zero elements of each image (containing data) are to be opaque. I think this should help in making the 3D reconstruction more recognizable in terms of the anatomical features of the patient.
This is my code (I am trying to write it myself before trying vol3d or other viewers):
for i = 1 : 39
I=imread(filename);
I(:,:,i) = imread(filename);
mask = imcomplement(I(:,:,i));
M(:,:,i) = I(:,:,i);
end
hs = slice(M,[],[],1:1:39, 'nearest');
colormap('gray');
shading interp
set(hs,'FaceAlpha',0.8);
Philip G
Philip G 2017-4-24
I found a nice way to use surf images ... see my anwer.

请先登录，再进行评论。

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

 采纳的回答

Philip G
Philip G 2017-4-24
编辑：Philip G 2017-4-24

2 个投票

Hi, if you do not want to use a simple isosurface or volume viewer as I suggested in my comment, every image can have an alphadata property (that can be a matrix and not only a single value), as explained here:
https://www.mathworks.com/company/newsletters/articles/image-overlay-using-transparency.html
You can extend this behavior towards 3D with the help of this:
https://de.mathworks.com/help/matlab/visualize/add-transparency-to-graphics-objects.html
For example by using surf.
One example:
img=peaks; % this is the color
t=abs(img); % this is the transparency
%t=img>1; % if you want to have a binary mask
for i=1:10
surf(img*0+i,p,'linestyle','none','alphadata',t,'facealpha','flat')
hold on
end
I hope you get the idea.

更多回答（1 个）

Bjorn Gustavsson
Bjorn Gustavsson 2017-4-21

1 个投票

If you change your colour-map from black-white to white-black then you might get some mileage out of experimentations with alpha-mapping, look at the help for alphamap and alim.
HTH

类别

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

产品

标签

Community Treasure Hunt

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

Start Hunting!

Translated by