matlab图像修复

基于距离向量矩阵，对有有污点的图像进行修复。

clear all
close all
I=imread（‘1.jpg‘）;
[mnhh]=size（I）;
A=rgb2gray（I）;
I1=I;
[rowindcolumnind]=find（A<15）;
pointnum=length（rowind）;
rowstart=rowind-20;rowend=rowind+20;
rowstart（rowstart<1）=1;rowend（rowend>m）=m;
columnstart=columnind-20;columnend=columnind+20;
columnstart（columnstart<1）=1;columnend（columnend>n）=n;

for num=1:pointnum
    x=rowind（num）;
    y=columnind（num）;   
    xstart=rowstart（num）;
    xend=rowend（num）;
    x1=x-xstart+1;
    ystart=columnstart（num）;
    yend=columnend（num）;
    y1=y-ystart+1;
    
    depthpos=A（xstart:xendystart:yend）>15;
    [locind1locind2]=find（depthpos）;
    W=（locind1-x1）.*（locind1-x1）+（locind2-y1）.*（locind2-y1）+1;
    
    TT=1./W;T=sum（TT）; 
    Aloc=I（xstart:xendystart:yend1）;
    Aloc=double（Aloc（depthpos））;
    R=Aloc.*TT;R=sum（R）;
    I1（xy1）=R/（T+0.1）;
    
    Aloc=I（xstart:xendystart:yend2）;
    Aloc=double（Aloc（depthpos））;
    R=Aloc.*TT;R=sum（R）;
    I1（xy2）=R/（T+0.1）;
    
    Aloc=I（xstart:xendystart:yend3）;
    Aloc=double（Aloc（depthpos））;
    R=Aloc.*TT;R=sum（R）;
    I1（xy3）=R/（T+0.1）;
end
figureimshow（I[]）;
title（‘污染图‘）;
figureimshow（I1[]）;
title（‘去污图‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1209  2012-05-31 12:29  xiufu.m

     文件     167213  2012-05-31 12:26  1.jpg

----------- ---------  ---------- -----  ----

               168422                    2