clear;clc;close?all
%處理第一張圖片（對(duì)應(yīng)距離為D）
I=imread（‘20180409-204942-278.jpg‘）;?????????%讀入圖片
figure
subplot（231）;
imshow（I）;
title（‘原圖‘）;???
I=rgb2gray（I）;?????????????%將RGB圖轉(zhuǎn)換成灰度圖
subplot（232）;
imshow（I）;
title（‘灰度圖‘）;
[mn]=size（I）;
for?K=1:10??????????????
????I2=wiener2（I[3?3]）;???%濾波10次，把多余的點(diǎn)濾去。
end
I=edge（I‘sobel‘）;?????????%集中邊緣檢測(cè)算子sobel、canny、log、roberts、prewitt。在此處sobel算子效果最好
subplot（233）;
imshow（I）;
title（‘邊緣提取圖‘）;
for?i=1:m??????????????????%y??????
????for?j=1:n??????????????%x
????????if???i<280?||?i>340||?j<430||?j>490?
???????????I（ij）=0;
????????end
????end
end
subplot（234）;
imshow（I）;
title（‘目標(biāo)物提取圖‘）;
???????????????????????????%形態(tài)學(xué)處理，填充輪廓內(nèi)像素
SE=strel（‘diamond‘2）;?????%?strel在操作結(jié)構(gòu)元素應(yīng)用，用于膨脹腐蝕及開(kāi)閉運(yùn)算等操作的結(jié)構(gòu)元素對(duì)象
???????????????????????????%具體用法：SE?=?strel（shapeparameters）
???????????????????????????%創(chuàng)建由指定形狀shape對(duì)應(yīng)的結(jié)構(gòu)元素。%參數(shù)parameters一般控制SE的大小。
I2=imdilate（ISE）;?????????%imdilate對(duì)圖像實(shí)現(xiàn)膨脹操作。
???????????????????????????%IM2?=?imdilate（IMSE）膨脹灰度，二值，壓縮二值圖像IM，返回IM2。
???????????????????????????%參數(shù)SE為由strel函數(shù)返回的結(jié)構(gòu)元素或者結(jié)構(gòu)元素對(duì)象組。
I2=imfill（I2‘holes‘）;?????%該函數(shù)用于填充二值圖像中的空洞區(qū)域。?如，?黑色的背景上有個(gè)白色的圓圈。?則這個(gè)圓圈內(nèi)區(qū)域?qū)⒈惶畛洹?br/>subplot（235）;
imshow（I2）;
title（‘填充圖‘）;
sum1=0;
for?i=1:m??????????????????%計(jì)算白像素的個(gè)數(shù)，即為圖形面積
????for?j=1:n
????????if?I（ij）==1;
????????????sum1=sum1+1;
????????end
????end
end
fprintf（‘在距離為D時(shí)目標(biāo)物的圖形面積（像素點(diǎn)個(gè)數(shù)）s1=%f\n‘sum1）;


%處理第二張圖片（對(duì)應(yīng)距離為D+d）
I=imread（‘20180409-205011-792.jpg‘）;????%讀入圖片
figure
subplot（231）;
imshow（I）;
title（‘原圖‘）;
I=rgb2gray（I）;??????