CUDA中FFT图像上采样输出异常问题排查求助

补零策略常见坑点 #

• 补零位置错误：FFT域补零必须在频谱的**高频区域（矩阵边缘）**操作，而非直接在矩阵末尾追加。比如2×2转4×4，正确做法是把原2×2频谱放在4×4矩阵的中心区域，周围补零；要是把原数据堆在左上角再补零到右下角，必然导致采样后图像出现周期性重复或相位偏移。
• 实信号对称性被破坏：如果输入是灰度图（实信号），正向FFT结果天然满足共轭对称性。补零时必须严格维持这个对称性，否则逆FFT会产生冗余虚部，最终图像出现伪影。比如N×N实信号FFT，只需处理前N/2+1列的频谱，补零后要镜像复制对称部分。

CUDA FFT执行环节问题 #

• 内存对齐不达标：CUDA FFT的cufftExecR2C/cufftExecC2R要求输入输出内存是16字节对齐的，要是用普通cudaMalloc分配内存没指定对齐，大概率会引发计算错误。建议改用cudaMallocPitch或cufftAlloc分配符合要求的内存。
• FFT计划维度不匹配：创建FFT计划时，正向FFT要对应原图像尺寸，逆FFT必须严格对应补零后的目标尺寸，尺寸不匹配直接导致输出乱码。

归一化环节易犯错误 #

• 归一化因子搞混：逆FFT后的结果必须除以补零后图像的总像素数（目标尺寸宽×高），不是原图像的像素数。用原尺寸做归一化会导致图像整体亮度异常（过亮或过暗）。
• 虚部处理遗漏：要是补零破坏了实信号的共轭对称性，逆FFT后会产生非零虚部，必须取实部再做归一化，不然会引入噪声或异常色块。

代码快速排查步骤 #

• 打印正向FFT后的频谱矩阵，检查是否符合实信号的共轭对称性；
• 打印补零后的频谱矩阵，确认补零位置和对称性是否正确；
• 验证逆FFT后的实部数据范围，若超出0~255（8位图像），需做截断或映射处理；
• 检查所有CUDA/FFT API的返回码（比如cufftPlan1d/cufftExecR2C的返回值），确认没有执行错误。

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