public function testSign() { $this->assertEquals('-', Time::sec(-10)->sign); $this->assertEquals('-', Time::min(-10)->sign); $this->assertEquals('-', Time::hours(-10)->sign); $this->assertEquals('-', Time::days(-10)->sign); $this->assertEquals('-', Time::weeks(-10)->sign); $this->assertEquals('-', Time::hms(-10, 0, 0, 0)->sign); $this->assertEquals('-', Time::hms(0, -10, 0, 0)->sign); $this->assertEquals('-', Time::hms(0, 0, -10, 0)->sign); $this->assertEquals('-', Time::hms(0, 0, 0, -10)->sign); $this->assertEquals('+', Time::hms(10, 0, 0, 0)->sign); $this->assertEquals('+', Time::hms(1, -10, 0, 0)->sign); $this->assertEquals('+', Time::hms(1, 0, -10, 0)->sign); $this->assertEquals('+', Time::hms(1, 0, 0, -10)->sign); $this->assertEquals('+', Time::sec(10)->sign); $this->assertEquals('+', Time::min(10)->sign); $this->assertEquals('+', Time::hours(10)->sign); $this->assertEquals('+', Time::days(10)->sign); $this->assertEquals('+', Time::weeks(10)->sign); }
/** * Calculates the nutation in Right Ascension (equation of the equinoxes) for * a given date * * @param AstroDate $date Date to calculate * @return Time Nutation in Right Ascension as time */ public static function inRA(AstroDate $date) { // Copy time as UTC $utc = $date->copy()->toUTC(); // Get the nutations and true obliquity of the date $n = static::find($utc); $ε = static::trueObliquity($utc); // Calculate the nutation in right ascension as seconds of time $nRA = rad2deg($n->long->rad * cos($ε->rad) / 15) * Time::SEC_IN_HOUR; return Time::sec($nRA)->round(4); }
public function test() { return; $frame = Frame::ICRF(); $epoch = Epoch::J(2000); echo $epoch->toDate()->toTDB(); echo "\n"; echo "\n"; exit; $ra = \Marando\Units\Time::hms(11, 16, 46.6); $dec = \Marando\Units\Angle::dms(5, 45, 32.5); $dist = Distance::au(5.8); echo "\n" . ($eq = new Equat($frame, $epoch, $ra, $dec, $dist)); echo "\n" . ($eq = (new Equat($frame, $epoch, $ra, $dec, $dist))->apparent()); echo "\n"; return; $frame = Frame::ICRF(); $epoch = AstroDate::jd(2457335.472615741)->toEpoch(); $ra = \Marando\Units\Time::hms(11, 16, 46.6); $dec = \Marando\Units\Angle::dms(5, 45, 32.5); $dist = Distance::au(5.8); echo "\n" . ($eq = new Equat($frame, $epoch, $ra, $dec, $dist)); echo "\n" . ($eq = (new Equat($frame, $epoch, $ra, $dec, $dist))->apparent()); $eq = new Equat($frame, $epoch, $ra, $dec, $dist); $eq->topo = Geo::deg(27, -82); echo "\n" . $eq->toHoriz(); echo "\n" . $eq->apparent(); echo "\n" . ($h = $eq->toHoriz()); echo "\n" . $h->az->deg . "\t" . $h->alt->deg; echo "\n" . $eq->apparent(\Marando\Units\Pressure::mbar(100), \Marando\Units\Temperature::F(90), 0.85); echo "\n" . $eq->toHoriz(); return; $ra = \Marando\Units\Time::hours(14.424354); $dec = \Marando\Units\Angle::deg(33.54366); $dist = Distance::au(1.5); $geo = Geo::deg(27, -82); $e = new Equat(Frame::ICRF(), Epoch::jd(2455586), $ra, $dec, $dist); $e->topo = $geo; echo "\n\n" . $e; echo "\n" . $e->apparent(); $e = new Equat(Frame::ICRF(), Epoch::jd(2455586), $ra, $dec, $dist); $e->topo = $geo; echo "\n" . $e->apparent(); $e = new Equat(Frame::ICRF(), Epoch::jd(2455586), $ra, $dec, $dist); $e->topo = $geo; echo "\n" . $e->toHoriz(); return; $e = new Equat(Frame::ICRF(), Epoch::jd(2455586), $ra, $dec, $dist); echo "\n\n" . $e; echo "\n" . $e->apparent($geo); echo "\n" . $e->apparent()->toHoriz($geo); $e = new Equat(Frame::ICRF(), Epoch::jd(2455586), $ra, $dec, $dist); echo "\n\n" . $e; echo "\n" . $e->apparent(); echo "\n" . $e->apparent()->toHoriz(); $e = new Equat(Frame::ICRF(), Epoch::jd(2451586), $ra, $dec, $dist); $e->topo = $geo; echo "\n\n" . $e; echo "\n" . $e->apparent(); echo "\n" . $e->apparent()->toHoriz(); return; // Position of Mercury $x = Distance::au(+1.18052679326447); $y = Distance::au(-0.3650485652522116); $z = Distance::au(-0.2123422968928603); // Velocity of Mercury $vx = Velocity::aud(0.02052320553396918); $vy = Velocity::aud(0.02828864020900079); $vz = Velocity::aud(0.01145246106968459); // Frame & epoch $frame = Frame::ICRF(); $epoch = AstroDate::parse('2015-Mar-20')->toEpoch(); $c = new Cartesian($frame, $epoch, $x, $y, $z, $vx, $vy, $vz); echo "\n\n"; echo $c; echo "\n\n"; echo $c->setUnit('km km/d'); echo "\n\n"; echo $c->setUnit('km km/s'); echo "\n\n"; echo $c->toEquat(); echo "\n\n"; echo $c->toEquat()->apparent(); echo "\n\n"; echo $c->toEquat()->apparent(Geo::deg(27, 278)); echo "\n\n"; echo $c->toEquat()->toHoriz(Geo::deg(27, 278)); echo "\n\n"; }
/** * @covers Marando\AstroDate\AstroDate::diff */ public function testDiff() { $tests = [[AstroDate::jd(2451545.5), AstroDate::jd(2451555.5), Time::days(10)], [AstroDate::jd(2451545.5), AstroDate::jd(2451556.0), Time::days(10.5)], [AstroDate::jd(2451545.5), AstroDate::jd(2451595.5), Time::days(50)]]; foreach ($tests as $t) { $dt1 = $t[0]; $dt2 = $t[1]; $diff = $t[2]; $this->assertEquals($diff, $dt1->diff($dt2)); } }
/** * Calculates the time required to travel the provided distance at the * velocity of this instance. * * @param Distance $dist * * @return Time */ public function time(Distance $dist) { return Time::sec($dist->m / $this->dist->m * $this->time->sec); }
/** * Performs a [IRCS -> observed] coordinate transformation for the parameters * of this instance * @param string $type Coordintate type, 'e' for equat 'h' for horiz * @return static|Horiz */ protected function ICRStoObserved($type = 'e', Pressure $pressure = null, Temperature $temp = null, $humidity = null) { // Instance initial properties $rc = $this->ra->toAngle()->rad; $dc = $this->dec->rad; $date1 = $this->epoch->toDate()->toTDB()->toJD(); $pr = 0; $pd = 0; $rv = 0; $px = $this->dist->au > 0 ? 8.794 / 3600 / $this->dist->au : 0; $utc1 = $this->epoch->toDate()->toUTC()->toJD(); $dut1 = IERS::jd($utc1)->dut1(); $elong = $this->topo ? $this->topo->lon->rad : 0; $phi = $this->topo ? $this->topo->lat->rad : 0; $hm = 0; //$this->obsrv->height->m; $xp = IERS::jd($utc1)->x() / 3600 * pi() / 180; $yp = IERS::jd($utc1)->y() / 3600 * pi() / 180; $phpa = $pressure ? $pressure->mbar : 0; $tc = $temp ? $temp->c : 0; $rh = $humidity ? $humidity : 0; $wl = 0.55; // ICRS -> CIRS (geocentric observer) IAU::Atci13($rc, $dc, $pr, $pd, $px, $rv, $date1, 0, $ri, $di, $eo); // CIRS -> ICRS (astrometric) //IAU::Atic13($ri, $di, $date1, 0, $rca, $dca, $eo); // ICRS (astrometric) -> CIRS (geocentric observer) //IAU::Atci13($rca, $dca, $pr, $pd, $px, $rv, $date1, 0, $ri, $di, $eo); // Apparent place ? //$ri = $ri - $eo; //$di = $di; // // CIRS -> topocentric IAU::Atio13($ri, $di, $utc1, 0, $dut1, $elong, $phi, $hm, $xp, $yp, $phpa, $tc, $rh, $wl, $aob, $zob, $hob, $dob, $rob); if ($type == 'e') { // Copy this instance, and override the apparent RA and Decl $topocentric = $this->copy(); $topocentric->ra = Angle::rad($rob)->toTime(); $topocentric->dec = Angle::rad($dob); $topocentric->apparent = true; // Return apparent coordinates return $topocentric; } else { // Prepare new horizontal instance $horiz = new Horiz(Angle::rad(deg2rad(90) - $zob), Angle::rad($aob), $this->dist); return $horiz; } }
public function test() { return; // Converts to J2000.0 echo "\n" . AstroDate::parse('2000-Jan-1 12:00:00 TT')->toEpoch(); // Isn't J2000.0 so represented as date echo "\n" . ($e = AstroDate::parse('2000-Jan-1 14:00:00 TT')->toEpoch()); echo "\n\n" . AstroDate::now()->toEpoch(); $d = AstroDate::now()->setTimezone('est'); echo "\n\n" . $d->sinceMidnight(); echo "\n\n"; return; echo "\n" . ($d = AstroDate::now()->setTimezone(TimeZone::UT(-7))); echo "\n" . TimeScale::TT()->abr; echo "\n" . TimeScale::TT()->name; echo "\n" . Epoch::J2000(); echo "\n" . Epoch::J1900(); echo "\n" . Epoch::B1950(); echo "\n" . Epoch::B1900(); echo "\n" . Epoch::B(1954.4423); echo "\n" . Epoch::jd(2455200.5); echo "\n\n" . Epoch::J2000()->toDate(); echo "\n" . Epoch::J1900()->toDate(); echo "\n" . Epoch::B1950()->toDate(); echo "\n" . Epoch::B1900()->toDate(); echo "\n" . Epoch::B(1954.4423)->toDate(); echo "\n\n" . AstroDate::now(); echo "\n" . ($d = AstroDate::now()->setTimezone(TimeZone::parse('EST'))); echo "\n" . $d->toTDB(); echo "\n" . $d->toUTC(); echo "\n" . ($d = AstroDate::now()->setTimezone('MST')); echo "\n" . ($d = AstroDate::now()->setTimezone('PST')); return; echo "\n" . AstroDate::parse('2015-Dec-10 6:00'); echo "\n" . AstroDate::parse('2015-Dec-10 6:00')->setTimezone(TimeZone::parse('EST')); echo "\n" . AstroDate::now(); echo "\n" . AstroDate::now()->setTimezone(TimeZone::parse('EST')); echo "\n" . AstroDate::now()->setTimezone(TimeZone::parse('PST')); echo "\n" . AstroDate::now()->setTimezone(TimeZone::UT(6)); echo "\n" . TimeZone::parse('est')->offset; echo "\n" . TimeZone::parse('est')->offset(2451545.5); echo "\n" . TimeZone::parse('est')->offset(2451589.5); return; echo "\n" . TimeZone::UT(-2); echo "\n" . TimeZone::parse('EST'); echo "\n" . AstroDate::now()->setTimezone(TimeZone::EST())->format(AstroDate::FORMAT_GOOGLE); echo "\n" . AstroDate::now()->setTimezone(TimeZone::UT(-10.5))->format(AstroDate::FORMAT_GOOGLE); echo "\n" . AstroDate::now()->setTimezone(TimeZone::UT(-10))->format(AstroDate::FORMAT_GOOGLE); return; $t = timezone_abbreviations_list(); $s = timezone_identifiers_list(); $e = timezone_name_from_abbr('est'); var_dump($t); return; $d = AstroDate::now()->setTimezone(TimeZone::EST()); echo "\n" . $d->format('r Y-M-c h:i:s.u A T'); echo "\n" . $d->sub(Time::days(15))->format(AstroDate::FORMAT_GENERIC); echo "\n" . $d->sidereal('a'); echo "\n" . $d->sidereal('m'); echo "\n" . $d->format(AstroDate::FORMAT_JPL_FRAC); echo "\n" . $d->sidereal('a', Angle::deg(-82.47)); echo "\n" . $d->sidereal('m', Angle::deg(-82.47)); echo "\n" . $d->sinceMidnight(); echo "\n" . $d->untilMidnight(); /** * TODO: * * - Figure out formula for DST start and end * - Figure out formula for Week # of year * * */ return; echo "\n" . $d->format(DateTime::RSS); echo "\n" . $d; echo "\n" . $d->format(AstroDate::FORMAT_GENERIC); echo "\n" . $d->format(DateTime::ISO8601); return; $str = '2016-Nov-14 17:07:07.120'; echo "\n" . ($d = AstroDate::parse($str)); echo "\n" . ($d = AstroDate::parse($str)->setTimezone(TimeZone::EST())); echo "\n" . $d->format('Z T P O e u s i h H G g a A y Y n M m F W z D Y-m-d j'); $str = '2016-Nov-14 17:07:07.120 TAI'; echo "\n" . ($d = AstroDate::parse($str)); echo "\n" . $d->format('O e u s i h H G g a A y Y n M m F W z D Y-m-d j'); $str = '2015-Nov-1 17:07:07.120 UTC'; $d = AstroDate::parse($str); echo "\n" . $d->format('z w S N l L D Y-m-d j'); return; echo "\n\n" . ($str = '2015-Nov-16 17:07:07.120 TT'); echo "\n" . AstroDate::parse($str); echo "\n\n" . ($str = '-1950-1-16 17:07:07 UTC'); echo "\n" . ($d = AstroDate::parse($str)); echo "\n" . $d->format('Y-m-d'); return; echo "\n" . ($d = AstroDate::now()); echo "\n" . $d->toUT1(); echo "\n" . $d->toTAI(); echo "\n\n" . ($d = AstroDate::now(TimeZone::EST())); echo "\n" . $d->toUT1(); echo "\n" . $d->toTAI(); return; $d = new AstroDate(2017, 11, 15, 7, 0, 0); echo "\n" . $d; echo "\n" . $d->setTimezone(TimeZone::EST()); //var_dump($d); $d = new AstroDate(2017, 6, 15, 7, 0, 0); echo "\n\n" . $d; echo "\n" . $d->setTimezone(TimeZone::EST()); //var_dump($d); $d = new AstroDate(2017, 11, 15, 7, 0, 0); echo "\n\n" . $d; echo "\n" . $d->setTimezone(TimeZone::EST()); //var_dump($d); $d = new AstroDate(2017, 6, 15, 6, 0, 0); echo "\n\n" . $d; echo "\n" . $d->setTimezone(TimeZone::EST()); //var_dump($d); return; $d = new AstroDate(2017, 11, 15, 0, 0, 0); echo "\n" . $d; echo "\n" . $d->jd(); echo "\n" . $d->jd(12); echo "\n" . $d->toTT(); echo "\n" . $d->toJD(); echo "\n" . $d->jd(12); echo "\n"; $d = new AstroDate(2015, 11, 15, 20, 23, 18, TimeZone::EST(), TimeScale::TT()); echo "\n" . $d; echo "\n" . $d->setTimezone(TimeZone::UTC()); echo "\n" . $d->toUT1(); echo "\n" . $d->setTimezone(TimeZone::EST()); echo "\n" . $d->jd(); echo "\n" . $d->jd(12); echo "\n" . $d->mjd(); echo "\n" . $d->mjd(12); echo "\n" . $d->monthName(); echo "\n" . $d->mjd(); echo "\n\n" . $d->toTAI(); echo "\n" . $d->year; echo "\n" . $d->month; echo "\n" . $d->day; echo "\n" . $d->hour; echo "\n" . $d->min; echo "\n" . $d->sec; echo "\n" . $d->micro; echo "\n" . $d->timezone; echo "\n" . $d->timescale; $d->year = 2020; echo "\n" . $d; $d->month = 12.13432; echo "\n" . $d; $d->sec = 12; echo "\n" . $d; echo "\n" . $d->add(Time::sec(0.5)); echo "\n" . $d->setDateTime(2017, 9, 4, 18, 34, 34.234); var_dump($d->isLeapYear()); echo "\n" . $d->setDateTime(2016, 9, 4, 18, 34, 34.234); var_dump($d->isLeapYear()); echo "\n" . $d->dayName(); $a = new AstroDate(2016, 12, 25, 19, 10, 2); $b = new AstroDate(2015, 12, 25, 19, 10, 2); echo "\n" . $a->diff($b); echo "\n" . ($a = new AstroDate(2016, 12, 31, 19, 10, 2)); echo "\n" . ($b = new AstroDate(2015, 12, 31, 19, 10, 2)); echo "\n" . $a->dayOfYear(); echo "\n" . $b->dayOfYear(); echo "\n\n" . $b; echo "\n" . $b->toTDB(); echo "\n" . $b->setTimezone(TimeZone::UTC()); echo "\n" . $b->toTDB(); echo "\n" . $b->setTimezone(TimeZone::EST()); echo "\n" . $b->toTDB(); echo "\n\n\n"; $d = new AstroDate(2015, 11, 15, 20, 23, 18.454334); echo "\n" . $d; echo "\n" . $d->setTimezone(TimeZone::EST()); echo "\n" . $d->toUTC(); echo "\n" . $d->toTAI(); echo "\n" . $d->setTimezone(TimeZone::EST()); echo "\n" . $d->toTT(); return; $d = new AstroDate(2015, 11, 15, 20, 23, 18.454334); echo "\n" . $d; $d = new AstroDate(2015, 11, 15, 20, 23, 18.454334, TimeZone::EST()); echo "\n" . $d; echo "\n" . $d->setDate(2020, 10, 2); echo "\n" . $d->setTime(23, 59, 1); echo "\n" . $d->setTimezone(TimeZone::UTC()); echo "\n" . $d->setTimezone(TimeZone::EST()); }
/** * Returns the time left until since midnight * @return Time */ public function untilMidnight() { return Time::days(1 - $this->dayFrac)->setUnit('hours'); }
/** * Calculates the position and velocity vector of a target with respect to a * given center accounting for light-travel time * * @param SSObj $target Target body * @param SSObj $center Center body * @param Time $lt Target -> center light travel time * * @return array Position/velocity vector in AU and AU/day */ public function observe(SSObj $target, SSObj $center, Time &$lt = null) { // Initial JDE $jde = $this->jde; // Light-travel time (days) and target distance (au) $τ = 0; $Δ = 0; // Maximum convergence iterations $maxIter = 100; // Converge on light travel time for ($i = 0; $i < $maxIter; $i++) { // Calculate position, distance and light travel time $pv = $this->jde($jde - $τ)->position($target, $center); $Δ0 = sqrt($pv[0] * $pv[0] + $pv[1] * $pv[1] + $pv[2] * $pv[2]); $τ = 0.0057755183 * $Δ0; // Check if distance equals last distance if ($Δ0 != $Δ) { // Doesn't equal push and re-iterate $Δ = $Δ0; } else { // Equal, so break $Δ = $Δ0; break; } } // Store light travel time as Time instance $lt = Time::days($τ); // Return the apparent position return $this->jde($jde - $τ)->position($target, $center); }
public function testDist() { $v = Velocity::mph(70); $this->assertEquals(35, $v->dist(Time::min(30))->mi); }
public function testTime() { $angle = Angle::time(Time::hours(12)); $this->assertEquals(180, $angle->deg); }