#?-*-?coding:?utf-8?-*-

import?numpy?as?np
import?time
from?scipy?import?optimize
from?matplotlib?import?pyplot?as?plt
from?mpl_toolkits.mplot3d?import?Axes3D
from?datetime?import?datetime?timedelta


#?把從Rinex讀取的字符串日期中的年，月，日。。提取出來
def?out_YMD（str）:
????Y?=?int（str[0:4]）
????M?=?int（str[5:7]）
????D?=?int（str[8:10]）
????hour?=?int（str[11:13]）
????minute?=?int（str[14:16]）
????sec?=?float（str[17:19]）
????return?（Y?M?D?hour?minute?sec）


#?格里高利里轉換成儒略日
def?d_to_mjd（tuple）:
????Y?M?D?hour?minute?sec?=?tuple
????D?+=?hour?/?24?+?minute?/?24?/?60?+?sec?/?24?/?3600
????if?M?>?2:
????????new_Y?=?Y
????????new_M?=?M
????else:
????????new_Y?=?Y?-?1
????????new_M?=?M?+?12
????A?=?int（new_Y?/?100）
????B?=?2?-?A?+?int（A?/?4）
????JD?=?int（365.25?*?（new_Y?+?4716））?+?int（30.6001?*?（new_M?+?1））?+?D?+?B?-?1524.5
????return?JD?-?2400000.5


#?由約化儒略日計算星期幾
def?mjd_wd（mjd）:
????wd?=?（（mjd?-?54899）?%?7）?+?1
????return?wd?-?7?if?（wd?>?7）?else?wd
#?應該再創建一個函數，把上面的幾個函數寫在一起，直接一步調用了

class?Renix_Read（object）:
????“““創建一個讀取文件的類，暫時先將讀取的各類值設為
????屬性#時間間隔單位為秒“““

????def?__init__（self?filename?interval）:
????????self.interval?=?interval
????????self.filename?=?filename
????????self.gravitational_constant?=?3.986005e14
????????self.We?=?7.29211567e-5
????????self.Prn?=?[]
????????self.Sqrt_a?=?[]
????????self.Delta_n?=?[]
????????self.a0?=?[]
????????self.a1?=?[]
????????self.a2?=?[]
????????self.Toc?=?[]
????????self.Toe?=?[]
????????self.M0?=?[]
????????self.Mk?=?[]
????????self.E?=?[]
????????self.W?=?[]
????????self.Cuc?=?[]
????????self.Crc?=?[]
????????self.Cic?=?[]
????????self.Cus?=?[]
????????self.Crs?=?[]
????????self.Cis?=?[]
????????self.I0?=?[]
????????self.Idot?=?[]
????????self.Omega?=?[]
????????self.Omega_dot?=?[]
????????self.time_observed?=?[]
????????self.length?=?0



????def?copy_n_times（self?list）:???#后來加上的，把數據重復n次
????????“““廣播星歷需要外推，根據要計算的時間間隔，那么就要重復多少次“““
????????new_list?=?[]
????????for?x?in?list:
????????????for?i?in?range（7200?//?self.interval）:
????????????????new_list.append（x）
????????return?new_list


????def?observe_time（self?alist）:
????????new_list?=?list（map（lambda?str_time:?datetime.strptime（str_time?‘%Y?%m?%d?%H?%M?%S‘）?-?timedelta（
????????????seconds=（3600?//?self.interval?-?1）?*?self.interval）?alist））
????????new_new_list?=?new_list[:]
????????for?j?in?range（self.length）:
????????????for?i?in?range（int（7200?//?self.interval））:
????????????????new_new_list[j?*?（7200?//?self.interval）?+?i]?=?new_list[j?*?（7200?//?self.interval）?+?i]?+?timedelta（
????????????????????seconds=self.interval?*?i）
????????str_new_list?=?list（map（lambda?time_obj:?datetime.strftime（time_obj?‘%Y?%m?%d?%H?%M?%S‘）?new_new_list））
????????return?str_new_list

????def?open_file（self）:
????????PRN_list?a0_list?a1_list?a2_list?Iode_list?Crs_list?=?[]?[]?[]?[]?[]?[]
????????Delta_n_list?M0_list?Cuc_list?E_list?Cus_list?=?[]?[]?[]?[]?[]