1 function [HDR] = ssave(FILENAME,DATA,TYPE,Fs,gdftyp)
8 % HDR = ssave(HDR,data);
9 % HDR = ssave(FILENAME,data,TYPE,Fs);
41         HDR = FILENAME;  variable
42         if isfield(HDR,'FileName'),
43                 FILENAME = HDR.FileName;
49         HDR.FileName = FILENAME;
50         HDR.SampleRate = Fs; 
52         HDR.TYPE = 'native';
56 	[HDR.SPR, HDR.NS] = size(DATA); HDR.NRec = 1; 
57 %	HDR.AS = rmfield(HDR.AS,'SPR'); 
58 …HDR.TYPE,'BDF') || strcmp(HDR.TYPE,'EDF') || strcmp(HDR.TYPE,'GDF')) && (~isfield(HDR,'DigMax') ||…
59 		HDR.PhysMax = max(DATA,[],1);
60 		HDR.PhysMin = min(DATA,[],1);
61 		ix = find(HDR.PhysMax == HDR.PhysMin);
62 		HDR.PhysMin(ix) = HDR.PhysMin(ix) - 1;
63 		if strcmp(HDR.TYPE,'BDF')
64 		   	[datatyp,HDR.THRESHOLD,datatypes,HDR.bits,HDR.GDFTYP] = gdfdatatype(511+24*ones(HDR.NS,1));
65 		   	HDR.DigMax = HDR.THRESHOLD(:,2)';
66 		   	HDR.DigMin = HDR.THRESHOLD(:,1)';
67 		elseif strcmp(HDR.TYPE,'EDF')
68 		   	[datatyp,HDR.THRESHOLD,datatypes,HDR.bits,HDR.GDFTYP] = gdfdatatype(3*ones(HDR.NS,1));
69 		   	HDR.DigMax = HDR.THRESHOLD(:,2)';
70 		   	HDR.DigMin = HDR.THRESHOLD(:,1)';
71 		elseif strcmp(HDR.TYPE,'GDF')
72 		   	[datatyp,HDR.THRESHOLD,datatypes,HDR.bits,HDR.GDFTYP] = gdfdatatype(16*ones(HDR.NS,1));
73 			HDR.DigMax = HDR.PhysMax;
74 			HDR.DigMin = HDR.PhysMin;
76 			HDR.DigMax = HDR.PhysMax;
77 			HDR.DigMin = HDR.PhysMin;
83         		HDR.TYPE = 'GDF';
84 	        	HDR.VERSION = 2;
85 			HDR.NS = 0; 
87         		HDR.TYPE = 'GDF';
88 	        	HDR.VERSION = 2;
90         		HDR.TYPE = 'GDF';
91       	 	 	HDR.VERSION = 1.25;
93 	        	HDR.TYPE = TYPE; 	% type of data format
96 	HDR = sopen(HDR,'w');  variable
98 		HDR = swrite(HDR,DATA);  variable
100 	HDR = sclose(HDR);  variable
106 if all([length(HDR.EVENT.POS), length(HDR.EVENT.TYP)]),
110 	HDR.EVENT.CodeDesc  = H.CodeDesc;
111 	HDR.EVENT.CodeIndex = H.CodeIndex;
112 	if isfield(HDR.EVENT,'DUR')
113 	        HDR.EVENT.POS = [HDR.EVENT.POS; HDR.EVENT.POS + HDR.EVENT.DUR];
114 	        HDR.EVENT.TYP = [HDR.EVENT.TYP; HDR.EVENT.TYP + hex2dec('8000')];
118 	[HDR.EVENT.POS, ix] = sort(HDR.EVENT.POS);
119 	HDR.EVENT.TYP       = HDR.EVENT.TYP(ix);
120 	[TYP, IX, IY]       = unique(HDR.EVENT.TYP);
123 	fid   = fopen(fullfile(HDR.FILE.Path,[HDR.FILE.Name,'.C07']),'w');
125 …fprintf(fid,'%2i %s (%s)\r\n', k, HDR.EVENT.CodeDesc{mod(TYP(k),2^15)==HDR.EVENT.CodeIndex}, OnOff…
130 	fid = fopen(fullfile(HDR.FILE.Path,[HDR.FILE.Name,'.007']),'w');
131 	fprintf(fid,'%8i %i\r\n', [round(HDR.EVENT.POS(:)),IY(:)]');