我期待将 3D 矩阵从 dicom 文件可视化到 matlab 中。由于我对matlab不太熟悉,我设法从this post获得帮助。 :
不同之处在于我的矩阵不是由 1 和 0 组成,而是由负数组成,这些负数通过 imshow(dimeread(dicomFile)) 正确显示为红色
如何通过 3D 渲染获得相同的对比度?
我的代码:
dicomFilesZm = dir(fullfile(myDcmFolder, 'SLZ-*.dcm')); %Get files name
dicomFilesZp = dir(fullfile(myDcmFolder, 'SLZ+*.dcm')); %~
Z = dicomFilesZm(end:-1:1); % sort
dicomFilesZ = [Z ; dicomFilesZp]; % recompose final array with files name
Iz1 = fullfile(myDcmFolder, dicomFilesZ(1).name);
v = NaN([size(dicomread(Iz1)) numel(dicomFilesZ)]); % creation of empty matrix with the good size
for i = 1 : numel(dicomFilesZ)
Iz = fullfile(myDcmFolder, dicomFilesZ(i).name);
v(:,:,i) = dicomread(Iz); % fill the matrix with each image
end
p = patch( isosurface(v,0) );
isonormals(v, p)
set(p, 'FaceColor','r', 'EdgeColor','none')
daspect([1 1 1])
感谢您的帮助。
最佳答案
解决问题的两种选择:
天真的,为什么不直接设置
v=v+min(v(:))
以便图像从零缩放到不同的最大值?为什么带有负 isovlaue 的
<isosurface(V,isovalue)
不能为您解决这个问题? 或者,如果您希望从 (-n:step_size:m) 循环遍历等值以获得所需的动态范围(使用一些alpha(0.2)
来查看所需的图层)
关于Matlab 3D 体积可视化 - dicom 文件,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/13084377/