高空探测数据处理

最新推荐文章于 2024-04-23 20:56:09 发布
最新推荐文章于 2024-04-23 20:56:09 发布
阅读量691 收藏 1
点赞数
文章标签： c++ 气象产品
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/zshluckydogs/article/details/128686458
版权

高空探测数据处理主要是针对探空观测数据，按照业务要求对数据进行质控、产品生成等。主要包括原始观测数据的初级质控、中心级质控、相关产品生成等。本文涉及到的主要算法实现基于C++11及其STL模版库。

一、数据初级质控

数据初级质控主要是对温度、湿度、气压、海拔高度、经度、纬度、东向风速、北向风速、升速等要素进行极值检查，在极值范围内的保留，超出极值的剔除。

1.1气压极值检查

气压正常范围0～1100hpa，超过此范围的气压数据直接剔除

int QCMethod::pressure_extreme_check(float p)
{
    
    if ((p<0)&&(p>1100))
        {
            qcode = 2;
        }
    else
        qcode = 0;
    return qcode;
}

范围一个标志位，2表示气压极值检查错误，0表示通过，1表示可疑。

1.2温度极值检查

温度受地理位置以及气压高度影响，所以做温度极值检查时同时需要站点纬度以及温度数据的气压高度。

int QCMethod::Tem_extreme_check(float tem,float station_lat,float p)
{
    if (station_lat < 45)
    {
        if ((p>925) && (p<=1000))
        {
            if ((tem < -55) || (tem>50))
                qcode = 2;
            else if ((tem < -40) || (tem > 45))
                qcode = 1;
            else
                qcode = 0;    
        }
        if ((p<=925) && (p>850))
        {
            if ((tem < -65) || (tem>40))
                qcode = 2;
            else if ((tem < -50) || (tem > 35))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=850) && (p>700))
        {
            if ((tem < -85) || (tem>30))
                qcode = 2;
            else if ((tem < -70) || (tem > 25))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=700) && (p>500))
        {
            if ((tem < -95) || (tem>15))
                qcode = 2;
            else if ((tem < -80) || (tem > 5))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=500) && (p>400))
        {
            if ((tem < -100) || (tem>0))
                qcode = 2;
            else if ((tem < -85) || (tem > -5))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=400) && (p>300))
        {
            if ((tem < -100) || (tem>-5))
                qcode = 2;
            else if ((tem < -85) || (tem > -10))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=300) && (p>250))
        {
            if ((tem < -100) || (tem>-5))
                qcode = 2;
            else if ((tem < -85) || (tem > -10))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=250) && (p>200))
        {
            if ((tem < -100) || (tem>-5))
                qcode = 2;
            else if ((tem < -85) || (tem > -10))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=200) && (p>150))
        {
            if ((tem < -100) || (tem>-5))
                qcode = 2;
            else if ((tem < -90) || (tem > -10))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=150) && (p>100))
        {
            if ((tem < -100) || (tem>-5))
                qcode = 2;
            else if ((tem < -90) || (tem > -10))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=100) && (p>70))
        {
            if ((tem < -100) || (tem>5))
                qcode = 2;
            else if ((tem < -90) || (tem > -5))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=70) && (p>50))
        {
            if ((tem < -100) || (tem>5))
                qcode = 2;
            else if ((tem < -85) || (tem > -5))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=50) && (p>30))
        {
            if ((tem < -100) || (tem>5))
                qcode = 2;
            else if ((tem < -85) || (tem > -5))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=30) && (p>20))
        {
            if ((tem < -100) || (tem>5))
                qcode = 2;
            else if ((tem < -85) || (tem > -5))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=20) && (p>10))
        {
            if ((tem < -100) || (tem>5))
                qcode = 2;
            else if ((tem < -85) || (tem > -5))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=10) && (p>7))
        {
            if ((tem < -90) || (tem>20))
                qcode = 2;
            else if ((tem < -80) || (tem > 10))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=7) && (p>5))
        {
            if ((tem < -80) || (tem>30))
                qcode = 2;
            else if ((tem < -70) || (tem > 15))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=5) && (p>3))
        {
            if ((tem < -70) || (tem>35))
                qcode = 2;
            else if ((tem < -60) || (tem > 25))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=3) && (p>2))
        {
            if ((tem < -70) || (tem>40))
                qcode = 2;
            else if ((tem < -60) || (tem > 30))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=2) && (p>1))
        {
            if ((tem < -70) || (tem>40))
                qcode = 2;
            else if ((tem < -60) || (tem > 30))
                qcode = 1;
            else
                qcode = 0;  
        }
    
    }
    else
    {
        if ((p>925) && (p<=1000))
        {
            if ((tem < -90) || (tem>45))
                qcode = 2;
            else if ((tem < -70) || (tem > 35))
                qcode = 1;
            else
                qcode = 0;    
        }
        if ((p<=925) && (p>850))
        {
            if ((tem < -90) || (tem>35))
                qcode = 2;
            else if ((tem < -70) || (tem > 30))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=850) && (p>700))
        {
            if ((tem < -90) || (tem>25))
                qcode = 2;
            else if ((tem < -70) || (tem > 20))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=700) && (p>500))
        {
            if ((tem < -100) || (tem>10))
                qcode = 2;
            else if ((tem < -80) || (tem > 5))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=500) && (p>400))
        {
            if ((tem < -100) || (tem>0))
                qcode = 2;
            else if ((tem < -85) || (tem > -5))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=400) && (p>300))
        {
            if ((tem < -100) || (tem>-5))
                qcode = 2;
            else if ((tem < -85) || (tem > -10))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=300) && (p>250))
        {
            if ((tem < -100) || (tem>-5))
                qcode = 2;
            else if ((tem < -85) || (tem > -10))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=250) && (p>200))
        {
            if ((tem < -100) || (tem>-5))
                qcode = 2;
            else if ((tem < -85) || (tem > -10))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=200) && (p>150))
        {
            if ((tem < -100) || (tem>-5))
                qcode = 2;
            else if ((tem < -90) || (tem > -10))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=150) && (p>100))
        {
            if ((tem < -100) || (tem>-5))
                qcode = 2;
            else if ((tem < -90) || (tem > -10))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=100) && (p>70))
        {
            if ((tem < -100) || (tem>5))
                qcode = 2;
            else if ((tem < -90) || (tem > -5))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=70) && (p>50))
        {
            if ((tem < -100) || (tem>5))
                qcode = 2;
            else if ((tem < -85) || (tem > -5))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=50) && (p>30))
        {
            if ((tem < -100) || (tem>5))
                qcode = 2;
            else if ((tem < -85) || (tem > -5))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=30) && (p>20))
        {
            if ((tem < -100) || (tem>5))
                qcode = 2;
            else if ((tem < -85) || (tem > -5))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=20) && (p>10))
        {
            if ((tem < -100) || (tem>5))
                qcode = 2;
            else if ((tem < -85) || (tem > -5))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=10) && (p>7))
        {
            if ((tem < -90) || (tem>20))
                qcode = 2;
            else if ((tem < -80) || (tem > 10))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=7) && (p>5))
        {
            if ((tem < -80) || (tem>30))
                qcode = 2;
            else if ((tem < -70) || (tem > 15))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=5) && (p>3))
        {
            if ((tem < -70) || (tem>35))
                qcode = 2;
            else if ((tem < -60) || (tem > 25))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=3) && (p>2))
        {
            if ((tem < -70) || (tem>40))
                qcode = 2;
            else if ((tem < -60) || (tem > 30))
                qcode = 1;
            else
                qcode = 0;  
        }
        if ((p<=2) && (p>1))
        {
            if ((tem < -70) || (tem>40))
                qcode = 2;
            else if ((tem < -60) || (tem > 30))
                qcode = 1;
            else
                qcode = 0;  
        }
    
    }

    return qcode;
}

