CDD

才用CDD对灰度或者彩色图像进行修复的一种方法

function [B] = CDDInpainting（）
tic;
%t0 = clock;
ImgOriginal=imread（‘Lena1.bmp‘）;
figure（1）;imshow（ImgOriginal）;
[widthheight] = size（ImgOriginal）;
Img = double（ImgOriginal）;
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+1） == 255）|（MASK（ij-1） == 255）|（MASK（i+1j） == 255）|（MASK（i-1j） == 255）
            if MASK（ij） == 0
                if MASK（ij） == 0
                %计算w1w2w3w4
                gridw2 = （V（ij）-V（i-1j））^2+（（V（i-1j-1）-V（i-1j+1））/2）^2;
                gride2 = （V（ij）-V（i+1j））^2+（（V（i+1j-1）-V（i+1j+1））/2）^2;
                grids2 = （V（ij）-V（ij-1））^2+（

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2369  2010-12-17 13:37  CDD\CDDgray.m

     文件       5723  2010-12-17 13:37  CDD\CDDRGB.m

     文件     263224  2010-12-17 13:37  CDD\Lena.bmp

     文件     263222  2010-12-17 13:37  CDD\Lena1.bmp

     文件     263222  2010-12-17 13:37  CDD\Lena2.bmp

     文件     263222  2010-12-17 13:37  CDD\Lena3.bmp

     文件     263222  2010-12-17 13:37  CDD\Lena4.bmp

     文件     263222  2010-12-17 13:37  CDD\Lena5.bmp

     文件     263222  2010-12-17 13:37  CDD\Lena6.bmp

     文件     274262  2010-12-17 13:37  CDD\Lincoln2.bmp

    ..A.SH.     42496  2010-12-17 14:04  CDD\Thumbs.db

     目录          0  2010-12-17 14:04  CDD

----------- ---------  ---------- -----  ----

              2167406                    12