方法一

主要函数：
graythresh() ; im2bw() ; bwperim() ;

代码如下

img=imread('ball.jpg'); %原图像为灰度图
subplot(1,2,1);
imshow(img);
title('原图像');

threshold= graythresh(img);  %使用最大类间方差法找到图片的一个合适的阈值(threshold=0.4745)
img_1=im2bw(img,threshold); %用这个阈值把灰度图像转换为二值图像(像素值>0.4745*255的变为255白，其他的变为0黑)

img_2=bwperim(img_1); %用bwperim函数查找二值图像的边缘

subplot(1,2,2);
imshow(img_2);
title('边缘提取后的图像')

运行结果

方法二

主要函数：
entropyfilt() ; bwareaopen() ; rangefilt() ;

代码如下

img=imread('ball.jpg'); 
subplot(2,3,1);
imshow(img);
title('原图像');

img_1 = entropyfilt(img);  %创建纹理图像
img_2 = mat2gray(img_1);  %转化为灰度图像
subplot(2,3,2)
imshow(img_2);  
title('纹理图像');

img_3 = im2bw(img_2,0.7);  %转化为二值图像
subplot(2,3,3);
imshow(img_3);  
title('二值图像');

img_4 = bwareaopen((1-img_3),3000,8); %提取大小球区域纹理
subplot(2,3,4);
imshow(img_4); 
title('不同区域纹理图像');

img_5= rangefilt (img,ones(7));  %rangefilt滤波
subplot(2,3,5); 
imshow(img_5); 
title('滤波后图像');

运行结果

对函数理解还不够深刻，如有错误望指正！