import copy
from datetime import datetime
from datetime import timedelta
from dateutil import parser
import matplotlib
matplotlib.use('TKAgg')
import matplotlib.pyplot as plt
import numpy as np
from numpy import mean
from tkinter import *
import time
import math
satellite=["sat0","sat1","sat2","sat3","sat4","sat5","sat6","sat7","sat8"]
address_name=["Abidjan","Accra","Asmara","Balikpapan","Bozeman","Crystal-Lake","Elk-City","Hanahan","Iquique",
"Lhasa","Madras","Manaus","Mangalore","Minsk","Munchen","Parole","Recife","Rio-de-Janeiro",
"Sarajevo","Wallingford-Center","Yaounde","北京","南宁","武汉"]
address_position={"Abidjan":[-4.02+360,5.32],"Accra":[-0.2+360,5.56],"Asmara":[15.33,38.97],"Balikpapan":[116.86,-1.21],
"Bozeman":[-111.04+360,45.68],"Crystal-Lake":[-88.33+360,42.23],"Elk-City":[-99.41+360,35.4],"Hanahan":[-80.01+360,32.93],
"Iquique":[-69.91+360,-20.26],"Lhasa":[91.13,29.65],"Madras":[80.25,13.06],"Manaus":[-60.01+360,-3.12],
"Mangalore":[75.16,12.95],"Minsk":[27.58,53.9],"Munchen":[11.54,48.14],"Parole":[-76.55+360,38.98],
"Recife":[-34.91+360,-8.09],"Rio-de-Janeiro":[-43.46+360,-22.72],"Sarajevo":[18.43,43.87],"Wallingford-Center":[-72.82+360,41.45],
"Yaounde":[11.51,3.87],"北京":[116.39,39.91],"南宁":[108.33,22.8],"武汉":[114.28,30.57]}
dict_circle={"sat0":[],"sat1":[],"sat2":[],"sat3":[],"sat4":[],"sat5":[],"sat6":[],"sat7":[],"sat8":[]}#每个卫星覆盖的圆，【时间，经度，纬度】
address_timewindow={"Abidjan":[],"Accra":[],"Asmara":[],"Balikpapan":[],"Bozeman":[],"Crystal-Lake":[],"Elk-City":[],"Hanahan":[],"Iquique":[],
"Lhasa":[],"Madras":[],"Manaus":[],"Mangalore":[],"Minsk":[],"Munchen":[],"Parole":[],"Recife":[],"Rio-de-Janeiro":[],"Sarajevo":[],"Wallingford-Center":[],"Yaounde":[],"北京":[],"南宁":[],"武汉":[]}#每个地点对应的窗口时间
address_single_timewindow={}
address_double_timewindow = {"Abidjan": [], "Accra": [], "Asmara": [], "Balikpapan": [], "Bozeman": [],
"Crystal-Lake": [], "Elk-City": [], "Hanahan": [], "Iquique": [],
"Lhasa": [], "Madras": [], "Manaus": [], "Mangalore": [], "Minsk": [], "Munchen": [],
"Parole": [], "Recife": [], "Rio-de-Janeiro": [], "Sarajevo": [],
"Wallingford-Center": [], "Yaounde": [], "北京": [], "南宁": [], "武汉": []}
address_freetime = {"Abidjan": [], "Accra": [], "Asmara": [], "Balikpapan": [], "Bozeman": [], "Crystal-Lake": [],
"Elk-City": [], "Hanahan": [], "Iquique": [],
"Lhasa": [], "Madras": [], "Manaus": [], "Mangalore": [], "Minsk": [], "Munchen": [],
"Parole": [], "Recife": [], "Rio-de-Janeiro": [], "Sarajevo": [], "Wallingford-Center": [],
"Yaounde": [], "北京": [], "南宁": [], "武汉": []}#间隙
max_freetime=[]#最大时间间隙
ave_freetime=[]#平均时间间隙
def getshijiancha(str1,str2):
#datetime.strptime是日期形式，相减之后变成datedelta时间差形式
#输入str，返回时间差
timetmp1=datetime.strptime(str(datetime.strptime(str1, '%Y/%m/%d,%H:%M:%S')-datetime.strptime(str2, '%Y/%m/%d,%H:%M:%S')),'%H:%M:%S')
return(timetmp1.hour * 3600 + timetmp1.minute * 60 + timetmp1.second)
def date2second(str):
str=datetime.strptime(str, '%Y/%m/%d,%H:%M:%S')
return (str.hour * 3600 + str.minute * 60 + str.second)
def dataprocessing():#传入相关的卫星和地点
global address_single_timewindow
for sat in dict_circle:
time = datetime(2022, 1, 1, 0, 0, 0) # 起始时间
file = open('Data/SatelliteInfo/SatCoverInfo_' + sat[3] + '.txt')
lines = file.readlines()
i = 0
while (i < len(lines)):
point = [] # 时间，圆心精度，圆心维度，半径
min1 = lines[i + 6].split() # 前为经度最小值，后为纬度
max1 = lines[i + 16].split() # 前为经度最大值，后为纬度
r = (float(max1[0]) - float(min1[0])) / 2
center = r + float(min1[0]) # 圆心的经度
point.append(time.strftime("%Y/%m/%d,%H:%M:%S"))
point.append(center)
point.append(float(min1[1]))
point.append(r)
dict_circle[sat].append(point)
time = time + timedelta(seconds=1) # 时间+1s
i += 22
# break#!!!!!!!!!!!!!!!!!!!!!!！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！1测试点，记得删掉
start = 0
end = 0
# 地点经度：address_position[address][0]，地点纬度：address_position[address][1]
# 卫星覆盖圆心：dict_circle[circle][t][1],dict_circle[circle][t][2]
# 半径：dict_circle[circle][t][3]
for address in address_position:
for circle in dict_circle:
for t in range(0, 86401):
if (start == 0): # 还没找到开始的时间窗口,则加入开始时间
if ((address_position[address][0] - dict_circle[circle][t][1]) ** 2 +
(address_position[address][1] - dict_circle[circle][t][2]) ** 2) <= 49.07:
start = dict_circle[circle][t][0]
elif ((start != 0) and (end == 0)): # 还没找到结束的时间窗口，则加入结束时间
if ((address_position[address][0] - dict_circle[circle][t][1]) ** 2 +
(address_position[address][1] - dict_circle[circle][t][2]) ** 2) > 49.07:
end = dict_circle[circle][t][0]
address_timewindow[address].append([circle, start, end])
start = end = 0
# print(address_timewindow["Abidjan"])
# 计算时间窗口
# repeat：重叠时间段
repeat = {"Abidjan": [], "Accra": [], "Asmara": [], "Balikpapan": [], "Bozeman": [], "Crystal-Lake": [],
"Elk-City": [], "Hanahan": [], "Iquique": [],
"Lhasa": [], "Madras": [], "Manaus": [], "Mangalore": [], "Minsk": [], "Munchen": [], "Parole": [],
"Recife": [], "Rio-de-Janeiro": [], "Sarajevo": [], "Wallingford-Center": [], "Yaounde": [], "北京": [],
"南宁": [], "武汉": []}
tmp = 0
for address in address_timewindow:
for i in range(0, len(address_timewindow[address])):
for t in range(i + 1, len(address_timewindow[address])):
if ((address_timewindow[address][t][2] <= address_timewindow[address][i][1]) or (
address_timewindow[address][t][1] >= address_timewindow[address][i][2])):
continue
repeat[address].append([i, t])#记录重叠的时间窗口的序号，以配对形式记录
# for address in address_timewindow:
# address_single_timewindow[address]=address_timewindow[address]
address_single_timewindow=copy.deepcopy(address_timewindow)
jiaoji = []
# for address in repeat:
# print(repeat[address])
for address in repeat:
if (repeat[address]):
# 按照重叠卫星的第二个降序排列，这样先删除后面的元素不会影响前面的元素
repeat[address] = sorted(repeat[address], key=lambda x: x[1], reverse=True)
for timedui in repeat[address]:
# print(type(address_single_timewindow[address][timedui[0]][1]))
# print(address_single_timewindow[address][timedui[0]][1])
# 求并集
address_single_timewindow[address][timedui[0]][1] = min(address_timewindow[address][timedui[0]][1],
