Modify the version of filter for tPAC (brainstorm-tools#479)

Co-authored-by: Niloofar Gharesi <niloofar.gharesi@aerial.ai>
mpompolas · Jan 12, 2022 · f278b97 · f278b97
1 parent 0b0682d
commit f278b97
Showing 1 changed file with 4 additions and 3 deletions.
diff --git a/toolbox/process/functions/process_pac_dynamic.m b/toolbox/process/functions/process_pac_dynamic.m
@@ -600,6 +600,7 @@
 bandNestingLen= max(2,1/(winLen+margin));  % Length of band nesting -- considering the resolution in FFT domain with available window length
 isMirror = 0;                              % Mirroring the data in filtering
 isRelax  = 1;                              % Attenuation of the filter in the stopband (1 => 40 dB, 0 => 60 dB)
+Method = 'bst-hfilter-2016'                % Version of the filter 
 minExtracFreq = max(1/winLen, fpBand(1));  % minimum frequency that could be extracted as nestingFreq
 doInterpolation = Options.doInterpolation; % Applying interpolation in frequency and time domain
 logCenters = Options.logCenters;           % Choose the center frequencies for f_A with log space in faBand
@@ -668,7 +669,7 @@
     bandNested = [nestedCenters(ifreq)-Fstep(ifreq),nestedCenters(ifreq)+Fstep(ifreq+1)];
 
     % Filtering in fA band
-    Xnested = bst_bandpass_hfilter(Xinput, sRate,bandNested(1), bandNested(2), isMirror, isRelax, [], fArolloff);    % Filtering
+    Xnested = bst_bandpass_hfilter(Xinput, sRate,bandNested(1), bandNested(2), isMirror, isRelax, [], fArolloff, Method);    % Filtering
     Xnested = Xnested(:,nMargin-nHilMar+1:end-nMargin+nHilMar);            % Removing part of the margin
 
     % Hilbert transform
@@ -767,11 +768,11 @@
 
         % Filtering in fP band
         if length(unique(bandNesting(:,1)))==1 && length(unique(bandNesting(:,2)))==1
-            Xnesting = bst_bandpass_hfilter(X, sRate,bandNesting(1,1), bandNesting(1,2), isMirror, isRelax, [], fProlloff);    % Filtering
+            Xnesting = bst_bandpass_hfilter(X, sRate,bandNesting(1,1), bandNesting(1,2), isMirror, isRelax, [], fProlloff, Method);    % Filtering
         else
             Xnesting = zeros(size(X));
             for i=1:length(isources)
-                Xnesting(i,:) = bst_bandpass_hfilter(X(i,:), sRate, bandNesting(i,1), bandNesting(i,2),isMirror, isRelax, [], fProlloff);    % Filtering
+                Xnesting(i,:) = bst_bandpass_hfilter(X(i,:), sRate, bandNesting(i,1), bandNesting(i,2),isMirror, isRelax, [], fProlloff, Method);    % Filtering
             end
         end        
         Xnesting = Xnesting(:,nMargin-nHilMar+1:fix((margin+winLen)*sRate)+nHilMar);              % Removing part of the margin