ge_hdr_read.html

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

0271 se_hdr = setfield(se_hdr, <span class="string">'se_no '</span>, fread(fid,1,<span class="string">'int16'</span>));              <span class="comment">%Series Number%</span>
0272 se_hdr = setfield(se_hdr, <span class="string">'se_datetime'</span>, fread(fid,1,<span class="string">'int32'</span>));        <span class="comment">%Allocation Series Data/Time stamp%</span>
0273 se_hdr = setfield(se_hdr, <span class="string">'se_actual_dt'</span>, fread(fid,1,<span class="string">'int32'</span>));       <span class="comment">%Actual Series Data/Time stamp%</span>
0274 se_hdr = setfield(se_hdr, <span class="string">'se_desc'</span>, fread(fid,30,<span class="string">'uchar'</span>));       <span class="comment">%Series Description%</span>
0275 se_hdr = setfield(se_hdr, <span class="string">'pr_sysid'</span>, fread(fid,9,<span class="string">'uchar'</span>));       <span class="comment">%Primary Receiver Suite and Host%</span>
0276 se_hdr = setfield(se_hdr, <span class="string">'pansysid'</span>, fread(fid,9,<span class="string">'uchar'</span>));       <span class="comment">%Archiver Suite and Host%</span>
0277 se_hdr = setfield(se_hdr, <span class="string">'se_typ '</span>, fread(fid,1,<span class="string">'int16'</span>));             <span class="comment">%Series Type%</span>
0278 se_hdr = setfield(se_hdr, <span class="string">'se_source '</span>, fread(fid,1,<span class="string">'int16'</span>));          <span class="comment">%Series from which prescribed%</span>
0279 se_hdr = setfield(se_hdr, <span class="string">'se_plane '</span>, fread(fid,1,<span class="string">'int16'</span>));           <span class="comment">%Most-like Plane (for L/S)%</span>
0280 se_hdr = setfield(se_hdr, <span class="string">'scan_type '</span>, fread(fid,1,<span class="string">'int16'</span>));          <span class="comment">%Scout or Axial (for CT)%</span>
0281 se_hdr = setfield(se_hdr, <span class="string">'position'</span>, fread(fid,1,<span class="string">'int32'</span>));           <span class="comment">%Patient Position%</span>
0282 se_hdr = setfield(se_hdr, <span class="string">'entry'</span>, fread(fid,1,<span class="string">'int32'</span>));              <span class="comment">%Patient Entry%</span>
0283 se_hdr = setfield(se_hdr, <span class="string">'anref'</span>, fread(fid,3,<span class="string">'uchar'</span>));          <span class="comment">%Anatomical reference%</span>
0284 fseek(fid,1,0); <span class="comment">% 16-bit alignment</span>
0285 se_hdr = setfield(se_hdr, <span class="string">'lmhor'</span>, fread(fid,1,<span class="string">'float32'</span>));              <span class="comment">%Horizontal Landmark%</span>
0286 se_hdr = setfield(se_hdr, <span class="string">'prtcl'</span>, fread(fid,25,<span class="string">'uchar'</span>));         <span class="comment">%Scan Protocol Name%</span>
0287 fseek(fid,1,0); <span class="comment">% 16-bit alignment</span>
0288 se_hdr = setfield(se_hdr, <span class="string">'se_contrast '</span>, fread(fid,1,<span class="string">'int16'</span>));        <span class="comment">%Non-zero if &gt; 0 image used contrast(L/S)%</span>
0289 se_hdr = setfield(se_hdr, <span class="string">'start_ras'</span>, fread(fid,1,<span class="string">'uchar'</span>));          <span class="comment">%RAS letter for first scan location (L/S)%</span>
0290 fseek(fid,1,0); <span class="comment">% 16-bit alignment</span>
0291 <span class="keyword">if</span> byte_align; fseek(fid,2,0); <span class="keyword">end</span> <span class="comment">% 32-bit alignment</span>
0292 se_hdr = setfield(se_hdr, <span class="string">'start_loc'</span>, fread(fid,1,<span class="string">'float32'</span>));          <span class="comment">%First scan location (L/S)%</span>
0293 se_hdr = setfield(se_hdr, <span class="string">'end_ras'</span>, fread(fid,1,<span class="string">'uchar'</span>));            <span class="comment">%RAS letter for last scan location (L/S)%</span>
0294 fseek(fid,1,0); <span class="comment">% 16-bit alignment</span>
0295 <span class="keyword">if</span> byte_align; fseek(fid,2,0); <span class="keyword">end</span> <span class="comment">% 32-bit alignment</span>
0296 se_hdr = setfield(se_hdr, <span class="string">'end_loc'</span>, fread(fid,1,<span class="string">'float32'</span>));            <span class="comment">%Last scan location (L/S)%</span>
0297 se_hdr = setfield(se_hdr, <span class="string">'se_pseq '</span>, fread(fid,1,<span class="string">'int16'</span>));            <span class="comment">%Last Pulse Sequence Used (L/S)%</span>
0298 se_hdr = setfield(se_hdr, <span class="string">'se_sortorder '</span>, fread(fid,1,<span class="string">'int16'</span>));       <span class="comment">%Image Sort Order (L/S)%</span>
0299 se_hdr = setfield(se_hdr, <span class="string">'se_lndmrkcnt'</span>, fread(fid,1,<span class="string">'int32'</span>));       <span class="comment">%Landmark Counter%</span>
0300 se_hdr = setfield(se_hdr, <span class="string">'se_nacq '</span>, fread(fid,1,<span class="string">'int16'</span>));            <span class="comment">%Number of Acquisitions%</span>
0301 se_hdr = setfield(se_hdr, <span class="string">'xbasest '</span>, fread(fid,1,<span class="string">'int16'</span>));            <span class="comment">%Starting number for baselines%</span>
0302 se_hdr = setfield(se_hdr, <span class="string">'xbaseend'</span>, fread(fid,1,<span class="string">'int16'</span>));            <span class="comment">%Ending number for baselines%</span>
0303 se_hdr = setfield(se_hdr, <span class="string">'xenhst'</span>, fread(fid,1,<span class="string">'int16'</span>));             <span class="comment">%Starting number for enhanced scans%</span>
0304 se_hdr = setfield(se_hdr, <span class="string">'xenhend'</span>, fread(fid,1,<span class="string">'int16'</span>));            <span class="comment">%Ending number for enhanced scans%</span>
0305 <span class="keyword">if</span> byte_align; fseek(fid,2,0); <span class="keyword">end</span> <span class="comment">% 32-bit alignment</span>
0306 se_hdr = setfield(se_hdr, <span class="string">'se_lastmod'</span>, fread(fid,1,<span class="string">'int32'</span>));         <span class="comment">%Date/Time of Last Change%</span>
0307 se_hdr = setfield(se_hdr, <span class="string">'se_alloc_key'</span>, fread(fid,13,<span class="string">'uchar'</span>));  <span class="comment">%Process that allocated this record%</span>
0308 fseek(fid,1,0); <span class="comment">% 16-bit alignment</span>
0309 <span class="keyword">if</span> byte_align; fseek(fid,2,0); <span class="keyword">end</span> <span class="comment">% 32-bit alignment</span>
0310 se_hdr = setfield(se_hdr, <span class="string">'se_delta_cnt'</span>, fread(fid,1,<span class="string">'int32'</span>));       <span class="comment">%Indicates number of updates to header%</span>
0311 se_hdr = setfield(se_hdr, <span class="string">'se_verscre'</span>, fread(fid,2,<span class="string">'uchar'</span>));     <span class="comment">%Genesis Version - Created%</span>
0312 se_hdr = setfield(se_hdr, <span class="string">'se_verscur'</span>, fread(fid,2,<span class="string">'uchar'</span>));     <span class="comment">%Genesis Version - Now%</span>
0313 se_hdr = setfield(se_hdr, <span class="string">'se_pds_a'</span>, fread(fid,1,<span class="string">'float32'</span>));           <span class="comment">%PixelData size - as stored%</span>
0314 se_hdr = setfield(se_hdr, <span class="string">'se_pds_c'</span>, fread(fid,1,<span class="string">'float32'</span>));           <span class="comment">%PixelData size - Compressed%</span>
0315 se_hdr = setfield(se_hdr, <span class="string">'se_pds_u'</span>, fread(fid,1,<span class="string">'float32'</span>));           <span class="comment">%PixelData size - UnCompressed%</span>
0316 se_hdr = setfield(se_hdr, <span class="string">'se_checksum'</span>, fread(fid,1,<span class="string">'uint32'</span>));        <span class="comment">%Series Record checksum%</span>
0317 se_hdr = setfield(se_hdr, <span class="string">'se_complete'</span>, fread(fid,1,<span class="string">'int32'</span>));        <span class="comment">%Series Complete Flag%</span>
0318 se_hdr = setfield(se_hdr, <span class="string">'se_numarch'</span>, fread(fid,1,<span class="string">'int32'</span>));         <span class="comment">%Number of Images Archived%</span>
0319 se_hdr = setfield(se_hdr, <span class="string">'se_imagect'</span>, fread(fid,1,<span class="string">'int32'</span>));         <span class="comment">%Last Image Number Used%</span>
0320 se_hdr = setfield(se_hdr, <span class="string">'se_numimages'</span>, fread(fid,1,<span class="string">'int32'</span>));       <span class="comment">%Number of Images Existing%</span>
0321 se_hdr = setfield(se_hdr, <span class="string">'se_images'</span>, struct(<span class="string">'length'</span>, fread(fid,1,<span class="string">'uint32'</span>), <span class="keyword">...</span>
0322     <span class="string">'data'</span>, fread(fid,1,<span class="string">'uint32'</span>))); <span class="comment">%Image Keys for this Series%</span>
0323 se_hdr = setfield(se_hdr, <span class="string">'se_numunimg'</span>, fread(fid,1,<span class="string">'int32'</span>));        <span class="comment">%Number of Unstored Images%</span>
0324 se_hdr = setfield(se_hdr, <span class="string">'se_unimages'</span>, struct(<span class="string">'length'</span>, fread(fid,1,<span class="string">'uint32'</span>), <span class="keyword">...</span>
0325     <span class="string">'data'</span>, fread(fid,1,<span class="string">'uint32'</span>))); <span class="comment">%Unstored Image Keys for this Series%</span>
0326 se_hdr = setfield(se_hdr, <span class="string">'se_toarchcnt'</span>, fread(fid,1,<span class="string">'int32'</span>));       <span class="comment">%Number of Unarchived Images%</span>
0327 se_hdr = setfield(se_hdr, <span class="string">'se_toarchive'</span>, struct(<span class="string">'length'</span>, fread(fid,1,<span class="string">'uint32'</span>), <span class="keyword">...</span>
0328     <span class="string">'data'</span>, fread(fid,1,<span class="string">'uint32'</span>))); <span class="comment">%Unarchived Image Keys for this Series%</span>
0329 se_hdr = setfield(se_hdr, <span class="string">'echo1_alpha'</span>, fread(fid,1,<span class="string">'float32'</span>));        <span class="comment">%Echo 1 Alpha Value%</span>
0330 se_hdr = setfield(se_hdr, <span class="string">'echo1_beta'</span>, fread(fid,1,<span class="string">'float32'</span>));         <span class="comment">%Echo 1 Beta Value%</span>
0331 se_hdr = setfield(se_hdr, <span class="string">'echo1_window'</span>, fread(fid,1,<span class="string">'uint16'</span>));       <span class="comment">%Echo 1 Window Value%</span>
0332 se_hdr = setfield(se_hdr, <span class="string">'echo1_level'</span>, fread(fid,1,<span class="string">'int16'</span>));        <span class="comment">%Echo 1 Level Value%</span>
0333 se_hdr = setfield(se_hdr, <span class="string">'echo2_alpha'</span>, fread(fid,1,<span class="string">'float32'</span>));        <span class="comment">%Echo 2 Alpha Value%</span>
0334 se_hdr = setfield(se_hdr, <span class="string">'echo2_beta'</span>, fread(fid,1,<span class="string">'float32'</span>));         <span class="comment">%Echo 2 Beta Value%</span>
0335 se_hdr = setfield(se_hdr, <span class="string">'echo2_window'</span>, fread(fid,1,<span class="string">'uint16'</span>));       <span class="comment">%Echo 2 Window Value%</span>
0336 se_hdr = setfield(se_hdr, <span class="string">'echo2_level'</span>, fread(fid,1,<span class="string">'int16'</span>));        <span class="comment">%Echo 2 Level Value%</span>
0337 se_hdr = setfield(se_hdr, <span class="string">'echo3_alpha'</span>, fread(fid,1,<span class="string">'float32'</span>));        <span class="comment">%Echo 3 Alpha Value%</span>
0338 se_hdr = setfield(se_hdr, <span class="string">'echo3_beta'</span>, fread(fid,1,<span class="string">'float32'</span>));         <span class="comment">%Echo 3 Beta Value%</span>
0339 se_hdr = setfield(se_hdr, <span class="string">'echo3_window'</span>, fread(fid,1,<span class="string">'uint16'</span>));       <span class="comment">%Echo 3 Window Value%</span>
0340 se_hdr = setfield(se_hdr, <span class="string">'echo3_level'</span>, fread(fid,1,<span class="string">'int16'</span>));        <span class="comment">%Echo 3 Level Value%</span>
0341 se_hdr = setfield(se_hdr, <span class="string">'echo4_alpha'</span>, fread(fid,1,<span class="string">'float32'</span>));        <span class="comment">%Echo 4 Alpha Value%</span>
0342 se_hdr = setfield(se_hdr, <span class="string">'echo4_beta'</span>, fread(fid,1,<span class="string">'float32'</span>));         <span class="comment">%Echo 4 Beta Value%</span>
0343 se_hdr = setfield(se_hdr, <span class="string">'echo4_window'</span>, fread(fid,1,<span class="string">'uint16'</span>));       <span class="comment">%Echo 4 Window Value%</span>
0344 se_hdr = setfield(se_hdr, <span class="string">'echo4_level'</span>, fread(fid,1,<span class="string">'int16'</span>));        <span class="comment">%Echo 4 Level Value%</span>
0345 se_hdr = setfield(se_hdr, <span class="string">'echo5_alpha'</span>, fread(fid,1,<span class="string">'float32'</span>));        <span class="comment">%Echo 5 Alpha Value%</span>
0346 se_hdr = setfield(se_hdr, <span class="string">'echo5_beta'</span>, fread(fid,1,<span class="string">'float32'</span>));         <span class="comment">%Echo 5 Beta Value%</span>
0347 se_hdr = setfield(se_hdr, <span class="string">'echo5_window'</span>, fread(fid,1,<span class="string">'uint16'</span>));       <span class="comment">%Echo 5 Window Value%</span>
0348 se_hdr = setfield(se_hdr, <span class="string">'echo5_level'</span>, fread(fid,1,<span class="string">'int16'</span>));        <span class="comment">%Echo 5 Level Value%</span>
0349 se_hdr = setfield(se_hdr, <span class="string">'echo6_alpha'</span>, fread(fid,1,<span class="string">'float32'</span>));        <span class="comment">%Echo 6 Alpha Value%</span>
0350 se_hdr = setfield(se_hdr, <span class="string">'echo6_beta'</span>, fread(fid,1,<span class="string">'float32'</span>));         <span class="comment">%Echo 6 Beta Value%</span>
0351 se_hdr = setfield(se_hdr, <span class="string">'echo6_window'</span>, fread(fid,1,<span class="string">'uint16'</span>));       <span class="comment">%Echo 6 Window Value%</span>
0352 se_hdr = setfield(se_hdr, <span class="string">'echo6_level'</span>, fread(fid,1,<span class="string">'int16'</span>));        <span class="comment">%Echo 6 Level Value%</span>
0353 se_hdr = setfield(se_hdr, <span class="string">'echo7_alpha'</span>, fread(fid,1,<span class="string">'float32'</span>));        <span class="comment">%Echo 7 Alpha Value%</span>
0354 se_hdr = setfield(se_hdr, <span class="string">'echo7_beta'</span>, fread(fid,1,<span class="string">'float32'</span>));         <span class="comment">%Echo 7 Beta Value%</span>
0355 se_hdr = setfield(se_hdr, <span class="string">'echo7_window'</span>, fread(fid,1,<span class="string">'uint16'</span>));       <span class="comment">%Echo 7 Window Value%</span>
0356 se_hdr = setfield(se_hdr, <span class="string">'echo7_level'</span>, fread(fid,1,<span class="string">'int16'</span>));        <span class="comment">%Echo 7 Level Value%</span>
0357 se_hdr = setfield(se_hdr, <span class="string">'echo8_alpha'</span>, fread(fid,1,<span class="string">'float32'</span>));        <span class="comment">%Echo 8 Alpha Value%</span>
0358 se_hdr = setfield(se_hdr, <span class="string">'echo8_beta'</span>, fread(fid,1,<span class="string">'float32'</span>));         <span class="comment">%Echo 8 Beta Value%</span>
0359 se_hdr = setfield(se_hdr, <span class="string">'echo8_window'</span>, fread(fid,1,<span class="string">'uint16'</span>));       <span class="comment">%Echo 8 Window Value%</span>
0360 se_hdr = setfield(se_hdr, <span class="string">'echo8_level'</span>, fread(fid,1,<span class="string">'int16'</span>));        <span class="comment">%Echo 8 Level Value%</span>
0361 se_hdr = setfield(se_hdr, <span class="string">'series_uid'</span>, fread(fid,32,<span class="string">'uchar'</span>));    <span class="comment">%Series Entity Unique ID%</span>
0362 se_hdr = setfield(se_hdr, <span class="string">'landmark_uid'</span>, fread(fid,32,<span class="string">'uchar'</span>));  <span class="comment">%Landmark Unique ID%</span>
0363 se_hdr = setfield(se_hdr, <span class="string">'equipmnt_uid'</span>, fread(fid,32,<span class="string">'uchar'</span>));  <span class="comment">%Equipment Unique ID%</span>
0364 se_hdr = setfield(se_hdr, <span class="string">'se_padding'</span>, fread(fid,588,<span class="string">'uchar'</span>));   <span class="comment">%Spare Space%</span>
0365 
0366 <span class="keyword">return</span>
0367 
0368 
0369 
0370 
0371 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span>
0372 
0373 <a name="_sub4" href="#_subfunctions" class="code">function im_hdr = ge_readHeaderImage(fid, byte_align)</a>
0374 
0375 <span class="comment">% returns the image header as a structure.</span>
0376 <span class="comment">% if byte_align = 1 then 32-bit alignment (SGI, LX2 format)</span>
0377 <span class="comment">% if byte_align = 0 then 16-bit alignment (Sun, 5.X format)</span>
0378 
0379 <span class="comment">% define the structure and read in the data</span>
0380 <span class="comment">% to overcome the byte alignment problems</span>
0381 <span class="comment">% break up the assignment into pieces using the setfield function</span>