clear,clc,close all;
f=imread('c2.jpg');
g=rgb2gray(f);figure,imshow(g);
gth=imtophat(g,strel('disk',4));
figure,imshow(gth);
bw=im2bw(gth,0.1);figure,imshow(bw);
bwao=bwareaopen(bw,20);
figure,imshow(cat(2,bw,bwao));
bwedge=edge(bwao);
figure,imshow(bwedge);
bwao=bwareaopen(bwedge,20);
figure,imshow(cat(2,bwedge,bwao));
bwerode=imerode(bwao,strel('line',5,90));
figure,imshow(bwerode);
bwdilate=imdilate(bwerode,strel('disk',25));
figure,imshow(bwdilate);
bwao=bwareaopen(bwdilate,20000);
figure,imshow(bwao);
[r,c]=find(bwao);
result=g(min(r):max(r),min(c):max(c));
figure,imshow(result);