【图像处理】FJFM 分数阶正交傅里叶矩图像重建附matlab代码

2026/6/16 13:17:43

✅作者简介：热爱科研的Matlab仿真开发者，擅长毕业设计辅导、数学建模、数据处理、程序设计科研仿真。

🍎完整代码获取定制创新论文复现点击：Matlab科研工作室

👇 关注我领取海量matlab电子书和数学建模资料

🍊个人信条：做科研，博学之、审问之、慎思之、明辨之、笃行之，是为：博学慎思，明辨笃行。

🔥 内容介绍

一、引言

在图像处理领域，图像重建及特征提取是关键任务。本文所提出的 FJFM（分数阶雅可比 - 傅里叶矩）凭借其独特的通用性质和时频分析能力，为图像重建和特征提取提供了新的视角。同时，基于 FJFM 的时频特性构建的 MLMF（混合低阶矩特征）框架，旨在提升图像全局表示的鲁棒性和判别能力。

二、FJFM 的特性

（一）通用性质

FJFM 是现有基于雅可比多项式的经典及分数阶圆盘基连续正交矩的通用形式。这一性质意义重大，它为分数阶正交矩的研究提供了统一的数学工具。以往不同形式的正交矩研究往往各自独立，而 FJFM 的出现使得可以在统一框架下对它们进行深入分析和比较，有助于发现不同正交矩之间的内在联系，推动分数阶正交矩理论的整体发展。

（二）时频分析能力

FJFM 能够通过改变分数参数的值来调整径向核的零分布。这种独特的特性在图像处理中具有重要应用价值。例如，在处理信息抑制问题时，可根据需要调整分数参数，使径向核的零分布能够更好地抑制图像中特定频率的信息，从而突出感兴趣的部分。在提取图像局部特征方面，通过合理选择分数参数，利用径向核零分布的变化，可以更准确地捕捉图像局部区域的细节特征，为后续的图像识别、分类等任务提供更具代表性的特征信息。

三、FJFM 计算的问题与解决策略

在实际实现中，基于多项式的矩计算通常涉及一些阶乘和 / 或伽马项，这可能会导致数值不稳定性和较高的计算成本。为解决这一问题，引入了一种新的递归策略。

该递归策略通过巧妙地利用多项式之间的递推关系，避免了直接计算阶乘和伽马项带来的数值问题。具体来说，它基于前一步的计算结果，通过简单的乘法和加法运算得到下一步的矩值，大大减少了计算量。这种递归方式不仅提高了计算效率，还增强了数值稳定性，使得 FJFM 在实际应用中能够更可靠地进行图像重建和特征提取。

四、MLMF 框架

（一）设计思路

MLMF 是一种基于 FJFM 时频特性的全新框架，旨在提升图像全局表示的鲁棒性和判别能力。其核心策略是将具有不同分数参数的低阶矩组合成单个特征向量，而非将它们作为单独的特征使用。

（二）优势分析

提升判别能力：分数参数与 FJFM 基函数的时域特性相关。不同分数参数的低阶矩能够捕捉图像不同方面的特征信息，将它们组合在一起，可以形成更具判别性的特征向量。例如，某些分数参数下的低阶矩可能对图像的纹理特征敏感，而另一些则对图像的边缘特征更有效，组合后能够全面地反映图像的特征，提高在图像分类等任务中的判别能力。
保证鲁棒性：只考虑低阶矩，即低频率分量。低频率分量对图像的噪声、光照变化等干扰因素相对不敏感，因此保证了特征的鲁棒性。在实际应用中，图像可能会受到各种噪声和环境因素的影响，基于低阶矩的 MLMF 框架能够在这些复杂情况下仍保持较好的性能，准确地表示图像的特征。

五、总结与展望

FJFM 的通用性质和时频分析能力为图像重建和特征提取提供了有力的理论支持，而新引入的递归策略有效解决了其计算过程中的难题，使得 FJFM 能够更高效、稳定地应用于实际。基于 FJFM 时频特性构建的 MLMF 框架，通过巧妙地组合低阶矩，提升了图像全局表示的鲁棒性和判别能力。

未来的研究可以进一步探索 FJFM 在更复杂图像场景下的应用，如在医学图像、遥感图像等领域的特征提取和图像重建。对于 MLMF 框架，可以研究如何自适应地选择分数参数和低阶矩的组合方式，以更好地适应不同类型图像的特点。同时，结合深度学习等新兴技术，将 FJFM 和 MLMF 与神经网络相结合，有望进一步提升图像处理的性能和效果。

⛳️ 运行结果

📣 部分代码

function [output] = getCHFM_RadialPoly(order,rho)% obtain the order and repetitionn = order;output = zeros(size(rho)); % initilization% compute the radial polynomialfor s = 0:floor(n/2)c = ((-1)^s)*factorial(n-s) / ...(factorial(s)*factorial(n-2*s));output = output + c * (4*rho-2) .^ (n-2*s);endoutput=output.*((64*(1-rho)./(pi*pi*rho)).^(1/4))*sqrt(1/(2*pi));end % end getRadialPoly method