Script PHP pour calculer la phase de la lune :)

[astro-fr.net]

Salut, comme promis, voici le code PHP pour calculer et afficher la phase de la Lune sous forme de petite image.

Vous pouvez voir ce que ça donne ici (http://www.astro-fr.net/phaselune.php) :) Et télécharger le script là (http://www.astro-fr.net/lune.zip).

J'espère qu'il n'y a pas trop de bugs ... en tout cas ça marche chez moi :D
Bon code à toutes et à tous :)



<?php

define("SYNODIC", 29.53058867); //constante pour la période synodique

define("MSPARJOUR", 24 * 60 * 60 * 1000); //constante pour le nombre de millisecondes par jour



function CalcPhase() {

$DateRef = mktime(8,53,0,7,29,2003); //date de référence : le 29/07/2003 à 6:53:00 GMT il y a eu une pleine lune

$Aujourdui = mktime(); //date du jour

$msDiff = ($Aujourdui - $DateRef)*1000 + MSPARJOUR; //on calcule la différence en millisecondes

$phase = ($msDiff * 100)/(SYNODIC * MSPARJOUR); //on calcule le pourcentage de la phase

while($phase>100) // tant que c'est supérieur à 100, on retire 100

{

$phase -= 100;

}

return $phase;

}



function ImgPhase($phase) {

$NumImage = round(29 * $phase / 100, 0); // on convertit le pourcentage en age de la lune (sur 29 jours)

if($NumImage==0) $NumImage==1;

$ImagePhase = 'day'.$NumImage.'.GIF';

return $ImagePhase;

}



function PartEntier($phase) {

if ( $phase <= 50.0 )

{

$plein = $phase * 2;

}

else

{

$plein = ( 100 - $phase ) * 2;

}

return $plein;

}



function NommerPhase($phase)

//convertit le pourcentage de lunaison en mots

{

if($phase >= 0 && $phase <= 2.49)

{

$NomPhase = "Nouvelle Lune";

}

else if($phase >= 2.5 && $phase <= 22.49)

{

$NomPhase = "Premier Croissant";

}

else if($phase >= 22.5 && $phase <= 27.49)

{

$NomPhase = "Premier Quartier";

}

else if($phase >= 27.5 && $phase <= 47.49)

{

$NomPhase = "Lune gibbeuse";

}

else if($phase >= 47.5 && $phase <= 52.49)

{

$NomPhase = "Pleine Lune";

}

else if($phase >= 52.5 && $phase <= 73.49)

{

$NomPhase = "Lune gibbeuse";

}

else if($phase >= 73.5 && $phase <= 77.49)

{

$NomPhase = "Dernier quartier";

}

else if($phase >= 77.5 && $phase <= 97.49)

{

$NomPhase = "Dernier croissant";

}

else

{

$NomPhase = "Nouvelle Lune";

}



return $NomPhase;

}



function JoursAvantNL($phase) { //calcule le nombre de jours avant la nouvelle Lune

$val = round((1 - $phase / 100) * SYNODIC,2);

return $val;

}



function JoursAvantPL($pourcentagePL) { //calcule le nombre de jours avant la pleine Lune

$val = round((1 - $pourcentagePL / 100) * SYNODIC,2);

return $val;

}

?>

[michel@astrocosmos.net]

Pas mal :D

[jfc_themia]

Salut la room.

je passais sur votre forum car je collecte différentes info concernant les lunaisons etc...

Voilà le prob.

Je viens d'integrer les fonctions dans un calendrier mensuel.
donc au moment de l'itération des dates pour le mois selectionné, j'envoie la date pointée afin de recupérer la n° du jour de lunaison.

Ca marche nickel, mais hélas ce n'est pas juste, pour le test du mois de juin 2005.
la nouvelle lune tombe juste mais pas la pleine lune(-1jr).

J'ai pensé que c'etait dû à l'heure passée, rien n'y fait, et pour le moment je traite les dates avec 12:00:00 Hr et ca ne colle toujours pas.

Donc c'est pas glop ! sorry.

Le systeme suivant :

http://perso.wanadoo.fr/pgj/phaselune.htm

fonctionne correctement. A étudier....

Cordialement.
Jean François.

[soulearth]

Moi je tente de programmer un script en php mais j'ai une erreur de 2h aujourd'hui et qui va jusqu'a 17h pour janvier 2012. Si vous pouviez m'aider; Ci dessous le code source un maximum commenté.



