📄 load_cor.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 load_cor</title> <meta name="keywords" content="load_cor"> <meta name="description" content=""> <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> > load_cor.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>load_cor</h1><h2><a name="_name"></a>PURPOSE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2><div class="box"><strong></strong></div><h2><a name="_synopsis"></a>SYNOPSIS <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2><div class="box"><strong>function [cor,M] = fmri_ldcor(subject,arg2,arg3) </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">
Loads the indicated corronal slices from
$SUBJECTS_DIR/subject/mri/volumeid
cor = fmri_ldcor(subject)
cor = fmri_ldcor(subject,volumeid)
cor = fmri_ldcor(subject,slices)
cor = fmri_ldcor(subject,volid,slices)
If unspecified, volumeid defaults to T1.
If unspecified, nslices defaults to [1:256].
$Id: load_cor.m,v 1.2 2004/02/07 01:41:51 psdlw Exp $</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 <a name="_sub0" href="#_subfunctions" class="code">function [cor,M] = fmri_ldcor(subject,arg2,arg3)</a>0002 <span class="comment">%</span>0003 <span class="comment">% Loads the indicated corronal slices from</span>0004 <span class="comment">% $SUBJECTS_DIR/subject/mri/volumeid</span>0005 <span class="comment">%</span>0006 <span class="comment">% cor = fmri_ldcor(subject)</span>0007 <span class="comment">% cor = fmri_ldcor(subject,volumeid)</span>0008 <span class="comment">% cor = fmri_ldcor(subject,slices)</span>0009 <span class="comment">% cor = fmri_ldcor(subject,volid,slices)</span>0010 <span class="comment">%</span>0011 <span class="comment">% If unspecified, volumeid defaults to T1.</span>0012 <span class="comment">% If unspecified, nslices defaults to [1:256].</span>0013 <span class="comment">%</span>0014 <span class="comment">% $Id: load_cor.m,v 1.2 2004/02/07 01:41:51 psdlw Exp $</span>0015 0016 <span class="keyword">if</span>(nargin < 1 | nargin > 3)0017 msg = <span class="string">'USAGE: cor = fmri_ldcor(subject,<volid>,<slices>)'</span>;0018 qoe(msg);error(msg);0019 <span class="keyword">end</span>0020 0021 volumeid = <span class="string">'T1'</span>;0022 slices = [1:256];0023 0024 <span class="keyword">if</span>(nargin == 2)0025 <span class="keyword">if</span>(ischar(arg2)) volumeid = arg2;0026 <span class="keyword">else</span> slices = arg2;0027 <span class="keyword">end</span>0028 <span class="keyword">end</span>0029 0030 <span class="keyword">if</span>(nargin == 3)0031 <span class="keyword">if</span>(ischar(arg2)) volumeid = arg2;0032 <span class="keyword">else</span> slices = arg2;0033 <span class="keyword">end</span>0034 <span class="keyword">if</span>(ischar(arg3)) volumeid = arg3;0035 <span class="keyword">else</span> slices = arg3;0036 <span class="keyword">end</span>0037 <span class="keyword">end</span>0038 0039 <span class="keyword">if</span>(min(slices) < 1)0040 msg = sprintf(<span class="string">'Min Slice No = %d, must be > 0'</span>,min(slices));0041 qoe(msg);error(msg);0042 <span class="keyword">end</span>0043 0044 <span class="keyword">if</span>(max(slices) > 256)0045 msg = sprintf(<span class="string">'Max Slice No = %d, must be <= 256'</span>,max(slices));0046 qoe(msg);error(msg);0047 <span class="keyword">end</span>0048 0049 0050 subject = deblank(subject);0051 volumeid = deblank(volumeid);0052 0053 SubjectsDir = deblank(getenv(<span class="string">'SUBJECTS_DIR'</span>));0054 <span class="keyword">if</span>(isempty(SubjectsDir))0055 msg = <span class="string">'Cannot find SUBJECTS_DIR environment variable'</span>;0056 qoe(msg);error(msg);0057 <span class="keyword">end</span>0058 0059 SubjDir = strcat(SubjectsDir,<span class="string">'/'</span>,subject);0060 CorDir = strcat(SubjDir,<span class="string">'/mri/'</span>,volumeid);0061 0062 d = dir(CorDir);0063 <span class="keyword">if</span>(isempty(d))0064 CorDir = deblank(subject);0065 d = dir(CorDir);0066 <span class="keyword">if</span>(isempty(d))0067 msg = sprintf(<span class="string">'Directory %s does not exist (DIR=%s, SD=%s)'</span>, <span class="keyword">...</span>0068 CorDir,volumeid, SubjectsDir);0069 qoe(msg);error(msg);0070 <span class="keyword">end</span>0071 <span class="keyword">end</span>0072 0073 nslices = length(slices);0074 cor = zeros(256,nslices,256);0075 0076 Endian = 0;0077 precision = <span class="string">'uint8'</span>;0078 Nv = 256*256;0079 0080 fprintf(1,<span class="string">'Loading coronals from subject %s, volume %s, %d slices ... \n'</span>,<span class="keyword">...</span>0081 subject,volumeid,nslices);0082 0083 <span class="keyword">for</span> s = 1:nslices,0084 0085 n = slices(s);0086 0087 corfile = sprintf(<span class="string">'%s/COR-%03d'</span>,CorDir,n);0088 d = dir(corfile);0089 <span class="keyword">if</span>(isempty(d))0090 msg = sprintf(<span class="string">'File %s does not exist'</span>,corfile);0091 qoe(msg);error(msg);0092 <span class="keyword">end</span>0093 0094 <span class="comment">%%%% Open the corfile %%%%%</span>0095 <span class="keyword">if</span>(Endian == 0) fid=fopen(corfile,<span class="string">'r'</span>,<span class="string">'b'</span>); <span class="comment">% Big-Endian</span>0096 <span class="keyword">else</span> fid=fopen(corfile,<span class="string">'r'</span>,<span class="string">'l'</span>); <span class="comment">% Little-Endian</span>0097 <span class="keyword">end</span>0098 <span class="keyword">if</span>(fid == -1)0099 msg = sprintf(<span class="string">'Could not open %s for reading.'</span>,corfile); 0100 qoe(msg); error(msg);0101 <span class="keyword">end</span>0102 0103 <span class="comment">%%% Read the file in corfile %%%</span>0104 z = fread(fid,Nv,precision);0105 cor(:,:,s) = reshape(z, [256 256])'; <span class="comment">%' transpose for row major</span>0106 fclose(fid); 0107 0108 <span class="keyword">end</span>0109 0110 <span class="comment">%cor = permute(cor, [3 2 1]);</span>0111 0112 M = [0,-1,0,129;0,0,1,-129;-1,0,0,129;0,0,0,1];0113 0114 <span class="keyword">return</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> © 2003</address></body></html>⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -