How to find which coordinate or pixel (x,y) contains which colour intensity?

I am writing a program where I am able to find the RGB values in the image using the below code

R = a(:,:,1); % Red color plane
G = a(:,:,2); % Green color plane
B = a(:,:,3); % Blue color plane

Now how I can find which coordinate or pixel (x,y) contains which type of colour intensity.

[José-Luis]

What do you mean by "which type of colour intensity"?

[Nimit Jain]

Means how to print the color present on each pixel. Means I need to print (0,0) has 150 value with R color.

[José-Luis]

Sorry but I still don't understand. RGB values go from 0 to 1 in Matlab.

[Guillaume]

By convention, RGB values go from 0 to 1 in matlab if the image is of type double.

The range will be [0 255] if the image is uint8, [0 65535] if the image is uint16, etc.

I too do not understand the question. The red value of pixel at row r and column c is

Rrc = a(r, c, 1);

But surely that is not the question.

[Nimit Jain]

Hi Guillaume Apologise for the confusion. Below is the image details

FileSize  20945
Format  jpg
FormatVersion  ''
Width  200
Height  200
BitDepth  24
ColorType  truecolor

Also I am applying Rrc = a(r,c,1); in my code it's giving me Index exceeds matrix dimensions error. Below is my code I am using

close all
clear all
a = imread('C:\Users\nimitja\Documents\MATLAB\example.jpg');
figure, imshow(a);
[r,c] = size(a);
totalSum = 0;
fileID = fopen('exp.txt','w');
fileID1 = fopen('totalValue.txt', 'w');

binaryImage = a >= 200;
numberOfWhitePixels = sum(binaryImage(:));

R = a(:,:,1); % Red color plane
rSize = size(R);
G = a(:,:,2); % Green color plane
B = a(:,:,3); % Blue color plane

check = a(r,c,1);

% Iterating pixel by pixel of the image and saving value in txt file
%for i=1:r
%    for j=1:c
%        value = double( a(i,j));
%        totalSum = totalSum + value;
%        fprintf(fileID,'%d\n',i,j, value);
%        fprintf(fileID1, '%d\n' , totalSum);
%        disp(a(i,j));
%    end
%end
fclose(fileID);
fclose(fileID1);

Answer 1

[Muhammad Usman Saleem]

On this link , read question and other answers..

find your solution to problem

Answer 2

[Image Analyst]

To interactively see the RGB values, use impixelinfo():

hp = impixelinfo();

Also, you're using size incorrectly:

[r,c] = size(a);

c is the product of the number of columns times the number of color channels, so is basically three times the number of columns and that's why you get an index out of bounds error. See Steve's blog for more info:

http://blogs.mathworks.com/steve/2011/03/22/too-much-information-about-the-size-function/

Don't call your image "a" - that is a not very descriptive name, and it seems like it might be susceptible for you to use "a" again later in the code, blowing away your image, because you forgot you used a. Call it rgbImage instead. It's much more descriptive. To use it correctly, do

[rows, columns, numberOfColorChannels] = size(rgbImage);

[Nimit Jain]

Thanks for the guidance. I have written the below code to iterate pixel by pixel in the image and also get the R, G, B value according to (x,y) position in a txt file. Below is the code. Can you review and tell me if there is anything wrong in it.

close all
clear all
rgbImage = imread('C:\Users\Desktop\Documents\MATLAB\example.jpg');
reducedA = zeros(size(rgbImage));
figure, imshow(rgbImage);
[r,c,numberOfColorChannels] = size(rgbImage);
totalSum = 0;
fileID = fopen('pixels.txt','w');
fileIDR = fopen('rColor.txt','w');
fileIDG = fopen('gColor.txt','w');
fileIDB = fopen('bColor.txt','w');
fileID1 = fopen('totalValue.txt', 'w');
fileIDCompressed = fopen('CompressedPixels.txt', 'w');

binaryImage = rgbImage >= 200;
numberOfWhitePixels = sum(binaryImage(:));

% Here we are calculating the color (RGB) planes
R = rgbImage(:,:,1); % Red color plane
G = rgbImage(:,:,2); % Green color plane
B = rgbImage(:,:,3); % Blue color plane

