材料科学中处理组织分布是否均匀是否存在聚集的matlab图像处理分析

%该脚本由BinarizedSTatic改编，添加循环，目的是分块图像，并统计每块图像中的黑色比例
%该m用于将RGB图片转化为灰度图以及二值化，二值化阈值最初采用由大津法自动确定，后根据实际情况自己确定

data=imread('E:\16image\even\xxxx.jpg'); %读取图像为矩阵
data;
subplot(221);
imshow(data);
gdata=rgb2gray(data);
subplot(222);
imshow(gdata);
nOtsu=graythresh(data);%OTSU方法得到最佳阈值

BW=im2bw(data,0.25); %#ok<IM2BW> %灰度转化为二进制图像BW
subplot(223);
imshow(BW);

[m,n]=size(BW);

Rowi=20;
Columni=m/Rowi;   %越大标准差越大,纵向像素数

Rowj=20;
Columnj=n/Rowj;   %越大标准差越大，横向像素数

%分割矩阵，mx中rowi*columni要等于行像素数 rowj*columnj要等于列像素数
mx=Rowi*ones(1,Columni);
ny=Rowj*ones(1,Columnj);
slicearray=mat2cell(BW,mx,ny);

%以下统计图像像素出现频率
tbl = tabulate(BW(:));
Cahcesum0=0;
Cahcesum1=0;
flag=0;
mflag=0;
Zeropercent=ones(1,Col