Discussion :

[Korithov]

Bonjour à tous!

J'ai un petit soucis avec mon code, je tente de trouver mon erreur mais impossible de mettre la main dessus. J'ai cette erreur qui survient :

Attempted to access v(2,:); index out of bounds because size(v)=[1,5].

Error in CreateAcq (line 47)
        xv = [xv v(i,:)]

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
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
 % We want to create an acquisition
% We request for  BCD
% We know A(x,y)
%  a and b side
%  alpha is random number
%
%
 
k=input('Write the number of satellites you wish : ');
sat=[];
for j=1:k
x=input('Write the first coordinate (West-Est) : ');
y=input('Write the second coordinate (North-South) : ');
 
p(k,:)=[x y];
 
if j==1
sat=[p(k,:)];
else
    sat=[sat;p(k,:)];
end
 
end
xv=[];
yv=[];
 
for i=1:size(sat,1)
p1=p(k,:);
a = 15;
b = 5;
alpha=-45;
% Formula to find B and C
 pa = [x+b*cos(alpha) y+b*sin(alpha)];
 pb = [x-a*sin(alpha) y+a*cos(alpha)];
% D = C + (B-A);
pc = [pb(1)+(pa(1)-p1(1)) pb(2)+(pa(2)-p1(2))];
 
% Creating the area
v=[p1(1) pa(1) pc(1) pb(1) p1(1)];
w=[p1(2) pa(2) pc(2) pb(2) p1(2)];
 
%M(i,:)=[v w]
 
 if i==size(sat,1)
        % if the last satellite do this
        xv = [xv v(i,:)]
        yv = [yv w(i,:)]
 
 
    else
        % else add
        xv = [xv v(i,:) NaN]
        yv = [yv w(i,:) NaN]
 
 end
 
end
 
[X,Y,xq,yq] = Init_Region;
% showing the acquisition
in = inpolygon(X,Y,xv,yv);
figure
axis equal
hold on
plot(xv,yv,'LineWidth',2) % acquisition
plot(X(in),Y(in),'rp') % points inside
plot(X(~in),Y(~in),'b+') % points outside
hold off

[Gooby]

Bonjour,

C'est très simple, tu tentes de récupérer la 2eme ligne d'un vecteur à une ligne et 5 colonnes.
remplaces par v(:,2)
(ce qui dans ton cas reviens à faire v(2))

Ou alors, tu veux réellement récupérer toute la ligne et dans ce cas je pense que tu as oublié de faire ceci:

Code :

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

v=cat(1,v,[p1(1) pa(1) pc(1) pb(1) p1(1)]);

[Gooby]

Un petit conseil de programmation: Lorsque tu fais face à une erreur de ce type, je te conseille de mettre un point d'arrêt à la ligne problématique et de vérifier le contenu des variables pour être sûr qu'elles contiennent bien ce que tu veux qu'elles contiennent. Si elles ne contiennent pas ce que tu désires, tu peux vérifier ligne par ligne ton code afin de toujours vérifier le contenu de tes variables.

[Korithov]

Merci à vous 2 pour vos réponses. Ce problème est résolu en changeant v(i;:) par v(:;i) ;)mais un problème survient ensuite, lors de l'affichage des rectangle j'ai qu'une droite qui apparait, en remontant l'historique, je vois que p1=p(k; :) pause problème donc je change k par i, là j'ai juste une grille

Je voudrais que p1 prenne les valeurs de chaque pk = [x y]

Désolé de vous retenir mais je suis en gros newbie ^^ il y a encore 2jours je ne savais à peine déclarer une variable. J'apprends sur le tas

[Gooby]

Bonjour,

Il va falloir être un peu plus précis que cela:

lors de l'affichage des rectangle j'ai qu'une droite qui apparait

De quels rectangles parles tu? Montre nous le code associé au tracé des rectangles.

je vois que p1=p(k, :) pause problème donc je change k par i, là j'ai juste une grille

Qu'est-ce qui pose problème exaactement? quel est l'erreur? Qu'est-ce que tu entends par "j'ai juste une grille" ?

PS: v(i,:) par v(:,i)

[Korithov]

A la fin du code, je souhaite faire apparaitre les zones crées pour chaque p(k,

Le soucis c'est que j'obtiens la grille de base dès que je met 2 satellites ou plus :



Avec 1 satellite, j'ai toujours la même figure :



Théoriquement, j'ai un(des) rectangle(s) et à l'intérieur de ce(s) rectangle(s) les croix deviennent rouges.


Je remet le code que j'ai actuellement :

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
% We want to create an acquisition
% We request for  BCD
% We know A(x,y)
%  a and b side
%  alpha is random number
%
%
 
k=input('Write the number of satellites you wish : ');
sat=[];
for j=1:k
x=input('Write the first coordinate (West-Est) : ');
y=input('Write the second coordinate (North-South) : ');
 
p(k,:)=[x y];
 
if j==1
sat=[p(k,:)];
else
    sat=[sat;p(k,:)];
end
 
end
fg =size(sat,1);
xv=[];
yv=[];
p1=[];
 
for i=1:fg
p1= p(i,:);
a = 15;
b = 5;
alpha=45;
% Formula to find B and C
 pa = [p1(1)+b*cos(alpha) p(2)+b*sin(alpha)];
 pb = [p1(1)-a*sin(alpha) p(2)+a*cos(alpha)];
% D = C + (B-A);
pc = [pb(1)+(pa(1)-p1(1)) pb(2)+(pa(2)-p1(2))];
 
% Creating the area
v(:,i)=[p1(1);pa(1);pc(1);pb(1);p1(1)];
w(:,i)=[p1(2);pa(2);pc(2);pb(2);p1(2)];
 
%M(i,:)=[v w]
 
 if i==fg
        % if the last satellite do this
        xv = [xv v(:,i)]
        yv = [yv w(:,i)]
 
 
    else
        % else add
        xv = [xv v(:,i) NaN]
        yv = [yv w(:,i) NaN]
 
 end
 
end
 
[X,Y,xq,yq] = Init_Region %Je charge la grille
% showing the acquisition
in = inpolygon(X,Y,xv,yv); %C'est l'ensemble des rectangles
figure
axis equal
hold on
plot(xv,yv,'LineWidth',2) % acquisition
plot(X(in),Y(in),'rp') % points inside
plot(X(~in),Y(~in),'b+') % points outside
hold off