matlab 中 FT 图像与 cuda 中的差异

Question

我正在尝试实现一个可以进行 2D 卷积的 matlab 代码。

matlab 代码首先将 3x3 内核置于图像大小的填充矩阵中。我在 C++ 代码中做了同样的事情。

我已将数组和 matlab 矩阵输出到 .csv 文件并确认它们是相同的。然后我对每一个都运行前向 FFT。在 matlab 中，内核的 FT 图像看起来就像您所期望的那样 - 基本上中间很强烈，呈圆形辐射。然而，CUDA FT 图像（我在将其导入为 csv 后在 matlab 中绘制）看起来像一个椭圆形。

可能是什么原因造成的？看起来几乎好像内核没有在图像中居中，但就像我说的，我将填充的内核数据转储到 csv 中，并使用 imagesc 在 matlab 中查看它，它看起来对我来说居中，事实上是与填充的 matlab 内核完全相同。

这是我用来将内核放在填充数组中心的代码：

kSize = 3;
halfl = 0.5*(kSize-1);

if(chipW%2 == 0)
    dcW = (.5*chipW) +1;
else
    dcW = round(chipW*.5);

if(chipH%2 == 0)
    dcH = (.5*chipH) +1;
else
    dcH = round(chipH*.5);

dcH--;
dcW--;


for(int i = dcH-halfl ; i <= dcH+halfl ; i++)

{
    for(int j = dcW -halfl ; j <= dcW+halfl ; j++)

    {
        h_PaddedKernel[i*chipW + j] = make_cuComplex(hp_kernel[(i-(dcH-halfl))*kSize + (j-(dcW-halfl))], 0.0);

    }


}

kSize 是内核一种大小的宽度，chipW 和 ChipH 是我尝试处理的图像的宽度和高度

Answer 1

有两件事：

• matlab 本身使用 double 进行计算。你的 CUDA 代码也在做同样的事情吗？
• 显示的 matlab 尝试在显示值之前缩放它们，图中有一些选项可以防止这种情况。