Discussion :

[Ciel-nocturne]

Bonjour,

Je souhaiterais représenter graphiquement, pour chacun des canaux R, G et B d'une image, l'intensité normalisée d'un pixel pour tous les pixels d'une même colonne de cette image (profil d'intensité). L'intensité normalisée est la quantité r=R/(R+G+B), g=G/(R+G+B) et b= B/(R+G+B). Le problème est que dans le cas où R=G=B = 0, les quantités r,g et b semblent ne pas être définies (0/0) si bien que les courbes obtenues ne sont pas exploitables.

1) Une commande existe t-elle pour extraire l'intensité normalisée ?
2) Sinon, comment modifier mon code pour faire cela ?

Voici mon code :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
clear all, close all
pkg load image
 
img = imread('C:/Users/Toto/Desktop\image.jpg'); % Read image
 
red = img(:,:,1);  
green = img(:,:,2); 
blue = img(:,:,3);  
 
red_curve = red(:,200);        % column 200
green_curve = green(:,200);  
blue_curve = blue(:,200);    
 
%Normalization
red_curve_norm = red_curve./(red_curve + green_curve + blue_curve);     
green_curve_norm = green_curve./(red_curve + green_curve + blue_curve); 
blue_curve_norm = blue_curve./(red_curve + green_curve + blue_curve);   
 
 
% FIGURES
figure
plot(red_curve_norm, 'r'), xlabel('N° pixel'), ylabel('Intensité'), title('Red channel')

[Jerome Briot]

Une idée pour la première composante :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
idx = (red_curve + green_curve + blue_curve)>0;
red_curve_norm = zeros(size(red_curve));
red_curve_norm(idx) = red_curve(idx)./(red_curve(idx) + green_curve(idx) + blue_curve(idx));

Voir la FAQ : Qu'est-ce que l'indexage logique ?

[Ciel-nocturne]

Bonjour Jerome

Je viens de tester ton code et il me donne exactement la même chose que le mien.

Par ailleurs, ce qui m'étonne c'est que la somme red_curve_norm(idx)+green_curve_norm(idx)+blue_curve_norm(idx) n'est pas égale à 1. Voici mon code :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
clear all, close all
pkg load image
 
img = imread('C:/Users/Toto/Desktop\image.jpg'); % Read image
 
red = img(:,:,1);   
green = img(:,:,2);
blue = img(:,:,3);  
 
red_curve = red(:,300);    
green_curve = green(:,300);
blue_curve = blue(:,300);    
 
idx = (red_curve + green_curve + blue_curve)> 0;
 
red_curve_norm = zeros(size(red_curve));
red_curve_norm(idx) = red_curve(idx)./(red_curve(idx) + green_curve(idx) + blue_curve(idx));
 
green_curve_norm = zeros(size(green_curve));
green_curve_norm(idx) = green_curve(idx)./(red_curve(idx) + green_curve(idx) + blue_curve(idx));
 
blue_curve_norm = zeros(size(blue_curve));
blue_curve_norm(idx) = blue_curve(idx)./(red_curve(idx) + green_curve(idx) + blue_curve(idx));
 
figure
plot(red_curve_norm(idx)+green_curve_norm(idx)+blue_curve_norm(idx) , 'r'), xlabel('N° pixel'), ylabel('Intensité'), title('sum channel')

[Jerome Briot]

Je viens de tester ton code et il me donne exactement la même chose que le mien.

Non, ton code renvoi NaN si les trois composantes sont nulles alors que

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
red_curve =   [ 0 100 0 255];
green_curve = [55  55 0 255];
blue_curve =  [ 0  10 0 255];
 
red_curve_norm = red_curve./(red_curve + green_curve + blue_curve)
 
idx = (red_curve + green_curve + blue_curve)>0;
red_curve_norm = zeros(size(red_curve));
red_curve_norm(idx) = red_curve(idx)./(red_curve(idx) + green_curve(idx) + blue_curve(idx))

Donne :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
red_curve_norm =
 
         0    0.6061       NaN    0.3333
 
 
red_curve_norm =
 
         0    0.6061         0    0.3333

Maintenant, si cela ne répond pas à ta question, donne nous plus de précisions

Par ailleurs, ce qui m'étonne c'est que la somme red_curve_norm(idx)+green_curve_norm(idx)+blue_curve_norm(idx) n'est pas égale à 1.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
red_curve =   [ 0 100 0 255];
green_curve = [55  55 0 255];
blue_curve =  [ 0  10 0 255];
 
idx = (red_curve + green_curve + blue_curve)>0;
 
red_curve_norm = zeros(size(red_curve));
green_curve_norm = zeros(size(red_curve));
blue_curve_norm = zeros(size(red_curve));
 
red_curve_norm(idx) = red_curve(idx)./(red_curve(idx) + green_curve(idx) + blue_curve(idx));
green_curve_norm(idx) = green_curve(idx)./(red_curve(idx) + green_curve(idx) + blue_curve(idx));
blue_curve_norm(idx) = blue_curve(idx)./(red_curve(idx) + green_curve(idx) + blue_curve(idx));
 
red_curve_norm+green_curve_norm+blue_curve_norm

Donne :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
ans =
 
     1     1     0     1

Peux-tu nous fournir l'image JPEG avec laquelle tu travailles ?

[phryte]

Bonjour,
J'ai essayé cela :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
img = imread('anon.jpg'); % Read image
img=im2double(img);

[Ciel-nocturne]

Bonjour phryte,

Merci pour ta contribution. Effectivement la commande que tu proposes est judicieuse et a permis de résoudre mon problème. Je vais me documenter sur celle-ci pour bien comprendre son intérêt.