nathan frad

Bonjour, Après avoir calculé les éphémérides des planètes selon le livre de Jean Meeus - "Calculs astronomiques à l'usage des amateurs" en utilisant la théorie VSOP87 avec sucées . Je souhaite calculer la position de la Lune ( éphémérides ), pour cela dans le livre ils utilisent la théorie ELP2000. Enfaîte mon problème , c'est que dans le livre pour le calcul de la position de la lune ils se limitent a une centaine de termes periodiqued ( les plus important selon eux ) de la théorie ELP2000 pour calculer sa longitude et latitude . Pour la longitude, il fait la somme des termes comme cela : Et moi du coup je voudrais faire la somme de la longitude mais avec tout les termes et non que les plus important. Mais dans la théorie elp200 , il y a plusieurs fichiers contenant le mot "longitude" dans le titre : Ci-dessous la liste des fichiers comprenant le mot "longitude" . ELP1 : Main problem. Longitude periodic terms (sine) ELP4 : Earth figure perturbations. Longitude ELP7 : Earth figure perturbations. Longitude/t ELP10 : Planetary perturbations. Table 1 Longitude ELP13 : Planetary perturbations. Table 1 Longitude/t ELP16 : Planetary perturbations. Table 2 Longitude ELP19 : Planetary perturbations. Table 2 Longitude/t ELP22 : Tidal effects. Longitude ELP25 : Tidal effects. Longitude/t ELP28 : Moon figure perturbations. Longitude ELP31 : Relativistic perturbations. Longitude ELP34 : Planetary perturbations - solar eccentricity. Longitude/t2 Il y a trop de fichier comportant le mot " longitude" , du coup je ne sais pas le qu’elle prendre pour faire la somme de tout les termes contenu dans les fichiers qui sont dans ce projet... Merci d'avance pour votre aide.

Ahah oui complétement ! Seul le résultat m'interrese dans mon cas ^^

Bonjour, Nous sommes trois étudiants qui sommes à la recherche d’aide dans le domaine de l’astronomie. Je m’explique : Depuis la rentrée 2017 nous travaillons dans une boite d’informatique en contrat de professionnalisation, afin d’avoir un diplôme de concepteur développeur informatique. Dans le cadre de ses études nous devons réaliser un projet sur deux ans. Celui-ci consiste à diriger un télescope Celestron via une raquette Nexstar. Lors de nos recherches nous nous sommes rendu compte qu’il fallait calculer les éphémérides des planètes, cependant nous n’avons pas les connaissances nécessaires dans ce domaine. C’est pourquoi nous sommes à la recherche de personnes qui souhaiteraient partager son savoir bénévolement. La difficultés que nous rencontrons est de comprendre le calcul des positions des planètes... Merci d'avance pour vos multiples réponses.

Holà tout le monde ! Voici le résultat des éphémérides des planètes ! : https://github.com/nathanfrad/PlanetaryEphemeris

Holà tout le monde ! Voici le résultat des éphémérides des planètes ! : https://github.com/nathanfrad/PlanetaryEphemeris

Merci beaucoup d'avoir pris le temps de répondre aussi précisément. Avec mon groupe d'étude nous sommes entrains de créer une application mobile qui permettra d'obtenir des clichés des planètes en temps réel. Mais du coup avec ce calcul, il y a une précision de 0.25 ° et cela correspond a un mouvement d'un peu près une demi-heure ( comme lu ci-dessus ) alors cette précision est largement suffisante pour mon projet ? Je ne sais pas de combien est la précision dans le livre de Jean Meeus. Je vais regarder. New: Je viens de regarder : pour le livre de Jean Meeus la précision est de 10 seconde d'arc sur la longitude et 4 seconde d'arc sur la latitude Donc je vais faire les calculs avec le livre de Meeus je pense. Merci

Oui tout a fait, c'est l'effet du Temps-lumière. Et pour la Lune j'ai tout les calculs, il me manque juste la somme des termes de la théorie ELP2000 . Merci d'avance . New: Je suis en entreprise jusqu'au 6 janvier (avant de retourner en formation) donc je vais continuer a rechercher des nouvelles informations. Et je vais peut être mettre un autre post sur le forum pour la théorie ELP2000.

Je viens de tout afficher correctement et cela est correcte pour toute les planètes, c’était un problème de variable non mis a jours. Merci beaucoup ! Exemple pour mars aujourd’hui : Julian day : 2458093.5606712964 t : 0.01792898198849127 Coordonnées écliptiques héliocentriques : Planète Mars : λ = 182.376633438691° β = 1.3600854103658517° ρ = 1.6488959513779677° Planète Terre : λo = 74.11666839850386° βo = 359.9999974415216° ρo = 0.9853179018538474° Coordonnées cartésiennes héliocentriques : Planète Mars : X = -1.9166746951200193 Y = -1.016056876472827 Z = 0.03913780831125706 Coordonnées écliptiques héliocentriques pour l'instant t-T (aberration) : Planète Mars : λ1 = 182.3710499309832° β1 = 1.360207558857363° ρ1 = 1.6489031342337546° Coordonnées cartésiennes héliocentriques pour l'instant t-T (aberration) : Planète Mars : X1 = -1.9166884400884336 Y1 = -1.0158966673340497 Z1 = 0.03914149309682227 Coordonnées écliptiques héliocentriques pour l'instant t-T : Planète Mars : λ1 = 207.9247984504529° β1 = 1.033710442737094° ρ1 = 1.6489031342337546° Coordonnées écliptiques héliocentriques pour l'instant t-T + abberation : Planète Mars : λ1 = 1.009780488577874° β1 = 1.0085960442814097° ρ1 = 1.6489031342337546° Coordonnées écliptiques héliocentriques pour l'instant t-T + abberation + nutation : Planète Mars : λ1 = 207.92138108076188° β1 = 1.0085960442814097° ρ1 = 1.6489031342337546° Coordonnées équatoriales moyennes pour l'équinoxe : : Planète Mars : α = 206.29608276909988° α heure = 13.759948720698961 h δ = -9.792715014994402°

