photo.m

MATLAB 基础编程 包含一些很基础的程序

function [f2,f3]=photo(ut)
%本例只考虑有一行横向贯穿破损，接下来思考存在
%多行横向贯穿破损、一行非贯穿破损以及不定方向破损的情况
f1=imread('lena256.bmp');
[m,n]=size(f1);

f2=f1;
f2(ut,:)=255;              %f2为破损图像

figure;
subplot(1,3,1);
imshow(f1);
xlabel('原始图像');
subplot(1,3,2);
imshow(f2);
xlabel('破损图像');
%考虑只有一行无数据的情况

f3=double(f2)+1;    %将uint8类型转换为双精度浮点类型
%经验公式，否则修补后f3矩阵与f1矩阵不吻合
ln=0;
mn=0;
ss=zeros(1,n);
st=0;
%mn=max(mean(f3'));   %找到每行平均值的最大值
for i=1:m
     ss=f3(i,:);
     st=sum(ss);
    if st/n>=max(max(f1))
    ln=i;
   end
   if ln~=0,break,end
   st=0;
end

%利用像素相关填补空白
if ln==1
   f3(ln,:)=f3(ln+1,:);
elseif ln==m
   f3(ln,:)=f3(ln-1,:);
else
   f3(ln,:)=(f3(ln-1,:)+f3(ln+1,:))./2; 
end

f3=uint8(f3-1);     %将双精度浮点类型转换为uint8类型
subplot(1,3,3);
imshow(f3);
xlabel('修补后图像');
imwrite(f2,'lena2.bmp');
imwrite(f3,'lena3.bmp');
error=abs(max(double(f1(ln,:))-double(f3(ln,:))));
bbt=['对应最大误差为',num2str(error)]
sop=0;