Description
December 17, 2011
The
following MATLAB code will
load the both outer and inner cortical surface mesh format
into MatLab. It is tested on MATLAB 2009 (version 7.9,
64-bit Mac intel). The complete data set is used in
studies [1] [2] [4]. If you want to use the data set for publication,
please contact Moo K. Chung.
load
autism.mesh.coord.mat
The data
set contains the both inner and outer cortical surfaces of 16
autistic and 11 control subjects. In the mat
file, the variable autism_inner is 3 dimensional array of size
16 x 3 x 40962. Each dimension indexces subject number,
(x,y,z) coordinates and coordinate values. To visualize
the inner surface of the 1st autistic subject, need to
reformat the data structure.
surf.vertices=squeeze(autism_inner(1,:,:))';
surf.faces=tri;
figure;
figure_wire(surf, 'yellow',
'white');
The outer
surface is visualized as follows:
surf.vertices=squeeze(autism_outer(1,:,:))';
surf.faces=tri;
figure;
figure_wire(surf, 'yellow', 'white');
To
visualize the y-coordinate values on the outer surface
figure;
figure_trimesh(surf,surf.vertices(:,2))
The
cortical thickness for the first autism subject is computed as
thickness=squeeze(L2norm(autism_outer(1,:,:)
- autism_inner(1,:,:)));
figure;
figure_trimesh(surf,thickness)
The
cortical thickness values are highly noise so we need to perform heat kernel
smoothing. See [1] and [2] for the detailed
implemenation of heat kernel smoothing. We can also extract
cortical thickness values at vertex 10 for all subjects
thickness1=L2norm(autism_outer(:,:,10) - autism_inner(:,:,10))
thickness2=L2norm(control_outer(:,:,10)
- control_inner(:,:,10))
figure;
plot(thickness1,'or');
hold
on;
plot(thickness2,'ob');
MNI surface mesh format
May 9
2006.
The following MATLAB
code will
load MNI cortical mesh format into MatLab. It only read
ASCII-based MNI
format. MNI surface mesh format always ends with *obj
extension. This is based on MATLAB
6.1. implementation. The MNI
cortical mesh datat set has
been used in the studies [1-4].
References
December 17,
2011