Merci , je remet de l'ordre dans mes résultat afin de pouvoir comparer mes résultat avec les tiens et je reviens vers toi en modifiant ce message. C'est bien Mars et non venus, juste une erreur d'affichage .

J'ai commencé a chercher les problème d'atan2, mais je continuerais demain. C'est bizarre car les longitudes de venus et saturne sont supérieur a 180°. ( venus longitude en deg : 233.40442810862078 , saturn longitude en deg : 269.7095372598329 ) ... C'est très gentil de proposer de comparer nos résultat : voici mes résultat pour mars : Un grand merci d'avance !! Debut constructeur greenwichPosition -- Julian day : 2458092.7674421296 t : 0.01792681024539234 heliocentrique L (deg) : 182.02328635700997 heliocentrique B (deg): 1.3677899490808634 heliocentrique R : 1.6493478011978169 heliocentrique terre Lo (deg) : 73.31189105773137 heliocentrique terre Bo (deg): 2.2184043058413104E-5 heliocentrique terre Ro : 0.9854272119530312 x : -1.9308268537879374 y : -1.0021378363422846 z : 0.03936983200723439 Vrai distance : 2.175757929092687 temps lumiere : 0.012566129735844917 Coordonnées heliocentrique de venus pour l'instant t-T en degres heliocentrique L1 (deg pour t-T): 182.01769026187776 heliocentrique B1 (deg pour t-T): 1.3679115509958741 heliocentrique R1 : 0.7242975917764743 x1 : -1.0066192270923175 y2 : -0.969416945334475 z3 : 0.01729026238490708 Vrai distance1 : 1.3976231378404187 temps lumiere1 : 0.00807199800910076 longitude en deg: 223.92143576175135 latitude en deg: 0.7088350970353545 l'exentricité de le l'orbite de la terre (e) : 0.016701094035012586 longitude du perihelie ( pi deg ) : 103.24561090717062 piRad : 1.8019758485630968 long geo moyenne du soleil ( LO) : 6734.256164048508 anomalie moyenne du soleil ( M) : 6811.010540692197 centre du soleil : 0.07360506468399858 LongVraiSoleil : 254.32976911319201 LongVraiSoleilRad : 4.438891856862179 K en radian: 9.936508497454117E-5 Aberration_Longitude_deg : -0.026625793469303344 Aberration_Latitude_deg : -0.014763885925450427 nutation_longitude en deg : -0.0034513864105610603 nutation_latitude : -7.477089645817657 Longitude_nutation : 223.91798437534078 rad en deg : 221.67866015827406 rad en heure : 14.785966632556878 dec en deg : -15.339496520557677 Fin constructeur greenwichPosition -- ******************************************* pi : 2.5849317577312025E-5 TSM :170.52953186631203 long geo moyenne du soleil ( L0) : 1.5764727194520578E9 anomalie moyenne du soleil ( M) : 118.87455932311805 nutation_longitude en deg : -0.0030281514152984018 nutation_latitude : -6.507131460476364 ObliquitM : 23.003626543129922 ObliquitVrai : 16.49649508265356 TSM : 170.51844129819844 TSM heure : 11.367896086546562 longitude : -0.5805 longitude heure : -0.0387 H_rad : -0.9030459014570468 u degres : 44.744230177359114 ra_position : 221.67951667644542 ra_position en heure : 14.786023762318909 dec_position : -15.33864000238631

Salut, Salut ! Bon bah me voila de retour . Après avoir testé pour chaque planètes leurs positions, je me rend compte que mon calcul me sort les bonnes positions ( RA , DEC ) que pour Saturne et Venus ^^... Bizarre le calcul étant pourtant universel " normalement" ... Je n'ai plus qu'a chercher l'erreur dans tout les calculs... Let's go !

