多媒体技术与应用之图像修复实验

普通的图像修补代码：

clear all 
close all 
flag = 20;%设置模板大小为41*41
I=imread('color.bmp'); 
subplot(121);imshow(I);title('原始图像 ') ;

[m,n,hh]=size(I); %长 宽 三层rgb
A=rgb2gray(I);%彩色图转换灰度图 
I1=I; 
[rowind,columnind]=find(A<4);
pointnum=length(rowind);%受污染像素的个数

rowstart=rowind-flag;
rowend=rowind+flag;%统计每个受污染像素修复模板的起止行 

rowstart(rowstart<1)=1;
rowend(rowend>m)=m;%起始行不能小于1，终止行不能大于m 
columnstart=columnind-flag;
columnend=columnind+flag;%统计模板的起止列 
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:xend,ystart:yend)>10 ; %参与加权的像素与基准点相差大于10的深度图像素点且不能是受污染点 
	
    
    [locind1,locind2]=find(depthpos); 
	W=(locind1-x1).^2+(locind2-y1).^2;% 权值为距离平方 的倒数 
	TT=1./W;
	T=sum(TT); 
    
	Aloc=I(xstart:xend,ystart:yend,1);%对红色通道进行处理 
	Aloc=double(Aloc(depthpos)); 
	R=Aloc.*TT;
	R=sum(R); 
	I1(x,y,1)=R/T;
    
	Aloc=I(xstart:xend,ystart:yend,2);%对绿色通道进行处理 
	Aloc=double(Aloc(depthpos)); 
	R=Aloc.*TT;
	R=sum(R); 
	I1(x,y,2)=R/T; 
    
	Aloc=I(xstart:xend,ystart:yend,3);%对蓝色通道进行处理 
	Aloc=double(Aloc(depthpos)); 
	R=Aloc.*TT;
	R=sum(R); 
    
	I1(x,y,3)=R/T; 
end 
subplot(122),imshow(I1),title(' 基于单幅图像修复 ')