<?php

$periode = (1/((1/27.322)-(1/365.25))) *24*60*60; //formule de copernic donnant la durée d'une periode lunaire ,environ 29.5j

//time de nouvelle lune 592500

//------------------------------------------------//
//Explication time Stamp
// Le Timestamp est un compteur servant de référence temporelle.
//Il a pour base le 01 Janvier 1970 à 00h00 UTC (soit GMT) et depuis, ce compteur s'incrémente de 1 toutes les secondes.
// La fonction time() renvoi donc le nombre de secondes qui se sont ecoulées depuis le 01/01/1970 à 00h00 GMT.

//Info: La première nouvelle lune suivant le déclenchement du timestamp était le 7 janvier 1970 à 20h35m UTC se qui
// donne 592500 seconde, 6 jour complets + 20h35min sois 6*24*3600+20*3600+35*60=252500
//Le timestamp de cette date est égale au nombre de secondes qui se sont déroulé depuis le 01/01/1970 sois 592500.



$tps_depuis_nouvelle_lune = (time() - 592500) % $periode ;

//Je calcul le temps qui s'est déroulé depuis la dernière nouvelle lune.
//Pour cela je prend le timestamp actuel ( fonction time() ) moins le temps qu'il s'est écoulé entre le déclanchement du timestamp et la première nouvelle lune car en 70 ils ont pas déclenché le timestamp au moment de la nouvelle lune et je fais modulo de la période, se qui me donne le reste en seconde. Avec ca j'ai le nombre de secondes qui s'est ecoulé depuis la dernière nouvelle lune.



echo "le temps passé depuis la dernière nouvelle lune est ". $tps_depuis_nouvelle_lune. " secondes";

$tps_restant_avant_prochaine = $periode - $tps_depuis_nouvelle_lune;
// le temps restant avant la prochaine nouvelle lune <=> temp restant = periode - temp ecoulé



echo "<br>";
echo "time stamp". time(); //J'affiche le timestamp actuel
echo "<br>";

echo " le temps qu'il reste avant la prochaine est ". $tps_restant_avant_prochaine. " secondes";
// J'affiche le temps restant avant la prochaine nouvelle lune.
echo "<br>";



echo "Date et heure de la prochaine nouvelle lune". date('d/m/Y H:i:s', time() + $tps_restant_avant_prochaine);
// J'affiche la date et l'heure de la prochaine nouvelle lune et la j'ai une erreur de 2 heures pour le moi actuelle et qui va en grandissant jusqu'a 17h d'erreur pour l'année 2012. C'est la que j'ai besoin de votre aide. Je peu avoir une erreur en prog comme en astro mais je vois pas ou.

//Site utiles:
//convertisseur timestamp: http://loranger.free.fr/convertisseur/index.php
// Calendrier de la lune : http://www.imcce.fr/page.php?nav=fr/ephemerides/astronomie/phases_lune/index.php


?>

[pejive]

bonjour,
je n'ai pas encore lu tout ton script, mais déjà je pense qu'il y a un problème au début: la durée du jour n'est pas de 24h mais 23h56mn04s soit 86164 s...
edit: en fait non c'est la période de rotation de la Terre; je ne sais pas s'il faut prendre 86164 ou 86400 ??
edit 2: il vaudrait mieux prendre 365.256 au lieu de 365.25 pour obtenir une période synodique plus précise

[soulearth]

Ta réponse sur la durée en seconde d'une journée est trés pertinante mais aprés verification j'ai vu que pour le timestamp une journée fait 86400 secondes donc le problème ne doit pas etre ca. Merci quand même d'essayer de m'aider.

[pejive]

Je pense que le problème vient de la précision sur la période en secondes.
Sur une durée de 37 ans un écart de 2 heures correspond à une erreur d'environ 15s sur la période soit 0,000 012 jour.
Il faudrait utiliser directement une valeur précise de la période synodique (29,53058868 jours) comme l'a fait astro-fr.net
Sinon, je ne vois pas d'erreur dans le programme

[Olivier_C8_SE]

ca ne peut pas être précis comme ça, car ça calcule une révolution sinodique comme si elle était toujours égale à "$periode = (1/((1/27.322)-(1/365.25))) *24*60*60;"... mais pas du tout... en fait, ce serait bien trop simple.
J'ai su faire en VB, je sais que certains s'en sont approché en javascript... en php ça me parait difficile (environ 3000 lignes de code)

