ge_hdr_read.html

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

0382 
0383 im_hdr = struct(<span class="string">'im_suid'</span>, fread(fid,4,<span class="string">'uchar'</span>)); <span class="comment">%Suite id for this image</span>
0384 im_hdr = setfield(im_hdr, <span class="string">'im_uniq'</span>, fread(fid,1,<span class="string">'int16'</span>));            <span class="comment">%The Make-Unique Flag</span>
0385 im_hdr = setfield(im_hdr, <span class="string">'im_diskid'</span>, fread(fid,1,<span class="string">'uchar'</span>));          <span class="comment">%Disk ID for this Image</span>
0386 fseek(fid, 1, 0);<span class="comment">% 16-bit alignment</span>
0387 im_hdr = setfield(im_hdr, <span class="string">'im_exno'</span>, fread(fid,1,<span class="string">'uint16'</span>));            <span class="comment">%Exam number for this image</span>
0388 im_hdr = setfield(im_hdr, <span class="string">'im_seno'</span>, fread(fid,1,<span class="string">'int16'</span>));            <span class="comment">%Series Number for this image</span>
0389 im_hdr = setfield(im_hdr, <span class="string">'im_no'</span>, fread(fid,1,<span class="string">'int16'</span>));              <span class="comment">%Image Number</span>
0390 <span class="keyword">if</span> byte_align; fseek(fid, 2, 0); <span class="keyword">end</span> <span class="comment">% 32-bit alignment</span>
0391 im_hdr = setfield(im_hdr, <span class="string">'im_datetime'</span>, fread(fid,1,<span class="string">'int32'</span>));        <span class="comment">%Allocation Image date/time stamp</span>
0392 im_hdr = setfield(im_hdr, <span class="string">'im_actual_dt'</span>, fread(fid,1,<span class="string">'int32'</span>));       <span class="comment">%Actual Image date/time stamp</span>
0393 im_hdr = setfield(im_hdr, <span class="string">'sctime'</span>, fread(fid,1,<span class="string">'float32'</span>));             <span class="comment">%Duration of scan</span>
0394 im_hdr = setfield(im_hdr, <span class="string">'slthick'</span>, fread(fid,1,<span class="string">'float32'</span>));            <span class="comment">%Slice Thickness (mm)</span>
0395 im_hdr = setfield(im_hdr, <span class="string">'imatrix_X'</span>, fread(fid,1,<span class="string">'int16'</span>));          <span class="comment">%Image matrix size - X</span>
0396 im_hdr = setfield(im_hdr, <span class="string">'imatrix_Y'</span>, fread(fid,1,<span class="string">'int16'</span>));          <span class="comment">%Image matrix size - Y</span>
0397 im_hdr = setfield(im_hdr, <span class="string">'dfov'</span>, fread(fid,1,<span class="string">'float32'</span>));               <span class="comment">%Display field of view - X (mm)</span>
0398 im_hdr = setfield(im_hdr, <span class="string">'dfov_rect'</span>, fread(fid,1,<span class="string">'float32'</span>));          <span class="comment">%Display field of view - Y (if different)</span>
0399 im_hdr = setfield(im_hdr, <span class="string">'dim_X'</span>, fread(fid,1,<span class="string">'float32'</span>));              <span class="comment">%Image dimension - X</span>
0400 im_hdr = setfield(im_hdr, <span class="string">'dim_Y'</span>, fread(fid,1,<span class="string">'float32'</span>));              <span class="comment">%Image dimension - Y</span>
0401 im_hdr = setfield(im_hdr, <span class="string">'pixsize_X'</span>, fread(fid,1,<span class="string">'float32'</span>));          <span class="comment">%Image pixel size - X</span>
0402 im_hdr = setfield(im_hdr, <span class="string">'pixsize_Y'</span>, fread(fid,1,<span class="string">'float32'</span>));          <span class="comment">%Image pixel size - Y</span>
0403 im_hdr = setfield(im_hdr, <span class="string">'pdid'</span>, fread(fid,14,<span class="string">'uchar'</span>));          <span class="comment">%Pixel Data ID</span>
0404 im_hdr = setfield(im_hdr, <span class="string">'contrastIV'</span>, fread(fid,17,<span class="string">'uchar'</span>));    <span class="comment">%IV Contrast Agent</span>
0405 im_hdr = setfield(im_hdr, <span class="string">'contrastOral'</span>, fread(fid,17,<span class="string">'uchar'</span>));  <span class="comment">%Oral Contrast Agent</span>
0406 im_hdr = setfield(im_hdr, <span class="string">'contmode'</span>, fread(fid,1,<span class="string">'int16'</span>));           <span class="comment">%Image Contrast Mode</span>
0407 im_hdr = setfield(im_hdr, <span class="string">'serrx'</span>, fread(fid,1,<span class="string">'int16'</span>));              <span class="comment">%Series from which prescribed</span>
0408 im_hdr = setfield(im_hdr, <span class="string">'imgrx'</span>, fread(fid,1,<span class="string">'int16'</span>));              <span class="comment">%Image from which prescribed</span>
0409 im_hdr = setfield(im_hdr, <span class="string">'screenformat'</span>, fread(fid,1,<span class="string">'int16'</span>));       <span class="comment">%Screen Format(8/16 bit)</span>
0410 im_hdr = setfield(im_hdr, <span class="string">'plane'</span>, fread(fid,1,<span class="string">'int16'</span>));              <span class="comment">%Plane Type</span>
0411 <span class="keyword">if</span> byte_align; fseek(fid, 2, 0); <span class="keyword">end</span> <span class="comment">% 32-bit alignment</span>
0412 im_hdr = setfield(im_hdr, <span class="string">'scanspacing'</span>, fread(fid,1,<span class="string">'float32'</span>));        <span class="comment">%Spacing between scans (mm?)</span>
0413 im_hdr = setfield(im_hdr, <span class="string">'im_compress'</span>, fread(fid,1,<span class="string">'int16'</span>));        <span class="comment">%Image compression type for allocation</span>
0414 im_hdr = setfield(im_hdr, <span class="string">'im_scouttype'</span>, fread(fid,1,<span class="string">'int16'</span>));       <span class="comment">%Scout Type (AP or lateral)</span>
0415 im_hdr = setfield(im_hdr, <span class="string">'loc_ras'</span>, fread(fid,1,<span class="string">'uchar'</span>));            <span class="comment">%RAS letter of image location</span>
0416 fseek(fid, 1, 0); <span class="comment">% 16-bit alignment</span>
0417 <span class="keyword">if</span> byte_align; fseek(fid, 2, 0); <span class="keyword">end</span> <span class="comment">% 32-bit alignment</span>
0418 im_hdr = setfield(im_hdr, <span class="string">'loc'</span>, fread(fid,1,<span class="string">'float32'</span>));                <span class="comment">%Image location</span>
0419 im_hdr = setfield(im_hdr, <span class="string">'ctr_R'</span>, fread(fid,1,<span class="string">'float32'</span>));              <span class="comment">%Center R coord of plane image</span>
0420 im_hdr = setfield(im_hdr, <span class="string">'ctr_A'</span>, fread(fid,1,<span class="string">'float32'</span>));              <span class="comment">%Center A coord of plane image</span>
0421 im_hdr = setfield(im_hdr, <span class="string">'ctr_S'</span>, fread(fid,1,<span class="string">'float32'</span>));              <span class="comment">%Center S coord of plane image</span>
0422 im_hdr = setfield(im_hdr, <span class="string">'norm_R'</span>, fread(fid,1,<span class="string">'float32'</span>));             <span class="comment">%Normal R coord</span>
0423 im_hdr = setfield(im_hdr, <span class="string">'norm_A'</span>, fread(fid,1,<span class="string">'float32'</span>));             <span class="comment">%Normal A coord</span>
0424 im_hdr = setfield(im_hdr, <span class="string">'norm_S'</span>, fread(fid,1,<span class="string">'float32'</span>));             <span class="comment">%Normal S coord</span>
0425 im_hdr = setfield(im_hdr, <span class="string">'tlhc_R'</span>, fread(fid,1,<span class="string">'float32'</span>));             <span class="comment">%R Coord of Top Left Hand Corner</span>
0426 im_hdr = setfield(im_hdr, <span class="string">'tlhc_A'</span>, fread(fid,1,<span class="string">'float32'</span>));             <span class="comment">%A Coord of Top Left Hand Corner</span>
0427 im_hdr = setfield(im_hdr, <span class="string">'tlhc_S'</span>, fread(fid,1,<span class="string">'float32'</span>));             <span class="comment">%S Coord of Top Left Hand Corner</span>
0428 im_hdr = setfield(im_hdr, <span class="string">'trhc_R'</span>, fread(fid,1,<span class="string">'float32'</span>));             <span class="comment">%R Coord of Top Right Hand Corner</span>
0429 im_hdr = setfield(im_hdr, <span class="string">'trhc_A'</span>, fread(fid,1,<span class="string">'float32'</span>));             <span class="comment">%A Coord of Top Right Hand Corner</span>
0430 im_hdr = setfield(im_hdr, <span class="string">'trhc_S'</span>, fread(fid,1,<span class="string">'float32'</span>));             <span class="comment">%S Coord of Top Right Hand Corner</span>
0431 im_hdr = setfield(im_hdr, <span class="string">'brhc_R'</span>, fread(fid,1,<span class="string">'float32'</span>));             <span class="comment">%R Coord of Bottom Right Hand Corner</span>
0432 im_hdr = setfield(im_hdr, <span class="string">'brhc_A'</span>, fread(fid,1,<span class="string">'float32'</span>));             <span class="comment">%A Coord of Bottom Right Hand Corner</span>
0433 im_hdr = setfield(im_hdr, <span class="string">'brhc_S'</span>, fread(fid,1,<span class="string">'float32'</span>));             <span class="comment">%S Coord of Bottom Right Hand Corner</span>
0434 im_hdr = setfield(im_hdr, <span class="string">'forimgrev'</span>, fread(fid,4,<span class="string">'uchar'</span>));      <span class="comment">%Foreign Image Revision</span>
0435 im_hdr = setfield(im_hdr, <span class="string">'tr'</span>, fread(fid,1,<span class="string">'int32'</span>));                 <span class="comment">%Pulse repetition time(usec)</span>
0436 im_hdr = setfield(im_hdr, <span class="string">'ti'</span>, fread(fid,1,<span class="string">'int32'</span>));                 <span class="comment">%Pulse inversion time(usec)</span>
0437 im_hdr = setfield(im_hdr, <span class="string">'te'</span>, fread(fid,1,<span class="string">'int32'</span>));                 <span class="comment">%Pulse echo time(usec)</span>
0438 im_hdr = setfield(im_hdr, <span class="string">'te2'</span>, fread(fid,1,<span class="string">'int32'</span>));                <span class="comment">%Second echo echo (usec)</span>
0439 im_hdr = setfield(im_hdr, <span class="string">'numecho'</span>, fread(fid,1,<span class="string">'int16'</span>));            <span class="comment">%Number of echoes</span>
0440 im_hdr = setfield(im_hdr, <span class="string">'echonum'</span>, fread(fid,1,<span class="string">'int16'</span>));            <span class="comment">%Echo Number</span>
0441 im_hdr = setfield(im_hdr, <span class="string">'tbldlta'</span>, fread(fid,1,<span class="string">'float32'</span>));            <span class="comment">%Table Delta</span>
0442 im_hdr = setfield(im_hdr, <span class="string">'nex'</span>, fread(fid,1,<span class="string">'float32'</span>));                <span class="comment">%Number of Excitations</span>
0443 im_hdr = setfield(im_hdr, <span class="string">'contig'</span>, fread(fid,1,<span class="string">'int16'</span>));             <span class="comment">%Continuous Slices Flag</span>
0444 im_hdr = setfield(im_hdr, <span class="string">'hrtrate'</span>, fread(fid,1,<span class="string">'int16'</span>));            <span class="comment">%Cardiac Heart Rate (bpm)</span>
0445 im_hdr = setfield(im_hdr, <span class="string">'tdel'</span>, fread(fid,1,<span class="string">'int32'</span>));               <span class="comment">%Delay time after trigger (msec)</span>
0446 im_hdr = setfield(im_hdr, <span class="string">'saravg'</span>, fread(fid,1,<span class="string">'float32'</span>));             <span class="comment">%Average SAR</span>
0447 im_hdr = setfield(im_hdr, <span class="string">'sarpeak'</span>, fread(fid,1,<span class="string">'float32'</span>));            <span class="comment">%Peak SAR</span>
0448 im_hdr = setfield(im_hdr, <span class="string">'monsar'</span>, fread(fid,1,<span class="string">'int16'</span>));             <span class="comment">%Monitor SAR flag</span>
0449 im_hdr = setfield(im_hdr, <span class="string">'trgwindow'</span>, fread(fid,1,<span class="string">'int16'</span>));          <span class="comment">%Trigger window (% of R-R interval)</span>
0450 im_hdr = setfield(im_hdr, <span class="string">'reptime'</span>, fread(fid,1,<span class="string">'float32'</span>));            <span class="comment">%Cardiac repetition time</span>
0451 im_hdr = setfield(im_hdr, <span class="string">'imgpcyc'</span>, fread(fid,1,<span class="string">'int16'</span>));            <span class="comment">%Images per cardiac cycle</span>
0452 im_hdr = setfield(im_hdr, <span class="string">'xmtgain'</span>, fread(fid,1,<span class="string">'int16'</span>));            <span class="comment">%Actual Transmit Gain (.1 db)</span>
0453 im_hdr = setfield(im_hdr, <span class="string">'rcvgain1'</span>, fread(fid,1,<span class="string">'int16'</span>));           <span class="comment">%Actual Receive Gain Analog (.1 db)</span>
0454 im_hdr = setfield(im_hdr, <span class="string">'rcvgain2'</span>, fread(fid,1,<span class="string">'int16'</span>));           <span class="comment">%Actual Receive Gain Digital (.1 db)</span>
0455 im_hdr = setfield(im_hdr, <span class="string">'mr_flip'</span>, fread(fid,1,<span class="string">'int16'</span>));            <span class="comment">%Flip Angle for GRASS scans (deg.)</span>
0456 <span class="keyword">if</span> byte_align; fseek(fid, 2, 0); <span class="keyword">end</span> <span class="comment">% byte alignment</span>
0457 im_hdr = setfield(im_hdr, <span class="string">'mindat'</span>, fread(fid,1,<span class="string">'int32'</span>));             <span class="comment">%Minimum Delay after Trigger (uSec)</span>
0458 im_hdr = setfield(im_hdr, <span class="string">'cphase'</span>, fread(fid,1,<span class="string">'int16'</span>));             <span class="comment">%Total Cardiac Phase prescribed</span>
0459 im_hdr = setfield(im_hdr, <span class="string">'swappf'</span>, fread(fid,1,<span class="string">'int16'</span>));             <span class="comment">%Swap Phase/Frequency Axis</span>
0460 im_hdr = setfield(im_hdr, <span class="string">'pauseint'</span>, fread(fid,1,<span class="string">'int16'</span>));           <span class="comment">%Pause Interval (slices)</span>
0461 <span class="keyword">if</span> byte_align; fseek(fid, 2, 0); <span class="keyword">end</span> <span class="comment">% 32-bit alignment</span>
0462 im_hdr = setfield(im_hdr, <span class="string">'pausetime'</span>, fread(fid,1,<span class="string">'float32'</span>));          <span class="comment">%Pause Time</span>
0463 im_hdr = setfield(im_hdr, <span class="string">'obplane'</span>, fread(fid,1,<span class="string">'int32'</span>));            <span class="comment">%Oblique Plane</span>
0464 im_hdr = setfield(im_hdr, <span class="string">'slocfov'</span>, fread(fid,1,<span class="string">'int32'</span>));            <span class="comment">%Slice Offsets on Freq axis</span>
0465 im_hdr = setfield(im_hdr, <span class="string">'xmtfreq'</span>, fread(fid,1,<span class="string">'int32'</span>));            <span class="comment">%Center Frequency (0.1 Hz)</span>
0466 im_hdr = setfield(im_hdr, <span class="string">'autoxmtfreq'</span>, fread(fid,1,<span class="string">'int32'</span>));        <span class="comment">%Auto Center Frequency (0.1 Hz)</span>
0467 im_hdr = setfield(im_hdr, <span class="string">'autoxmtgain'</span>, fread(fid,1,<span class="string">'int16'</span>));        <span class="comment">%Auto Transmit Gain (0.1 dB)</span>
0468 im_hdr = setfield(im_hdr, <span class="string">'prescan_r1'</span>, fread(fid,1,<span class="string">'int16'</span>));         <span class="comment">%PreScan R1 - Analog</span>
0469 im_hdr = setfield(im_hdr, <span class="string">'prescan_r2'</span>, fread(fid,1,<span class="string">'int16'</span>));         <span class="comment">%PreScan R2 - Digital</span>
0470 <span class="keyword">if</span> byte_align; fseek(fid, 2, 0); <span class="keyword">end</span> <span class="comment">% 32-bit alignment</span>
0471 im_hdr = setfield(im_hdr, <span class="string">'user_bitmap'</span>, fread(fid,1,<span class="string">'int32'</span>));        <span class="comment">%Bitmap defining user CVs</span>
0472 im_hdr = setfield(im_hdr, <span class="string">'cenfreq'</span>, fread(fid,1,<span class="string">'int16'</span>));            <span class="comment">%Center Frequency Method</span>
0473 im_hdr = setfield(im_hdr, <span class="string">'imode'</span>, fread(fid,1,<span class="string">'int16'</span>));              <span class="comment">%Imaging Mode</span>
0474 im_hdr = setfield(im_hdr, <span class="string">'iopt'</span>, fread(fid,1,<span class="string">'int32'</span>));               <span class="comment">%Imaging Options</span>
0475 im_hdr = setfield(im_hdr, <span class="string">'pseq'</span>, fread(fid,1,<span class="string">'int16'</span>));               <span class="comment">%Pulse Sequence</span>
0476 im_hdr = setfield(im_hdr, <span class="string">'pseqmode'</span>, fread(fid,1,<span class="string">'int16'</span>));           <span class="comment">%Pulse Sequence Mode</span>
0477 im_hdr = setfield(im_hdr, <span class="string">'psdname'</span>, fread(fid,33,<span class="string">'uchar'</span>));       <span class="comment">%Pulse Sequence Name</span>
0478 fseek(fid,1,0); <span class="comment">% 16-bit alignment</span>
0479 <span class="keyword">if</span> byte_align; fseek(fid, 2, 0); <span class="keyword">end</span> <span class="comment">% byte alignment</span>
0480 im_hdr = setfield(im_hdr, <span class="string">'psd_datetime'</span>, fread(fid,1,<span class="string">'int32'</span>));       <span class="comment">%PSD Creation Date and Time</span>
0481 im_hdr = setfield(im_hdr, <span class="string">'psd_iname'</span>, fread(fid,13,<span class="string">'uchar'</span>));     <span class="comment">%PSD name from inside PSD</span>
0482 fseek(fid, 1, 0); <span class="comment">% 16-bit alignment</span>
0483 im_hdr = setfield(im_hdr, <span class="string">'ctyp'</span>, fread(fid,1,<span class="string">'int16'</span>));               <span class="comment">%Coil Type</span>
0484 im_hdr = setfield(im_hdr, <span class="string">'cname'</span>, fread(fid,17,<span class="string">'uchar'</span>));         <span class="comment">%Coil Name</span>
0485 fseek(fid, 1, 0); <span class="comment">% 16-bit alignment</span>
0486 im_hdr = setfield(im_hdr, <span class="string">'surfctyp'</span>, fread(fid,1,<span class="string">'int16'</span>));           <span class="comment">%Surface Coil Type</span>
0487 im_hdr = setfield(im_hdr, <span class="string">'surfcext'</span>, fread(fid,1,<span class="string">'int16'</span>));           <span class="comment">%Extremity Coil Flag</span>
0488 <span class="keyword">if</span> byte_align; fseek(fid, 2, 0); <span class="keyword">end</span> <span class="comment">% 32-bit alignment</span>
0489 im_hdr = setfield(im_hdr, <span class="string">'rawrunnum'</