function nbrhood=cosmo_interval_neighborhood(ds, label, varargin)
% compute neighborhoods stretching intervals
%
% nbrhood=cosmo_interval_neighborhood(ds, label, 'radius', radius)
%
% Inputs:
% ds dataset struct
% label dimension label in ds.a.fdim.labels
% 'radius', r neighborhood radius
%
% Returns:
% nbrhood struct with fields:
% .a } dataset and feature attributes of neighborhood
% .fa }
% .neighbors If ds has N values in the dimension label, then
% neighbors is a Nx1 cell with indices of features
% where the indices differ at most radius from each other.
%
% Examples
% % Illustrate the 'time' dimension in MEEG time-frequency dataset,
% ds=cosmo_synthetic_dataset('type','timefreq','size','big');
% %
% % neighborhoods with bins 5 (=2*1+1) frequency bands wide
% % (every neighborhood contains all channels and time points)
% nbrhood=cosmo_interval_neighborhood(ds,'freq','radius',2);
% cosmo_disp(nbrhood.a.fdim)
% %|| .labels
% %|| { 'freq' }
% %|| .values
% %|| { [ 2 4 6 ... 10 12 14 ]@1x7 }
% cosmo_disp(nbrhood.fa.freq)
% %|| [ 1 2 3 ... 5 6 7 ]@1x7
% cosmo_disp(nbrhood.neighbors)
% %|| { [ 1 2 3 ... 9.48e+03 9.48e+03 9.49e+03 ]@1x4590
% %|| [ 1 2 3 ... 9.79e+03 9.79e+03 9.79e+03 ]@1x6120
% %|| [ 1 2 3 ... 1.01e+04 1.01e+04 1.01e+04 ]@1x7650
% %|| :
% %|| [ 613 614 615 ... 1.07e+04 1.07e+04 1.07e+04 ]@1x7650
% %|| [ 919 920 921 ... 1.07e+04 1.07e+04 1.07e+04 ]@1x6120
% %|| [ 1.22e+03 1.23e+03 1.23e+03 ... 1.07e+04 1.07e+04 1.07e+04 ]@1x4590 }@7x1
%
% % ds is an MEEG dataset with a time dimension
% ds=cosmo_synthetic_dataset('type','timelock','size','big');
% %
% % Neighborhoods just the frequency bin itself
% % (every neighborhood contains all channels)
% nbrhood=cosmo_interval_neighborhood(ds,'time','radius',2);
% cosmo_disp(nbrhood.a.fdim)
% %|| .labels
% %|| { 'time' }
% %|| .values
% %|| { [ -0.2 -0.15 -0.1 ... 0 0.05 0.1 ]@1x7 }
% cosmo_disp(nbrhood.fa.time)
% %|| [ 1 2 3 ... 5 6 7 ]@1x7
% cosmo_disp(nbrhood.neighbors)
% %|| { [ 1 2 3 ... 916 917 918 ]@1x918
% %|| [ 1 2 3 ... 1.22e+03 1.22e+03 1.22e+03 ]@1x1224
% %|| [ 1 2 3 ... 1.53e+03 1.53e+03 1.53e+03 ]@1x1530
% %|| :
% %|| [ 613 614 615 ... 2.14e+03 2.14e+03 2.14e+03 ]@1x1530
% %|| [ 919 920 921 ... 2.14e+03 2.14e+03 2.14e+03 ]@1x1224
% %|| [ 1.22e+03 1.23e+03 1.23e+03 ... 2.14e+03 2.14e+03 2.14e+03 ]@1x918 }@7x1
%
%
% Notes:
% - to combine neighborhoods from different dimensions (such as
% time, freq, chan, use cosmo_neighborhood
% - the output can be used for a searchlight using cosmo_searchlight
%
% See also: cosmo_neighborhood, cosmo_searchlight
%
% # For CoSMoMVPA's copyright information and license terms, #
% # see the COPYING file distributed with CoSMoMVPA. #
