Discussion :

[Kratein]

Bonjour,

Je cherche à calculer l'intégral d'une courbe ( son aire) entre deux points.

J'obtiens la courbe suivant à l'aide de la fonction improfile, pour obtenir le profil d'une image en nuance de gris.


Je recherche à positionner des droites (je sais comment faire) afin de ne calculer l'aire que pour le pic.

Je ne trouve pas et ne comprends pas comment faire le calcul de l'intégral sur cette courbe...
Pouvez-vous m'aider?

Merci d'avance!

[phryte]

Bonjour,

Tu peux utiliser :

Code :

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

trapz(x,y)

Exemple :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
figure
x=0:0.01:100;
y=-0.01*x.^2+x+1;
plot(x,y)
grid
line([10 10],[0 10],'Color','r')
line([90 90],[0 10],'Color','r')
intz=trapz(x(1001:9001),y(1001:9001))
%Vérification
syms xx yy
v=eval(int(-0.01*xx.^2+xx+1,10,90))

[le fab]

salut

il te faut commencer par récuper les vecteur vX et vX de ton tracé (y=f(x))
et le calcul s'opère comme ceci

Code :

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

aire = trapz(vX,vY);

et si veux restreindre le calcul entre x = 25 et x = 75 par exemple :

Code :

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

1
2
idx = (vX >= 25) & (vX <= 75);
aire = trapz(vX(idx),vY(idx));

à ce propos je te suggère un petit tour vers la FAQ

enfin pour tracer les trait verticaux à x = 25 et x = 75

Code :

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

1
2
3
hold on % pour superposer les courbes
plot([25 25], [min(vY) max(vY)],'r-.');
plot([75 75], [min(vY) max(vY)],'r-.');

Fabien

[Kratein]

Merci de votre aide!

Du coup, j'ai récupéré les Y que me fourni improfile, et les X sont la taille de mon image.
Voici le code:

Code :

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

1
2
3
4
5
6
7
        largeur=size(image,2);
        Y=improfile(rectangeGris,x,y);
        vX=[1:largeur];
        vY=Y;
        idx = (vX >= 25) & (vX <= 100);
        %Calcul de l'aire
        aire = trapz(vX(idx),vY(idx));

Merci encore!