高空探测数据处理--中心级质控（三）

最新推荐文章于 2024-04-26 19:59:23 发布

nobrody

最新推荐文章于 2024-04-26 19:59:23 发布

阅读量150

点赞数

文章标签：算法 c++ 开发语言

本文链接：https://blog.csdn.net/zshluckydogs/article/details/128753462

版权

接上。

4.8东向风质控

东向风质控包括东向风速离群值检查。

vector<qcFlag> espQC(vector<float> esp,int length,int outside,int inside){
    qcFlag temp;
    QCMethod qcm;
    vector<int> code;
    vector<qcFlag> result;
    int m2;
    qcm.N = length;
    code = qcm.ESPD_Double_weighted_outlier_check(esp);
    for (int i=0;i<length;i++){
        
        m2 = code[i];
        if (m2 > 1)
            temp.code = "2";
        else if (m2 == 1)
            temp.code = "1";
        else
            temp.code = "0";
        //cout<<m2<<" ";
        
        switch (m2)
        {
        case 0:{
            temp.label = "";
            break;
        }
        case 1:{
            if (i<outside)
                temp.label += "A16_ESPD_1,";
            if (i>=outside && i< inside)
                temp.label += "B16_ESPD_1,";
            if (i>inside)
                temp.label += "C16_ESPD_1,";
            break;
            }
        case 2:{
            if (i<outside)
                temp.label += "A16_ESPD_2,";
            if (i>=outside && i< inside)
                temp.label += "B16_ESPD_2,";
            if (i>inside)
                temp.label += "C16_ESPD_2,";
            break;
            }
        }
        result.push_back(temp);  
    }
    return result;
}

4.9垂直速度检查

垂直速度检查主要包括双权重离群值检查

vector<qcFlag> vspQC(vector<float> vsp,int length,int outside,int inside){
    qcFlag temp;
    QCMethod qcm;
    vector<int> code;
    vector<qcFlag> result;
    int m2;
    qcm.N = length;
    code = qcm.Vsp_Double_weighted_outlier_check(vsp);
    for (int i=0;i<length;i++){
        
        m2 = code[i];
        if (m2 > 1)
            temp.code = "2";
        else if (m2 == 1)
            temp.code = "1";
        else
            temp.code = "0";
        //cout<<m2<<" ";
        
        switch (m2)
        {
        case 0:{
            temp.label = "";
            break;
        }
        case 1:{
            if (i<outside)
                temp.label += "A23_VSPD_1,";
            if (i>=outside && i< inside)
                temp.label += "B23_VSPD_1,";
            if (i>inside)
                temp.label += "C23_VSPD_1,";
            break;
            }
        case 2:{
            if (i<outside)
                temp.label += "A23_VSPD_2,";
            if (i>=outside && i< inside)
                temp.label += "B23_VSPD_2,";
            if (i>inside)
                temp.label += "C23_VSPD_2,";
            break;
            }
        }
        result.push_back(temp);  
    }
    return result;
}

4.10风速/风向质控

风速/风向质控包括风速极值检查、风速风向一致性检查、风切变检查等。

vector<qcFlag> windSpeedQC(vector<float> wsd,vector<float> wd,vector<float> pre,int length,int outside,int inside){
    qcFlag temp;
    QCMethod qcm;
    vector<int> code;
    vector<qcFlag> result;
    int m1,m2,m3;
    qcm.N = length;
    
    for (int i=0;i<length-1;i++){
        m1 = qcm.WindSpeed_extreme_check(wsd[i],pre[i]);
        m2 = qcm.Windspeed_WindDirection_consistency(wsd[i],wd[i]);
        m3 = qcm.WindSheer(wsd[i],wsd[i+1],wd[i]);
        
        if (m1+m2+m3 > 1)
            temp.code = "22";
        else if (m1+m2+m3 == 1)
            temp.code = "11";
        else
            temp.code = "00";
        //cout<<m2<<" ";
        if (m1==0){
            temp.label = "";
            }
        if (m1==1){
            if (i<outside)
                temp.label = "A7_SPD_1,A7_DIR_1,";
            if (i>=outside && i< inside)
                temp.label = "B7_SPD_1,B7_DIR_1,";
            if (i>inside)
                temp.label = "C7_SPD_1,C7_DIR_1,";
            }
        if(m1==2){
            if (i<outside)
                temp.label = "A7_SPD_2,A7_DIR_2,";
            if (i>=outside && i< inside)
                temp.label = "B7_SPD_2,B7_DIR_2,";
            if (i>inside)
                temp.label = "C7_SPD_2,C7_DIR_2,";
            
            }
        switch (m2)
        {
        case 0:{
            temp.label += "";
            break;
        }
        case 1:{
            if (i<outside)
                temp.label += "A10_SPD_1,A10_DIR_1,";
            if (i>=outside && i< inside)
                temp.label += "B10_SPD_1,B10_DIR_1,";
            if (i>inside)
                temp.label += "C10_SPD_1,C10_DIR_1,";
            break;
            }
        case 2:{
            if (i<outside)
                temp.label += "A10_SPD_2,A10_DIR_2,";
            if (i>=outside && i< inside)
                temp.label += "B10_SPD_2,B10_DIR_2,";
            if (i>inside)
                temp.label += "C10_SPD_2,C10_DIR_2,";
            break;
            }
        }
        switch (m3)
        {
        case 0:{
            temp.label += "";
            break;
        }
        case 1:{
            if (i<outside)
                temp.label += "A12_SPD_1,A12_DIR_1,";
            if (i>=outside && i< inside)
                temp.label += "B12_SPD_1,B12_DIR_1,";
            if (i>inside)
                temp.label += "C12_SPD_1,C12_DIR_1,";
            break;
            }
        case 2:{
            if (i<outside)
                temp.label += "A12_SPD_2,A12_DIR_2,";
            if (i>=outside && i< inside)
                temp.label += "B12_SPD_2,B12_DIR_2,";
            if (i>inside)
                temp.label += "C12_SPD_2,C12_DIR_2,";
            break;
            }
        }
        result.push_back(temp);  
    }
    temp.label = "";
    temp.code = "00";
    result.push_back(temp);
    return result;
}

