Question

我正在尝试创建一个非常简单的C ++程序，该程序在范围[0-100]中给出一个参数，将低通滤波器应用于灰度图像，该图像应该以预定方式“压缩”到给定参数的值。我正在使用FFTW库。

我对如何定义频率阈值cut有疑问。有没有更有效的方法来定义这样的价值？

//fftw_complex *fft
//double[] magnitude
// . . . 

int percent = 100;
    if (percent < 0 || percent > 100) {
        cerr << "Compression rate must be a value between 0 and 100." << endl;
        return -1;
    }

double cut =(double)(w*h) * ((double)percent / (double)100);
    for (i = 0; i < (w * h); i++) {
        magnitude[i] = sqrt(pow(fft[i][0], 2.0) + pow(fft[i][1], 2.0));
        if (magnitude[i] < cut) {
            fft[i][0] = 0.0;
            fft[i][1] = 0.0;
        }
    }

UPDATE1：

我已将代码更改为此，但我不确定这是过滤频率的正确方法。图像肯定会被压缩，但是非方形图像会混乱并且将压缩设置为100％并不是真正的最大可用压缩（我可以达到~140％）。 Here你可以找到我现在看到的图像。

int cX = w/2;
int cY = h/2;
cout<<"TEST "<<((double)percent/(double)100)*h<<endl;
for(i = 0; i<(w*h);i++){
    int row = i/s;
    int col = i%s;
    int distance = sqrt((col-cX)*(col-cX)+(row-cY)*(row-cY));
    if(distance<((double)percent/(double)100)*min(cX,cY)){
        fft[i][0] = 0.0;
        fft[i][1] = 0.0;
    }
}

Answer 1

这根本不是低通滤波器。低通滤波器通过低频，即它去除精细细节（模糊）。你显然需要2D FFT。

这段代码基本上只删除了随机位。

[编辑] 新代码看起来更像是低通滤波器。预期141％的设置：正方形的对角线是sqrt（2）=其侧面的1.41倍。将索引转换为行/列对应使用图像宽度，而不是一些随机无法解释的s。

我不知道零频率的位置。这应该很容易发现（最大值）但它可能在（0,0）而不是（w / 2，h / 2）

用FFTW在C中实现FFT低通滤波器

1 个答案: