计算粒度分布

%计算粒度分布
%-----------------------------------------------------------------
clear all; clc; close all; %清除内存及数据、变量等;清除命令窗口；关闭所有图形窗口
disp（‘提示：图像转换中，需数分钟，请稍候...‘）

i=3;%处理的图像
dirIn=[‘image\‘]; %数据文件输入路径（目录）


I=imread（[dirIn‘DSC_‘num2str（i）‘.JPG‘]）;


for counter=0:24
    remain=imopen（Istrel（‘disk‘counter））;
    intensity_area（counter+1）=sum（remain（:））;
end
subplot（131）;plot（intensity_area‘m-*‘）