Providing you have the image processing toolbox, you can follow the suggestions of @rayryeng and @Shai and do something like this:
Morphological filtering:
BW1=bwareaopen(YourOriginalBWImage, 100);
Filtered image:
Hough transform (adapted from this):
[H,theta,rho] = hough(BW1);
P = houghpeaks(H,7,'threshold',ceil(0.3*max(H(:))));
lines = houghlines(BW1,theta,rho,P,'FillGap',5,'MinLength',3);
figure, imshow(YourOriginalBWImage), hold on
for k = 1:length(lines)
xy = [lines(k).point1; lines(k).point2];
plot(xy(:,1),xy(:,2),'LineWidth',2,'Color','green');
end
Result: