基于分数阶滤波器的遥感图像处理.pdfVIP

第 1期，总第 88期 国 土 资 源 遥 感 No．1．2011 2011年3月 15日 REMOTE SENSING FOR LAND ＆RESOURCES Mar．．2011 基于分数阶滤波器的遥感图像处理 郭 丁 ，顾行发 ，余 涛 ，赵 辉 ，马红涛 ， (1．电子科技大学 自动化X-程学院，成都 610054；2．中国科学院遥感应用研究所 遥感科学国家重点实验室，北京 100101；3．国家航天局航天遥感论证中心，北京 100101) 摘要：在研究分数阶微积分的基础上，提出了一种基于WLS误差准则的滤波方法，该滤波器的系数可通过调节其 阶次来改变；给出了该滤波器的原理，分析了其幅频特性。通过用Riemann—Liouville定义验证，并将此滤波器运 用于遥感图像处理 ，结果表明，该方法在抗干扰和提取更多细节方面取得了很好的平衡 ，具有较好的滤波效果。 关键词 ：分数阶微积分；FIR；数字滤波器；WLS误差准则 ；遥感 中图法分类号：TP751．1 文献标识码：A 文章编号：1001—070X(2011)01—0048—04 FOI)问题 …。 0 引言 现实中的时间序列在频域上通常是带限的。可 分数阶微积分 (FractionalOrderCalculus，FOC) 以设计出这样的FIR滤波器，使其在时间序列的有 概念 的提 出可 以上溯 300a，早在整数阶微积分 效频带内能很好地逼近 FOC算子，那么通过该滤波 (IOC)产生时，FOC就产生了；在 19世纪，其理论 器的输入、输出处理，就可以在保证有限精度的前提 体系就已基本完善 ¨J。FOC作为一个纯数学分支 下实现时间序列的FOC运算。 已逐渐发展成体系，但其物理意义不明确。近年来， 与泰勒展开法的出发点不同，本文充分考虑了 多种信息提取与处理领域的应用研究 表明，作 时问序列的频域带限性质，它能够有效地控制滤波 为一种新颖而有效的工具，FOC拥有 IOC所不能替 器频率响应带来的计算误差，而且无论对于FOD还 代的作用。 是FOI都具有普适性。 很多数学家从不同角度给 FOC作出了多种不 同的定义 ，主要有Grumivald—Letinikov(G—L)、 l 滤波器原理与设计实现 Caputo(C)及 Riemann—Liouville(R—L)等，其中 R—L定义最具有普适性 J。然而，所有 FOC定义 1．1 时间序列FOC运算的FIR滤波器实现原理 的表达式都是复杂的，实现其计算比实现 IOC复杂 对于实数功率时间序列 (n)，其 IOC记为 得多 。可以说，探索FOC有效的计算方法是深化 (n)，P为整数，当P0时表示微分 ，P0时 和拓展其应用的关键。 表示积分。根据傅立叶分析理论， (n)的傅立 离散化的时问序列是表达诸多处理对象信息的 叶变换为 重要 形式。由于有 限冲激 响应滤 波器 (Finite ’(n)]=( ) ( ) (1) ImpulseResponse，FIR)的实现便利及稳定性，借助 式中，( )为 (n)的傅立叶变换， ∈[一盯， 其实现时间序列的FOC运算是优先考虑的途径之 盯]。由式(1)可以定义 IOC数字滤波器的频率响 一 。 有人使用了基于泰勒级数展开设计分数阶微分 应 ，即 (FractionalOrderDerivative，FOD)滤波 器 的方 ／4(P)(