如何查找MATLAB freqz2在Python中的等效图像处理库？

• 首选工具：scipy.signal.freqz2
  该工具是官方提供的二维频率响应计算接口，功能和MATLAB的freqz2完全对齐，专门支持二维滤波器的频率响应计算，完全适配图像处理场景，无需自己额外封装逻辑。基础用法示例如下：

import numpy as np
from scipy import signal

# 示例输入：3x3的二维高斯滤波器核
filter_kernel = np.array([
    [1, 2, 1],
    [2, 4, 2],
    [1, 2, 1]
]) / 16

# 计算64x64点的二维频率响应
freq_grid, freq_resp = signal.freqz2(filter_kernel, fs=[64, 64])

返回值中freq_grid是归一化频率网格，freq_resp是对应位置的复数频率响应，输出逻辑和MATLAB freqz2完全一致。

• 自定义实现方案：基于numpy.fft.fft2
  如果有定制化需求，也可以手动通过二维FFT实现等效功能，核心逻辑是对补零到指定尺寸的滤波器核做二维FFT，再调整零频率位置即可，示例代码如下：

import numpy as np

def custom_freqz2(kernel, fft_size=(128, 128)):
    # 把滤波器核补零到指定FFT尺寸
    padded = np.zeros(fft_size, dtype=np.float64)
    padded[:kernel.shape[0], :kernel.shape[1]] = kernel
    # 做二维FFT并将零频率点移动到数组中心
    freq_response = np.fft.fftshift(np.fft.fft2(padded))
    return freq_response

如果需要和MATLAB输出完全匹配，调整FFT的归一化参数、补零对齐规则即可。

内容的提问来源于stack exchange，提问作者Facundo Farall