WAFS BUFR 航 空 数 据 的 解 码 及 可 视 化 研 究 欧 壮 杰 ( 民 航 中 南 空 管 局 通 信 网 络 管 理 中 心, 广 州,510405) 摘 要 本 文 简 要 地 介 绍 了 BUFR 数 据 的 格 式 特 点, 对 WAFS BUFR 航 空 数 据 的 引 入 原 由 做 了 分 析 ; 对 BUFR 文 件 格 式 描 述 符 及 码 表 进 行 了 详 细 的 解 剖, 给 出 了 解 码 公 式 和 解 码 描 述 符 例 子 本 文 还 通 过 解 码 后 的 数 据, 进 行 可 视 化 的 探 讨 研 究, 给 出 了 BUFR 数 据 可 视 化 的 内 容, 包 括 底 图 投 影 曲 线 生 成 图 形 缩 放 编 辑 的 方 法 公 式 和 实 现 程 序 本 文 对 BUFR 数 据 的 解 码 及 可 视 化 工 作 有 着 十 分 重 要 的 参 考 价 值 关 键 词 :BUFR, 解 码, 描 述 符, 码 表, 可 视 化, 地 球 投 影, 平 滑 曲 线 1 BUFR 数 据 格 式 简 介 BUFR 是 Binary Universal Form for the Representation of meteorological data 的 缩 写, 即 描 述 气 象 数 据 的 二 进 制 统 一 格 式 BUFR 最 早 于 1988 年 开 始 应 用 于 描 述 数 据, 那 时,BUFR 数 据 格 式 应 用 于 卫 星 航 空 器 风 洞 和 热 带 气 旋 的 观 测 2002 年,WMO( 世 界 气 象 组 织 ) 发 布 了 FM94 BUFR 数 据 格 式 指 导 文 档, 对 BUFR 数 据 的 格 式 编 码 / 解 码 方 法 做 了 详 细 的 描 述 此 后, 大 量 的 气 象 数 据 采 用 BUFR 格 式 来 描 述 存 储 和 交 换 BUFR 是 WMO 用 于 描 述 和 交 换 观 测 数 据 的 唯 一 推 荐 编 码 格 式, 被 建 议 在 现 在 或 未 来 的 WMO 应 用 上 使 用 BUFR 是 一 种 描 述 表 驱 动 的 编 码 格 式, 其 特 点 是 : BUFR 是 一 种 自 描 述 的 数 据 格 式, 自 描 述 是 指 在 BUFR 信 息 中 的 格 式 和 数 据 的 内 容 就 是 由 BUFR 信 息 自 己 来 描 述 BUFR 采 用 描 述 符 和 码 表 相 结 合 的 方 式 来 描 述 数 据, 描 述 符 由 公 用 描 述 符 和 扩 充 描 述 符 组 成, 便 于 BUFR 信 息 的 扩 充 BUFR 数 据 由 二 进 制 组 成, 信 息 项 所 占 位 数 (bits) 可 不 相 同 1
BUFR 采 用 描 述 符 重 复 循 环 的 方 式 来 用 尽 量 少 的 节 字 数 描 述 尽 量 多 的 内 容 2 BUFR 格 式 引 入 航 空 气 象 中 2.1 WAFS 简 介 WAFS(World Area Forecast System) 即 世 界 区 域 预 报 系 统 它 是 一 套 由 国 际 民 航 组 织 (ICAO) 认 可 的 国 际 航 空 气 象 业 务 系 统 全 球 有 两 大 世 界 区 域 预 报 中 心 (WAFC) 建 立 了 WAFS 系 统, 它 们 分 别 是 伦 敦 中 心 和 华 盛 顿 中 心 的 WAFS 系 统, 两 个 系 统 互 为 备 份 现 在 北 京 上 海 和 广 州 安 装 了 WAFS 系 统, 其 接 收 的 世 界 区 域 预 报 产 品 为 国 内 的 民 航 飞 行 提 供 了 强 有 力 的 气 象 保 障, 在 国 际 航 线 上 的 航 班 飞 行, 这 些 产 品 更 是 唯 一 的 必 不 可 少 的 飞 行 气 象 保 障 资 料 2.2 WAFS 中 引 入 BUFR 格 式 2.2.1 T4 格 式 的 局 限 T4 是 G3 传 真 机 采 用 的 图 像 格 式, 采 用 CCITT 建 议 的 改 进 霍 夫 曼 码 编 码 (MH) 方 案 进 行 压 缩 由 于 传 真 图 为 二 像 素 ( 黑 和 白 ) 位 图 格 式, 采 用 MH 算 法 压 缩 后, 传 输 量 十 分 小, 因 此 早 期 的 气 象 无 线 传 真 广 播 系 统 普 遍 采 用 T4 格 式 来 传 输 图 形 WAFS 原 来 也 采 用 T4 格 式 来 保 存 和 广 播 航 路 重 要 天 气 预 告 图, 但 T4 格 式 存 在 无 法 编 辑 打 印 与 发 布 效 果 差 传 输 变 形 等 诸 多 局 限 2.2.2 WAFS 引 入 BUFR 格 式 由 于 T4 格 式 存 在 种 种 的 局 限, 因 此,WAFC 开 始 试 验 采 用 BUFR 格 式 来 保 存 和 交 换 航 空 气 象 图 形, 也 就 是 开 展 航 空 气 象 图 形 的 数 字 化 矢 量 化 工 作 WAFS 伦 敦 中 心 分 别 于 2001 年 和 2002 年 完 成 了 高 空 中 空 重 要 天 气 预 告 图 的 BUFR 码 化 WAFS 的 BUFR 格 式 现 已 发 展 到 3.0 版 本, 原 来 的 T4 格 式 也 已 被 淘 汰 使 用 2.2.3 BUFR 格 式 的 优 点 BUFR 格 式 相 对 于 T4 格 式 而 言, 有 其 特 定 的 优 势 : 2
资 料 传 输 量 缩 减 为 原 所 有 T4 图 的 一 半, 各 文 件 之 间 没 有 重 复 交 叉 的 内 容 用 户 可 以 根 据 需 要, 选 择 自 己 感 兴 趣 的 区 域, 定 做 重 要 天 气 预 告 图, 或 选 特 定 航 线 供 飞 行 计 划 之 用 可 以 实 现 矢 量 打 印, 可 以 对 图 中 的 气 象 元 素 进 行 编 辑 可 以 实 现 彩 色 图 表 显 示, 上 网 效 果 较 好 但 是, 使 用 BUFR 数 据 的 条 件 比 较 高, 收 到 BUFR 数 据 后 必 须 对 数 据 进 行 解 码, 然 后 必 须 编 写 专 门 的 图 形 软 件 进 行 可 视 化 处 理, 这 也 是 本 文 要 探 讨 的 重 点 3 WAFS BUFR 数 据 的 解 码 3.1 WAFS BUFR 数 据 文 件 WAFS BUFR 数 据 文 件 包 含 的 内 容 是 全 球 范 围 的 航 路 重 要 天 气 的 元 素,WAFS 一 天 发 布 四 个 时 次, 每 一 时 次 有 11 个 文 件, 每 一 个 文 件 包 括 一 到 几 种 气 象 元 素 这 些 文 件 包 含 了 高 层 (SWH) 和 中 层 (SWM) 的 重 要 天 气 预 告 图 的 内 容 3.2 WAFS BUFR 文 件 格 式 WAFS BUFR 文 件 和 通 用 的 BUFR 文 件 格 式 是 一 样 的, 由 六 个 部 分 组 成, 如 下 表 所 示 : 部 分 序 号 部 分 名 字 内 容 0 指 示 部 分 BUFR 标 识, 信 息 长 度 和 版 本 号 1 信 息 头 数 据 时 间 发 布 中 心 等 2 可 选 部 分 数 据 处 理 中 心 的 局 地 化 信 息 3 数 据 描 述 部 分 数 据 描 述 符 F + X + Y 4 数 据 部 分 真 正 的 数 据 5 结 束 标 记 7777 ASCII 码 3
WAFS BUFR 文 件 没 有 可 选 部 分 的 内 容 其 它 每 一 部 分 由 一 连 串 的 字 节 构 成, 每 一 单 独 部 分 总 是 由 偶 数 字 节 组 成, 需 要 时 将 补 上 0 值 来 凑 上 偶 数 字 节 任 何 一 个 BUFR 文 件 最 少 应 包 括 46 个 字 节 数 据 描 述 和 数 据 部 分 可 以 包 括 多 种 气 象 图 形 元 素, 使 用 顺 序 方 式 依 次 存 放 3.3 数 据 描 述 符 能 够 理 解 数 据 描 述 符 概 念, 也 就 基 本 上 能 够 理 解 BUFR 数 据 的 存 放 格 式 BUFR 信 息 是 一 种 自 描 述 的 信 息, 其 数 据 项 的 内 容 由 数 据 描 述 符 来 描 述 描 述 符 部 分 和 数 据 部 分 是 分 开 的, 但 描 述 符 会 和 数 据 相 对 应 一 个 描 述 符 总 共 占 两 个 字 节, 分 为 三 个 部 分, 如 表 1 所 示 : 表 1 FXY 描 述 符 F X Y 2 BITS 6 BITS 8 BITS F 表 示 描 述 符 本 身 的 内 容, 当 F 为 表 入 口 时, 会 到 相 应 的 表 去 检 索, 查 出 本 描 述 符 所 描 述 的 数 据 内 容 描 述 符 的 X 部 分 表 示 航 空 气 象 所 有 可 能 的 类 别, 如 温 度 风 等 参 数, 而 Y 部 分 则 是 X 部 分 参 数 的 入 口 值, 可 以 理 解 为 子 类 型 3.4 码 表 码 表 由 描 述 符 和 其 所 对 应 用 的 参 数 组 成 在 WMO FM94 BUFR 文 档 中, 码 表 分 为 多 个 表, 可 由 描 述 符 的 F 部 分 来 定 位 哪 一 个 表, 但 在 进 行 WAFS BUFR 的 解 码 时, 我 们 把 所 有 重 要 天 气 预 告 图 要 用 到 的 码 表 的 条 目 集 中 在 一 个 表 中, 这 样 做 的 好 处 是 表 中 的 所 有 项 都 有 用, 加 快 了 描 述 符 的 定 位 表 2 为 码 表 部 分 条 目 示 例 : 表 2 码 表 部 分 条 目 FXY 单 位 比 例 因 参 考 值 数 据 宽 度 说 明 子 (bit) 4
031001 数 值 0 0 8 有 多 少 个 描 述 符 将 被 复 制 001031 码 表 0 0 16 发 布 中 心 代 码 001022 CCITTIA5 0 0 224 天 气 现 象 名 字 004001 年 0 0 12 时 间 : 年 005002 度 2-9000 15 粗 精 度 纬 度, 负 数 为 南 半 球 008011 码 表 0 0 6 气 象 特 征 : 颠 簸 019006 米 /S 2 0 14 移 速 下 面 为 一 个 WAFS BUFR 的 重 要 云 区 元 素 的 描 述 符 组 成 : 112000 F=1, 重 复 标 志 ;X=12, 表 示 接 下 来 有 12 个 描 述 符 被 循 环, 外 框 重 复 031001 重 复 次 数, 从 数 据 区 获, 得 即 是 有 多 少 个 云 区 008011 外 框 开 始 重 复 数 据, 气 象 特 征,12:CLOUD 008007 空 间 ( 维 数 ) 指 示 码 1: LINE 007002 云 底 高 007002 云 顶 高 102000 F=1, 重 复 标 志 ;X=02, 接 下 来 有 02 个 描 述 符 被 循 环, 内 框 重 复 031001 重 复 次 数, 从 数 据 区 获 得, 即 是 云 区 的 点 数 005002 内 框 开 始 重 复 数 据 纬 度 006002 内 框 结 束 重 复 数 据 经 度 020008 云 散 程 序 020012 云 类 型, 如 CB 等 008007 008011 外 框 结 束 重 复 数 据 3.5 解 码 BUFR 数 据 的 解 码 由 描 述 符 码 表 和 取 出 的 数 据 共 同 完 成, 在 描 述 符 部 分 取 出 一 个 描 述 5
符, 比 如 取 出 的 FXY 为 005002, 到 码 表 中 定 位 005002, 根 据 表 2 可 知 这 是 个 粗 精 度 纬 度, 根 据 数 据 的 宽 度 15, 在 数 据 部 分 取 出 15 位 (bit) 的 数 值 (Value), 然 后 根 据 比 例 因 子 (Scale) 参 考 值 (Ref) 算 出 真 实 的 数 值 (Real Value) 计 算 公 式 如 下 : Real Value = (Value + Ref) / Power(10, Scale) 在 这 里 Power 为 幂 运 算, 假 设 在 数 值 部 分 (15bits) 取 出 的 值 为 12140, 则 真 实 的 值 : Real Value = (12140 + (-9000)) / Power(10, 2) = 31.4 得 出 北 纬 31.4 度 的 真 实 值 4 WAFS BUFR 数 据 的 可 视 化 4.1 可 视 化 的 基 本 功 能 WAFS BUFR 数 据 的 可 视 化, 也 即 是 要 把 WAFS BUFR 数 据 转 换 为 航 路 重 要 天 气 图 的 图 表, 是 一 项 非 常 复 杂 艰 巨 的 工 程, 它 涉 及 航 空 气 象 GIS 学 科 计 算 图 形 学 以 及 软 件 工 程 等 诸 多 知 识 我 们 经 过 多 年 潜 心 研 究,WAFS BUFR 数 据 的 可 视 化 功 能 已 基 本 上 实 现 其 主 要 功 能 必 须 包 括 : BUFR 重 要 天 气 预 告 图 底 图, 包 括 Mercator Lambert 和 极 射 赤 面 三 种 投 影 的 底 图 定 义, 范 围 选 取 和 底 图 其 它 格 式 定 制 将 解 码 出 来 的 数 据 按 气 象 元 素 的 格 式 转 换 为 可 视 图 形 画 到 底 图 上, 这 是 BUFR 可 视 化 必 有 的 基 本 功 能 可 对 气 象 元 素 进 行 编 辑 操 作 可 以 自 动 以 位 图 格 式 输 出, 也 可 以 按 矢 量 方 式 打 印 输 出 图 形 可 视 化 后 重 要 天 气 预 告 图 的 样 式 如 图 1 所 示 : 6
图 1 BUFR 数 据 可 视 化 图 形 4.2 底 图 地 球 投 影 的 生 成 底 图 的 地 球 投 影, 就 是 地 球 椭 圆 球 坐 标 系 和 平 面 坐 标 系 的 相 互 转 换 为 和 国 外 制 作 的 图 形 相 一 致,WAFS BUFR 的 底 图 以 WGS-84 坐 标 系 为 基 准 面 本 文 仅 介 绍 Mercator 投 影, 其 它 投 影 请 参 考 列 出 的 文 献 Mercator 投 影 是 一 种 等 角 正 轴 圆 柱 投 影,Mercator 投 影 没 有 角 度 变 形, 由 每 一 点 向 各 方 向 的 长 度 比 相 等, 它 的 经 纬 线 都 是 平 行 直 线, 且 相 交 成 直 角, 经 线 间 隔 相 等, 纬 线 间 隔 从 标 准 纬 线 向 两 极 逐 渐 增 大 所 以 纬 度 越 高, 变 形 越 大, 只 适 合 于 制 作 中 低 纬 度 的 底 图 下 面 给 出 Mercator 投 影 的 正 反 解 公 式 : 正 解 公 式 :(B,L) (X,Y), 标 准 纬 度 B0, 原 点 纬 度 0, 原 点 经 度 L0: X N π B 1 esin B = K ln[ tg( + )*( ) 4 2 1+ esin B e 2 ] YE = K( L L0) 7
2 a / b K = NB0*cos( B0) = *cos( B0) 2 2 1+ e *cos ( B0) 反 解 公 式 :(X,Y) (B,L), 标 准 纬 度 B0, 原 点 纬 度 0, 原 点 经 度 L0( 图 3 所 示 ) XN π ( ) ( )ln( K 2 1+ B = 2arctg( EXP * EXP 2 YE L = + L0 K e 1 esin B ) esin B ) 在 公 式 中,EXP 为 自 然 对 数 底 数,e 为 地 球 偏 心 率,a b 为 地 球 椭 球 长 短 半 轴 长,N B0 为 在 标 准 纬 度 B0 处 的 卯 酉 圈 曲 率 半 径,e 为 地 球 第 二 偏 心 率 反 解 公 式 中,B 的 计 算 通 过 迭 代 算 法, 可 以 很 快 收 敛 根 据 正 解 公 式, 给 出 Mercator 投 影 的 标 准 纬 度 原 始 经 度, 读 出 全 球 海 岸 线 的 纬 度 经 度, 则 可 算 出 在 底 图 范 围 内 的 海 岸 线 每 一 点 的 XY 坐 标, 再 经 过 平 移 和 缩 小, 画 在 底 图 上, 再 画 上 经 纬 线 和 标 注, 底 图 就 生 成 了 4.3 BUFR 气 象 要 素 可 视 化 及 编 辑 根 据 BUFR 数 据 解 码 出 来 的 经 纬 度 及 其 它 气 象 元 素, 通 过 投 影 坐 标 变 换 和 计 算 机 图 形 学 的 方 法, 将 其 绘 制 到 底 图 上, 这 就 是 可 视 化 ; 绘 制 出 的 气 象 图 形 还 应 可 通 过 人 工 操 作 的 方 式 进 行 编 辑 下 面 探 讨 几 个 关 键 的 技 术 4.3.1 平 滑 曲 线 的 实 现 平 滑 曲 线 是 绘 制 锋 面 急 流 轴 等 线 型 天 气 系 统 的 关 键 技 术 平 滑 曲 线 必 须 通 过 给 定 的 点, 插 值 加 入 其 它 的 点, 是 光 滑 无 拐 点 的 线 型 好 的 算 法 必 须 有 较 高 的 计 算 效 率, 而 且 改 变 一 个 点 位 置 时, 不 会 对 其 它 的 点 之 间 的 线 段 造 成 较 大 的 影 响 4.3.1.1 Akima 算 法 Akima 算 法 是 一 种 5 点 平 滑 二 次 多 项 式 插 值 算 法, 可 以 由 前 面 的 点 追 到 下 一 个 点 而 点 与 点 之 间 通 过 插 值 的 方 式 得 出 其 它 点, 再 通 过 连 线 实 现 平 滑 当 插 值 步 长 过 于 大 时, 曲 线 无 法 光 顺, 当 插 值 步 长 过 小 时, 造 成 计 算 量 增 大, 而 且 步 长 过 小 时, 反 而 会 出 现 曲 线 粗 糙 现 象 8
因 此, 必 须 根 据 线 长, 动 态 调 整 插 值 步 长 Akima 的 追 踪 算 法, 十 分 适 合 于 线 中 间 要 插 入 其 它 符 号 的 急 流 轴 绘 制, 当 追 踪 到 某 位 置 时, 可 先 画 出 其 它 气 象 符 号, 比 如 风 标 符 号, 再 接 着 跳 过 符 号 画 出 余 下 部 分 图 2 为 通 过 四 个 关 键 点, 采 用 Akima 算 法 画 出 急 流 轴, 并 加 入 风 标 的 例 子 图 2 Akima 算 法 画 出 的 急 流 轴 通 过 我 们 的 实 践 表 明,Akima 曲 线 算 法 能 很 好 地 应 用 于 曲 线 形 天 气 系 统 的 绘 制 4.3.1.2 Spline 算 法 样 条 (Spline) 算 法 是 一 种 多 次 多 项 式 的 平 滑 算 法, 该 算 法 比 较 常 用, 也 是 WAFC 在 指 导 BUFR 可 视 化 时 所 推 荐 的 方 法 之 一 该 算 法 虽 能 画 出 平 滑 的 曲 线, 但 画 出 的 曲 线 通 常 不 过 中 间 的 关 键 点, 经 我 们 多 次 研 究 后, 写 出 了 通 过 中 间 关 键 点 的 样 条 曲 线 程 序 图 3 为 用 该 算 法 画 出 的 样 条 曲 线 : 9
图 3 Spline 算 法 画 出 的 平 滑 曲 线 4.3.1.3 Bezier 曲 线 Bezier 曲 线 是 一 种 简 单 的 多 项 式 曲 线 算 法, 由 四 个 关 键 点 插 值 出 曲 线, Bezier 曲 线 不 适 合 于 画 多 点 平 滑 的 曲 线, 但 Windows 系 统 的 GDI 应 用 接 口 可 支 持 直 接 画 出, 编 写 程 序 简 单 但 由 Windows GDI 画 出 的 Bezier 曲 线 不 通 过 中 间 两 点, 应 用 性 差 我 们 利 用 给 定 的 中 间 点, 来 反 求 新 中 间 点 的 坐 标, 再 由 新 中 间 点 来 计 算 新 曲 线, 则 曲 线 可 通 过 给 定 的 中 间 点 在 画 重 要 云 区 曲 线 时, 可 以 采 用 Akima 曲 线 叠 加 Bezier 曲 线 的 方 法 来 实 现, 图 4 为 画 出 的 示 例 : 图 4 Akima + Bezier 算 法 画 出 的 云 区 曲 线 10
4.3.2 气 象 符 号 的 放 缩 及 旋 转 气 象 符 号 的 放 缩 和 旋 转, 可 以 达 到 手 工 调 整 气 象 符 号 的 目 的 气 象 符 号 应 该 可 以 定 义 大 小 长 宽 比 例 角 度 等 图 5 显 示 气 象 符 号 随 经 纬 度 旋 转 的 例 子 图 5 气 象 图 形 随 经 纬 度 旋 转 除 了 通 过 操 作 气 象 符 号 上 的 方 块 手 柄, 可 实 现 气 象 符 号 的 缩 放 外, 还 应 可 通 过 设 置 图 形 的 样 式 来 精 确 设 置 图 形 的 大 小 颜 色 和 旋 转 等 属 性 4.3.3 其 它 图 形 学 问 题 为 了 实 现 整 个 BUFR 数 据 的 可 视 化, 还 必 须 解 决 如 下 问 题 : 必 须 实 现 文 字 的 旋 转 通 过 图 形 分 层 的 原 理 和 底 图 缓 冲 的 方 式, 实 现 图 形 互 相 覆 盖 时 的 删 除 移 动 问 题 实 现 图 上 气 象 元 素 的 操 作 感 知 功 能 4.4 图 形 的 输 出 11
做 好 的 图 形 必 须 能 够 输 出 和 发 布, 这 是 BUFR 数 据 可 视 化 的 目 的, 通 常 有 两 种 输 出 方 式, 一 是 将 做 好 的 图 形 存 为 位 图 格 式, 发 布 到 网 上, 供 用 户 浏 览 ; 一 是 将 图 形 以 矢 量 的 方 式 输 出 到 打 印 机 上 5 BUFR 数 据 可 视 化 水 平 展 望 全 自 动 生 成 重 要 天 气 图 图 形,BUFR 航 空 数 据 的 可 视 化 工 作 还 存 在 一 些 难 题 : 注 释 方 框 会 出 现 相 互 重 叠 现 象, 即 使 是 由 WAFC 生 成 的 图 也 会 存 在 此 问 题 ; 注 释 和 天 气 现 象 的 结 合 程 度 比 较 松 散 由 于 BUFR 数 据 属 于 全 球 范 围 预 报 数 据, 而 制 作 图 形 时, 则 只 有 一 定 的 范 围, 必 然 存 在 气 象 符 号 或 天 气 系 统 符 号 在 图 边 被 裁 剪 问 题, 使 到 无 法 完 整 了 解 图 边 被 裁 剪 符 号 的 所 有 属 性 上 述 问 题, 只 有 经 过 逐 渐 提 高 BUFR 数 据 的 可 视 化 算 法, 才 能 更 加 完 美 地 实 现 BUFR 数 据 的 可 视 化 6 总 结 WAFS BUFR 航 空 数 据 的 解 码 和 可 视 化 是 近 年 来 航 空 气 象 的 一 个 重 要 研 究 课 题, 如 何 更 好 地 应 用 这 一 全 球 数 字 化 模 式 预 报 的 重 要 天 气 产 品, 来 为 日 益 扩 大 范 围 的 国 际 航 空 飞 行 服 务, 迫 切 地 摆 在 了 航 空 气 象 科 研 人 员 的 面 前, 本 文 通 过 对 WAFS BUFR 航 空 数 据 的 解 码 和 可 视 化 这 一 复 杂 课 题 进 行 探 讨, 以 期 达 到 抛 砖 引 玉 的 作 用, 切 实 提 高 应 用 WAFS 产 品 的 能 力 参 考 文 献 : (1) Guide to WMO Table Driven Code Forms:FM 94 BUFR and FM95 CREX Layer 1 2 3. (2) World Area Forecast Centres (WAFCs) London and Washington. Representing WAFS Significant Weather (SIGWX) Data in BUFR Version 3.0. (3) European Centre for Medium-Range Weather Forecasts. BUFR reference manual (2006) (4) A guide to the encoding into BUFR of Significant Weather Features extracted from a SIGWX chart, $Revision: 1.9 $ 12
(5) 孙 家 广. 计 算 机 图 形 学. 北 京 : 清 华 大 学 出 版 社,1986 年 12 月 (6) http://www.ossim.org/ 网 站 之 投 影 部 分 内 容 (7) International Association of Oil & Gas Producers(OGP). Surveying and Positioning Guidance. (8) 戴 勤 奋. 常 用 地 图 投 影 转 换 公 式. 青 岛 海 洋 地 质 研 究 所 13