计算fitEllipse返回的旋转椭圆上点的切线时部分出错如何解决？

核心错误原因 #

• 你硬编码将角度减90度、固定取height为长轴、width为短轴的逻辑仅适用于角度接近90°的场景，当椭圆角度较小时（比如你报错的6.27°场景），长轴的角度和a/b的对应关系完全错位，所以计算结果错误。
• OpenCV的fitEllipse返回的RotatedRect参数有特殊约定：angle是椭圆主轴顺时针旋转的角度，取值范围为[0,90)，size的两个分量分别对应两个主轴的长度，没有固定谁是长轴谁是短轴。
• 原实现用斜率表示直线还存在潜在问题：当切线为竖直线时斜率无穷大，会导致数值计算异常。

修复后的代码 #

// 优先返回直线一般式，避免斜率无穷大的兼容问题
struct LineGeneral {
    double A, B, C; // 满足 Ax + By + C = 0
};

static LineGeneral getTangent(const RotatedRect& ell, const Point& p) {
    // 统一长半轴、短半轴定义
    double a = max(ell.size.width, ell.size.height) / 2.0;
    double b = min(ell.size.width, ell.size.height) / 2.0;
    // OpenCV的angle是顺时针旋转角度，转换为标准逆时针旋转弧度
    double rad = -ell.angle * CV_PI / 180.0;
    double cr = cos(rad);
    double sr = sin(rad);
    // 椭圆中心坐标
    double s = ell.center.x;
    double t = ell.center.y;
    // 将点p转换到椭圆局部坐标系（中心在原点，长轴沿x轴）
    double x_local = (p.x - s) * cr + (p.y - t) * sr;
    double y_local = -(p.x - s) * sr + (p.y - t) * cr;
    // 局部坐标系下的切线方程：x_local*x/a² + y_local*y/b² = 1，转换回全局坐标系一般式
    double A = x_local * cr / (a*a) - y_local * sr / (b*b);
    double B = x_local * sr / (a*a) + y_local * cr / (b*b);
    double C = -1.0 - A * s - B * t;
    return {A, B, C};
}

// 如果必须保留原返回斜率的Line接口，加一层转换即可
static Line getTangentWithSlope(const RotatedRect& ell, const Point& p) {
    LineGeneral general = getTangent(ell, p);
    if (fabs(general.B) < 1e-6) {
        // 竖直线场景，根据你的Line类定义适配无穷大斜率的表示即可
        return Line(INFINITY, p);
    }
    double k = -general.A / general.B;
    return Line(k, p);
}

如果不想改动原逻辑的计算方式，仅替换原代码中rad、a、b的赋值逻辑也可以解决当前的报错问题。

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