%計算粒度分布
%-----------------------------------------------------------------
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-*‘）