| 查看: 981 | 回复: 2 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
图片分形维数计算已有2人参与
|
|||
|
如何用盒子法计算图片分形维数？具体流程？ 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
第一篇论文投稿进程记录
已经有11人回复
-
博后网上办公系统维护
已经有22人回复
-
暑假孩子一个人在家可以干什么
已经有7人回复
-
专科院校申请国自然
已经有21人回复
-
offer选择
已经有3人回复
-
夏至，要不要硕博联谊
已经有8人回复
-
基金评审时，评审专家还回去看申请人代表作之外的文章吗？
已经有12人回复
-
CSC博士联培对将来就业有用吗
已经有5人回复
-
对象没有，还非常想读博，难以抉择
已经有52人回复
-
自然基金变化？？
已经有25人回复
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
老师，我有个分形维数的程序有些地方不明白，想请教一下您！小弟不胜感激！ A=imread('e:\1.jpg');% 读入图像，并转化为矩阵 t=graythresh(A);% 计算二值化阈值 B=im2bw(A,t);% 输出图像的二值化矩阵% image(B); [x,y]=size(B); u=1;% 图像尺寸 V=nonzeros(B); Area=sum(V)/(x*y); for side_length=2:110% 定义观测尺度，即正方形(盒子)边长 Hang=mod(x,side_length); Lie=mod(y,side_length); C=B(1:x-Hang,1:y-Lie);% 若图像大小“x×y”不能被正方形%“side_length ×side_length”整除，则舍去多余的“边幅”部分； [m,n]=size(C); X=reshape(C,side_length, numel(C)/side_length);interim1=sum(X);% 计数有像素(“1”)的盒子数Y=reshape(interim1,side_length,numel(interim1)/side_length); interim2=sum(Y); Number=numel(nonzeros(interim2)); interim=sum(interim2');%当盒子边长为side_length时，统计出各个 % 盒子内的“1”值之和“interim2”，以及所有盒子“1”值总和“interim”， W(u,1)= Number; u=u+1; endy=log(W); x=log(2:110); plot(x,y,'o'); Dbox=polyfit(x',y,1)% 线性方程斜率的绝对值即为盒子维数 B=nonzeros(B); Area % 附带给出图形(像素为1)，占整个图幅的面积比例 上面是分形维数的程序，但是有几个地方不明白，side_length=2:110% 定义观测尺度，即正方形(盒子)边长，这个110是什么意思？？？？？还有2应该是盒子边长吧？ |
3楼2018-09-27 19:18:20
5