close?all;
clc;
clear;
%--------------------------讀入圖像------------------------------------------
imgoriginal=imread（‘c1003.bmp‘）;?????????????????%待處理的圖片轉(zhuǎn)化為矩陣，黑白2維矩陣，彩色三維，例a=imread（xxx.jpg）R=a（::1）;?G=a（::2）;?B=a（::3）;
figure（1）;???????????????????????????????????????%figure創(chuàng)建圖像，imshow顯示圖片，一起為避免覆蓋
imshow（imgoriginal）;?????????????????????????????%顯示待處理的圖片
[widthheight]?=?size（imgoriginal）;??????????????%size（），獲取圖片矩陣的行數(shù)為width、列數(shù)為height
img=?double（imgoriginal）;????????????????????????%轉(zhuǎn)換成double精度
%--------------------------獲取掩碼圖像（mask遮擋物）---------------------------------------------
mask?=?zeros（widthheight/33）;?
for?j?=?1:height/3
????for?i?=?1:width
????????if?（（imgoriginal（ij1）?>220）&&（imgoriginal（ij2）?>220）&&（imgoriginal（ij3）?>220））?
????????????mask（ij1）?=?255;???????????????????%灰度255=白?，0=黑
????????????mask（ij2）?=?255;?
????????????mask（ij3）?=?255;?
????????else?
????????????mask（ij1）?=?0;?
????????????mask（ij2）?=?0;
????????????mask（ij3）?=?0;?
????????end?
????end?
end?
figure（2）;
imshow（mask）;????????????????????????????????????%顯示遮擋物
%--------------------------更新原圖修復(fù)區(qū)域內(nèi)每點(diǎn)值------------------------------------------
a=zeros（widthheight）;???????????????????????????%生成全0矩陣
I=cat（3a2*a3*a）;
J=cat（3a2*a3*a）;
n?=?1;?
itertimes=1000;??????????????????????????????????%迭代時(shí)間?iteration迭代
tic;?????????????????????????????????????????????%tic與toc，記錄命令執(zhí)行的時(shí)間
while?n?<=?itertimes
????for?i?=?2:width-1
????????for?j?=?2:height/3-1
???????????if?（mask（ij+11）?==?255）||（mask（ij-11）?==?255）||（mask（i+1j1）?==?255）||（mask（i-1j1）?==?255）
????????????????for?k=1:3?
????????????????grid_w（k）?=?（img（ijk）-img（i-1jk））^2+（1.0/16）*（img（i-1j+1k）+img（ij+1k）-img（i-1j-1k）-img（ij-1k））^2;
????????????????grid_e（k）?=?（img（ijk）-img（i+1jk））^2+（1.0/16）*（img（ij+1k）+img（i+1j+1k）-img（ij-1k）-img（i+1j-1k））^2;
????????????????grid_s（k）?=?（img（ijk）-img（ij-1k））^2+（1.0/16）*（img（i+1jk）+img（i+1j-1k）-img（i-1jk）-img（i-1j-1k））^2;
????????????????grid_n（k）?=?（img（ijk）-img（ij+1k））^2+（1.0/16）*（img（i+1jk）+img（i+1j+1k）-img（i-1jk）-img（i-1j+1k））^2;
????????????????I（ijk）=0.5*（img（i+1jk）-img（i-1jk））/sqrt（0.25*（img（i+1jk）-img（i-1jk））^2+0.25*（img（ij+1k）-img（ij-1k））^2+1）;
????????????????J（ijk）=0.5*（img（ij+1k）-img（ij-1k））/sqrt（0.25*（img（i+1jk）-img（i-1jk））^2+0.25*（img（ij+1k）-img（ij-1k））^2+1）;
????????????????Kw（k）=sqrt（（I（ijk）-I（i-1jk）+（I（i-1j+1k）+I（ij+1k）-I（i-1j-1k）-I（ij-1k））/2）^2+（J（ijk）-J（i-1jk）+（J（i-1j+1k）+J（ij+1k）-J（i-1j-1k）-J（ij-1k））/2）^2）;
????????????????Ke（k）=sqrt（（I（i+1jk）-I（ijk）+（I（ij+1k）+I（i+1j+1k）-I（ij-1k）-I（i+1j-1k））/2）^2+（J（i+1jk）-J（ijk）+（J（ij+1k）+J（i+1j+1k）-J（ij-1k）-J（i+1j-1k））/2）^2）;
????????????????Ks（k）=sqrt（（I（ijk）-I（ij-1k）+（I（i+1jk）+I（i+1j-1k）-I（i-1jk）-I（i-1j-1k））/2）^2+（J（ijk）-J（ij-1k）+（J（i+1jk）+J（i+1j-1k）-J（i-1jk）-J（i-1j-1k））/2）^2）;
????????????????

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-10-21?17:37??CDD\
?????文件??????196662??2012-10-04?21:21??CDD\C1.bmp
?????文件??????196662??2012-10-04?22:35??CDD\C1001.bmp
?????文件??????196662??2012-10-04?23:30??CDD\C1003.bmp
?????文件??????196662??2012-10-04?23:34??CDD\C1004.bmp
?????文件??????196662??2012-04-20?20:58??CDD\C13.bmp
?????文件??????196662??2012-04-20?20:55??CDD\C13x.bmp
?????文件????????3960??2018-10-21?17:36??CDD\CDD.m
?????文件??????196662??2012-05-01?10:26??CDD\CDD_n100.bmp
?????文件??????196662??2012-10-04?23:32??CDD\CDD_n1000.bmp
?????文件??????196662??2012-10-04?23:36??CDD\CDD_n10004.bmp
?????文件??????196662??2012-05-01?10:27??CDD\CDD_n200.bmp
?????文件??????196662??2012-05-01?10:23??CDD\CDD_n50.bmp
?????文件??????196662??2018-10-21?15:28??CDD\CDD修復(fù)后.bmp
?????文件???????66614??2012-10-07?20:44??CDD\TV1.bmp
?????文件???????66614??2012-04-14?18:39??CDD\bc.bmp
?????文件???????66614??2012-11-06?23:20??CDD\cdd.bmp
?????文件???????66614??2012-10-04?20:24??CDD\cdd2_abcd.bmp
?????文件??????786486??2012-04-20?22:03??CDD\dt.bmp
?????文件??????786486??2012-04-20?22:03??CDD\dtx.bmp
?????文件??????126849??2016-10-27?13:52??CDD\yuan.png
?????文件??????113992??2018-10-05?22:04??CDD\圖片2.png