loaddta.m

A MATLAB tool for analysis of Transient Otoacoustic Emission signals

function [r,stim,name,m_case,ear,m_date,rfft,noise] = loaddta(dtafile)

%LOADDTA Load .DTA- data files created by ILO-88 
%     [r,stim,name,m_case,ear,m_date,rfft,noise] = loaddta(dtafile)
%
%   in: dtafile     name of the ilo88 data file
%
%  out: r           responses 512 x 2 matrix
%       stim        stimulus
%       name        patient's name
%       m_case        
%       ear         which ear, left or right
%       m_date      date of measurement (dd.mm.-yr)
%       rfft        fft of the response
%       noise       fft of the noise

% (c) Pekka Kumpulainen 21.6.1993 (P.K. 6.8.1993)

%currdir = pwd;
if nargin<1
    [w p] = uigetfile('*.dta'); dtafile = [p w];
    %eval(['cd ' currdir]);
end

if isempty(findstr(dtafile,'.'))
    dtafile = [dtafile,'.dta'];
end

fid = fopen(dtafile,'rb');
if fid ~= -1 
    
    if nargout < 3; 
        fseek(fid,34,-1); % skip patient information
    
    else
        fread(fid,1,'uchar');
        NAME = setstr(fread(fid,20,'uchar')');
        ind = abs(NAME)==0;
        NAME = setstr(NAME+ ' '*ind);
        NAME = dos2ansi(abs(NAME));
        m_case = NAME(18:20);name = NAME(1:17);
        N = fread(fid,1,'uchar');
        EAR = setstr(fread(fid,5,'uchar')');
        if N ~= 5; 
            m_case = [m_case EAR(2:N)];
        end
        if EAR(1) == 'l' | EAR(1) == 'L'; ear = 'Left '; 
            elseif  EAR(1) == 'r' | EAR(1) == 'R'; ear = 'Right';
            else ear = EAR;
        end    
        fread(fid,1,'uchar');
        d = setstr(fread(fid,6,'uchar')'); 
        m_date = [d([5 6]) '.' d([3 4]) '.-' d([1 2])];
    end
    
    % read data
    r = zeros(512,2);
    r(:,1) = 0.238e-6*(fread(fid,512,'short'));
    r(:,2) = 0.238e-6*(fread(fid,512,'short'));
    stim = 0.00996*(fread(fid,128,'short')-128);
    
    %read fft's if necessary
    if nargout>5; 
        rfft = fread(fid,128,'short');
        noise = fread(fid,128,'short');
    end    
end

fclose(fid);

if fid == -1
    error('Can''t open .DTA file for input!');
end;