Coordonnées GPS, mille nautique, cercle.. problème réel urgent*!

[invitecd41c164]

Bonjour à tous,

Je suis marin à bord d'un «*navire de garde*»*; nous sommes à poste 24h/24 7j/7 pendant un mois pour nous assurer qu'aucun navire ne viendra perturber des travaux sous-marin (en l’occurrence la réparation d'un câble)

Pour se faire nous utilisons une cartographie numérique couplée au radar qui nous permet de connaître la distance entre les navires de passage et notre zone de garde. Hier un de nos ordinateurs a lâché, nous laissant notre seul ordinateur de réserve pour assurer la surveillance. Il est équipé d'un logiciel de cartographie très basique, qui n'autorise rien de plus que le tracé d'une (unique) «*route*» à partir de coordonnées de points.

Voici notre problème*:

Notre zone garde est circulaire, de rayon 0,5 mille nautique
Nous connaissons les coordonnées GPS du point central de ce cercle, définit en degré, minute, décimale*: 47°09.78'N – 002°25.34'W

Il me faut afficher la zone de garde sur notre interface cartographique. Comment pourrait-on trouver les coordonnées des points du cercle*? J'aimerais pouvoir tracer grossièrement ce cercle par une «*route*» qui lirait un certain nombre de ses points.

Je suis à l'écoute de toutes vos idées même si ce ne sont que des pistes, les travaux débutent le 26*! Merci
-----

[obi76]

Bonjour,

allez, je vais être sympa :

 lat :    47.1635170     lon :   -2.41415215    
 lat :    47.1645584     lon :   -2.40602112    
 lat :    47.1661072     lon :   -2.39797235    
 lat :    47.1681709     lon :   -2.39003801    
 lat :    47.1707268     lon :   -2.38225007    
 lat :    47.1737747     lon :   -2.37464023    
 lat :    47.1772919     lon :   -2.36723876    
 lat :    47.1812859     lon :   -2.36007571    
 lat :    47.1857147     lon :   -2.35317993    
 lat :    47.1905708     lon :   -2.34657907    
 lat :    47.1958427     lon :   -2.34029984    
 lat :    47.2015038     lon :   -2.33436728    
 lat :    47.2075233     lon :   -2.32880545    
 lat :    47.2138824     lon :   -2.32363677    
 lat :    47.2205582     lon :   -2.31888199    
 lat :    47.2275238     lon :   -2.31456041    
 lat :    47.2347488     lon :   -2.31068921    
 lat :    47.2422028     lon :   -2.30728412    
 lat :    47.2498589     lon :   -2.30435872    
 lat :    47.2576942     lon :   -2.30192471    
 lat :    47.2656555     lon :   -2.29999256    
 lat :    47.2737312     lon :   -2.29856896    
 lat :    47.2818756     lon :   -2.29766059    
 lat :    47.2900658     lon :   -2.29727077    
 lat :    47.2982635     lon :   -2.29740071    
 lat :    47.3064346     lon :   -2.29805017    
 lat :    47.3145485     lon :   -2.29921675    
 lat :    47.3225708     lon :   -2.30089569    
 lat :    47.3304672     lon :   -2.30308008    
 lat :    47.3382187     lon :   -2.30576110    
 lat :    47.3457794     lon :   -2.30892801    
 lat :    47.3531227     lon :   -2.31256819    
 lat :    47.3602180     lon :   -2.31666660    
 lat :    47.3670425     lon :   -2.32120752    
 lat :    47.3735657     lon :   -2.32617140    
 lat :    47.3797569     lon :   -2.33153939    
 lat :    47.3856010     lon :   -2.33728981    
 lat :    47.3910637     lon :   -2.34339881    
 lat :    47.3961372     lon :   -2.34984207    
 lat :    47.4007797     lon :   -2.35659409    
 lat :    47.4049950     lon :   -2.36362696    
 lat :    47.4087486     lon :   -2.37091279    
 lat :    47.4120293     lon :   -2.37842250    
 lat :    47.4148369     lon :   -2.38612509    
 lat :    47.4171448     lon :   -2.39399028    
 lat :    47.4189529     lon :   -2.40198588    
 lat :    47.4202461     lon :   -2.41008019    
 lat :    47.4210320     lon :   -2.41824007    
 lat :    47.4212875     lon :   -2.42643309    
 lat :    47.4210320     lon :   -2.43462610    
 lat :    47.4202461     lon :   -2.44278598    
 lat :    47.4189529     lon :   -2.45088005    
 lat :    47.4171448     lon :   -2.45887589    
 lat :    47.4148369     lon :   -2.46674109    
 lat :    47.4120293     lon :   -2.47444415    
 lat :    47.4087486     lon :   -2.48195338    
 lat :    47.4049950     lon :   -2.48923922    
 lat :    47.4007797     lon :   -2.49627209    
 lat :    47.3961372     lon :   -2.50302410    
 lat :    47.3910637     lon :   -2.50946760    
 lat :    47.3856010     lon :   -2.51557660    
 lat :    47.3797569     lon :   -2.52132678    
 lat :    47.3735657     lon :   -2.52669477    
 lat :    47.3670425     lon :   -2.53165913    
 lat :    47.3602180     lon :   -2.53619957    
 lat :    47.3531227     lon :   -2.54029799    
 lat :    47.3457794     lon :   -2.54393816    
 lat :    47.3382187     lon :   -2.54710507    
 lat :    47.3304672     lon :   -2.54978609    
 lat :    47.3225708     lon :   -2.55197048    
 lat :    47.3145485     lon :   -2.55364943    
 lat :    47.3064346     lon :   -2.55481601    
 lat :    47.2982635     lon :   -2.55546546    
 lat :    47.2900658     lon :   -2.55559540    
 lat :    47.2818756     lon :   -2.55520558    
 lat :    47.2737312     lon :   -2.55429697    
 lat :    47.2656555     lon :   -2.55287337    
 lat :    47.2576942     lon :   -2.55094123    
 lat :    47.2498589     lon :   -2.54850721    
 lat :    47.2422028     lon :   -2.54558182    
 lat :    47.2347488     lon :   -2.54217672    
 lat :    47.2275238     lon :   -2.53830552    
 lat :    47.2205582     lon :   -2.53398395    
 lat :    47.2138824     lon :   -2.52922916    
 lat :    47.2075233     lon :   -2.52406049    
 lat :    47.2015038     lon :   -2.51849866    
 lat :    47.1958427     lon :   -2.51256609    
 lat :    47.1905708     lon :   -2.50628686    
 lat :    47.1857147     lon :   -2.49968600    
 lat :    47.1812859     lon :   -2.49279022    
 lat :    47.1772919     lon :   -2.48562717    
 lat :    47.1737747     lon :   -2.47822571    
 lat :    47.1707268     lon :   -2.47061586    
 lat :    47.1681709     lon :   -2.46282792    
 lat :    47.1661072     lon :   -2.45489359    
 lat :    47.1645584     lon :   -2.44684482    
 lat :    47.1635170     lon :   -2.43871379    
 lat :    47.1629982     lon :   -2.43053293    
 lat :    47.1629982     lon :   -2.42233562    
 lat :    47.1635170     lon :   -2.41415501

Vérifiez quand même mais ça doit être bon.

Sinon carte + compas ça se fait immédiatement...

[invitecd41c164]

Merci beaucoup. Comment avez-vous fait ? Y a t-il une formule que je pourrais appliquer directement à partir des coordonnées du centre du cercle et de son rayon ? (car la zone est amenée à changer)

[obi76]

Vous prenez le centre (latitude / longitude), vous convertissez en km (y = latitude*rayon de le terre*2 pi / 360 ; x = longitude*rayon de la terre*2 pi / 360)
Vous prenez un angle de 0 à 360 °, vous ajoutez le cercle : (y + rayon*sin(angle); x + rayon*cos(angle)) ; avec le rayon de votre zone en [km] (pas en miles)
Et vous reconvertissez chaque point obtenu en angle :
lat = y*360/(2*pi*rayon de la terre)
lon = x*360/(2*pi*rayon de la terre)

Mais je ne comprends pas pourquoi vous ne le faites pas à l'ancienne, c'est quand même nettement plus simple : vous prenez votre carte, vous faites un cercle, et vous prenez plusieurs points sur ce cercle. Ca sera nettement plus rapide que de vous amuser à calculer les points avec cette méthode un par un (à moins que quelqu'un chez vous sache programmer, dans ce cas c'est une option aussi).

[A voir en vidéo sur Futura]

[Biname]

Salut,

Mais je ne comprends pas pourquoi vous ne le faites pas à l'ancienne, c'est quand même nettement plus simple : vous prenez votre carte, vous faites un cercle, et vous prenez plusieurs points sur ce cercle. Ca sera nettement plus rapide que de vous amuser à calculer les points avec cette méthode un par un (à moins que quelqu'un chez vous sache programmer, dans ce cas c'est une option aussi).

Je me suis posé la même question !

Les navires n'apparaissent pas sur une carte papier , faut reporter les coordonnées radar et/ou AIS à la main, ils n'en ont pas le temps. 1/2 miles = 90 secondes à 20 noeuds.
N'importe quel logiciel est préférable. Les logiciels proposent des fonctionnalités traçage de route basées sur une série de points. ...

Biname

[Biname]

Salut,

Il me semble qu'une zone d'exclusion est possible sur QTVLM ?????
https://www.meltemus.com/index.php/f...ne-d-exclusion

Les logiciels gratuits OPENCPN, QTVLM(français, contacter Maitai), ...

Biname

[obi76]

Les navires n'apparaissent pas sur une carte papier , faut reporter les coordonnées radar et/ou AIS à la main, ils n'en ont pas le temps. 1/2 miles = 90 secondes à 20 noeuds.
N'importe quel logiciel est préférable. Les logiciels proposent des fonctionnalités traçage de route basées sur une série de points. ...

Ha ceci expliquant cela, merci

[jiherve]

bonjour,
savez vous quel est le système de projection utilisé et s'il y en a un?
Où votre navire se situe t il sur l'image ?
Y a t il une fonction de zoom ?
Les plots radar apparaissent ils sur l'écran ?
Si la réponse à la dernière question est oui alors il y a de fortes chance que votre navire soit au centre et qu'il n'y ait aucune projection ou une projection sympa type Lambert glissant et si vous connaissez le facteur d’échelle alors un bête bout de transparent sur lequel vous aurez tracé un cercle placé devant l’écran fera le boulot, on peut même tracer directement sur l'écran avec des feutres si possible effaçables.
S'il y a différents zoom discrets alors il faudra faire plusieurs cercles de couleurs ou autre moyen typo pour matérialiser votre zone de veille.
JR

[invitecd41c164]

Merci à tous pour vos super réponses ! Je vous réponds en détail à mon prochain quart en milieu de nuit. Bonne soirée à vous

[jiherve]

bonsoir
bon quart je pense que nous serons couchés!

Il y a trois sortes d'hommes : les vivants, les morts, et ceux qui vont sur la mer.

Aristote apocryphe
JR

[SK69202]

Mais je ne comprends pas pourquoi vous ne le faites pas à l'ancienne,

On celle ante natum de cocolastic, qu'on pratiquait dans la tempête sur le rail au large d'Ouessant, avec le crayon gras, la règle à CPA et le tracé sur le plot anti-parallaxe du scope à tube cathodique ?
"Une main pour toi, une main pour le bateau, et si le bateau souffre rajoute y trois doigts"....

[invite12e4e957]

Bonjour,
Je ne trouve pas les mêmes résultats qu' Obi76 :
- le centre du cercle devrait normalement se situer aux coordonnées décimales [47.16300, -2.42233] ;
- vu qu' on n'est pas sur l'équateur, 1 degré en latitude et 1 degré en longitude ne donnent pas la même distance. Et en France, ce n'est pas anodin. Il faut donc entrer l' équation d'une ellipse, pas celle d' un cercle. Explications :
https://fr.wikipedia.org/wiki/Coordo...%A9ographiques
- Autre souci, 0.5 mille correspond à 1/120e de degré en latitude, soit 0.00833 degré. Là, les coordonnées des points décrivent un cercle bien plus grand.

Perso, je trouverais plutôt cela :

47.17133 , -2.42233
47.1712 , -2.42021
47.17083 , -2.41816
47.17021 , -2.41623
47.16938 , -2.41448
47.16835 , -2.41298
47.16717 , -2.41176
47.16585 , -2.41086
47.16445 , -2.41031
47.163 , -2.41012
47.16155 , -2.41031
47.16015 , -2.41086
47.15884 , -2.41176
47.15765 , -2.41298
47.15662 , -2.41448
47.15579 , -2.41623
47.15517 , -2.41816
47.1548 , -2.42021
47.15467 , -2.42233
47.1548 , -2.42445
47.15517 , -2.42651
47.15579 , -2.42844
47.15662 , -2.43018
47.15765 , -2.43169
47.15884 , -2.43291
47.16015 , -2.43381
47.16155 , -2.43436
47.163 , -2.43454
47.16445 , -2.43436
47.16585 , -2.43381
47.16717 , -2.43291
47.16835 , -2.43169
47.16938 , -2.43018
47.17021 , -2.42844
47.17083 , -2.42651
47.1712 , -2.42445

Cocolastic, je te laisse tracer ces points sur ton PC. Dis-moi si souci.

[Biname]

Salut les gras de la narine,
0.5 mille nautique = 1/2 minute d'arc au centre de la Terre
Le diamètre de ce cercle doit faire en lattitude une minute d'arc

OBI76 : Max_Lat 47.4212875 - Min_Lat 47.1635170 = 0.2577 soit 15.46 minutes ou rayon = 7.7 milles

Centre du cercle : 47°09.78'N – 002°25.34'W ou 47.163°N et 2.422°W
Point Lat Max Nord : 47°10.28
Point Lat Min Nord : 47°09.28

Par 47.163N la circonférence de la Terre vaut http://www.easysurf.cc/circle.htm#cetol1 16961.4 milles
1/2 mille en minute de longitude par 47.163 vaut 360 * 60/(16961.4 * 2) = 0.636 minutes d'arc

Point Max Ouest 002°25.34' + 0.636 = 002°25.967' (carte 2°26.08')
Point Min Ouest 002°25.34' - 0.636 = 002°24.704' (carte 2°24.60')

Ah, tu es repéré ! Une carte avec un cercle ... very simple

Voilà mission accomplie

Biname

[obi76]

- vu qu' on n'est pas sur l'équateur, 1 degré en latitude et 1 degré en longitude ne donnent pas la même distance

en effet, erreur bête... :-/

[SK69202]

1/2 mille en minute de longitude par 47.163 vaut 360 * 60/(16961.4 * 2) = 0.636 minutes d'arc

Si ma mémoire est bonne, 1 nautique = 1 minute d'arc de méridien et c'est constant sur toute la sphère terrestre.
Chadok d'ailleurs le précise.

0.5 mille correspond à 1/120e de degré en latitude

[invitecd41c164]

Déjà merci à tous pour vos réponses, je ne m'attendais vraiment pas à en recevoir autant ! Notre sistership nous a envoyé le tracé de la zone -donc il n'y pas d'urgence- mais je cherche tout de même à comprendre vos méthodes pour pouvoir me débrouiller à l'avenir

« Mais je ne comprends pas pourquoi vous ne le faites pas à l'ancienne, c'est quand même nettement plus simple : vous prenez votre carte, vous faites un cercle, et vous prenez plusieurs points sur ce cercle. Ca sera nettement plus rapide que de vous amuser à calculer les points avec cette méthode un par un (à moins que quelqu'un chez vous sache programmer, dans ce cas c'est une option aussi). »

obi76

« et si vous connaissez le facteur d’échelle alors un bête bout de transparent sur lequel vous aurez tracé un cercle placé devant l’écran fera le boulot, on peut même tracer directement sur l'écran avec des feutres si possible effaçables.
S'il y a différents zoom discrets alors il faudra faire plusieurs cercles de couleurs ou autre moyen typo pour matérialiser votre zone de veille.
JR »

jiherve

En fait nous n'avons ni cartes (il paraît que certains dans la marine marchande naviguent toujours avec des cartes physiques, mais ce n'est pas notre cas) ni feutres, ni papier calque.. on doit faire avec les moyens du bord c'est Apollo 13 notre affaire !

Et juste pour être un peu plus clair dans l’énoncé de mon problème : notre interface cartographique est connecté au GPS et au radar, ce qui permet d'afficher en temps réel la position des navires sur la carte, leur route et leur vitesse. Notre seul soucis est que contrairement à la majorité des programmes, celui-ci ne permet pas de tracer des cercles à partir d'un centre. Il autorise uniquement le tracé d'une route reliant les points dont l'on rentre les coordonnées. Avec un certain nombre de points je pourrais donc tracer un polygone régulier (d'une grosse dizaine de côtés par exemple) qui pourrait s'apparenter à un cercle.

« On celle ante natum de cocolastic, qu'on pratiquait dans la tempête sur le rail au large d'Ouessant, avec le crayon gras, la règle à CPA et le tracé sur le plot anti-parallaxe du scope à tube cathodique ?
"Une main pour toi, une main pour le bateau, et si le bateau souffre rajoute y trois doigts".... »

SK69202

Tout a changé, tout est dématérialisé... ton pays s'est endormi

« Il faut donc entrer l' équation d'une ellipse, pas celle d' un cercle. »

Chadok

Ces points correspondent parfaitement au tracé de la zone ! J'ai essayé d'aller voir du côté des ellipses pour ajuster les équations d'Obi (un grand merci au passage pour le détail de tes calculs) ça m'a refroidit ! Pourrais-tu m'expliquer comment tu es arrivé à ces résultats ? J'aimerais beaucoup pouvoir refaire ces calculs avec un point central et un rayon différent à l'avenir.

« Ah, tu es repéré ! Une carte avec un cercle ... very simple »

Biname

Merci pour tout ! j'ai téléchargé OpenCPM, je vais creuser les fonctionnalités

[SK69202]

Tout a changé, tout est dématérialisé... ton pays s'est endormi

Oui, je sais, j'ai vu la transition, entre ceux qui ont besoin de l'électricité partout sur Terre pour naviguer et ceux déjà pas nombreux, qui savaient où on était dans le Pacifique avec une crépu, une droite de Lune et quelques tables abimées.

[Biname]

Salut,

Si ma mémoire est bonne, 1 nautique = 1 minute d'arc de méridien et c'est constant sur toute la sphère terrestre.
Chadok d'ailleurs le précise.

Oui mais les méridiens sont gradués en degré de latitude et vice versa ! Mon texte manque de rigueur mais mes quatre valeurs N,W,S,E sont 'exacte'.
Biname

[jiherve]

Bonjour,
c'est dommage que vous n'ayez pas de feutre car c'est un instrument basique à avoir si l'on dispose d'un écran couplé à un radar car votre problème était facilement résolu.
Attention cela n'est possible que s'il n'y a aucune déformation induite par une projection carto et que la zone affichée ne couvre que qqs miles.
JR

[Biname]

Salut,
J'ai chercher un calculateur en ligne mais n'ai rien trouvé ????

Pour le calcul (approximation) :
P_lat = latitude du centre
P_lon = longitude centre
R_lat = rayon du cercle en latitude (ici 0.5')
R_lon = rayon du cercle en longitude (ici 0.636' mon calcul ci-dessus, rayon terre à latitude ...)
alpha l'azimut géographique (N=0°, E=90°, S=180°, W=270°)

X un point du cercle
Lat_X = P_lat + R_lat * cos (alpha)
Lon_X = P_lon - R_lon * sin (alpha)

Faudrait vérifier les signes ... laisser un peu de boulot pour les autres et un petit javascript ?

Biname

[invite12e4e957]

Biname,
Tu as fait une erreur de lecture en entrant la latitude 47.163 sur ton site http://www.easysurf.cc/circle.htm#cetol1 :
tu as lu les miles US, et non les milles nautiques. Par 47.163N la circonférence de la Terre vaut 14739 milles nautiques, et non 16961.
Si je reprends la suite,
1/2 mille nautique en minute de longitude par 47.163 vaut 360 * 60/(14739 * 2) = 0.7327 minutes d'arc, et non 0.636.

Mon approximation d'hier est à 0.733, donc Cocolastic peut faire son quart tranquille.
Mais merci pour ton lien easysurf, Biname. Je ne connaissais pas ce site, il y a un tas de choses pratiques

[invitecd41c164]

J'ai chercher un calculateur en ligne mais n'ai rien trouvé ????

Biname

Vassillia du forum "https://www.maths-forum.com/" https://www.maths-forum.com/enigmes/....html#p1457793

a proposé le lien https://www.lexilogos.com/calcul_distances.htm pour calculer les distances entre deux points.

D'ailleurs, je vous prie de m'excusez pour le multi-post, ça vous fait cogiter sans profiter des avancées des autres, je n'avais pas réfléchi à ça, je n'aurais pas dû.
Voici les liens des forums où j'ai posé ma question :

https://www.maths-forum.com/enigmes/....html#p1457793
http://www.bibmath.net/forums/viewto...d=93520#p93520

Je me suis engagé à faire un topo complet de toutes les méthodes proposées pour résoudre ce problème de façon à ce que tout le monde puisse en profiter.

Pour ma part j'utilise ce site pour convertir les unité de mes coordonnées : https://www.lecampingsauvage.fr/gps-convertisseur
En passant une petite anecdote - dans les eaux françaises il y a deux standards en vigueur : le degré, minute, décimal de minute mais aussi le degré, minute, seconde (parfois même utilisé par des officiels). Ca rend mes collègues britanniques fou car si le système d'unité n'est pas précisé il est très facile de faire une erreur de position, les coordonnées s'énonçant de la même façon

[invitecd41c164]

savez vous quel est le système de projection utilisé et s'il y en a un?
Où votre navire se situe t il sur l'image ?
Y a t il une fonction de zoom ?
Les plots radar apparaissent ils sur l'écran ?

Attention cela n'est possible que s'il n'y a aucune déformation induite par une projection carto et que la zone affichée ne couvre que qqs miles.

Jiherve

Notre interface cartographique se présente comme Google Map, il est possible de zoomer et de se déplacer indépendamment de la position du navire même si il y a une option de "centrage" forcé. La projection est celle du standard des cartes marines, le WGS84 :

https://en.wikipedia.org/wiki/World_Geodetic_System

Notre interface radar elle est un peu à l'ancienne, il est possible de faire des zoom mais elle est centrée sur notre navire. Je ne sais pas de quelle projection il s'agit. Je suppose que c'est également une WGS84 car elle semble couplé à notre système carto sans correction (visible du moins)

[invite12e4e957]

Comme tu cherches une méthode pour le futur, je te laisse ce code. Ça liste les 36 points que je t'ai donnés hier.
Je ne pratique que Python, donc c'est du Python, mais tu pourras facilement l'adapter dans un autre langage. Et si jamais tu ne programmes pas, c'est franchement facile à reproduire sur une feuille Excel.
NB : c'est une méthde simplifiée, ça devient approximatif pour des cercles qui font des milliers de km

# Cercle centré en latitude et longitude
# Traçage d'une ellipse, angulaire sur parallèle <> angulaire sur méridien

from math import *

lat_centre_deg = 47.163
long_centre_deg = -2.42233
rayon_mille = 0.5
nb_points = 36

rayon_km = rayon_mille*1.852
diam_terre_km = 40005
parallele_km = diam_terre_km*cos(lat_centre_deg*pi/180)

angle_meridien = rayon_km*360/diam_terre_km
angle_parallele = rayon_km*360/parallele_km
print("Angle sur méridien [deg] :", round(angle_meridien, 5))
print("Angle sur parallèle [deg] :", round(angle_parallele, 5))

step = int(360/nb_points)
for theta in range(0, 360, step):
    lat_point = lat_centre_deg+angle_meridien*cos(theta*pi/180)
    long_point = long_centre_deg+angle_parallele*sin(theta*pi/180)
    print(theta, " ", round(lat_point, 5), " ", round(long_point, 5))

[Biname]

@Chadok : oui erreur de ma part miles US et pas milles nautiques !

@Cocolastic : Je pensais à un modèle qui donnant des points d'un cercle polaire basé sur la trigonométrie sphérique.

@Tous : On trouve des millions de modèles, le truc est d'ajouter 'online calculator'

Je viens de m'amuser deux heures sur un modèle Javascript qui tient dans un HTML !
ChapmanCoordCerclePolaire.zip sans garantie

Biname

[jiherve]

bonsoir
Le WGS84 ne défini pas une projection mais le géoïde de référence et le système de coordonnées affèrent.
Google utilise une projection type Mercator(standard des cartes marines) à partir de coordonnées WGS84 modifiées, la projection Mercator ne conserve que les angles (c'est donc une projection conforme), c'est utile pour la navigation, mais surfaces et distances représentées sont dépendantes de la latitude(voir indicatrice de Tissot), a noter qu'aucun système de projection ne permet de conserver les distances pour des représentations à petite échelle.
Le radar travaille fondamentalement en coordonnées circulaires qu'il est assez facile de convertir en coordonnées cartésiennes dont l'origine du repère instantané est défini par les coordonnées du porteur ; Donc si les plots radars se superposent à une carte alors il y a deux possibilités:
1:La carte est générée de façon telle qu'elle corresponde à un projection présentant un minimum d'erreur (angles, distance)pour la représentation locale (Mercator out) et les données radar sont plotées telles quelles.
2:Les plot radars sont transposés dans le système de projection choisi.

La cartographie c'est compliqué.
JR

[invitecd41c164]

Je vous répondrais plus en détail quand j'aurais un peu de temps mais déjà :

@Chadok : merci ! je ne savais pas du tout ce qu'était Python, je télécharge Anaconda dés que notre parabole arrêtera de pédaler dans la semoule et je m'y colle.

@Biname : ça fonctionne parfaitement ! merci pour ce gros boulot. Je ne m'attendais vraiment pas à tout ça. Vous êtes une communauté très surprenante. Merci encore !

Merci à tous ! Mes collègues sont comme des fous "bloody hell !" "bloody hell !"

@Jiherve : Je ne comprends pas certains mots techniques que tu utilises, je te réponds très vite quand j'aurai un peu creusé ! à première vue je dirais l'option 1 car par exemple le profil de la côte en nuage de points que l'on peut observer sur le radar semble correspondre parfaitement au profil de la côte sur la carte (mais peut-être y a t'il une faible distorsion que je ne parviens pas à discerner ? je ne me rends pas compte du niveau de distorsion auquel il faudrait s'attendre)

Belle journée à tous !

[Biname]

Salut,

On s'amuse ! Et avec toi, on voyage dans l'estuaire de la Loire, plein de vahinés en bigouden.

On trouve en ligne des interpréteurs et compilateurs pour tous les langages.

Ici un interpréteur Python où il suffit de copier/coller le code Python de Chadok msg #24 et cliquer sur RUN pour obtenir les points,
aucune nécessité d'installer Python.

https://www.programiz.com/python-pro...line-compiler/

Et pour vérifier l'adage 'forgeant ... forgeron', nos codes Javascript et Python donnent les mêmes résultats, image attachée.
@Chadok : on est trop fort

Biname
Résultats Javascript à Gauche et Python à droite pour 12 points

[SK69202]

le profil de la côte en nuage de points que l'on peut observer sur le radar semble correspondre parfaitement au profil de la côte sur la carte (mais peut-être y a t'il une faible distorsion que je ne parviens pas à discerner ? je ne me rends pas compte du niveau de distorsion auquel il faudrait s'attendre)

Ce que l'on voit sur l'écran, ce n'est pas exactement la côte, c'est ce que voit le radar en fonction de sa hauteur, de l'altitude et du relief de la côte, avec la distance la coïncidence image radar, trait de côte diminue.
Le profil de la côte a aussi son mot à dire, entre les falaises du Pas de Calais et les terres basses de Belgique ou de Hollande, ce n'est pas la même chose.
Un autre soucis qui peut apparaitre en fonction de la saison et de la météo, on voit parfois la terre de très loin ou vachement plus que d'habitude, c'est souvent couplé avec les effets optiques des mirages et là il peut y avoir de sacrés écarts.
Il y a aussi d'autres contraintes venant de la technologie du radar.

C'est suffisant pour naviguer, mais au temps où on faisait la police, à Ouessant, au Casquets ou Pas de Calais, jamais le PV n'était établi sur les seuls informations radar.

"La terre la plus proche est sous la quille".

[Biname]

Salut,
Il est aussi intéressant dans le code calculant les coordonnées de points du cercle polaire d'ajouter l'erreur en calculant la distance sphérique entre le centre du cercle polaire et les points calculés (idéalement cette valeur est égale au rayon de l'énoncé).

En python, c'est plus court ! Je reprends donc le code de Chadok et j'ajoute le calcul de la distance

# Cercle centré en latitude et longitude
# Traçage d'une ellipse, angulaire sur parallèle <> angulaire sur méridien

from math import *

lat_centre_deg = 47.163
long_centre_deg = -2.42233
rayon_mille = 0.5
nb_points = 12
rayon_km = rayon_mille*1.852
diam_terre_km = 40005
parallele_km = diam_terre_km*cos(lat_centre_deg*pi/180)

angle_meridien = rayon_km*360/diam_terre_km
angle_parallele = rayon_km*360/parallele_km
print("Angle sur méridien [deg] :", round(angle_meridien, 5))
print("Angle sur parallèle [deg] :", round(angle_parallele, 5))

k=pi/180
step = int(360/nb_points)
for theta in range(0, 360, step):
    lat_point = lat_centre_deg+angle_meridien*cos(theta*pi/180)
    long_point = long_centre_deg+angle_parallele*sin(theta*pi/180)
    milCalcRayon = (21600/(2*pi))*acos(sin(lat_centre_deg*k)*sin(lat_point*k) + cos(lat_centre_deg*k)*cos(lat_point*k)*cos((long_centre_deg*k) - (long_point*k)))

    print(theta, " ", round(lat_point, 5), " ", round(long_point, 5), "  Rayon Calculé=", round(milCalcRayon,5))

Copier-Coller et RUN à https://www.programiz.com/python-pro...line-compiler/


Résultats pour 12 points et R=0.5 milles

    0   47.17133   -2.42233   Rayon Calculé= 0.49998
  30   47.17022   -2.4162     Rayon Calculé= 0.49997
  60   47.16717   -2.41172   Rayon Calculé= 0.49996
  90   47.163       -2.41007   Rayon Calculé= 0.49998
120   47.15883   -2.41172   Rayon Calculé= 0.49999
150   47.15578   -2.4162     Rayon Calculé= 0.49999
180   47.15467   -2.42233   Rayon Calculé= 0.49998
210   47.15578   -2.42846   Rayon Calculé= 0.49999
240   47.15883   -2.43294   Rayon Calculé= 0.49999
270   47.163       -2.43459   Rayon Calculé= 0.49998
300   47.16717   -2.43294   Rayon Calculé= 0.49996
330   47.17022   -2.42846   Rayon Calculé= 0.49997

Résultats pour 12 points et R=120 milles (2 degrés)

    0   49.16291   -2.42233   Rayon Calculé= 119.9946
  30   48.89497   -0.95162   Rayon Calculé= 119.49957
  60   48.16296    0.12501   Rayon Calculé= 119.13601
  90   47.163       0.51908   Rayon Calculé= 119.98751
120   46.16304    0.12501   Rayon Calculé= 120.83011
150   45.43103   -0.95162   Rayon Calculé= 120.47762
180   45.16309   -2.42233   Rayon Calculé= 119.9946
210   45.43103   -3.89304   Rayon Calculé= 120.47762
240   46.16304   -4.96967   Rayon Calculé= 120.83011
270   47.163       -5.36374   Rayon Calculé= 119.98751
300   48.16296   -4.96967   Rayon Calculé= 119.13601
330   48.89497   -3.89304   Rayon Calculé= 119.49957

Résultats pour 12 points et R=600 milles (10 degrés)

   0   57.16255   -2.42233   Rayon Calculé= 599.973
  30   55.82286    4.9312    Rayon Calculé= 586.98273
  60   52.16278  10.31436   Rayon Calculé= 577.31138
  90   47.163      12.28473   Rayon Calculé= 599.08501
120   42.16322   10.31436   Rayon Calculé= 619.7486
150   38.50314    4.9312     Rayon Calculé= 611.46081
180   37.16345   -2.42233   Rayon Calculé= 599.973
210   38.50314   -9.77586   Rayon Calculé= 611.46081
240   42.16322  -15.15902   Rayon Calculé= 619.7486
270   47.163     -17.12939   Rayon Calculé= 599.08501
300   52.16278   -15.15902  Rayon Calculé= 577.31138
330   55.82286   -9.77586   Rayon Calculé= 586.98273

Une solution exacte en trigonométrie sphérique est possible ? Ca vous dit ?

Biname