围绕某点旋转点时COS/SIN计算异常，求排查解决

问题描述 #

实现点绕中心点旋转功能时，参考的代码返回结果完全不符合预期。例如将点{6.13,0}绕{0,0}旋转90度，预期结果为{0,6.13}，但实际得到{3.75341847446559E-16,3.75341847446559E-16}。

原无效代码如下：

static Point RotatePoint(Point pointToRotate, Point centerPoint, double angleInDegrees)
{
    double angleInRadians = angleInDegrees * (Math.PI / 180);
    double cosTheta = Math.Cos(angleInRadians);
    double sinTheta = Math.Sin(angleInRadians);

    // 将中心点转为原点
    pointToRotate.X -= centerPoint.X;
    pointToRotate.Y -= centerPoint.Y;

    // 旋转计算（错误逻辑）
    pointToRotate.X = pointToRotate.X * cosTheta - pointToRotate.Y * sinTheta;
    pointToRotate.Y = pointToRotate.X * sinTheta + pointToRotate.Y * cosTheta;

    // 转回原坐标系
    pointToRotate.X += centerPoint.X;
    pointToRotate.Y += centerPoint.Y;

    return pointToRotate;
}

核心错误原因 #

代码中直接覆盖了偏移后的X值，导致计算Y坐标时使用的是已经修改过的旋转后X，而非原始偏移X。这是典型的变量复用错误——旋转公式需要基于同一个原始偏移坐标进行计算，不能中途修改其中一个值后再用它计算另一个。

以90度旋转的测试点为例：

1. 偏移后X=6.13，Y=0
2. 计算新X：6.13*0 - 0*1 = 0（此时X被改为0）
3. 计算新Y：0*1 + 0*0 = 0（错误地使用了修改后的X=0，而非原始的6.13）

修复后的代码 #

将偏移后的X、Y保存到临时变量，用临时变量完成旋转计算，避免覆盖原始值：

static Point RotatePoint(Point pointToRotate, Point centerPoint, double angleInDegrees)
{
    double angleInRadians = angleInDegrees * (Math.PI / 180);
    double cosTheta = Math.Cos(angleInRadians);
    double sinTheta = Math.Sin(angleInRadians);

    // 计算偏移量并保存到临时变量
    double offsetX = pointToRotate.X - centerPoint.X;
    double offsetY = pointToRotate.Y - centerPoint.Y;

    // 使用临时变量进行旋转计算
    double rotatedX = offsetX * cosTheta - offsetY * sinTheta;
    double rotatedY = offsetX * sinTheta + offsetY * cosTheta;

    // 转回原坐标系并返回新Point
    return new Point(rotatedX + centerPoint.X, rotatedY + centerPoint.Y);
}

测试验证 #

用你的测试数据验证：

• 旋转90度，中心点{0,0}，点1{6.13,0}：
  计算得rotatedX=6.13*0 - 0*1=0，rotatedY=6.13*1 +0*0=6.13，结果为{0,6.13}，符合预期。
• 点2{0.13,4.25}：
  rotatedX=0.13*0 -4.25*1= -4.25，rotatedY=0.13*1 +4.25*0=0.13，结果为{-4.25,0.13}，符合逆时针旋转90度的预期。

补充说明 #

1. 坐标系方向：上述代码基于数学坐标系（Y轴向上）的逆时针旋转。如果是屏幕坐标系（Y轴向下），旋转方向会相反，此时需要将sinTheta取反，或者调整旋转公式为顺时针版本：

   double rotatedX = offsetX * cosTheta + offsetY * sinTheta;
   double rotatedY = -offsetX * sinTheta + offsetY * cosTheta;
   

2. 浮点数精度：测试结果中的3.75341847446559E-16是浮点数计算的精度误差，实际可近似为0，属于正常现象。

内容的提问来源于stack exchange，提问作者ds900