Discover MakerZone

MATLAB and Simulink resources for Arduino, LEGO, and Raspberry Pi

Learn more

Discover what MATLAB® can do for your career.

Opportunities for recent engineering grads.

Apply Today

Thread Subject:
Import Tucker Davis Technology (TDT) Data

Subject: Import Tucker Davis Technology (TDT) Data

From: Weimin

Date: 28 Apr, 2009 00:42:02

Message: 1 of 4

Does anyone know any or can nayone write MatLab scripts to import data from TDT Tank to sigTOOL or other data analysis programs?

Subject: Import Tucker Davis Technology (TDT) Data

From: Malcolm Lidierth

Date: 12 May, 2009 12:28:01

Message: 2 of 4

This should be possible using the TDT tank and neuroshare compliant dll provided by Tucker Davis.

Subject: Import Tucker Davis Technology (TDT) Data

From: Weimin

Date: 12 May, 2009 18:07:02

Message: 3 of 4

Thanks, Malcolm. I will give a try by using neuroshare.
"Malcolm Lidierth" <> wrote in message <gubq0h$1h6$>...
> This should be possible using the TDT tank and neuroshare compliant dll provided by Tucker Davis.

Subject: Import Tucker Davis Technology (TDT) Data

From: kai yu

Date: 2 Apr, 2011 04:54:04

Message: 4 of 4

"Weimin" <> wrote in message <gt5jcq$ri3$>...
> Does anyone know any or can nayone write MatLab scripts to import data from TDT Tank to sigTOOL or other data analysis programs?

i found this code, hope it is useful to you!
function sampFreq = ecogTDTData2Matlab(rawDataPath,blockName,chans2Convert,dataTag,outputPath,outputPrefix)
% sampFreq = ecogTDTData2Matlab(rawDataPath,blockName,chans2Convert,dataTag,outputPath,outputPrefix) Covert a set of channels recorded with a TDT sytem to matlab
% PURPOSE: Transform raw TDT data to matlab files
% Each channel is written into a separate mat-file
% Precision is single.
% rawDataPath: Where the raw data live
% e.g. D:\Data\Ecog\bp21FingerMov090127\GP21_B12
% blockName: The data block of interest e.g. GP21_B12. Typically
% this is the last part of the pathname and is included in
% the name of the raw files. An error is thron and the
% availabe blocks are listed if the requested block does
% not exist in the tank.
% chans2Convert: A list of channel to convert.
% 1:64 will convert the first 64 channels. Currently ther
% is no checking if the channels exist. You will just
% find very small files.
% dataTag: A case sensitive four letter tag defining what we read
% Most likely 'ECOG' or 'Wave' if you want to convert ecog
% channels.
% Most likely 'ANIN' if you want to convert analogue
% channels.
% Check 'ecogTDTGetDataTags()' if you are not sure
% outputPath: The path where the exported dat are written.
% Typically the same as rawDataPath
% outputPrefix: The prefix for the output files. Helps to discrimintae
% the contents.
% For ECoG use 'gdat_'
% For analogue channels use 'analog_'
% OUTPUT: The data are written to mat-files. The file contains
% data: A vector containing the time series
% params: A structure holding descriptive information
% tagName: The tag name of the channel
% chanNum: The number of the current channel
% sampFreq: The sampling frequency
% nChannels: The total number of channels with
% that tag
% sampFreq: The sampling frequency used to collect the data
% REQUIREMENTS: This function requires Microsoft Windows, the
% TDT drivers, and OpenDeveloper (the latter requires openEx)
% installed
% Get the software from OpenEX and OpenDeveloper
% are password protected.
% ecogTDTData2Matlab('D:\Data\Berkeley\Ecog\bp22090323\resultsFile\GP22_B28','GP22_B28',1:80,'HECG','D:\Data\Berkeley\Ecog\bp22090323\resultsFile\GP22_B28','gdat_')
% reads 80 channels of ECoG-data from tag 'HECG' and writes each of them to
% a file starting with name 'gdat_[1:80]'. Returns the sampling frequency
% Use

