使用MyRobustCrust使用Matlab从MNI脑模板创建3D线框皮质模型

时间:2015-03-24 20:57:08

标签: matlab wireframe neuroscience nifti

我正在尝试使用matlab上的MyRobustCrust从MNI脑模板创建线框3d模型。我的目标是能够确定皮质表面上的MNI坐标。

为此,请按照此网站上列出的说明进行操作:http://soundsfromthenorth.com/?p=3

不幸的是,管理员似乎不再活跃了。

按照网站的说明,我已成功从FSL中提取数据集(" MNI152_T1_2mm_brain_mask.nii")。但是,当我去运行MyRobustCrust时,使用以下输入(按照站点的指示),我得到一个空白图。

代码输入

>> [t,tnorm]=MyRobustCrust(p);
figure(1)
hold on
axis equal
h=trisurf(t,p(:,1),p(:,2),p(:,3),'facealpha',0.4,...
...'facecolor',[1 1 1],'edgecolor',[0.8 0.8 0.8]);
view(3);
axis vis3d

代码输出:

Added Shield: 0.0559 s
Delaunay Triangulation Time: 0.6515 s
Connectivity Time: 0.0796 s
Circumcenters Time: 0.0168 s
Warning: Brute continuation necessary 
> In MyRobustCrust>Marking (line 380)
  In MyRobustCrust (line 106) 
Warning: 1000 th level was reached\n 
> In MyRobustCrust>Marking (line 412)
  In MyRobustCrust (line 106) 
99.9941 % of Tetraedroms were checked
Walking Time: 0.2084 s
Manifold extraction Time: 0.3782 s
Total Time: 1.4809 s
 view(3);
        |
Error: Unbalanced or unexpected parenthesis or bracket.

不幸的是,我对神经影像学和Matlab都很陌生。有没有人对我在这里做错了什么有任何建议?错误信号是否表示MyRobustCrust或数据代码中的不平衡或意外括号?

1 个答案:

答案 0 :(得分:1)

我希望我能正确理解你,但在我看来,你试图将横向平面(轴面)标记带入MNI坐标。

如果是这样,那么也许你应该尝试xjView(注意:需要SPM)。这是一个图形界面,主要用于覆盖模板上的对比,其他功能的音调(创建蒙版,生成统计表等)。在您的情况下,您只需加载默认的MNI图像(来自SPM的avg152T1或T2)并单击右上角的3轴视图; MNI坐标位于左下方。

enter image description here