extractresults2freesurfer.html

mri_toolbox是一个工具用来MRI. 来自于SourceForge, 我上传这个软件,希望能结识对医疗软件感兴趣的兄弟.

0124                 <span class="keyword">for</span> t=1:length(sa),
0125                     
0126                     <span class="comment">% find column of ImageGridAmp corresponding to wm(t)</span>
0127                     col = find(ImageGridTime == sa(t)/1000);
0128                     <span class="keyword">if</span> isempty(col), error(<span class="string">'Cannot find time column in ImageGridTime\n'</span>); <span class="keyword">end</span>
0129                     <span class="keyword">if</span> length(col) &gt; 1, error(<span class="string">'&gt; 1 matching time column in ImageGridTime\n'</span>); <span class="keyword">end</span>
0130                     
0131                     p.mesh.data.Cdata{p.mesh.current} = abs(ImageGridAmp(:,col) .* scale);
0132                     Dif.oac(:,t) = p.mesh.data.Cdata{p.mesh.current};
0133                     
0134                     p.mesh.file = sprintf(<span class="string">'rh.sa_%04d_oac_float'</span>,sa(t));
0135                     mesh_write(p);
0136                     
0137                 <span class="keyword">end</span>
0138                 
0139                 clear ImageGridAmp ImageGridTime; pack;
0140                 
0141                 load(oucfile,<span class="string">'ImageGridAmp'</span>,<span class="string">'ImageGridTime'</span>); <span class="comment">% Time is in sec here</span>
0142                 <span class="keyword">for</span> t=1:length(sa),
0143                     
0144                     <span class="comment">% find column of ImageGridAmp corresponding to sa(t)</span>
0145                     col = find(ImageGridTime == sa(t)/1000);
0146                     <span class="keyword">if</span> isempty(col), error(<span class="string">'Cannot find time column in ImageGridTime\n'</span>); <span class="keyword">end</span>
0147                     <span class="keyword">if</span> length(col) &gt; 1, error(<span class="string">'&gt; 1 matching time column in ImageGridTime\n'</span>); <span class="keyword">end</span>
0148                     
0149                     p.mesh.data.Cdata{p.mesh.current} = abs(ImageGridAmp(:,col) .* scale);
0150                     Dif.ouc(:,t) = p.mesh.data.Cdata{p.mesh.current};
0151                     
0152                     p.mesh.file = sprintf(<span class="string">'rh.sa_%04d_ouc_float'</span>,sa(t));
0153                     mesh_write(p);
0154                 <span class="keyword">end</span>
0155                 
0156                 clear ImageGridAmp ImageGridTime; pack;
0157                 
0158                 <span class="keyword">for</span> t=1:length(sa),
0159                     
0160                     p.mesh.data.Cdata{p.mesh.current} = abs(Dif.oac(:,t) - Dif.ouc(:,t));
0161                     p.mesh.file = sprintf(<span class="string">'rh.sa_%04d_dif_float'</span>,sa(t));
0162                     mesh_write(p);
0163                     
0164                 <span class="keyword">end</span>
0165                 
0166                 clear Dif; pack;
0167                 
0168                 
0169                 
0170             <span class="keyword">case</span> <span class="string">'wm'</span>,
0171                 
0172                 <span class="comment">%---------------------------------------------------------------------</span>
0173                 <span class="comment">% WM condition</span>
0174                 
0175                 vert = size(p.mesh.data.vertices{p.mesh.current},1);
0176                 times = length(wm);
0177                 
0178                 Dif.tac = zeros(vert,times);
0179                 Dif.oac = Dif.tac;
0180                 
0181                 load(tacfile,<span class="string">'ImageGridAmp'</span>,<span class="string">'ImageGridTime'</span>); <span class="comment">% Time is in sec here</span>
0182                 <span class="keyword">for</span> t=1:length(wm),
0183                     
0184                     <span class="comment">% find column of ImageGridAmp corresponding to wm(t)</span>
0185                     col = find(ImageGridTime == wm(t)/1000);
0186                     <span class="keyword">if</span> isempty(col), error(<span class="string">'Cannot find time column in ImageGridTime\n'</span>); <span class="keyword">end</span>
0187                     <span class="keyword">if</span> length(col) &gt; 1, error(<span class="string">'&gt; 1 matching time column in ImageGridTime\n'</span>); <span class="keyword">end</span>
0188                     
0189                     p.mesh.data.Cdata{p.mesh.current} = abs(ImageGridAmp(:,col) .* scale);
0190                     Dif.tac(:,t) = p.mesh.data.Cdata{p.mesh.current};
0191                     
0192                     p.mesh.file = sprintf(<span class="string">'rh.wm_%04d_tac_float'</span>,wm(t));
0193                     mesh_write(p);
0194                     
0195                 <span class="keyword">end</span>
0196                 
0197                 clear ImageGridAmp ImageGridTime; pack;
0198                 
0199                 load(oacfile,<span class="string">'ImageGridAmp'</span>,<span class="string">'ImageGridTime'</span>); <span class="comment">% Time is in sec here</span>
0200                 <span class="keyword">for</span> t=1:length(wm),
0201                     
0202                     <span class="comment">% find column of ImageGridAmp corresponding to sa(t)</span>
0203                     col = find(ImageGridTime == wm(t)/1000);
0204                     <span class="keyword">if</span> isempty(col), error(<span class="string">'Cannot find time column in ImageGridTime\n'</span>); <span class="keyword">end</span>
0205                     <span class="keyword">if</span> length(col) &gt; 1, error(<span class="string">'&gt; 1 matching time column in ImageGridTime\n'</span>); <span class="keyword">end</span>
0206                     
0207                     p.mesh.data.Cdata{p.mesh.current} = abs(ImageGridAmp(:,col) .* scale);
0208                     Dif.oac(:,t) = p.mesh.data.Cdata{p.mesh.current};
0209                     
0210                     p.mesh.file = sprintf(<span class="string">'rh.wm_%04d_oac_float'</span>,wm(t));
0211                     mesh_write(p);
0212                 <span class="keyword">end</span>
0213                 
0214                 clear ImageGridAmp ImageGridTime; pack;
0215                 
0216                 <span class="keyword">for</span> t=1:length(wm),
0217                     
0218                     p.mesh.data.Cdata{p.mesh.current} = abs(Dif.tac(:,t) - Dif.oac(:,t));
0219                     p.mesh.file = sprintf(<span class="string">'rh.wm_%04d_dif_float'</span>,wm(t));
0220                     mesh_write(p);
0221                     
0222                 <span class="keyword">end</span>
0223                 
0224                 clear Dif; pack;
0225                 
0226                 
0227                 
0228             <span class="keyword">case</span> <span class="string">'ea'</span>,
0229                 
0230                 <span class="comment">%---------------------------------------------------------------------</span>
0231                 <span class="comment">% EA condition</span>
0232                 
0233                 vert = size(p.mesh.data.vertices{p.mesh.current},1);
0234                 times = length(ea);
0235                 
0236                 Dif.oat = zeros(vert,times);
0237                 Dif.oac = Dif.oat;
0238                 
0239                 load(oatfile,<span class="string">'ImageGridAmp'</span>,<span class="string">'ImageGridTime'</span>); <span class="comment">% Time is in sec here</span>
0240                 <span class="keyword">for</span> t=1:length(ea),
0241                     
0242                     <span class="comment">% find column of ImageGridAmp corresponding to sa(t)</span>
0243                     col = find(ImageGridTime == ea(t)/1000);
0244                     <span class="keyword">if</span> isempty(col), error(<span class="string">'Cannot find time column in ImageGridTime\n'</span>); <span class="keyword">end</span>
0245                     <span class="keyword">if</span> length(col) &gt; 1, error(<span class="string">'&gt; 1 matching time column in ImageGridTime\n'</span>); <span class="keyword">end</span>
0246                     
0247                     p.mesh.data.Cdata{p.mesh.current} = abs(ImageGridAmp(:,col) .* scale);
0248                     Dif.oat(:,t) = p.mesh.data.Cdata{p.mesh.current};
0249                     
0250                     p.mesh.file = sprintf(<span class="string">'rh.ea_%04d_oat_float'</span>,ea(t));
0251                     mesh_write(p);
0252                     
0253                 <span class="keyword">end</span>
0254                 
0255                 clear ImageGridAmp ImageGridTime; pack;
0256                 
0257                 load(oacfile,<span class="string">'ImageGridAmp'</span>,<span class="string">'ImageGridTime'</span>); <span class="comment">% Time is in sec here</span>
0258                 <span class="keyword">for</span> t=1:length(ea),
0259                     
0260                     <span class="comment">% find column of ImageGridAmp corresponding to sa(t)</span>
0261                     col = find(ImageGridTime == ea(t)/1000);
0262                     <span class="keyword">if</span> isempty(col), error(<span class="string">'Cannot find time column in ImageGridTime\n'</span>); <span class="keyword">end</span>
0263                     <span class="keyword">if</span> length(col) &gt; 1, error(<span class="string">'&gt; 1 matching time column in ImageGridTime\n'</span>); <span class="keyword">end</span>
0264                     
0265                     p.mesh.data.Cdata{p.mesh.current} = abs(ImageGridAmp(:,col) .* scale);
0266                     Dif.oac(:,t) = p.mesh.data.Cdata{p.mesh.current};
0267                     
0268                     p.mesh.file = sprintf(<span class="string">'rh.ea_%04d_oac_float'</span>,ea(t));
0269                     mesh_write(p);
0270                 <span class="keyword">end</span>
0271                 
0272                 clear ImageGridAmp ImageGridTime; pack;
0273                 
0274                 <span class="keyword">for</span> t=1:length(ea),
0275                     
0276                     p.mesh.data.Cdata{p.mesh.current} = abs(Dif.oat(:,t) - Dif.oac(:,t));
0277                     p.mesh.file = sprintf(<span class="string">'rh.ea_%04d_dif_float'</span>,ea(t));
0278                     mesh_write(p);
0279                     
0280                 <span class="keyword">end</span>
0281                 
0282                 clear Dif; pack;
0283             <span class="keyword">end</span>
0284         <span class="keyword">end</span>
0285     <span class="keyword">end</span>
0286 <span class="keyword">end</span></pre></div>
<hr><address>Generated on Fri 21-May-2004 12:38:21 by <strong><a href="http://www.artefact.tk/software/matlab/m2html/">m2html</a></strong> &copy; 2003</address>
</body>
</html>