[rR, cR, zR] = size(R);
[rG, cG, zG] = size(G);
[rB, cB, zB] = size(B);

% try to print all Red value present in the image
for iR = 1:rR
    for jR = 1:cR
        valueR = double(R(iR,jR));
        fprintf(fileIDR,'%d %d %d\n',iR,jR, valueR);
    end
end

% try to print all Green value present in the image
for iG = 1:rG
    for jG = 1:cG
        valueG = double(G(iG,jG));
        fprintf(fileIDG,'%d %d %d\n',iG,jG, valueG);
    end
end

% try to print all Blue value present in the image
for iB = 1:rB
    for jB = 1:cB
        valueB = double(B(iB,jB));
        fprintf(fileIDB,'%d %d %d\n',iB,jB, valueB);
    end
end

% Iterating pixel by pixel of the image and saving value in txt file
for i=1:r
    for j=1:c
        value = double(rgbImage(i,j));
        compressedValue = value - 1;

          totalSum = totalSum + value;
          fprintf(fileIDCompressed, '%d %d %d \n',i,j, compressedValue);
          fprintf(fileID,'%d %d %d \n',i,j, value);
          fprintf(fileID1, '%d\n' , totalSum);
  %        disp(a(i,j));
      end
  end

fclose(fileID);
fclose(fileID1);
fclose(fileIDR);
fclose(fileIDG);
fclose(fileIDB);
fclose(fileIDCompressed);

[Guillaume]

There are a few things in your code that show a lack of understanding of how images are represented and of how matlab works. I would suggest grabbing a book on image processing in your favorite library:

[rR, cR, zR] = size(R);
[rG, cG, zG] = size(G);
[rB, cB, zB] = size(B);

R, G, and B are the three colour planes of your images. The z* is always going to be 1, there's no point asking for it. The size of the colour planes is going to be the same as the size of the image, so rR == rG == rB == r, same for c*. In other words, the above three queries are completely unnecessary. You've already got the information.

valueR = double(something uint);
%then simply print value

There's absolutely no point in converting to double. The exact same number will be printed than if you hadn't bothered.

value = double(rgbImage(i,j));
%...
fprintf(format, value)

Note that rgbImage is an r x c x 3 matrix. You haven't specified the 3rd dimension index in the above, so due to the way matlab indexing work, it's simply 1. Therefore that loop is only going over the red pixels, same as your first loop.

Note that to make it easier to spot bugs, I would move the initialisation of totalsum to 0 just before the loop starts.

Answer 3

[Guillaume]

A much simpler way to save your pixels to text files would be:

