ICode9

I=imread('附件B.PNG'); f=rgb2gray(I);%灰度化 BW=edge(f,'canny'); %edge:以灰度图像为输入，'canny'为边缘检测算子 % 输出BW为二值图像，边缘处为白（255）其余部分为黑（0） [row,col]=size(BW); rhomax=round((row*row+col*col)^0.5);%计算对角线距离 A=zero