F=imread('e:\毕业论文\分水岭\汽车.jpg');
I=imread('e:\毕业论文\分水岭\背景.jpg');
bw=im2bw(F,graythresh(F)); % 转换为黑白二值图像
bwc=~bw; % 图像反色
dst=bwdist(bwc); % 图像距离
ws=watershed(-dst);
w=ws==0;
rf=bw&~w;
figure(1),imshow(F);title('原图像'); % 显示原图像
figure(2),subplot(1,2,1),imshow(F1);title('灰度图'); % 显示灰度图像
figure(2),subplot(1,2,2),imhist(F1);title('灰度图直方图'); % 显示灰度直方图图像
figure(3),imshow(F1);title('运动车辆灰度图'); % 显示处理后的图像
figure(4),imshow(I1);title('背景灰度图'); % 显示处理后的图像
h=fspecial('sobel'); title('滤波器'); % 设置滤波器
fd=im2double(F); % 数据类型转换
sq=sqrt(imfilter(fd,h,'replicate').^2+imfilter(fd,h','replicate').^2); %距离变换
sqoc=imclose(imopen(sq,ones(3,3)),ones(3,3)); % 图像闭合运算
wsd=watershed(sqoc);
wg=wsd==0;
rfg=F;
rfg(wg)=255;
im=imextendedmin(F,20);
Lim=watershed(bwdist(im));
em=Lim==0;
rfmin=imimposemin(sq,im|em);
wsdmin=watershed(rfmin);
rfgm=F;
rfgm(wsdmin==0)=255;