tdclient.m

toolbox of BVQX, This is the access between BV and matlab. It will help you to analysis data from BV

function [tdreply, st] = tdclient(varargin)
% tdclient  - send request to internet TDdatabase / local TD
%
%       FORMAT:     tdreply = tdclient( x,y,z  [, kind [, cubesize, local]])
%       alt.        tdreply = tdclient([x,y,z] [, kind [, cubesize, local]])
%
%       FORMAT:     tdclient('format', formstring, formsep [, extraline])
%
% Input fields:
%       x           x coordinate of request
%       y           y coordinate of request
%       z           z coordinate of request
%       kind        string, 'SPMap', {'TalLabel'}, 'Cube', 'NGM', 'RNGM'
%       cubesize    size of voxelcube to search, {7} (from 3 through 11)
%       local       override tddolocal setting (s.b.)
%
%       formstring  specifies the output format for RNGM call
%       formsep     sets the field seperator
%       extraline   enable/disable extra line feed in output
%
%   NOTE:           in the alternative calling form, a ?x3 matrix can be
%                   processed, then a cell array will be returned!
%
% functionality range (kind) depends on setting of tddolocal variable
% at the first line of code in this file (tdclient.m). when set to any
% true value (1) tdclient will ask tdlocal(2).m with the request; then you
% only have TalLabel, Cube, and NGM available. alternatively, when set to
% a false value (0) then tdclient will try to reach the TDserver set in
% variables tdsrv/tdport with over the internet; hence SPMap will be
% available while NGM won't be!
%
% when calling with kind='Cube' or kind='NGM' a cubesize of 7mm is
% assumed if not specified in the argument list; in NGM mode, only those
% points will be returned that lie within a radius of cubesize !
%
% as a special feature, NGM offers a 1x2 cubesize argument, where the
% second value must be greater than 1 to show the first n occurences of
% the search, in command mode, seperate with a '#' (hash)
%
% tdclient can be called as a MATLAB command with strings as input!
% here are some examples for this kind of request:
%
% tdclient -21 -24  12          returns what lies at -21,-24,12
% tdclient   5  17  -5 Cube 9   returns the number of voxels for each label
%                                 found in a 9x9x9 cube around 5,17,-5
% tdclient -22 -21  22 NGM 5    finds the nearest Gray Matter from
%                                 -22,-21,22; with a radius of 5mm
% tdclient -22 -21  22 NGM 7#5  find the 5 nearest distinct Gray Matter
%                                 occurences from the same point, radius 7mm
%
% when run in command mode, the return value will always be formatted
%
% due to the fact that extraline is a numeric value, the formatting call
% must be made as a function:
%
% tdclient('format','$xc;$yc;$zc;$diff;$xd;$yd;$zd;$t1;$t3;$t5;',';',0);
%
% where $xc, $yc, $zc are replaced with the actual coordinates of the
% gray matter found, $diff is the calculated distance between the
% requested and the actual coordinates, and $xd, $yd, $zd are the
% differences for each coordinate. $t1 through $t5 are the text patterns
% for each group of location descriptors.
%
% See also tdlocal, tdlocal2.

% Version:  v0.6f
% Build:    7061812
% Date:     Jun-18 2007, 12:42 PM CET
% Author:   Jochen Weber, Brain Innovation, B.V., Maastricht, NL
% URL/Info: http://wiki.brainvoyager.com/BVQXtools

% persistent variable
persistent td_iset;
if isempty(td_iset)
    td_iset.gui = [];
end

% enough arguments ?
if nargin < 1
    error( ...
        'BVQXtools:TooFewArguments',...
        'Too few arguments.' ...
    );
end

% preset output
if nargout > 0
    tdreply = '';
end

% settings
tdsrv     = 'biad73.uthscsa.edu';
tdport    = '1600';
tddecimal = '.';
tdformat  = '$xc;$yc;$zc;$diff;$xd;$yd;$zd;$t1;$t2;$t3;$t4;$t5;';
tdforsep  = ';';
tdforxln  = false;
tddolocal = 2;
tdsilent  = false;
cmdpre    = [tdsrv ' ' tdport ' '];
indent    = '                 ';
lfc       = char(10);

% test tdlocal()
if tddolocal > 0
   try
       if tddolocal == 1
           tdlocal(2, 0, 0, 0);
       else
           tdlocal2(2, 0, 0, 0);
       end
   catch
       error( ...
           'BVQXtools:InternalError',...
           'Call to function tdlocal/tdlocal2(...) failed.' ...
       );
   end
end

