close?all;
clc;
clear;
%--------------------------讀入圖像------------------------------------------
imgoriginal=imread（‘yuan.png‘）;??
figure（1）;
imshow（imgoriginal）;?
[widthheight]?=?size（imgoriginal）;?
img=?double（imgoriginal）;?
%--------------------------獲取掩碼圖像---------------------------------------------
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;?
????????????mask（ij2）?=?255;?
????????????mask（ij3）?=?255;?
????????else?
????????????mask（ij1）?=?0;?
????????????mask（ij2）?=?0;
????????????mask（ij3）?=?0;?
????????end?
????end?
end?
figure（2）;
imshow（mask）;
%--------------------------更新原圖修復區域內每點值------------------------------------------
a=zeros（widthheight）;
I=cat（3a2*a3*a）;
J=cat（3a2*a3*a）;
n?=?1;?
itertimes=1000;?
tic;
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）;
????????????????Kn（k）=sqrt（（I（ij+1k）-I（ijk）+（I（i+1jk）+I（i+1j+1k）-I（i-1jk）-I（i-1j+1k））/2）^2+（J（ij+1k）-J（ijk）+（I（i+1jk）+J（i+1j+1k）-J（i-1jk）-J（i-1j+1k））/2）^2）;
????????????????w1（k）?=?Kw（k）/sqrt（1+grid_w（k））+1;
????????????????w2（k）?=?Ke（k）/sqrt（1+grid_e（k））+1;
????????????????w3（k）?=?Ks（k）/sqrt（1+grid_s（k））+1;
????????????????w4（k）?=?Kn（k）/sqrt（1+grid_n（k））+1;
????????????????img（ijk）?=（w1（k）*img（i-1jk）+w2（k）*img（i+1jk）+w3（k）*img（ij-1k）+w4（k）*img（ij+1k））/（w1（k）+w2（k）+w3（k）+w4（k））;
???

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????66614??2012-04-14?18:39??CDD\bc.bmp

?????文件?????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

?????文件??????66614??2012-11-06?23:20??CDD\cdd.bmp

?????文件???????3281??2016-10-27?15:02??CDD\CDD.m

?????文件??????66614??2012-10-04?20:24??CDD\cdd2_abcd.bmp

?????文件?????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

?????文件?????504456??2016-10-27?15:05??CDD\CDD修復后.bmp

?????文件?????786486??2012-04-20?22:03??CDD\dt.bmp

?????文件?????786486??2012-04-20?22:03??CDD\dtx.bmp

?????文件??????66614??2012-10-07?20:44??CDD\TV1.bmp

?????文件?????126849??2016-10-27?13:52??CDD\yuan.png

?????文件?????786486??2012-03-18?22:16??TV\lena.bmp

?????文件????1710666??2012-06-06?16:42??TV\lena1.bmp

?????文件?????786486??2012-06-06?16:42??TV\lena_修復后.bmp

?????文件?????786486??2012-06-06?15:05??TV\lena_原圖.bmp

?????文件?????786486??2012-06-06?16:42??TV\lena_掩碼圖像.bmp

?????文件???????2189??2016-10-27?15:10??TV\TV.m

?????文件????2270846??2012-06-06?16:47??TV\男兵.bmp

?????文件????1044989??2012-06-06?16:47??TV\男兵_修復后.png

?????文件?????136150??2012-06-06?15:24??TV\男兵_原圖.jpg

?????文件??????35655??2012-06-06?16:43??TV\男兵_掩碼圖像.png

............此處省略5個文件信息