comparison aux/autoalign.m @ 3:1f7b986dab05

Added autoalign.m: aligns audio by zero-padding before the signal to align the peak of the crosscorrelation function

-function autoalign(folderName, masterIndex)
-% AUTOALIGN aligns audio in folder with selected file (index 'masterIndex')
-% by inserting zeros at the start of the other files. 
+function autoalign(foldername)
+% AUTOALIGN aligns audio in folder with selected file by inserting zeros
+%  at the start of the files. 
 %
 % by Brecht De Man at Centre for Digital Music on 17 November 2014
 
-    error('Work in progress, function not finished');
+    % error('Work in progress, function not finished'); % while unfinished
     
     tic;                                % start measuring time
-    list = dir([foldername '\*.wav']);  % get names of files
+    list = dir([foldername '/*.wav']);  % get names of files
     
     % remove hidden files from list
     % see http://www.mathworks.co.uk/matlabcentral/newsreader/view_thread/258220
     for k = length(list):-1:1
         fname = list(k).name;
         if fname(1) == '.'
             list(k) = [ ];
         end
     end
-
-    % load master file
-    masterAudioChannels = audioread([foldername '/' list(masterIndex).name]); 
+    
+    shiftvec = zeros(length(list), 1); 
+    
+    % first file: 
+    masterAudioChannels = audioread([foldername '/' list(1).name]); % read file
     masterAudio         = sum(masterAudioChannels, 2); % sum to mono
-    disp(['Aligning audio files with ' list(masterIndex).name '...']);
     
     % for each other file:
     for i = 1:length(list)
-        
         slaveAudioChannels = audioread([foldername '/' list(i).name]); % read file
         slaveAudio         = sum(slaveAudioChannels, 2); % sum to mono
         
-        % check sampling rate is the same
+        % check sampling rate is the same (or: express in terms of seconds)
         % ... 
-
+        
         % crosscorrelation function
-        xcorr(masterAudio, slaveAudio); 
+        [xc, timeaxis] = xcorr(masterAudio, slaveAudio);
+        [value, lag] = max(abs(xc));
+        shiftvec(i) = timeaxis(lag);
         
-        % add zeros to beginning 
+        % debugging plots
+%         figure;
+%         title(['Crosscorrelation file 1 vs ' num2str(i)]);
+%         subplot(2,1,1);
+%         plot(timeaxis, xc);
+%         xlabel('Samples');
+%         ylabel('Crosscorrelation');
+%         subplot(2,1,2);
+%         plot(timeaxis, abs(xc));
+%         xlabel('Samples');
+%         ylabel('Absolute value of crosscorrelation');
         
         % monitoring/debugging
         
     end
+    
+    % add absolute of largest negative lag worth of zeros to all; 
+    % then add 'shiftvec' value to that
+    maxlag = min(shiftvec); 
+    for i = 1:length(list)
+        [audio, fs] = audioread([foldername '/' list(i).name]);
+        audiowrite([foldername '/SHIFT' list(i).name], ... % remove SHIFT
+            [zeros(-maxlag + shiftvec(i), size(audio, 2)) ; audio], ...
+            fs, 'BitsPerSample', 24);
+    end
 
     elapsedTime = toc;  % stop measuring time
-    disp(['autoalign took ' elapsedTime ' seconds to align the files in folder ' ...
-        folderName ' with master file ' list(masterIndex).name]);
+    disp(['autoalign took ' num2str(elapsedTime) ...
+        ' seconds to align the files in folder ' foldername]);
 
 end