图像复原最新课件

《图像复原最新课件》由会员分享，可在线阅读，更多相关《图像复原最新课件（98页珍藏版）》请在装配图网上搜索。

1、图像复原 最新 第5讲 图像复原力求保持图像的本来面目，以保真原则为前提，找出图像降质的原因，描述其物理过程，提出数学模型。根据该模型重建或恢复被退化的图像。不考虑图像是如何退化的，而是试图采用各种技术来增强图像的视觉效果。因此，图像增强可以不顾增强后的图像是否失真，只要看得舒服就行。n图像增强主要是主观过程，而图像复原主要是客观过程n图像的钝化或平滑可看作是图像复原技术（去噪）图像复原 最新第5讲 图像复原n5.1 图像退化/复原过程的模型n5.2 噪声模型n5.3 噪声存在下的图像复原n5.4 估计退化函数n5.5 图像复原的滤波方法n5.6 图像复原的非线性方法n5.7 图像复原的盲反卷

2、积方法n5.8 几何校正图像复原 最新进化与退化图像复原 最新图像退化因素图像退化指由场景得到的图像没能完全地反映场景的真实内容，产生了失真等问题w 透镜象差/色差,聚焦不准等 w 噪声（是一个统计过程）w 照相机/摄像机与目标的相对移动，如抖动（机械、电子）图像复原 最新图像退化模型H：退化过程n(x,y)：加性噪声（统计特性已知）退化过程描述：x,yH()x,yf()n()x,yg),(),(),(yxnyxfHyxg图像复原 最新n假设退化函数H为线性、位置不变性过程，则图像退化过程可在空间域表示为：h(x,y)是退化函数的空间描述n空间域上的卷积等同于频域上的乘积，则频域的等价描述为：

3、,*,g x yh x yf x yn x y vuNvuFvuHvuG,h(x,y)和 n(x,y)图像恢复：对原始图像作出尽可能好的估计。已知退化图像，要作这种估计，须知道退化参数的有关知识。图像复原 最新第5讲 图像复原n5.1 图像退化/复原过程的模型n5.2 噪声模型n5.3 噪声存在下的图像复原n5.4 估计退化函数n5.5 图像复原的滤波方法n5.6 图像复原的非线性方法n5.7 图像复原的盲反卷积方法n5.8 几何校正图像复原 最新n图像获取过程的噪声n如，使用CCD摄像机获取图像，光照水平和传感器温度是生成图像中产生大量噪声的主要因素n图像传输过程的噪声n如，通过无线电网络传

4、输的图像可能会因为光或其他大气因素的干扰被污染 图像噪声的来源图像复原 最新噪声的空间和频率特性n噪声与图像的相关性n相关乘性噪声n不相关加性噪声n图像平面上不同点的噪声是不相关的，其谱密度为常数。n一般假设图像上的噪声是白噪声。n实用上，只要噪声带宽远大于图像带宽，就可把它当作白噪声。图像复原 最新一些重要噪声的概率密度函数n高斯噪声n瑞利噪声n伽马（爱尔兰）噪声n指数分布噪声n均匀分布噪声n脉冲噪声（椒盐噪声）图像复原 最新噪声概率密度函数,1、高斯噪声、高斯噪声222)(exp21)(zzp高斯函数的曲线如左图所示。当z服从高斯分布时，其值有70%落在范围 内，且 有 9 5%落 在 范

5、围 。2,2图像复原 最新瑞利噪声n瑞利噪声的概率密度函数：n概率密度的均值和方差：azeazbazbazzp2204ba442b图像复原 最新伽马（爱尔兰）噪声n伽马噪声PDF：n其中，a0，b为正整数且“！”表示阶乘。其密度的均值和方差为：101!00bbaza zezbzpzab22ab图像复原 最新指数分布噪声n指数噪声的PDF：n其中，a0。概率密度函数的期望值和方差：n注意，指数分布的概率密度函数是当b=1时爱尔兰概率分布的特殊情况。000azzaezp za1221a图像复原 最新均匀分布噪声n均匀分布噪声的概率密度：n概率密度函数的期望值和方差是：bzaabzp其他102ba

6、1222ab 图像复原 最新脉冲（椒盐噪声）噪声n脉冲噪声的PDF是：n如果ba，灰度值b在图像中将显示为一个亮点，相反，a的值将显示为一个暗点。若 或 为零，则脉冲噪声称为单级脉冲。如果 和 均不可能为零，尤其是他们近似相等时，脉冲噪声值将类似于随机分布在图像上的胡椒和盐粉微粒。azbzPPabzp其他0aPbPaPbP图像复原 最新样本噪声图像和它们的直方图图像复原 最新各种噪声的产生n 高斯噪声：电子噪声、弱光照/温度条件下的传感器噪声n 瑞利分布：深度成像、超声波图像n 指数和Gamma分布：激光成像n 椒盐噪声：快速瞬变、误切换图像复原 最新周期噪声产生原因：图像采集过程中的电子或电

7、磁干扰消除方法：频率域滤波Butterworth带阻滤波器图像复原 最新周期噪声参数的估计n 一般可以通过图像的频谱进行估计n 特殊情况下可以直接从图像中噪声分量的周期性进行推断（简单情形）n 自动分析在下列情形可行：噪声信号非常明确、对噪声频率成分的位置有一定的先验知识图像复原 最新一般噪声参数的估计n 可以根据所采用的传感器类型进行噪声分布的部分推断；通常通过特定的成像安排进行估计n 当只有已采集到的图像时，一般通过图像中的平滑区域进行PDF参数的估计图像复原 最新噪声的确定 不同的恢复方法需要关于噪声的不同的数字特征。eg:维纳滤波要知道噪声的谱密度。约束最小平方滤波要知道噪声的方差。图

8、像复原 最新第5讲 图像复原n5.1 图像退化/复原过程的模型n5.2 噪声模型n5.3 噪声存在下的图像复原n5.4 估计退化函数n5.5 图像复原的滤波方法n5.6 图像复原的非线性方法n5.7 图像复原的盲反卷积方法n5.8 几何校正图像复原 最新n当在图像中唯一存在的退化是噪声时，则退化方程为：n噪声项是未知的，从g(x,y)或G(u,v)中减去它们不是一个现实的选择。,g x yf x yn x yvuNvuFvuG,图像复原 最新5.3.1 均值滤波器n算术均值滤波器n几何均值滤波器n谐波均值滤波器n逆谐波均值滤波器图像复原 最新n算术均值滤波器xyStstsgmnyxf,1,这个

9、操作可以用系数为 1/mn 的卷积模板来实现。n几何均值滤波器1,mns tSf x yg s t几何均值滤波器所达到的平滑度可以与算术均值滤波器相比，但在滤波过程中丢失更少的图像细节。图像复原 最新n谐波均值滤波器xyStstsgmnyxf,1,xyxyStsQStsQtsgtsgyxf,1,谐波均值滤波器对于“盐”噪声效果比较好，但是不适用于“胡椒”噪声。它善于处理像高斯噪声那样的其他噪声。n逆谐波均值滤波器当Q是正数时，用于消除“胡椒”噪声；当Q是负数时，用于消除“盐”噪声。但不能同时消除这两种噪声。Q的选择很重要！图像复原 最新5.3.2 顺序统计滤波器n中值滤波器n最大值和最小值滤波

10、器n中点滤波器n修正后的阿尔法均值滤波器图像复原 最新n中值滤波器tsgyxfxySts,max,tsgyxfxySts,min,用该像素相邻像素的灰度中值来代替该像素的值。适用于单级或双级脉冲噪声。n最大值滤波器这种滤波器在发现图像中的最亮点时非常有用。n最小值滤波器n这种滤波器在发现图像中的最暗点时非常有用。图像复原 最新n中点滤波器tsgtsgyxfxyxyStsSts,21,minmax,这种滤波器结合了顺序统计和求均匀，对于高斯和均匀随机分布噪声有最好的效果。n修正后的阿尔法均值滤波器图像复原 最新修正后的阿尔法均值滤波器n假设在邻域内去掉g(s,t)最高灰度值的d/2和最低灰度值的

11、d/2。由这些剩余后的像素点的平均值形成的滤波器n当d=0时，退变为算术均值滤波器；当 d=mn-1时，退变为中值滤波器；当d为其他值时，修正后的阿尔法均值滤波器在包括多种噪声的情况下非常适用。xyStsrtsgdmnyxf,1,图像复原 最新5.3.3 自适应滤波器n自适应滤波器考虑了图像中的像素点与其他像素点的特征的差异性，其行为变化基于由mn矩形窗口定义的区域内图像的统计特征。n自适应滤波器要优于迄今为止讨论过的所有滤波器的性能。但自适应滤波器的复杂度提高了。n两种简单的自适应滤波器n自适应、局部噪声消除滤波器n自适应中值滤波器图像复原 最新自适应、局部噪声消除滤波器n均值和方差是自适应

12、滤波器的基础22,nLLfx yg x yg x ym唯一需要知道或估计的量就是全部噪声的方差n2。其他参数需要从中各个坐标处的像素计算出来，在该位置滤波器窗口已被中心化。图像复原 最新自适应中值滤波器n自适应中值滤波器可以处理具有更大概率的冲激噪声。它的另一个优点是平滑非冲激噪声时可以保留细节。n工作在两层（参考书上）n主要目的n除去“椒盐”噪声n平滑其他非冲激噪声n并减少诸如物体边界细化或粗化等失真。图像复原 最新带阻滤波器n理想带阻滤波器nn阶巴特沃思带阻滤波器的表达式：2,2,22,1,0,1,0000WDvuDWDvuDWDWDvuDvuHnDvuDWvuDvuH2202,11,5.

13、3.4 频域滤波消减周期噪声图像复原 最新n高斯带阻滤波器表达式：n三个滤波器的透视图。2202,211,WvuDDvuDevuH带阻滤波器的主要应用之一是，在频率域噪声分量的一般位置近似已知的应用中消除噪声。图像复原 最新Butterworth带阻滤波器图像复原 最新带通滤波器n带通滤波器执行与带阻滤波器想反的操作。vuHvuHbrbp,1,n通常不在在一幅图像上直接执行带通滤波器，这通常会消除太多图像细节。n带通滤波器可以获取噪声模式，从而帮助屏蔽噪声模式。获取的噪声模式图像复原 最新陷波滤波器n陷波滤波器阻止（或通过）事先定义的中心频率邻域内的频率。图像复原 最新n理想陷波带阻滤波器其他

14、或,1,0,0201DvuDDvuDvuHnn阶的巴特沃思陷波带阻滤波器n高斯陷波带阻滤波器nvuDvuDDvuH,11,21202021,211,DvuDvuDevuH问题：当u0=v0=0时，以上三个滤波器都是什么滤波器？图像复原 最新n陷波滤波器vuHvuHnrnp,1,n最佳陷波滤波器n存在多种干扰最佳陷波滤波器n第一步屏蔽干扰的主要因素，第二步从被干扰的图像中减去一个可变的加权部分。n（具体内容见书）n小波分析方法更好！图像复原 最新第5讲 图像复原n5.1 图像退化/复原过程的模型n5.2 噪声模型n5.3 噪声存在下的图像复原n5.4 估计退化函数n5.5 图像复原的滤波方法n5

15、.6 图像复原的非线性方法n5.7 图像复原的盲反卷积方法n5.8 几何校正图像复原 最新退化退化H的性质的性质(1)线性：(2)相加性（k1=k2=1）：(3)一致性（f2(x,y)=0）：(4)位置（空间）不变性：),(),(),(),(22112211yxfHkyxfHkyxfkyxfkH),(),(),(),(2121yxfHyxfHyxfyxfH),(),(1111yxfHkyxfkH),(),(byaxgbyaxfH图像复原 最新点扩展函数的确定 大气湍流 光学系统散焦 照相机与景物相对运动 根据导致模糊的物理过程（先验知识）来确定h(x,y)或H(u,v)。图像复原 最新（1）长

16、时间曝光下大气湍流造成的传递函数5/622(,)expHu vkuvk是与湍流性质有关的常数。图像复原 最新图像复原 最新（2）光学散焦 d是散焦点扩展函数的直径,J1()是第一类贝塞尔函数。2/1221)()(),(vuddJvuH图像复原 最新图像复原 最新（3）照相机与景物相对运动 设T为快门时间，x0(t)，y0(t)是位移的x分量和y分量 dttvytuxjvuHT)()(2exp),(000图像复原 最新图像复原 最新 从退化图像本身来估计h(x,y)。（1）若有把握断定原始景物某部位有一个清晰的点，于是那个点再退回图像的模糊图像就是h(x,y)。（2）原景物含有明显的直线，从这些

17、线条的退化图像得出h(x,y)。（3）有明显的界限 可以证明：界线的退化图像的导数平行与该界线的线源的退化图像。图像复原 最新第5讲 图像复原n5.1 图像退化/复原过程的模型n5.2 噪声模型n5.3 噪声存在下的图像复原n5.4 估计退化函数n5.5 图像复原的滤波方法n5.6 图像复原的非线性方法n5.7 图像复原的盲反卷积方法n5.8 几何校正图像复原 最新5.5 图像复原的滤波方法n5.5.1 逆滤波n5.5.2 维纳滤波n5.5.3 等功率谱滤波n5.5.4 约束最小二乘方滤波图像复原 最新逆滤波 假定退化图像遵从以下模型),(/),(),(),(/),(),(1vuHvuGyxf

18、vuHvuGvuFF在不考虑噪声的情况下 则),(),(*),(),(yxnyxhyxfyxg),(),(),(vuHvuFvuG图像复原 最新函数。代表恢复滤波器的转移之逆，它是令),(,),(1),(vuHvuHvuP逆滤波图像复原 最新实际应用时的缺点n无噪声情况n若在频谱平面对图像信号有决定影响的点或区域上，H(u,v)的值为零，那么G(u,v)的值也为零，故不能确定这些频率处的F(u,v)值，也就难以恢复原始图像f(x,y)。n有噪声情况nG(u,v)=F(u,v)H(u,v)+N(u,v)n仍采用逆滤波器P(u,v)=1/H(u,v)nH(u,v)=0,没有定义。nH(u,v)=0

19、附近，H(u,v)较小，N(u,v)/H(u,v)会非常大，结果 与 F(u,v)大不相同，就不再象f(x,y),(vuF),(vuF(,)f x y图像复原 最新 在H(u,v）=0及其附近，人为地仔细设置H-1(u,v)的值，使N(u,v)*H-1(u,v)不会对F（u，v）产生太大影响。下图给出了H(u,v)、H-1(u,v)同改进的滤波特性HI(u,v)的一维波形，从中可看出与正常的滤波的差别。使H-1(u,v)具有低通滤波性质001DD 0DD ),(1),(vuHvuH改进的方法图像复原 最新逆滤波情况下的图像退化与复原记M(u,v)为复原传递函数，并不正好是1/H(u,v)图像退

20、化和复原模型图像退化和复原模型除去除去H(u,v)为零的点为零的点 减少振铃效应减少振铃效应k和d均为小于1的常数 fM()(x,yn()x,yH()u,vf()x,yg()x,yu,v20222022 1 ),(1),(wvuwvuvuHvuM如如其它如),(1),(),(vuHdvuHkvuM图像复原 最新模糊点源以获得转移函数模糊点源以获得转移函数将点源图像看做单位脉冲函数（F (x,y)=1）的近似则有G(u,v)=H(u,v)F(u,v)H(u,v)图像退化和恢复示例图像退化和恢复示例 退化图 滤波器 除去零点 减少振铃图像复原 最新匀速直线运动匀速直线运动 消除匀速直线运动模糊 )

21、,(),()d()(j2exp),(d dd )(j2)exp(),(-d)d(j2)exp,(),(000000vuFvuHttvytuxvuFtyxvyuxtyytxxfyxvyuxyxgvuGTT Tttyytxxfyxg000d )(),(),(T:采集时间长度采集时间长度x方向运动分量 y方向运动分量 图像复原 最新水平方向匀速直线运动水平方向匀速直线运动x0(t)=ct/T，y0(t)=0 当n为整数时，H在u=n/c处为零当 f(x,y)在区间0 x L之外为零或已知时 ucucucTtTctuvuHTjexp)sin(d j2exp)(0,LxftTctxfxgxcxT0)d(

22、d )(0图像复原 最新图像复原 最新5.5 图像复原的滤波方法n5.5.1 逆滤波n5.5.2 维纳滤波n5.5.3 等功率谱滤波n5.5.4 约束最小二乘方滤波图像复原 最新维纳滤波(,)(,)f x yf x y采用一准则使与之间的均方误差最小。求准则：),(),(),(:min22yxfyxfyxfEe图像复原 最新可推出 2222(,)(,)(,)(,)1(,)(,)(,)(,)(,)1(,)(,)fgnfS u vP u vS u vH u vS u vH u vH u vS u vH u vH u vH u vK图像复原 最新*11(,)(,)(2)(,)(,)(,)(,)(,)

23、fnP u vH u vH u vSu vP u vHu vSu v讨论：（）无噪声：。有噪声：在很小时，图像复原 最新图像复原 最新图像复原 最新5.5 图像复原的滤波方法n5.5.1 逆滤波n5.5.2 维纳滤波n5.5.3 等功率谱滤波n5.5.4 约束最小二乘方滤波图像复原 最新等功率谱滤波 方法：（1）退化模型 ),(),(),(),(yxnddyxhfyxg2(,)(,)(,)(,)(,)fnfSu vSu vS u vP u vH u v（）用、和表示。图像复原 最新等功率谱滤波 1/221/223(,)(,)(,)(,)(,)(,)(,)1(,)(,)(,)ffffnnfSu

24、vSu vSu vP u vSu v H u vS u vS u vH u vSu v（）由图像复原 最新1/21(,)0,(,)2(,)(,)(,)(,)nfnS u vP u vH u vSu vP u vS u v讨论：（）无噪声时：为逆滤波器。（）有噪声时：当很小或者零时，。图像复原 最新几何均值滤波器12*2*),(),(),(),(),(),(),(vuSvuSvuHvuHvuHvuHvuPfnn=1，逆滤波；n=0，参数维纳滤波；n=1，标准维纳滤波；n=1/2，普通逆滤波器和维纳滤波器的几何平均。n=1/2，=1，等功率谱均衡滤波器。图像复原 最新5.5 图像复原的滤波方法n5

25、.5.1 逆滤波n5.5.2 维纳滤波n5.5.3 等功率谱滤波n5.5.4 约束最小二乘方滤波图像复原 最新n约束最小二乘方滤波仅要求噪声方差和均值的知识。n滤波器特点n对于处理的每一幅图像都能产生最优结果n滤波器算法的选择往往由结果图像的感知视觉质量决定。n最优准则22min(,)f x yn约束条件22 nfHg图像复原 最新n 这个最佳化问题的频域解决方法),(),(),(),(*),(22vuGvuPvuHvuHvuF010(,)141010p x y图像复原 最新维纳滤波约束最小二乘方滤波图像复原 最新第5讲 图像复原n5.1 图像退化/复原过程的模型n5.2 噪声模型n5.3 噪

26、声存在下的图像复原n5.4 估计退化函数n5.5 图像复原的滤波方法n5.6 图像复原的非线性方法n5.7 图像复原的盲反卷积方法n5.8 几何校正图像复原 最新图像恢复的非线性方法 两种有代表性的非线性方法：（一）最大后验法：考虑了图像记录过程所具有的非线性并且把图像各点的集合平均看作是非平稳的随机过程。（二）最大熵法：考虑图像形成时的非线性，并且能保证图像函数有非负值。图像复原 最新最大后验法 若把图像记录的非线性考虑进去，图像退化模型为(,)(,)(,)(,)(,)(,)()g x yS b x yn x yb x yh x yfd dS gbnbHf|ln|m axpfgpfgff使

27、后 验 条 件 密 度 函 数取 最 大值 的 时 候，的 估 值 为，使。图像复原 最新最大后验法 NbnbTffbuuuSbuuuSbuuuSS)()()(211SHfgCSHCmf最终解得图像复原 最新最大熵方法 xyyxfE),(xyfyxfyxfH),(ln),(xyxyTTnByxnByxnyxnyxnH),(ln(),(),(ln),(假设图像函数具有非负值，定义一幅图像的总能量为图像的熵为噪声熵为最大熵复原就是在满足上式和图像退化模型的约束条件下使得复原后图像的图像熵和噪声熵最大。B为最小的噪声负值，以使上式中的对数有意义（令0ln0=0）图像复原 最新第5讲 图像复原n5.1

28、 图像退化/复原过程的模型n5.2 噪声模型n5.3 噪声存在下的图像复原n5.4 估计退化函数n5.5 图像复原的滤波方法n5.6 图像复原的非线性方法n5.7 图像复原的盲反卷积方法n5.8 几何校正图像复原 最新图像复原的盲反卷积方法n图像复原可以看作为一个反卷积的过程。n图像和退化系统先验知识非常有限甚至没 盲反卷积方法n盲反卷积的求解过程是相当困难的，具有病态性 具体问题具体对待n两大类方法n从退化图像中估计点扩散函数，然后采用合适的复原方法进行图像复原。（较容易实现）n在复原过程中同时估计点扩散函数和未降质图像（复杂）图像复原 最新迭代盲反卷积算法（IBD）n主要思想：利用真实图像

29、和PSF的某些先验知识，在随机选取真实图像的初始估计后，采用某种变换将图像在时域和频率域之间反复地变换，直到执行一定的迭代次数或算法最终收敛。基于FFT变换的IBD算法结构图像复原 最新图像复原与考古图像复原 最新第5讲 图像复原n5.1 图像退化/复原过程的模型n5.2 噪声模型n5.3 噪声存在下的图像复原n5.4 估计退化函数n5.5 图像复原的滤波方法n5.6 图像复原的非线性方法n5.7 图像复原的盲反卷积方法n5.8 几何校正图像复原 最新几何失真n几何失真原因n图像在获取过程中，由于成像系统本身具有非线性、拍摄角度等因素的影响，会使获得的图像产生几何失真。n几何失真n系统失真是有

30、规律的、能预测的；n非系统失真则是随机的。n几何校正的必要性n当对图像作定量分析时，就要对失真的图像先进行精确的几何校正，以免影响定量分析的精度。图像复原 最新几何校正方法n基本方法：先建立几何校正的数学模型；其次利用已知条件确定模型参数；最后根据模型对图像进行几何校正。n图像空间坐标变换n确定各像素的灰度值（灰度内插）。图像复原 最新设图像f(x,y)是经过几何失真产生图像g(x,y)表示，则变换可以表示为：),(1yxhx),(2yxhy 空间坐标变换空间坐标变换通常h1(x，y)和h2(x，y)可用多项式来近似 niinjjiijyxax00niinjjiijyxby00图像复原 最新当

31、n=2时，畸变关系式为20211220011000yaxyaxayaxaax20211220011000ybxybxbybxbby当n=1时，畸变关系为线性变换，yaxaax011000ybxbby011000问题：求解上式需要已知几个像素点？问题：求解上式需要已知几个像素点？图像复原 最新像素灰度内插方法像素灰度内插方法n作用：处理空间变换后图像中像素灰度级的赋值n常用的像素灰度内插法n最近邻元法n双线性内插法n三次内插法图像复原 最新最近邻元法n在待求点的四邻像素中，将距离这点最近的相邻像素灰度赋给该待求点。n方法最简单，效果尚佳，但校正后的图像有明显锯齿状，即存在灰度不连续性。图像复原

32、最新对于(i，j+v)有f(i，j+v)=f(i，j+1)-f(i，j)v +f(i，j)对于(i+1，j+v)有f(i+1，j+v)=f(i+1，j+1)-f(i+1，j)v+f(i+1，j)双线性内插法 n双线性内插法是利用待求点四个邻像素的灰度在两个方向上作线性内插。图像复原 最新 对于(i+u,j+v)有f(i+u,j+v)=f(i+1,j+v)-f(i,j+v)u+f(i,j+v)=)1,1(),1()1()1,()1(),()1)(1(jiuvfjifvujivfujifvu该方法要比最近邻元法复杂，计算量大。但没有灰度不连续性的缺点，结果令人满意。它具有低通滤波性质，使高频分量受

33、损，图像轮廓有一定模糊。图像复原 最新2|02|1|5|841|0|21)(3232xxxxxxxxxS(i-1,j-1)(i-1,j+2)(i+2,j-1)(i+2,j+2)(x,y)u v三次内插法n该方法利用三次多项式S(x)来逼近理论上的最佳插值函数sin(x)/x。其数学表达式为图像复原 最新其中A=s(1+v)s(v)s(1-v)s(2-v)(1,1)(1,)(1,1)(1,2)(,1)(,)(,1)(,2)(1,1)(1,)(1,1)(1,2)(2,1)(2,)(2,1)(2,2)f ijf ijf ijf ijf i jf i jf i jf i jBf ijf ijf ijf ijf ijf ijf ijf ijc=s(1+u)s(u)s(1-u)s(2-u)T该算法计算量最大，但内插效果最好，精度最高。待求像素(x,y)的灰度值由其周围十六个点的灰度值加权内插得到。可推导出待求像素的灰度计算式如下：f(x,y)=A B C图像复原 最新原始影像灰度表面 最近邻内插法双线性内插法 三次内插法像素灰度内插法效果比较