00001 function se_hdr = GE_readHeaderSeries(fid, byte_align)
00002 %
00003 % se_hdr = read_series_header(fid, byte_align)
00004 %
00005 % Loads the series header from a file with filed id fid
00006 % and returns it as a structure.
00007 % if byte_align = 1 then 32-bit alignment (SGI, LX2 format)
00008 % if byte_align = 0 then 16-bit alignment (Sun, 5.X format)
00009 %
00010 %
00011 % Souheil J. Inati
00012 % Dartmouth College
00013 % May 2000
00014 % souheil.inati@dartmouth.edu
00015 %
00016
00017 % define the structure and read in the data
00018 % to overcome the byte alignment problems
00019 % break up the assignment into pieces using the setfield function
00020 se_hdr = struct('se_suid', fread(fid,4,'uchar')); %Suite ID for this Series%
00021 se_hdr = setfield(se_hdr, 'se_uniq', fread(fid,1,'int16')); %The Make-Unique Flag%
00022 se_hdr = setfield(se_hdr, 'se_diskid', fread(fid,1,'uchar')); %Disk ID for this Series%
00023 fseek(fid,1,0); % 16-bit alignment
00024 se_hdr = setfield(se_hdr, 'se_exno', fread(fid,1,'uint16')); %Exam Number%
00025 se_hdr = setfield(se_hdr, 'se_no ', fread(fid,1,'int16')); %Series Number%
00026 se_hdr = setfield(se_hdr, 'se_datetime', fread(fid,1,'int32')); %Allocation Series Data/Time stamp%
00027 se_hdr = setfield(se_hdr, 'se_actual_dt', fread(fid,1,'int32')); %Actual Series Data/Time stamp%
00028 se_hdr = setfield(se_hdr, 'se_desc', fread(fid,30,'uchar')); %Series Description%
00029 se_hdr = setfield(se_hdr, 'pr_sysid', fread(fid,9,'uchar')); %Primary Receiver Suite and Host%
00030 se_hdr = setfield(se_hdr, 'pansysid', fread(fid,9,'uchar')); %Archiver Suite and Host%
00031 se_hdr = setfield(se_hdr, 'se_typ ', fread(fid,1,'int16')); %Series Type%
00032 se_hdr = setfield(se_hdr, 'se_source ', fread(fid,1,'int16')); %Series from which prescribed%
00033 se_hdr = setfield(se_hdr, 'se_plane ', fread(fid,1,'int16')); %Most-like Plane (for L/S)%
00034 se_hdr = setfield(se_hdr, 'scan_type ', fread(fid,1,'int16')); %Scout or Axial (for CT)%
00035 se_hdr = setfield(se_hdr, 'position', fread(fid,1,'int32')); %Patient Position%
00036 se_hdr = setfield(se_hdr, 'entry', fread(fid,1,'int32')); %Patient Entry%
00037 se_hdr = setfield(se_hdr, 'anref', fread(fid,3,'uchar')); %Anatomical reference%
00038 fseek(fid,1,0); % 16-bit alignment
00039 se_hdr = setfield(se_hdr, 'lmhor', fread(fid,1,'float32')); %Horizontal Landmark%
00040 se_hdr = setfield(se_hdr, 'prtcl', fread(fid,25,'uchar')); %Scan Protocol Name%
00041 fseek(fid,1,0); % 16-bit alignment
00042 se_hdr = setfield(se_hdr, 'se_contrast ', fread(fid,1,'int16')); %Non-zero if > 0 image used contrast(L/S)%
00043 se_hdr = setfield(se_hdr, 'start_ras', fread(fid,1,'uchar')); %RAS letter for first scan location (L/S)%
00044 fseek(fid,1,0); % 16-bit alignment
00045 if byte_align; fseek(fid,2,0); end % 32-bit alignment
00046 se_hdr = setfield(se_hdr, 'start_loc', fread(fid,1,'float32')); %First scan location (L/S)%
00047 se_hdr = setfield(se_hdr, 'end_ras', fread(fid,1,'uchar')); %RAS letter for last scan location (L/S)%
00048 fseek(fid,1,0); % 16-bit alignment
00049 if byte_align; fseek(fid,2,0); end % 32-bit alignment
00050 se_hdr = setfield(se_hdr, 'end_loc', fread(fid,1,'float32')); %Last scan location (L/S)%
00051 se_hdr = setfield(se_hdr, 'se_pseq ', fread(fid,1,'int16')); %Last Pulse Sequence Used (L/S)%
00052 se_hdr = setfield(se_hdr, 'se_sortorder ', fread(fid,1,'int16')); %Image Sort Order (L/S)%
00053 se_hdr = setfield(se_hdr, 'se_lndmrkcnt', fread(fid,1,'int32')); %Landmark Counter%
00054 se_hdr = setfield(se_hdr, 'se_nacq ', fread(fid,1,'int16')); %Number of Acquisitions%
00055 se_hdr = setfield(se_hdr, 'xbasest ', fread(fid,1,'int16')); %Starting number for baselines%
00056 se_hdr = setfield(se_hdr, 'xbaseend', fread(fid,1,'int16')); %Ending number for baselines%
00057 se_hdr = setfield(se_hdr, 'xenhst', fread(fid,1,'int16')); %Starting number for enhanced scans%
00058 se_hdr = setfield(se_hdr, 'xenhend', fread(fid,1,'int16')); %Ending number for enhanced scans%
00059 if byte_align; fseek(fid,2,0); end % 32-bit alignment
00060 se_hdr = setfield(se_hdr, 'se_lastmod', fread(fid,1,'int32')); %Date/Time of Last Change%
00061 se_hdr = setfield(se_hdr, 'se_alloc_key', fread(fid,13,'uchar')); %Process that allocated this record%
00062 fseek(fid,1,0); % 16-bit alignment
00063 if byte_align; fseek(fid,2,0); end % 32-bit alignment
00064 se_hdr = setfield(se_hdr, 'se_delta_cnt', fread(fid,1,'int32')); %Indicates number of updates to header%
00065 se_hdr = setfield(se_hdr, 'se_verscre', fread(fid,2,'uchar')); %Genesis Version - Created%
00066 se_hdr = setfield(se_hdr, 'se_verscur', fread(fid,2,'uchar')); %Genesis Version - Now%
00067 se_hdr = setfield(se_hdr, 'se_pds_a', fread(fid,1,'float32')); %PixelData size - as stored%
00068 se_hdr = setfield(se_hdr, 'se_pds_c', fread(fid,1,'float32')); %PixelData size - Compressed%
00069 se_hdr = setfield(se_hdr, 'se_pds_u', fread(fid,1,'float32')); %PixelData size - UnCompressed%
00070 se_hdr = setfield(se_hdr, 'se_checksum', fread(fid,1,'uint32')); %Series Record checksum%
00071 se_hdr = setfield(se_hdr, 'se_complete', fread(fid,1,'int32')); %Series Complete Flag%
00072 se_hdr = setfield(se_hdr, 'se_numarch', fread(fid,1,'int32')); %Number of Images Archived%
00073 se_hdr = setfield(se_hdr, 'se_imagect', fread(fid,1,'int32')); %Last Image Number Used%
00074 se_hdr = setfield(se_hdr, 'se_numimages', fread(fid,1,'int32')); %Number of Images Existing%
00075 se_hdr = setfield(se_hdr, 'se_images', struct('length', fread(fid,1,'uint32'), ...
00076 'data', fread(fid,1,'uint32'))); %Image Keys for this Series%
00077 se_hdr = setfield(se_hdr, 'se_numunimg', fread(fid,1,'int32')); %Number of Unstored Images%
00078 se_hdr = setfield(se_hdr, 'se_unimages', struct('length', fread(fid,1,'uint32'), ...
00079 'data', fread(fid,1,'uint32'))); %Unstored Image Keys for this Series%
00080 se_hdr = setfield(se_hdr, 'se_toarchcnt', fread(fid,1,'int32')); %Number of Unarchived Images%
00081 se_hdr = setfield(se_hdr, 'se_toarchive', struct('length', fread(fid,1,'uint32'), ...
00082 'data', fread(fid,1,'uint32'))); %Unarchived Image Keys for this Series%
00083 se_hdr = setfield(se_hdr, 'echo1_alpha', fread(fid,1,'float32')); %Echo 1 Alpha Value%
00084 se_hdr = setfield(se_hdr, 'echo1_beta', fread(fid,1,'float32')); %Echo 1 Beta Value%
00085 se_hdr = setfield(se_hdr, 'echo1_window', fread(fid,1,'uint16')); %Echo 1 Window Value%
00086 se_hdr = setfield(se_hdr, 'echo1_level', fread(fid,1,'int16')); %Echo 1 Level Value%
00087 se_hdr = setfield(se_hdr, 'echo2_alpha', fread(fid,1,'float32')); %Echo 2 Alpha Value%
00088 se_hdr = setfield(se_hdr, 'echo2_beta', fread(fid,1,'float32')); %Echo 2 Beta Value%
00089 se_hdr = setfield(se_hdr, 'echo2_window', fread(fid,1,'uint16')); %Echo 2 Window Value%
00090 se_hdr = setfield(se_hdr, 'echo2_level', fread(fid,1,'int16')); %Echo 2 Level Value%
00091 se_hdr = setfield(se_hdr, 'echo3_alpha', fread(fid,1,'float32')); %Echo 3 Alpha Value%
00092 se_hdr = setfield(se_hdr, 'echo3_beta', fread(fid,1,'float32')); %Echo 3 Beta Value%
00093 se_hdr = setfield(se_hdr, 'echo3_window', fread(fid,1,'uint16')); %Echo 3 Window Value%
00094 se_hdr = setfield(se_hdr, 'echo3_level', fread(fid,1,'int16')); %Echo 3 Level Value%
00095 se_hdr = setfield(se_hdr, 'echo4_alpha', fread(fid,1,'float32')); %Echo 4 Alpha Value%
00096 se_hdr = setfield(se_hdr, 'echo4_beta', fread(fid,1,'float32')); %Echo 4 Beta Value%
00097 se_hdr = setfield(se_hdr, 'echo4_window', fread(fid,1,'uint16')); %Echo 4 Window Value%
00098 se_hdr = setfield(se_hdr, 'echo4_level', fread(fid,1,'int16')); %Echo 4 Level Value%
00099 se_hdr = setfield(se_hdr, 'echo5_alpha', fread(fid,1,'float32')); %Echo 5 Alpha Value%
00100 se_hdr = setfield(se_hdr, 'echo5_beta', fread(fid,1,'float32')); %Echo 5 Beta Value%
00101 se_hdr = setfield(se_hdr, 'echo5_window', fread(fid,1,'uint16')); %Echo 5 Window Value%
00102 se_hdr = setfield(se_hdr, 'echo5_level', fread(fid,1,'int16')); %Echo 5 Level Value%
00103 se_hdr = setfield(se_hdr, 'echo6_alpha', fread(fid,1,'float32')); %Echo 6 Alpha Value%
00104 se_hdr = setfield(se_hdr, 'echo6_beta', fread(fid,1,'float32')); %Echo 6 Beta Value%
00105 se_hdr = setfield(se_hdr, 'echo6_window', fread(fid,1,'uint16')); %Echo 6 Window Value%
00106 se_hdr = setfield(se_hdr, 'echo6_level', fread(fid,1,'int16')); %Echo 6 Level Value%
00107 se_hdr = setfield(se_hdr, 'echo7_alpha', fread(fid,1,'float32')); %Echo 7 Alpha Value%
00108 se_hdr = setfield(se_hdr, 'echo7_beta', fread(fid,1,'float32')); %Echo 7 Beta Value%
00109 se_hdr = setfield(se_hdr, 'echo7_window', fread(fid,1,'uint16')); %Echo 7 Window Value%
00110 se_hdr = setfield(se_hdr, 'echo7_level', fread(fid,1,'int16')); %Echo 7 Level Value%
00111 se_hdr = setfield(se_hdr, 'echo8_alpha', fread(fid,1,'float32')); %Echo 8 Alpha Value%
00112 se_hdr = setfield(se_hdr, 'echo8_beta', fread(fid,1,'float32')); %Echo 8 Beta Value%
00113 se_hdr = setfield(se_hdr, 'echo8_window', fread(fid,1,'uint16')); %Echo 8 Window Value%
00114 se_hdr = setfield(se_hdr, 'echo8_level', fread(fid,1,'int16')); %Echo 8 Level Value%
00115 se_hdr = setfield(se_hdr, 'series_uid', fread(fid,32,'uchar')); %Series Entity Unique ID%
00116 se_hdr = setfield(se_hdr, 'landmark_uid', fread(fid,32,'uchar')); %Landmark Unique ID%
00117 se_hdr = setfield(se_hdr, 'equipmnt_uid', fread(fid,32,'uchar')); %Equipment Unique ID%
00118 se_hdr = setfield(se_hdr, 'se_padding', fread(fid,588,'uchar')); %Spare Space%
00119
00120 return
