Bonjour à tous,
Grâce à JCR, j'ai découvert le RGC (
http://professionnels.ign.fr/ficheProduitCMS.do?idDoc=5323862 )
J'ai téléchargé les 2 fichiers :
RGC® sur la France métropolitaine
http://professionnels.ign.fr/DISPLAY/000/526/741/5267412/RGC.zip
et
RGC® sur les DOM
http://professionnels.ign.fr/DISPLAY/000/526/741/5267413/RGC_DOM.zip
La spécification du contenu est dans ce doc PDF :
http://professionnels.ign.fr/DISPLAY/000/526/601/5266016/DC_RGC.pdf.
J'aimerai faire de l'orthodromie entre la ville de Rennes ( INSEE 35238 ) et
Saint Denis à la réunion ( INSEE 974011 ).
Formule :
Ortho(A,B)=6371 x acos[cos(LatA) x cos(LatB) x cos(LongB-LongA)+sin(LatA) x
sin(LatB)]
Source :
http://www.lion1906.com/Pages/francais/infos/orthodromie_loxodromie_navigation_cap.html
Ci dessous l'extrait complet du RGC pour ces deux points.
ADMI --- 2
ARRD --- 3
CANT --- 98
CARTE --- 1218O
COM --- 238
DEP --- 35
LATI_DMS --- 480651
LATI_GRD --- 53460
LONGI_DMS --- -14051
LONGI_GRD --- -4463
NOM --- RENNES
POPU --- 2062
SURFACE --- 5039
XLAMB2 --- 3010
XLAMB93 --- 3519
XLAMBZ --- 3010
YLAMB2 --- 23538
YLAMB93 --- 67896
YLAMBZ --- 23538
ZCHL --- 39
ZMAX --- 74
ZMIN --- 20
---------------------------------------------
DEP --- 974
COM --- 11
ARRD --- 1
CANT --- 93
ADMI --- 3
POPU --- 1316
SURFACE --- 14279
NOM --- SAINT-DENIS
XPROJ_LOC --- 1518
YPROJ_LOC --- 777
LONGI_DMS_LOC --- +552715
LATI_DMS_LOC --- -205159
XPROJ_MOND --- 3386
YPROJ_MOND --- 76905
LONGI_DMS_MOND --- +552653
LATI_DMS_MOND --- -205244
ZMIN --- 0
ZMAX --- 2276
ZCHL --- 23
CARTE --- 4402RT
Comment feriez vous?
sinon tu as tout ce qu'il faut dans ce que tu donnes (les lati, les
longi, et la formule).
"François Meignien" <preno...@yahoo.fr> a écrit dans le message de news:
49246f66$0$13869$426a...@news.free.fr...
Bonjour,
Voici un calcul fait avec Mathematica,
que vous pourrez peut-être adapter:
In[1]:= rad[{deg_, mn_, sec_}] =
(Pi*(3600*deg + 60*mn + sec))/648000;
In[2]:= Ortho[{LatA_,LongA_},{LatB_,LongB_}]:=
(20000/Pi)*ArcCos[Cos[LatA]*Cos[LatB]*Cos[LongA-LongB]+
Sin[LatA]*Sin[LatB]];
Vérification du degré de latitude:
In[3]:= Ortho[{rad[{45,0,0}],rad[{01,0,0}]},
{rad[{46,0,0}],rad[{01,0,0}]}]//N
Out[3]= 111.111
Votre calcul:
In[4]:= Ortho[{rad[{48,06,51}],-rad[{01,40,51}]},
{-rad[{20,52,44}],rad[{55,26,53}]}]//N
Out[4]= 9533.29
V.Astanoff
p.s. Je suis bien conscient de l'anomalie que
présente mon rayon de la terre, mais j'ai la
satisfaction incomparable de me conformer
aux spécifications de la glorieuse Révolution Française !
Bonjour,