计算将灰度图像转化为二值图像所需的门限BWim2bwI1threshold将灰度图像转化为二值图像figureimshowBWtitle二值图像dimsizeBWcolrou
ddim2290计算起始点列坐标rowfi
dBWcol1计算起始点行坐标co
ectivity8
um_poi
ts180co
tourbwtracebou
daryBWrowcolNco
ectivity
um_poi
ts提取边界figureimshowI1holdo
plotco
tour2co
tour1gLi
eWidth2title边界跟踪图像12Hough变换Iimreadxia
bmprotIrgb2grayIsubplot221imshowrotItitle灰度图像axis5025050200grido
axiso
BWedgerotIprewittsubplot222imshowBWtitleprewitt算子边缘检测后图像axis5025050200grido
axiso
HTRhoughBWsubplot223imshowHXDataTYDataRI
itialMag
ificatio
fittitle霍夫变换图xlabelthetaylabelrhoaxiso
axis
ormalholdo
PhoughpeaksH5thresholdceil03maxH
fxTP2yRP1plotxyscolorwhiteli
eshoughli
esBWTRPFillGap5Mi
Le
gth7subplot224imshowrotItitle霍夫变换图像检测axis5025050200grido
axiso
holdo
max_le
0fork1le
gthli
esxyli
eskpoi
t1li
eskpoi
t2plotxy1xy2Li
eWidth2Colorgree
plotxy11xy12xLi
eWidth2Coloryellowplotxy21xy22xLi
eWidth2Colorredle
ormli
eskpoi
t1li
eskpoi
t2ifle
max_le
max_le
le
xy_lo
gxye
de
dplotxy_lo
g1xy_lo
g2Li
eWidth2Colorcya
13直方图阈值法用MATLAB实现直方图阈值法：Iimreadxia
bmpI1rgb2grayIfiguresubplot221imshowI1title灰度图像axis5025050200grido
显示网格线axiso
显示坐标系m
sizeI1测量图像尺寸参数GPzeros1256预创建存放灰度出现概率的向量fork0255GPk1le
gthfi
dI1km
计算每级灰度出现的概率，将其存入GP中相应位置e
dsubplot222bar0255GPg绘制直方图title灰度直方图xlabel灰度值ylabel出现概率I2im2bwI150255
fsubplot223imshowI2title阈值150的分割图像axis5025050200grido
显示网格线axiso
显示坐标系I3im2bwI200255subplot224imshowI3title阈值200的分割图像axis5025050200grido
显示网格线axiso
显示坐标系14自动阈值法：Otsu法用MATLAB实现Otsu算法：clcclearallIimreadxia
bmpsubplot121imshowItitle原始图像axis5025050200grido
显示网格线axiso
显示坐标系levelgraythreshI确定灰度阈值BWim2bwIlevelsubplot122imshowBWtitleOtsu法阈值分割图像axis5025050200grido
显示网格线axiso
显示坐标系15膨胀操作Iimreadxia
r