Enfaîte mon problème , c'est que dans le livre pour le calcul de la position de la lune il se limite aux termes periodique les plus important de la théorie ELP2000 pour calculer sa longitude et latitude . Pour la longitude, il fait la somme des termes cela donne ça : Et moi du coup je voudrais faire la somme de la longitude ( par example ) mais avec tout les termes et non que les plus important. Et dans la theorie elp200 il y a plusieurs fichier contenant des termes longitude : ELP1 : Main problem. Longitude periodic terms (sine) ELP4 : Earth figure perturbations. Longitude ELP7 : Earth figure perturbations. Longitude/t ELP10 : Planetary perturbations. Table 1 Longitude ELP13 : Planetary perturbations. Table 1 Longitude/t ELP16 : Planetary perturbations. Table 2 Longitude ELP19 : Planetary perturbations. Table 2 Longitude/t ELP22 : Tidal effects. Longitude ELP25 : Tidal effects. Longitude/t ELP28 : Moon figure perturbations. Longitude ELP31 : Relativistic perturbations. Longitude ELP34 : Planetary perturbations - solar eccentricity. Longitude/t2 Il y a trop de fichier comportant des " longitude" , du coup je ne sais pas le qu’elle prendre pour faire la somme de tout les termes qui sont dans ce projet...

Ah super ! J'ai trouvé une autre source : http://cdsarc.u-strasbg.fr/viz-bin/Cat?cat=VI%2F79&target=http& Dans celui ci il y a un fichier exemple tout en bas . Et il y aussi un fichier nommé elp82b.ps.gz ( vers le bas ) ( je te le met aussi en fichier joint), c’est un pdf explicatif mais j'avoue que je ne comprend pas bien tout ^^ elp82b.ps.gz Modification : J'ai trouvé le calcul du livre de Meeus sur internet. Les terme sont presque identique. ( http://www.stargazing.net/kepler/ch30.txt )

Je viens de me rendre contre que dans la théorie vsop87, il n'y a pas la Lune... Du coup je vais calculer la position de la lune avec la méthode décrite dans le livre de Jean MEEUS. Dans le livre, il utilise la théorie ELP2000 de M.Chapront et H.Chapront. Ce pendant, pour l'instant, je ne retrouve pas sur internet, les termes de la théorie ELP2000 présente dans le livre ... Exemple de source trouvé : ftp://cdsarc.u-strasbg.fr/pub/cats/VI/79/ ftp://ftp.imcce.fr/pub/ephem/moon/elp82b/ http://astronomie.baillet.org/logiciels/heclium/telechargement.php .... Si vous avez des informations sur le sujet n’hésiter pas. Merci d'avance

Oui c'est bien ça, dans mon cas aussi je les ai additionné.

oh la la ! Merci beaucoup! C'était bien ça l'erreur. J'ai remplacé double longitude = Math.atan(tan_longitude); par double longitude = atan2(Y1, X1); Et voici que le programme marche pour tout les JJ ! Merci grandement pour votre aide. J'optimise le code puis l'adapte pour toute les planètes et je le mettrais ensuite en Open Source sur github. Je vous enverrai le liens github quand je l'aurais fait . Et je vais devoir l'adapter pour le lieux du télescope, car pour l’instant je calcul le temps sidéral à Greenwich mais le télescope ne se trouve pas forcément à Greenwich ... Je fais des recherches et des tests à ce sujet et je reviens vers vous si besoin. Encore Merci ! PS : J'essayerai plus tard de calculer la position des étoiles mais c'est pas gagné

