[发明专利]ATB格式气象数据的解析方法

权利要求书

1.ATB格式气象数据的解析方法，其特征在于，包括以下步骤：

(1)读取和解析ATB数据文件，生成浮点数矩阵；

(2)按ATB数据格式对步骤(1)的浮点数矩阵进行拆分提取，得到给定经纬度点所在的格点数据；

(3)根据步骤(2)的格点数据，生成任意经纬度点相关格点、固定气压层、特定时刻数据对象；

所述步骤(1)包括以下分步骤：

(1.1)开始读取ATB数据文件；

(1.2)计算数据文件包含的字节数；

(1.3)声明大小等于数据文件包含字节数的字节类型数组，用于存放ATB二进制读取结果；

(1.4)将ATB数据文件读入该字节类型数组；

(1.5)从字节数据的第一个字节开始读取数据；

(1.6)判断是否有剩余字节存在，如果是执行步骤(1.7)，否则执行步骤(1.11)；

(1.7)从字节类型数组中顺序读出四个字节的数据，即四个8位二进制数；

(1.8)将四个8位二进制数转化为浮点数，得到了十进制的浮点数；

(1.9)储存生成的浮点数；

(1.10)当前数组指针后移4位，执行步骤(1.6)；

(1.11)ATB二进制数据读取与解析结束；

所述步骤(2)包括以下分步骤：

(2.1)开始提取格点数据，读取准确的经度x、纬度y；

(2.2)经度数据位置判定，声明经度位置结果变量i，如果是东经经度执行步骤(2.3)，如果是西经经度执行步骤(2.4)，0度执行步骤(2.7)，180度执行步骤(2.8)；

(2.3)i＝x/1.25+24，转步骤(2.9)；

(2.4)判断经度是否小于西经30度，如果是执行步骤(2.5)，否则执行步骤(2.6)；

(2.5)i＝(x-30)/(-1.25)，转步骤(2.9)；

(2.6)i＝(x-178.5)/(-1.25)+168，转步骤(2.9)；

(2.7)i＝24，转步骤(2.9)；

(2.8)i＝168，转步骤(2.9)；

(2.9)对i进行向下取整运算，得到经度对应得标准位置；

(2.10)纬度数据位置判定，声明纬度位置结果变量j，如果是北纬纬度执行步骤(2.11)，如果是南纬纬度执行步骤(2.12)；

(2.11)j＝(y+90)/1.25，转步骤(2.13)；

(2.12)j＝(90-y)/1.25，转步骤(2.13)；

(2.13)对j进行向下取整运算，得到纬度对应得标准位置；

(2.14)确定格点数据：将(i，j)位置对应数据附值给(2，1)位置数据，(i，j+1)位置对应数据附值给(2，2)位置数据，(i+1，j)位置对应数据附值给(1，1)位置数据，(i+1，j+1)位置对应数据附值给(1，2)位置数据；

(2.15)浮点数矩阵进行拆分提取结束；

所述步骤(3)包括以下分步骤：

(3.1)开始提取数据，传入坐标点i，j的位置；

(3.2)从第一个温度和高度数据块开始处理，声明数据块标识变量k＝0；

(3.3)判断是否已经读取过最后一个数据块，如果是执行步骤(3.11)，否则执行步骤(3.4)；

(3.4)从当前数据块的第一个气压层开始处理，声明气压层标识变量l＝0；

(3.5)判断是否已经读取过最后一个气压层，如果是执行步骤(3.10)，否则执行步骤(3.6)；

(3.6)计算当前格点数据标识位置；

(3.7)按照步骤(2.14)的方法取当前格点数据；

(3.8)将当前的读取结果写入储存空间；

(3.9)指向下一个气压层，执行步骤(3.5)；

(3.10)指向下一个温度和高度数据块，执行步骤(3.3)；

(3.11)数据对象生成结束。

2.根据权利要求1所述的ATB格式气象数据的解析方法，其特征在于，所述步骤(1.8)包括以下分步骤：

(1.8.1)开始将四阶二进制数组转化为浮点数；

(1.8.2)判断当前数组是否包含全部四阶二进制数，如果是执行步骤(1.8.3)，否则执行步骤(1.8.8)；

(1.8.3)将数组[3]号元素左移24位；

(1.8.4)将数组[2]号元素左移16位；

(1.8.5)将数组[1]号元素左移8位；

(1.8.6)将数组[0]、[1]、[2]和[3]号元素合并；

(1.8.7)根据IEEE 754规则要求转化为浮点数，转步骤(1.8.9)；

(1.8.8)提示当前数组有误；

(1.8.9)浮点数转化结束。