24小时热门版块排行榜

返回列表

ppt1210

至尊木虫 (著名写手)

应助: 1 (幼儿园)
金币: 17232.6
散金: 848
红花: 24
帖子: 1848
在线: 323.5小时
虫号: 257328
注册: 2006-06-04
性别: GG
专业: 生物化工与食品化工

[交流] 请问多糖的TEM图可否用分形的概念来解析呢？有图已有11人参与

本人是学生物化学的，实验过程中所纯化的样品，在透射电镜下得到了如下的图，苦于不知如何解析，偶然了解到分形的理论，但是对这方面实在是认识浅薄，特发此贴求教专业人士，不知我的这图能否用分形的理论来解析？能得到什么样的结果？

望请不吝赐教，非常感激！

SHW50-25K_副本.jpg

SHW50H2O-4_副本.jpg

回复此楼

» 猜你喜欢

写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有3人回复
需要合成515-64-0，50g，能接单的留言已经有3人回复
最近几年招的学生写论文不引自己组发的文章已经有10人回复
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）已经有3人回复
临港实验室与上科大联培博士招生1名已经有8人回复
26申博自荐已经有7人回复
想换工作。大多数高校都是评职称时认可5年内在原单位取得的成果吗？已经有4人回复
带资进组求博导收留已经有9人回复
求助大佬们，伤口沾上了乙腈已经有6人回复

我对人生无悔，奈何心中有愧。。。

1楼 2014-05-07 08:30:52

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

回帖支持 ( 显示支持度最高的前 50 名 )

feixiaolin

荣誉版主 (文坛精英)

专家经验: +518
应助: 942 (博后)
贵宾: 1.275
金币: 3430
散金: 58785
红花: 532
沙发: 11
帖子: 24215
在线: 2601.8小时
虫号: 2139575
注册: 2012-11-21
专业: 光学信息获取与处理
管辖: 数学

★
小木虫: 金币+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);