diff --git a/CellsortApplyFilter.m b/CellsortApplyFilter.m index 13663d1..825db2c 100644 --- a/CellsortApplyFilter.m +++ b/CellsortApplyFilter.m @@ -44,7 +44,7 @@ fprintf('Loading %5g frames for %d ROIs.\n', nt, size(ica_segments,1)) tic -while k 0).*tempImage(:,:,i); -end -%scale tempImage -% maximum = max(tempImage(:)); -% minimum = min(tempImage(:)); -% tempImage = uint16(65535*(tempImage - minimum)/(maximum-minimum)); -tempImage = uint16(tempImage*65535); - -%make background white -% blackMask = uint16(all((tempImage < 1),3)); -% for i = 1:3 -% tempImage(:,:,i) = blackMask * 65535 + tempImage(:,:,i); -% end -meanImage = mean(data_r,3); -maximum = max(meanImage(:)); -minimum = min(meanImage(:)); -meanImage = (meanImage - minimum)/(maximum-minimum); -meanImage = uint16(imadjust(meanImage, [0, 1], [0, 1], 0.4)*65535); -%copy meanImage in every colorchannel -for i = 1:3 - tempImage(:,:,i) = tempImage(:,:,i) + meanImage; -end - -subplot(1,2,1); -imshow(tempImage); -%draw number in each roi -for drawNumberCt = 1:size(cluster{i}) - number = cluster{i}(drawNumberCt); -end -subplot(1,2,2); -hold on -cell_sigPlot = cell_sig(sortedIdx,:); - - - -signal = zeros(size(cluster,2),size(cell_sig,2)); -for i = 1:size(cluster, 2) - cell_sig_merged{i} = mean(cell_sig_merged{i}(cluster{i},:), 1); - for clusterIdx = cluster{i} - %subtract minimum value - minimum = min(cell_sigPlot(clusterIdx,:)); - cell_sigPlot(clusterIdx,:) = cell_sigPlot(clusterIdx,:) - minimum; - %norm cell signals - integral = median(cell_sigPlot(clusterIdx,:)); - cell_sigPlot(clusterIdx,:)=cell_sigPlot(clusterIdx,:)/abs(integral); - end - %mean value of all cluster signals - signal(i,:) = mean(cell_sigPlot(cluster{i},:), 1); - %add maximum position from signal from before - if i ~= 1 - %maximum = max(signal(i-1,:)); - maximum = max(max(cell_sigPlot(cluster{i-1},:))); - else - maximum = 0; - end - signal(i,:) = signal(i,:) + maximum; - rgb = cm(round(64/subPlotCols * (i-1)+1),:); - plot(signal(i,:), 'Color', rgb); - text(size(signal,2), maximum, num2str(i)) - %plot single components dotted - for clusterIdx = cluster{i} - cell_sigPlot(clusterIdx,:) = cell_sigPlot(clusterIdx,:) + maximum; - plot(cell_sigPlot(clusterIdx, :),':', 'Color', rgb); - end -end -xlim([0,size(signal,2)]); -%maximum = max(cell_sigPlot(:)); -ylim auto;%([0, maximum]); -end \ No newline at end of file diff --git a/alternativePlot.m b/alternativePlot.m index ebfcf7a..a20388e 100644 --- a/alternativePlot.m +++ b/alternativePlot.m @@ -55,7 +55,7 @@ %get position to draw number in mask = squeeze(ica_segments(number,:,:)); cog = centerofgravity(mask); - text(cog(2), cog(1), num2str(i), 'Color', 'white') + text(cog(2), cog(1), num2str(i), 'Color', 'white', 'FontSize', 15) end end