图像分36块后，求局部熵

clc
clear all
img = imread('all.bmp');
%% resize the image into the new size with 500x*500y
[row, col] = size(img);
new_row = ceil(row/46) * 46;
new_col = ceil(col/46) * 46;
new_r_img = imresize(img, [new_row new_col], 'bilinear');
new_img = new_r_img;
[y_row y_col dim] = size(new_img);
row_blk_num = y_row/46;  % 46
col_blk_num = y_col/46;  % 46
%no overlap
blocks={};
k=1;
for i = 1:row_blk_num
    for j = 1:col_blk_num
        disp(blocks);
        block{k} = new_img((i - 1) * 46 + 1 : i * 46, (j - 1) * 46 + 1 : j * 46, :);
        k=k+1;
    end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:36
     p(i)=entropy(block{i});
end
sum=0;
for i=1:36
    sum=sum+(p(i)/30);
end