📄 ctf_make_mri.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 ctf_make_mri</title> <meta name="keywords" content="ctf_make_mri"> <meta name="description" content="ctf_make_mri - create a CTF mri struct with zero image"> <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> > ctf_make_mri.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>ctf_make_mri</h1><h2><a name="_name"></a>PURPOSE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2><div class="box"><strong>ctf_make_mri - create a CTF mri struct with zero image</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 mri = ctf_make_mri </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"> ctf_make_mri - create a CTF mri struct with zero image
mri = ctf_make_mri
mri is a data struct similar to those returned from ctf_read_mri,
although this one return a zero matrix.
see also: <a href="ctf_read_mri.html" class="code" title="function mri = ctf_read_mri(file)">ctf_read_mri</a>, <a href="ctf_write_mri.html" class="code" title="function ctf_write_mri(mri, fileName, force)">ctf_write_mri</a> <>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> % < > % < DISCLAIMER: > % < > % < THIS PROGRAM IS INTENDED FOR RESEARCH PURPOSES ONLY. > % < THIS PROGRAM IS IN NO WAY INTENDED FOR CLINICAL OR > % < OFFICIAL USE. > % < > % <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<> %</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)"><li><a href="avw2ctf.html" class="code" title="function mri = avw2ctf(avw)">avw2ctf</a> mri = avw2ctf(avw)</li></ul><!-- crossreference --><h2><a name="_subfunctions"></a>SUBFUNCTIONS <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2><ul style="list-style-image:url(../matlabicon.gif)"><li><a href="#_sub1" class="code">function Version_2_Header = Version_2_Header_make,</a></li><li><a href="#_sub2" class="code">function HeadModel_Info = headModel_make,</a></li><li><a href="#_sub3" class="code">function Image_Info = imageInfo_make,</a></li></ul><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 mri = ctf_make_mri</a>0002 0003 <span class="comment">% ctf_make_mri - create a CTF mri struct with zero image</span>0004 <span class="comment">%</span>0005 <span class="comment">% mri = ctf_make_mri</span>0006 <span class="comment">%</span>0007 <span class="comment">% mri is a data struct similar to those returned from ctf_read_mri,</span>0008 <span class="comment">% although this one return a zero matrix.</span>0009 <span class="comment">%</span>0010 <span class="comment">% see also: ctf_read_mri, ctf_write_mri</span>0011 <span class="comment">%</span>0012 <span class="comment">% <>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> %</span>0013 <span class="comment">% < > %</span>0014 <span class="comment">% < DISCLAIMER: > %</span>0015 <span class="comment">% < > %</span>0016 <span class="comment">% < THIS PROGRAM IS INTENDED FOR RESEARCH PURPOSES ONLY. > %</span>0017 <span class="comment">% < THIS PROGRAM IS IN NO WAY INTENDED FOR CLINICAL OR > %</span>0018 <span class="comment">% < OFFICIAL USE. > %</span>0019 <span class="comment">% < > %</span>0020 <span class="comment">% <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<> %</span>0021 <span class="comment">%</span>0022 0023 0024 0025 <span class="comment">% $Revision: 1.1 $ $Date: 2004/05/21 18:58:11 $</span>0026 0027 <span class="comment">% Licence: GNU GPL, no implied or express warranties</span>0028 <span class="comment">% History: 08/2003, Darren.Weber_at_radiology.ucsf.edu</span>0029 <span class="comment">% - adapted from an appendex to CTF document</span>0030 <span class="comment">% MRIConverter.pdf, which is copied at the end of this</span>0031 <span class="comment">% function.</span>0032 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>0033 0034 0035 ver = <span class="string">'[$Revision: 1.1 $]'</span>;0036 fprintf(<span class="string">'\nCTF_MAKE_MRI [v%s]\n'</span>,ver(12:16)); tic;0037 0038 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>0039 <span class="comment">% create the header</span>0040 fprintf(<span class="string">'...creating header\n'</span>);0041 mri.hdr = Version_2_Header_make;0042 0043 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>0044 <span class="comment">% ensure the data is 16 bits (it can be 8 or 16)</span>0045 mri.hdr.dataSize = 2;0046 0047 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>0048 <span class="comment">% create a zero image matrix, CTF mri files are always 256 x 256 x 256</span>0049 fprintf(<span class="string">'...creating zero image\n'</span>);0050 0051 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>0052 <span class="comment">% imageOrientation, 0 = left on left, 1 = left on right</span>0053 mri.hdr.imageOrientation = 0;0054 fprintf(<span class="string">'...creating sagittal slices in neurological orientation (left is on the left)\n'</span>);0055 fprintf(<span class="string">'...+X left to right, +Y anterior to posterior, +Z superior to inferior\n'</span>);0056 0057 0058 mri.img = zeros(256,256,256);0059 0060 t=toc; fprintf(<span class="string">'...done (%5.2f sec).\n\n'</span>,t);0061 0062 <span class="keyword">return</span>0063 0064 0065 0066 0067 0068 0069 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>0070 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>0071 <a name="_sub1" href="#_subfunctions" class="code">function Version_2_Header = Version_2_Header_make,</a>0072 0073 Version_2_Header.identifierString = <span class="string">'CTF_MRI_FORMAT VER 2.2 '</span>;0074 0075 0076 Version_2_Header.imageSize = 256; <span class="comment">% always = 256</span>0077 Version_2_Header.dataSize = 2; <span class="comment">% 1 or 2 (bytes), 8 or 16 bits</span>0078 Version_2_Header.clippingRange = 65536; <span class="comment">% max. integer value of data</span>0079 Version_2_Header.imageOrientation = 1; <span class="comment">% eg., 0 = left on left, 1 = left on right</span>0080 0081 <span class="comment">% voxel dimensions in mm</span>0082 Version_2_Header.mmPerPixel_sagittal = 1.0;0083 Version_2_Header.mmPerPixel_coronal = 1.0;0084 Version_2_Header.mmPerPixel_axial = 1.0;0085 0086 Version_2_Header.HeadModel_Info = headModel_make; <span class="comment">% defined below...</span>0087 Version_2_Header.Image_Info = imageInfo_make; <span class="comment">% defined below...</span>0088 0089 <span class="comment">% voxel location of head origin</span>0090 Version_2_Header.headOrigin_sagittal = 128;0091 Version_2_Header.headOrigin_coronal = 128;0092 Version_2_Header.headOrigin_axial = 128;0093 0094 <span class="comment">% euler angles to align MR to head coordinate system (angles in degrees!)</span>0095 <span class="comment">% 1. rotate in coronal plane by this angle</span>0096 <span class="comment">% 2. rotate in sagittal plane by this angle</span>0097 <span class="comment">% 3. rotate in axial plane by this angle</span>0098 Version_2_Header.rotate_coronal = 0;0099 Version_2_Header.rotate_sagittal = 0;0100 Version_2_Header.rotate_axial = 0;0101 0102 Version_2_Header.orthogonalFlag = 0; <span class="comment">% if set then image is orthogonal</span>0103 Version_2_Header.interpolatedFlag = 0; <span class="comment">% if set than image was interpolated</span>0104 0105 <span class="comment">% original spacing between slices before interpolation to CTF format</span>0106 Version_2_Header.originalSliceThickness = 0;0107 0108 <span class="comment">% transformation matrix head->MRI [column][row]</span>0109 Version_2_Header.transformMatrix = eye(4);0110 0111 <span class="keyword">return</span>0112 0113 0114 0115 0116 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>0117 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>0118 <a name="_sub2" href="#_subfunctions" class="code">function HeadModel_Info = headModel_make,</a>0119 0120 <span class="comment">% this function is called from Version_2_Header_make</span>0121 0122 <span class="comment">% fid. point coordinate (in voxels)</span>0123 HeadModel_Info.Nasion_Sag = 0; <span class="comment">% nasion - sagittal</span>0124 HeadModel_Info.Nasion_Cor = 0; <span class="comment">% nasion - coronal</span>0125 HeadModel_Info.Nasion_Axi = 0; <span class="comment">% nasion - axial</span>0126 HeadModel_Info.LeftEar_Sag = 0; <span class="comment">% left ear - sagittal</span>0127 HeadModel_Info.LeftEar_Cor = 0; <span class="comment">% left ear - coronal</span>0128 HeadModel_Info.LeftEar_Axi = 0; <span class="comment">% left ear - axial</span>0129 HeadModel_Info.RightEar_Sag = 0; <span class="comment">% right ear - sagittal</span>0130 HeadModel_Info.RightEar_Cor = 0; <span class="comment">% right ear - coronal</span>0131 HeadModel_Info.RightEar_Axi = 0; <span class="comment">% right ear - axial</span>0132 0133 <span class="comment">% fid. point coordinate (in voxels)</span>0134 <span class="comment">% HeadModel_Info.Nasion_Sag = 128; % nasion - sagittal</span>0135 <span class="comment">% HeadModel_Info.Nasion_Cor = 64; % nasion - coronal</span>0136 <span class="comment">% HeadModel_Info.Nasion_Axi = 128; % nasion - axial</span>0137 <span class="comment">% HeadModel_Info.LeftEar_Sag = 64; % left ear - sagittal</span>0138 <span class="comment">% HeadModel_Info.LeftEar_Cor = 128; % left ear - coronal</span>0139 <span class="comment">% HeadModel_Info.LeftEar_Axi = 128; % left ear - axial</span>0140 <span class="comment">% HeadModel_Info.RightEar_Sag = 192; % right ear - sagittal</span>0141 <span class="comment">% HeadModel_Info.RightEar_Cor = 128; % right ear - coronal</span>0142 <span class="comment">% HeadModel_Info.RightEar_Axi = 128; % right ear - axial</span>0143 0144 <span class="comment">% default sphere origin</span>0145 HeadModel_Info.defaultSphereX = 0.0; <span class="comment">% sphere origin x coordinate ( in mm )</span>0146 HeadModel_Info.defaultSphereY = 0.0; <span class="comment">% sphere origin y coordinate ( in mm )</span>0147 HeadModel_Info.defaultSphereZ = 0.0; <span class="comment">% sphere origin z coordinate ( in mm )</span>0148 HeadModel_Info.defaultSphereRadius = 75.0; <span class="comment">% default sphere radius ( in mm )</span>0149 0150 <span class="keyword">return</span>0151 0152 0153 0154 0155 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>0156 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>0157 <a name="_sub3" href="#_subfunctions" class="code">function Image_Info = imageInfo_make,</a>0158 0159 <span class="comment">% this function is called from Version_2_Header_make</span>0160 0161 Image_Info.modality = 0; <span class="comment">% 0 = MRI, 1 = CT, 2 = PET, 3 = SPECT, 4 = OTHER</span>0162 Image_Info.manufacturerName = char(repmat(double(<span class="string">' '</span>),1,64));0163 Image_Info.instituteName = char(repmat(double(<span class="string">' '</span>),1,64));0164 Image_Info.patientID = char(repmat(double(<span class="string">' '</span>),1,32));0165 Image_Info.dateAndTime = char(repmat(double(<span class="string">' '</span>),1,32));0166 Image_Info.scanType = char(repmat(double(<span class="string">' '</span>),1,32));0167 Image_Info.contrastAgent = char(repmat(double(<span class="string">' '</span>),1,32));0168 Image_Info.imagedNucleus = char(repmat(double(<span class="string">' '</span>),1,32));0169 0170 Image_Info.Frequency = 0;0171 Image_Info.FieldStrength = 0;0172 Image_Info.EchoTime = 0;0173 Image_Info.RepetitionTime = 0;0174 Image_Info.InversionTime = 0;0175 Image_Info.FlipAngle = 0;0176 Image_Info.NoExcitations = 0;0177 Image_Info.NoAcquisitions = 0;0178 0179 Image_Info.commentString = char(repmat(double(<span class="string">' '</span>),1,256));0180 Image_Info.forFutureUse = char(repmat(double(<span class="string">' '</span>),1,64));0181 0182 <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 + -