Bonjour , Je reviens vers vous que maintenant car j’étais un peu perdu du coup j'ai décidé me relancer dans les calculs des éphémérides depuis le début (en utilisant la théorie vsop87 toujours ) et en suivant le livre je Jean Meeus - Calcul astronomiques à l'usage des amateurs et en suivant les conseils ci-dessus de bruno ... Afin de réaliser les calculs et de pouvoir les vérifier, j'ai réalisé le même calcul que dans leurs exemple. C'est a dire que j'ai calculé l' ascension droite et la déclinaison ( la position apparente comme dit dans le livre ) de Venus pour le 20 décembre 1992 à 0h Temps dynamique = JJ 2448976.5 Après deux jours de calcul et de test j'arrive au bon résultat, c'est a dire : rad en deg : 316.1767716972735 dec en deg : -18.88643717742539 Cependant mon calcul marche pour certaine date mais pas toutes ... et ça je ne comprend pas pourquoi ... Je vous met mon calcul entiers et commenté avec les résultat à la fin . Merci d'avance aux personnes qui prendrons le temps de se pencher sur le sujet, je leurs suis très reconnaissant ! Et pour les plus courageux je leurs met le fichier java en fichier joint si vous voulez tester. ( il suffit juste de l’exécuté dans votre éditeur ) fichier : CompletVenusNow.java Merci beaucoup !!! Le code / Le calcul est le suivant : package ephemeride; public class CompletVenusNow { public static void main(String[] args) { CompletVenusNow tt = new CompletVenusNow(); double tVenus; tVenus = (( 2448976.5 - 2451545) / 365250); /* ******** Coordonnées heliocentrique de venus en deg ********** */ double L = (tt.Venus_L0(tVenus)+ tt.Venus_L1(tVenus) + tt.Venus_L2(tVenus) + tt.Venus_L3(tVenus) + tt.Venus_L4(tVenus)+ tt.Venus_L5(tVenus))*(180/Math.PI); if ( L < -360 ) { L = L+ (360 * Math.floor(L / 360)*-1) ; } else if (L > 0 ) { L = L- (360 * Math.floor(L / 360)) ; } System.out.println("\n heliocentrique L (deg) : " + L ); double B = (tt.Venus_B0(tVenus)+ tt.Venus_B1(tVenus) + tt.Venus_B2(tVenus) + tt.Venus_B3(tVenus) + tt.Venus_B4(tVenus)+ tt.Venus_B5(tVenus))*(180/Math.PI); if ( B < -360 ) { B = B+ (360 * Math.floor(B / 360)*-1) ; } else if (L > 360 ) { B = B- (360 * Math.floor(B / 360)) ; } System.out.println(" heliocentrique B (deg): " + B); double R = (tt.Venus_R0(tVenus)+ tt.Venus_R1(tVenus) + tt.Venus_R2(tVenus) + tt.Venus_R3(tVenus) + tt.Venus_R4(tVenus)+ tt.Venus_R5(tVenus)); System.out.println(" heliocentrique R : " + R); /* ******** Coordonnées heliocentrique de venus en rad ********** */ L = L *Math.PI / 180; // 0,45378560552 rad B = B *Math.PI / 180; // -0,03490658504 rad /* ******** Coordonnées heliocentrique de la terre en deg ********** */ double Lo = (tt.Earth_L0(tVenus)+ tt.Earth_L1(tVenus) + tt.Earth_L2(tVenus) + tt.Earth_L3(tVenus) + tt.Earth_L4(tVenus)+ tt.Earth_L5(tVenus))*(180/Math.PI); if ( Lo < -360 ) { Lo = Lo+ (360 * Math.floor(Lo / 360)*-1) ; } else if (Lo > 0 ) { Lo = Lo- (360 * Math.floor(Lo / 360)) ; } System.out.println("\n heliocentrique terre Lo (deg) : " + Lo ); double Bo = (tt.Earth_B0(tVenus)+ tt.Earth_B1(tVenus) + tt.Earth_B2(tVenus) + tt.Earth_B3(tVenus) + tt.Earth_B4(tVenus)+ tt.Earth_B5(tVenus))*(180/Math.PI); if ( Bo < -360 ) { Bo = Bo+ (360 * Math.floor(Bo / 360)*-1) ; } else if (Bo > 360 ) { Bo = Bo- (360 * Math.floor(Bo / 360)) ; } System.out.println(" heliocentrique terre Bo (deg): " + Bo); double Ro = (tt.Earth_R0(tVenus)+ tt.Earth_R1(tVenus) + tt.Earth_R2(tVenus) + tt.Earth_R3(tVenus) + tt.Earth_R4(tVenus)+ tt.Earth_R5(tVenus)); System.out.println(" heliocentrique terre Ro : " + Ro); /* ******** Coordonnées heliocentrique de la terre en rad ********** */ Lo = Lo *Math.PI / 180; // 1,53588974175 Bo = Bo *Math.PI / 180; Ro = Ro; /* ******************** Calcul ********************/ double X = (R * Math.cos(B) * Math.cos(L)) - (Ro * Math.cos(Bo) * Math.cos(Lo)); double Y = (R * Math.cos(B) * Math.sin(L)) - (Ro * Math.cos(Bo) * Math.sin(Lo)); double Z = (R * Math.sin(B)) - (Ro * Math.sin(Bo)); System.out.println("x : " + X); System.out.println("y : " + Y); System.out.println("z : " + Z); /* ******** Vrai distance terre - venus ********** */ double VraiDis = Math.sqrt((X * X) + (Y * Y) + (Z * Z)); System.out.println("Vrai distance : " + VraiDis); /* ******** Temps lumière ( T) ********** */ double TempLum = 0.0057755183 * VraiDis; System.out.println("temps lumiere : " + TempLum); /* ******** Coordonnées heliocentrique de venus pour l'instant t-T en degres ********** */ tVenus = (( (2448976.5 - TempLum ) - 2451545) / 365250); double L1 = (tt.Venus_L0(tVenus)+ tt.Venus_L1(tVenus) + tt.Venus_L2(tVenus) + tt.Venus_L3(tVenus) + tt.Venus_L4(tVenus)+ tt.Venus_L5(tVenus))*(180/Math.PI); if ( L1 < -360 ) { L1 = L1+ (360 * Math.floor(L1 / 360)*-1) ; } else if (L1 > 0 ) { L1 = L1- (360 * Math.floor(L1 / 360)) ; } System.out.println("\n heliocentrique L1 (deg pour t-T): " + L1 ); double B1 = (tt.Venus_B0(tVenus)+ tt.Venus_B1(tVenus) + tt.Venus_B2(tVenus) + tt.Venus_B3(tVenus) + tt.Venus_B4(tVenus)+ tt.Venus_B5(tVenus))*(180/Math.PI); if ( B1 < -360 ) { B1 = B1+ (360 * Math.floor(B1 / 360)*-1) ; } else if (L > 360 ) { B1 = B1- (360 * Math.floor(B1 / 360)) ; } System.out.println(" heliocentrique B1 (deg pour t-T): " + B1); double R1 = (tt.Venus_R0(tVenus)+ tt.Venus_R1(tVenus) + tt.Venus_R2(tVenus) + tt.Venus_R3(tVenus) + tt.Venus_R4(tVenus)+ tt.Venus_R5(tVenus)); System.out.println(" heliocentrique R1 : " + R1); /* ******** Coordonnées heliocentrique de venus pour l'instant t-T en rad ********** */ L1 = L1 *Math.PI / 180; // 0,45378560552 rad B1 = B1 *Math.PI / 180; // -0,03490658504 rad double X1 = (R1 * Math.cos(B1) * Math.cos(L1)) - (Ro * Math.cos(Bo) * Math.cos(Lo)); double Y1 = (R1 * Math.cos(B1) * Math.sin(L1)) - (Ro * Math.cos(Bo) * Math.sin(Lo)); double Z1 = (R1 * Math.sin(B1)) - (Ro * Math.sin(Bo)); System.out.println("\nx1 : " + X1); System.out.println("y2 : " + Y1); System.out.println("z3 : " + Z1); /* ******** Vrai distance terre - venus ( avec t-T )********** */ double VraiDis1 = Math.sqrt((X1 * X1) + (Y1 * Y1) + (Z1 * Z1)); System.out.println("Vrai distance1 : " + VraiDis1); /* ******** Temps lumière (T) ( avec t-T ) ********** */ double TempLum1 = 0.0057755183 * VraiDis1; System.out.println("temps lumiere1 : " + TempLum1); /* ******** calcul longitude de venus ********** */ double tan_longitude = Y1 / X1; double longitude = Math.atan(tan_longitude); double longitude_deg = (longitude * (180 / Math.PI)) + 360; System.out.println("\nlongitude en deg: " + longitude_deg ); /* ******** calcul latitude de venus ********** */ double tan_latitude = Z1 / (Math.sqrt((X1 * X1) + (Y1 * Y1))); double latitude = Math.atan(tan_latitude); double latitude_deg = latitude * (180 / Math.PI); System.out.println("latitude en deg: " + latitude_deg); double t = (((2448976.5 ) - 2451545) / 36525); /* ******** l'exentricité de le l'orbite de la terre (e) ********** */ double e = 0.016708634 - 0.000042037 * t - 0.0000001267* (t*t) ; System.out.println("l'exentricité de le l'orbite de la terre (e) : " + e); /* ******** longitude du perihelie ( pi ) ********** */ double pi = 102.93735 + 1.71947*t + 0.00046* (t*t) ; System.out.println("longitude du perihelie ( pi deg ) : " + pi); double piRad = pi * Math.PI / 180; System.out.println("piRad : " + piRad); /* ******** longitude geocentrique moyenne du soleil ( LO) ********** */ double L0 = 280.46646 + 36000.76983 * t + 0.0003032 * (t * t); System.out.println("\nlong geo moyenne du soleil ( LO) : " + L0); /* ******** anomalie moyenne du soleil (M) ********** */ double M = 357.52911 * Math.PI / 180 + (35999.05029 * Math.PI) / 180 * t - 0.0001537 * Math.PI / 180 * (t * t); System.out.println("anomalie moyenne du soleil ( M) : " + M); /* ******** centre du soleil (C) ********** */ double C = (1.914602 - 0.004817 * t - 0.000014 * (t * t)) * Math.sin(M) + (0.019993 - 0.000101 * t) * Math.sin(2 * M) + 0.000289 * Math.sin(3 * M); System.out.println("centre du soleil : " + C); /* ******** Longitude Vrai du Soleil ********** */ double LongVraiSoleil = (L0 + C) - 360 * Math.floor((L0 + C) / 360); System.out.println("LongVraiSoleil : " + LongVraiSoleil); double LongVraiSoleilRad = LongVraiSoleil * Math.PI / 180; System.out.println("LongVraiSoleilRad : " + LongVraiSoleilRad); /* ******** constante de l'abberation ********** */ double K = ((20.49552 / 3600)*Math.PI)/180 ; System.out.println("K : " + K); /* ******************* **************************** Calcul Aberration *********************************************************** */ double Aberration_Longitude = (-K * Math.cos(LongVraiSoleilRad - longitude) + e * K * Math.cos(piRad - longitude)) / Math.cos(latitude) ; double Aberration_Longitude_deg = Aberration_Longitude* ( 180 / Math.PI) ; System.out.println("Aberration_Longitude_deg : " + Aberration_Longitude_deg ); double Aberration_Latitude = -K * Math.sin(latitude)* (Math.sin(LongVraiSoleilRad - longitude) - e * Math.sin(piRad - longitude)); double Aberration_Latitude_deg = (Aberration_Latitude * ( 180 / Math.PI)) ; System.out.println("Aberration_Latitude_deg : " + Aberration_Latitude_deg); /* ******************* **************************** Longitude et latitude moins l'Aberration *********************************************************** */ double latitude_aberration = latitude_deg - Aberration_Longitude_deg ; double longitude_aberration = latitude_deg - Aberration_Latitude_deg ; /* *********** LongitudeMoyenneSoleil = L0 ***************** */ L0 = L0 * (Math.PI / 180 ) ; /* *********** LongitudeMoyenneLune ***************** */ double LLune = (218.3164 + 481267.8812 * t + 0.001599* (t*t)) * (Math.PI / 180 ); /* *********** anomalie moyenne du soleil = M ***************** */ M = M * (Math.PI / 180 ); /* *********** anomalie moyenne de la lune ***************** */ double MLune = (134.9634 + 477198.8675 * t + 0.008721 * (t*t)) * (Math.PI / 180 ) ; /* *********** anomalie moyenne de la lune ***************** */ double omega = (125.0443 - 1934.1363 * t + 0.002075 * (t*t)) * (Math.PI / 180 ) ; /* *********** nutation de la longitude ***************** */ double nutation_longitude = 0 ; nutation_longitude += -(17.1996 + 0.01742 * t)* Math.sin(omega); nutation_longitude += -(1.3187 + 0.00016 *t) * Math.sin(2*L0) ; nutation_longitude += -0.2274 * Math.sin(2*LLune); nutation_longitude += +0.2062 * Math.sin(2*omega) ; nutation_longitude += +(0.1426 - 0.00034 *t) * Math.sin(M); nutation_longitude += +0.0712 * Math.sin(MLune); nutation_longitude += -(0.0517 - 0.00012 * t) * Math.sin(2 * L0 + M); nutation_longitude +=-0.0386 * Math.sin(2*LLune-omega); nutation_longitude += -0.0301 * Math.sin(2*LLune+MLune); nutation_longitude += +0.0217 * Math.sin(2*L0-M) ; nutation_longitude += -0.0158 * Math.sin(2*L0-2*LLune + MLune) ; nutation_longitude += +0.0129 * Math.sin(2*L0 - omega) ; nutation_longitude += +0.0123 * Math.sin(2*LLune - MLune) ; nutation_longitude = nutation_longitude / 3600; System.out.println("\nnutation_longitude en deg : " + nutation_longitude); // 0,00444444444 /* *********** nutation de la latitude ***************** */ double nutation_latitude = 0 ; nutation_latitude += +(9.2025 + 0.00089*t) * Math.cos(omega) ; nutation_latitude += +(0.5736 - 0.00031* t) * Math.cos(2 * L0) ; nutation_latitude += +0.0977 * Math.cos(2 * LLune) ; nutation_latitude += -0.0895 * Math.cos(2* omega) ; nutation_latitude += +0.0224 * Math.cos(2*L0 + M); nutation_latitude += + 0.0200 * Math.cos(2*LLune-omega) ; nutation_latitude += +0.0129 * Math.cos(2*LLune + MLune) ; nutation_latitude += -0.0095 * Math.cos(2*L0 - M); nutation_latitude += -0.0070 * Math.cos(2 * L0 - omega); System.out.println("nutation_latitude : " + nutation_latitude); // deg -0,016666667 /* *********** obliquité moyenne ( qui est non utilisé dans mon cas ) ***************** */ double ObliquitM= (23+ 26/60 + 21.448/3600 )- ( 46.8150 / 3600 ) *t - (0.00059/3600) *(t*t) + (0.001813/3600) *(t*t*t) ; System.out.println("ObliquitM : " + ObliquitM); /* *********** longitude corrigé pour la nutation ***************** */ double Longitude_nutation = longitude_deg + nutation_longitude ; System.out.println( "Longitude_nutation : " + Longitude_nutation); /* *********** obliquité de l'écliptique ***************** */ double E = (23.4392911 * Math.PI / 180); /* *********** Conversion de Longitude_nutation en radian ***************** */ Longitude_nutation = Longitude_nutation * Math.PI / 180 ; /* *********** Calcul de l'Ascension droite ***************** */ double F; F = Math.cos(Longitude_nutation) * Math.cos(latitude); double G; G = Math.cos(E) * Math.sin(Longitude_nutation) * Math.cos(latitude) - Math.sin(E) * Math.sin(latitude); double ra_rad; ra_rad = Math.atan2(G, F); // dans cette ordre G,F System.out.println("rad en deg : "+ (ra_rad* ( 180 / Math.PI)+360 )); ; /* *********** Calcul de la Déclinaison ***************** */ double sin_dec; sin_dec = ((Math.sin(E)) * Math.sin(Longitude_nutation) * Math.cos(latitude) + (Math.cos(E) * Math.sin(latitude))); double dec_rad; dec_rad = Math.asin(sin_dec); System.out.println("dec en deg : " + dec_rad * ( 180 / Math.PI)); } Le résultat est le suivant : heliocentrique L (deg) : 26.11411959514089 heliocentrique B (deg): -2.6206031211961918 heliocentrique R : 0.7246016759555223 heliocentrique terre Lo (deg) : 88.3570045917113 heliocentrique terre Bo (deg): 1.662695729726217E-4 heliocentrique terre Ro : 0.9838241642711414 x : 0.6217453110684215 y : -0.6648122386978192 z : -0.03313324305849642 Vrai distance : 0.9108459564347573 temps lumiere : 0.005260607489869944 heliocentrique L1 (deg pour t-T): 26.10571812103899 heliocentrique B1 (deg pour t-T): -2.620919380428602 heliocentrique R1 : 0.7246023679452037 x1 : 0.6217924790901884 y2 : -0.6649073232679508 z3 : -0.0331372701491552 Vrai distance1 : 0.9109477011669834 temps lumiere1 : 0.005261195118432844 longitude en deg: 313.0808463015071 latitude en deg: -2.0846912526550136 l'exentricité de le l'orbite de la terre (e) : 0.01671158948664706 longitude du perihelie ( pi deg ) : 102.81643622562049 piRad : 1.794485337303849 long geo moyenne du soleil ( LO) : -2251.1687830825636 anomalie moyenne du soleil ( M) : -37.94319988043709 centre du soleil : -0.47236015117494784 LongVraiSoleil : 268.3588567662614 LongVraiSoleilRad : 4.683745627459124 K : 9.936508497454117E-5 Aberration_Longitude_deg : -0.004130095207340988 Aberration_Latitude_deg : -1.4747370015569854E-4 nutation_longitude en deg : 0.004629498315663484 nutation_latitude : -1.9219445605245453 ObliquitM : 23.006872251866344 Longitude_nutation : 313.08547579982275 rad en deg : 316.1767716972735 dec en deg : -18.88643717742539