温度数据跟观测站点纬度以及数据的气压高度关系较为密切，无法脱离这两个要素单独判断。返回的标志位含义同上。

1.3湿度极值检查

由于实际业务中发现湿度传感器异常情况较多，所以对湿度数据不进行初级质控。

1.4海拔高度极值检查

海拔高度极值检查与观测站点的海拔有密切关系，不能低于观测站海拔，不高于大气层高度。

int QCMethod::Height_extreme_check(float hg,float station_altitude)
{
    if ((hg > station_altitude) && (hg < 40000.0))
        qcode = 0;
    else
        qcode = 2;
    return qcode;
}

1.5经纬度极值检查

观测数据的经度和纬度应在站点经度的10度范围内，超过这个距离则认为是异常数据。

int QCMethod::Lon_extreme_check(float lon,float station_lon)
{
    if (abs(lon-station_lon) < 10)
        qcode = 0;
    else
        qcode = 2;
    return qcode;
}

int QCMethod::Lat_extreme_check(float lat,float station_lat)
{
    if (abs(lat-station_lat) < 10)
        qcode = 0;
    else
        qcode = 2;
    return qcode;
}

1.6东向风速北向风速极值检查

东向风速和北向风速不应超过100m/s。

int QCMethod::esp_extreme_check(float esp){
    if (fabsf(esp)>=100)
        return 2;
    else
        return 0;
}
int QCMethod::nsp_extreme_check(float nsp){
    if (fabsf(nsp)>=100)
        return 2;
    else
        return 0;
}