rgbImage = imread('C:\Users\Desktop\Documents\MATLAB\example.jpg');
[height, width, ~] = size(rgbImage); %height and width are more meaningful than r and c.
RedChannel = rgbImage(:, :, 1);
GreenChannel = rgbImage(:, :, 2);
BlueChannel = rgbImage(:, :, 3);
[rows, columns] = ndgrid(1:height, 1:column);
%because you save by rows and matlab is column based, we need to transpose all the arrays before reshaping them into one column
%it can then be written as one matrix 
rows = reshape(rows.', [], 1);
columns = reshape(columns.', [], 1);
dlmwrite('rColor.txt', [rows, columns, reshape(RedChannel.', [], 1)], ' ');
dlmwrite('gColor.txt', [rows, columns, reshape(GreenChannel.', [], 1)], ' ');
dlmwrite('bColor.txt', [rows, columns, reshape(BlueChannel.', [], 1)], ' ');

I have no idea what you're trying to do with your last loop, but I'm certain you don't need a loop.

[Nimit Jain]

I did the transpose. In the last loop I am trying to make a copy of the image and showing it. Then I want to remove the high intensity colours. How I will remove the high intensity colours. This is the reason I was thinking to get all the R,G,B and then accordingly operate on the x,y coordinates.

I have also attached the image I am using. Can you please see the code once again.

close all
clear all
rgbImage = imread('C:\Users\nimitja\Documents\MATLAB\example.jpg');
reducedImage = zeros(size(rgbImage));
figure, imshow(rgbImage);
[height,width,~] = size(rgbImage);

[rows, columns] = ndgrid(1:height, 1:width);

rows = reshape(rows.', [], 1);
columns = reshape(columns.', [], 1);

totalSum = 0;
fileID = fopen('pixels.txt','w');
fileIDR = fopen('rColor.txt','w');
fileIDG = fopen('gColor.txt','w');
fileIDB = fopen('bColor.txt','w');
fileID1 = fopen('totalValue.txt', 'w');
fileIDCompressed = fopen('CompressedPixels.txt', 'w');

binaryImage = rgbImage >= 200;
numberOfWhitePixels = sum(binaryImage(:));

RedChannel = rgbImage(:, :, 1);
GreenChannel = rgbImage(:, :, 2);
BlueChannel = rgbImage(:, :, 3);
[rows, columns] = ndgrid(1:height, 1:width);
%because you save by rows and matlab is column based, we need to transpose all the arrays before reshaping them into one column
%it can then be written as one matrix 
rows = reshape(rows.', [], 1);
columns = reshape(columns.', [], 1);

dlmwrite('rColor.txt', [rows, columns, reshape(RedChannel.', [], 1)], ' ');
dlmwrite('gColor.txt', [rows, columns, reshape(GreenChannel.', [], 1)], ' ');
dlmwrite('bColor.txt', [rows, columns, reshape(BlueChannel.', [], 1)], ' ');

% Iterating pixel by pixel of the image and saving value in txt file
for i=1:height
    for j=1:width
        for d=1:3
            value = double(rgbImage(i,j,d));

              compressedValue = rgbImage(i,j,d);
              reducedImage(i,j,d) = rgbImage(i,j,d);
              totalSum = totalSum + value;
              fprintf(fileIDCompressed, '%d %d %d %d \n',i,j,d, compressedValue);
              fprintf(fileID,'%d %d %d %d \n',i,j,d, value);
              fprintf(fileID1, '%d\n' , totalSum);
          end

      end

end

meanIntensityValue = mean(rgbImage(:));
figure, imshow(uint8(reducedImage)), title('Reduced Image');
imwrite((uint8(reducedImage)), 'compressedImage.jpg');
fclose(fileID);
fclose(fileID1);
fclose(fileIDR);
fclose(fileIDG);
fclose(fileIDB);
fclose(fileIDCompressed);

[Guillaume]

Please see my comments to the code you've posted in Image Analyst's answer. They apply to the code above, in addition you now have lines duplicated in two separate places. You need to think a bit more about what you are writing.

You'll note that there is absolutely no fopen whatsoever in my answer. dlmwrite takes care of opening and closing the file for you. In fact, I would suspect that with your code dlmwrite will fail with a permission error because the file is already open.

I still have no idea what you're trying to do with your loop. At present, all it does is equivalent to:

reducedImage = rgbImage; 
totalSum = sum(double(rgbImage(:))); %conversion to double is required here to avoid overflowing uint8

I don't know what you mean by "I want to remove the high intensity colours". What do you call high intensity? And what do you mean by remove? You can't remove pixels from an image (unless you crop it into a smaller rectangle), you can only change their colour.

[Nimit Jain]

I apologise for the confusion. Actually I am want to reduce the size of the image for that I am trying to get the high color density values in an image and then either reduce it or remove it. Sorry for the bad language that created confusion. Also I copy paste my code without filtered it but I corrected that and also remove the extra lines.

Please suggest how to reduce or remove the high color density values/bits.

[Image Analyst]

We're still confused and your most recent explanation did not answer any of Guillaume's questions. Perhaps you want the color frequency image: http://www.mathworks.com/matlabcentral/fileexchange/28164-color-frequency-image

or perhaps you just want imresize(). I have no since your wording doesn't correlate with your code where all you do is copy and image and take it's mean and sum. Sure, the mean or sum is a reduced set of data compared to the original image, but is that what you want or need? Still, I have no idea.

[Guillaume]

"I want to reduce the size of the image" Physical size (i.e. imresize the image)?, the memory footprint but not size, maybe by reducing the number colours and converting to indexed with rgb2ind?

"I am trying to get the high color density values" What does that mean? What is the density of a colour?

"either reduce it or remove it" What does it refer to?

" high color density values/bits" Again what does that mean? Why are you suddenly talking about bits?