当前位置：首页 > wzjs >正文

网站开发技术考试题目长沙网络营销群

wzjs 2025/9/22 22:22:13

网站开发技术考试题目,长沙网络营销群,网站开发多长时间,中国企业信用信息公示网查询系统对CH报文数据进行每小时插值原始数据文件数据三小时一次的报文数据需求按小时补齐热带气旋路径信息插值后数据效果如下： 插值代码 # 对ch文件插值import pandas as pd import datetime import osdef interpolate_ch_one_hour (file_name):new_file_name…

对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)