求助!!!重建ECT 算法
本帖最后由 冷凌川L 于 2012-6-19 19:04 编辑根据ECT数据重建图像。实验对象是一个6cm柱形容器内靠近管壁相对放置的两个有机玻璃棒,玻璃棒直径分别为2.5cm和1.5cm。使用12电极ECT传感器测量。给出灵敏度矩阵S,满场数据P2,空场数据P0,测量数据P1。另外nodeeit.m文件为重建图像剖分网格函数。可生成方形剖分网格数812个。要求是:根据所给数据,用一种图像重建算法重建ECT图像,使重建图像质量最好。数据见群共享压缩文件“ECT数据.rar”。
S:
P0:
P1:
P2:
function =nodeeit();
aa=32;
x=(-1+1/aa:2/aa:1-1/aa);
y=(-1+1/aa:2/aa:1-1/aa);
for i=1:aa
for j=1:aa
node(aa*(i-1)+j,1)=x(i);
node(aa*(i-1)+j,2)=y(j);
end
end
j=1;
for i=1:aa*aa
if node(i,1)^2+node(i,2)^2<=1
newnode(j,1)=node(i,1);
newnode(j,2)=node(i,2);
nodenum(j)=i;
j=j+1;
end
end
就这些资料。。。
页:
[1]