ge_readrawheaderrdbrec.m

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

rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ps_command', fread(fid,1,'int32'));	 % PS internal use only	%
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ps_mps_r1', fread(fid,1,'int32'));	 % PS MPS R1 setting		1 - 7 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ps_mps_r2', fread(fid,1,'int32'));	 % PS MPS R2 setting		1 - 30 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ps_mps_tg', fread(fid,1,'int32'));	 % PS MPS Transmit gain setting	0 - 200%
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ps_mps_freq', fread(fid,1,'int32')); % PS MPS Center frequency hz	+-3.40282e38 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ps_aps_r1', fread(fid,1,'int32'));	 % PS APS R1 setting		1 - 7 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ps_aps_r2', fread(fid,1,'int32'));	 % PS APS R2 setting		1 - 30 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ps_aps_tg', fread(fid,1,'int32'));	 % PS APS Transmit gain setting	0 - 200%
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ps_aps_freq', fread(fid,1,'int32')); % PS APS Center frequency hz	+-3.40282e38 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ps_scalei', fread(fid,1,'float32')); % PS rational scaling 		+-3.40282e38 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ps_scaleq', fread(fid,1,'float32')); % PS unused %			
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ps_snr_warning', fread(fid,1,'int32')); % PS noise test 0=16 1=32 bits	0, 1 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ps_aps_or_mps', fread(fid,1,'int32')); % PS prescan order logic	0 - 5 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ps_mps_bitmap', fread(fid,1,'int32')); % PS bit mask			4 bits%
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ps_powerspec', fread(fid,256,'char')); % PS %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ps_filler1', fread(fid,1,'int32'));	 % PS filler %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ps_filler2', fread(fid,1,'int32'));	 % PS filler %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_rec_noise_mean', fread(fid,16,'float32')); % PS mean noise each receiver +-3.40282e38 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_rec_noise_std', fread(fid,16,'float32')); % PS noise calc for muti rec	+-3.40282e38 %
rdb_hdr = setfield(rdb_hdr, 'halfecho', fread(fid,1,'int16')); % spectro full, half echo 0, 1 %
% 858 bytes %
% New fields 02-19-92 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_im_size_y', fread(fid,1,'int16')); % rh???? 			0 - 512 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_data_collect_type1', fread(fid,1,'int32')); % rh???? bit mask		31 bits %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_freq_scale', fread(fid,1,'float32')); % rh???? freq k-space step +-3.40282e38 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_phase_scale', fread(fid,1,'float32')); % rh???? freq k-space step +-3.40282e38 %
% 14 bytes %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ovl', fread(fid,1,'int16')); % rhovl - overlaps for MOTSA % 
% Phase Correction Control Param. %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pclin', fread(fid,1,'int16')); % Linear Corr. 0:off, 1:linear, 2:polynomial %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pclinnpts', fread(fid,1,'int16')); % fit number of points %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pclinorder', fread(fid,1,'int16')); % fit order %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pclinavg', fread(fid,1,'int16')); % linear phase corr avg 0:off, 1:on %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pccon', fread(fid,1,'int16')); % Const Corr. 0:off, 1:Ky spec., 2:polyfit(2/ilv), 3:polyfit(1/ilv) %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcconnpts', fread(fid,1,'int16')); % fit number of points %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcconorder', fread(fid,1,'int16')); % fit order %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcextcorr', fread(fid,1,'int16')); % external correction file 0:don't use, 1: use %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcgraph', fread(fid,1,'int16')); % Phase Correction coef. image 0:off, 1:linear & constant %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcileave', fread(fid,1,'int16')); % Interleaves to use for correction: 0=all, 1=only first %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_hdbestky', fread(fid,1,'int16')); % bestky view for fractional Ky scan %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcctrl', fread(fid,1,'int16')); % phase correction research control %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcthrespts', fread(fid,1,'int16')); % 2..512 adjacent points %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcdiscbeg', fread(fid,1,'int16')); % 0..512 beginning point to discard %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcdiscmid', fread(fid,1,'int16')); % 0..512 middle point to discard %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcdiscend', fread(fid,1,'int16')); % 0..512 ending point to discard %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcthrespct', fread(fid,1,'int16')); % Threshold percentage %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcspacial', fread(fid,1,'int16')); % Spacial best ref scan index 0..512 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pctemporal', fread(fid,1,'int16')); % Temporal best ref scan index 0..512 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcspare', fread(fid,1,'int16')); % spare for phase correction %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ileaves', fread(fid,1,'int16')); % Nunmber of interleaves %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_kydir', fread(fid,1,'int16')); % Ky traversal dircetion 0: top-down, 1:center out %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_alt', fread(fid,1,'int16')); % Alt read sign 0=no, 1=odd/even, 2=pairs %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_reps', fread(fid,1,'int16')); % Number of scan repetitions %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ref', fread(fid,1,'int16')); % Ref Scan 0: off 1: on %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcconnorm', fread(fid,1,'float32')); % Constant S term normalization factor %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcconfitwt', fread(fid,1,'float32')); % Constant polyfit weighting factor %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pclinnorm', fread(fid,1,'float32')); % Linear S term normalization factor %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pclinfitwt', fread(fid,1,'float32')); % Linear polyfit weighting factor %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcbestky', fread(fid,1,'float32')); % Best Ky location %
% VRG Filter param %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_vrgf', fread(fid,1,'int32')); % control word for VRG filter %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_vrgfxres', fread(fid,1,'int32')); % control word for VRGF final x resolution %
% Bandpass AsymmetryCorrection Param. %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_bp_corr', fread(fid,1,'int32')); % control word for bandpass asymmetry %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_recv_freq_s', fread(fid,1,'float32')); % starting frequency (+62.5) %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_recv_freq_e', fread(fid,1,'float32')); % ending frequency (-62.5) %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_hniter', fread(fid,1,'int32')); % Selects the number of
% iterations used in homodyne processing %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_fast_rec', fread(fid,1,'int32')); % Added for homodyne II, tells if
% the fast receiver is being used
% and the lpf setting of the fast
				% receiver, 0: fast receiver off,
				% 1 - 5: lpf settings %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_refframes', fread(fid,1,'int32')); % total # of frames for ref scan %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_refframep', fread(fid,1,'int32')); % # of frames per pass for a ref scan %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_scnframe', fread(fid,1,'int32')); % total # of frames for a entire scan %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pasframe', fread(fid,1,'int32')); % # of frames per pass %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user_usage_tag', fread(fid,1,'uint32')); % for spectro %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user_fill_mapMSW', fread(fid,1,'uint32')); % for spectro %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user_fill_mapLSW', fread(fid,1,'uint32')); % for Spectro %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user20', fread(fid,1,'float32')); % all following usercv are for spectro %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user21', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user22', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user23', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user24', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user25', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user26', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user27', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user28', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user29', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user30', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user31', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user32', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user33', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user34', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user35', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user36', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user37', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user38', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user39', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user40', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user41', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user42', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user43', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user44', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user45', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user46', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user47', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_user48', fread(fid,1,'float32'));
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcfitorig', fread(fid,1,'int16')); % Adjust view indexes if set so bestky view = 0 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcshotfirst', fread(fid,1,'int16')); % First view within an echo group used for fit%
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcshotlast', fread(fid,1,'int16')); % Last view within an echo group used for fit %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcmultegrp', fread(fid,1,'int16')); % If = 1, force pts from other egrps to be used %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pclinfix', fread(fid,1,'int16')); % If = 2, force slope to be set to pclinslope %
 % If = 1, neg readout slope = pos readout slope %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcconfix', fread(fid,1,'int16')); % If = 2, force slope to be set to pcconslope %
 % If = 1, neg readout slope = pos readout slope %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pclinslope', fread(fid,1,'float32')); % Value to set lin slope to if forced %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pcconslope', fread(fid,1,'float32')); % Value to set con slope to if forced %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_pccoil', fread(fid,1,'int16')); % If 1,2,3,4, use that coil's results for all %
