应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 数学综合 » 请问多糖的TEM图可否用分形的概念来解析呢？有图
412/5返回列表上一页12345下一页
查看: 4669  |  回复: 40
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

feixiaolin

荣誉版主 (文坛精英)

优秀版主
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
10楼: Originally posted by ppt1210 at 2014-05-08 17:11:45
已赞助，谢谢赐教。
可否告知，所得信息能说明什么问题？以及如果解析得到，谢谢！...

C code
//  定义待用变量
#define int i, j, k, i1, j1;
                int kk[6], ll[6], gmax, gmin;
                int s[6]={2,3,4,6,8,12};
                double  M=24, D[6], Nr[6];

// 算法部分
for(k=0; k<6; k++)
{
        kk[k]=0;
        ll[k]=0;
        Nr[k]=0;
   
        // read image
        for(i=0; i<nWidth-1-s[k], i+=s[k])
                for(j=0; j<nHight-1-s[k], j+=s[k])
                {
                        max=0;
                        min=255;
                        for(i1=i; i1<i+s[k], i1+=s[k])
                                for(j1=j; j1<j+s[k], j1+=s[k])
                                {
                                        if(gray[i1, j1]<min)
                                                min=gray[i1, j1];
                                        if(gray[i1, j1]>max)
                                                max=gray[i1, j1];
                                }
                        kk[k]+=floor((float)min/s[k]);
                        ll[k]+=floor((float)max/s[k] + 0.5);
                        Nr[k]+=ll[k] - kk[k] + 1;                               
                }
        D[k]=log(N[k])/log(M/s[k]);
}

for(k=1; k<6; k++)
        D[0]+=D[k];
D[0]/=6.0;

Matlab code
% 该程序只能计算size为NxN（N=2^n）的图形。
function fd=box_frac_dem(X);
% 差分盒维数
% Example:
%    X=double(imread('rice.tif'));
%    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;
   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);
赞一下(10人) 回复此楼

» 本帖附件资源列表
11楼2014-05-08 17:17:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ppt1210

至尊木虫 (著名写手)
引用回帖:
11楼: Originally posted by feixiaolin at 2014-05-08 17:17:46
code
//  定义待用变量
#define int i, j, k, i1, j1;
                int kk, ll, gmax, gmin;
                int s={2,3,4,6,8,12};
                double  M=24, D, Nr;

// 算法部分
for(k=0; k<6; k++)
{
        kk=0;
        ll=0;
        Nr= ...

非常感谢，这个是用什么软件来处理吗？很抱歉，我对这方面的了解几乎为0
赞一下 回复此楼
我对人生无悔，奈何心中有愧。。。
12楼2014-05-09 09:06:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
12楼: Originally posted by ppt1210 at 2014-05-09 09:06:14
非常感谢，这个是用什么软件来处理吗？很抱歉，我对这方面的了解几乎为0...

图像处理软件做的。
赞一下 回复此楼
13楼2014-05-09 09:33:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ppt1210

至尊木虫 (著名写手)
引用回帖:
13楼: Originally posted by feixiaolin at 2014-05-09 09:33:25
图像处理软件做的。...

是matlab吗？
回复此楼
我对人生无悔，奈何心中有愧。。。
14楼2014-05-09 10:09:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
引用回帖:
14楼: Originally posted by ppt1210 at 2014-05-09 10:09:19
是matlab吗？...

看11楼
回复此楼
15楼2014-05-09 10:19:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ppt1210

至尊木虫 (著名写手)
引用回帖:
15楼: Originally posted by feixiaolin at 2014-05-09 10:19:53
看11楼

好的，谢谢，只是，我不知道，如果我要写到论文里面，可以写到什么结论呢？
赞一下 回复此楼
我对人生无悔，奈何心中有愧。。。
16楼2014-05-09 10:27:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
16楼: Originally posted by ppt1210 at 2014-05-09 10:27:36
好的，谢谢，只是，我不知道，如果我要写到论文里面，可以写到什么结论呢？...

经过分析，分形维数是多少云云。
赞一下 回复此楼
17楼2014-05-09 11:00:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ppt1210

至尊木虫 (著名写手)
引用回帖:
17楼: Originally posted by feixiaolin at 2014-05-09 11:00:10
经过分析，分形维数是多少云云。...

OK，多谢您非常耐心的帮助
赞一下 回复此楼
我对人生无悔，奈何心中有愧。。。
18楼2014-05-09 11:09:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

门前绝壁

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
3楼: Originally posted by feixiaolin at 2014-05-07 16:09:58
完全行。

层主你好，我想做电镜扫描图片的分形维数，由灰度图片计算记盒维数，向层主求救，维数不是应该在1.几迈，我看这个维数达到了十几π_π

[ 发自小木虫客户端 ]
赞一下(2人) 回复此楼
19楼2014-11-11 17:46:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖仅楼主可见
赞一下
20楼2015-02-01 22:02:04
已阅   申请数学EPI   回复此楼   编辑   查看我的主页
相关版块跳转 我要订阅楼主 ppt1210 的主题更新
412/5返回列表上一页12345下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 10/500 2026-08-21 17:32 by chengyan1220
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] filecode，4个jtjc了 +11 ziyangfang 2026-08-19 14/700 2026-08-21 16:51 by applepetty
[基金申请] 让我中一个面上吧！ +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 人气不行了 +7 fansofjerry 2026-08-21 7/350 2026-08-21 14:13 by Vivilian
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +15 医学老男孩 2026-08-20 18/900 2026-08-21 13:46 by 医学老男孩
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 只有每年这种时候来逛逛小木虫 +23 yaoyewhu2008 2026-08-20 25/1250 2026-08-21 12:19 by Poppy1104
[基金申请] 什么时候开奖？ +7 CrisMessi 2026-08-18 7/350 2026-08-21 12:15 by weiyin
[基金申请] 今天基金会出结果吗？20260819 +15 kkkl_v 2026-08-19 16/800 2026-08-21 11:22 by 365372687
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
信息提示
关闭
请填处理意见
关闭