Discussion :

[lunalutcha]

Bonjour,

voila c'est un programme basé sur les descripteurs de Fourier (issu du livre Feature Extraction and Image Processing que mon encadreur m'a donné)

Sur ce dernier il y a un programme que j'ai du mal à compiler. En fait c'est une fonction créée pour laquelle je dois connaitre le paramètre qui est "curve"

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
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
%Angular function
function AngFuncDescrp(curve)
%Function
X=curve(1,:);Y=curve(2,:);
M=size(X,2);%number points
%Arc length
S=zeros(1,m);
S(1)=sqrt((X(1)-X(m))^2+(Y(1)-Y(m))^2);
for i=2:m
    S(i)=S(i-1)+sqrt((X(i)-X(i-1))^2+(Y(i)-Y(i-1))^2);
end
L=S(m);
%Normalised Parameter
t=(2*pi*S)/L;
%Graph of the curve
subplot(3,3,1);
plot(X,Y);
mx=max(max(X),max(Y))+10;
axis([0,mx,0,mx]); axis square; %Aspect ratio
%Graph of the angular function y’/x’
avrg=10;
A=zeros(1,m);
for i=1:m
    x1=0; x2=0; y1=0; y2=0;
    for j=1:avrg
        pa=i-j; pb=i+j;
        if(pa<1) pa=m+pa;
        end
        if(pb>m) pb=pb-m;
        end
        x1=x1+X(pa); y1=y1+Y(pa);
        x2=x2+X(pb); y2=y2+Y(pb);
    end
    x1=x1/avrg; y1=y1/avrg;
    x2=x2/avrg; y2=y2/avrg;
    dx=x2-x1; dy=y2-y1;
    if(dx==0) dx=.00001;
    end
    if dx>0 & dy>0
        A(i)=atan(dy/dx);
    elseif dx>0 & dy<0
        A(i)=atan(dy/dx)+2*pi;
    else
        A(i)=atan(dy/dx)+pi;
    end
end
subplot(3,3,2);
plot(S,A);
axis([0,S(m),-1,2*pi+1]);
%Cumulative angular G(s)=-2pi
G=zeros(1,m);
for i=2:m
    d=min(abs(A(i)-A(i-1)),abs(abs(A(i)-A(i-1))-2*pi));
    if d>.5
        G(i)=G(i-1);
    elseif (A(i)-A(i-1))<-pi
        G(i)=G(i-1)-(A(i)-A(i-1)+2*pi);
    elseif (A(i)-A(i-1))>pi
        G(i)=G(i-1)-(A(i)-A(i-1)-2*pi);
    else
        G(i)=G(i-1)-(A(i)-A(i-1));
    end
end
subplot(3,3,3);
plot(S,G);
axis([0,S(m),-2*pi-1,1]);
%Cumulative angular Normalised
F=G+t;
subplot(3,3,4);
plot(t,F);
axis([0,2*pi,-2*pi,2*pi]);

Je voudrais aussi pouvoir créer une courbe en forme d'avion si c'est possible

voila si vous pouvez m'aider à comprendre et surtout quels sont les paramètres à mettre lors de l'appel de cette fonction

merci

[Jerome Briot]

paramètre qui est "curve"

curve semble être simplement un tableau de dimension 2xn où n serait le nombre de points du profil, non ?

Je voudrais aussi pouvoir créer une courbe en forme d'avion si c'est possible

Tu veux dire, tracer toi-même le profile sous MATLAB avec la souris ?
Dans ce cas, tu peux utiliser la fonction GINPUT

[lunalutcha]

j'ai essayé pour la fonction ginput mais sa reste comme même vague j'arrive toujours pas a tracer ma figure
pour les paramètre curve j'ai essayé d'introduire aussi pour n=50 si j'ai bien compris
il trouve le x et le x parcontre arriver a S il affiche cette erreur je crois que c m qui est indefinie

Code :

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

1
2
Error in ==> AngFuncDescrp at 8
S=zeros(1,m);

voila svp si vous pouvez m'aider merci

[Jerome Briot]

Il faut remplacer la ligne suivante :

Code :

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

M=size(X,2);%number points

par

Code :

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

m=size(X,2);%number points

[lunalutcha]

voila j'ai déjà eu cette idée de le remplacer
mais bon j'ai le premier graph
il reste le problem de x bon l'erreur est

Code :

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

1
2
3
4
??? Attempted to access X(0); index must be a positive integer or logical.
 
Error in ==> AngFuncDescrp at 32
x1=x1+X(pa); y1=y1+Y(pa);

voila j'ai essayé de modifier la boucle enfin la valeur de pa sa ne marche pas

[Jerome Briot]

Voir la : Subscript indices must either be real positive integers or logicals.

[lunalutcha]

ben donc le problème je l'ai compris et je trouve sa logique
enfin mais je dois donc changer le boucle

Code :

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

1
2
3
4
if(pa<1) pa=m+pa 
end
if(pb>m) pb=pb-m
end

de maniere a se que ''pa'' ne prenne pas la valeur nul vu que quand j'ai compiler elle prend deux valeur

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
pa =
 
     1
 
 
pb =
 
     1
 
 
pa =
 
     0
 
 
pb =
 
     2

donc voila j'ai changer la valeur de curve d'éviter cette erreur mais c'est pas les graphe que je voulais avoir enfaite voila les graphes