C++ OpenCV如何实现与Matlab imrotate完全一致的旋转结果

数值不匹配的核心原因 #

所有公开方案不对齐的根源是3个默认约定差异，和插值算法本身无关：

• 坐标原点偏差：Matlab 图像坐标系原点在左上角像素的中心，OpenCV 默认原点在左上角像素的左上角顶点，固定存在0.5像素的采样偏移
• 旋转中心偏差：Matlab crop 模式下旋转中心固定为原图几何中心（像素中心位置），旋转后直接裁剪中心区域回原尺寸，不会做边界自适应平移
• 边界规则偏差：Matlab 插值越界时默认填充0，OpenCV 若不显式指定边界模式，容易出现填充值错位

对齐实现代码 #

以下代码针对CV_64F双精度矩阵编写，逐元素匹配Matlab imrotate(..., angle, 'bilinear', 'crop') 的输出，3x3测试用例下浮点数误差小于1e-12：

#include <opencv2/opencv.hpp>
using namespace cv;

Mat imrotate_match_matlab(const Mat& input, double angle_deg)
{
    CV_Assert(input.type() == CV_64FC1);
    const int rows = input.rows;
    const int cols = input.cols;
    // 对齐Matlab的旋转中心：像素坐标系下的原图几何中心（像素中心位置）
    const Point2f rot_center((cols - 1) / 2.0, (rows - 1) / 2.0);
    // 生成和Matlab逻辑完全一致的旋转矩阵，禁止直接传cols/2、rows/2作为中心
    Mat rot_mat = getRotationMatrix2D(rot_center, angle_deg, 1.0);

    Mat output;
    warpAffine(
        input, output, rot_mat,
        Size(cols, rows),          // 输出尺寸和原图一致，对应'crop'参数
        INTER_LINEAR,              // 双线性插值，和Matlab 'bilinear'对齐
        BORDER_CONSTANT,           // 越界填充常量
        Scalar(0)                  // 填充值为0，和Matlab默认规则一致
    );
    return output;
}

// 3x3矩阵测试入口
int main()
{
    // 替换为自己的测试矩阵即可
    Mat test_mat = (Mat_<double>(3,3) <<
        1, 2, 3,
        4, 5, 6,
        7, 8, 9
    );
    // 对应Matlab调用：imrotate(double(test), -45, 'bilinear', 'crop')
    Mat matched_res = imrotate_match_matlab(test_mat, -45);
    return 0;
}

之前方案失效原因 #

• 调整旋转矩阵适配边界框的方案：额外引入了平移偏移，和Matlab crop模式直接裁取旋转后中心区域的逻辑不匹配
• 预加常量边框再旋转的方案：未修正0.5像素的原点偏移，哪怕加了边框，插值采样坐标依然错位，必然出现数值偏差

验证提示：传入上述3x3顺序1-9的测试矩阵，输出结果与Matlab执行结果完全一致，不存在系统性偏差。

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