CN115201937B - 一种下投探空阶段自动划分方法、装置、介质和设备 - Google Patents
一种下投探空阶段自动划分方法、装置、介质和设备 Download PDFInfo
- Publication number
- CN115201937B CN115201937B CN202210829076.XA CN202210829076A CN115201937B CN 115201937 B CN115201937 B CN 115201937B CN 202210829076 A CN202210829076 A CN 202210829076A CN 115201937 B CN115201937 B CN 115201937B
- Authority
- CN
- China
- Prior art keywords
- time
- sounding data
- group
- phase
- corrected
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000012163 sequencing technique Methods 0.000 claims abstract description 7
- 230000001133 acceleration Effects 0.000 claims description 42
- 230000006641 stabilisation Effects 0.000 claims description 13
- 238000011105 stabilization Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 7
- 239000000523 sample Substances 0.000 claims description 3
- 230000007812 deficiency Effects 0.000 claims description 2
- 238000005266 casting Methods 0.000 abstract description 12
- 238000012545 processing Methods 0.000 description 19
- 238000001514 detection method Methods 0.000 description 8
- 230000009467 reduction Effects 0.000 description 4
- 238000002872 Statistical quality control Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 239000010410 layer Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000005526 G1 to G0 transition Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/08—Adaptations of balloons, missiles, or aircraft for meteorological purposes; Radiosondes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Environmental & Geological Engineering (AREA)
- Computational Mathematics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Environmental Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Atmospheric Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Biology (AREA)
- Algebra (AREA)
- Aviation & Aerospace Engineering (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Ecology (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
- Electric Clocks (AREA)
Abstract
本发明涉及一种下投探空阶段自动划分方法,包括获取下投式探空仪在一次探测过程中的N组下投探空数据,对N组所述下投探空数据按照时间序列从早到晚进行排序;根据N组所述下投探空数据和气压‑高度模型,得到下投探空数据修正后的降速值和修正后的降加速值,根据所有下投探空数据的修正后的降速值和修正后的降加速值,划分下投探空的各个阶段,本申请基于气压、高度、降速、降加速等要素,提出下投探空阶段自动划分方法,可作为数据统计分析和质量控制的基础。本发明还涉及一种下投探空阶段自动划分装置、介质和设备。
Description
技术领域
本发明涉及气象探测技术领域,尤其涉及一种下投探空阶段自动划分方法、装置、介质和设备。
背景技术
下投式探空仪是指从飞机、气球等高空平台施放的,用于探测空中温度、湿度、气压、风向和风速等气象要素垂直分布的仪器。下投式探空仪可弥补常规高空探测在时间和空间上的不足,在观测交互科学试验、台风监测预报、特殊任务执行等方面发挥着重要的作用。
与常规业务用的探空仪相比,下投式探空仪的数据处理方法具有特殊性,整个探测过程可分为不同的阶段,每个阶段的探空数据具有不同的特点,为了进行后续的数据处理工作,需要人工划分各个阶段之间的界限。
发明内容
本发明所要解决的技术问题是针对现有技术的不足,提供一种下投探空阶段自动划分方法、装置、介质和设备。
本发明解决上述技术问题的技术方案如下:
一种下投探空阶段自动划分方法,所述方法包括2
S1、获取一枚下投式探空仪在一次探测过程中的N组下投探空数据,对N组所述下投探空数据按照时间序列进行排序,将所述下投探空数据中的最早时间作为所述下投探空仪在本次探测过程中飞行阶段的开始时间ts1;
S2、根据所述下投探空数据、气压-高度模型、第一阈值和第二阈值,计算每组所述下投探空数据修正后的海拔高度h′i、修正后的降速v′i和修正后的降加速a′i,i=1,2,..N;
S3、根据所有所述修正后的降加速a′i、第一预设数量n1和第三阈值,确定所述飞行阶段的结束时间te1和投放阶段的开始时间ts2;
S4、根据所有所述修正后的降速v′i、所有所述修正后的降加速a′i、第二预设数量n2、第四阈值和第五阈值,判断是否开伞成功,若是,转入步骤S5,否则,确定所述投放阶段的结束时间te2、开伞阶段的开始时间ts3、开伞阶段的结束时间te3和平稳阶段的开始时间ts4,转入步骤S7;
S5、根据所有所述修正后的降加速a′i、第三预设数量n3和第六阈值,确定所述投放阶段的结束时间te2和开伞阶段的开始时间ts3;
S6、根据所有所述修正后的降加速a′i、第四预设数量n4和第七阈值,确定所述开伞阶段的结束时间te3和平稳阶段的开始时间ts4;
S7、根据所有所述修正后的降速v′i、第五预设数量n5和第八阈值,确定所述平稳阶段的结束时间te4、着地阶段的开始时间ts5和着地阶段的结束时间te5。
本方法发明的有益效果是:提出了一种下投探空阶段自动划分方法,包括获取下投式探空仪在一次探测过程中的N组下投探空数据,对N组所述下投探空数据按照时间序列从早到晚进行排序;根据N组所述下投探空数据和气压-高度模型,得到下投探空数据的降速值和降加速值,根据所有下投探空数据的降速值和所有降加速值,划分下投探空的各个阶段,本申请基于气压、高度、降速、降加速等要素,提出下投探空阶段自动划分方法,可作为数据统计分析和质量控制的基础。
本发明解决上述技术问题的另一技术方案如下:
一种下投探空阶段自动划分装置,所述装置包括:
第一处理模块,用于获取一枚下投式探空仪在一次探测过程中的N组下投探空数据,对N组所述下投探空数据按照时间序列进行排序,将所述下投探空数据中的最早时间作为所述下投探空仪在本次探测过程中飞行阶段的开始时间ts1;
第二处理模块,用于根据所述下投探空数据、气压-高度模型、第一阈值和第二阈值,计算每组所述下投探空数据修正后的海拔高度h′i、修正后的降速v′i和修正后的降加速a′i,i=1,2,..N;
第三处理模块,用于根据所有所述修正后的降加速a′i、第一预设数量n1和第三阈值,确定所述飞行阶段的结束时间te1和投放阶段的开始时间ts2;
第四处理模块,用于根据所有所述修正后的降速v′i、所有所述修正后的降加速a′i、第二预设数量n2、第四阈值和第五阈值,判断是否开伞成功,若是,转入所述第五处理模块,否则,确定所述投放阶段的结束时间te2、开伞阶段的开始时间ts3、开伞阶段的结束时间te3和平稳阶段的开始时间ts4,转入第七处理模块
所述第五处理模块,用于根据所有所述修正后的降加速a′i、第三预设数量n3和第六阈值,确定所述投放阶段的结束时间te2和开伞阶段的开始时间ts3;
第六处理模块,用于根据所有所述修正后的降加速a′i、第四预设数量n4和第七阈值,确定所述开伞阶段的结束时间te3和平稳阶段的开始时间ts4;
所述第七处理模块,用于根据所有所述修正后的降速v′i、第五预设数量n5和第八阈值,确定所述平稳阶段的结束时间te4、着地阶段的开始时间ts5和着地阶段的结束时间te5。
本申请还提供一种计算机可读存储介质,其上存储有计算机程序指令,该程序指令被处理器执行时实现上述技术方案中任一项所述方法的步骤。
本申请提供一种电子设备,包括:存储器,其上存储有计算机程序;处理器,用于执行所述存储器中的所述计算机程序,以实现上述技术方案中任一项所述方法的步骤。
本发明附加的方面的优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明实践了解到。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对本发明实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面所描述的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例所示的一种下投探空阶段自动划分方法的流程示意图;
图2为本发明另一实施例所示的一种下投探空阶段自动划分装置的模块示意图;
图3为本发明另一实施例所示的一种下投探空阶段自动划分方法的流程示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明的一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都应属于本发明保护的范围。
如图1所示,本发明实施例所述的一种下投探空阶段自动划分方法,方法包括:
110、获取一枚下投式探空仪在一次探测过程中的N组下投探空数据,对N组下投探空数据按照时间序列进行排序,将下投探空数据中的最早时间作为下投探空仪在本次探测过程中飞行阶段的开始时间ts1。
120、根据下投探空数据、气压-高度模型、第一阈值和第二阈值,计算每组下投探空数据修正后的海拔高度h′i、修正后的降速v′i和修正后的降加速a′i,i=1,2,..N。
130、根据所有修正后的降加速a′i、第一预设数量n1和第三阈值,确定飞行阶段的结束时间te1和投放阶段的开始时间ts2。
140、根据所有修正后的降速v′i、所有修正后的降加速a′i、第二预设数量n2、第四阈值和第五阈值,判断是否开伞成功,若是,转入步骤150,否则,确定投放阶段的结束时间te2、开伞阶段的开始时间ts3、开伞阶段的结束时间te3和平稳阶段的开始时间ts4,转入步骤170。
150、根据所有修正后的降加速a′i、第三预设数量n3和第六阈值,确定投放阶段的结束时间te2和开伞阶段的开始时间ts3。
160、根据所有修正后的降加速a′i、第四预设数量n4和第七阈值,确定开伞阶段的结束时间te3和平稳阶段的开始时间ts4。
170、根据所有修正后的降速v′i、第五预设数量n5和第八阈值,确定平稳阶段的结束时间te4、着地阶段的开始时间ts5和着地阶段的结束时间te5。
进一步地,步骤120具体包括:
121、计算气压海拔高度hp:
对从第2组至第N-1组气压海拔高度hpi,分别使用第i-1组下投探空数据和第i+1组下投探空数据,根据气压-高度模型分别计算得到的高度值后,再计算所得到的两个高度值的平均值;
第1组下投探空数据的气压海拔高度hp1使用第2组下投探空数据输入下式计算;
第N组下投探空数据的气压海拔高度hpN使用第N-1组下投探空数据输入下式计算;
计算公式为:
当i=2,3,...N-1时
其中,Rd是干空气的比气体常数,g是重力加速度,Tv是虚温,hi是第i组下投探空数据的海拔高度,p′是各高度层的气压值;
122、计算修正后的海拔高度h′:
将原始海拔高度h与气压海拔高度hpi进行对比,对于第i组下投探空数据,若hi缺测或hi与hpi之差的绝对值大于所述第一阈值,则第i组修正后的海拔高度hi′=hpi,否则h′i=hi;
123、计算修正后的降速v′:
将第2组至第N组所述修正后的海拔高度h′i和时间ti输入下式进行计算,得到第2组至第N组气压降速vpi,
计算公式为:
当i=2,3,...N时
将原始降速v与气压降速vpi进行对比,对于从第2组至第N组下投探空数据,若第i组下投探空数据中的降速vi缺测或vi与vpi之差的绝对值大于所述第二阈值,则从第2组至第N组修正后的降速v′i=vpi,否则v′i=vi;
124、计算修正后的降加速a′:
将第2组至第N组修正后的降速和时间ti输入公式
当i=2,3,...N时,得到第2组至第N组修正后的降加速a′i。
进一步地,步骤130中根据所有修正后的降加速a′i、第一预设数量n1和第三阈值,确定飞行阶段的结束时间te1和投放阶段的开始时间ts2,具体包括:
按照时间序列遍历下投探空数据。
当首次出现连续第一预设数量n1组修正后的降加速a′i的值均大于第三阈值时,则投放阶段的开始时间ts2是连续n1组下投探空数据中的第一组下投探空数据的时间。
飞行阶段的结束时间te1是连续n1组下投探空数据的前一组下投探空数据的时间。其中,n1是正整数。
进一步地,步骤140中根据所有修正后的降速v′i、所有修正后的降加速a′i、第二预设数量n2、第四阈值和第五阈值,判断是否开伞成功,具体包括:
按照时间序列从投放阶段的开始时间ts2对应的下投探空数据的组号开始遍历。
当首次出现连续第二预设数量n2组修正后的降加速a′i的绝对值均小于第四阈值,且连续第二预设数量n2组下投探空数据的修正后的降速v′i的平均值小于第五阈值时,则开伞成功;
否则,开伞失败。其中,n2是正整数。
进一步地,步骤140中确定投放阶段的结束时间te2、开伞阶段的开始时间ts3、开伞阶段的结束时间te3和平稳阶段的开始时间ts4,具体包括:
平稳阶段的开始时间ts4是连续n2组下投探空数据中的第一组下投探空数据的时间。
投放阶段的结束时间te2是连续n2组下投探空数据的前一组下投探空数据的时间。
开伞阶段的开始时间ts3和结束时间te3均为空值。
进一步地,步骤150中根据所有修正后的降加速a′i、第三预设数量n3和第六阈值,确定投放阶段的结束时间te2和开伞阶段的开始时间ts3,具体包括:
按照时间序列从投放阶段的开始时间te2对应的下投探空数据的组号开始遍历。
当首次出现连续第三预设数量n3组修正后的降加速a′i的值均小于第六阈值时,则开伞阶段的开始时间ts3是连续n3组下投探空数据中的第一组下投探空数据的时间。
投放阶段的结束时间te2是连续n3组下投探空数据的前一组下投探空数据的时间。其中,n3是正整数。
进一步地,步骤160中根据所有修正后的降加速a′i、第四预设数量n4和第七阈值,确定开伞阶段的结束时间te3和平稳阶段的开始时间ts4,具体包括:
按照时间序列从投放阶段的开始时间ts3对应的下投探空数据的组号开始遍历。
当首次出现连续第四预设数量n4组修正后的降加速a′i的绝对值均小于第七阈值时,则平稳阶段的开始时间ts4是连续n4组下投探空数据中的第一组下投探空数据的时间。
开伞阶段的结束时间te3是连续n4组下投探空数据的前一组下投探空数据的时间。其中,n4是正整数。
进一步地,步骤170中根据所有修正后的降速v′i、第五预设数量n5和第八阈值,确定平稳阶段的结束时间te4、着地阶段的开始时间ts5和着地阶段的结束时间te5,具体包括:
按照时间序列从平稳阶段的开始时间ts4对应的下投探空数据的组号开始遍历。
当首次出现连续第五预设数量n5组修正后的降速v′i的值均小于第八阈值时,则着地阶段的开始时间ts5是连续n5组下投探空数据中的第一组下投探空数据的时间。
平稳阶段的结束时间te4是连续n5组下投探空数据的前一组下投探空数据的时间。其中,n5是整数。
着地阶段的结束时间te5是第N组下投探空数据的时间。
应理解,本申请基于气压、海拔高度、降速、降加速等要素,提出下投探空阶段自动划分方法,可作为数据统计分析和质量控制的基础。
一、各阶段定义
根据降速和姿态的不同,将下投式探空仪的探测过程分为5个阶段:
(一)飞行阶段
探空仪装载在吊舱内,温度、湿度、气压传感器测量的是吊舱内的环境状态,风向、风速数据表征飞行平台的运动状态。记飞行阶段的开始和结束时间分别为ts1和te1。
(二)投放阶段
探空仪从吊舱弹出,连续翻滚,姿态不稳定,降速不断增大。记投放阶段的开始和结束时间分别为ts2和te2。
(三)开伞阶段
降落伞打开,探空仪姿态逐渐稳定,降速继续增加至最大值,然后开始减小。记开伞阶段的开始和结束时间分别为ts3和te3。
(四)平稳阶段
此阶段的时间最长,探空仪的降速趋于稳定,风向、风速仅在此阶段开伞成功的情况下表征大气运动状态。记平稳阶段的开始和结束时间分别为ts4和te4。
(五)着地阶段
探空仪降落到地面,降速迅速变为零。传感器可能由于损坏而导致数据无效,也可能因地形遮蔽而导致无数据。记着地阶段的开始和结束时间分别为ts5和te5。
二、阶段划分流程
假设一枚下投式探空仪共探测到N组数据,每组数据包括时间ti、海拔高度hi、气压pi、温度Ti、相对湿度Ui、风向θi、风速ri、降速vi,i=1,2,..N。本申请提出的下投探空阶段划分流程如所图3所示,具体步骤如下:
(一)初始化
读取原始下投探空数据,将ts1赋值为t1。
(二)静力学检查
首先计算气压海拔高度hp。对第i组海拔高度hpi(i=2,3,..N-1),分别使用第i-1组和第i+1组探空数据,根据气压-高度模型计算高度,再求出二者的平均值;对于第1组和第N组,分别使用第2组和第N-1组探空数据单独计算。计算方法为:
当i=2,3,...N-1时
其中,Rd是干空气的比气体常数,在近地面层的值为287.05J/(kg·K);g是重力加速度,与纬度和高度有关;Tv是虚温,可由各层的温度、湿度、气压探测数据计算得出。
其次,计算修正后的海拔高度h′。将原始海拔高度h与气压海拔高度hp进行对比,对于第i组数据(i=1,2,..N),若hi缺测或hi与hpi之差的绝对值大于阈值一,则h′i=hpi,否则h′i=hi。
再次,计算修正后的降速v′。根据修正后的海拔高度h′和时间,计算气压降速vp。计算方法为:
将原始降速v与气压降速vp进行对比,对于第i组数据(i=2,3,..N),若vi缺测或vi与vpi之差的绝对值大于阈值二,则v′i=vpi,否则v′i=vi。
最后,计算修正后的降加速a′。计算方法为:
(三)划分飞行和投放阶段的界限
从第2组探空数据开始向后遍历,当首次出现连续n1个降加速值大于阈值三的情况时,将ts2赋值为连续n1组探空数据的开始时间,te1赋值为ts2的前一组探空数据时间。
(四)判断开伞是否成功
从ts2开始向后遍历,当首次出现连续n2个降加速的绝对值小于阈值四的情况时,若连续n2组探空数据的降速平均值小于阈值五,则开伞成功,转入步骤(五);否则开伞失败,将ts4赋值为连续n2组探空数据的开始时间,te2赋值为ts4的前一组探空数据时间,ts3、te3均赋值为空值,转入步骤(七)。
(五)划分投放和开伞阶段的界限
从ts2开始向后遍历,当首次出现连续n3个降加速值小于阈值六的情况时,将ts3赋值为连续n3组探空数据的开始时间,te2赋值为ts3的前一组探空数据时间。
(六)划分开伞和平稳阶段的界限
从ts3开始向后遍历,当首次出现连续n4个降加速的绝对值小于阈值七的情况时,将ts4赋值为连续n4组的前一组探空数据时间,te3赋值为ts4的前一组探空数据时间。
(七)划分平稳和着地阶段的界限
从ts4开始向后遍历,若出现连续n5个降速值小于阈值八的情况时,则将te4赋值为连续n5组的前一组探空数据时间,ts5赋值为首个连续n5组探空数据的开始时间,te5赋值为tN;若未出现,则将te4赋值为tN,ts5、te5均赋值为空值。
三、结果个例
选取2枚某型探空仪的试验数据,对本专利的方法加以说明。阈值一取值为10m,阈值二取值为5m/s,阈值三取值为2m/s2,阈值四取值为3m/s2,阈值五取值为15m/s,阈值六取值为2m/s2,阈值七取值为1m/s2,阈值八取值为3m/s,n1取值为3,n2取值为10,n3取值为3,n4取值为5,n5取值为5。
表1探空仪1部分数据
探空仪1为开伞成功的个例,表1显示了部分数据。根据前述方法,划分出飞行阶段为14:08:16至14:09:06,投放阶段14:09:07为至14:09:13,开伞阶段为14:09:14至14:09:20,平稳阶段为14:09:21至14:13:26,着地阶段为14:13:27至14:13:31。
表2探空仪2部分数据
探空仪2为开伞失败的个例,表2显示了部分数据。根据前述方法,划分出飞行阶段为14:16:22至14:17:24,投放阶段14:17:25为至14:17:30,平稳阶段为14:17:31至14:19:22,无开伞阶段和着地阶段。
如图2所示,一种下投探空阶段自动划分装置,装置包括:
第一处理模块,用于获取一枚下投式探空仪在一次探测过程中的N组下投探空数据,对N组下投探空数据按照时间序列进行排序,将下投探空数据中的最早时间作为下投探空仪在本次探测过程中飞行阶段的开始时间ts1;
第二处理模块,用于根据下投探空数据、气压-高度模型、第一阈值和第二阈值,计算每组下投探空数据修正后的海拔高度h′i、修正后的降速v′i和修正后的降加速a′i,i=1,2,..N;
第三处理模块,用于根据所有修正后的降加速a′i、第一预设数量n1和第三阈值,确定飞行阶段的结束时间te1和投放阶段的开始时间ts2;
第四处理模块,用于根据所有修正后的降速v′i、所有修正后的降加速a′i、第二预设数量n2、第四阈值和第五阈值,判断是否开伞成功,若是,转入所述第五处理模块,否则,确定投放阶段的结束时间te2、开伞阶段的开始时间ts3、开伞阶段的结束时间te3和平稳阶段的开始时间ts4,转入第七处理模块
所述第五处理模块,用于根据所有修正后的降加速a′i、第三预设数量n3和第六阈值,确定投放阶段的结束时间te2和开伞阶段的开始时间ts3;
第六处理模块,用于根据所有修正后的降加速a′i、第四预设数量n4和第七阈值,确定开伞阶段的结束时间te3和平稳阶段的开始时间ts4;
所述第七处理模块,用于根据所有修正后的降速v′i、第五预设数量n5和第八阈值,确定平稳阶段的结束时间te4、着地阶段的开始时间ts5和着地阶段的结束时间te5。
本申请还提供一种计算机可读存储介质,其上存储有计算机程序指令,该程序指令被处理器执行时实现上述技术方案中任一项所述方法的步骤。
本申请提供一种电子设备,包括:存储器,其上存储有计算机程序;处理器,用于执行所述存储器中的所述计算机程序,以实现上述技术方案中任一项所述方法的步骤。
以上,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到各种等效的修改或替换,这些修改或替换都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。
Claims (8)
1.一种下投探空阶段自动划分方法,其特征在于,所述方法包括:
S1、获取一枚下投式探空仪在一次探测过程中的N组下投探空数据,对N组所述下投探空数据按照时间序列进行排序,将所述下投探空数据中的最早时间作为所述下投探空仪在本次探测过程中飞行阶段的开始时间ts1;
S2、根据所述下投探空数据、气压-高度模型、第一阈值和第二阈值,计算每组所述下投探空数据修正后的海拔高度hi′、修正后的降速vi′和修正后的降加速ai′,i=1,2,..N;
S3、根据所有所述修正后的降加速ai′、第一预设数量n1和第三阈值,确定所述飞行阶段的结束时间te1和投放阶段的开始时间ts2;
S4、根据所有所述修正后的降速vi′、所有所述修正后的降加速ai′、第二预设数量n2、第四阈值和第五阈值,判断是否开伞成功,若是,转入步骤S5,否则,确定所述投放阶段的结束时间te2、开伞阶段的开始时间ts3、开伞阶段的结束时间te3和平稳阶段的开始时间ts4,转入步骤S7;
S5、根据所有所述修正后的降加速ai′、第三预设数量n3和第六阈值,确定所述投放阶段的结束时间te2和开伞阶段的开始时间ts3;
S6、根据所有所述修正后的降加速ai′、第四预设数量n4和第七阈值,确定所述开伞阶段的结束时间te3和平稳阶段的开始时间ts4;
S7、根据所有所述修正后的降速vi′、第五预设数量n5和第八阈值,确定所述平稳阶段的结束时间te4、着地阶段的开始时间ts5和着地阶段的结束时间te5;
所述S2中具体包括:
S21、计算气压海拔高度hp:
对从第2组至第N-1组气压海拔高度hpi,分别使用第i-1组所述下投探空数据和第i+1组所述下投探空数据,根据气压-高度模型分别计算得到的高度值后,再计算所得到的两个高度值的平均值;
第1组所述下投探空数据的气压海拔高度hp1使用第2组所述下投探空数据输入下式计算;
第N组所述下投探空数据的气压海拔高度hpN使用第N-1组所述下投探空数据输入下式计算;
计算公式为:
其中,Rd是干空气的比气体常数,g是重力加速度,Tv是虚温,hi是第i组所述下投探空数据的海拔高度,p′是各高度层的气压值;
S22、计算修正后的海拔高度h′:
将原始海拔高度h与所述气压海拔高度hpi进行对比,对于第i组所述下投探空数据,若hi缺测或hi与hpi之差的绝对值大于所述第一阈值,则第i组修正后的海拔高度hi′=hpi,否则hi′=hi;
S23、计算修正后的降速v′:
将第2组至第N组所述修正后的海拔高度hi′和时间ti输入下式进行计算,得到第2组至第N组所述气压降速vpi,
计算公式为:
当i=2,3,...N时
将原始降速v与所述气压降速vpi进行对比,对于从第2组至第N组所述下投探空数据,若第i组所述下投探空数据中的降速vi缺测或vi与所述vpi之差的绝对值大于所述第二阈值,则从第2组至第N组所述修正后的降速vi′=vpi,否则vi′=vi;
S24、计算修正后的降加速a′:
将第2组至第N组所述修正后的降速和时间ti输入公式
当i=2,3,...N时,得到第2组至第N组所述修正后的降加速ai′。
2.根据权利要求1所述的下投探空阶段自动划分方法,其特征在于,所述S3具体包括:
按照时间序列遍历所述下投探空数据;
当首次出现连续所述第一预设数量n1组所述修正后的降加速ai′的值均大于所述第三阈值时,则所述投放阶段的开始时间ts2是连续n1组所述下投探空数据中的第一组所述下投探空数据的时间;
所述飞行阶段的结束时间te1是连续n1组所述下投探空数据的前一组所述下投探空数据的时间,n1是正整数。
3.根据权利要求1所述的下投探空阶段自动划分方法,其特征在于,所述步骤S4中根据所有修正后的降速vi′、所有修正后的降加速ai′、第二预设数量n2、第四阈值和第五阈值,判断是否开伞成功,具体包括:
按照时间序列从所述投放阶段的开始时间ts2对应的所述下投探空数据的组号开始遍历;
当首次出现连续所述第二预设数量n2组所述修正后的降加速ai′的绝对值均小于所述第四阈值,且连续所述第二预设数量n2组所述下投探空数据的修正后的降速vi′的平均值小于所述第五阈值时,则开伞成功;
否则,开伞失败,n2是正整数。
4.根据权利要求1所述的下投探空阶段自动划分方法,其特征在于,所述S4中确定投放阶段的结束时间te2、开伞阶段的开始时间ts3、开伞阶段的结束时间te3和平稳阶段的开始时间ts4,具体包括:
所述平稳阶段的开始时间ts4是连续n2组所述下投探空数据中的第一组所述下投探空数据的时间;
所述投放阶段的结束时间te2是连续n2组所述下投探空数据的前一组所述下投探空数据的时间;
所述开伞阶段的开始时间ts3和结束时间te3均为空值。
5.根据权利要求1所述的下投探空阶段自动划分方法,其特征在于,所述S5具体包括:
按照时间序列从投放阶段的开始时间te2对应的下投探空数据的组号开始遍历;
当首次出现连续第三预设数量n3组修正后的降加速ai′的值均小于第六阈值时,则开伞阶段的开始时间ts3是连续n3组下投探空数据中的第一组下投探空数据的时间;
投放阶段的结束时间te2是连续n3组下投探空数据的前一组所述下投探空数据的时间,其中,n3是正整数;
所述S6具体包括:
按照时间序列从投放阶段的开始时间ts3对应的下投探空数据的组号开始遍历;
当首次出现连续第四预设数量n4组修正后的降加速ai′的绝对值均小于第七阈值时,则平稳阶段的开始时间ts4是连续n4组下投探空数据中的第一组下投探空数据的时间;
开伞阶段的结束时间te3是连续n4组下投探空数据的前一组下投探空数据的时间,其中,n4是正整数。
6.根据权利要求1所述的下投探空阶段自动划分方法,其特征在于,所述S7中根据所有修正后的降速vi′、第五预设数量n5和第八阈值,确定平稳阶段的结束时间te4、着地阶段的开始时间ts5和着地阶段的结束时间te5,具体包括:
按照时间序列从所述平稳阶段的开始时间ts4对应的所述下投探空数据的组号开始遍历;
当首次出现连续所述第五预设数量n5组所述修正后的降速vi′的值均小于所述第八阈值时,则所述着地阶段的开始时间ts5是连续n5组所述下投探空数据中的第一组所述下投探空数据的时间;
所述平稳阶段的结束时间te4是连续n5组所述下投探空数据的前一组所述下投探空数据的时间,其中,n5是整数;
所述着地阶段的结束时间te5是第N组所述下投探空数据的时间。
7.一种计算机可读存储介质,其上存储有计算机程序指令,其特征在于,该程序指令被处理器执行时实现权利要求1至6中任一项所述方法的步骤。
8.一种电子设备,其特征在于,包括:存储器,其上存储有计算机程序;处理器,用于执行所述存储器中的所述计算机程序,以实现权利要求1至6中任一项所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210829076.XA CN115201937B (zh) | 2022-07-15 | 2022-07-15 | 一种下投探空阶段自动划分方法、装置、介质和设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210829076.XA CN115201937B (zh) | 2022-07-15 | 2022-07-15 | 一种下投探空阶段自动划分方法、装置、介质和设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115201937A CN115201937A (zh) | 2022-10-18 |
CN115201937B true CN115201937B (zh) | 2024-01-26 |
Family
ID=83582078
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210829076.XA Active CN115201937B (zh) | 2022-07-15 | 2022-07-15 | 一种下投探空阶段自动划分方法、装置、介质和设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115201937B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB191505750A (en) * | 1915-04-16 | 1919-03-06 | Leonard Eugene Cowey | Improvements in or relating to Sighting Devices for use in Dropping Bombs and other Missiles from Aircraft. |
CN102183962A (zh) * | 2011-03-11 | 2011-09-14 | 中国气象局气象探测中心 | 无人机下投探空仪及其释放控制方法 |
CN104417759A (zh) * | 2013-08-28 | 2015-03-18 | 无锡慧思顿科技有限公司 | 一种基于多参数多传感的智能降落伞 |
CN110576972A (zh) * | 2019-08-30 | 2019-12-17 | 中国气象局气象探测中心 | 一种多筒式下投探空仪弹射装置、投放装置及投放方法 |
CN209946419U (zh) * | 2019-06-26 | 2020-01-14 | 中国科学院大气物理研究所 | 基于高空气球飞行平台的全天候下投气象探空设备 |
CN110927826A (zh) * | 2019-12-06 | 2020-03-27 | 航天新气象科技有限公司 | 下投式探空仪及其投放方法和系统 |
WO2021169318A1 (zh) * | 2020-02-25 | 2021-09-02 | 东南大学 | 一种基于抛物线的区域对流层湿延迟计算方法 |
CN214375343U (zh) * | 2021-03-22 | 2021-10-08 | 中国气象局气象探测中心 | 一种浮标型下投式探空装置 |
-
2022
- 2022-07-15 CN CN202210829076.XA patent/CN115201937B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB191505750A (en) * | 1915-04-16 | 1919-03-06 | Leonard Eugene Cowey | Improvements in or relating to Sighting Devices for use in Dropping Bombs and other Missiles from Aircraft. |
CN102183962A (zh) * | 2011-03-11 | 2011-09-14 | 中国气象局气象探测中心 | 无人机下投探空仪及其释放控制方法 |
CN104417759A (zh) * | 2013-08-28 | 2015-03-18 | 无锡慧思顿科技有限公司 | 一种基于多参数多传感的智能降落伞 |
CN209946419U (zh) * | 2019-06-26 | 2020-01-14 | 中国科学院大气物理研究所 | 基于高空气球飞行平台的全天候下投气象探空设备 |
CN110576972A (zh) * | 2019-08-30 | 2019-12-17 | 中国气象局气象探测中心 | 一种多筒式下投探空仪弹射装置、投放装置及投放方法 |
CN110927826A (zh) * | 2019-12-06 | 2020-03-27 | 航天新气象科技有限公司 | 下投式探空仪及其投放方法和系统 |
WO2021169318A1 (zh) * | 2020-02-25 | 2021-09-02 | 东南大学 | 一种基于抛物线的区域对流层湿延迟计算方法 |
CN214375343U (zh) * | 2021-03-22 | 2021-10-08 | 中国气象局气象探测中心 | 一种浮标型下投式探空装置 |
Non-Patent Citations (4)
Title |
---|
火箭弹下投探测台风气象参数新技术及初步试验;雷小途;雷明;赵兵科;王江;汤杰;高辉文;曾智华;李国兰;罗建平;吴敏;余晖;叶雷;方平治;杨道安;张帅;邹勇波;周鄂炯;林立旻;陈平;钟业诚;石保平;何海龙;李运生;杨晓晖;赵涛;鲍旭炜;刘卓;陈国民;邬华明;吴丹;黄飞;王建军;罗珊;伊蕾;;科学通报(第32期);全文 * |
热带气旋科学观测试验及研究进展概况;廖菲;邓华;赵中阔;刘春霞;梁建茵;;热带气象学报(04);全文 * |
钱山 ; 孙守明 ; 李恒年 ; 王海涛 ; 马宏 ; .基于气象风修正的返回舱落点快速预报方法.飞行器测控学报.2013,(03),全文. * |
飞机下投探空在台风探测中的应用;李杨;马舒庆;官福顺;;气象科技(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN115201937A (zh) | 2022-10-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109800749A (zh) | 一种文字识别方法及装置 | |
CN115201937B (zh) | 一种下投探空阶段自动划分方法、装置、介质和设备 | |
CN113607205B (zh) | 一种航空发动机传感器故障检测方法及装置 | |
CN110816874B (zh) | 一种通过地面滑行试验辨识双发飞机平衡极曲线的方法 | |
WO2022145759A1 (ko) | 비행시험 데이터를 이용한 6 자유도 비행 시뮬레이션 정확도 향상 방법 및 장치 | |
CN112836581A (zh) | 一种基于相关性分析的敏感故障特征提取方法及装置 | |
CN117037006B (zh) | 一种高续航能力的无人机跟踪方法 | |
CN115828416A (zh) | 一种两级vtvl运载火箭点对点运输全程弹道设计方法 | |
CN113906353A (zh) | 用于对风力涡轮机的阻力系数进行计算机实现确定的方法 | |
CN116702439B (zh) | 基于再入落点预测模型的仿真用飞行器再入落点预测方法 | |
CN114492176A (zh) | 一种基于生成对抗网络的动力学模型参数辨识方法及系统 | |
CN112763138A (zh) | 一种飞越塔静压校准方法 | |
CN110826189B (zh) | 一种飞行器缩比模型实验系统的确定方法 | |
CN113837254B (zh) | 一种计算机视觉中图像目标的检测方法 | |
CN117113870A (zh) | 一种考虑风干扰的固定翼飞行器气动力建模方法 | |
CN110411707B (zh) | 串联飞行器级间分离气动特性干扰量预测方法 | |
CN114324975A (zh) | 一种风向确定方法、装置、设备和介质 | |
CN113221719B (zh) | 飞行器故障诊断方法、装置和电子设备 | |
Huang et al. | Research on recognition and location method of insulator in infrared image based on deep learning | |
CN117390899B (zh) | 一种运输类飞机副翼最大铰链力矩确定方法 | |
Bhat | Near surface atmospheric characteristics over the North Bay of Bengal during the Indian summer monsoon | |
CN113532445A (zh) | 一种卷帘曝光星敏感器的高动态快速自主捕获方法 | |
CN114757332A (zh) | 一种航空器气动传感器故障智能检测方法 | |
CN114781281B (zh) | 机弹分离投放试验垂直加速度模拟不足修正方法及装置 | |
CN112182857A (zh) | 火箭子级残骸落点预示方法、设备及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |