adopt code to support SIESTA usecase 2.3 (fieldtrip#2416)

* do not use compression when saving

* also detect senstype from data.label, not only from grad/elec/opto

* allow for BIDS inheritance when searching for JSON sidecars

- the JSON sidecar might have fewer entities (e.g., no run-xxx)
- the JSON sidecar might be at a higher directory

fileio/private/bids_datafile.m

@@ -7,7 +7,7 @@
 %
 % See also BIDS_SIDECAR, BIDS_TSV, EVENTS_TSV
 
-% Copyright (C) 2020, Robert Oostenveld
+% Copyright (C) 2020-2024, Robert Oostenveld
 %
 % This file is part of FieldTrip, see http://www.fieldtriptoolbox.org
 % for the documentation and details.
@@ -52,11 +52,81 @@
 % this is the reverse procedure than the one implemented in BIDS_SIDECAR
 for i=1:numel(modality)
   f = [sprintf('%s_', entities{:}), modality{i}];
+
   % we don't know the file extension, hence we have to search with a wildcard
   d = dir(fullfile(p, [f '.*']));
   if ~isempty(d)
     datafile = fullfile(p, d.name);
-    jsonfile = fullfile(p, [sprintf('%s_', entities{:}), modality{i} '.json']);
-    break
-  end
+
+    % the json file can have fewer entities and/or can be at a higher level due to the inheritance principle
+    % we search for the one that is the "closest by"
+    foundjson = false;
+
+    if ~foundjson
+      % find the corresponding json file in the same directory as the data
+      alloptions = binoall(numel(entities));
+      for j=1:size(alloptions,1)
+        jsonfile = fullfile(p, [sprintf('%s_', entities{alloptions(j,:)}), modality{i} '.json']);
+        if exist(jsonfile, 'file')
+          foundjson = true;
+          break
+        end
+      end
+    end % if ~foundjson
+    if ~foundjson
+      % find the corresponding json file one directory up
+      alloptions = binoall(numel(entities));
+      p = fileparts(p); % remove the last directory
+      for j=1:size(alloptions,1)
+        jsonfile = fullfile(p, [sprintf('%s_', entities{alloptions(j,:)}), modality{i} '.json']);
+        if exist(jsonfile, 'file')
+          foundjson = true;
+          break
+        end
+      end
+    end % if ~foundjson
+    if ~foundjson
+      % find the corresponding json file two directories up
+      alloptions = binoall(numel(entities));
+      p = fileparts(p); % remove the last directory
+      for j=1:size(alloptions,1)
+        jsonfile = fullfile(p, [sprintf('%s_', entities{alloptions(j,:)}), modality{i} '.json']);
+        if exist(jsonfile, 'file')
+          foundjson = true;
+          break
+        end
+      end
+    end % if ~foundjson
+    if ~foundjson
+      % find the corresponding json file three directories up
+      alloptions = binoall(numel(entities));
+      p = fileparts(p); % remove the last directory
+      for j=1:size(alloptions,1)
+        jsonfile = fullfile(p, [sprintf('%s_', entities{alloptions(j,:)}), modality{i} '.json']);
+        if exist(jsonfile, 'file')
+          foundjson = true;
+          break
+        end
+      end
+    end % if ~foundjson
+
+    if ~foundjson
+      error('cannot find JSON sidecar for %s', datafile);
+    end
+  end % ~isempty(d)
+end % for each modality
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% SUBFUNCTION
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+function r =  binoall(n)
+% this returns an 2^n by n matrix with all possible binomial options
+% See also BINORND
+
+r = zeros(2^n, n, 'logical');
+for i=0:(2^n-1)
+  r(i+1,:) = bitget(uint32(i), 1:n);
 end
+r = fliplr(r);
+r = flipud(r);
+

private/savevar.m

@@ -24,7 +24,7 @@ function savevar(filename, varname, value, hashfile)
 
 % if variable < ~500 MB, store it in old (uncompressed) format, which is faster
 if (s.bytes < 500000000)
-  save(filename, varname, '-v6');
+  save(filename, varname, '-v7', '-nocompression');
 else
   save(filename, varname, '-v7.3');
 end

utilities/ft_checkdata.m

@@ -566,17 +566,12 @@
     stype = {stype};
   end
 
-  if isfield(data, 'grad') || isfield(data, 'elec') || isfield(data, 'opto')
-    if any(strcmp(ft_senstype(data), stype))
-      okflag = 1;
-    elseif any(cellfun(@ft_senstype, repmat({data}, size(stype)), stype))
-      % this is required to detect more general types, such as "meg" or "ctf" rather than "ctf275"
-      okflag = 1;
-    else
-      okflag = 0;
-    end
+  if any(strcmp(ft_senstype(data), stype))
+    okflag = 1;
+  elseif any(cellfun(@ft_senstype, repmat({data}, size(stype)), stype))
+    % this is required to detect more general types, such as "meg" or "ctf" rather than "ctf275"
+    okflag = 1;
   else
-    % the data does not contain a sensor array
     okflag = 0;
   end