nathan frad

Membre
 Afficher le profil  Parcourir son contenu

  • Compteur de contenus

    24

  • Inscription

  • Dernière visite

  • Last Connexion

    Soon available - 79390

Messages posté(e)s par nathan frad

  4. Posté(e) · Modifié par nathan frad

    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

  5. Posté(e)

    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  :20171130_105554.jpg.d80dc594900ab7848df2aaf88ac8509f.thumb.jpg.d3cd6896117b0246f0d2259c45b9e325.jpg

    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.

  6. Posté(e) · Modifié par nathan frad

    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.

  7. Posté(e)

    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°

  9. Posté(e)

    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

     

     

  10. Posté(e)

    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 ! 

  11. Posté(e)

    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 : 20171130_105554.thumb.jpg.e27787d12f63495523930527df0e2e46.jpg

    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...

  12. Posté(e) · Modifié par nathan frad
    Apres de longue recherche

    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 )

  13. Posté(e)

    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

  15. Posté(e) · Modifié par nathan frad

    :D 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 ! xD

     

    PS : J'essayerai plus tard de calculer la position des étoiles mais c'est pas gagné ^_^

  16. Posté(e) · Modifié par nathan frad

    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

     

     

     

  17. Posté(e)

    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 ...

     

     

     

     

     

     

     

  18. Posté(e)

    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;
       

     

     

     

  19. Posté(e)

    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 )

     

     

  20. Posté(e) · Modifié par nathan frad

    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)

    3ffa24d1f3f58aa66a4874b1b2d12b271da49ab6

     

    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 :D

     

  21. Posté(e)

    Le 09/10/2017 à 19:32, Yves BACHELET a dit :

    Il y a le livre de Jean Meeus qui explique comment calculer les positions des planètes :

    https://www.maison-astronomie.com/calculs/2950-calculs-astronomiques-a-l-usage-des-amateurs.html

     

    A mon avis, on peut le trouver dans toutes les bibliothèques universitaires ou des grandes écoles. Ça peut être ardu comme calcul !

     

    Merci, sur vos conseils nous avons commandé ce livre.

     

    Le 09/10/2017 à 21:42, bricodob300 a dit :

    Salut Nathan et bienvenue.

    Si tu veux pas te cogner les qques milliers de terme des formules : http://neoprogrammics.com/vsop87/

    tu choisis les astres, le référentiel, les coordonnées et le langage  et voilou.

    JP

     

     

    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.

  22. Posté(e)

    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.  ;)

  24. Posté(e)

    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. :)