@@ -36,6 +36,10 @@ discard block |
||
| 36 | 36 | /* the year 1999 assuming that two-digit years in the range 00-56 */ |
| 37 | 37 | /* correspond to 2000-2056. Until the two-line element set format */ |
| 38 | 38 | /* is changed, it is only valid for dates through 2056 December 31. */ |
| 39 | + |
|
| 40 | + /** |
|
| 41 | + * @param double $epoch |
|
| 42 | + */ |
|
| 39 | 43 | public static function Julian_Date_of_Epoch($epoch) |
| 40 | 44 | { |
| 41 | 45 | $year = 0; |
@@ -54,6 +58,11 @@ discard block |
||
| 54 | 58 | } |
| 55 | 59 | |
| 56 | 60 | /* Equivalent to the C modf function */ |
| 61 | + |
|
| 62 | + /** |
|
| 63 | + * @param double $x |
|
| 64 | + * @param integer $ipart |
|
| 65 | + */ |
|
| 57 | 66 | public static function modf($x, &$ipart) { |
| 58 | 67 | $ipart = (int)$x; |
| 59 | 68 | return $x - $ipart; |
@@ -63,6 +72,10 @@ discard block |
||
| 63 | 72 | /* of Day 0.0 of {year}. This function is used to calculate the */ |
| 64 | 73 | /* Julian Date of any date by using Julian_Date_of_Year, DOY, */ |
| 65 | 74 | /* and Fraction_of_Day. */ |
| 75 | + |
|
| 76 | + /** |
|
| 77 | + * @param integer $year |
|
| 78 | + */ |
|
| 66 | 79 | public static function Julian_Date_of_Year($year) |
| 67 | 80 | { |
| 68 | 81 | /* Astronomical Formulae for Calculators, Jean Meeus, */ |
@@ -85,6 +98,10 @@ discard block |
||
| 85 | 98 | /* 1999, as described above. The function ThetaG_JD provides the */ |
| 86 | 99 | /* same calculation except that it is based on an input in the */ |
| 87 | 100 | /* form of a Julian Date. */ |
| 101 | + |
|
| 102 | + /** |
|
| 103 | + * @param double $epoch |
|
| 104 | + */ |
|
| 88 | 105 | public static function ThetaG($epoch, Predict_DeepArg $deep_arg) |
| 89 | 106 | { |
| 90 | 107 | /* Reference: The 1992 Astronomical Almanac, page B6. */ |
@@ -157,6 +174,10 @@ discard block |
||
| 157 | 174 | /* (approximately the same as UTC) and ET (now referred to as TDT).*/ |
| 158 | 175 | /* This function is based on a least squares fit of data from 1950 */ |
| 159 | 176 | /* to 1991 and will need to be updated periodically. */ |
| 177 | + |
|
| 178 | + /** |
|
| 179 | + * @param double $year |
|
| 180 | + */ |
|
| 160 | 181 | public static function Delta_ET($year) |
| 161 | 182 | { |
| 162 | 183 | /* Values determined using data from 1950-1991 in the 1990 |
@@ -204,7 +225,7 @@ discard block |
||
| 204 | 225 | * |
| 205 | 226 | * @param Predict_TLE $tle The TLE object |
| 206 | 227 | * |
| 207 | - * @return int |
|
| 228 | + * @return double |
|
| 208 | 229 | */ |
| 209 | 230 | public static function getEpochTimeStamp(Predict_TLE $tle) |
| 210 | 231 | { |
@@ -30,194 +30,194 @@ |
||
| 30 | 30 | */ |
| 31 | 31 | class Predict_Time |
| 32 | 32 | { |
| 33 | - /* The function Julian_Date_of_Epoch returns the Julian Date of */ |
|
| 34 | - /* an epoch specified in the format used in the NORAD two-line */ |
|
| 35 | - /* element sets. It has been modified to support dates beyond */ |
|
| 36 | - /* the year 1999 assuming that two-digit years in the range 00-56 */ |
|
| 37 | - /* correspond to 2000-2056. Until the two-line element set format */ |
|
| 38 | - /* is changed, it is only valid for dates through 2056 December 31. */ |
|
| 39 | - public static function Julian_Date_of_Epoch($epoch) |
|
| 40 | - { |
|
| 41 | - $year = 0; |
|
| 42 | - |
|
| 43 | - /* Modification to support Y2K */ |
|
| 44 | - /* Valid 1957 through 2056 */ |
|
| 45 | - $day = self::modf($epoch * 1E-3, $year) * 1E3; |
|
| 46 | - if ($year < 57) { |
|
| 47 | - $year = $year + 2000; |
|
| 48 | - } else { |
|
| 49 | - $year = $year + 1900; |
|
| 50 | - } |
|
| 51 | - /* End modification */ |
|
| 52 | - |
|
| 53 | - return self::Julian_Date_of_Year($year) + $day; |
|
| 54 | - } |
|
| 55 | - |
|
| 56 | - /* Equivalent to the C modf function */ |
|
| 57 | - public static function modf($x, &$ipart) { |
|
| 58 | - $ipart = (int)$x; |
|
| 59 | - return $x - $ipart; |
|
| 60 | - } |
|
| 61 | - |
|
| 62 | - /* The function Julian_Date_of_Year calculates the Julian Date */ |
|
| 63 | - /* of Day 0.0 of {year}. This function is used to calculate the */ |
|
| 64 | - /* Julian Date of any date by using Julian_Date_of_Year, DOY, */ |
|
| 65 | - /* and Fraction_of_Day. */ |
|
| 66 | - public static function Julian_Date_of_Year($year) |
|
| 67 | - { |
|
| 68 | - /* Astronomical Formulae for Calculators, Jean Meeus, */ |
|
| 69 | - /* pages 23-25. Calculate Julian Date of 0.0 Jan year */ |
|
| 70 | - $year = $year - 1; |
|
| 71 | - $i = (int) ($year / 100); |
|
| 72 | - $A = $i; |
|
| 73 | - $i = (int) ($A / 4); |
|
| 74 | - $B = (int) (2 - $A + $i); |
|
| 75 | - $i = (int) (365.25 * $year); |
|
| 76 | - $i += (int) (30.6001 * 14); |
|
| 77 | - $jdoy = $i + 1720994.5 + $B; |
|
| 78 | - |
|
| 79 | - return $jdoy; |
|
| 80 | - } |
|
| 81 | - |
|
| 82 | - /* The function ThetaG calculates the Greenwich Mean Sidereal Time */ |
|
| 83 | - /* for an epoch specified in the format used in the NORAD two-line */ |
|
| 84 | - /* element sets. It has now been adapted for dates beyond the year */ |
|
| 85 | - /* 1999, as described above. The function ThetaG_JD provides the */ |
|
| 86 | - /* same calculation except that it is based on an input in the */ |
|
| 87 | - /* form of a Julian Date. */ |
|
| 88 | - public static function ThetaG($epoch, Predict_DeepArg $deep_arg) |
|
| 89 | - { |
|
| 90 | - /* Reference: The 1992 Astronomical Almanac, page B6. */ |
|
| 91 | - // double year,day,UT,jd,TU,GMST,_ThetaG; |
|
| 92 | - |
|
| 93 | - /* Modification to support Y2K */ |
|
| 94 | - /* Valid 1957 through 2056 */ |
|
| 95 | - $year = 0; |
|
| 96 | - $day = self::modf($epoch * 1E-3, $year) * 1E3; |
|
| 97 | - |
|
| 98 | - if ($year < 57) { |
|
| 99 | - $year += 2000; |
|
| 100 | - } else { |
|
| 101 | - $year += 1900; |
|
| 102 | - } |
|
| 103 | - /* End modification */ |
|
| 104 | - |
|
| 105 | - $UT = fmod($day, $day); |
|
| 106 | - $jd = self::Julian_Date_of_Year($year) + $day; |
|
| 107 | - $TU = ($jd - 2451545.0) / 36525; |
|
| 108 | - $GMST = 24110.54841 + $TU * (8640184.812866 + $TU * (0.093104 - $TU * 6.2E-6)); |
|
| 109 | - $GMST = Predict_Math::Modulus($GMST + Predict::secday * Predict::omega_E * $UT, Predict::secday); |
|
| 110 | - $deep_arg->ds50 = $jd - 2433281.5 + $UT; |
|
| 111 | - |
|
| 112 | - return Predict_Math::FMod2p(6.3003880987 * $deep_arg->ds50 + 1.72944494); |
|
| 113 | - } |
|
| 114 | - |
|
| 115 | - /* See the ThetaG doc block above */ |
|
| 116 | - public static function ThetaG_JD($jd) |
|
| 117 | - { |
|
| 118 | - /* Reference: The 1992 Astronomical Almanac, page B6. */ |
|
| 119 | - $UT = Predict_Math::Frac($jd + 0.5); |
|
| 120 | - $jd = $jd - $UT; |
|
| 121 | - $TU = ($jd - 2451545.0) / 36525; |
|
| 122 | - $GMST = 24110.54841 + $TU * (8640184.812866 + $TU * (0.093104 - $TU * 6.2E-6)); |
|
| 123 | - $GMST = Predict_Math::Modulus($GMST + Predict::secday * Predict::omega_E * $UT, Predict::secday); |
|
| 124 | - |
|
| 125 | - return Predict::twopi * $GMST / Predict::secday; |
|
| 126 | - } |
|
| 127 | - |
|
| 128 | - /** |
|
| 129 | - * Read the system clock and return the current Julian day. From phpPredict |
|
| 130 | - * |
|
| 131 | - * @return float |
|
| 132 | - */ |
|
| 133 | - public static function get_current_daynum() { |
|
| 134 | - // Gets the current decimal day number from microtime |
|
| 135 | - |
|
| 136 | - list($usec, $sec) = explode(' ', microtime()); |
|
| 137 | - return self::unix2daynum($sec, $usec); |
|
| 138 | - } |
|
| 139 | - |
|
| 140 | - /** |
|
| 141 | - * Converts a standard unix timestamp and optional |
|
| 142 | - * milliseconds to a daynum |
|
| 143 | - * |
|
| 144 | - * @param int $sec Seconds from the unix epoch |
|
| 145 | - * @param int $usec Optional milliseconds |
|
| 146 | - * |
|
| 147 | - * @return float |
|
| 148 | - */ |
|
| 149 | - public static function unix2daynum($sec, $usec = 0) |
|
| 150 | - { |
|
| 151 | - $time = ((($sec + $usec) / 86400.0) - 3651.0); |
|
| 152 | - return $time + 2444238.5; |
|
| 153 | - } |
|
| 154 | - |
|
| 155 | - /* The function Delta_ET has been added to allow calculations on */ |
|
| 156 | - /* the position of the sun. It provides the difference between UT */ |
|
| 157 | - /* (approximately the same as UTC) and ET (now referred to as TDT).*/ |
|
| 158 | - /* This function is based on a least squares fit of data from 1950 */ |
|
| 159 | - /* to 1991 and will need to be updated periodically. */ |
|
| 160 | - public static function Delta_ET($year) |
|
| 161 | - { |
|
| 162 | - /* Values determined using data from 1950-1991 in the 1990 |
|
| 33 | + /* The function Julian_Date_of_Epoch returns the Julian Date of */ |
|
| 34 | + /* an epoch specified in the format used in the NORAD two-line */ |
|
| 35 | + /* element sets. It has been modified to support dates beyond */ |
|
| 36 | + /* the year 1999 assuming that two-digit years in the range 00-56 */ |
|
| 37 | + /* correspond to 2000-2056. Until the two-line element set format */ |
|
| 38 | + /* is changed, it is only valid for dates through 2056 December 31. */ |
|
| 39 | + public static function Julian_Date_of_Epoch($epoch) |
|
| 40 | + { |
|
| 41 | + $year = 0; |
|
| 42 | + |
|
| 43 | + /* Modification to support Y2K */ |
|
| 44 | + /* Valid 1957 through 2056 */ |
|
| 45 | + $day = self::modf($epoch * 1E-3, $year) * 1E3; |
|
| 46 | + if ($year < 57) { |
|
| 47 | + $year = $year + 2000; |
|
| 48 | + } else { |
|
| 49 | + $year = $year + 1900; |
|
| 50 | + } |
|
| 51 | + /* End modification */ |
|
| 52 | + |
|
| 53 | + return self::Julian_Date_of_Year($year) + $day; |
|
| 54 | + } |
|
| 55 | + |
|
| 56 | + /* Equivalent to the C modf function */ |
|
| 57 | + public static function modf($x, &$ipart) { |
|
| 58 | + $ipart = (int)$x; |
|
| 59 | + return $x - $ipart; |
|
| 60 | + } |
|
| 61 | + |
|
| 62 | + /* The function Julian_Date_of_Year calculates the Julian Date */ |
|
| 63 | + /* of Day 0.0 of {year}. This function is used to calculate the */ |
|
| 64 | + /* Julian Date of any date by using Julian_Date_of_Year, DOY, */ |
|
| 65 | + /* and Fraction_of_Day. */ |
|
| 66 | + public static function Julian_Date_of_Year($year) |
|
| 67 | + { |
|
| 68 | + /* Astronomical Formulae for Calculators, Jean Meeus, */ |
|
| 69 | + /* pages 23-25. Calculate Julian Date of 0.0 Jan year */ |
|
| 70 | + $year = $year - 1; |
|
| 71 | + $i = (int) ($year / 100); |
|
| 72 | + $A = $i; |
|
| 73 | + $i = (int) ($A / 4); |
|
| 74 | + $B = (int) (2 - $A + $i); |
|
| 75 | + $i = (int) (365.25 * $year); |
|
| 76 | + $i += (int) (30.6001 * 14); |
|
| 77 | + $jdoy = $i + 1720994.5 + $B; |
|
| 78 | + |
|
| 79 | + return $jdoy; |
|
| 80 | + } |
|
| 81 | + |
|
| 82 | + /* The function ThetaG calculates the Greenwich Mean Sidereal Time */ |
|
| 83 | + /* for an epoch specified in the format used in the NORAD two-line */ |
|
| 84 | + /* element sets. It has now been adapted for dates beyond the year */ |
|
| 85 | + /* 1999, as described above. The function ThetaG_JD provides the */ |
|
| 86 | + /* same calculation except that it is based on an input in the */ |
|
| 87 | + /* form of a Julian Date. */ |
|
| 88 | + public static function ThetaG($epoch, Predict_DeepArg $deep_arg) |
|
| 89 | + { |
|
| 90 | + /* Reference: The 1992 Astronomical Almanac, page B6. */ |
|
| 91 | + // double year,day,UT,jd,TU,GMST,_ThetaG; |
|
| 92 | + |
|
| 93 | + /* Modification to support Y2K */ |
|
| 94 | + /* Valid 1957 through 2056 */ |
|
| 95 | + $year = 0; |
|
| 96 | + $day = self::modf($epoch * 1E-3, $year) * 1E3; |
|
| 97 | + |
|
| 98 | + if ($year < 57) { |
|
| 99 | + $year += 2000; |
|
| 100 | + } else { |
|
| 101 | + $year += 1900; |
|
| 102 | + } |
|
| 103 | + /* End modification */ |
|
| 104 | + |
|
| 105 | + $UT = fmod($day, $day); |
|
| 106 | + $jd = self::Julian_Date_of_Year($year) + $day; |
|
| 107 | + $TU = ($jd - 2451545.0) / 36525; |
|
| 108 | + $GMST = 24110.54841 + $TU * (8640184.812866 + $TU * (0.093104 - $TU * 6.2E-6)); |
|
| 109 | + $GMST = Predict_Math::Modulus($GMST + Predict::secday * Predict::omega_E * $UT, Predict::secday); |
|
| 110 | + $deep_arg->ds50 = $jd - 2433281.5 + $UT; |
|
| 111 | + |
|
| 112 | + return Predict_Math::FMod2p(6.3003880987 * $deep_arg->ds50 + 1.72944494); |
|
| 113 | + } |
|
| 114 | + |
|
| 115 | + /* See the ThetaG doc block above */ |
|
| 116 | + public static function ThetaG_JD($jd) |
|
| 117 | + { |
|
| 118 | + /* Reference: The 1992 Astronomical Almanac, page B6. */ |
|
| 119 | + $UT = Predict_Math::Frac($jd + 0.5); |
|
| 120 | + $jd = $jd - $UT; |
|
| 121 | + $TU = ($jd - 2451545.0) / 36525; |
|
| 122 | + $GMST = 24110.54841 + $TU * (8640184.812866 + $TU * (0.093104 - $TU * 6.2E-6)); |
|
| 123 | + $GMST = Predict_Math::Modulus($GMST + Predict::secday * Predict::omega_E * $UT, Predict::secday); |
|
| 124 | + |
|
| 125 | + return Predict::twopi * $GMST / Predict::secday; |
|
| 126 | + } |
|
| 127 | + |
|
| 128 | + /** |
|
| 129 | + * Read the system clock and return the current Julian day. From phpPredict |
|
| 130 | + * |
|
| 131 | + * @return float |
|
| 132 | + */ |
|
| 133 | + public static function get_current_daynum() { |
|
| 134 | + // Gets the current decimal day number from microtime |
|
| 135 | + |
|
| 136 | + list($usec, $sec) = explode(' ', microtime()); |
|
| 137 | + return self::unix2daynum($sec, $usec); |
|
| 138 | + } |
|
| 139 | + |
|
| 140 | + /** |
|
| 141 | + * Converts a standard unix timestamp and optional |
|
| 142 | + * milliseconds to a daynum |
|
| 143 | + * |
|
| 144 | + * @param int $sec Seconds from the unix epoch |
|
| 145 | + * @param int $usec Optional milliseconds |
|
| 146 | + * |
|
| 147 | + * @return float |
|
| 148 | + */ |
|
| 149 | + public static function unix2daynum($sec, $usec = 0) |
|
| 150 | + { |
|
| 151 | + $time = ((($sec + $usec) / 86400.0) - 3651.0); |
|
| 152 | + return $time + 2444238.5; |
|
| 153 | + } |
|
| 154 | + |
|
| 155 | + /* The function Delta_ET has been added to allow calculations on */ |
|
| 156 | + /* the position of the sun. It provides the difference between UT */ |
|
| 157 | + /* (approximately the same as UTC) and ET (now referred to as TDT).*/ |
|
| 158 | + /* This function is based on a least squares fit of data from 1950 */ |
|
| 159 | + /* to 1991 and will need to be updated periodically. */ |
|
| 160 | + public static function Delta_ET($year) |
|
| 161 | + { |
|
| 162 | + /* Values determined using data from 1950-1991 in the 1990 |
|
| 163 | 163 | Astronomical Almanac. See DELTA_ET.WQ1 for details. */ |
| 164 | 164 | |
| 165 | - $delta_et = 26.465 + 0.747622 * ($year - 1950) + |
|
| 166 | - 1.886913 * sin(Predict::twopi * ($year - 1975) / 33); |
|
| 167 | - |
|
| 168 | - return $delta_et; |
|
| 169 | - } |
|
| 170 | - |
|
| 171 | - /** |
|
| 172 | - * Converts a daynum to a unix timestamp. From phpPredict. |
|
| 173 | - * |
|
| 174 | - * @param float $dn Julian Daynum |
|
| 175 | - * |
|
| 176 | - * @return float |
|
| 177 | - */ |
|
| 178 | - public static function daynum2unix($dn) { |
|
| 179 | - // Converts a daynum to a UNIX timestamp |
|
| 180 | - |
|
| 181 | - return (86400.0 * ($dn - 2444238.5 + 3651.0)); |
|
| 182 | - } |
|
| 183 | - |
|
| 184 | - /** |
|
| 185 | - * Converts a daynum to a readable time format. |
|
| 186 | - * |
|
| 187 | - * @param float $dn The julian date |
|
| 188 | - * @param string $zone The zone string, defaults to America/Los_Angeles |
|
| 189 | - * @param string $format The date() function's format string. Defaults to m-d-Y H:i:s |
|
| 190 | - * |
|
| 191 | - * @return string |
|
| 192 | - */ |
|
| 193 | - public static function daynum2readable($dn, $zone = 'America/Los_Angeles', $format = 'm-d-Y H:i:s') |
|
| 194 | - { |
|
| 195 | - $unix = self::daynum2unix($dn); |
|
| 196 | - $date = new DateTime("@" . round($unix)); |
|
| 197 | - $dateTimezone = new DateTimezone($zone); |
|
| 198 | - $date->setTimezone($dateTimezone); |
|
| 199 | - return $date->format($format); |
|
| 200 | - } |
|
| 201 | - |
|
| 202 | - /** |
|
| 203 | - * Returns the unix timestamp of a TLE's epoch |
|
| 204 | - * |
|
| 205 | - * @param Predict_TLE $tle The TLE object |
|
| 206 | - * |
|
| 207 | - * @return int |
|
| 208 | - */ |
|
| 209 | - public static function getEpochTimeStamp(Predict_TLE $tle) |
|
| 210 | - { |
|
| 211 | - $year = $tle->epoch_year; |
|
| 212 | - $day = $tle->epoch_day; |
|
| 213 | - $sec = round(86400 * $tle->epoch_fod); |
|
| 214 | - |
|
| 215 | - $zone = new DateTimeZone('GMT'); |
|
| 216 | - $date = new DateTime(); |
|
| 217 | - $date->setTimezone($zone); |
|
| 218 | - $date->setDate($year, 1, 1); |
|
| 219 | - $date->setTime(0, 0, 0); |
|
| 220 | - |
|
| 221 | - return $date->format('U') + (86400 * $day) + $sec - 86400; |
|
| 222 | - } |
|
| 165 | + $delta_et = 26.465 + 0.747622 * ($year - 1950) + |
|
| 166 | + 1.886913 * sin(Predict::twopi * ($year - 1975) / 33); |
|
| 167 | + |
|
| 168 | + return $delta_et; |
|
| 169 | + } |
|
| 170 | + |
|
| 171 | + /** |
|
| 172 | + * Converts a daynum to a unix timestamp. From phpPredict. |
|
| 173 | + * |
|
| 174 | + * @param float $dn Julian Daynum |
|
| 175 | + * |
|
| 176 | + * @return float |
|
| 177 | + */ |
|
| 178 | + public static function daynum2unix($dn) { |
|
| 179 | + // Converts a daynum to a UNIX timestamp |
|
| 180 | + |
|
| 181 | + return (86400.0 * ($dn - 2444238.5 + 3651.0)); |
|
| 182 | + } |
|
| 183 | + |
|
| 184 | + /** |
|
| 185 | + * Converts a daynum to a readable time format. |
|
| 186 | + * |
|
| 187 | + * @param float $dn The julian date |
|
| 188 | + * @param string $zone The zone string, defaults to America/Los_Angeles |
|
| 189 | + * @param string $format The date() function's format string. Defaults to m-d-Y H:i:s |
|
| 190 | + * |
|
| 191 | + * @return string |
|
| 192 | + */ |
|
| 193 | + public static function daynum2readable($dn, $zone = 'America/Los_Angeles', $format = 'm-d-Y H:i:s') |
|
| 194 | + { |
|
| 195 | + $unix = self::daynum2unix($dn); |
|
| 196 | + $date = new DateTime("@" . round($unix)); |
|
| 197 | + $dateTimezone = new DateTimezone($zone); |
|
| 198 | + $date->setTimezone($dateTimezone); |
|
| 199 | + return $date->format($format); |
|
| 200 | + } |
|
| 201 | + |
|
| 202 | + /** |
|
| 203 | + * Returns the unix timestamp of a TLE's epoch |
|
| 204 | + * |
|
| 205 | + * @param Predict_TLE $tle The TLE object |
|
| 206 | + * |
|
| 207 | + * @return int |
|
| 208 | + */ |
|
| 209 | + public static function getEpochTimeStamp(Predict_TLE $tle) |
|
| 210 | + { |
|
| 211 | + $year = $tle->epoch_year; |
|
| 212 | + $day = $tle->epoch_day; |
|
| 213 | + $sec = round(86400 * $tle->epoch_fod); |
|
| 214 | + |
|
| 215 | + $zone = new DateTimeZone('GMT'); |
|
| 216 | + $date = new DateTime(); |
|
| 217 | + $date->setTimezone($zone); |
|
| 218 | + $date->setDate($year, 1, 1); |
|
| 219 | + $date->setTime(0, 0, 0); |
|
| 220 | + |
|
| 221 | + return $date->format('U') + (86400 * $day) + $sec - 86400; |
|
| 222 | + } |
|
| 223 | 223 | } |
@@ -42,7 +42,7 @@ discard block |
||
| 42 | 42 | |
| 43 | 43 | /* Modification to support Y2K */ |
| 44 | 44 | /* Valid 1957 through 2056 */ |
| 45 | - $day = self::modf($epoch * 1E-3, $year) * 1E3; |
|
| 45 | + $day = self::modf($epoch*1E-3, $year)*1E3; |
|
| 46 | 46 | if ($year < 57) { |
| 47 | 47 | $year = $year + 2000; |
| 48 | 48 | } else { |
@@ -55,7 +55,7 @@ discard block |
||
| 55 | 55 | |
| 56 | 56 | /* Equivalent to the C modf function */ |
| 57 | 57 | public static function modf($x, &$ipart) { |
| 58 | - $ipart = (int)$x; |
|
| 58 | + $ipart = (int) $x; |
|
| 59 | 59 | return $x - $ipart; |
| 60 | 60 | } |
| 61 | 61 | |
@@ -68,12 +68,12 @@ discard block |
||
| 68 | 68 | /* Astronomical Formulae for Calculators, Jean Meeus, */ |
| 69 | 69 | /* pages 23-25. Calculate Julian Date of 0.0 Jan year */ |
| 70 | 70 | $year = $year - 1; |
| 71 | - $i = (int) ($year / 100); |
|
| 71 | + $i = (int) ($year/100); |
|
| 72 | 72 | $A = $i; |
| 73 | - $i = (int) ($A / 4); |
|
| 73 | + $i = (int) ($A/4); |
|
| 74 | 74 | $B = (int) (2 - $A + $i); |
| 75 | - $i = (int) (365.25 * $year); |
|
| 76 | - $i += (int) (30.6001 * 14); |
|
| 75 | + $i = (int) (365.25*$year); |
|
| 76 | + $i += (int) (30.6001*14); |
|
| 77 | 77 | $jdoy = $i + 1720994.5 + $B; |
| 78 | 78 | |
| 79 | 79 | return $jdoy; |
@@ -93,7 +93,7 @@ discard block |
||
| 93 | 93 | /* Modification to support Y2K */ |
| 94 | 94 | /* Valid 1957 through 2056 */ |
| 95 | 95 | $year = 0; |
| 96 | - $day = self::modf($epoch * 1E-3, $year) * 1E3; |
|
| 96 | + $day = self::modf($epoch*1E-3, $year)*1E3; |
|
| 97 | 97 | |
| 98 | 98 | if ($year < 57) { |
| 99 | 99 | $year += 2000; |
@@ -104,12 +104,12 @@ discard block |
||
| 104 | 104 | |
| 105 | 105 | $UT = fmod($day, $day); |
| 106 | 106 | $jd = self::Julian_Date_of_Year($year) + $day; |
| 107 | - $TU = ($jd - 2451545.0) / 36525; |
|
| 108 | - $GMST = 24110.54841 + $TU * (8640184.812866 + $TU * (0.093104 - $TU * 6.2E-6)); |
|
| 109 | - $GMST = Predict_Math::Modulus($GMST + Predict::secday * Predict::omega_E * $UT, Predict::secday); |
|
| 107 | + $TU = ($jd - 2451545.0)/36525; |
|
| 108 | + $GMST = 24110.54841 + $TU*(8640184.812866 + $TU*(0.093104 - $TU*6.2E-6)); |
|
| 109 | + $GMST = Predict_Math::Modulus($GMST + Predict::secday*Predict::omega_E*$UT, Predict::secday); |
|
| 110 | 110 | $deep_arg->ds50 = $jd - 2433281.5 + $UT; |
| 111 | 111 | |
| 112 | - return Predict_Math::FMod2p(6.3003880987 * $deep_arg->ds50 + 1.72944494); |
|
| 112 | + return Predict_Math::FMod2p(6.3003880987*$deep_arg->ds50 + 1.72944494); |
|
| 113 | 113 | } |
| 114 | 114 | |
| 115 | 115 | /* See the ThetaG doc block above */ |
@@ -118,11 +118,11 @@ discard block |
||
| 118 | 118 | /* Reference: The 1992 Astronomical Almanac, page B6. */ |
| 119 | 119 | $UT = Predict_Math::Frac($jd + 0.5); |
| 120 | 120 | $jd = $jd - $UT; |
| 121 | - $TU = ($jd - 2451545.0) / 36525; |
|
| 122 | - $GMST = 24110.54841 + $TU * (8640184.812866 + $TU * (0.093104 - $TU * 6.2E-6)); |
|
| 123 | - $GMST = Predict_Math::Modulus($GMST + Predict::secday * Predict::omega_E * $UT, Predict::secday); |
|
| 121 | + $TU = ($jd - 2451545.0)/36525; |
|
| 122 | + $GMST = 24110.54841 + $TU*(8640184.812866 + $TU*(0.093104 - $TU*6.2E-6)); |
|
| 123 | + $GMST = Predict_Math::Modulus($GMST + Predict::secday*Predict::omega_E*$UT, Predict::secday); |
|
| 124 | 124 | |
| 125 | - return Predict::twopi * $GMST / Predict::secday; |
|
| 125 | + return Predict::twopi*$GMST/Predict::secday; |
|
| 126 | 126 | } |
| 127 | 127 | |
| 128 | 128 | /** |
@@ -148,7 +148,7 @@ discard block |
||
| 148 | 148 | */ |
| 149 | 149 | public static function unix2daynum($sec, $usec = 0) |
| 150 | 150 | { |
| 151 | - $time = ((($sec + $usec) / 86400.0) - 3651.0); |
|
| 151 | + $time = ((($sec + $usec)/86400.0) - 3651.0); |
|
| 152 | 152 | return $time + 2444238.5; |
| 153 | 153 | } |
| 154 | 154 | |
@@ -162,8 +162,8 @@ discard block |
||
| 162 | 162 | /* Values determined using data from 1950-1991 in the 1990 |
| 163 | 163 | Astronomical Almanac. See DELTA_ET.WQ1 for details. */ |
| 164 | 164 | |
| 165 | - $delta_et = 26.465 + 0.747622 * ($year - 1950) + |
|
| 166 | - 1.886913 * sin(Predict::twopi * ($year - 1975) / 33); |
|
| 165 | + $delta_et = 26.465 + 0.747622*($year - 1950) + |
|
| 166 | + 1.886913*sin(Predict::twopi*($year - 1975)/33); |
|
| 167 | 167 | |
| 168 | 168 | return $delta_et; |
| 169 | 169 | } |
@@ -178,7 +178,7 @@ discard block |
||
| 178 | 178 | public static function daynum2unix($dn) { |
| 179 | 179 | // Converts a daynum to a UNIX timestamp |
| 180 | 180 | |
| 181 | - return (86400.0 * ($dn - 2444238.5 + 3651.0)); |
|
| 181 | + return (86400.0*($dn - 2444238.5 + 3651.0)); |
|
| 182 | 182 | } |
| 183 | 183 | |
| 184 | 184 | /** |
@@ -193,7 +193,7 @@ discard block |
||
| 193 | 193 | public static function daynum2readable($dn, $zone = 'America/Los_Angeles', $format = 'm-d-Y H:i:s') |
| 194 | 194 | { |
| 195 | 195 | $unix = self::daynum2unix($dn); |
| 196 | - $date = new DateTime("@" . round($unix)); |
|
| 196 | + $date = new DateTime("@".round($unix)); |
|
| 197 | 197 | $dateTimezone = new DateTimezone($zone); |
| 198 | 198 | $date->setTimezone($dateTimezone); |
| 199 | 199 | return $date->format($format); |
@@ -210,7 +210,7 @@ discard block |
||
| 210 | 210 | { |
| 211 | 211 | $year = $tle->epoch_year; |
| 212 | 212 | $day = $tle->epoch_day; |
| 213 | - $sec = round(86400 * $tle->epoch_fod); |
|
| 213 | + $sec = round(86400*$tle->epoch_fod); |
|
| 214 | 214 | |
| 215 | 215 | $zone = new DateTimeZone('GMT'); |
| 216 | 216 | $date = new DateTime(); |
@@ -218,6 +218,6 @@ discard block |
||
| 218 | 218 | $date->setDate($year, 1, 1); |
| 219 | 219 | $date->setTime(0, 0, 0); |
| 220 | 220 | |
| 221 | - return $date->format('U') + (86400 * $day) + $sec - 86400; |
|
| 221 | + return $date->format('U') + (86400*$day) + $sec - 86400; |
|
| 222 | 222 | } |
| 223 | 223 | } |
@@ -6,30 +6,30 @@ |
||
| 6 | 6 | /* Common arguments between deep-space functions */ |
| 7 | 7 | class Predict_DeepArg |
| 8 | 8 | { |
| 9 | - /* Used by dpinit part of Deep() */ |
|
| 10 | - public $eosq; |
|
| 11 | - public $sinio; |
|
| 12 | - public $cosio; |
|
| 13 | - public $betao; |
|
| 14 | - public $aodp; |
|
| 15 | - public $theta2; |
|
| 16 | - public $sing; |
|
| 17 | - public $cosg; |
|
| 18 | - public $betao2; |
|
| 19 | - public $xmdot; |
|
| 20 | - public $omgdot; |
|
| 21 | - public $xnodot; |
|
| 22 | - public $xnodp; |
|
| 9 | + /* Used by dpinit part of Deep() */ |
|
| 10 | + public $eosq; |
|
| 11 | + public $sinio; |
|
| 12 | + public $cosio; |
|
| 13 | + public $betao; |
|
| 14 | + public $aodp; |
|
| 15 | + public $theta2; |
|
| 16 | + public $sing; |
|
| 17 | + public $cosg; |
|
| 18 | + public $betao2; |
|
| 19 | + public $xmdot; |
|
| 20 | + public $omgdot; |
|
| 21 | + public $xnodot; |
|
| 22 | + public $xnodp; |
|
| 23 | 23 | |
| 24 | - /* Used by dpsec and dpper parts of Deep() */ |
|
| 25 | - public $xll; |
|
| 26 | - public $omgadf; |
|
| 27 | - public $xnode; |
|
| 28 | - public $em; |
|
| 29 | - public $xinc; |
|
| 30 | - public $xn; |
|
| 31 | - public $t; |
|
| 24 | + /* Used by dpsec and dpper parts of Deep() */ |
|
| 25 | + public $xll; |
|
| 26 | + public $omgadf; |
|
| 27 | + public $xnode; |
|
| 28 | + public $em; |
|
| 29 | + public $xinc; |
|
| 30 | + public $xn; |
|
| 31 | + public $t; |
|
| 32 | 32 | |
| 33 | - /* Used by thetg and Deep() */ |
|
| 34 | - public $ds50; |
|
| 33 | + /* Used by thetg and Deep() */ |
|
| 34 | + public $ds50; |
|
| 35 | 35 | } |
@@ -3,16 +3,16 @@ |
||
| 3 | 3 | /** Brief satellite pass info. */ |
| 4 | 4 | class Predict_Pass |
| 5 | 5 | { |
| 6 | - public $satname; /*!< satellite name */ |
|
| 7 | - public $aos; /*!< AOS time in "jul_utc" */ |
|
| 8 | - public $tca; /*!< TCA time in "jul_utc" */ |
|
| 9 | - public $los; /*!< LOS time in "jul_utc" */ |
|
| 10 | - public $max_el; /*!< Maximum elevation during pass */ |
|
| 11 | - public $aos_az; /*!< Azimuth at AOS */ |
|
| 12 | - public $los_az; /*!< Azimuth at LOS */ |
|
| 13 | - public $orbit; /*!< Orbit number */ |
|
| 14 | - public $maxel_az; /*!< Azimuth at maximum elevation */ |
|
| 15 | - public $vis; /*!< Visibility string, e.g. VSE, -S-, V-- */ |
|
| 16 | - public $details = array(); /*!< List of pass_detail_t entries */ |
|
| 17 | - public $max_apparent_magnitude = null; /* maximum apparent magnitude, experimental */ |
|
| 6 | + public $satname; /*!< satellite name */ |
|
| 7 | + public $aos; /*!< AOS time in "jul_utc" */ |
|
| 8 | + public $tca; /*!< TCA time in "jul_utc" */ |
|
| 9 | + public $los; /*!< LOS time in "jul_utc" */ |
|
| 10 | + public $max_el; /*!< Maximum elevation during pass */ |
|
| 11 | + public $aos_az; /*!< Azimuth at AOS */ |
|
| 12 | + public $los_az; /*!< Azimuth at LOS */ |
|
| 13 | + public $orbit; /*!< Orbit number */ |
|
| 14 | + public $maxel_az; /*!< Azimuth at maximum elevation */ |
|
| 15 | + public $vis; /*!< Visibility string, e.g. VSE, -S-, V-- */ |
|
| 16 | + public $details = array(); /*!< List of pass_detail_t entries */ |
|
| 17 | + public $max_apparent_magnitude = null; /* maximum apparent magnitude, experimental */ |
|
| 18 | 18 | } |
@@ -3,16 +3,16 @@ |
||
| 3 | 3 | /** Brief satellite pass info. */ |
| 4 | 4 | class Predict_Pass |
| 5 | 5 | { |
| 6 | - public $satname; /*!< satellite name */ |
|
| 7 | - public $aos; /*!< AOS time in "jul_utc" */ |
|
| 8 | - public $tca; /*!< TCA time in "jul_utc" */ |
|
| 9 | - public $los; /*!< LOS time in "jul_utc" */ |
|
| 10 | - public $max_el; /*!< Maximum elevation during pass */ |
|
| 11 | - public $aos_az; /*!< Azimuth at AOS */ |
|
| 12 | - public $los_az; /*!< Azimuth at LOS */ |
|
| 13 | - public $orbit; /*!< Orbit number */ |
|
| 6 | + public $satname; /*!< satellite name */ |
|
| 7 | + public $aos; /*!< AOS time in "jul_utc" */ |
|
| 8 | + public $tca; /*!< TCA time in "jul_utc" */ |
|
| 9 | + public $los; /*!< LOS time in "jul_utc" */ |
|
| 10 | + public $max_el; /*!< Maximum elevation during pass */ |
|
| 11 | + public $aos_az; /*!< Azimuth at AOS */ |
|
| 12 | + public $los_az; /*!< Azimuth at LOS */ |
|
| 13 | + public $orbit; /*!< Orbit number */ |
|
| 14 | 14 | public $maxel_az; /*!< Azimuth at maximum elevation */ |
| 15 | - public $vis; /*!< Visibility string, e.g. VSE, -S-, V-- */ |
|
| 16 | - public $details = array(); /*!< List of pass_detail_t entries */ |
|
| 15 | + public $vis; /*!< Visibility string, e.g. VSE, -S-, V-- */ |
|
| 16 | + public $details = array(); /*!< List of pass_detail_t entries */ |
|
| 17 | 17 | public $max_apparent_magnitude = null; /* maximum apparent magnitude, experimental */ |
| 18 | 18 | } |
@@ -6,78 +6,78 @@ |
||
| 6 | 6 | /* static data for DEEP */ |
| 7 | 7 | class Predict_DeepStatic |
| 8 | 8 | { |
| 9 | - public $thgr; |
|
| 10 | - public $xnq; |
|
| 11 | - public $xqncl; |
|
| 12 | - public $omegaq; |
|
| 13 | - public $zmol; |
|
| 14 | - public $zmos; |
|
| 15 | - public $savtsn; |
|
| 16 | - public $ee2; |
|
| 17 | - public $e3; |
|
| 18 | - public $xi2; |
|
| 19 | - public $xl2; |
|
| 20 | - public $xl3; |
|
| 21 | - public $xl4; |
|
| 22 | - public $xgh2; |
|
| 23 | - public $xgh3; |
|
| 24 | - public $xgh4; |
|
| 25 | - public $xh2; |
|
| 26 | - public $xh3; |
|
| 27 | - public $sse; |
|
| 28 | - public $ssi; |
|
| 29 | - public $ssg; |
|
| 30 | - public $xi3; |
|
| 31 | - public $se2; |
|
| 32 | - public $si2; |
|
| 33 | - public $sl2; |
|
| 34 | - public $sgh2; |
|
| 35 | - public $sh2; |
|
| 36 | - public $se3; |
|
| 37 | - public $si3; |
|
| 38 | - public $sl3; |
|
| 39 | - public $sgh3; |
|
| 40 | - public $sh3; |
|
| 41 | - public $sl4; |
|
| 42 | - public $sgh4; |
|
| 43 | - public $ssl; |
|
| 44 | - public $ssh; |
|
| 45 | - public $d3210; |
|
| 46 | - public $d3222; |
|
| 47 | - public $d4410; |
|
| 48 | - public $d4422; |
|
| 49 | - public $d5220; |
|
| 50 | - public $d5232; |
|
| 51 | - public $d5421; |
|
| 52 | - public $d5433; |
|
| 53 | - public $del1; |
|
| 54 | - public $del2; |
|
| 55 | - public $del3; |
|
| 56 | - public $fasx2; |
|
| 57 | - public $fasx4; |
|
| 58 | - public $fasx6; |
|
| 59 | - public $xlamo; |
|
| 60 | - public $xfact; |
|
| 61 | - public $xni; |
|
| 62 | - public $atime; |
|
| 63 | - public $stepp; |
|
| 64 | - public $stepn; |
|
| 65 | - public $step2; |
|
| 66 | - public $preep; |
|
| 67 | - public $pl; |
|
| 68 | - public $sghs; |
|
| 69 | - public $xli; |
|
| 70 | - public $d2201; |
|
| 71 | - public $d2211; |
|
| 72 | - public $sghl; |
|
| 73 | - public $sh1; |
|
| 74 | - public $pinc; |
|
| 75 | - public $pe; |
|
| 76 | - public $shs; |
|
| 77 | - public $zsingl; |
|
| 78 | - public $zcosgl; |
|
| 79 | - public $zsinhl; |
|
| 80 | - public $zcoshl; |
|
| 81 | - public $zsinil; |
|
| 82 | - public $zcosil; |
|
| 9 | + public $thgr; |
|
| 10 | + public $xnq; |
|
| 11 | + public $xqncl; |
|
| 12 | + public $omegaq; |
|
| 13 | + public $zmol; |
|
| 14 | + public $zmos; |
|
| 15 | + public $savtsn; |
|
| 16 | + public $ee2; |
|
| 17 | + public $e3; |
|
| 18 | + public $xi2; |
|
| 19 | + public $xl2; |
|
| 20 | + public $xl3; |
|
| 21 | + public $xl4; |
|
| 22 | + public $xgh2; |
|
| 23 | + public $xgh3; |
|
| 24 | + public $xgh4; |
|
| 25 | + public $xh2; |
|
| 26 | + public $xh3; |
|
| 27 | + public $sse; |
|
| 28 | + public $ssi; |
|
| 29 | + public $ssg; |
|
| 30 | + public $xi3; |
|
| 31 | + public $se2; |
|
| 32 | + public $si2; |
|
| 33 | + public $sl2; |
|
| 34 | + public $sgh2; |
|
| 35 | + public $sh2; |
|
| 36 | + public $se3; |
|
| 37 | + public $si3; |
|
| 38 | + public $sl3; |
|
| 39 | + public $sgh3; |
|
| 40 | + public $sh3; |
|
| 41 | + public $sl4; |
|
| 42 | + public $sgh4; |
|
| 43 | + public $ssl; |
|
| 44 | + public $ssh; |
|
| 45 | + public $d3210; |
|
| 46 | + public $d3222; |
|
| 47 | + public $d4410; |
|
| 48 | + public $d4422; |
|
| 49 | + public $d5220; |
|
| 50 | + public $d5232; |
|
| 51 | + public $d5421; |
|
| 52 | + public $d5433; |
|
| 53 | + public $del1; |
|
| 54 | + public $del2; |
|
| 55 | + public $del3; |
|
| 56 | + public $fasx2; |
|
| 57 | + public $fasx4; |
|
| 58 | + public $fasx6; |
|
| 59 | + public $xlamo; |
|
| 60 | + public $xfact; |
|
| 61 | + public $xni; |
|
| 62 | + public $atime; |
|
| 63 | + public $stepp; |
|
| 64 | + public $stepn; |
|
| 65 | + public $step2; |
|
| 66 | + public $preep; |
|
| 67 | + public $pl; |
|
| 68 | + public $sghs; |
|
| 69 | + public $xli; |
|
| 70 | + public $d2201; |
|
| 71 | + public $d2211; |
|
| 72 | + public $sghl; |
|
| 73 | + public $sh1; |
|
| 74 | + public $pinc; |
|
| 75 | + public $pe; |
|
| 76 | + public $shs; |
|
| 77 | + public $zsingl; |
|
| 78 | + public $zcosgl; |
|
| 79 | + public $zsinhl; |
|
| 80 | + public $zcoshl; |
|
| 81 | + public $zsinil; |
|
| 82 | + public $zcosil; |
|
| 83 | 83 | } |
@@ -12,26 +12,26 @@ |
||
| 12 | 12 | */ |
| 13 | 13 | class Predict_PassDetail |
| 14 | 14 | { |
| 15 | - public $time; /*!< time in "jul_utc" */ |
|
| 16 | - public $pos; /*!< Raw unprocessed position at time */ |
|
| 17 | - public $vel; /*!< Raw unprocessed velocity at time */ |
|
| 18 | - public $velo; |
|
| 19 | - public $az; |
|
| 20 | - public $el; |
|
| 21 | - public $range; |
|
| 22 | - public $range_rate; |
|
| 23 | - public $lat; |
|
| 24 | - public $lon; |
|
| 25 | - public $alt; |
|
| 26 | - public $ma; |
|
| 27 | - public $phase; |
|
| 28 | - public $footprint; |
|
| 29 | - public $vis; |
|
| 30 | - public $orbit; |
|
| 15 | + public $time; /*!< time in "jul_utc" */ |
|
| 16 | + public $pos; /*!< Raw unprocessed position at time */ |
|
| 17 | + public $vel; /*!< Raw unprocessed velocity at time */ |
|
| 18 | + public $velo; |
|
| 19 | + public $az; |
|
| 20 | + public $el; |
|
| 21 | + public $range; |
|
| 22 | + public $range_rate; |
|
| 23 | + public $lat; |
|
| 24 | + public $lon; |
|
| 25 | + public $alt; |
|
| 26 | + public $ma; |
|
| 27 | + public $phase; |
|
| 28 | + public $footprint; |
|
| 29 | + public $vis; |
|
| 30 | + public $orbit; |
|
| 31 | 31 | |
| 32 | - public function __construct() |
|
| 33 | - { |
|
| 34 | - $this->pos = new Predict_Vector(); |
|
| 35 | - $this->vel = new Predict_Vector(); |
|
| 36 | - } |
|
| 32 | + public function __construct() |
|
| 33 | + { |
|
| 34 | + $this->pos = new Predict_Vector(); |
|
| 35 | + $this->vel = new Predict_Vector(); |
|
| 36 | + } |
|
| 37 | 37 | } |
@@ -12,9 +12,9 @@ |
||
| 12 | 12 | */ |
| 13 | 13 | class Predict_PassDetail |
| 14 | 14 | { |
| 15 | - public $time; /*!< time in "jul_utc" */ |
|
| 16 | - public $pos; /*!< Raw unprocessed position at time */ |
|
| 17 | - public $vel; /*!< Raw unprocessed velocity at time */ |
|
| 15 | + public $time; /*!< time in "jul_utc" */ |
|
| 16 | + public $pos; /*!< Raw unprocessed position at time */ |
|
| 17 | + public $vel; /*!< Raw unprocessed velocity at time */ |
|
| 18 | 18 | public $velo; |
| 19 | 19 | public $az; |
| 20 | 20 | public $el; |
@@ -5,8 +5,8 @@ |
||
| 5 | 5 | */ |
| 6 | 6 | class Predict_ObsSet |
| 7 | 7 | { |
| 8 | - public $az = 0.0; /*!< Azimuth [deg] */ |
|
| 9 | - public $el = 0.0; /*!< Elevation [deg] */ |
|
| 10 | - public $range = 0.0; /*!< Range [km] */ |
|
| 11 | - public $range_rate = 0.0; /*!< Velocity [km/sec] */ |
|
| 8 | + public $az = 0.0; /*!< Azimuth [deg] */ |
|
| 9 | + public $el = 0.0; /*!< Elevation [deg] */ |
|
| 10 | + public $range = 0.0; /*!< Range [km] */ |
|
| 11 | + public $range_rate = 0.0; /*!< Velocity [km/sec] */ |
|
| 12 | 12 | } |
@@ -5,8 +5,8 @@ |
||
| 5 | 5 | */ |
| 6 | 6 | class Predict_ObsSet |
| 7 | 7 | { |
| 8 | - public $az = 0.0; /*!< Azimuth [deg] */ |
|
| 9 | - public $el = 0.0; /*!< Elevation [deg] */ |
|
| 10 | - public $range = 0.0; /*!< Range [km] */ |
|
| 11 | - public $range_rate = 0.0; /*!< Velocity [km/sec] */ |
|
| 8 | + public $az = 0.0; /*!< Azimuth [deg] */ |
|
| 9 | + public $el = 0.0; /*!< Elevation [deg] */ |
|
| 10 | + public $range = 0.0; /*!< Range [km] */ |
|
| 11 | + public $range_rate = 0.0; /*!< Velocity [km/sec] */ |
|
| 12 | 12 | } |
@@ -12,221 +12,221 @@ |
||
| 12 | 12 | */ |
| 13 | 13 | class Predict_TLE |
| 14 | 14 | { |
| 15 | - public $header; /* Header line of TLE file */ |
|
| 16 | - public $line1; /* Line 1 of TLE */ |
|
| 17 | - public $line2; /* Line 2 of TLE */ |
|
| 18 | - public $epoch; /*!< Epoch Time in NORAD TLE format YYDDD.FFFFFFFF */ |
|
| 19 | - public $epoch_year; /*!< Epoch: year */ |
|
| 20 | - public $epoch_day; /*!< Epoch: day of year */ |
|
| 21 | - public $epoch_fod; /*!< Epoch: Fraction of day. */ |
|
| 22 | - public $xndt2o; /*!< 1. time derivative of mean motion */ |
|
| 23 | - public $xndd6o; /*!< 2. time derivative of mean motion */ |
|
| 24 | - public $bstar; /*!< Bstar drag coefficient. */ |
|
| 25 | - public $xincl; /*!< Inclination */ |
|
| 26 | - public $xnodeo; /*!< R.A.A.N. */ |
|
| 27 | - public $eo; /*!< Eccentricity */ |
|
| 28 | - public $omegao; /*!< argument of perigee */ |
|
| 29 | - public $xmo; /*!< mean anomaly */ |
|
| 30 | - public $xno; /*!< mean motion */ |
|
| 31 | - |
|
| 32 | - public $catnr; /*!< Catalogue Number. */ |
|
| 33 | - public $elset; /*!< Element Set number. */ |
|
| 34 | - public $revnum; /*!< Revolution Number at epoch. */ |
|
| 35 | - |
|
| 36 | - public $sat_name; /*!< Satellite name string. */ |
|
| 37 | - public $idesg; /*!< International Designator. */ |
|
| 38 | - public $status; /*!< Operational status. */ |
|
| 39 | - |
|
| 40 | - /* values needed for squint calculations */ |
|
| 41 | - public $xincl1; |
|
| 42 | - public $xnodeo1; |
|
| 43 | - public $omegao1; |
|
| 44 | - |
|
| 45 | - |
|
| 46 | - /* Converts the strings in a raw two-line element set */ |
|
| 47 | - /* to their intended numerical values. No processing */ |
|
| 48 | - /* of these values is done, e.g. from deg to rads etc. */ |
|
| 49 | - /* This is done in the select_ephemeris() function. */ |
|
| 50 | - public function __construct($header, $line1, $line2) |
|
| 51 | - { |
|
| 52 | - if (!$this->Good_Elements($line1, $line2)) { |
|
| 53 | - throw new Predict_Exception('Invalid TLE contents'); |
|
| 54 | - } |
|
| 55 | - |
|
| 56 | - $this->header = $header; |
|
| 57 | - $this->line1 = $line1; |
|
| 58 | - $this->line2 = $line2; |
|
| 59 | - |
|
| 60 | - /** Decode Card 1 **/ |
|
| 61 | - /* Satellite's catalogue number */ |
|
| 62 | - $this->catnr = (int) substr($line1, 2, 5); |
|
| 63 | - |
|
| 64 | - /* International Designator for satellite */ |
|
| 65 | - $this->idesg = substr($line1, 9, 8); |
|
| 66 | - |
|
| 67 | - /* Epoch time; this is the complete, unconverted epoch. */ |
|
| 68 | - /* Replace spaces with 0 before casting, as leading spaces are allowed */ |
|
| 69 | - $this->epoch = (float) str_replace(' ', '0', substr($line1, 18, 14)); |
|
| 70 | - |
|
| 71 | - /* Now, convert the epoch time into year, day |
|
| 15 | + public $header; /* Header line of TLE file */ |
|
| 16 | + public $line1; /* Line 1 of TLE */ |
|
| 17 | + public $line2; /* Line 2 of TLE */ |
|
| 18 | + public $epoch; /*!< Epoch Time in NORAD TLE format YYDDD.FFFFFFFF */ |
|
| 19 | + public $epoch_year; /*!< Epoch: year */ |
|
| 20 | + public $epoch_day; /*!< Epoch: day of year */ |
|
| 21 | + public $epoch_fod; /*!< Epoch: Fraction of day. */ |
|
| 22 | + public $xndt2o; /*!< 1. time derivative of mean motion */ |
|
| 23 | + public $xndd6o; /*!< 2. time derivative of mean motion */ |
|
| 24 | + public $bstar; /*!< Bstar drag coefficient. */ |
|
| 25 | + public $xincl; /*!< Inclination */ |
|
| 26 | + public $xnodeo; /*!< R.A.A.N. */ |
|
| 27 | + public $eo; /*!< Eccentricity */ |
|
| 28 | + public $omegao; /*!< argument of perigee */ |
|
| 29 | + public $xmo; /*!< mean anomaly */ |
|
| 30 | + public $xno; /*!< mean motion */ |
|
| 31 | + |
|
| 32 | + public $catnr; /*!< Catalogue Number. */ |
|
| 33 | + public $elset; /*!< Element Set number. */ |
|
| 34 | + public $revnum; /*!< Revolution Number at epoch. */ |
|
| 35 | + |
|
| 36 | + public $sat_name; /*!< Satellite name string. */ |
|
| 37 | + public $idesg; /*!< International Designator. */ |
|
| 38 | + public $status; /*!< Operational status. */ |
|
| 39 | + |
|
| 40 | + /* values needed for squint calculations */ |
|
| 41 | + public $xincl1; |
|
| 42 | + public $xnodeo1; |
|
| 43 | + public $omegao1; |
|
| 44 | + |
|
| 45 | + |
|
| 46 | + /* Converts the strings in a raw two-line element set */ |
|
| 47 | + /* to their intended numerical values. No processing */ |
|
| 48 | + /* of these values is done, e.g. from deg to rads etc. */ |
|
| 49 | + /* This is done in the select_ephemeris() function. */ |
|
| 50 | + public function __construct($header, $line1, $line2) |
|
| 51 | + { |
|
| 52 | + if (!$this->Good_Elements($line1, $line2)) { |
|
| 53 | + throw new Predict_Exception('Invalid TLE contents'); |
|
| 54 | + } |
|
| 55 | + |
|
| 56 | + $this->header = $header; |
|
| 57 | + $this->line1 = $line1; |
|
| 58 | + $this->line2 = $line2; |
|
| 59 | + |
|
| 60 | + /** Decode Card 1 **/ |
|
| 61 | + /* Satellite's catalogue number */ |
|
| 62 | + $this->catnr = (int) substr($line1, 2, 5); |
|
| 63 | + |
|
| 64 | + /* International Designator for satellite */ |
|
| 65 | + $this->idesg = substr($line1, 9, 8); |
|
| 66 | + |
|
| 67 | + /* Epoch time; this is the complete, unconverted epoch. */ |
|
| 68 | + /* Replace spaces with 0 before casting, as leading spaces are allowed */ |
|
| 69 | + $this->epoch = (float) str_replace(' ', '0', substr($line1, 18, 14)); |
|
| 70 | + |
|
| 71 | + /* Now, convert the epoch time into year, day |
|
| 72 | 72 | and fraction of day, according to: |
| 73 | 73 | |
| 74 | 74 | YYDDD.FFFFFFFF |
| 75 | 75 | */ |
| 76 | 76 | |
| 77 | - // Adjust for 2 digit year through 2056 |
|
| 78 | - $this->epoch_year = (int) substr($line1, 18, 2); |
|
| 79 | - if ($this->epoch_year > 56) { |
|
| 80 | - $this->epoch_year = $this->epoch_year + 1900; |
|
| 81 | - } else { |
|
| 82 | - $this->epoch_year = $this->epoch_year + 2000; |
|
| 83 | - } |
|
| 77 | + // Adjust for 2 digit year through 2056 |
|
| 78 | + $this->epoch_year = (int) substr($line1, 18, 2); |
|
| 79 | + if ($this->epoch_year > 56) { |
|
| 80 | + $this->epoch_year = $this->epoch_year + 1900; |
|
| 81 | + } else { |
|
| 82 | + $this->epoch_year = $this->epoch_year + 2000; |
|
| 83 | + } |
|
| 84 | 84 | |
| 85 | - /* Epoch day */ |
|
| 86 | - $this->epoch_day = (int) substr($line1, 20, 3); |
|
| 85 | + /* Epoch day */ |
|
| 86 | + $this->epoch_day = (int) substr($line1, 20, 3); |
|
| 87 | 87 | |
| 88 | - /* Epoch fraction of day */ |
|
| 89 | - $this->epoch_fod = (float) substr($line1, 23, 9); |
|
| 88 | + /* Epoch fraction of day */ |
|
| 89 | + $this->epoch_fod = (float) substr($line1, 23, 9); |
|
| 90 | 90 | |
| 91 | 91 | |
| 92 | - /* Satellite's First Time Derivative */ |
|
| 93 | - $this->xndt2o = (float) substr($line1, 33, 10); |
|
| 92 | + /* Satellite's First Time Derivative */ |
|
| 93 | + $this->xndt2o = (float) substr($line1, 33, 10); |
|
| 94 | 94 | |
| 95 | - /* Satellite's Second Time Derivative */ |
|
| 96 | - $this->xndd6o = (float) (substr($line1, 44, 1) . '.' . substr($line1, 45, 5) . 'E' . substr($line1, 50, 2)); |
|
| 95 | + /* Satellite's Second Time Derivative */ |
|
| 96 | + $this->xndd6o = (float) (substr($line1, 44, 1) . '.' . substr($line1, 45, 5) . 'E' . substr($line1, 50, 2)); |
|
| 97 | 97 | |
| 98 | - /* Satellite's bstar drag term |
|
| 98 | + /* Satellite's bstar drag term |
|
| 99 | 99 | FIXME: How about buff[0] ???? |
| 100 | 100 | */ |
| 101 | - $this->bstar = (float) (substr($line1, 53, 1) . '.' . substr($line1, 54, 5) . 'E' . substr($line1, 59, 2)); |
|
| 101 | + $this->bstar = (float) (substr($line1, 53, 1) . '.' . substr($line1, 54, 5) . 'E' . substr($line1, 59, 2)); |
|
| 102 | 102 | |
| 103 | - /* Element Number */ |
|
| 104 | - $this->elset = (int) substr($line1, 64, 4); |
|
| 103 | + /* Element Number */ |
|
| 104 | + $this->elset = (int) substr($line1, 64, 4); |
|
| 105 | 105 | |
| 106 | - /** Decode Card 2 **/ |
|
| 107 | - /* Satellite's Orbital Inclination (degrees) */ |
|
| 108 | - $this->xincl = (float) substr($line2, 8, 8); |
|
| 106 | + /** Decode Card 2 **/ |
|
| 107 | + /* Satellite's Orbital Inclination (degrees) */ |
|
| 108 | + $this->xincl = (float) substr($line2, 8, 8); |
|
| 109 | 109 | |
| 110 | - /* Satellite's RAAN (degrees) */ |
|
| 111 | - $this->xnodeo = (float) substr($line2, 17, 8); |
|
| 110 | + /* Satellite's RAAN (degrees) */ |
|
| 111 | + $this->xnodeo = (float) substr($line2, 17, 8); |
|
| 112 | 112 | |
| 113 | - /* Satellite's Orbital Eccentricity */ |
|
| 114 | - $this->eo = (float) ('.' . substr($line2, 26, 7)); |
|
| 113 | + /* Satellite's Orbital Eccentricity */ |
|
| 114 | + $this->eo = (float) ('.' . substr($line2, 26, 7)); |
|
| 115 | 115 | |
| 116 | - /* Satellite's Argument of Perigee (degrees) */ |
|
| 117 | - $this->omegao = (float) substr($line2, 34, 8); |
|
| 116 | + /* Satellite's Argument of Perigee (degrees) */ |
|
| 117 | + $this->omegao = (float) substr($line2, 34, 8); |
|
| 118 | 118 | |
| 119 | - /* Satellite's Mean Anomaly of Orbit (degrees) */ |
|
| 120 | - $this->xmo = (float) substr($line2, 43, 8); |
|
| 119 | + /* Satellite's Mean Anomaly of Orbit (degrees) */ |
|
| 120 | + $this->xmo = (float) substr($line2, 43, 8); |
|
| 121 | 121 | |
| 122 | - /* Satellite's Mean Motion (rev/day) */ |
|
| 123 | - $this->xno = (float) substr($line2, 52, 11); |
|
| 122 | + /* Satellite's Mean Motion (rev/day) */ |
|
| 123 | + $this->xno = (float) substr($line2, 52, 11); |
|
| 124 | 124 | |
| 125 | - /* Satellite's Revolution number at epoch */ |
|
| 126 | - $this->revnum = (float) substr($line2, 63, 5); |
|
| 127 | - } |
|
| 125 | + /* Satellite's Revolution number at epoch */ |
|
| 126 | + $this->revnum = (float) substr($line2, 63, 5); |
|
| 127 | + } |
|
| 128 | 128 | |
| 129 | - /* Calculates the checksum mod 10 of a line from a TLE set and */ |
|
| 130 | - /* returns true if it compares with checksum in column 68, else false.*/ |
|
| 131 | - /* tle_set is a character string holding the two lines read */ |
|
| 132 | - /* from a text file containing NASA format Keplerian elements. */ |
|
| 133 | - /* NOTE!!! The stuff about two lines is not quite true. |
|
| 129 | + /* Calculates the checksum mod 10 of a line from a TLE set and */ |
|
| 130 | + /* returns true if it compares with checksum in column 68, else false.*/ |
|
| 131 | + /* tle_set is a character string holding the two lines read */ |
|
| 132 | + /* from a text file containing NASA format Keplerian elements. */ |
|
| 133 | + /* NOTE!!! The stuff about two lines is not quite true. |
|
| 134 | 134 | The function assumes that tle_set[0] is the begining |
| 135 | 135 | of the line and that there are 68 elements - see the consumer |
| 136 | 136 | */ |
| 137 | - public function Checksum_Good($tle_set) |
|
| 138 | - { |
|
| 139 | - if (strlen($tle_set) < 69) { |
|
| 140 | - return false; |
|
| 141 | - } |
|
| 142 | - |
|
| 143 | - $checksum = 0; |
|
| 144 | - |
|
| 145 | - for ($i = 0; $i < 68; $i++) { |
|
| 146 | - if (($tle_set[$i] >= '0') && ($tle_set[$i] <= '9')) { |
|
| 147 | - $value = $tle_set[$i] - '0'; |
|
| 148 | - } else if ($tle_set[$i] == '-' ) { |
|
| 149 | - $value = 1; |
|
| 150 | - } else { |
|
| 151 | - $value = 0; |
|
| 152 | - } |
|
| 153 | - |
|
| 154 | - $checksum += $value; |
|
| 155 | - } |
|
| 156 | - |
|
| 157 | - $checksum %= 10; |
|
| 158 | - $check_digit = $tle_set[68] - '0'; |
|
| 159 | - |
|
| 160 | - return $checksum == $check_digit; |
|
| 161 | - } |
|
| 162 | - |
|
| 163 | - /* Carries out various checks on a TLE set to verify its validity */ |
|
| 164 | - /* $line1 is the first line of the TLE, $line2 is the second line */ |
|
| 165 | - /* from a text file containing NASA format Keplerian elements. */ |
|
| 166 | - public function Good_Elements($line1, $line2) |
|
| 167 | - { |
|
| 168 | - /* Verify checksum of both lines of a TLE set */ |
|
| 169 | - if (!$this->Checksum_Good($line1) || !$this->Checksum_Good($line2)) { |
|
| 170 | - return false; |
|
| 171 | - } |
|
| 172 | - |
|
| 173 | - /* Check the line number of each line */ |
|
| 174 | - if (($line1[0] != '1') || ($line2[0] != '2')) { |
|
| 175 | - return false; |
|
| 176 | - } |
|
| 177 | - |
|
| 178 | - /* Verify that Satellite Number is same in both lines */ |
|
| 179 | - if (strncmp($line1[2], $line2[2], 5) != 0) { |
|
| 180 | - return false; |
|
| 181 | - } |
|
| 182 | - |
|
| 183 | - /* Check that various elements are in the right place */ |
|
| 184 | - if (($line1[23] != '.') || |
|
| 185 | - ($line1[34] != '.') || |
|
| 186 | - ($line2[11] != '.') || |
|
| 187 | - ($line2[20] != '.') || |
|
| 188 | - ($line2[37] != '.') || |
|
| 189 | - ($line2[46] != '.') || |
|
| 190 | - ($line2[54] != '.') || |
|
| 191 | - (strncmp(substr($line1, 61), ' 0 ', 3) != 0)) { |
|
| 192 | - |
|
| 193 | - return false; |
|
| 194 | - } |
|
| 195 | - |
|
| 196 | - return true; |
|
| 197 | - } |
|
| 198 | - |
|
| 199 | - /** |
|
| 200 | - * A function to allow checksum creation of a line. This is driven by |
|
| 201 | - * the fact that some TLEs from SpaceTrack are missing checksum numbers. |
|
| 202 | - * You can use this to create a checksum for a line, but you should |
|
| 203 | - * probably have confidence that the TLE data itself is good. YMMV. |
|
| 204 | - * |
|
| 205 | - * @throws Predict_Exception if the line is not exactly 68 chars |
|
| 206 | - * @return string |
|
| 207 | - */ |
|
| 208 | - static public function createChecksum($line) |
|
| 209 | - { |
|
| 210 | - if (strlen($line) != 68) { |
|
| 211 | - throw Predict_Exception('Invalid line, needs to e 68 chars'); |
|
| 212 | - } |
|
| 213 | - |
|
| 214 | - $checksum = 0; |
|
| 215 | - |
|
| 216 | - for ($i = 0; $i < 68; $i++) { |
|
| 217 | - if (($line[$i] >= '0') && ($line[$i] <= '9')) { |
|
| 218 | - $value = (int) $line[$i]; |
|
| 219 | - } else if ($line[$i] == '-' ) { |
|
| 220 | - $value = 1; |
|
| 221 | - } else { |
|
| 222 | - $value = 0; |
|
| 223 | - } |
|
| 224 | - |
|
| 225 | - $checksum += $value; |
|
| 226 | - } |
|
| 227 | - |
|
| 228 | - $checksum %= 10; |
|
| 229 | - |
|
| 230 | - return $checksum; |
|
| 231 | - } |
|
| 137 | + public function Checksum_Good($tle_set) |
|
| 138 | + { |
|
| 139 | + if (strlen($tle_set) < 69) { |
|
| 140 | + return false; |
|
| 141 | + } |
|
| 142 | + |
|
| 143 | + $checksum = 0; |
|
| 144 | + |
|
| 145 | + for ($i = 0; $i < 68; $i++) { |
|
| 146 | + if (($tle_set[$i] >= '0') && ($tle_set[$i] <= '9')) { |
|
| 147 | + $value = $tle_set[$i] - '0'; |
|
| 148 | + } else if ($tle_set[$i] == '-' ) { |
|
| 149 | + $value = 1; |
|
| 150 | + } else { |
|
| 151 | + $value = 0; |
|
| 152 | + } |
|
| 153 | + |
|
| 154 | + $checksum += $value; |
|
| 155 | + } |
|
| 156 | + |
|
| 157 | + $checksum %= 10; |
|
| 158 | + $check_digit = $tle_set[68] - '0'; |
|
| 159 | + |
|
| 160 | + return $checksum == $check_digit; |
|
| 161 | + } |
|
| 162 | + |
|
| 163 | + /* Carries out various checks on a TLE set to verify its validity */ |
|
| 164 | + /* $line1 is the first line of the TLE, $line2 is the second line */ |
|
| 165 | + /* from a text file containing NASA format Keplerian elements. */ |
|
| 166 | + public function Good_Elements($line1, $line2) |
|
| 167 | + { |
|
| 168 | + /* Verify checksum of both lines of a TLE set */ |
|
| 169 | + if (!$this->Checksum_Good($line1) || !$this->Checksum_Good($line2)) { |
|
| 170 | + return false; |
|
| 171 | + } |
|
| 172 | + |
|
| 173 | + /* Check the line number of each line */ |
|
| 174 | + if (($line1[0] != '1') || ($line2[0] != '2')) { |
|
| 175 | + return false; |
|
| 176 | + } |
|
| 177 | + |
|
| 178 | + /* Verify that Satellite Number is same in both lines */ |
|
| 179 | + if (strncmp($line1[2], $line2[2], 5) != 0) { |
|
| 180 | + return false; |
|
| 181 | + } |
|
| 182 | + |
|
| 183 | + /* Check that various elements are in the right place */ |
|
| 184 | + if (($line1[23] != '.') || |
|
| 185 | + ($line1[34] != '.') || |
|
| 186 | + ($line2[11] != '.') || |
|
| 187 | + ($line2[20] != '.') || |
|
| 188 | + ($line2[37] != '.') || |
|
| 189 | + ($line2[46] != '.') || |
|
| 190 | + ($line2[54] != '.') || |
|
| 191 | + (strncmp(substr($line1, 61), ' 0 ', 3) != 0)) { |
|
| 192 | + |
|
| 193 | + return false; |
|
| 194 | + } |
|
| 195 | + |
|
| 196 | + return true; |
|
| 197 | + } |
|
| 198 | + |
|
| 199 | + /** |
|
| 200 | + * A function to allow checksum creation of a line. This is driven by |
|
| 201 | + * the fact that some TLEs from SpaceTrack are missing checksum numbers. |
|
| 202 | + * You can use this to create a checksum for a line, but you should |
|
| 203 | + * probably have confidence that the TLE data itself is good. YMMV. |
|
| 204 | + * |
|
| 205 | + * @throws Predict_Exception if the line is not exactly 68 chars |
|
| 206 | + * @return string |
|
| 207 | + */ |
|
| 208 | + static public function createChecksum($line) |
|
| 209 | + { |
|
| 210 | + if (strlen($line) != 68) { |
|
| 211 | + throw Predict_Exception('Invalid line, needs to e 68 chars'); |
|
| 212 | + } |
|
| 213 | + |
|
| 214 | + $checksum = 0; |
|
| 215 | + |
|
| 216 | + for ($i = 0; $i < 68; $i++) { |
|
| 217 | + if (($line[$i] >= '0') && ($line[$i] <= '9')) { |
|
| 218 | + $value = (int) $line[$i]; |
|
| 219 | + } else if ($line[$i] == '-' ) { |
|
| 220 | + $value = 1; |
|
| 221 | + } else { |
|
| 222 | + $value = 0; |
|
| 223 | + } |
|
| 224 | + |
|
| 225 | + $checksum += $value; |
|
| 226 | + } |
|
| 227 | + |
|
| 228 | + $checksum %= 10; |
|
| 229 | + |
|
| 230 | + return $checksum; |
|
| 231 | + } |
|
| 232 | 232 | } |
@@ -12,30 +12,30 @@ discard block |
||
| 12 | 12 | */ |
| 13 | 13 | class Predict_TLE |
| 14 | 14 | { |
| 15 | - public $header; /* Header line of TLE file */ |
|
| 16 | - public $line1; /* Line 1 of TLE */ |
|
| 17 | - public $line2; /* Line 2 of TLE */ |
|
| 18 | - public $epoch; /*!< Epoch Time in NORAD TLE format YYDDD.FFFFFFFF */ |
|
| 15 | + public $header; /* Header line of TLE file */ |
|
| 16 | + public $line1; /* Line 1 of TLE */ |
|
| 17 | + public $line2; /* Line 2 of TLE */ |
|
| 18 | + public $epoch; /*!< Epoch Time in NORAD TLE format YYDDD.FFFFFFFF */ |
|
| 19 | 19 | public $epoch_year; /*!< Epoch: year */ |
| 20 | - public $epoch_day; /*!< Epoch: day of year */ |
|
| 21 | - public $epoch_fod; /*!< Epoch: Fraction of day. */ |
|
| 22 | - public $xndt2o; /*!< 1. time derivative of mean motion */ |
|
| 23 | - public $xndd6o; /*!< 2. time derivative of mean motion */ |
|
| 24 | - public $bstar; /*!< Bstar drag coefficient. */ |
|
| 25 | - public $xincl; /*!< Inclination */ |
|
| 26 | - public $xnodeo; /*!< R.A.A.N. */ |
|
| 27 | - public $eo; /*!< Eccentricity */ |
|
| 28 | - public $omegao; /*!< argument of perigee */ |
|
| 29 | - public $xmo; /*!< mean anomaly */ |
|
| 30 | - public $xno; /*!< mean motion */ |
|
| 31 | - |
|
| 32 | - public $catnr; /*!< Catalogue Number. */ |
|
| 33 | - public $elset; /*!< Element Set number. */ |
|
| 34 | - public $revnum; /*!< Revolution Number at epoch. */ |
|
| 35 | - |
|
| 36 | - public $sat_name; /*!< Satellite name string. */ |
|
| 37 | - public $idesg; /*!< International Designator. */ |
|
| 38 | - public $status; /*!< Operational status. */ |
|
| 20 | + public $epoch_day; /*!< Epoch: day of year */ |
|
| 21 | + public $epoch_fod; /*!< Epoch: Fraction of day. */ |
|
| 22 | + public $xndt2o; /*!< 1. time derivative of mean motion */ |
|
| 23 | + public $xndd6o; /*!< 2. time derivative of mean motion */ |
|
| 24 | + public $bstar; /*!< Bstar drag coefficient. */ |
|
| 25 | + public $xincl; /*!< Inclination */ |
|
| 26 | + public $xnodeo; /*!< R.A.A.N. */ |
|
| 27 | + public $eo; /*!< Eccentricity */ |
|
| 28 | + public $omegao; /*!< argument of perigee */ |
|
| 29 | + public $xmo; /*!< mean anomaly */ |
|
| 30 | + public $xno; /*!< mean motion */ |
|
| 31 | + |
|
| 32 | + public $catnr; /*!< Catalogue Number. */ |
|
| 33 | + public $elset; /*!< Element Set number. */ |
|
| 34 | + public $revnum; /*!< Revolution Number at epoch. */ |
|
| 35 | + |
|
| 36 | + public $sat_name; /*!< Satellite name string. */ |
|
| 37 | + public $idesg; /*!< International Designator. */ |
|
| 38 | + public $status; /*!< Operational status. */ |
|
| 39 | 39 | |
| 40 | 40 | /* values needed for squint calculations */ |
| 41 | 41 | public $xincl1; |
@@ -93,12 +93,12 @@ discard block |
||
| 93 | 93 | $this->xndt2o = (float) substr($line1, 33, 10); |
| 94 | 94 | |
| 95 | 95 | /* Satellite's Second Time Derivative */ |
| 96 | - $this->xndd6o = (float) (substr($line1, 44, 1) . '.' . substr($line1, 45, 5) . 'E' . substr($line1, 50, 2)); |
|
| 96 | + $this->xndd6o = (float) (substr($line1, 44, 1).'.'.substr($line1, 45, 5).'E'.substr($line1, 50, 2)); |
|
| 97 | 97 | |
| 98 | 98 | /* Satellite's bstar drag term |
| 99 | 99 | FIXME: How about buff[0] ???? |
| 100 | 100 | */ |
| 101 | - $this->bstar = (float) (substr($line1, 53, 1) . '.' . substr($line1, 54, 5) . 'E' . substr($line1, 59, 2)); |
|
| 101 | + $this->bstar = (float) (substr($line1, 53, 1).'.'.substr($line1, 54, 5).'E'.substr($line1, 59, 2)); |
|
| 102 | 102 | |
| 103 | 103 | /* Element Number */ |
| 104 | 104 | $this->elset = (int) substr($line1, 64, 4); |
@@ -111,7 +111,7 @@ discard block |
||
| 111 | 111 | $this->xnodeo = (float) substr($line2, 17, 8); |
| 112 | 112 | |
| 113 | 113 | /* Satellite's Orbital Eccentricity */ |
| 114 | - $this->eo = (float) ('.' . substr($line2, 26, 7)); |
|
| 114 | + $this->eo = (float) ('.'.substr($line2, 26, 7)); |
|
| 115 | 115 | |
| 116 | 116 | /* Satellite's Argument of Perigee (degrees) */ |
| 117 | 117 | $this->omegao = (float) substr($line2, 34, 8); |
@@ -145,7 +145,7 @@ discard block |
||
| 145 | 145 | for ($i = 0; $i < 68; $i++) { |
| 146 | 146 | if (($tle_set[$i] >= '0') && ($tle_set[$i] <= '9')) { |
| 147 | 147 | $value = $tle_set[$i] - '0'; |
| 148 | - } else if ($tle_set[$i] == '-' ) { |
|
| 148 | + } else if ($tle_set[$i] == '-') { |
|
| 149 | 149 | $value = 1; |
| 150 | 150 | } else { |
| 151 | 151 | $value = 0; |
@@ -216,7 +216,7 @@ discard block |
||
| 216 | 216 | for ($i = 0; $i < 68; $i++) { |
| 217 | 217 | if (($line[$i] >= '0') && ($line[$i] <= '9')) { |
| 218 | 218 | $value = (int) $line[$i]; |
| 219 | - } else if ($line[$i] == '-' ) { |
|
| 219 | + } else if ($line[$i] == '-') { |
|
| 220 | 220 | $value = 1; |
| 221 | 221 | } else { |
| 222 | 222 | $value = 0; |
@@ -6,8 +6,8 @@ |
||
| 6 | 6 | */ |
| 7 | 7 | class Predict_Geodetic |
| 8 | 8 | { |
| 9 | - public $lat; /*!< Lattitude [rad] */ |
|
| 10 | - public $lon; /*!< Longitude [rad] */ |
|
| 11 | - public $alt; /*!< Altitude [km] */ |
|
| 12 | - public $theta; |
|
| 9 | + public $lat; /*!< Lattitude [rad] */ |
|
| 10 | + public $lon; /*!< Longitude [rad] */ |
|
| 11 | + public $alt; /*!< Altitude [km] */ |
|
| 12 | + public $theta; |
|
| 13 | 13 | } |
@@ -8,32 +8,32 @@ |
||
| 8 | 8 | */ |
| 9 | 9 | class Predict_SGSDPStatic |
| 10 | 10 | { |
| 11 | - public $aodp; |
|
| 12 | - public $aycof; |
|
| 13 | - public $c1; |
|
| 14 | - public $c4; |
|
| 15 | - public $c5; |
|
| 16 | - public $cosio; |
|
| 17 | - public $d2; |
|
| 18 | - public $d3; |
|
| 19 | - public $d4; |
|
| 20 | - public $delmo; |
|
| 21 | - public $omgcof; |
|
| 22 | - public $eta; |
|
| 23 | - public $omgdot; |
|
| 24 | - public $sinio; |
|
| 25 | - public $xnodp; |
|
| 26 | - public $sinmo; |
|
| 27 | - public $t2cof; |
|
| 28 | - public $t3cof; |
|
| 29 | - public $t4cof; |
|
| 30 | - public $t5cof; |
|
| 31 | - public $x1mth2; |
|
| 32 | - public $x3thm1; |
|
| 33 | - public $x7thm1; |
|
| 34 | - public $xmcof; |
|
| 35 | - public $xmdot; |
|
| 36 | - public $xnodcf; |
|
| 37 | - public $xnodot; |
|
| 38 | - public $xlcof; |
|
| 11 | + public $aodp; |
|
| 12 | + public $aycof; |
|
| 13 | + public $c1; |
|
| 14 | + public $c4; |
|
| 15 | + public $c5; |
|
| 16 | + public $cosio; |
|
| 17 | + public $d2; |
|
| 18 | + public $d3; |
|
| 19 | + public $d4; |
|
| 20 | + public $delmo; |
|
| 21 | + public $omgcof; |
|
| 22 | + public $eta; |
|
| 23 | + public $omgdot; |
|
| 24 | + public $sinio; |
|
| 25 | + public $xnodp; |
|
| 26 | + public $sinmo; |
|
| 27 | + public $t2cof; |
|
| 28 | + public $t3cof; |
|
| 29 | + public $t4cof; |
|
| 30 | + public $t5cof; |
|
| 31 | + public $x1mth2; |
|
| 32 | + public $x3thm1; |
|
| 33 | + public $x7thm1; |
|
| 34 | + public $xmcof; |
|
| 35 | + public $xmdot; |
|
| 36 | + public $xnodcf; |
|
| 37 | + public $xnodot; |
|
| 38 | + public $xlcof; |
|
| 39 | 39 | } |