% Variable View Sharing %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_vvsmode', fread(fid,1,'int16')); % Variable view sharing mode %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_vvsaimgs', fread(fid,1,'int16')); % number of original images %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_vvstr', fread(fid,1,'int16')); % TR in microseconds %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_vvsgender', fread(fid,1,'int16')); % gender: male or female %
% 3D Slice ZIP %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_zip_factor', fread(fid,1,'int16')); % Slice ZIP factor: 0=OFF, 2, or 4 %
% Maxwell Term Correction Coefficients %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_maxcoef1a', fread(fid,1,'float32')); % Coefficient A for flow image 1 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_maxcoef1b', fread(fid,1,'float32')); % Coefficient B for flow image 1 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_maxcoef1c', fread(fid,1,'float32')); % Coefficient C for flow image 1 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_maxcoef1d', fread(fid,1,'float32')); % Coefficient D for flow image 1 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_maxcoef2a', fread(fid,1,'float32')); % Coefficient A for flow image 2 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_maxcoef2b', fread(fid,1,'float32')); % Coefficient B for flow image 2 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_maxcoef2c', fread(fid,1,'float32')); % Coefficient C for flow image 2 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_maxcoef2d', fread(fid,1,'float32')); % Coefficient D for flow image 2 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_maxcoef3a', fread(fid,1,'float32')); % Coefficient A for flow image 3 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_maxcoef3b', fread(fid,1,'float32')); % Coefficient B for flow image 3 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_maxcoef3c', fread(fid,1,'float32')); % Coefficient C for flow image 3 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_maxcoef3d', fread(fid,1,'float32')); % Coefficient D for flow image 3 %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_ut_ctrl', fread(fid,1,'int32')); % System utility control variable %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_dp_type', fread(fid,1,'int16')); % EPI II diffusion control cv %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_arw', fread(fid,1,'int16')); % Arrhythmia rejection window(percentage:1-100)%
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_vps', fread(fid,1,'int16'));	% View Per Segment for FastCine % 
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_mcReconEnable', fread(fid,1,'int16')); % N-Coil recon map %
rdb_hdr = setfield(rdb_hdr, 'rdb_hdr_excess', fread(fid,420,'int16')); % free space for later expansion %

return;