quality plot for shift correction

functions/CellSortLoadTrialsShiftCorrOnTrials.m

@@ -51,4 +51,13 @@
 end
 
 disp('Finished Loading Trial Files!');
+%calculate quality plot
+reference = mean(data_r_array, 3);
+qualityPlot = zeros(1,size(data_r_array, 3));
+parfor i = 1:size(data_r_array,3)
+    qualityPlot(i) = squeeze(mean(mean(data_r_array(:,:,i) - reference,1),2));
+end
+qualityPlot = qualityPlot.^2;
+figure;
+plot(qualityPlot);
 end

main.m

@@ -11,7 +11,7 @@
 CellSortXCorr(inFolder, outFolder, nSubIm, binningFactor);
 %% load files
 %go to corr folder
-outFolder = 'F:\Arbeit\testData\PCAICA';
+outFolder = 'C:\Users\ftvol\Desktop\data';
 %[data_g, data_r, frames, fn_green, fn_red, fileNamesGreen, fileNamesRed] = CellSortLoadTrials(outFolder);
 [data_g, data_r, frames, fn_green, fn_red, fileNamesGreen, fileNamesRed] = CellSortLoadTrialsShiftCorrOnTrials(outFolder);
 %If you want to load single files, use the functions below.