24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wyjie2004

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 99.1
散金: 5
帖子: 21
在线: 7.4小时
虫号: 1032565
注册: 2010-05-31
性别: MM
专业: 环境工程

[求助] matlab计盒维数程序问题

function fd=box_frac_dem(X);
% 差分盒维数
% Example:
% X=double(imread('1-1-4.jpg));
% fd=box_frac_dem(X);
% Author email of the program:% zjliu2001@163.com
%
% Reference:
% Sarkar N,Chaudhuri B B. An efficient approach to estimate
% fractal dimensionof textural images [J].Pattern Recognition,
% 1992,25(9):1035-1041

if size(X,1)~=size(X,2);
error('The size of X must NxN.');
end
B=size(X,1);
if mod(log2(B),1)>0;
error('The size of X must 2^n');
end
t=log2(B);
s=2.^(1:t);
Nr=zeros(1,t);
for k=1:t;
d=s(k);
h=256/d; % Let B instead of 256 for the size [2^nX2^n]
for m=1:h;
for n=1:h;
A=X(d*(m-1)+[1:d],d*(n-1)+[1:d]);
mn=min(A(1:end));
mx=max(A(1:end));
nr=fix(mx/d)-fix(mn/d)+1;
Nr(k)=Nr(k)+nr;
end
end
end
r=B./s;
p=polyfit(log10(r),log10(Nr),1);
fd=p(1);

出现错误

回复此楼