【Cesium风场】如何解决风向角度定义不一致导致的风向计算误差

如何解决风向角度定义不一致导致的风向计算误差

背景

在处理风场数据时，我们通常需要将风速和风向转换为 U 和 V 分量，用于进一步的计算和可视化。然而，在这一过程中，风向角度的定义方式可能导致计算结果出现严重误差。本文将详细介绍这一问题的根源，并提供解决方案，确保风向计算的准确性。

问题描述

在气象学中，风向角度是按照从北向顺时针旋转定义的。例如：

• 0° 表示北风（从北吹来）
• 90° 表示东风（从东吹来）
• 180° 表示南风（从南吹来）
• 270° 表示西风（从西吹来）

然而，在数学中，角度的定义方式不同：

• 0° 通常表示正东方向（x轴正方向）
• 角度按逆时针旋转计算
• 90° 表示北方，180° 表示西方，270° 表示南方

在使用三角函数（sin 和 cos）计算 U 和 V 分量时，如果直接使用气象学定义的风向角度，而不进行转换，可能会导致严重的计算误差。

风场计算公式

在风场计算中，我们使用以下公式将风速和风向转换为 U 和 V 分量：

• U 分量计算：表示东西方向的风速（东为正方向，西为负方向）：
  $$U=-V_{\text{风速}}\times \sin ({\theta }_{\text{风向}})$$

• V 分量计算：表示南北方向的风速（北为正方向，南为负方向）：
  $$V=-V_{\text{风速}}\times \cos ({\theta }_{\text{风向}})$$

其中，θ风向是气象学中定义的风向角度，需要在使用前转换为数学角度。

误差的表现

以下是一个错误的风向转换导致的结果示例：

// 错误的 U 和 V 分量计算（直接使用气象学角度）
const u = -windSpeed * Math.sin(windDirectionRad); 
const v = -windSpeed * Math.cos(windDirectionRad);

假设原始风向数据的加权平均值为 206.57°，使用错误的公式计算得到的结果却是 62.28°，风向完全偏离实际值。这是由于没有考虑到气象学角度和数学角度的定义差异。

解决方案

为了正确地进行 U 和 V 分量的计算，我们需要将气象学角度转换为数学角度。转换公式如下：

const mathDirection = (270 - windDirection + 360) % 360;

该公式将气象角度转换为数学角度，从而确保三角函数使用的角度与风向定义保持一致。

修改前后的代码对比

错误的代码示例：

// 错误的风向转换，未考虑角度定义差异
const u = -windSpeed * Math.sin((windDirection * Math.PI) / 180);
const v = -windSpeed * Math.cos((windDirection * Math.PI) / 180);

正确的代码示例：

// 将气象学角度转换为数学角度
const convertToMathDirection = (degree: number): number => {
   
    return (270 - degree + 360)<