clear?all
close?all

I=imread（‘color.bmp‘）;
B=imread（‘depth.bmp‘）;
B=rgb2gray（B）;
[mnhh]=size（I）;
A=rgb2gray（I）;%彩色圖轉換灰度圖
I1=I;
figure（1）imshow（I）;title（‘原始受污染的圖像‘）;
figure（2）imshow（B）;title（‘深度圖‘）;
figure（3）imshow（A）;title（‘灰度圖‘）;

[rowindcolumnind]=find（A<4）;%灰度圖中小于4灰度值的認為是受到污染的像素統計受污染像素的位置，rowind和columnind均是向量，并且兩者一一對應組成位置坐標
pointnum=length（rowind）;%受污染像素的個數

rowstart=rowind-20;rowend=rowind+20;%統計每個受污染像素修復模板的起止行
rowstart（rowstart<1）=1;rowend（rowend>m）=m;%起始行不能小于1，終止行不能大于m

columnstart=columnind-20;columnend=columnind+20;%統計模板的起，止列
columnstart（columnstart<1）=1;columnend（columnend>n）=n;%以基準點（受污染點）為中心41*41模板

for?num=1:pointnum%對每個受污染的像素進行修復
????x=rowind（num）;
????y=columnind（num）;
????basedepth=B（xy）;%深度圖基準點灰度值
????
????xstart=rowstart（num）;
????xend=r

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

?????文件????2359350??2011-10-21?08:00??深度圖像修補\color.bmp

?????文件????2359350??2011-10-21?08:00??深度圖像修補\depth.bmp

?????文件???????1903??2012-06-12?22:51??深度圖像修補\impainting.m

?????目錄??????????0??2012-06-12?22:50??深度圖像修補

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

??????????????4720603????????????????????4