发明名称 高空格点气象数据的四维插值方法
摘要 公开了一种高空格点气象数据的四维插值方法,包括以下步骤:(1)对等气压层间实际放行高度进行插值计算;(2)进行指定经纬度格点间的双线性插值计算;(3)进行指定时间点的线性插值计算。本方法通过将高空格点气象数据的四维插值分析和计算,能够得到标准大气条件下规定的高空任意飞行高度层、任意经纬度和任意时间点在实际大气环境中的气压高度数据,为后续的数据处理提供了有效的数据源。
申请公布号 CN101655568B 申请公布日期 2011.05.04
申请号 CN200910088653.9 申请日期 2009.07.06
申请人 民航数据通信有限责任公司;北京航空航天大学 发明人 朱衍波;金开研;张军;唐金翔;许有臣;兆珺;吕嘉川
分类号 G01W1/00(2006.01)I 主分类号 G01W1/00(2006.01)I
代理机构 北京中北知识产权代理有限公司 11253 代理人 冯梦洪
主权项 1.高空格点气象数据的四维插值方法,其特征在于,包括以下步骤:(1)对等气压层间实际放行高度进行插值计算;(2)进行指定经纬度格点间的双线性插值计算;(3)进行指定时间点的线性插值计算;其中所述步骤(1)包括以下分步骤:(1.1)开始插值计算分配高度层的真实高度;(1.2)根据FL的数值取定各已知参数数值,FL是标准百英尺高度层;(1.3)计算第二辅助参数K<sub>2</sub>,其中<img file="DEST_PATH_FSB00000439240600011.GIF" wi="300" he="115" />t<sub>1</sub>和t<sub>2</sub>是许可飞行高度层所对应上下mb气压层的真实温度,H<sub>1</sub>和H<sub>2</sub>是标准大气环境下许可飞行高度层所对应上下mb气压层的标准高度;(1.4)判断当前高度层是否高于FL370,如果是执行步骤(1.6),否则执行步骤(1.5);(1.5)按较低空域计算公式,计算<img file="DEST_PATH_FSB00000439240600012.GIF" wi="483" he="107" />较低空域包含FL280至FL370共10个高度层,执行步骤(1.7);(1.6)按较高空域计算公式,计算<img file="DEST_PATH_FSB00000439240600013.GIF" wi="490" he="105" />较高空域包含FL380至FL450共8个高度层,执行步骤(1.7);(1.7)计算第一辅助参数K<sub>1</sub>,其中<img file="DEST_PATH_FSB00000439240600014.GIF" wi="564" he="150" />P<sub>1</sub>和P<sub>2</sub>是许可飞行高度层所对应上下mb气压层的mb值,p<sub>fl</sub>是标准大气环境下高度层 对应的标准大气压强,<img file="FSB00000372134200021.GIF" wi="30" he="41" />是许可飞行高度层所对应上下mb气压层的真实平均温度,<img file="FSB00000372134200022.GIF" wi="219" he="99" /><img file="FSB00000372134200023.GIF" wi="37" he="58" />是由相关上下气压层和百英尺高度层确定的平均环境温度;(1.8)判断当前高度层是否高于FL370,如果是执行步骤(1.10),否则执行步骤(1.9);(1.9)按较低空域计算公式,计算H<sub>fl</sub>=h<sub>1</sub>+K<sub>1</sub>·(h<sub>2</sub>-h<sub>1</sub>),h<sub>1</sub>和h<sub>2</sub>是许可飞行高度层所对应上下mb气压层的真实高度,H<sub>fl</sub>是差值公式最终结果,实际大气环境下的放行高度,较低空域包含FL280至FL370共10个高度层,执行步骤(1.11);(1.10)按较高空域计算公式,计算H<sub>fl</sub>=h<sub>2</sub>+K<sub>1</sub>·(h<sub>2</sub>-h<sub>1</sub>),较高空域包含FL380至FL450共8个高度层,执行步骤(1.11);(1.11)返回计算结果H<sub>fl</sub>值,为高空分配高度层的真实高度数据;(1.12)等气压层间实际放行高度插值计算结束;所述步骤(2)包括以下分步骤:(2.1)开始双线性插值计算;(2.2)给各项变量附初值,令P(1,1)点的经纬度分别为x<sub>1</sub>、y<sub>1</sub>,数据值为f(1,1),P(1,2)点的经纬度分别为x<sub>1</sub>、y<sub>2</sub>,数据值为f(1,2),P(2,1)点的经纬度分别为x<sub>2</sub>、y<sub>1</sub>,数据值为f(2,1),P(2,2)点的经纬度分别为x<sub>2</sub>、y<sub>2</sub>,数据值为f(2,2),M(x,y)点的经纬度为x、y,数据值为f(x,y);(2.3)计算Q<sub>1</sub>(x,y<sub>1</sub>)点的数据值,即<img file="FSB00000372134200024.GIF" wi="924" he="111" />(2.4)计算Q<sub>2</sub>(x,y<sub>2</sub>)点的数据值,即<img file="FSB00000372134200025.GIF" wi="941" he="111" />(2.5)计算M(x,y)点的数据值,即<img file="FSB00000372134200031.GIF" wi="873" he="110" />返回f(x,y)数据值;(2.6)双线性插值过程结束;所述步骤(3)包括以下分步骤:(3.1)开始线性插值计算;(3.2)给各项变量附初值,令T时刻点的时间为t<sub>1</sub>,数据值为f(1),T+6时刻点的时间为t<sub>2</sub>,数据值为f(2),待求解时刻的时间为t,数据值为f(t);(3.3)计算t时刻的数据值,即<img file="FSB00000372134200032.GIF" wi="638" he="129" />返回f(t)数据值;(3.4)线性插值过程结束。 
地址 100083 北京市海淀区北四环中路238号柏彦大厦16层