| 查看: 1106 | 回复: 2 | |||
[交流]
图片分形维数计算 已有2人参与
|
|
如何用盒子法计算图片分形维数？具体流程？ 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
材料调剂
已经有4人回复
-
面上模板改不了页边距吧？
已经有6人回复
-
307求调剂
已经有6人回复
-
304求调剂
已经有5人回复
-
317一志愿华南理工电气工程求调剂
已经有8人回复
-
272求调剂
已经有3人回复
-
化工专硕348，一志愿985求调剂
已经有6人回复
-
292求调剂
已经有3人回复
-
290求调剂
已经有6人回复
-
295求调剂
已经有5人回复
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
2楼2018-05-05 11:02:21
★
小木虫: 金币+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
