1.一种对加密自动气象站实测温度数值进行质量检查的方法,其特征在于所述方法包括如下步骤:
步骤1)构建矩形地理区域Area_AWS,将Area_AWS区域所包含的自动站作为处理对象;
步骤2)以GI(t)为格距,将Area_AWS划分为若干个子网格,GI(t)∈{GI1,GI2,…,GIz},且GI1<GI2<…<GIz,t表示子网划分的组数,z表示子网划分的最大组数,t∈{1,2,…z};
每个子网格用Sub_Gird(GI(t),x,y)表示,其中,x、y分别对应地表示各个子网格的横向坐标与纵向坐标;
步骤3)定义一个数据集DataSet_Sub[x,y,Count_Station(GI(t),x,y),Sum_Temp(GI(t),x,y),Avg_Temp(GI(t),x,y),SD_Temp(GI(t),x,y)]用于记录各个子网格内包含的自动站数量和温度情况,其中,Count_Station(GI(t),x,y)表示某一子网格Sub_Gird(GI(t),x,y)中的自动站数量,Sum_Temp(GI(t),x,y)表示子网格Sub_Gird(GI(t),x,y)中所有自动站实测温度的累加和,Avg_Temp(GI(t),x,y)表示子网格Sub_Gird(GI(t),x,y)中所有自动站实测温度的算术平均值,SD_Temp(GI(t),x,y)表示子网格Sub_Gird(GI(t),x,y)中所有自动站实测温度的标准差,遍历数据集DataSet_Sub中的每条记录,将各个自动站根据自动站对应的经纬度归并到相应的子网格Sub_Gird(GI(t),x,y)中,并计算不同格矩GI(t)的偏差指数CI(k);
步骤4)以子网格中包含的自动站数量为自变量,以偏差指数CI(k)为因变量构建一个离散数据集,对该离散数据集进行函数拟合,得到一个反映子网格自动站数量与偏差指数变化关系的函数C(z),通过函数C(z),得到单位子网格内只有1个自动站时的偏差指数,即C(1);
步骤5)由所述CI(1)乘以子网格内自动站的平均温度,即得到反映该子网格内所有自动站观测温度的平均误差的参量,将所述参量引入原有的温度数值的判定式形成改进的温度数值的判定式,如式(1),根据改进的温度数值的判定式对自动站实测温度进行空间一致性检查:Va‑m×(TH+Va×C(1))≤Vt≤Va+n×(TH+Va×C(1)) (1)
其中,Vt表示自动站StationX的实测值,Va表示与StationX地理位置相邻的若干个自动站的实测记录平均值,TH表示判定阈值,m、n都是自然数;当Vt满足式(1)的条件时,判定该数值可信或正确,反之,则判定Vt可疑或错误;
所述步骤3)中计算不同格矩GI(t)的偏差指数具体包括如下步骤:
步骤3‑1)根据式(2)、式(3)逐一计算各个子网格的平均温度Avg_Temp(GI(t),其中,AWS_Temp(GI(t),x,y,i)表示在GI(t)为格距划分的子网格Sub_Gird(GI(t),x,y)中某一自动站i实测的温度数值,i∈[1,Count_Station(GI(t),x,y)],且为整数;步骤3‑
2)根据式(4)逐一计算各个子网格温度的标准差SD_Temp(GI(t),x,y),
步骤3‑3)根据所述数据集DataSet_Sub中的数值通过式(5)逐一计算各个子网格的温度偏差系数,
步骤3‑4)对所有子网格,以子网格包含的自动站数量Count_Station为归并基准,根据式(6)计算不同站点数量情况下,子网格的温度偏差系数平均值Avg_BS(GI(t),k),其中,KK=Max(Count_Station(GI(t),x,y)),Max()表示求最大值的函数,k表示子网格中包含的自动站数量,且k∈[0,KK];
步骤3‑5)以各子网格包含的自动站数量为归并基准,计算不同格矩GI(t)的偏差指数CI(k),计算方法为:
2.根据权利要求1所述的对加密自动气象站实测温度数值进行质量检查的方法,其特征在于所述步骤1)包括如下步骤:步骤1‑1)读取自动站数据文件AWS_File,从AWS_File中筛选出对温度数据进行质量控制需要的信息,形成数据集DataSet[SID,Lo,La,Temp],其中SID表示自动站编号,Lo表示自动站经度,La表示自动站纬度,Temp表示自动站实测温度;
步骤1‑2)从数据集DataSet中提取自动站经度最小值Lo_Min、自动站经度最大值Lo_Max、自动站纬度最小值La_Min以及自动站纬度最大值La_Max,根据所述自动站经度最小值、自动站经度最大值、自动站纬度最小值和自动站纬度最大值构建矩形地理区域Area_AWS。
3.根据权利要求2所述的对加密自动气象站实测温度数值进行质量检查的方法,其特征在于所述步骤2)中,x和y都是整数,且满足: