# Create a map with 3 different variables

19 views (last 30 days)
Nicholas Engle on 13 Apr 2020
Edited: Adam Danz on 14 Apr 2020
I am trying to create a table that has bins of RPM on the x-axis, bins of throttle position on the Y-axis, and filled in the middle based on these valuse is the fuel used (z-axis). The fuel used vairable needs to be the average of all fuel values that fall in the RPM range and throttle range it falls in and place that average value in the corresping cell.
The code I have so far to try and do this is attached as well as the data file I am using. Also a screenshot of the table I am trying to recreate is attached.

Adam Danz on 13 Apr 2020
Edited: Adam Danz on 13 Apr 2020
Plotting bivariate binned means using heatmap
% Compute binned mean
[N,~,~,binX,binY] = histcounts2(RPM, Throttle, binEdgesRPM, binEdgesThrottle)
ind = sub2ind(size(N),binX,binY);
binMeans = reshape(arrayfun(@(i)mean(Fuel(ind==i)), 1:numel(N)),size(N));
% Compute bin centers (for heatmap)
binCentersRPM = binEdgesRPM(2:end)-(binEdgesRPM(2)-binEdgesRPM(1))/2;
binCentersThrottle = binEdgesThrottle(2:end)-(binEdgesThrottle(2)-binEdgesThrottle(1))/2;
% Plot heatmap
heatmap(binCentersRPM, binCentersThrottle, binMeans.')
... using imagesc
imagesc(binEdgesRPM, binEdgesThrottle, binMeans')
Plotting bivariates binned counts
histogram2() or histcounts2() + heatmap()
filename = 'Scarface_tune_run2.csv';
%% Throttle Base Injection Qty Map
RPM = Raw(:,4);
Throttle= Raw(:,5);
Fuel = Raw(:,11);
Map1 = cat(2,RPM,Throttle,Fuel);
topEdge = 4800;
botEdge = 0;
numBins = 24;
binEdgesRPM = linspace(botEdge, topEdge, numBins+1);
binEdgesThrottle = linspace(0, 100, 11);
% THIS LINE BELOW PRODUCES THE PLOT
histogram2(RPM, Throttle, binEdgesRPM, binEdgesThrottle, 'DisplayStyle','tile','ShowEmptyBins','on')
Or
N = histcounts2(RPM, Throttle, binEdgesRPM, binEdgesThrottle)
heatmap(binEdgesRPM, binEdgesThrottle,N)
Adam Danz on 14 Apr 2020
Add 1 column to the end showing the row-means; and 1 row to the end showing the column means. After plotting the heatmap, change the last x and y label to "mean".
binMeans(end+1,:) = mean(binMeans,1);
binMeans(:,end+1) = mean(binMeans,2);
binMeans(end,end) = NaN;
binCentersRPM(end+1) = binCentersRPM(end) + binEdgesRPM(2)-binEdgesRPM(1);
binCentersThrottle(end+1) = binEdgesThrottle(end) + binEdgesThrottle(2)-binEdgesThrottle(1);
hm = heatmap(binCentersRPM, binCentersThrottle, binMeans.');
hm.XDisplayLabels(end) = {'mean'};
hm.YDisplayLabels(end) = {'mean'};