基于OpenCV的C++旋转后自定义双线性插值出现伪影问题求助

嘿，我明白自己实现双线性插值时遇到伪影有多头疼——尤其是看着Lena图出现奇怪的色块，肯定特别挠头。结合你提到的用了cvRound()而非ceil()、手动处理RGB分量的情况，我来梳理几个最可能踩坑的点，以及对应的排查和修复方法：

1. 先确认坐标映射方向是否搞反了 #

这是旋转插值里最容易犯的低级错误：很多人会下意识遍历输入图像像素，再映射到输出图像，但正确的做法应该是遍历输出图像的每个像素，反向计算它在原始输入图像中的对应坐标。如果搞反方向，会出现大量输出像素未被覆盖或重复覆盖的情况，直接导致错位、色块伪影。

举个简单逻辑：假设要旋转θ角，先拿输出图像的像素坐标(x_out, y_out)，先平移到以输出中心为原点的坐标系，做逆旋转（-θ角），再平移回原始图像的中心坐标系，得到浮点型的(x_in, y_in)，最后对这个浮点坐标做双线性插值。

2. cvRound()的使用是否导致邻域像素选择错误 #

双线性插值需要的是浮点坐标周围的四个整数邻点（比如floor(x_in)、ceil(x_in)对应的四个点），但cvRound()是四舍五入逻辑，这可能让你选到错误的邻域点，进而打乱权重计算。

比如x_in=2.4时，cvRound()得到2（和cvFloor()结果一致），但x_in=2.6时cvRound()得到3，这时候如果用cvRound(x_in)作为邻点，会混淆左右/上下的权重对应关系。正确的邻点获取方式应该是：

float x_in = ...; // 反向映射得到的输入图像浮点坐标
float y_in = ...;
int x1 = cvFloor(x_in);
int x2 = x1 + 1;
int y1 = cvFloor(y_in);
int y2 = y1 + 1;

再计算对应权重：

float dx = x_in - x1; // 水平方向距离权重
float dy = y_in - y1; // 垂直方向距离权重

3. 边界处理：是否忽略了坐标越界的情况 #

当反向映射得到的x_in或y_in超出原始图像范围（比如小于0，或大于等于图像宽/高），直接访问这些坐标会读取到内存中的垃圾数据，表现为随机的异常色块。

解决方法是对邻点坐标做边界钳位，确保它们在合法范围内：

int img_width = src.cols;
int img_height = src.rows;
// 钳位到0~宽/高-1的合法区间
x1 = max(0, min(x1, img_width - 1));
x2 = max(0, min(x2, img_width - 1));
y1 = max(0, min(y1, img_height - 1));
y2 = max(0, min(y2, img_height - 1));

你也可以选择镜像、填充等边界处理方式，但钳位是最简单直接避免非法内存访问的手段。

4. RGB分量处理是否存在精度丢失或计算错误 #

你提到对每个RGB分量单独处理，这里要注意：如果输入图像是CV_8UC3类型，每个分量是0-255的uchar，直接做浮点计算前必须先转换成float，否则会出现整数溢出或截断错误。

正确的分量处理逻辑示例：

// 先把四个邻点的RGB分量转成float类型，避免整数运算溢出
float b1 = src.at<Vec3b>(y1, x1)[0];
float g1 = src.at<Vec3b>(y1, x1)[1];
float r1 = src.at<Vec3b>(y1, x1)[2];

float b2 = src.at<Vec3b>(y1, x2)[0];
float g2 = src.at<Vec3b>(y1, x2)[1];
float r2 = src.at<Vec3b>(y1, x2)[2];

float b3 = src.at<Vec3b>(y2, x1)[0];
float g3 = src.at<Vec3b>(y2, x1)[1];
float r3 = src.at<Vec3b>(y2, x1)[2];

float b4 = src.at<Vec3b>(y2, x2)[0];
float g4 = src.at<Vec3b>(y2, x2)[1];
float r4 = src.at<Vec3b>(y2, x2)[2];

// 计算双线性插值后的分量
float b_interp = b1*(1-dx)*(1-dy) + b2*dx*(1-dy) + b3*(1-dx)*dy + b4*dx*dy;
float g_interp = g1*(1-dx)*(1-dy) + g2*dx*(1-dy) + g3*(1-dx)*dy + g4*dx*dy;
float r_interp = r1*(1-dx)*(1-dy) + r2*dx*(1-dy) + r3*(1-dx)*dy + r4*dx*dy;

// 转换回uchar前，钳位到0-255范围避免溢出
dst.at<Vec3b>(y_out, x_out)[0] = cvRound(clamp(b_interp, 0.0f, 255.0f));
dst.at<Vec3b>(y_out, x_out)[1] = cvRound(clamp(g_interp, 0.0f, 255.0f));
dst.at<Vec3b>(y_out, x_out)[2] = cvRound(clamp(r_interp, 0.0f, 255.0f));

如果直接用uchar做乘法加法，很容易因为溢出导致颜色完全异常——比如两个200的uchar相加会变成400，截断后变成144，完全偏离插值预期。

最后检查旋转矩阵的逆变换是否正确 #

旋转的坐标变换很容易在符号、平移步骤上出错。正确的逆变换步骤应该是：

1. 将输出像素坐标(x_out, y_out)平移到以输出图像中心为原点的坐标系：x_out_centered = x_out - dst_center_x; y_out_centered = y_out - dst_center_y;
2. 做逆旋转（旋转-θ角）：x_in_centered = x_out_centered * cosθ + y_out_centered * sinθ; y_in_centered = -x_out_centered * sinθ + y_out_centered * cosθ;
3. 平移回原始图像的中心坐标系：x_in = x_in_centered + src_center_x; y_in = y_in_centered + src_center_y;
   如果这里的旋转方向、符号搞反，不仅会导致图像旋转方向错误，还可能间接引发坐标越界或插值错误，产生伪影。

你可以对照上面的几点逐一排查，尤其是坐标映射方向和边界处理这两点，大概率是伪影的根源。

内容的提问来源于stack exchange，提问作者Štěpán vrša