|
|
【答案】应助回帖
感谢参与，应助指数 +1
%% Write an ECG header in MIT format
% This function receives a "header" struct with ECG signal propeties and
% writes [header.name '.hea'] file into "header_path".
%
% Example
%
% MIT_filename = 'your_filename';
% MIT_path = ['.' filesep];
% ECG_header.recname = MIT_filename;
% ECG_header.nsig = 8;
% ECG_header.nsamp = 9978;
% ECG_header.freq = 1000; %Hz
% ECG_header.desc = char( {'I','II','V1','V2','V3','V4','V5','V6'} );
% ECG_header.adczero = zeros( ECG_header.nsig,1 );
% ECG_header.gain = ones( ECG_header.nsig,1 );
% ECG_header.units = repmat( 'uV',ECG_header.nsig,1 );
% ECG_header.btime = '00:00:00';
% ECG_header.bdate = '01/01/2000';
%
% writeheader( MIT_path, ECG_header );
%
% See also read_ECG, ECGwrapper
%
% Author: Salvador Olmos, Juan Pablo Mart穩nez and/or Juan Bolea
% adapted to ECGkit by: Mariano Llamedo Soria
% <matlab:web( 'mailto:llamedom@electron.frba.utn.edu.ar','-browser' ) ( email )>
% Version: 0.1 beta
% Birthdate: 5/01/2014
% Last update: 19/11/2014
%
function writeheader( header_path, header )
%
% WRITHEAD function writes header file for signal data struct in directory header_path
% Input parameters.
% header_path: directory of work
% header: struct of header to write
if( header_path( end ) ~= filesep )
header_path = [header_path filesep];
end
% Opening ASCII header file
if ~isempty( header.recname )
fid=fopen( [header_path header.recname '.hea'], 'w' );
header.fname = repmat( [ header.recname '.dat'], header.nsig,1 );
else
fid = fopen( [header_path header.fname( 1,: ) '.hea'], 'w' );
end
if ( fid<0 )
error( ['Can not write ' [header_path header.recname '.hea'] '\n'] );
end
try
% writing first line of record_name, # signals and so on
fprintf( fid,'%s %d %d %d\n',header.recname,header.nsig );
if isfield( header,'freq' )
fprintf( fid,'%d ',header.freq );
end
if isfield( header,'nsamp' )
fprintf( fid,'%d ',header.nsamp );
end
if isfield( header,'btime' )
fprintf( fid,'%s ',header.btime );
end
if isfield( header,'bdate' );
fprintf( fid,'%s ',header.bdate );
end
fmt = 16;
% writing the rest of lines
for ii=1:header.nsig
if isfield( header,'spf' ) || isfield( header,'skew' ) || isfield( header,'offset' )
fprintf( fid,'\n%s %d',header.fname( ii,: ),fmt );
if isfield( header,'spf' ) && ~isnan( header.spf( ii ) )
fprintf( fid,'x%d',header.spf( ii ) );
end
if isfield( header,'skew' ) && ~isnan( header.skew( ii ) )
fprintf( fid,':%d',header.skew( ii ) );
end
if isfield( header,'offset' ) && ~isnan( header.offset( ii ) )
fprintf( fid,'+%d',header.offset( ii ) );
end
fprintf( fid,' ' );
else
fprintf( fid,'\n%s %d ',header.fname( ii,: ),fmt );
end
if isfield( header,'gain' )
fprintf( fid,'%f',header.gain( ii ) );
else
fprintf( fid,'%f',200 );
end
if( isfield( header,'baseline' ) && ~isnan( header.baseline( ii ) ) )
fprintf( fid,'( %d )',header.baseline( ii ) );
end
if ( isfield( header,'units' ) )
fprintf( fid,'/%s ',header.units( ii,: ) );
else
fprintf( fid,' ' );
end
if isfield( header,'adcres' )
fprintf( fid,'%d ',header.adcres( ii ) );
else
fprintf( fid,'%d ', 12 );
end
if isfield( header,'adczero' )
fprintf( fid,'%d ',header.adczero( ii ) );
else
fprintf( fid,'%d ', 0 );
end
if isfield( header,'initval' )
fprintf( fid,'%d ',header.adczero( ii ) );
else
fprintf( fid,'%d ', 0 );
end
if isfield( header,'cksum' )
fprintf( fid,'%d ',header.cksum( ii ) );
else
fprintf( fid,'%d ', 0 );
end
if isfield( header,'bsize' )
fprintf( fid,'%d ',header.bsize( ii ) );
else
fprintf( fid,'%d ', 0 );
end
if isfield( header,'desc' )
fprintf( fid,'%s',header.desc( ii,: ) );
end
end
fprintf( fid,'\n' );
fclose( fid );
catch MEE
fclose( fid );
rethrow( MEE );
end
參考:
https://groups.google.com/forum/ ... U0xsMY/CMXs2ZOvNoEJ |
|
回复此楼
