clc
clear?all;
close?all;
OriginalImg?=?rgb2gray（imread（‘zhiwen.jpg‘））;
figure（1）;imshow（OriginalImg）;title（‘原圖像‘）;
DirtyImg=rgb2gray（imread（‘zhiwen1.jpg‘））;
figure（2）;imshow（DirtyImg）;title（‘破損后圖像‘）;
OriginalImg=double（OriginalImg）;
DirtyImg=double（DirtyImg）;
%獲取圖像寬度和高度
[widthheight]?=?size（DirtyImg）;

Img=DirtyImg;

U?=?Img;
V?=?Img;
IterTimes=100;
%利用邏輯運算提取掩模
%?MASK=~（Img）;?%人為選擇閾值s
MASK=zeros（widthheight）;
for?j=1:（height-1）
???for??i=1:（width-1）
?????if?Img（ij）<15
?????????MASK（ij）=0;
?????else
?????????MASK（ij）=255;
?????end
???end
end
%?MASK?=?（Img?>?251）;%lena;集市圖像
MASK=floor（MASK）
imwrite（MASK‘mask1.bmp‘）;
figure（3）;imshow（‘mask1.bmp‘）;
title（‘受損區域‘）;
%計算受損率



%?%將二值掩模中間的空洞填充，使之成為純粹的二值掩模
%?noborder=imclearborder（MASK4）;
%?%?figureim

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????153078??2012-05-07?09:03??畢設\inpainted.bmp

?????文件?????153078??2012-05-07?09:03??畢設\mask1.bmp

?????文件??????13216??2012-05-05?09:04??畢設\sheying.JPG

?????文件??????13597??2012-05-05?09:00??畢設\sheying1.JPG

?????文件???????2971??2012-05-05?09:14??畢設\tv.m

?????文件??????23836??2012-05-05?08:59??畢設\zhiwen.jpg

?????文件??????15246??2012-05-05?08:35??畢設\zhiwen1.jpg

?????目錄??????????0??2012-05-07?09:37??畢設

-----------?---------??----------?-----??----

???????????????375022????????????????????8