Hello Nikolay,
The code below computes the density using histcounts2() and then plots the results using imagesc(). You'll have to play around with the 'nBins' parameter to control the density within each bin. It works on your data from Trial1hEM.mat.
nBins = 50;
minData = min([Trial1hEM;Trial1vEM]);
maxData = max([Trial1hEM;Trial1vEM]);
bins = linspace(minData, maxData, nBins);
N = histcounts2(Trial1vEM(:), Trial1hEM(:), bins, bins);
figure;
subplot(1,2,1)
scatter(Trial1hEM, Trial1vEM, 'b.');
axis equal;
xlim([minData,maxData])
ylim([minData,maxData])
% Plot heatmap:
sh = subplot(1,2,2);
imagesc(bins, bins, N);
axis equal;
xlim([minData,maxData])
ylim([minData,maxData])
set(sh, 'YDir', 'normal')
Uspeh! Adam
