为什么我无法从原始图像

Question

我有以下代码：

close all;
star = imread('/Users/name/Desktop/folder/pics/OnTheBeach.png');

blrtype = fspecial('average',[3 3]);
blurred = imfilter(star, blrtype);

[rows,cols,planes] = size(star);
R = star(:,:,1); G = star(:,:,2); B = star(:,:,3);
starS = zeros(rows,cols);
ind = find(R > 190 & R < 240 & G > 100 & G < 170 & B > 20 & B < 160);
starS(ind) = 1;

K = imfill(starS,'holes');

stats = regionprops(logical(K), 'Area', 'Solidity');
ind = ([stats.Area] > 250 & [stats.Solidity] > 0.1);
L = bwlabel(K);
result = ismember(L,find(ind));

到目前为止，我加载图像，模糊以过滤掉一些噪点，进行颜色分割以找到落在该范围内的特定对象，然后为对象创建值为1的二进制图像＆＃ 39; s颜色，0用于所有其他东西。最后，我进行了区域过滤以消除图像中留下的任何杂乱，因此我只留下了我正在寻找的对象。

现在我想根据分割蒙版重新着色原始图像，以改变海星的颜色。我想创建红色，绿色，蓝色通道，为它们赋值，然后在图像上放置遮罩。 （例如，有红色的海星）

red = star;
red(starS) = starS(:,:,255);
green = star;
green(starS) = starS(:,:,0);
blue = star;
blue(starS) = star(:,:,0);
out = cat(3, red, green, blue);
imshow(out);

这给了我一个错误：索引超出矩阵维度。

Project4中的错误（第28行） red（starS）= starS（：，：，255）;

我目前的做法有什么问题？

Answer 1

你的代码有点令人困惑......我不明白你想要使用的面具是starS还是result，因为它们看起来都像是2d索引器。在您使用starS的第二个代码段中，您在问题中发布的掩码为result。

无论如何，无论你想要的面具是什么，你所要做的就是使用imoverlay function。以下是基于您的代码的小示例：

out = imoverlay(star,result,[1 0 0]);
imshow(out);

这是输出：

Answer 2

如果Tommaso建议的imoverlay的不透明遮罩不是您所追求的，则可以修改输入的RGB值，以在所选像素上投射色调而不会使其饱和。它只涉及更多。

I = find(result);

为您提供2D图像中像素的索引。但是，star是3D。这些索引将指向相同的像素，但仅指向第一个2D切片。也就是说，如果I指向像素(x,y)，则它等效地指向像素(x,y,1)。那是像素的红色成分。要对(x,y,2)和(x,y,2)（绿色和蓝色组件）进行索引，您需要将I增加numel(result)和2*numel(result)。也就是说，star(I)访问所选像素的红色成分，star(I+numel(result))访问绿色成分，star(I+2*numel(result))访问蓝色成分。

既然我们可以访问这些值，我们如何修改它们的颜色？

这是imoverlay的作用：

I = find(result);
out = star;
out(I) = 255; % red channel
I = I + numel(result);
out(I) = 0; % green channel
I = I + numel(result);
out(I) = 0; % blue channel

相反，您可以按比例增加红色的亮度，并减少绿色和蓝色。这将改变色调，增加饱和度，并保持恒星内强度的变化。我建议使用伽马函数，因为它不会导致强烈的饱和度伪影：

I = find(result);
out = double(star)/255;
out(I) = out(I).^0.5; % red channel
I = I + numel(result);
out(I) = out(I).^1.5; % green channel
I = I + numel(result);
out(I) = out(I).^1.5; % blue channel
imshow(out)

通过增加1.5并减少0.5，您可以使效果更强。