用快速傅立叶变换加卡尔曼滤波方法进行图象恢复.pdfVIP

第五卷第三期 蒜AL号PR Es理1Na 989 l989年 9月 8I 长 用快速傅立叶变换加卡尔曼滤波 的方法进行图象的恢复 曹 尔 强 (长春邮电学院) (搞謇】本文叔述对变模糊井附加了噪声的二堆攀宁化胃錾 行俊复处 一种方法。它是在行方向 上进行一维的快速博立叶变换，而在列方向上进行一维卡焦曼i鼠逵．然后再复原的方法。由于告尔曼谗渡 的 “动态 性优点．使它从理论上解决了实时图象灶理的闸题。 一 ， 引 言 观测臻I或记录到的图象不仪因为各种随机观测噪声 (如胶卷微粒不匀、电子散弹噪声， 热噪声，量化噪声等)变坏，而且也因相机移动、散焦、大气拢动等造成模糊。利用退化模 型的某种先验的知识，把已经退化了的图象加 以恢复或重建，具有重要的现实意义。对静止 图象的恢复，已经研究的较多，而且 已经取得了很多喜人 的成果。但在某些应 用场合，则需 要对实时输入的数字化图象进行实时处理。美国科学家J．W ．W oods ’等人开发的二维卡 尔曼滤波，荷兰科学家J．Biemond_2_等人，f发ff々、本文在这里引用的方法，就是在进行这 种尝试。笔者根据J．Biemond等人的论文，进行了计算机上模拟处理，取得 了初 步 的成 果．这种方法可 用硬件 电路实现， 达到 “宴时处理”的 目的。因此，具有重大的现实意 艾 。 二，数学模型和原理框 图 假设 图象为零均值，均匀的离散Ⅳ ×N随机场 则可用下述表达式表示 (t，J)= ∑ 4(，q)x(t— ， g)+(i，J) (1) P，q∈ 其中， ： {，，glO≤，≤，3’一g≤‘g≤g2且 (p，g)≠ (O，0)} uCi，J)是噪声输入项。 如果成象系统是线性和空不变怕，并且假设噪声被 口于系统 i输 出端而点 散画数议波 眨周围少量的点，则被观察或记录到的图象能被表达成t kl 12 Y(I，J)= ∑ ∑ g(，z)(一 ，J—)¨ (，J) L2) k=kJl=ll 其中，∞(t，j)是与 (，i)无关的霉均值 白噪声； 第三期 用快速傅立叶嗖换加卡身：曼滤涟的方法进行目象的协复 ’l5 ‘ g(k， )是定义区间内∞脉冲响应或造成模糊的 扩散函数，并假设这种模糊波及的 城大大小于图象的尺寸。 我们知道，摄象机 一实际扫描情况是每次扫描一整行，如果我们定义I Ⅱ(0，0)：一1} ● (i)： ( (t，1)，#(，2)， … … ， (‘，Ⅳ)] ‘ Y()= (( 1)， (}，2)， … … ， ，(t，N )] U ()= [ (， 。)， (‘，2)， … … ， ( N )) W ()‘： (】(1，1)， (t，2)， … … ，∞ (， Ⅳ )] 并设置边界条件为0，则 (I) ， 2‘)式可改写成1 ．4o ()： 一∑ A， (1一 )+ ()‘ P =I (3) k 一 y()二∑ X(i )十W(t)