📄 extractresults2freesurfer.html
字号:
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" "http://www.w3.org/TR/REC-html40/loose.dtd"><html><head> <title>Description of extractresults2freesurfer</title> <meta name="keywords" content="extractresults2freesurfer"> <meta name="description" content="Loop over all controls and patients"> <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1"> <meta name="generator" content="m2html © 2003 Guillaume Flandin"> <meta name="robots" content="index, follow"> <link type="text/css" rel="stylesheet" href="../m2html.css"></head><body><a name="_top"></a><div><a href="../index.html">Home</a> > <a href="index.html">mri_toolbox</a> > extractresults2freesurfer.m</div><!--<table width="100%"><tr><td align="left"><a href="../index.html"><img alt="<" border="0" src="../left.png"> Master index</a></td><td align="right"><a href="index.html">Index for mri_toolbox <img alt=">" border="0" src="../right.png"></a></td></tr></table>--><h1>extractresults2freesurfer</h1><h2><a name="_name"></a>PURPOSE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2><div class="box"><strong>Loop over all controls and patients</strong></div><h2><a name="_synopsis"></a>SYNOPSIS <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2><div class="box"><strong>This is a script file. </strong></div><h2><a name="_description"></a>DESCRIPTION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2><div class="fragment"><pre class="comment"> Loop over all controls and patients</pre></div><!-- crossreference --><h2><a name="_cross"></a>CROSS-REFERENCE INFORMATION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2>This function calls:<ul style="list-style-image:url(../matlabicon.gif)"></ul>This function is called by:<ul style="list-style-image:url(../matlabicon.gif)"></ul><!-- crossreference --><h2><a name="_source"></a>SOURCE CODE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2><div class="fragment"><pre>0001 0002 <span class="comment">% Loop over all controls and patients</span>0003 0004 <span class="comment">% Search ImageGridTime for selected time point</span>0005 <span class="comment">% Extract Cdata from ImageGridAmp at selected time (column)</span>0006 0007 <span class="comment">%----------------------------------------------</span>0008 <span class="comment">% Define component timing</span>0009 0010 sa=[100 <span class="comment">% first positive peak P100 (OT)</span>0011 150 <span class="comment">% positive & negative peak P150 (SF) / N150 (OT)</span>0012 180 <span class="comment">% negative N180 @ PT,AT and SP,SC</span>0013 250 <span class="comment">% positive P250 (SP)</span>0014 400 ]; <span class="comment">% P400 (SC,SP)</span>0015 0016 wm=[ 800017 90 0018 1500019 2600020 3000021 4000022 550 ];0023 0024 ea=[ 80 <span class="comment">% N80 @ SF , P80 @ OT/PT</span>0025 145 <span class="comment">% N150 @ OT/PT , P150 @ SPF (SC)</span>0026 240 <span class="comment">% , P240 @ SPF/IPF & OT/SP</span>0027 300 <span class="comment">% N300 @ L_PT/OT/IP,</span>0028 350 <span class="comment">% Series: P350 @ Frontal (IPF), P500 @ SP</span>0029 450 <span class="comment">% Series: P350 @ Frontal (IPF), P500 @ SP</span>0030 550 ];<span class="comment">% Series: P350 @ Frontal (IPF), P500 @ SP; Also N550 @ left frontal</span>0031 0032 <span class="comment">%----------------------------------------------</span>0033 <span class="comment">% Define groups</span>0034 groups = {<span class="string">'c'</span>,<span class="string">'p'</span>};0035 0036 <span class="comment">% cd E:\matlab\brainstorm_v1\studies; extractresults2freesurfer</span>0037 0038 <span class="comment">%----------------------------------------------</span>0039 <span class="comment">% Define subjects</span>0040 subs = [1:10];0041 0042 cond = {<span class="string">'ea'</span>,<span class="string">'sa'</span>,<span class="string">'wm'</span>};0043 0044 <span class="keyword">for</span> g=1:length(groups),0045 <span class="keyword">for</span> s=1:length(subs),0046 0047 sub = sprintf(<span class="string">'%s%02d'</span>,groups{g},subs(s));0048 0049 0050 <span class="comment">% temporarily skip subject c04</span>0051 <span class="keyword">if</span> strmatch(sub,<span class="string">'c04'</span>), <span class="keyword">continue</span>; <span class="keyword">end</span>0052 0053 0054 0055 <span class="comment">%----------------------------------------------------</span>0056 <span class="comment">% Load cortical tesselation from subjecttess.mat</span>0057 0058 subtessfile = [<span class="string">'E:\matlab\brainstorm_v1\subjects\'</span>,sub,filesep,sub,<span class="string">'_subjecttess.mat'</span>];0059 0060 <span class="keyword">if</span> exist(subtessfile) ~= 2,0061 <span class="comment">% skip this group/subject, as no data available</span>0062 fprintf(<span class="string">'\nskipping subject: %s\n\n'</span>,sub);0063 <span class="keyword">continue</span>;0064 <span class="keyword">else</span>0065 p.mesh.path = [<span class="string">'E:\matlab\brainstorm_v1\subjects\'</span>,sub,filesep];0066 p.mesh.file = [sub,<span class="string">'_subjecttess.mat'</span>];0067 p.mesh.type = <span class="string">'BrainStorm'</span>;0068 p = mesh_open(p);0069 0070 p.mesh.current = mesh_check(p,<span class="string">'cortex'</span>);0071 0072 <span class="comment">% effectively remove all extraneous surfaces for mesh_write_freesurfer below</span>0073 indices = strcmp(p.mesh.data.meshtype,<span class="string">'cortex'</span>);0074 exclude = find(indices == 0);0075 p.mesh.data.meshtype(exclude) = [];0076 0077 <span class="keyword">end</span>0078 0079 0080 <span class="comment">%----------------------------------------------------</span>0081 <span class="comment">% define cortical source activity matrix file names</span>0082 0083 <span class="comment">%c01_oac_volts_data_results_image_svd12.mat</span>0084 <span class="comment">%c01_ouc_volts_data_results_image_svd12.mat</span>0085 oucfile = [<span class="string">'E:\matlab\brainstorm_v1\studies\'</span>,sub,filesep,sub,<span class="string">'_ouc_volts_data_results_image_svd12.mat'</span>];0086 oacfile = [<span class="string">'E:\matlab\brainstorm_v1\studies\'</span>,sub,filesep,sub,<span class="string">'_oac_volts_data_results_image_svd12.mat'</span>];0087 oatfile = [<span class="string">'E:\matlab\brainstorm_v1\studies\'</span>,sub,filesep,sub,<span class="string">'_oat_volts_data_results_image_svd12.mat'</span>];0088 tacfile = [<span class="string">'E:\matlab\brainstorm_v1\studies\'</span>,sub,filesep,sub,<span class="string">'_tac_volts_data_results_image_svd12.mat'</span>];0089 0090 <span class="keyword">if</span> exist(oucfile) ~= 2, error([<span class="string">'no '</span>,oucfile]); <span class="keyword">end</span>;0091 <span class="keyword">if</span> exist(oacfile) ~= 2, error([<span class="string">'no '</span>,oacfile]); <span class="keyword">end</span>;0092 <span class="keyword">if</span> exist(oatfile) ~= 2, error([<span class="string">'no '</span>,oatfile]); <span class="keyword">end</span>;0093 <span class="keyword">if</span> exist(tacfile) ~= 2, error([<span class="string">'no '</span>,tacfile]); <span class="keyword">end</span>;0094 0095 <span class="comment">%----------------------------------------------------</span>0096 <span class="comment">% Save Cdata to freesurfer curvature files</span>0097 0098 <span class="comment">% D:\freesurfer\subjects\ptsdpet-c01\surf</span>0099 p.mesh.path = sprintf(<span class="string">'D:\\freesurfer\\subjects\\ptsdpet-%s\\surf\\'</span>,sub);0100 0101 p.mesh.type = <span class="string">'FS_curv'</span>; <span class="comment">% output freesurfer binary curvature files</span>0102 0103 <span class="comment">% cortical values are scaled into nano A.m (dipole moment)</span>0104 scale = 10^9;0105 0106 0107 <span class="keyword">for</span> c = 1:length(cond),0108 0109 0110 <span class="keyword">switch</span> cond{c},0111 0112 <span class="keyword">case</span> <span class="string">'sa'</span>,0113 0114 <span class="comment">%---------------------------------------------------------------------</span>0115 <span class="comment">% SA condition</span>0116 0117 vert = size(p.mesh.data.vertices{p.mesh.current},1);0118 times = length(sa);0119 0120 Dif.oac = zeros(vert,times);0121 Dif.ouc = Dif.oac;0122 0123 load(oacfile,<span class="string">'ImageGridAmp'</span>,<span class="string">'ImageGridTime'</span>); <span class="comment">% Time is in sec here</span>⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -