Retrieving the lung pixels through matrix operations with images

2 次查看（过去 30 天）

Sivaramakrishnan Rajaraman 2018-8-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/414381-retrieving-the-lung-pixels-through-matrix-operations-with-images

评论： Sivaramakrishnan Rajaraman 2018-8-13

I have two images, an original CXR, and its corresponding lung mask. I would like to perform an operation with these two images, to generate only the lung ROI from the original CXR and make the rest of the background black. The resultant image should be in RGB as the original CXR and mask, and not grayscale or binary. I tried performing bitwise multiplications and X-OR operations but doesn't work. Can you suggest me the code to do this?

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

KALYAN ACHARJYA 2018-8-11

I have edited the answer, Have you check? Please confirm?

请先登录，再进行评论。

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

采纳的回答

KALYAN ACHARJYA 2018-8-10

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/414381-retrieving-the-lung-pixels-through-matrix-operations-with-images#answer_332317

编辑：KALYAN ACHARJYA 2018-8-11

在 MATLAB Online 中打开

%This can do without for loop also, recomemded to follow the without loop  
[rows colm]=size(rgb_image);
      for i=1:rows
          for j=1:colm
              if (mask(i,j)==0) %Mask should be binary Image, if not convert it.
                  rgb_image(i,j,1)=0;
                  rgb_image(i,j,2)=0;
                  rgb_image(i,j,3)=0;
              end
          end 
      end 
      figure,imshow(rgb_image);

9 个评论
显示 7更早的评论隐藏 7更早的评论

Image Analyst 2018-8-11

在 MATLAB Online 中打开

Hmmmm. Not sure why. Here's a fully working example.

originalImage = imread('peppers.png');
subplot(2, 2, 1);
imshow(originalImage);
axis('on', 'image');
[rows, columns, numberOfColorChannels] = size(originalImage);
% Make mask
mask = false(rows, columns);
mask(100:300, 100:400) = true;
subplot(2, 2, 2);
imshow(mask);
axis('on', 'image');
maskedRgbImage = originalImage; % Initialize
% Mask the image using bsxfun() function to multiply the mask by each channel individually.
% Will work regardless if it's grayscale or color.
maskedRgbImage = bsxfun(@times, maskedRgbImage, cast(mask, 'like', maskedRgbImage));
subplot(2, 2, 3);
imshow(maskedRgbImage);
axis('on', 'image');

Sivaramakrishnan Rajaraman 2018-8-13

It works pretty well for me. Many thanks.

请先登录，再进行评论。

更多回答（1 个）

Image Analyst 2018-8-10

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/414381-retrieving-the-lung-pixels-through-matrix-operations-with-images#answer_332377

编辑：Image Analyst 2018-8-10

在 MATLAB Online 中打开

This should work

maskedRgbImage = originalImage; % Initialize
% Convert to RGB, if needed.
[rows, columns, numberOfColorChannels] = size(maskedRgbImage);
if numberOfColorChannels == 1
    maskedImage = cat(3, maskedRgbImage, maskedRgbImage, maskedRgbImage);
end
% Mask the image using bsxfun() function to multiply the mask by each channel individually.
% Will work regardless if it's grayscale or color.
maskedRgbImage = bsxfun(@times, maskedRgbImage, cast(mask, 'like', maskedRgbImage));