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
????%更新修復(fù)區(qū)域內(nèi)每點(diǎn)值
????for?i?=?2:width-1
????????for?j?=?2:height-1
????????????

????????????if?MASK（ij）?==?0
??????????????
????????????????%計(jì)算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）;
??????????????
????????????????%計(jì)算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

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????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——網(wǎng)格.jpg
?????文件???????64636??2011-11-15?18:46??CDD程序\7——網(wǎng)格.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