我必须计算图像内一个区域（不总是矩形）的“低频带功率”。我必须这样做，因为我必须使用 Matlab 实现本文中解释的算法（参见“火灾检测算法”段落中的第 2 页）。所以我的问题是：
如何计算低频段功率？
我们的想法是遵循以下三个步骤：
- FI = 使用 fft2 计算二维傅里叶变换
- LBFI = 使圆外的所有值都为 0（低带通）
- PSD = 功率计算
在代码中
FI = fftshift(fft2(img));
LBFI = % not yet done
PSD = sum(sum(abs(LBFI).^2))
有没有办法让这个算法适应一个非矩形区域？
我认为没有办法做出类似的事情，但我希望是错的