[Olivier_C8_SE]

Tiens, d'ailleurs c'est amusant, j'en ai un qui doit arriver au même résultat approché... visible ici: (au milieu de la page) Scripts PHP (http://olravet.fr/pagephp.php)

[Cristian]

Lorsque je chercher un script en php pour calculer les phases de la Lune, je suis tombé sur ce forum, il est vrai que le scripte proposé plus haut n'est pas très précis, surtout pour les dates s'écartant de la date de référence.

Voici un script beaucoup plus précis.

C'est le même scripte qui est utilisé par IMCCE (http://www.imcce.fr/page.php?nav=fr/ephemerides/astronomie/phases_lune/index.php)

<?
/*
Adaptation en php du fameux et excellent scripte Astro-MoonPhase de Brett Hamilton écrit en Perl.
http://search.cpan.org/~brett/Astro-MoonPhase-0.60/

Ce Scripte vous permettra de connaître, à une date donnée, l'illumination de la Lune, son age,
sa distance en km par rapport à la Terre, son angle en degrés, sa distance par rapport au soleil,
et son angle par rapport au soleil.

*/
class Moon
{
function phase($Year, $Month, $Day, $Hour, $Minutes, $Seconds)
{
$DateSec = mktime($Hour, $Minutes, $Seconds, $Month, $Day, $Year, 0);

ini_set(precision, "20"); //Defini la precision des calcules

# Astronomical constants.
$Epoch = 2444238.5; # 1980 January 0.0

# Constants defining the Sun's apparent orbit.
$Elonge = 278.833540; # ecliptic longitude of the Sun at epoch 1980.0
$Elongp = 282.596403; # ecliptic longitude of the Sun at perigee
$Eccent = 0.016718; # eccentricity of Earth's orbit
$Sunsmax = 1.495985e8; # semi-major axis of Earth's orbit, km
$Sunangsiz = 0.533128; # sun's angular size, degrees, at semi-major axis distance

# Elements of the Moon's orbit, epoch 1980.0.
$Mmlong = 64.975464; # moon's mean longitude at the epoch
$Mmlongp = 349.383063; # mean longitude of the perigee at the epoch
$Mlnode = 151.950429; # mean longitude of the node at the epoch
$Minc = 5.145396; # inclination of the Moon's orbit
$Mecc = 0.054900; # eccentricity of the Moon's orbit
$Mangsiz = 0.5181; # moon's angular size at distance a from Earth
$Msmax = 384401.0; # semi-major axis of Moon's orbit in km
$Mparallax = 0.9507; # parallax at distance a from Earth
$Synmonth = 29.53058868; # synodic month (new Moon to new Moon)

$pdate = Moon::jtime($DateSec);

$pphase; # illuminated fraction
$mage; # age of moon in days
$dist; # distance in kilometres
$angdia; # angular diameter in degrees
$sudist; # distance to Sun
$suangdia; # sun's angular diameter


# Calculation of the Sun's position.

$Day = $pdate - $Epoch; # date within epoch
$N = Moon::fixangle((360 / 365.2422) * $Day); # mean anomaly of the Sun
$M = Moon::fixangle($N + $Elonge - $Elongp); # convert from perigee
# co-ordinates to epoch 1980.0
$Ec = Moon::kepler($M, $Eccent); # solve equation of Kepler
$Ec = sqrt((1 + $Eccent) / (1 - $Eccent)) * tan($Ec / 2);
$Ec = 2 * Moon::todeg(atan($Ec)); # true anomaly
$Lambdasun = Moon::fixangle($Ec + $Elongp); # Sun's geocentric ecliptic
# longitude
# Orbital distance factor.
$F = ((1 + $Eccent * cos(Moon::torad($Ec))) / (1 - $Eccent * $Eccent));
$SunDist = $Sunsmax / $F; # distance to Sun in km
$SunAng = $F * $Sunangsiz; # Sun's angular size in degrees


# Calculation of the Moon's position.

# Moon's mean longitude.
$ml = Moon::fixangle(13.1763966 * $Day + $Mmlong);

# Moon's mean anomaly.
$MM = Moon::fixangle($ml - 0.1114041 * $Day - $Mmlongp);

# Moon's ascending node mean longitude.
$MN = Moon::fixangle($Mlnode - 0.0529539 * $Day);

# Evection.
$Ev = 1.2739 * sin(Moon::torad(2 * ($ml - $Lambdasun) - $MM));

# Annual equation.
$Ae = 0.1858 * sin(Moon::torad($M));

# Correction term.
$A3 = 0.37 * sin(Moon::torad($M));

# Corrected anomaly.
$MmP = $MM + $Ev - $Ae - $A3;

# Correction for the equation of the centre.
$mEc = 6.2886 * sin(Moon::torad($MmP));

# Another correction term.
$A4 = 0.214 * sin(Moon::torad(2 * $MmP));

# Corrected longitude.
$lP = $ml + $Ev + $mEc - $Ae + $A4;

# Variation.
$V = 0.6583 * sin(Moon::torad(2 * ($lP - $Lambdasun)));

# True longitude.
$lPP = $lP + $V;

# Corrected longitude of the node.
$NP = $MN - 0.16 * sin(Moon::torad($M));

# Y inclination coordinate.
$y = sin(Moon::torad($lPP - $NP)) * cos(Moon::torad($Minc));

# X inclination coordinate.
$x = cos(Moon::torad($lPP - $NP));

# Ecliptic longitude.
$Lambdamoon = Moon::todeg(atan2($y, $x));
$Lambdamoon += $NP;

# Ecliptic latitude.
$BetaM = Moon::todeg(asin(sin(Moon::torad($lPP - $NP)) * sin(Moon::torad($Minc))));

# Calculation of the phase of the Moon.

# Age of the Moon in degrees.
$MoonAge = $lPP - $Lambdasun;

# Phase of the Moon.
$MoonPhase = (1 - cos(Moon::torad($MoonAge))) / 2;

# Calculate distance of moon from the centre of the Earth.

$MoonDist = ($Msmax * (1 - $Mecc * $Mecc)) /
(1 + $Mecc * cos(Moon::torad($MmP + $mEc)));

# Calculate Moon's angular diameter.

$MoonDFrac = $MoonDist / $Msmax;
$MoonAng = $Mangsiz / $MoonDFrac;

# Calculate Moon's parallax.

$MoonPar = $Mparallax / $MoonDFrac;

$pphase = $MoonPhase; # illuminated fraction
$mage = $Synmonth * (Moon::fixangle($MoonAge) / 360.0); # age of moon in days
$dist = $MoonDist; # distance in kilometres
$angdia = $MoonAng; # angular diameter in degrees
$sudist = $SunDist; # distance to Sun
$suangdia = $SunAng; # sun's angular diameter
$mpfrac = Moon::fixangle($MoonAge) / 360.0;
return array( $pphase, $mage, $dist, $angdia, $sudist, $suangdia, $mpfrac, $mpfrac );
}

function fixangle($x) { return ($x - 360.0 * (floor($x / 360.0))); } # fix angle
function torad($x) { return ($x * (M_PI / 180.0)); } # deg->rad
function todeg($x) { return ($x * (180.0 / M_PI)); } # rad->deg

function jtime($t)
{
$julian = ($t / 86400) + 2440587.5; # (seconds /(seconds per day)) + julian date of epoch 2440587.5 / 86400 = 28,24753472222 Days
return ($julian);
}

function kepler($m, $ecc)
{
$EPSILON = 1e-6;

$m = Moon::torad($m);
$e = $m;
while (abs($delta) > $EPSILON)
{
$delta = $e - $ecc * sin($e) - $m;
$e -= $delta / (1 - $ecc * cos($e));
}
return ($e);
}

}


//Exemple d'utilisation :

//Pour le 11 Avril 2009 à 00h00
list($MoonPhase, $MoonAge, $MoonDist, $MoonAng, $SunDist, $SunAng, $mpfrac) = Moon::phase(2009, 04, 11, 00, 00, 01);
echo "La Lune est éclairée à ".number_format($MoonPhase*100, 2, ',', '')."%"."<br>";
echo "Son age est de ".number_format($MoonAge, 0, ',', '')." jours"."<br>";
echo "Et elle se situe à une distance de ".number_format($MoonDist, 0, ',', '')." km par rapport à la Terre."."<br>";
?>

Vous pouvez voir le scripte en fonctionnement sur www.CropCirclesDataBase.org dans la rubrique calendrier.

Au plaisir
Cristian