En faite je pense que mon problème vient de mon calcul pour la longitude et latitude : pour JD : 2451545 Exemple pour longitude : JD_t = ((JD - 2451545) / 365250); JD_t = ((2451545 - 2451545) / 365250); JD_t = 0.017876814586660537; Après en utilisant les méthodes du programme fourni par vsop87 ( Mercury_L0(JD_t) ) cela me donne : L0 = 4.4293481036110025 L1 =0.0 L2=0.0 L3 =0.0 L4 =-0.0 L5 =-0.0 Du conformément a la règle j’additionne tout les L ( ici ce n'a pas grand intérêt car les autres valeurs sont null): cela donne Longitude in radians : 4.035154791673375 Longitude in degrés : 231.19733924487534 Alors que ça devrais me donner λ = 351,033 si j'avais le bon calcul ... J'ai vue sur internet : Il faut ajouter 180° à λ si x est négatif. Mais au pire des cas ça donne 411,197339245 degrés dans mon cas. Ce qui n'est pas bon non plus ...

C'est vrai que j'ai arrondi très rapidement disons ^^ Super ça fonctionne d'enfer ! Mais j'ai encore quelque interrogation... Exemple pour le 1 janvier 2000 a 12 h 00 .00 pour Mercure ( JD : 2451545) : Le convertisseur en ligne donne : (http://pgj.pagesperso-orange.fr/position-planetes.htm) RA= 18h08m19.7s = 18,138610778 heures = 272.08 Degrés DEC = -24°25'11.34'' Résultat de mon calcul avec vsop87 mercure version D donne : Mercury Longitude in radians : 4.045026248245587 Longitude in degres : 231.76293204410973 Mercury Latitude in radians : -0.12081520741805907 Latitude in degres : -6.92220148605242 RA_rad : -2.3124504136790285 DEC_rad : -0.4286088586741546 RA_deg : 227.50635096291037 DEC_deg : -24.557478663948224 On remarque que le RA calculé ( 227.5 ) et différent du RA donné en ligne ( 272.08 ) D’où vient cette différence... ?? Voir mes calculs en fin de message ... J'ai l'impression que les longitudes et latitudes donné par le convertisseur en ligne sont différentes de ce donné par vsop87 Mercure D ... Je me suis dit je vais essayé de calcul les coordonnées de Venus avec la même technique que Mercure mais avec le vsop87 de venus version D bien sur, mais la le résultat et encore plus différent que les résultats donné par le convertisseur en ligne .... Exemple pour le 1 janvier 2000 a 12 h 00 .00 pour Venus ( JD : 2451545) : Le convertisseur en ligne donne : (http://pgj.pagesperso-orange.fr/position-planetes.htm) RA= 15h59m34.77s =15,992777444 heures = 239.89 Degrés DEC = -18°26'59.96'' Résultat de mon calcul avec vsop87 Venus version D donne : Venus Longitude in radians : 3.187022183287107 Longitude in degres : 182.60292031692032 Venus Latitude in radians : 0.05697828490452141 Latitude in degres : 3.2646152489230458 dec_deg : 1.98882550755495 ra_deg : 183.67090718660307 Ici on remarque bien la différence entre les deux résultat ..... ( DEC : 239.89 vs 1.98882550755495 ) Après du coup je me suis mis a chercher dans le livre calculs astronomiques de Jean MEEUS et j'ai trouver ce calcul pour obtenir la longitude héliocentriques Ecliptique a l'aide de vsop87 : L = (tt.Venus_L0(t)+ tt.Venus_L1(t)*(t) + tt.Venus_L2(t)*Math.pow(t, 2) + tt.Venus_L3(t)*Math.pow(t, 3) + tt.Venus_L4(t)*Math.pow(t, 4)+ tt.Venus_L5(t)*Math.pow(t, 5))/Math.pow(10, 8); Ce qui me donne : 3.1870221832871074E-8 Cela ne m'a pas l'aire bon non plus malheureusement ...... Voila, je vous remercie d'avance pour votre aide et j'espère avoir explicité clairement mes problèmes. ****************************************************************************************************************************************************************************************************************************------------------ CALCUL RA ET DEC -------------------------------------****************************************************************************** // calcul de la Acension Droite ( RA) en radian static double RA_rad(double e, double longitude, double latitude) { double X; X = Math.cos(longitude) * Math.cos(latitude); double Y; Y = Math.cos(e) * Math.sin(longitude) * Math.cos(latitude) - Math.sin(e) * Math.sin(latitude); double ra_rad; ra_rad = Math.atan2(Y, X); // dans cette ordre X,Y return ra_rad ; } // calcul de la Declinaison en radian static double DEC_rad(double e, double longitude, double latitude) { double sin_dec; sin_dec = ((Math.sin(e)) * Math.sin(longitude) * Math.cos(latitude) + (Math.cos(e) * Math.sin(latitude))); double dec_rad; dec_rad = Math.asin(sin_dec); return dec_rad;

Merci beaucoup Bruno ! J'ai réussi a calculer la déclinaison (DEC) avec votre aide mais maintenant je n’obtient pas le bon RA ( Right acccension ) .... voici mon calcul ci-dessous ... je pense avoir fait attention aux unités pourtant... // calcul RA en utilisant la deuxième équation *************************************************** cos_ra = (Math.cos(longitude)* Math.cos(latitude)/ Math.cos(dec_rad); cos_ra = (Math.cos(4.4293481036110025)* Math.cos(-0.05275734092008231)/ Math.cos(-0.444270950257026); cos_ra = -0.30887225220868564 ra_rad = Math.acos(cos_ra); ra_rad = Math.acos(-0.30887225220868564); ra_rad = 1.8848034049959312 ra_deg = ra_rad*180/Math.PI; ra_deg = 1.8848034049959312 *180/Math.PI; ra_deg = 107.99128031815367 Ce qui donne 107 alors que ça devrais donner 252 ..... ( source de 252 : nasa )

Bonjour, je reviens vers vous pour plus ample information. Je viens de tester le code généré par vsop87 ( http://neoprogrammics.com/vsop87/) qui me renvoie soit : L = Longitude in radians B = Latitude in radians R = Radius vector in AU (Astronomical Units) ( qui sont les coordonnées écliptiques si j'ai bien compris ) soit x, y , z . ( HELIOCENTRIC DYNAMICAL ECLIPTIC AND EQUINOX OF THE DATE Rectangular (X,Y,Z) Coordinates in AU (Astronomical Units) Toutes ces informations renvoyé par le programme sont correct, j'ai vérifié avec des logiciels de calcul des éphéméride en ligne. Cependant nous devons envoyé au télescope les donnes sous forme RA DEC ( coordonnées équatoriales). J'ai cherché comment convertir les coordonnées écliptiques vers les coordonnées équatoriales . J'ai trouvé ce calcul : ( https://fr.wikipedia.org/wiki/Système_de_coordonnées_célestes#Des_coordonn.C3.A9es_horizontales_aux_coordonn.C3.A9es_horaires) ce qui me donne ce calcul ( j'ai essayé juste pour DEC pour l'instant mais le résultat n'est pas bon comparé au résultat donné par les logiciels de calcul des éphémérides en ligne ) : dec = (Math.sin(23.439281)* Math.sin(longitude)* Math.cos(latitude)+ (Math.cos(23.439281)* Math.sin(latitude); exemple pour le 1 janvier 2000 a 12 h 00 .00 : dec = (Math.sin(23.439281)* Math.sin(4.4293481036110025)* Math.cos(-0.05275734092008231)+ (Math.cos(23.439281)* Math.sin(-0.05275734092008231); dec = 0.9581220776863877 dec2 =Math.asin(dec) 1.2803706926523695 Mais devrais donné : 1.711989 selon le site de la nasa ( https://heasarc.gsfc.nasa.gov/cgi-bin/Tools/convcoord/convcoord.pl?CoordVal=4.4293481036110025+%2C+-0.05275734092008231&CoordType=Ecliptic&Resolver=GRB%2FSIMBAD%2BSesame%2FNED&NoCache=on&Epoch=2000 ) Si vous avez une idée pour que je trouve les coordonnées équatorial des planètes, je vous en serait très reconnaissant. Merci d'avance

Merci, sur vos conseils nous avons commandé ce livre. Ce lien est très intéressant. Nous allons nous penché sur le sujet et nous vous tenons au courant. Merci beaucoup. Tout vos messages sont très important pour l'avancer de notre projet.

merci de votre réponse. Cependant nous savions déjà cela mais nous voulons l’adapter a un système embarquer! C'est pour cela que nous avons besoin du code en "dur" . Nous avons déjà examiner les logiciels que vous avez cité mais sans aboutissement. Nous souhaitons juste le calcul mathématique , même si il est très complexe. C'est pour cela que nous souhaitons un partage de connaissance dans ce domaine. Nous sommes de plus très curieux dans ce domaine.

Nous pensions aussi, au début, pouvoir utiliser la basse de donnes de la raquette en envoyant une commande externe ( ex : "venus") mais ce n'est pas le cas. Nous pouvons envoyer que des commandes en hexadécimale ( ra/dec converti en hexadecimal ).