Identify which color zone a x - y point is in

4 vues (au cours des 30 derniers jours)
Mikkel Ibsen
Mikkel Ibsen le 5 Sep 2017
Hi
I have this image where I have identified the different zones using RBG.
After that I stored the different colorareas from zone 1-9.
Ive placed the image inside a x - y axis.
xlim([0 10]) & ylim([0 10])
Now I have a point value of: point = ([5 2]);
This when i plot this, it is inside the "blue" area stored as zone_6.
The RGB color is: ([0 0 204])
How can I make a code that, when I plot a value it gives me which zone it is in?
Ive placed my code in a zip file, so you don't have to identify the different zones, Ive already done that. But I cant solve how to find the zones from a x-y coordinate.
  2 commentaires
José-Luis
José-Luis le 6 Sep 2017
Modifié(e) : José-Luis le 6 Sep 2017
Two quick comments about your code:
  • Using numbered variables is horrible code™. Best not to do it or it is a near certainty you'll regret it later.
  • You pass nine (numbered) variables to your function and proceed to immediately delete everything. What's the point of that? To add insult to injury, you also hard-code them after having destroyed them.
I don't have the Image Processing Toolbox so it is tricky for me to help you further.
Stephen23
Stephen23 le 6 Sep 2017
"After that I stored the different colorareas from zone 1-9."
As everyone has already told you, creating and accessing numbered variables is invariably a bad way to write code. Read this to know why:

Connectez-vous pour commenter.

Réponse acceptée

Guillaume
Guillaume le 6 Sep 2017
Modifié(e) : Guillaume le 6 Sep 2017
As per José-Luis' comment, using numbered variables is a bad idea. In your case, it's going to make the code harder to expand as well. What if you want 10 zones? Ok, you can copy/paste two more lines. What if you want 100 zones? copy/paste is going to be a nightmare.
You could make your zone creation more clever (and shorter), which at the same time would help with your question, identifying which zone a pixel belongs to.
let's start by defining your zones colours:
zonecolours = [128 128 0 %rgb triplet
193 221 198
126 0 1
127 127 127
205 204 0
0 0 204
229 126 127
0 229 0
185 122 87];
Note that I use the power of matlab: matrices, to just have everything in one variable. It's trivial to add more zones: just add rows.
Now let's find out which zone a pixel belongs to:
%img: a rgb colour image of arbitrary size
[~, pixelzone] = ismember(reshape(img, [], 3), zonecolours, 'rows');
pixelzone = reshape(pixelzone, size(img, 1), size(img, 2));
What I've done here is reshape the image into just rows of pixels where each column is an RGB triplet. Then the second output of ismember tells me which row of zonecolour each pixel belongs to (it'll be 0 if the pixel doesn't belong to any zone). I then reshape that output back into the image shape.
Now you can simply use that pixelzone matrix for your imfill:
for zoneindex = 1:size(zonecolours, 1)
imfill(pixelzone == zoneindex, 'holes');
end
And finding which zone a pixel belongs to is trivial: it's the pixelzone matrix, pixelzone(row, col) tells you which zone the pixel at (row, col) belongs to. If you work in axes coordinates, you just have to convert these coordinates to pixel coordinates.
  3 commentaires
Guillaume
Guillaume le 6 Sep 2017
There are many ways you could do the conversion. I've just discovered the function axes2pix (by doing a search for matlab axes coordinates to image pixel), which you can use like so:
%...
himg = imagesc([0 10], [0 10], img); %get a handle to the image
%...
x = 5; y = 2;
pixelrow = axes2pix(size(img, 1), himg.YData, y);
pixelcol = axes2pix(size(img, 2), himg.XData, x);
Or you could calculate it yourself since you know the mapping:
pixelrow = round(interp1(himg.YData, [1 size(img, 1)], y);
pixelcol = round(interp1(himg.XData, [1 size(img, 2)], x);
Mikkel Ibsen
Mikkel Ibsen le 6 Sep 2017
Modifié(e) : Mikkel Ibsen le 6 Sep 2017
Thanks Guillaume :) just what I was looking for. If I wanted to make the x and y axis logarithmic why doesnt this work?
himg = imagesc(loglog([0.1 10], [1 1000]), img); %get a handle to the image
Is there a way that I can make the picture in logarithmic scale so when I go in with [0.1 600] it says its in zone 7 and not 5? Its a small tester for the code.

Connectez-vous pour commenter.

Plus de réponses (2)

Floriane Madeleine Schreiber
Hello Mikkel,
I think, the first thing to do, is to find the value in your matrix corresponding to your coordinates. For example : x = 5. in your matrix the length of x is 2000. You use a cross multiplication, and you find the index_x = 1000. for y=2, you do the same, but now you find a decimal number index_y=337.8000. You cab round it (index_y=338) or look the value before (index_y=337) and after (index_y=338); Now, i have the two index that i can put in my matrix img; like : a=img(338,1000,:);
and a will be :
a(:,:,1) = 0 a(:,:,2) = 0 a(:,:,3) = 204
so you have your RGB code for your x and y coordinates.
now, you can add some if-condition with the number like :
b=[a(:,:,1) a(:,:,2) a(:,:,3)];
if b(1)==0 && b(2)==0 && b(3)==204 sprintf('it is blue and it is in the zone_6') zone='zone_6'; end
with your different RGB codes !
Hope it helps you and answers at your question. Don't hesitate to ask me if it is not clear or not complete.
  2 commentaires
Mikkel Ibsen
Mikkel Ibsen le 6 Sep 2017
Hey Floriane
Im a bit confused on how you use cross multiplication to get the x value in the image?
Can you write the code for me to find the x and y component in the image? (x = 5) == (index_x = 1000) and (y = 2) == (index_y = 337.8)
Guillaume
Guillaume le 6 Sep 2017
Modifié(e) : Guillaume le 6 Sep 2017
@Floriane,
Please use the formatting tools to make your post more readable. See the help link that you can click each time you write/edit a comment/answer.

Connectez-vous pour commenter.


Kevin Krüger
Kevin Krüger le 6 Sep 2017
If I understand this correctly, you already specified each Zone by the points which belongs to it. Now to solve your Problem, you have to go the other direction. With a given Point you can determine the corresponding color. And with the color you can determine the corresponding Zone. I think this is what Floriane is doing.

Catégories

En savoir plus sur Modify Image Colors dans Help Center et File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by