Matlab 3D 体积可视化 - dicom 文件

标签 matlab 3d dicom

我期待将 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])

感谢您的帮助。

最佳答案

解决问题的两种选择:

  1. 天真的,为什么不直接设置 v=v+min(v(:)) 以便图像从零缩放到不同的最大值?

  2. 为什么带有负 isovlaue 的 isosurface(V,isovalue) 不能为您解决这个问题? 或者,如果您希望从 (-n:step_size:m) 循环遍历等值以获得所需的动态范围(使用一些 alpha(0.2) 来查看所需的图层)

    <

关于Matlab 3D 体积可视化 - dicom 文件,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/13084377/

相关文章:

image - 如何在 MATLAB 中从视频中选择特定帧?

matlab - 如何创建带有轮廓边框的矩形?

c++ - 给定 3D 旋转矩阵和 3D 方向 vector 求角度差

dicom - 显示相同的 Dicom 图像标题

image-processing - 如何以可靠的方式获取 DICOM 像素数据的颜色深度?

matlab - 计算单元格数组matlab中的单词

matlab - 如何在 MATLAB 中将 theta 限制设置在 -90 到 90 范围内?

3d - COLLADA 文件查看器

c++ - 具有连续结果的夹角

dicom - 在覆盖 DICOM 中的像素数据时,应替换除 UID 之外的哪些 DICOM 标签?