%計(jì)算粒度分布
%-----------------------------------------------------------------
clear?all;?clc;?close?all;?%清除內(nèi)存及數(shù)據(jù)、變量等;清除命令窗口；關(guān)閉所有圖形窗口
disp（‘提示：圖像轉(zhuǎn)換中，需數(shù)分鐘，請(qǐng)稍候...‘）

i=3;%處理的圖像
dirIn=[‘image\‘];?%數(shù)據(jù)文件輸入路徑（目錄）


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-*‘）