Project01| 气象站点数据处理|降水量缺测值插补替代

需要注意年份的索引以及闰年平年的天数…
1981年至2010年共10957天

import numpy as np
import pandas as pd
import time
import math
import warnings
warnings.filterwarnings("ignore")

pre_txt=np.loadtxt(r"F:\09_climate_site_data\03_日降水集中指数\02_data\pre-day.txt")
print(pre_txt)
pre_txt=np.array(pre_txt)

print(pre_txt.shape)
pre_txt1=[]
for i in range(pre_txt.shape[0]):
    for j in range(pre_txt.shape[1]):
        if pre_txt[i,j]!=32766:
            pre_txt1.append(pre_txt[i,j]*0.1)
        else:
            pre_txt1.append(32766)
pre_txt1=np.array(pre_txt1)
pre_txt2=pre_txt1.reshape(1672,19359).T
print(pre_txt2.shape)
pre_cols_txt=np.loadtxt(r"F:\09_climate_site_data\id-value.txt")
pre_cols=np.array