做网站组织架构站长工具pr值查询

对CH报文数据进行每小时插值

原始数据文件

数据

三小时一次的报文数据

需求

按小时补齐热带气旋路径信息
插值后数据效果如下：

插值代码

# 对ch文件插值import pandas as pd
import datetime
import osdef interpolate_ch_one_hour (file_name):new_file_name=file_name.split('.')[0]+'_new.txt'with open(file_name,'r') as f:content=f.readlines()def write_line(line):if not line.endswith('\n'):line+='\n'try:with open(new_file_name,'r') as f:content=f.readlines()if line == content[-1]:print(line+' is already in the file!')return  # 重复行不写入except Exception as e:pass# print('开始写入数据:',line)with open(new_file_name,'a') as f:f.write(line)# 每个台风路径插值for i in range(len(content)):if content[i].startswith('66666'):write_line(content[i])  # 写入第一行for j in range(i+1,len(content)-1):if content[j+1].startswith('66666'):write_line(content[j])  # 写入中间行breakelse:st_line=content[j].split(' ')et_line=content[j+1].split(' ')if st_line[2]=='':st_line=st_line[:2]+st_line[3:]if et_line[2]=='':et_line=et_line[:2]+et_line[3:]if st_line[3]=='':st_line=st_line[:3]+st_line[4:]if et_line[3]=='':et_line=et_line[:3]+et_line[4:]# 输入前后时间和经纬度，等小时间隔插值st_time = st_line[0]st_lat = float(st_line[2])st_lon = float(st_line[3])et_time = et_line[0]et_lat = float(et_line[2])try:et_lon = float(et_line[3])except Exception as e:print(e)print(et_line)print(et_line[3])print(j)time_interval = (datetime.datetime.strptime(et_time,'%Y%m%d%H')-datetime.datetime.strptime(st_time,'%Y%m%d%H')).seconds/60/60lat_interval = (et_lat-st_lat)/time_intervallon_interval = (et_lon-st_lon)/time_intervalfor k in range(int(time_interval)):time_now = datetime.datetime.strptime(st_time,'%Y%m%d%H')+pd.Timedelta(minutes=k*60)lat_now = st_lat+lat_interval*klon_now = st_lon+lon_interval*kline_now = time_now.strftime('%Y%m%d%H')+' '+st_line[1]+' '+str(int(lat_now))+' '+str(int(lon_now))+' '+' '.join(st_line[4:])+''write_line(line_now)write_line(content[-1])if __name__ == '__main__':list_file=['CH2021BST.txt']  #　修改成需要插值的文件名，支持多个文件# list_file=['CH2023BST.txt','CH2022BST.txt','CH2021BST.txt']for file_name in list_file:interpolate_ch_one_hour(file_name)