3

我会尽量做到准确和简短。

我有一个体积（128x128x128）和一个掩码（与 [0|1|2] 值相同的大小）

我想使 3D 体积矩阵成为具有 RGB 的 3D 图像，并将掩码中标记的点存储在每个通道（红色、绿色、蓝色）中。

这是通过获取该 3D 立方体的切片来使用 2D 表示，而不是一遍又一遍地计算它以使事情变得更快（在我的项目中非常重要），所以实际上，3D 体积 + rgb 就像一个商店用于 128 张 2D 图像。

问题是，我必须采取哪些步骤以及如何完成所有这些： - 创建一个 128x128x128x3 的卷？- 定义一个新的颜色图（原来是灰色的）？- 加入每个频道？- 如何使用 imagesc/whatever 来显示该立方体的一个切片，其颜色点与掩码中标记的颜色相同（例如：imageRGB(:,:,64)）？

这只是我的猜测，但我什至不知道如何正确地做到这一点......我有点迷茫，希望你能帮助我，这是一段可能是错误的代码，但可能会帮助你

% Create the matrix 4D
ovImg = zeros(size(volImg,1),size(volImg,2),size(volImg,3),3);    % 128x128x128x3
% Store in each channel the points marked as groups
ovImg(:,:,:,1) = volImg .* (mask==1);
ovImg(:,:,:,2) = volImg .* (mask==2);
ovImg(:,:,:,3) = volImg .* (mask==3);

非常感谢！！

更新：

我在透明度和颜色图方面遇到了一些问题，这就是我所做的。

% Create the matrix 4D
ovImg = zeros(size(volImg,1),size(volImg,2),size(volImg,3),3);
% Store in each channel the points marked as groups
ovImg(:,:,:,1) = imaNorm.*(mask==1);
ovImg(:,:,:,2) = imaNorm.*(mask==2);
ovImg(:,:,:,3) = imaNorm.*(mask==3);

[X,Y,Z] = meshgrid(1:128,1:128,1:128);
imaNorm = volImg - min(volImg(:));   
maxval = max(imaNorm(:));        
ovImg  = imaNorm + mask * maxval;


N= ceil(maxval);
c = [linspace(0,1,N)' zeros(N,2)];
my_colormap = [c(:,[1 2 3]) ; c(:,[3 1 2]) ; c(:,[2 3 1])];

figure;
imshow(squeeze(ovImg(:,:,64)),my_colormap);

figure;
imagesc(squeeze(mask(:,:,64)));

结果（叠加图像/蒙版）有 叠加和蒙版图片 什么想法吗？再次感谢大家

最后更新：使用 Gunther Struyf 建议的另一种方法，我得到了我想要的。谢谢朋友，我真的很感激，希望这对其他人也有帮助。

4

1 回答 1

2

您可以使用带有颜色图的imshow从灰度图像（您拥有）中“伪造”RGB 图像。对于比例，我不会将其相乘，而是为该值添加一个偏移量，因此每个蒙版都是颜色图中的不同范围。

要绘制 3d 矩阵的切片，您可以对其进行索引，然后对其进行挤压以删除生成的单例维度：

例子：

[X,Y,Z]=meshgrid(1:128,1:128,1:128);
volImg =5*sin(X/3)+13*cos(Y/5)+8*sin(Z/10);
volImg=volImg-min(volImg(:));
mask = repmat(floor(linspace(0,3-2*eps,128))',[1 128 128]);

maxval=max(volImg(:));
ovImg=volImg+mask*maxval;
imshow(squeeze(ovImg(:,:,1)),jet(ceil(max(ovImg(:)))));

未屏蔽的原始图像 ( imshow(squeeze(volImg(:,:,1)),jet(ceil(maxval))))

在此处输入图像描述

结果与掩码（上面的代码块）：

在此处输入图像描述

对于不同的颜色图，请参见此处，或创建您自己的颜色图。例如，你的面具有三个值，所以让我们将它们与 R、G 和 B 匹配：

N = ceil(maxval);
c = [linspace(0,1,N)' zeros(N,2)];
my_colormap = [c(:,[1 2 3]) ; c(:,[3 1 2]) ; c(:,[2 3 1])];
figure
imshow(squeeze(ovImg(:,:,1)),my_colormap);

这使：

在此处输入图像描述

其他方法：

现在我理解了您的问题，我看到您从一开始就完全正确，您只需要将变量重新调整为 0 到 1 之间的值，因为来自imshow：

颜色强度可以在 0.0 到 1.0 的区间内指定。

您可以使用：

minval=min(volImg(:));
maxval=max(volImg(:));
volImg=(volImg-minval)/(maxval-minval);

接下来是您的代码：

ovImg = zeros([size(volImg),3]);
ovImg(:,:,:,1) = volImg .* (mask==1);
ovImg(:,:,:,2) = volImg .* (mask==2);
ovImg(:,:,:,3) = volImg .* (mask==3);

你现在只需要绘制它：

imshow(squeeze(ovImg(:,:,64,:)))
于 2012-08-18T07:53:36.393 回答