实现图像处理方法

rgb = imread(‘s1.bmp’);
I = rgb2gray(rgb);
[n1,n2] = size(I);
%灰度图
figure,imshow(I),title(‘灰度图’)
tic
h = ones(9)/81;
I = uint8(conv2(I,h));
figure,imshow(I),title(‘线性均值滤波’)
BW = imbinarize(I);
figure,imshow(BW),title(‘二值化’)
B = ones(21);%结构元素
BW = -imerode(BW,B) + BW;
figure,imshow(BW),title(‘形态学边界提取’)
BW = bwmorph(BW, ‘thicken’);
figure,imshow(BW),title(‘加粗边界’)
BW = not(bwareaopen(not(BW), 300));
figure,imshow(BW),title(‘把空洞填了’)
%进行形态学运算
B = strel(‘line’50,90);
BW = imdilate(BW,B);
BW = imerode(BW,B);
figure,imshow(BW),title(‘再闭操作之后’)
B = strel('line',10,0);
BW = imerode(BW,B);
figure,imshow(BW),title(‘闭操作之后再腐蚀’)
%%BW = gpuArray(BW);%%将数据载入gpu可以加速，电脑不一定支持
 
%最小化背景
%细分
div =