Discussion :

[procrastination]

Bonjour,
J'ai des mesures de terrain avec données GPS et autres grandeurs. J'ai en outre point des repères sur mon parcours GPS que je voudrais retrouver.

Pour résumer, j'ai mes mesures GPS:

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
Latitude 	Longitude
4,02059	49,282936
4,020588	49,282933
4,020587	49,282931
4,020585	49,282929
4,020584	49,282927
4,020582	49,282925
4,020581	49,282923
4,02058	49,282921
4,020578	49,282919
4,020577	49,282916
4,020575	49,282914
4,020574	49,282912
4,020572	49,28291
4,020571	49,282908
4,02057	49,282906
4,020569	49,282903
4,020568	49,282901
4,020567	49,282899
4,020567	49,282897
4,020566	49,282894
4,020565	49,282892
4,020564	49,28289
4,020563	49,282887
4,020562	49,282885
4,020562	49,282883
4,020561	49,282881
4,02056	49,282878
4,020559	49,282876
4,020558	49,282874
4,020558	49,282872
4,020557	49,282869
4,020556	49,282867
4,020555	49,282865
4,020554	49,282862
4,020554	49,28286
4,020553	49,282858
4,020553	49,282856
4,020553	49,282853
4,020553	49,282851
4,020553	49,282849
4,020553	49,282846

J'ai un point de repère :
Latitude Longitude
4,020593 49,28294 Station 5.B

Je souhaiterai retrouver mon point de repère dans mes coordonnées GPS mesurées afin d'effectuer un recalage de mes mesures.

Dois-je calculer la distance euclidienne (ou une autre distance) pour chaque point ? Existe-t-il une méthode rapide et fiable pour effectuer cela ?

Merci de votre réponse.

Mehdi

[Gooby]

Bonjour,

Tu peux éventuellement commencer par un

Code :

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

1
2
toto = ismember(coordonnées_gps,point_de_repere,'rows')
ligne_point_de_repere = find(toto);

[le fab]

Salut

Dois-je calculer la distance euclidienne (ou une autre distance) pour chaque point ?

comme je suppose peu probable que ton point repère soit dans ton jeux de mesures, je dirais oui
sinon si tu es sur que ton point repère est exactement dans ton jeux de mesures, alors cf méthode expliqué dans le post ci dessus

Existe-t-il une méthode rapide et fiable pour effectuer cela ?

comme ceci :

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
mes = [4.02059	49.282936
4.020588	49.282933
4.020587	49.282931
4.020585	49.282929
4.020584	49.282927
4.020582	49.282925
4.020581	49.282923
4.02058	49.282921
4.020578	49.282919
4.020577	49.282916
4.020575	49.282914
4.020574	49.282912
4.020572	49.28291
4.020571	49.282908
4.02057	49.282906
4.020569	49.282903
4.020568	49.282901
4.020567	49.282899
4.020567	49.282897
4.020566	49.282894
4.020565	49.282892
4.020564	49.28289
4.020563	49.282887
4.020562	49.282885
4.020562	49.282883
4.020561	49.282881
4.02056	49.282878
4.020559	49.282876
4.020558	49.282874
4.020558	49.282872
4.020557	49.282869
4.020556	49.282867
4.020555	49.282865
4.020554	49.282862
4.020554	49.28286
4.020553	49.282858
4.020553	49.282856
4.020553	49.282853
4.020553	49.282851
4.020553	49.282849
4.020553	49.282846];
rep = [4.020593	49.28294];
>> dist = ((mes(:,1)-rep(1)).^2 + (mes(:,2)-rep(2)).^2).^0.5;

ensuite utilise la fonction min (avec la deuxième sortie) pour trouver le(s) point(s) le plus proche de ton repère reference

[procrastination]

Merci pour votre réponse ;-)