image-inpainting 图像修复

图像修复，使用曲率方法，热传递。边缘修复非常好，参考某博士论文写。

clc
clear
R=imread（‘L4_6.jpg‘）; %读入待修复图像
R=double（R）;        %把图像值uint8型转为double型
[nrow ncol]=size（R）;  %读图像大小，nrow为图像高度，ncol为图像宽度
% R_mask=imread（‘mask.jpg‘）;%读模版
mask=double（R>80 & R<150）;  %把模版转换成0、1阵列，待修复为1，其他为0=1;
dt=0.2;%时间步长
I=R;
tic   %计时开始
for n=1:2000  %迭代次数，其值可以任选
    for i=3:nrow-2
        for j=3:ncol-2  %选择开始点和终点，以消除在图像边缘，进行中心差分法            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