% 090130 got the original code from LS
% 090130 JR included comments, OS-check, try catch, removed repeated calls and modified path handling
% 090324 JR rewritten to read only data associated with one tag
% tag is passed as an argument
% fileprefix is passed as an argument
% channels are passed as an argument
% 090324 JR The time series is now named 'data' and a params structure is
% written together with the timeseries.
% 090321 JR Tank name is now found in the directory. No more assumptions
% about specific prefixes are made
% If a block name is not available, all available names are listed before throwing an error
% 090812 JR Use dataTag as outputprefix if non was passed instead of chrashing

switch computer
    case 'PCWIN'
    case 'PCWIN64'
        error('Requires Microsoft Windows. Uses ActiveX controls.')

TT = actxcontrol('TTank.X');
if (e==0)
    error('Cannot connect to server');

%e=invoke(TT,'OpenTank',[rawDataPath '\Ecog2_' blockName],'R');%PLACE PATH OF TANK TO COLLECT DATA FROM
% % the sumption here is that there is only one tank in a directory and
    % that the extension is tev
 e=invoke(TT,'OpenTank',[rawDataPath filesep tankName],'R');%PLACE PATH OF TANK TO COLLECT DATA FROM
 if (e==0)
    %error(['Cannot open data tank: ' rawDataPath '\Ecog2_' blockName]);
    error(['Cannot open data tank: ' rawDataPath filesep blockName ' No *.tev file there']);
e=invoke(TT,'SelectBlock',blockName);%PLACE NAME OF BLOCK
if e==0
    display('Available block names are:')
    while (~isempty(TT.QueryBlockName(count)))
    error(['Cannot find block name: ' blockName]);
totaltime = floor(invoke(TT,'CurBlockStopTime')- invoke(TT,'CurBlockStartTime'));%Gets total tank duration.
display(['Found ' num2str(totaltime) ' second of data']);



for k=1:length(chans2Convert)
    e=invoke (TT,'SetGlobals', ['Options = ALL; WavesMemLimit = 134217728;Channel=' num2str(chans2Convert(k)) '; T1=0; T2=' num2str(totaltime)]);
    data=invoke(TT, 'ReadWavesV', dataTag)';%read data
    if ~isempty(data)
        if ~exist('outputPrefix','var') || isempty(outputPrefix)
    disp(['saving ' outputPrefix num2str(chans2Convert(k)) ' of ' num2str(chans2Convert(end))]);
    %eval([outputPrefix num2str(chans2Convert(k)) ' = data;']);
    %eval(['save ' outputPrefix num2str(chans2Convert(k)) ' ' outputPrefix num2str(chans2Convert(k)) 'params;']);
    save([rawDataPath filesep outputPrefix num2str(chans2Convert(k))],'data','params');
        disp(['No data found for channeL: ' num2str(chans2Convert(k)) '. Skipping!']);

%clear gdat*
if 0
% anin channels
for ch = nchans_anin
    e=invoke (TT,'SetGlobals', ['Options = ALL; WavesMemLimit = 134217728; Channel=' num2str(ch) '; T1=0; T2=' num2str(totaltime)]);
    gdat=invoke(TT, 'ReadWavesV', 'ANIN')';%grabs the file
    disp(['saving analog ' num2str(ch) ' of ' num2str(nchans_anin(end))]);
    eval(['analog_' num2str(ch) ' = gdat;']);
    eval(['save analog_' num2str(ch) ' analog_' num2str(ch) ';']);

invoke(TT,'CloseTank')%closes the server.
% got to where we started

catch %Make sure the server is closed if anything chrashes
    invoke(TT,'CloseTank')%closes the server.
    display(lasterr) %Tell the user the problem
% cd(oldPath)

<html xmlns=""><body xmlns=""><div xmlns=""><img width='77' height='5' src='chrome://livemargins/skin/monitor-background-horizontal.png' style='position: absolute;left:-77px;top:-5px'/> <img src='chrome://livemargins/skin/monitor-background-vertical.png' style='position: absolute;left:0;top:-5px;'/> <img id='monitor-play-button' src='chrome://livemargins/skin/monitor-play-button.png' onmouseover='' onmouseout='' style='position: absolute;left:1px;top:0;opacity:0.5;cursor:pointer'/>

Tags for this Thread

What are tags?

A tag is like a keyword or category label associated with each thread. Tags make it easier for you to find threads of interest.

Anyone can tag a thread. Tags are public and visible to everyone.

Contact us