deal with RawSources for surface derivatives

Author: Remi-Gau

src/utils/identify_rawsources.m

@@ -3,6 +3,8 @@
     % find the most likely files in the raw dataset
     % that was used to create this derivatives
     %
+    % Assumption: the raw file was zipped
+    %
     % USAGE::
     %
     %   rawsource = identify_rawsources(derivatives)
@@ -43,6 +45,7 @@
 
     % - remove prefix
     % - remove eventual derivatives entities
+    % - change surface extension to a volume one
     % - use only .nii.gz
     bf = bids.File(derivatives, 'verbose', verbose, 'use_schema', false);
 
@@ -53,6 +56,11 @@
     for i = 1:numel(idx)
       bf.entities.(entities{idx(i)}) = '';
     end
+    
+    if strcmp(bf.extension, '.surf.gii')
+      bf.extension = '.nii';
+    end
+    
     if strcmp(bf.extension, '.nii')
       bf.extension = '.nii.gz';
     end

tests/test_identify_rawsources.m

@@ -8,6 +8,24 @@
     initTestSuite;
 end
 
+function test_identify_rawsources_surface()
+
+    input_output = {'c1sub-01_T1w.surf.gii', 'sub-01_T1w.nii.gz'};
+
+    verbose = false;
+    
+    map = default_mapping();
+
+    for i = 1:size(input_output, 1)
+
+        rawsource = identify_rawsources(input_output{i, 1}, map, verbose);
+
+        assertEqual(rawsource, {['sub-01/' input_output{i, 2}]});
+
+    end
+
+end
+
 function test_identify_rawsources_when_der_entities()
 
     input_output = {'sub-01_desc-skullstripped_T1w.nii', 'sub-01_T1w.nii.gz'};