selecting and removing all grey pixels in an RGB image

34 次查看（过去 30 天）
Sara
Sara 2012-1-5
评论： Image Analyst 2017-5-3
I have a set of fused optical and infrared images I need to analyse. The optical part of the image is visually greyscale and the infrared part is represented by colored blobs. The size, location, range and intensity of each blob changes with each image. Each image is an RGB 480 * 640 * 3 matrix of rgb values. I want to isolate the thermal section from each image and turn the optical portion to black.
Obviously, I can't apply a simple threshold to this because the "grey" part still ranges from 0:255. Edge detection and normal image techniques aren't working well because the grey part of the image is very noisy and grainy. I don't want to convert the whole image to grey because then I'd lose the colour mapping of the thermal portion.
I thought the best way would be masking out the grey images was to find where each pixel that is the same in each colour dimension.
for i = 0:255
A = im(:,:,1) & im(:,:,2) & im(:,:,3) == i;
end
That is obviously just for one image. I thought that if I could get a matrix of A which shows all the pixel locations where rgb is the same in all dimensions, for each value of i, I could use that to mask the non-thermal parts of my photo.
But the code above only identifies about 20 pixels, and not all of them are in the right locations. I'm not a programmer so none of this is obvious to me.
Where am I going wrong?

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

采纳的回答

Titus Edelhofer
Titus Edelhofer 2012-1-5
Hi,
what about looking for all pixels that are somewhat "greyish"? Instead of R, G and B exactly the same you could test
idx = abs(im(:,:,1)-im(:,:,2))<=tresh & abs(im(:,:,2)-im(:,:,3))<=thresh;
Thresh could be 1 or 2 (give it a try) so a pixel [50 51 50] would be considered grey as well.
Titus
  2 个评论
Sara
Sara 2012-1-5
Thank you! That worked perfectly, setting the threshold to 19 actually stripped out virtually all of the noise. Only a few tiny blobs of noise remain and I can filter them out very easily.
Here's a question though, and maybe I'm still being a bit dumb -- how do I use that to mask my original image and leave the rgb values in place? Do I have to apply it to each channel? Now that I'm trying it, most of the roi functions I knew of seem to require greyscale images. I'm still quite new to Matlab and completely new to image processing.
Sara
Sara 2012-1-5
Ignore that, it's working fine now. Thanks again for your help, it is brilliant.

请先登录，再进行评论。

更多回答（4 个）

Image Analyst
Image Analyst 2012-1-5
You can't deal with the pseudocolored image. You need to obtain the actual grayscale infrared image. If that image is just one of the color planes and the optical image is in the other two color channels, then you're in luck. If not, you're going to have to find out how to obtain the thermal image by itself. I've dealt with thermal cameras so I know it can be done.
  1 个评论
Sara
Sara 2012-1-5
Thank you. Nope, the IR isn't in one channel, which is why I had to search for shades of grey across each channel. I did manage to pull out those interest regions perfectly this time, but that is probably due more to luck from the way I set up the original IR fusion. I wouldn't do it the same way next time, it would probably be safer to take thermal and digital photos separately then fuse them AFTER selecting the interest features.

请先登录，再进行评论。

Sara
Sara 2012-1-5
For info, this is what I used to mask through the code and create a corrected image with only the thermal features and everything else blacked out:
idx = abs(I(:,:,1)-I(:,:,2))<=thresh & abs(I(:,:,2)-I(:,:,3))<=thresh;
figure(2);imshow(idx)
[x1 y1 z1] = size(I);
I2 = I;
for x = 1:x1
for y = 1:y1
for z = 1:z1
if idx(x,y) == 1
I2(x,y,z) = 0;
else I2(x,y,z) = I2(x,y,z);
end
end
end
end
It might not have been necessary loop through the whole image like that but I couldn't think of a better way.
  1 个评论
Titus Edelhofer
Titus Edelhofer 2012-1-5
Hi Sara,
for sake of completeness: in your solution at least the most inner loop could be saved: inside the y-loop it would have been sufficient to write
if idx(x,y)==1
I2(x,y,:) = 0;
end
Titus

请先登录，再进行评论。

Titus Edelhofer
Titus Edelhofer 2012-1-5
Hi Sara,
the trick is to change the image from NxMx3 to (N*M)x3:
% store the size
s = size(I);
% make I a matrix with 3 columns and 480*640 rows:
I = reshape(I, s(1)*s(2), 3);
% apply the "filter:
I(idx(:),:) = 0;
% restore the original size
I = reshape(I, s);
Titus
  4 个评论
显示 2更早的评论
Abdoo
Abdoo 2015-4-8
编辑：Abdoo 2015-4-8
I dont have Idea to how can select and convert the value of the range color from RGB image, for example the below image show the select color and the range of this color want to convert it to black and all the remaining colors convert to white.
Image Analyst
Image Analyst 2015-4-8
Abdoo, I think you missed my comment, between your two, about posting a new, separate question in a separate thread.
In the meantime, see my File Exchange: http://www.mathworks.com/matlabcentral/fileexchange/?term=authorid%3A31862. My demos there are very close to what you want to do.

请先登录，再进行评论。

Gebra maryam Alehegn
How to read multiple image from folder and extract RGB color features????

类别

在 Help Center 和 File Exchange 中查找有关 Modify Image Colors 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by