图像修复cdd程序

图像修复 基于matlab的cdd技术修复代码 附有修复图片和详细过程

function [B] = CDDInpainting（）
tic;
%t0 = clock;
ImgOriginal=imread（‘F:\wall\CDD程序\Lena1.bmp‘）;
%ImgOriginal=rgb2gray（ImgOriginal）;
figure（1）imshow（ImgOriginal）;
Img = double（ImgOriginal）;
[widthheight]=size（Img）;
U = Img;
V = Img;
n = 1;
a=0.0001;
IterTimes=150;
MASK=~（Img>=251）; %人为选择阈值
figure（2）;imshow（MASK）;



while n <= IterTimes
    %更新修复区域内每点值
    for i = 2:width-1
        for j = 2:height-1
            

            if MASK（ij） == 0
              
                %计算w1w2w3w4
                gridw2 = （（V（ij）-V（i-1j））^2+（（V（i-1j-1）-V（i-1j+1））/2）^2）^0.25;
                gride2 = （（V（ij）-V（i+1j））^2+（（V（i+1j-1）-V（i+1j+1））/2）^2）^0.25;
                grids2 = （（V（ij）-V（ij-1））^2+（（V（i-1j-1）-V（i+1j-1））/2）^2）^0.25;
                gridn2 = （（V（ij）-V（ij+1））^2+（（V（i-1j+1）-V（i+1j+1））/2）^2）^0.25;
                a1 = 1/sqrt（gridw2+a）;
                a2 = 1/sqrt（gride2+a）;
                a3 = 1/sqrt（grids2+a）;
                a4 = 1/sqrt（gridn2+a）;
          k = （a1*V（ij+1）+a2*V（ij-1）+a3*V（i-1j）+a4*V（i+1j））/（a1+a2+a3+a4）;
              
                %计算w1w2w3w4
                gridUw2 = （V（ij）-V（i-1j））^2+（（V（i-1j-1）-V（i-1j+1））/2）^2;
                gridUe2 = （V（ij）-V（i+1j））^2+（（V（i+1j-1）-V（i+1j+1））/2）^2;
                gridUs2 = （V（ij）-V（ij-1））^2+（（V（i-1j-1）-V（i+1j-1））/2）^2;
                gridUn2 = （V（ij）-V（ij+1））^2+（（V（i-1j+1）-V（i+1j+1））/2）^2;
                w1 = k/sqrt（gridUw2+a）;
                w2 = k/sqrt（gridUe2+a）;
                w3 = k/sqrt（gridUs2+a）;
                w4 = k/sqrt（gridUn2+a）;
               U（ij） =（w1*V（ij+1）+w2*V（ij-1）+w3*V（i-1j）+w4*V（i+1j））/（w1+w2+w3+w4）;
            end
        end
    end
    n = n+1;
    V = U;
end
D = floor（V）;
figure（3）;imshow（D[]）;
toc

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2012-10-12 16:12  CDD程序\
     文件       68957  2011-11-15 18:14  CDD程序\1——建大图像.jpg
     文件       65478  2011-11-15 18:11  CDD程序\2——建大图像.jpg
     文件       67743  2011-11-15 18:40  CDD程序\3——去除.jpg
     文件       68122  2011-11-15 18:39  CDD程序\4——去除.jpg
     文件       66155  2011-11-15 18:30  CDD程序\5 tv1——建大图像.jpg
     文件       69006  2011-11-15 18:47  CDD程序\6——网格.jpg
     文件       64636  2011-11-15 18:46  CDD程序\7——网格.jpg
     文件        1853  2012-05-20 17:43  CDD程序\CDDgrayimproved.asv
     文件        1852  2012-05-02 11:19  CDD程序\CDDgrayimproved.m
     文件      263222  2010-11-22 22:01  CDD程序\Lena1.bmp
     文件      263222  2010-11-22 22:01  CDD程序\Lena6.bmp
     文件         766  2011-11-22 11:47  CDD程序\MATLAB 7.0.1.lnk
     文件        1197  2012-05-22 15:00  CDD程序\tv.asv
     文件        1202  2012-05-02 11:18  CDD程序\tv.m
     文件      263222  2011-11-14 17:38  CDD程序\原图.jpg