How to display a 2D square as a 3D cube

Question

Elysi Cochin 2022-3-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1675244-how-to-display-a-2d-square-as-a-3d-cube

编辑： Elysi Cochin 2022-3-23

采纳的回答： Abolfazl Chaman Motlagh

IMG-0004-00079.zip

How can i display the attached 2D square input as a 3D cube from?

I used the same image in all slice. But the display is not clear.

Please can someone help me with displaying

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Abolfazl Chaman Motlagh 2022-3-19

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1675244-how-to-display-a-2d-square-as-a-3d-cube#answer_921519

编辑：Abolfazl Chaman Motlagh 2022-3-19

Use volshow function :

thisSlice = dicomread('IMG-0004-00079.dcm');
thisSlice = double(thisSlice);
thisSlice = (thisSlice-min(thisSlice(:)))/(max(thisSlice(:))-min(thisSlice(:)));
thisSlice = uint8(255*thisSlice);
volSlice = repmat(thisSlice,[1,1,100]);
queryPoints = 0:255;
alpha = [0 0.01 0.02 .05 0.9];
intensity = linspace(0,255,5);
alphamap = interp1(intensity,alpha,queryPoints)';
volshow(volSlice,'Alphamap',alphamap)

change Alphamap and Colormap option of volshow to have better representation for your application. there is a simple example in documentation of dicomreadvolume in MATLAB documentation.

the image in final step :

3D volumetric plot :

4 个评论
显示 2更早的评论隐藏 2更早的评论

Elysi Cochin 2022-3-20

编辑：Elysi Cochin 2022-3-20

Sir, is it possible to view the 3D image as the same texture of 2D image (that is as the image you have displayed under the image in final step)

Is it possible to do so?

Abolfazl Chaman Motlagh 2022-3-20

as i said earlier, you can change alphamap and colormap to see what is best illustration for your work.

for example if you don't use any option even dark pixels have intensity in volumentric figure:

volshow(volSlice)

but for example a little more precise number would be this: (i see that most pixel in box have intensity grater than 130 in 8bit image so i did this:)

queryPoints = 0:255;
alpha = [0 0.01 0.02 .6 1 1];
intensity = [0 50 100 130 200 255];
alphamap = interp1(intensity,alpha,queryPoints)';
volshow(volSlice,'Alphamap',alphamap)