function Divided_image = Divide_into_block(I,m,n)
%把二维图像矩阵I分割成m*n的矩阵块，当边界区域不足m*n时，用0填充
Divided_image_new = Divided_image_judge(I,m,n);
[line,row] = size(Divided_image_new);
Divided_image = mat2cell(Divided_image_new,ones(line/m,1)*m,ones(row/n,1)*n);
end

function Divided_image_new = Divided_image_judge(I,m,n)
%函数用来判断I矩阵是否可以被m*n整数划分，如果不能，则把矩阵边界扩充，扩充部分用0填充
%I表示输入矩阵；m*n表示分块大小
[line,row] = size(I);
rem_m = rem(line,m);
rem_n = rem(row,n);
if rem_m == 0 & rem_n == 0
    Divided_image_new = I;
else
    if rem_m > 0
        line_new =line+m-rem_m;
    end
    if rem_n > 0
        row_new = row+n-rem_n;
    end
    Divided_image_new = zeros(line_new,row_new);
    Divided_image_new(1:line,1:row) = I;
end
end