/**
* Finds the current location of the sun based on the observer location
*
* @param Predict_QTH $qth The observer location
* @param int $daynum The daynum or null to use the current daynum
*
* @return Predict_ObsSet
*/
public static function FindSun(Predict_QTH $qth, $daynum = null)
{
if ($daynum === null) {
$daynum = Predict_Time::get_current_daynum();
}
$obs_geodetic = new Predict_Geodetic();
$obs_geodetic->lon = $qth->lon * Predict::de2ra;
$obs_geodetic->lat = $qth->lat * Predict::de2ra;
$obs_geodetic->alt = $qth->alt / 1000.0;
$obs_geodetic->theta = 0;
$solar_vector = new Predict_Vector();
$zero_vector = new Predict_Vector();
$solar_set = new Predict_ObsSet();
self::Calculate_Solar_Position($daynum, $solar_vector);
Predict_SGPObs::Calculate_Obs($daynum, $solar_vector, $zero_vector, $obs_geodetic, $solar_set);
$solar_set->az = Predict_Math::Degrees($solar_set->az);
$solar_set->el = Predict_Math::Degrees($solar_set->el);
return $solar_set;
}
public static function Calculate_LatLonAlt($_time, Predict_Vector $pos, Predict_Geodetic $geodetic)
{
/* Reference: The 1992 Astronomical Almanac, page K12. */
/* double r,e2,phi,c; */
$geodetic->theta = Predict_Math::AcTan($pos->y, $pos->x);
/*radians*/
$geodetic->lon = Predict_Math::FMod2p($geodetic->theta - Predict_Time::ThetaG_JD($_time));
/*radians*/
$r = sqrt($pos->x * $pos->x + $pos->y * $pos->y);
$e2 = Predict::__f * (2 - Predict::__f);
$geodetic->lat = Predict_Math::AcTan($pos->z, $r);
/*radians*/
do {
$phi = $geodetic->lat;
$sinPhi = sin($phi);
$c = 1 / sqrt(1 - $e2 * ($sinPhi * $sinPhi));
$geodetic->lat = Predict_Math::AcTan($pos->z + Predict::xkmper * $c * $e2 * $sinPhi, $r);
} while (abs($geodetic->lat - $phi) >= 1.0E-10);
$geodetic->alt = $r / cos($geodetic->lat) - Predict::xkmper * $c;
/*kilometers*/
if ($geodetic->lat > Predict::pio2) {
$geodetic->lat -= Predict::twopi;
}
}
require_once 'Predict/Time.php';
// Set to true to profile
$profile = true;
// Example check at the top of your application for enabling profiling
if ($profile && extension_loaded('xhprof')) {
$xproflib = './xhprof_lib';
include_once $xproflib . '/utils/xhprof_lib.php';
include_once $xproflib . '/utils/xhprof_runs.php';
xhprof_enable(XHPROF_FLAGS_CPU + XHPROF_FLAGS_MEMORY);
function stop_xhprof_profiling()
{
$xhprofData = xhprof_disable();
$xhprofNameSpace = 'predict';
$xhprofRuns = new XHProfRuns_Default();
$xhprofRunID = $xhprofRuns->save_run($xhprofData, $xhprofNameSpace);
}
register_shutdown_function('stop_xhprof_profiling');
}
$start = microtime(true);
$predict = new Predict();
$qth = new Predict_QTH();
$qth->lat = 37.655;
$qth->lon = -122.407;
$qth->alt = 0;
$tleFile = file('examples/iss.tle');
$tle = new Predict_TLE($tleFile[0], $tleFile[1], $tleFile[2]);
$sat = new Predict_Sat($tle);
$now = Predict_Time::get_current_daynum();
$results = $predict->get_passes($sat, $qth, $now, 10);
echo "Execution time: " . number_format((microtime(true) - $start) * 1000, 2) . "ms\n";
exit;
public function Deep($ientry, Predict_Sat $sat)
{
switch ($ientry) {
case self::dpinit:
/* Entrance for deep space initialization */
$sat->dps->thgr = Predict_Time::ThetaG($sat->tle->epoch, $sat->deep_arg);
$eq = $sat->tle->eo;
$sat->dps->xnq = $sat->deep_arg->xnodp;
$aqnv = 1.0 / $sat->deep_arg->aodp;
$sat->dps->xqncl = $sat->tle->xincl;
$xmao = $sat->tle->xmo;
$xpidot = $sat->deep_arg->omgdot + $sat->deep_arg->xnodot;
$sinq = sin($sat->tle->xnodeo);
$cosq = cos($sat->tle->xnodeo);
$sat->dps->omegaq = $sat->tle->omegao;
$sat->dps->preep = 0;
/* Initialize lunar solar terms */
$day = $sat->deep_arg->ds50 + 18261.5;
/* Days since 1900 Jan 0.5 */
if ($day != $sat->dps->preep) {
$sat->dps->preep = $day;
$xnodce = 4.523602 - 0.00092422029 * $day;
$stem = sin($xnodce);
$ctem = cos($xnodce);
$sat->dps->zcosil = 0.91375164 - 0.03568096 * $ctem;
$sat->dps->zsinil = sqrt(1.0 - $sat->dps->zcosil * $sat->dps->zcosil);
$sat->dps->zsinhl = 0.08968351099999999 * $stem / $sat->dps->zsinil;
$sat->dps->zcoshl = sqrt(1.0 - $sat->dps->zsinhl * $sat->dps->zsinhl);
$c = 4.7199672 + 0.2299715 * $day;
$gam = 5.8351514 + 0.001944368 * $day;
$sat->dps->zmol = Predict_Math::FMod2p($c - $gam);
$zx = 0.39785416 * $stem / $sat->dps->zsinil;
$zy = $sat->dps->zcoshl * $ctem + 0.91744867 * $sat->dps->zsinhl * $stem;
$zx = Predict_Math::AcTan($zx, $zy);
$zx = $gam + $zx - $xnodce;
$sat->dps->zcosgl = cos($zx);
$sat->dps->zsingl = sin($zx);
$sat->dps->zmos = 6.2565837 + 0.017201977 * $day;
$sat->dps->zmos = Predict_Math::FMod2p($sat->dps->zmos);
}
/* End if(day != preep) */
/* Do solar terms */
$sat->dps->savtsn = 1.0E+20;
$zcosg = Predict::zcosgs;
$zsing = Predict::zsings;
$zcosi = Predict::zcosis;
$zsini = Predict::zsinis;
$zcosh = $cosq;
$zsinh = $sinq;
$cc = Predict::c1ss;
$zn = Predict::zns;
$ze = Predict::zes;
$zmo = $sat->dps->zmos;
$xnoi = 1.0 / $sat->dps->xnq;
/* Loop breaks when Solar terms are done a second */
/* time, after Lunar terms are initialized */
for (;;) {
/* Solar terms done again after Lunar terms are done */
$a1 = $zcosg * $zcosh + $zsing * $zcosi * $zsinh;
$a3 = -$zsing * $zcosh + $zcosg * $zcosi * $zsinh;
$a7 = -$zcosg * $zsinh + $zsing * $zcosi * $zcosh;
$a8 = $zsing * $zsini;
$a9 = $zsing * $zsinh + $zcosg * $zcosi * $zcosh;
$a10 = $zcosg * $zsini;
$a2 = $sat->deep_arg->cosio * $a7 + $sat->deep_arg->sinio * $a8;
$a4 = $sat->deep_arg->cosio * $a9 + $sat->deep_arg->sinio * $a10;
$a5 = -$sat->deep_arg->sinio * $a7 + $sat->deep_arg->cosio * $a8;
$a6 = -$sat->deep_arg->sinio * $a9 + $sat->deep_arg->cosio * $a10;
$x1 = $a1 * $sat->deep_arg->cosg + $a2 * $sat->deep_arg->sing;
$x2 = $a3 * $sat->deep_arg->cosg + $a4 * $sat->deep_arg->sing;
$x3 = -$a1 * $sat->deep_arg->sing + $a2 * $sat->deep_arg->cosg;
$x4 = -$a3 * $sat->deep_arg->sing + $a4 * $sat->deep_arg->cosg;
$x5 = $a5 * $sat->deep_arg->sing;
$x6 = $a6 * $sat->deep_arg->sing;
$x7 = $a5 * $sat->deep_arg->cosg;
$x8 = $a6 * $sat->deep_arg->cosg;
$z31 = 12 * $x1 * $x1 - 3 * $x3 * $x3;
$z32 = 24 * $x1 * $x2 - 6 * $x3 * $x4;
$z33 = 12 * $x2 * $x2 - 3 * $x4 * $x4;
$z1 = 3 * ($a1 * $a1 + $a2 * $a2) + $z31 * $sat->deep_arg->eosq;
$z2 = 6 * ($a1 * $a3 + $a2 * $a4) + $z32 * $sat->deep_arg->eosq;
$z3 = 3 * ($a3 * $a3 + $a4 * $a4) + $z33 * $sat->deep_arg->eosq;
$z11 = -6 * $a1 * $a5 + $sat->deep_arg->eosq * (-24 * $x1 * $x7 - 6 * $x3 * $x5);
$z12 = -6 * ($a1 * $a6 + $a3 * $a5) + $sat->deep_arg->eosq * (-24 * ($x2 * $x7 + $x1 * $x8) - 6 * ($x3 * $x6 + $x4 * $x5));
$z13 = -6 * $a3 * $a6 + $sat->deep_arg->eosq * (-24 * $x2 * $x8 - 6 * $x4 * $x6);
$z21 = 6 * $a2 * $a5 + $sat->deep_arg->eosq * (24 * $x1 * $x5 - 6 * $x3 * $x7);
$z22 = 6 * ($a4 * $a5 + $a2 * $a6) + $sat->deep_arg->eosq * (24 * ($x2 * $x5 + $x1 * $x6) - 6 * ($x4 * $x7 + $x3 * $x8));
$z23 = 6 * $a4 * $a6 + $sat->deep_arg->eosq * (24 * $x2 * $x6 - 6 * $x4 * $x8);
$z1 = $z1 + $z1 + $sat->deep_arg->betao2 * $z31;
$z2 = $z2 + $z2 + $sat->deep_arg->betao2 * $z32;
$z3 = $z3 + $z3 + $sat->deep_arg->betao2 * $z33;
$s3 = $cc * $xnoi;
$s2 = -0.5 * $s3 / $sat->deep_arg->betao;
$s4 = $s3 * $sat->deep_arg->betao;
$s1 = -15 * $eq * $s4;
$s5 = $x1 * $x3 + $x2 * $x4;
$s6 = $x2 * $x3 + $x1 * $x4;
$s7 = $x2 * $x4 - $x1 * $x3;
$se = $s1 * $zn * $s5;
$si = $s2 * $zn * ($z11 + $z13);
$sl = -$zn * $s3 * ($z1 + $z3 - 14 - 6 * $sat->deep_arg->eosq);
$sgh = $s4 * $zn * ($z31 + $z33 - 6);
$sh = -$zn * $s2 * ($z21 + $z23);
if ($sat->dps->xqncl < 0.052359877) {
$sh = 0;
}
$sat->dps->ee2 = 2 * $s1 * $s6;
$sat->dps->e3 = 2 * $s1 * $s7;
$sat->dps->xi2 = 2 * $s2 * $z12;
$sat->dps->xi3 = 2 * $s2 * ($z13 - $z11);
$sat->dps->xl2 = -2 * $s3 * $z2;
$sat->dps->xl3 = -2 * $s3 * ($z3 - $z1);
$sat->dps->xl4 = -2 * $s3 * (-21 - 9 * $sat->deep_arg->eosq) * $ze;
$sat->dps->xgh2 = 2 * $s4 * $z32;
$sat->dps->xgh3 = 2 * $s4 * ($z33 - $z31);
$sat->dps->xgh4 = -18 * $s4 * $ze;
$sat->dps->xh2 = -2 * $s2 * $z22;
$sat->dps->xh3 = -2 * $s2 * ($z23 - $z21);
if ($sat->flags & self::LUNAR_TERMS_DONE_FLAG) {
break;
}
/* Do lunar terms */
$sat->dps->sse = $se;
$sat->dps->ssi = $si;
$sat->dps->ssl = $sl;
$sat->dps->ssh = $sh / $sat->deep_arg->sinio;
$sat->dps->ssg = $sgh - $sat->deep_arg->cosio * $sat->dps->ssh;
$sat->dps->se2 = $sat->dps->ee2;
$sat->dps->si2 = $sat->dps->xi2;
$sat->dps->sl2 = $sat->dps->xl2;
$sat->dps->sgh2 = $sat->dps->xgh2;
$sat->dps->sh2 = $sat->dps->xh2;
$sat->dps->se3 = $sat->dps->e3;
$sat->dps->si3 = $sat->dps->xi3;
$sat->dps->sl3 = $sat->dps->xl3;
$sat->dps->sgh3 = $sat->dps->xgh3;
$sat->dps->sh3 = $sat->dps->xh3;
$sat->dps->sl4 = $sat->dps->xl4;
$sat->dps->sgh4 = $sat->dps->xgh4;
$zcosg = $sat->dps->zcosgl;
$zsing = $sat->dps->zsingl;
$zcosi = $sat->dps->zcosil;
$zsini = $sat->dps->zsinil;
$zcosh = $sat->dps->zcoshl * $cosq + $sat->dps->zsinhl * $sinq;
$zsinh = $sinq * $sat->dps->zcoshl - $cosq * $sat->dps->zsinhl;
$zn = Predict::znl;
$cc = Predict::c1l;
$ze = Predict::zel;
$zmo = $sat->dps->zmol;
$sat->flags |= self::LUNAR_TERMS_DONE_FLAG;
}
/* End of for(;;) */
$sat->dps->sse = $sat->dps->sse + $se;
$sat->dps->ssi = $sat->dps->ssi + $si;
$sat->dps->ssl = $sat->dps->ssl + $sl;
$sat->dps->ssg = $sat->dps->ssg + $sgh - $sat->deep_arg->cosio / $sat->deep_arg->sinio * $sh;
$sat->dps->ssh = $sat->dps->ssh + $sh / $sat->deep_arg->sinio;
/* Geopotential resonance initialization for 12 hour orbits */
$sat->flags &= ~self::RESONANCE_FLAG;
$sat->flags &= ~self::SYNCHRONOUS_FLAG;
if (!($sat->dps->xnq < 0.0052359877 && $sat->dps->xnq > 0.0034906585)) {
if ($sat->dps->xnq < 0.00826 || $sat->dps->xnq > 0.00924) {
return;
}
if ($eq < 0.5) {
return;
}
$sat->flags |= self::RESONANCE_FLAG;
$eoc = $eq * $sat->deep_arg->eosq;
$g201 = -0.306 - ($eq - 0.64) * 0.44;
if ($eq <= 0.65) {
$g211 = 3.616 - 13.247 * $eq + 16.29 * $sat->deep_arg->eosq;
$g310 = -19.302 + 117.39 * $eq - 228.419 * $sat->deep_arg->eosq + 156.591 * $eoc;
$g322 = -18.9068 + 109.7927 * $eq - 214.6334 * $sat->deep_arg->eosq + 146.5816 * $eoc;
$g410 = -41.122 + 242.694 * $eq - 471.094 * $sat->deep_arg->eosq + 313.953 * $eoc;
$g422 = -146.407 + 841.88 * $eq - 1629.014 * $sat->deep_arg->eosq + 1083.435 * $eoc;
$g520 = -532.114 + 3017.977 * $eq - 5740 * $sat->deep_arg->eosq + 3708.276 * $eoc;
} else {
$g211 = -72.099 + 331.819 * $eq - 508.738 * $sat->deep_arg->eosq + 266.724 * $eoc;
$g310 = -346.844 + 1582.851 * $eq - 2415.925 * $sat->deep_arg->eosq + 1246.113 * $eoc;
$g322 = -342.585 + 1554.908 * $eq - 2366.899 * $sat->deep_arg->eosq + 1215.972 * $eoc;
$g410 = -1052.797 + 4758.686 * $eq - 7193.992 * $sat->deep_arg->eosq + 3651.957 * $eoc;
$g422 = -3581.69 + 16178.11 * $eq - 24462.77 * $sat->deep_arg->eosq + 12422.52 * $eoc;
if ($eq <= 0.715) {
$g520 = 1464.74 - 4664.75 * $eq + 3763.64 * $sat->deep_arg->eosq;
} else {
$g520 = -5149.66 + 29936.92 * $eq - 54087.36 * $sat->deep_arg->eosq + 31324.56 * $eoc;
}
}
/* End if (eq <= 0.65) */
if ($eq < 0.7) {
$g533 = -919.2277 + 4988.61 * $eq - 9064.77 * $sat->deep_arg->eosq + 5542.21 * $eoc;
$g521 = -822.71072 + 4568.6173 * $eq - 8491.4146 * $sat->deep_arg->eosq + 5337.524 * $eoc;
$g532 = -853.6660000000001 + 4690.25 * $eq - 8624.77 * $sat->deep_arg->eosq + 5341.4 * $eoc;
} else {
$g533 = -37995.78 + 161616.52 * $eq - 229838.2 * $sat->deep_arg->eosq + 109377.94 * $eoc;
$g521 = -51752.104 + 218913.95 * $eq - 309468.16 * $sat->deep_arg->eosq + 146349.42 * $eoc;
$g532 = -40023.88 + 170470.89 * $eq - 242699.48 * $sat->deep_arg->eosq + 115605.82 * $eoc;
}
/* End if (eq <= 0.7) */
$sini2 = $sat->deep_arg->sinio * $sat->deep_arg->sinio;
$f220 = 0.75 * (1 + 2 * $sat->deep_arg->cosio + $sat->deep_arg->theta2);
$f221 = 1.5 * $sini2;
$f321 = 1.875 * $sat->deep_arg->sinio * (1 - 2 * $sat->deep_arg->cosio - 3 * $sat->deep_arg->theta2);
$f322 = -1.875 * $sat->deep_arg->sinio * (1 + 2 * $sat->deep_arg->cosio - 3 * $sat->deep_arg->theta2);
$f441 = 35 * $sini2 * $f220;
$f442 = 39.375 * $sini2 * $sini2;
$f522 = 9.84375 * $sat->deep_arg->sinio * ($sini2 * (1 - 2 * $sat->deep_arg->cosio - 5 * $sat->deep_arg->theta2) + 0.33333333 * (-2 + 4 * $sat->deep_arg->cosio + 6 * $sat->deep_arg->theta2));
$f523 = $sat->deep_arg->sinio * (4.92187512 * $sini2 * (-2 - 4 * $sat->deep_arg->cosio + 10 * $sat->deep_arg->theta2) + 6.56250012 * (1 + 2 * $sat->deep_arg->cosio - 3 * $sat->deep_arg->theta2));
$f542 = 29.53125 * $sat->deep_arg->sinio * (2 - 8 * $sat->deep_arg->cosio + $sat->deep_arg->theta2 * (-12 + 8 * $sat->deep_arg->cosio + 10 * $sat->deep_arg->theta2));
$f543 = 29.53125 * $sat->deep_arg->sinio * (-2 - 8 * $sat->deep_arg->cosio + $sat->deep_arg->theta2 * (12 + 8 * $sat->deep_arg->cosio - 10 * $sat->deep_arg->theta2));
$xno2 = $sat->dps->xnq * $sat->dps->xnq;
$ainv2 = $aqnv * $aqnv;
$temp1 = 3 * $xno2 * $ainv2;
$temp = $temp1 * Predict::root22;
$sat->dps->d2201 = $temp * $f220 * $g201;
$sat->dps->d2211 = $temp * $f221 * $g211;
$temp1 = $temp1 * $aqnv;
$temp = $temp1 * Predict::root32;
$sat->dps->d3210 = $temp * $f321 * $g310;
$sat->dps->d3222 = $temp * $f322 * $g322;
$temp1 = $temp1 * $aqnv;
$temp = 2 * $temp1 * Predict::root44;
$sat->dps->d4410 = $temp * $f441 * $g410;
$sat->dps->d4422 = $temp * $f442 * $g422;
$temp1 = $temp1 * $aqnv;
$temp = $temp1 * Predict::root52;
$sat->dps->d5220 = $temp * $f522 * $g520;
$sat->dps->d5232 = $temp * $f523 * $g532;
$temp = 2 * $temp1 * Predict::root54;
$sat->dps->d5421 = $temp * $f542 * $g521;
$sat->dps->d5433 = $temp * $f543 * $g533;
$sat->dps->xlamo = $xmao + $sat->tle->xnodeo + $sat->tle->xnodeo - $sat->dps->thgr - $sat->dps->thgr;
$bfact = $sat->deep_arg->xmdot + $sat->deep_arg->xnodot + $sat->deep_arg->xnodot - Predict::thdt - Predict::thdt;
$bfact = $bfact + $sat->dps->ssl + $sat->dps->ssh + $sat->dps->ssh;
} else {
$sat->flags |= self::RESONANCE_FLAG;
$sat->flags |= self::SYNCHRONOUS_FLAG;
/* Synchronous resonance terms initialization */
$g200 = 1 + $sat->deep_arg->eosq * (-2.5 + 0.8125 * $sat->deep_arg->eosq);
$g310 = 1 + 2 * $sat->deep_arg->eosq;
$g300 = 1 + $sat->deep_arg->eosq * (-6 + 6.60937 * $sat->deep_arg->eosq);
$f220 = 0.75 * (1 + $sat->deep_arg->cosio) * (1 + $sat->deep_arg->cosio);
$f311 = 0.9375 * $sat->deep_arg->sinio * $sat->deep_arg->sinio * (1 + 3 * $sat->deep_arg->cosio) - 0.75 * (1 + $sat->deep_arg->cosio);
$f330 = 1 + $sat->deep_arg->cosio;
$f330 = 1.875 * $f330 * $f330 * $f330;
$sat->dps->del1 = 3 * $sat->dps->xnq * $sat->dps->xnq * $aqnv * $aqnv;
$sat->dps->del2 = 2 * $sat->dps->del1 * $f220 * $g200 * Predict::q22;
$sat->dps->del3 = 3 * $sat->dps->del1 * $f330 * $g300 * Predict::q33 * $aqnv;
$sat->dps->del1 = $sat->dps->del1 * $f311 * $g310 * Predict::q31 * $aqnv;
$sat->dps->fasx2 = 0.13130908;
$sat->dps->fasx4 = 2.8843198;
$sat->dps->fasx6 = 0.37448087;
$sat->dps->xlamo = $xmao + $sat->tle->xnodeo + $sat->tle->omegao - $sat->dps->thgr;
$bfact = $sat->deep_arg->xmdot + $xpidot - Predict::thdt;
$bfact = $bfact + $sat->dps->ssl + $sat->dps->ssg + $sat->dps->ssh;
}
/* End if( !(xnq < 0.0052359877) && (xnq > 0.0034906585) ) */
$sat->dps->xfact = $bfact - $sat->dps->xnq;
/* Initialize integrator */
$sat->dps->xli = $sat->dps->xlamo;
$sat->dps->xni = $sat->dps->xnq;
$sat->dps->atime = 0;
$sat->dps->stepp = 720;
$sat->dps->stepn = -720;
$sat->dps->step2 = 259200;
/* End case self::dpinit: */
return;
case self::dpsec:
/* Entrance for deep space secular effects */
$sat->deep_arg->xll = $sat->deep_arg->xll + $sat->dps->ssl * $sat->deep_arg->t;
$sat->deep_arg->omgadf = $sat->deep_arg->omgadf + $sat->dps->ssg * $sat->deep_arg->t;
$sat->deep_arg->xnode = $sat->deep_arg->xnode + $sat->dps->ssh * $sat->deep_arg->t;
$sat->deep_arg->em = $sat->tle->eo + $sat->dps->sse * $sat->deep_arg->t;
$sat->deep_arg->xinc = $sat->tle->xincl + $sat->dps->ssi * $sat->deep_arg->t;
if ($sat->deep_arg->xinc < 0) {
$sat->deep_arg->xinc = -$sat->deep_arg->xinc;
$sat->deep_arg->xnode = $sat->deep_arg->xnode + Predict::pi;
$sat->deep_arg->omgadf = $sat->deep_arg->omgadf - Predict::pi;
}
if (~$sat->flags & self::RESONANCE_FLAG) {
return;
}
do {
if ($sat->dps->atime == 0 || $sat->deep_arg->t >= 0 && $sat->dps->atime < 0 || $sat->deep_arg->t < 0 && $sat->dps->atime >= 0) {
/* Epoch restart */
if ($sat->deep_arg->t >= 0) {
$delt = $sat->dps->stepp;
} else {
$delt = $sat->dps->stepn;
}
$sat->dps->atime = 0;
$sat->dps->xni = $sat->dps->xnq;
$sat->dps->xli = $sat->dps->xlamo;
} else {
if (abs($sat->deep_arg->t) >= abs($sat->dps->atime)) {
if ($sat->deep_arg->t > 0) {
$delt = $sat->dps->stepp;
} else {
$delt = $sat->dps->stepn;
}
}
}
do {
if (abs($sat->deep_arg->t - $sat->dps->atime) >= $sat->dps->stepp) {
$sat->flags |= self::DO_LOOP_FLAG;
$sat->flags &= ~self::EPOCH_RESTART_FLAG;
} else {
$ft = $sat->deep_arg->t - $sat->dps->atime;
$sat->flags &= ~self::DO_LOOP_FLAG;
}
if (abs($sat->deep_arg->t) < abs($sat->dps->atime)) {
if ($sat->deep_arg->t >= 0) {
$delt = $sat->dps->stepn;
} else {
$delt = $sat->dps->stepp;
}
$sat->flags |= self::DO_LOOP_FLAG | self::EPOCH_RESTART_FLAG;
}
/* Dot terms calculated */
if ($sat->flags & self::SYNCHRONOUS_FLAG) {
$xndot = $sat->dps->del1 * sin($sat->dps->xli - $sat->dps->fasx2) + $sat->dps->del2 * sin(2 * ($sat->dps->xli - $sat->dps->fasx4)) + $sat->dps->del3 * sin(3 * ($sat->dps->xli - $sat->dps->fasx6));
$xnddt = $sat->dps->del1 * cos($sat->dps->xli - $sat->dps->fasx2) + 2 * $sat->dps->del2 * cos(2 * ($sat->dps->xli - $sat->dps->fasx4)) + 3 * $sat->dps->del3 * cos(3 * ($sat->dps->xli - $sat->dps->fasx6));
} else {
$xomi = $sat->dps->omegaq + $sat->deep_arg->omgdot * $sat->dps->atime;
$x2omi = $xomi + $xomi;
$x2li = $sat->dps->xli + $sat->dps->xli;
$xndot = $sat->dps->d2201 * sin($x2omi + $sat->dps->xli - Predict::g22) + $sat->dps->d2211 * sin($sat->dps->xli - Predict::g22) + $sat->dps->d3210 * sin($xomi + $sat->dps->xli - Predict::g32) + $sat->dps->d3222 * sin(-$xomi + $sat->dps->xli - Predict::g32) + $sat->dps->d4410 * sin($x2omi + $x2li - Predict::g44) + $sat->dps->d4422 * sin($x2li - Predict::g44) + $sat->dps->d5220 * sin($xomi + $sat->dps->xli - Predict::g52) + $sat->dps->d5232 * sin(-$xomi + $sat->dps->xli - Predict::g52) + $sat->dps->d5421 * sin($xomi + $x2li - Predict::g54) + $sat->dps->d5433 * sin(-$xomi + $x2li - Predict::g54);
$xnddt = $sat->dps->d2201 * cos($x2omi + $sat->dps->xli - Predict::g22) + $sat->dps->d2211 * cos($sat->dps->xli - Predict::g22) + $sat->dps->d3210 * cos($xomi + $sat->dps->xli - Predict::g32) + $sat->dps->d3222 * cos(-$xomi + $sat->dps->xli - Predict::g32) + $sat->dps->d5220 * cos($xomi + $sat->dps->xli - Predict::g52) + $sat->dps->d5232 * cos(-$xomi + $sat->dps->xli - Predict::g52) + 2 * ($sat->dps->d4410 * cos($x2omi + $x2li - Predict::g44) + $sat->dps->d4422 * cos($x2li - Predict::g44) + $sat->dps->d5421 * cos($xomi + $x2li - Predict::g54) + $sat->dps->d5433 * cos(-$xomi + $x2li - Predict::g54));
}
/* End of if (isFlagSet(SYNCHRONOUS_FLAG)) */
$xldot = $sat->dps->xni + $sat->dps->xfact;
$xnddt = $xnddt * $xldot;
if ($sat->flags & self::DO_LOOP_FLAG) {
$sat->dps->xli = $sat->dps->xli + $xldot * $delt + $xndot * $sat->dps->step2;
$sat->dps->xni = $sat->dps->xni + $xndot * $delt + $xnddt * $sat->dps->step2;
$sat->dps->atime = $sat->dps->atime + $delt;
}
} while ($sat->flags & self::DO_LOOP_FLAG && ~$sat->flags & self::EPOCH_RESTART_FLAG);
} while ($sat->flags & self::DO_LOOP_FLAG && $sat->flags & self::EPOCH_RESTART_FLAG);
$sat->deep_arg->xn = $sat->dps->xni + $xndot * $ft + $xnddt * $ft * $ft * 0.5;
$xl = $sat->dps->xli + $xldot * $ft + $xndot * $ft * $ft * 0.5;
$temp = -$sat->deep_arg->xnode + $sat->dps->thgr + $sat->deep_arg->t * Predict::thdt;
if (~$sat->flags & self::SYNCHRONOUS_FLAG) {
$sat->deep_arg->xll = $xl + $temp + $temp;
} else {
$sat->deep_arg->xll = $xl - $sat->deep_arg->omgadf + $temp;
}
return;
/* End case dpsec: */
/* End case dpsec: */
case self::dpper:
/* Entrance for lunar-solar periodics */
$sinis = sin($sat->deep_arg->xinc);
$cosis = cos($sat->deep_arg->xinc);
if (abs($sat->dps->savtsn - $sat->deep_arg->t) >= 30) {
$sat->dps->savtsn = $sat->deep_arg->t;
$zm = $sat->dps->zmos + Predict::zns * $sat->deep_arg->t;
$zf = $zm + 2 * Predict::zes * sin($zm);
$sinzf = sin($zf);
$f2 = 0.5 * $sinzf * $sinzf - 0.25;
$f3 = -0.5 * $sinzf * cos($zf);
$ses = $sat->dps->se2 * $f2 + $sat->dps->se3 * $f3;
$sis = $sat->dps->si2 * $f2 + $sat->dps->si3 * $f3;
$sls = $sat->dps->sl2 * $f2 + $sat->dps->sl3 * $f3 + $sat->dps->sl4 * $sinzf;
$sat->dps->sghs = $sat->dps->sgh2 * $f2 + $sat->dps->sgh3 * $f3 + $sat->dps->sgh4 * $sinzf;
$sat->dps->shs = $sat->dps->sh2 * $f2 + $sat->dps->sh3 * $f3;
$zm = $sat->dps->zmol + Predict::znl * $sat->deep_arg->t;
$zf = $zm + 2 * Predict::zel * sin($zm);
$sinzf = sin($zf);
$f2 = 0.5 * $sinzf * $sinzf - 0.25;
$f3 = -0.5 * $sinzf * cos($zf);
$sel = $sat->dps->ee2 * $f2 + $sat->dps->e3 * $f3;
$sil = $sat->dps->xi2 * $f2 + $sat->dps->xi3 * $f3;
$sll = $sat->dps->xl2 * $f2 + $sat->dps->xl3 * $f3 + $sat->dps->xl4 * $sinzf;
$sat->dps->sghl = $sat->dps->xgh2 * $f2 + $sat->dps->xgh3 * $f3 + $sat->dps->xgh4 * $sinzf;
$sat->dps->sh1 = $sat->dps->xh2 * $f2 + $sat->dps->xh3 * $f3;
$sat->dps->pe = $ses + $sel;
$sat->dps->pinc = $sis + $sil;
$sat->dps->pl = $sls + $sll;
}
$pgh = $sat->dps->sghs + $sat->dps->sghl;
$ph = $sat->dps->shs + $sat->dps->sh1;
$sat->deep_arg->xinc = $sat->deep_arg->xinc + $sat->dps->pinc;
$sat->deep_arg->em = $sat->deep_arg->em + $sat->dps->pe;
if ($sat->dps->xqncl >= 0.2) {
/* Apply periodics directly */
$ph = $ph / $sat->deep_arg->sinio;
$pgh = $pgh - $sat->deep_arg->cosio * $ph;
$sat->deep_arg->omgadf = $sat->deep_arg->omgadf + $pgh;
$sat->deep_arg->xnode = $sat->deep_arg->xnode + $ph;
$sat->deep_arg->xll = $sat->deep_arg->xll + $sat->dps->pl;
} else {
/* Apply periodics with Lyddane modification */
$sinok = sin($sat->deep_arg->xnode);
$cosok = cos($sat->deep_arg->xnode);
$alfdp = $sinis * $sinok;
$betdp = $sinis * $cosok;
$dalf = $ph * $cosok + $sat->dps->pinc * $cosis * $sinok;
$dbet = -$ph * $sinok + $sat->dps->pinc * $cosis * $cosok;
$alfdp = $alfdp + $dalf;
$betdp = $betdp + $dbet;
$sat->deep_arg->xnode = Predict_Math::FMod2p($sat->deep_arg->xnode);
$xls = $sat->deep_arg->xll + $sat->deep_arg->omgadf + $cosis * $sat->deep_arg->xnode;
$dls = $sat->dps->pl + $pgh - $sat->dps->pinc * $sat->deep_arg->xnode * $sinis;
$xls = $xls + $dls;
$xnoh = $sat->deep_arg->xnode;
$sat->deep_arg->xnode = Predict_Math::AcTan($alfdp, $betdp);
/* This is a patch to Lyddane modification */
/* suggested by Rob Matson. */
if (abs($xnoh - $sat->deep_arg->xnode) > Predict::pi) {
if ($sat->deep_arg->xnode < $xnoh) {
$sat->deep_arg->xnode += Predict::twopi;
} else {
$sat->deep_arg->xnode -= Predict::twopi;
}
}
$sat->deep_arg->xll = $sat->deep_arg->xll + $sat->dps->pl;
$sat->deep_arg->omgadf = $xls - $sat->deep_arg->xll - cos($sat->deep_arg->xinc) * $sat->deep_arg->xnode;
}
/* End case dpper: */
return;
}
/* End switch(ientry) */
}
<?php
/**
* An example for looking up the current solar position in the sky from a given
* location (elevation/azimuth)
*/
require_once 'Predict/Solar.php';
require_once 'Predict/QTH.php';
// Use current time in the form of a daynum
$time = time();
$daynum = Predict_Time::unix2daynum($time);
// Set up the observer position on the ground
$qth = new Predict_QTH();
$qth->lat = 37.786759;
$qth->lon = -122.405162;
$qth->alt = 10;
// Altitude above sea level in meters
$sunInfo = Predict_Solar::FindSun($qth, $daynum);
$output = array('elevation' => $sunInfo->el, 'azimuth' => $sunInfo->az, 'timestamp' => $time);
// output results
echo json_encode($output);
// You can modify some preferences in Predict(), the defaults are below
//
// $predict->minEle = 10; // Minimum elevation for a pass
// $predict->timeRes = 10; // Pass details: time resolution in seconds
// $predict->numEntries = 20; // Pass details: number of entries per pass
// $predict->threshold = -6; // Twilight threshold (sun must be at this lat or lower)
// Get the passes and filter visible only, takes about 4 seconds for 10 days
$results = $predict->get_passes($sat, $qth, $now, 10);
$filtered = $predict->filterVisiblePasses($results);
$zone = 'America/Los_Angeles';
// Pacific time zone
$format = 'm-d-Y H:i:s';
// Time format from PHP's date() function
// Format the output similar to the heavens-above.com website
foreach ($filtered as $pass) {
echo "AOS Daynum: " . $pass->visible_aos . "\n";
echo "AOS Time: " . Predict_Time::daynum2readable($pass->visible_aos, $zone, $format) . "\n";
echo "AOS Az: " . $predict->azDegreesToDirection($pass->visible_aos_az) . "\n";
echo "AOS El: " . round($pass->visible_aos_el) . "\n";
echo "Max Time: " . Predict_Time::daynum2readable($pass->visible_tca, $zone, $format) . "\n";
echo "Max Az: " . $predict->azDegreesToDirection($pass->visible_max_el_az) . "\n";
echo "Max El: " . round($pass->visible_max_el) . "\n";
echo "LOS Time: " . Predict_Time::daynum2readable($pass->visible_los, $zone, $format) . "\n";
echo "LOS Az: " . $predict->azDegreesToDirection($pass->visible_los_az) . "\n";
echo "LOS El: " . round($pass->visible_los_el) . "\n";
echo "Magnitude: " . number_format($pass->max_apparent_magnitude, 1) . "\n";
echo "\n";
}
// How long did this take?
echo "Execution time: " . number_format((microtime(true) - $start) * 1000, 2) . "ms\n";
exit;
/**
* Predict the next pass.
*
* This function simply wraps the get_pass function using the current time
* as parameter.
*
* Note: the data in sat will be corrupt (future) and must be refreshed
* by the caller, if the caller will need it later on (eg. if the caller
* is GtkSatList).
*
* @param Predict_Sat $sat The satellite data.
* @param Predict_QTH $qth The observer data.
* @param int $maxdt The maximum number of days to look ahead.
*
* @return Predict_Pass Pointer instance or NULL if no pass can be
* found.
*/
public function get_next_pass(Predict_Sat $sat, Predict_QTH $qth, $maxdt)
{
/* get the current time and call the get_pass function */
$now = Predict_Time::get_current_daynum();
return $this->get_pass($sat, $qth, $now, $maxdt);
}
/** Initialise satellite data.
* @param sat The satellite to initialise.
* @param qth Optional QTH info, use (0,0) if NULL.
*
* This function calculates the satellite data at t = 0, ie. epoch time
* The function is called automatically by gtk_sat_data_read_sat.
*/
public function sat_data_init_sat(Predict_Sat $sat, Predict_QTH $qth = null)
{
$obs_geodetic = new Predict_Geodetic();
$obs_set = new Predict_ObsSet();
$sat_geodetic = new Predict_Geodetic();
/* double jul_utc, age; */
$jul_utc = Predict_Time::Julian_Date_of_Epoch($sat->tle->epoch);
// => tsince = 0.0
$sat->jul_epoch = $jul_utc;
/* initialise observer location */
if ($qth != null) {
$obs_geodetic->lon = $qth->lon * Predict::de2ra;
$obs_geodetic->lat = $qth->lat * Predict::de2ra;
$obs_geodetic->alt = $qth->alt / 1000.0;
$obs_geodetic->theta = 0;
} else {
$obs_geodetic->lon = 0.0;
$obs_geodetic->lat = 0.0;
$obs_geodetic->alt = 0.0;
$obs_geodetic->theta = 0;
}
/* execute computations */
$sdpsgp = Predict_SGPSDP::getInstance($sat);
if ($sat->flags & Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG) {
$sdpsgp->SDP4($sat, 0.0);
} else {
$sdpsgp->SGP4($sat, 0.0);
}
/* scale position and velocity to km and km/sec */
Predict_Math::Convert_Sat_State($sat->pos, $sat->vel);
/* get the velocity of the satellite */
$sat->vel->w = sqrt($sat->vel->x * $sat->vel->x + $sat->vel->y * $sat->vel->y + $sat->vel->z * $sat->vel->z);
$sat->velo = $sat->vel->w;
Predict_SGPObs::Calculate_Obs($jul_utc, $sat->pos, $sat->vel, $obs_geodetic, $obs_set);
Predict_SGPObs::Calculate_LatLonAlt($jul_utc, $sat->pos, $sat_geodetic);
while ($sat_geodetic->lon < -Predict::pi) {
$sat_geodetic->lon += Predict::twopi;
}
while ($sat_geodetic->lon > Predict::pi) {
$sat_geodetic->lon -= Predict::twopi;
}
$sat->az = Predict_Math::Degrees($obs_set->az);
$sat->el = Predict_Math::Degrees($obs_set->el);
$sat->range = $obs_set->range;
$sat->range_rate = $obs_set->range_rate;
$sat->ssplat = Predict_Math::Degrees($sat_geodetic->lat);
$sat->ssplon = Predict_Math::Degrees($sat_geodetic->lon);
$sat->alt = $sat_geodetic->alt;
$sat->ma = Predict_Math::Degrees($sat->phase);
$sat->ma *= 256.0 / 360.0;
$sat->footprint = 2.0 * Predict::xkmper * acos(Predict::xkmper / $sat->pos->w);
$age = 0.0;
$sat->orbit = floor(($sat->tle->xno * Predict::xmnpda / Predict::twopi + $age * $sat->tle->bstar * Predict::ae) * $age + $sat->tle->xmo / Predict::twopi) + $sat->tle->revnum - 1;
/* orbit type */
$sat->otype = $sat->get_orbit_type($sat);
}