% - set initial values, settings...
kinds  = struct(   ...
    'spmap',    1, ...
    'tallabel', 2, ...
    'cube',     3, ...
    'ngm',      5, ...
    'rngm',    -5  ...
);

if ischar(varargin{1}) && ...
    strcmpi('visual', varargin{1})
    if nargin == 1 && ...
        isempty(td_iset.gui) 
        try
            tfgp = bvqxtools_path('tfg');
            tfgc = [];
            if exist([tfgp '/tdclient_tfg.mat'], 'file') == 2
                tfgc = [tfgp '/tdclient_tfg.mat'];
                tfga = {'IsStruct'};
            elseif exist([tfgp '/tdclient.tfg'], 'file') == 2
                tfgc = [tfgp '/tdclient.tfg'];
                tfga = {};
            else
                error( ...
                    'BVQXtools:InternalError', ...
                    'Error parsing tdclient GUI window TFG.' ...
                );
            end
            hWnd = BVQXfigure(tfgc, tfga{:});
            td_iset.gui.hWnd = hWnd;
            td_iset.gui.hnds = hWnd.TagStruct;
            tdclient('visual', 'clear');
        catch
            warning( ...
                'BVQXtools:BVQXfigureError', ...
                'Couldn''t open GUI' ...
            );
            return;
        end
    elseif nargin == 1
        return;
    elseif ischar(varargin{2})
        switch (lower(varargin{2}))
            case {'clear'}
                try
                    td_iset.gui.hnds.ED_tdclient_input.String = '';
                    td_iset.gui.hnds.ED_tdclient_result.String = '';
                    td_iset.gui.hnds.RB_tdclient_TALLabel.DoCallback;
                    td_iset.gui.hnds.ED_tdclient_SrcRange.String = '7';
                    td_iset.gui.hnds.ED_tdclient_NGMoccur.String = '5';
                catch
                    % do nothing
                end
            case {'close'}
                try
                    td_iset.gui = [];
                    hWnd = BVQXfigure('Wnd_BVQXtools_tdclient');
                    hWnd.Delete;
                catch
                    % do nothing
                end
            case {'send'}
                try
                    try
                        input = ...
                            evalin('base', ...
                            ['[' td_iset.gui.hnds.ED_tdclient_input.String ']']);
                    catch
                        try
                            input = ...
                                eval(['[' ...
                                td_iset.gui.hnds.ED_tdclient_input.String ']'], '[]');
                        catch
                            input = '';
                        end
                    end
                    if isempty(input)
                        return;
                    end
                    if td_iset.gui.hnds.RB_tdclient_TALLabel.IsActive
                        retval = tdclient(input);
                    elseif td_iset.gui.hnds.RB_tdclient_CubeSearch.IsActive
                        try
                            srcrange = round(str2double( ...
                                td_iset.gui.hnds.ED_tdclient_SrcRange.String));
                        catch
                            srcrange = 7;
                        end
                        retval = tdclient(input, 'Cube', srcrange);
                    else
                        try
                            srcrange = round(str2double( ...
                                td_iset.gui.hnds.ED_tdclient_SrcRange.String));
                            numhits = round(str2double( ...
                                td_iset.gui.hnds.ED_tdclient_NGMoccur.String));
                        catch
                            srcrange = 7;
                            numhits = 5;
                        end
                        retval = tdclient(input, 'NGM', [srcrange, numhits]);
                    end
                    if iscell(retval)
                        retval = gluetostring(retval);
                    end
                    td_iset.gui.hnds.ED_tdclient_result.String = retval;
                catch
                    % do nothing
                end
        end
    end
	return;
end

% argument check
if ischar(varargin{1}) && ....
    any(varargin{1}(:) == ',')
    if tddolocal == 0
        [st, tdreply] = i_tdclient([cmdpre varargin{1}]);
    else
        st = 0;
        tdargs = splittocell(varargin{1}, ',', 1);
        if any(tdargs{1} == ':')
            tdargsx = splittocell(tdargs{1}, ':', 1);
            tdargs{1} = tdargsx{1};
            tdargs{end+1} = tdargsx{2};
        end
        if tddolocal == 1
            tdreply = tdlocal(tdargs{1:end});
        else
            tdreply = tdlocal2(tdargs{1:end});
        end
    end
    tdreply = [varargin{1} ': ' tdreply];

else
    if ischar(varargin{1})
        if nargin < 3
            error( ...
                'BVQXtools:TooFewArguments',...
                'Too few arguments. Try ''help %s''.',...
                mfilename ...
            );
        end
        varargin{1} = str2double(varargin{1});
        if isempty(varargin{1})
            error( ...
                'BVQXtools:BadArgument',...
                'Bad argument.' ...
            );
        end