4.11经度质控

经度质控包括经度极值检查。

vector<qcFlag> lonQC(vector<float> lon,float sta_lon,int length,int outside,int inside){
    qcFlag temp;
    QCMethod qcm;
    vector<int> code;
    vector<qcFlag> result;
    int m1;
    qcm.N = length;
    
    for (int i=0;i<length;i++){
        m1 = qcm.Lon_extreme_check(lon[i],sta_lon);
        
        if (m1 > 1)
            temp.code = "2";
        else if (m1 == 1)
            temp.code = "1";
        else
            temp.code = "0";
        //cout<<m2<<" ";
        if (m1==0){
            temp.label = "";
            }
        if (m1==1){
            if (i<outside)
                temp.label = "A4_LNG_1,";
            if (i>=outside && i< inside)
                temp.label = "B4_LNG_1,";
            if (i>inside)
                temp.label = "C4_LNG_1";
            }
        if(m1==2){
            if (i<outside)
                temp.label = "A4_LNG_2,";
            if (i>=outside && i< inside)
                temp.label = "B4_LNG_2,";
            if (i>inside)
                temp.label = "C4_LNG_2,";
            
            }
        
        result.push_back(temp);  
    }
    return result;
}

4.12纬度质控

纬度质控主要是包括纬度的极值检查

vector<qcFlag> latQC(vector<float> lat,float sta_lat,int length,int outside,int inside){
    qcFlag temp;
    QCMethod qcm;
    vector<int> code;
    vector<qcFlag> result;
    int m1;
    qcm.N = length;
    
    for (int i=0;i<length;i++){
        m1 = qcm.Lat_extreme_check(lat[i],sta_lat);
        
        if (m1 > 1)
            temp.code = "2";
        else if (m1 == 1)
            temp.code = "1";
        else
            temp.code = "0";
        //cout<<m2<<" ";
        if (m1==0){
            temp.label = "";
            }
        if (m1==1){
            if (i<outside)
                temp.label = "A5_LAT_1,";
            if (i>=outside && i< inside)
                temp.label = "B5_LAT_1,";
            if (i>inside)
                temp.label = "C5_LAT_1";
            }
        if(m1==2){
            if (i<outside)
                temp.label = "A5_LAT_2,";
            if (i>=outside && i< inside)
                temp.label = "B5_LAT_2,";
            if (i>inside)
                temp.label = "C5_LAT_2,";
            
            }
        
        result.push_back(temp);  
    }
    return result;
}

4.13海拔高度质控

海拔高度质控主要是包括海拔极值检查等。

vector<qcFlag> altQC(vector<float> alt,float sta_alt,int length,int outside,int inside){
    qcFlag temp;
    QCMethod qcm;
    vector<int> code;
    vector<qcFlag> result;
    int m1;
    qcm.N = length;
    
    for (int i=0;i<length;i++){
        m1 = qcm.Height_extreme_check(alt[i],sta_alt);
        
        if (m1 > 1)
            temp.code = "2";
        else if (m1 == 1)
            temp.code = "1";
        else
            temp.code = "0";
        //cout<<m2<<" ";
        if (m1==0){
            temp.label = "";
            }
        if (m1==1){
            if (i<outside)
                temp.label = "A6_ALT_1,";
            if (i>=outside && i< inside)
                temp.label = "B6_ALT_1,";
            if (i>inside)
                temp.label = "C6_ALT_1";
            }
        if(m1==2){
            if (i<outside)
                temp.label = "A6_ALT_2,";
            if (i>=outside && i< inside)
                temp.label = "B6_ALT_2,";
            if (i>inside)
                temp.label = "C6_ALT_2,";
            
            }
        
        result.push_back(temp);  
    }
    return result;
}

基本要素质控方法以及参考实现就到这里了，可以留言或者私信讨论问题。