1.7升速极值检查

升速不应超过1000m/s

int QCMethod::vsp_extreme_check(float vsp){
    if (fabsf(vsp)>1000)
        return 2;
    else 
        return 0;
}

解释一下，由于都是使用的float类型的数据，所以在对数据取绝对值时使用的是fabsf()方法，abs()方法是对整形数据取绝对值，这两个方法都在C++标准库中。

下一篇文章介绍比较重要，同时也是比较复杂的中心级质控。

nobrody
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
高空气象观测秒数据的处理应用研究
08-21
针对如何对L波段探空雷达系统秒级高精度观测数据加以有效利用的问题，从测站上空高空温湿特性层识别以及整层可降水量计算这两个方面，探讨了探空秒数据序列客观分析和自动处理技术，提出了能够适应天气分析和预报服务业务需要、有明确物理依据的筛选高空温湿特性层的“窗口缩放分析算法”、以及具有更高准确度的大气可降水量“秒数据直接积分算法”。结果表明：直接对探空温、压、湿（相对湿度）观测值进行秒数据积分来求取大气可降水量，与GPS反演降水量值吻合度更好、比通常仅使用规定层探空资料求取可降水量的做法准确度更高。采用“窗口缩放分析算法”能够快速、有效地定位高空气象观测秒数据序列中有明确物理意义的“温、湿显著变化点”，避免规范方法有可能出现的漏选情况，其筛选结果可提供比业务算法更丰富的温湿显著变化细节，更好地满足天气分析和预报服务业务对大气层结垂直剖面精细化分析和绘图需求；将该算法输出的特性层资料与常规等压面资料相结合进行PW积分计算，甚至可获得比探空秒数据直接积分计算PW更接近于GPS反演PW的结果。
《失控》
热门推荐
每一个不曾起舞的日子，都是对人生的辜负。
04-01 3万+
人造与天生 新生物文明 造化所生的自然王国和人类建造的人造国度正在融为一体。机器，正在生物化；而生物，正在工程化。这种趋势正验证着某些古老的隐喻——将机器比喻为生物，将生物比喻为机器。如今，那些久远的隐喻不再只是诗意的遐想，它们正在变成现实。 由于我们自己创造的这个世界变得过于复杂，我们不得不求助于自然世界以了解管理它的方法。这也就意味着，想要保证一切正常运转，我们最终制造出来的环境越机械化，可能...
大气模型软件：WRF、CMAQ、SMOKE、MCM、CAMx、Calpuff、人工智能气象、WRFchem、PMF、FLEXPART拉格朗日粒子扩散、WRF-UCM、EKMA
weixin_49857990的博客
04-20 2444
WRF模式、FLEXPART拉格朗日粒子扩散、Calpuff大气污染扩散、WRFDA资料同化、PMF源解析、MCM箱模型、EKMA曲线、WRF-Chem、CMAQ空气质量、SMOKE排放清单、WRF-Hydro、WRF-UCM、python-WRF融合技术、camx空气质量、CMIP6数据处理、人工智能气象等...
高空探测数据处理--对流层顶选取
zshluckydogs的博客
01-14 1225
在第一对流层顶以上存在一个厚度至少达1km、平均温度垂直递减率大于3摄氏度/km的气层，在该气层以上又出现温度垂直递减率小于等于2摄氏度/km的最低高度，假如此高度以上2km及其以内的任何高度与此高度间的平均温度垂直递减率也都小于等于2摄氏度/km，则此最低高度也应选为第二对流层顶。就纬度而言，对流层上界在低纬度地区平均为16－18 km，在中纬度的地区则为9－12 km，而在高纬度地区只有7－8 km。风向和风速经常变化；第一对流层顶只能有一个，如有几个气层都符合第一对流层顶条件，则选取高度最低的一个。
探空仪测得数据按照分层报表要求进行处理（程序）
算法小站
09-25 1362
clear  clc  load D:\MATLAB7\work\qiu\snd1.tpu %存放有时间、温度、湿度、气压、高度  load D:\MATLAB7\work\qiu\snd2.wnd %存放有风向、风速  temp=snd1(:,2);%温度  shi=snd1(:,3);%湿度  pres=snd1(:,4);%压强  h=snd1(:,6);%高度  fx=
气象-探空资料绘图
qq_41627642的博客
08-11 1750
''' Author: your name Date: 2021-08-10 12:00:04 LastEditTime: 2021-08-11 16:22:53 LastEditors: Please set LastEditors Description: In User Settings Edit FilePath: \Test\QHTJ\Python气象数据处理与绘图\TanKong.py ''' import matplotlib.pyplot as plt from mpl_toolkits..
高空探测数据处理--中心级质控（三）
zshluckydogs的博客
01-23 150
接上，补充剩余中心级质控方法。
高空气象探测系统数据处理与误差分析.pdf
08-15
#资源达人分享计划#
高空气象探测数据处理软件的国产化设计和实现.pdf
08-15
#资源达人分享计划#
高空气象探测微机数据处理系统软件存在的问题及注意事项.pdf
08-15
#资源达人分享计划#
读取探空数据
01-09
读取不规则排列的探空数据，写成规则排列的数据格式，或者转成其他所需的二进制格式的数据；
高空气象观测秒数据的处理应用研究-(资料数据细节）
08-21
高空气象观测秒数据 PW比较分析细节
高空探测程序“59-701”数据处理系统使用中特殊情况的处理.pdf
08-15
#资源达人分享计划#
高空气象探测“59—701”微机数据处理系统特殊情况处理的探讨.pdf
08-15
#资源达人分享计划#
【C】 初识C语言
努力学习C/C++
04-21 759
从本篇开始，带你从不会C语言到精通C语言。这篇是让大家对C语言有个初步的了解。下面的内容会拆分不同篇幅单独进行详细讲叙。数据类型变量、常量字符串转义字符注释选择语句循环语句函数数组操作符常见关键字define定义常量和宏指针结构体1.什么是C语言？C语言是一门面向过程的计算机编程语言，与C++、C#、Java等面向对象编程语言有所不同。C语言的设计目标是提供一种能以简易的方式编译、处理低级存储器、仅产生少量的机器码以及不需要任何运行环境支持便能运行的编程语言。
C++ WebServer的细节理解
我的博客
04-22 594
在非阻塞模式下，如果 I/O 操作不能立即完成（例如，因为没有数据可读 或 写缓冲区满），系统不会让调用进程阻塞等待，而是立即返回一个错误（通常为 EAGAIN 或 EWOULDBLOCK）。先查询文件描述符 fd 当前的标志，然后将 O_NONBLOCK 标志加入其中，并通过 F_SETFL 来更新文件描述符，最终实现将该文件描述符切换为非阻塞模式。这意味着，如果没有待处理的连接请求（即没有客户端尝试连接到服务器），accept() 会一直阻塞，直到有新的连接请求到达或发生其他特定条件（如超时）为止。
C++——初识模板
qq_74245477的博客
04-22 884
之前实现的数据结构，比如栈，我们为了便于修改存储的数据类型，通常使用typedef来重定义数据类型，但是这样有一个缺陷，会导致我们定义的栈存储的都是一种数据，如果同时需要两个栈，一个存储int，一个存储自定义类型，那么typedef就不管用了。比如：当用int类型使用函数模板时，编译器通过对实参类型的推演，将T确定为int类型，然后产生一份专门处理int类型的代码，对于字符类型也是如此。概念：函数模板代表了一个函数家族，该函数模板与类型无关，在使用时被参数化，根据实参类型产生函数的特定类型的版本。
C语言：贪吃蛇游戏（万字解读超详细）
2301_80655639的博客
04-20 1045
贪吃蛇这个游戏应该我们小时候都玩过吧我相信能够点进这个博客的应该都很熟悉这个游戏了如果对这个游戏不熟悉的可以网上查查怎么玩，这样才能对这个游戏有个清晰的认知该篇博客需要拥有C语言基础（对指针熟悉、会链表的使用）下面我们需要先学习一些C语言之外的东西认真坚持看完，我相信你一定能完成这个接近500行的C语言贪吃蛇项目该类型内部封装了两个成员，一个x，一个y，代表坐标我们可以将我们的控制台理解为一个平面坐标图横向正方向是x轴正方向，竖向负方向是y轴正方向。
C++中的抽象类
最新发布
hh121386的博客
04-23 182
C++中的抽象类
动力声呐探测设备数据处理方法
06-11
动力声呐探测设备采集到的数据需要进行处理和分析，以确定管道的状态和问题。一般来说，动力声呐探测设备的数据处理方法如下： 1. 数据预处理：对采集到的原始数据进行预处理，包括滤波、去噪、校准等，以保证数据...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

nobrody

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值