Matlab实现Sandbox方法计算分形维数

具体的算法讲解可以去这里看：Sandbox法计算分形维数

function [r,average_num,DIMENSION,R_2] = sdbox_dim(varargin)
%% [scale,num,dim,r^2] = sandbox_dim(Image, NPoints)
% 该程序的功能是计算SandBox分形维数
%% 参数输入：1.二值图像矩阵；2.不同尺度的点数（并非最终点数，后面会调整）
Image = varargin{1};
NPoints = varargin{2};
if NPoints~=fix(NPoints)
    error('NPoints应为一个整数')
end
%% 计算质心，以及距质心最近边界的距离
[row,column] = size(Image);
points = find(Image==1);
[temp1,temp2] = ind2sub(size(Image), points);
Index_all = cat(2,temp1,temp2);
mass_center_r = fix(sum(Index_all(:,1))/size(Index_all,1));
mass_center_c = fix(sum(Index_all(:,2))/size(Index_all,1));
L_max = min(min([mass_center_r-1,row-mass_center_r]),min([mass_center_c-1,column-mass_center_c]));
%% 获得距离各个边界L_max长度的区域内的所有非零点的坐标
roi = Image(L_max+1:row-L_max,L_max+1:column-L_max);
Origins = find(roi=