clc
clear
R=imread（‘L4_6.jpg‘）;?%讀入待修復(fù)圖像
R=double（R）;????????%把圖像值uint8型轉(zhuǎn)為double型
[nrow?ncol]=size（R）;??%讀圖像大小，nrow為圖像高度，ncol為圖像寬度
%?R_mask=imread（‘mask.jpg‘）;%讀模版
mask=double（R>80?&?R<150）;??%把模版轉(zhuǎn)換成0、1陣列，待修復(fù)為1，其他為0=1;
dt=0.2;%時間步長
I=R;
tic???%計時開始
for?n=1:2000??%迭代次數(shù)，其值可以任選
????for?i=3:nrow-2
????????for?j=3:ncol-2??%選擇開始點和終點，以消除在圖像邊緣，進(jìn)行中心差分法????????????if?mask（ij）==1
????????????%求J_e^1
????????????deta_I1=（（I（ij+1）-I（ij））^2+（0.25*（I（i-1j+1）+I（i-1j）-I（i+1j+1）-I（i+1j）））^2+eps）^（1/2）;?
????????????k1_x=（I（ij+2）+I（ij）-2*I（ij+1））/（（I（ij+2）+I（ij）-2*I（ij+1））^2+（I（i-1j+1）+I（i+1j+1）-2*I（ij+1））^2+eps）^（1/2）...
????????????????-（（I（ij+1）+I（ij-1）-2*I（ij））/（（I（ij+1）+I（ij-1）-2*I（ij））^2+（I（i-1j）+I（i+1j）-2*I（ij））^2+eps）^（1/2））;
????????????k1_y=0.25*（（I（i-2j+1

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2011-11-25?15:23??image?inpainting\
?????文件??????229888??2011-11-03?16:29??image?inpainting\CDD.doc
?????文件????????4972??2011-11-04?10:29??image?inpainting\CDD.m
?????文件????????1313??2010-04-05?20:37??image?inpainting\L4_6.jpg
?????目錄???????????0??2011-11-25?15:23??image?inpainting\RRR\
?????文件???????21078??2011-11-03?16:10??image?inpainting\RRR\0.bmp