0001 function vb_job_curry2matlab
0002
0003
0004
0005
0006 global define;
0007 define = vb_define_set;
0008
0009
0010
0011
0012 parm.dir.brain = './';
0013 parm.file.brain = '100008d.brain.mat';
0014
0015
0016
0017
0018 parm.dir.mri = ...
0019 '/home/vbmeg2/taku-y/realdata/data/SBIdata-new/mri/100008d/';
0020 parm.file.dicom = '100008d.tbl.mat';
0021
0022
0023
0024
0025 parm.dir.curry_srf = './';
0026 parm.file.curry_srf = '100008d.s01';
0027
0028
0029
0030
0031 parm.dir.head = './';
0032 parm.file.head = '100008d.s01.head.mat';
0033
0034
0035
0036
0037 parm.brain_parm.ignore_offset = 0;
0038
0039 [V,F,xx] = vb_read_curry_surface(parm);
0040 Vhead = V;
0041 Fhead = F;
0042 XXhead = xx;
0043
0044 patch('Vertices',V,'Faces',F,'FaceColor','none','EdgeColor','g');
0045 hold on;
0046 axis equal;
0047 quiver3(V(:,1),V(:,2),V(:,3),xx(:,1),xx(:,2),xx(:,3),1.0);
0048
0049 load([parm.dir.brain parm.file.brain],'V','F');
0050 patch('Vertices',V,'Faces',F.F3,'FaceColor','none',...
0051 'EdgeColor','r');
0052
0053 vb_fsave([parm.dir.head parm.file.head],'Vhead','Fhead','XXhead');
0054