How to make a Heatmap Calendar plot

Question

Emanuel Valdes 2018-8-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/416583-how-to-make-a-heatmap-calendar-plot

编辑： CdC 2021-4-26

data.mat

I'm trying to make a Heatmap Calendar for MP2.5 like the one made for R's "Open Air Project", using date and concentration (data attached)

http://www.openair-project.org/PDF/OpenAir_Manual.pdf

It should look like this:

I know there's a calendar function, but I have no idea of how to plot it. I would really appreciate any help you could give me. Thanks a lot!

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

jonas 2018-8-28

Should be possible but probably requires a lot of code unless you can find something on FEX. you may be able to overlay a heatmap on the calender function.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

jonas 2018-8-28

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/416583-how-to-make-a-heatmap-calendar-plot#answer_334460

编辑：jonas 2018-8-29

在 MATLAB Online 中打开

I made something similar to what you were asking for. It's just something I put together for fun, so there are definitely some shortcomings and lots of shortcuts. However, if you don't find anything better on FEX, then you can develop this further. You have to play around with the figure size, markersize and fontsize to make it look good. If you change one you have to adjust the others.

% Load data
data=load('cal.mat')
Date=data.Date;
MP25=data.MP25;
MP25(isnan(MP25))=-100;
% Group data
x=weekday(Date);
y=week(Date);
[G,m]=findgroups(month(Date));
% Figure starts here!
figure;hold on
for ii=1:12
% Extract subgroup of monthly data and name of month
xm=x(G==ii);
ym=y(G==ii);
zm=MP25(G==ii);
Mname = unique(month(Date(G==ii),'name'));
% Initiate subplot ii
ax(ii) = subplot(3,4,ii);hold on
% Write out day of month on grid
str=cellstr(num2str(day(Date(G==ii),'dayofmonth')));
text(xm,ym,ones(size(xm)).*100,str,'fontsize',6,...
    'fontweight','bold',...
    'verticalalignment','middle',...
    'horizontalalignment','center')
zm(zm<0)=NaN;
% Plot color scale
H(ii) = scatter3(xm,ym,zm,150,zm,...
    's','filled','markeredgecolor','none');
% Fix some properties that are annoying to change after plotting
title(Mname)
set(gca,'ylim',[min(ym)-.5 min(ym)+5.5],...
    'ytick',[min(ym)-.5:max(ym)+.5])
grid on
end
% Add colorbar
colorbar('west','position',[0.93 0.30 0.03 0.4])
% Set up properties for all axes
set(ax,'clim',[0 90],...
    'xcolor','k',...
    'ycolor','k',...
    'xtick',1.5:8,...
    'box','on',...
    'xticklabels',{'S     ','M     ','T     ','W     ','T     ','F     ','S     '},...
    'yticklabels',{},...
    'ticklength',[0 0],...
    'ydir','reverse',...
    'xlim',[0.5 7.5],...
    'view',[0 90],...
    'layer','top');
% Adjust figure size to fit with the font/markersize
set(gcf,'color',[1 1 1],...
    'position',[955   140   640   428])
colormap(flipud(autumn))

14 个评论
显示 12更早的评论隐藏 12更早的评论

jonas 2018-8-29

编辑：jonas 2018-8-29

在 MATLAB Online 中打开

Replace the previous code after the for-loop with this, and it should give you what you want.

I would try a different colormap with more distinct colors, perhaps 'jet'.

% Add colorbar
cb=colorbar('west','position',[0.85 0.30 0.02 0.4])
colors=flipud(jet(5));
cmap=[repmat(colors(1,:),5,1);repmat(colors(2,:),3,1);repmat(colors(3,:),3,1);repmat(colors(4,:),6,1);repmat(colors(5,:),3,1)]
colormap(cmap)
cb.Ticks=[50/2 50+30/2 80+30/2 110+60/2 180+20/2];
cb.TickLabels={'Good','Regular','Alert','Pre-Emergency','Emergency'};
% Set up properties for all axes
set(ax,'clim',[0 200],...
    'xcolor',[.7 .7 .7],...
    'ycolor',[.7 .7 .7],...
    'xtick',1.5:8,...
    'xticklabels',{'S     ','M     ','T     ','W     ','T     ','F     ','S     '},...
    'yticklabels',{},...
    'ticklength',[0 0],...
    'ydir','reverse',...
    'xlim',[0.5 7.5],...
    'view',[0 90],...
    'layer','top',...
    'box','on');
% Adjust figure size to fit with the font/markersize
set(gcf,'color',[1 1 1],...
    'position',[955   140   640   428])
for i=1:12;
    old_pos=ax(i).Position;
    ax(i).Position=old_pos-[0.1 0 0 0];
end

For segments of the colorbar are not equal in length, as you may want. This is more difficult to achieve. If it's important, I would probably make a 'fake' colorbar with patch objects, and hide the real one.

Emanuel Valdes 2018-8-29

编辑：Emanuel Valdes 2018-8-29

在 MATLAB Online 中打开

cal2.mat

I think it's more a inssue about how to set the intervals for each color of the colormap.

I've been trying to apply this example to a new dataset (cal2.mat): https://www.mathworks.com/matlabcentral/answers/275085-how-to-have-a-discrete-color-bar-for-scatter3

yet I'm having some problems...

 % Load data
 data=load('cal2.mat')
 Date=data.Date;
 MP25=data.MP25;
 MP25(isnan(MP25))=-100;
 % Group data
 x=weekday(Date);
 y=week(Date);
 [G,m]=findgroups(month(Date));
 % Figure starts here!
 figure;hold on
 for ii=1:12
 % Extract subgroup of monthly data and name of month
 xm=x(G==ii);
 ym=y(G==ii);
 zm=MP25(G==ii);
 Mname = unique(month(Date(G==ii),'name'));
 % Initiate subplot ii
 ax(ii) = subplot(3,4,ii);hold on
 % Write out day of month on grid
 str=cellstr(num2str(day(Date(G==ii),'dayofmonth')));
 text(xm,ym,ones(size(xm)).*100,str,'fontsize',6,...
     'fontweight','bold',...
     'verticalalignment','middle',...
     'horizontalalignment','center')
 zm(zm<0)=NaN;
 % Plot color scale
 H(ii) = scatter3(xm,ym,zm,150,zm,...
     's','filled','markeredgecolor','none');
 % Fix some properties that are annoying to change after plotting
 title(Mname)
 set(gca,'ylim',[min(ym)-.5 min(ym)+5.5],...
     'ytick',[min(ym)-.5:max(ym)+.5])
 grid on
 end
 % Add colorbar
 colorbar('west','position',[0.93 0.30 0.03 0.4],...
          'Ticks',[30,60,90,120,180],...
          'TickLabels',{'Good','Regular','Alert','Pre- Emercengy','Emergency'})
 % Set up properties for all axes
 set(ax,'clim',[0 90],...
     'xcolor','k',...
     'ycolor','k',...
     'xtick',1.5:8,...
     'box','on',...
     'xticklabels',{'M     ','T     ','W     ','T     ','F     ','S     
 ','S     '},...
     'yticklabels',{},...
     'ticklength',[0 0],...
     'ydir','reverse',...
     'xlim',[0.5 7.5],...
     'view',[0 90],...
     'layer','top');
 % Adjust figure size to fit with the font/markersize
 set(gcf,'color',[1 1 1],...
     'position',[955   140   640   428]); 
 %%colormap([0 1 0; 1 1 0; 1.0000 0.6000 0; 1 0 0; 0.8000 0 1.0000])
 crange = caxis;
 t = linspace(crange(1),crange(2),64);
 cmap = zeros(64,3);
 %Good
 mask = t<=50;
 cmap(mask,:) = repmat([0 1 0],[sum(mask), 1]);
 %Regular
 mask = t>=51 & t<=79;
 cmap(mask,:) = repmat([1 1 0],[sum(mask), 1]);
 %Alert
 mask = t>=80 & t<=109;
 cmap(mask,:) = repmat([1.0000 0.6000 0],[sum(mask), 1]);
 %Pre-Emergency
 mask = t>=110 & t<=169;
 cmap(mask,:) = repmat([1 0 0],[sum(mask), 1]);
 %Emergency
 mask = t>=170;
 cmap(mask,:) = repmat([0.8000 0 1.0000],[sum(mask), 1]);
 colormap(cmap); 
 clear all

jonas 2018-8-29

编辑：jonas 2018-8-29

在 MATLAB Online 中打开

matlabhelp3.m

Other than the fact that some numbers seem to be missing, I don't see any major problem. You have to change the 'clim' of the axes (currently [0 90]) if you want to see the entire range. I will take a look.

Also, it seems you are still using US/UK week days, but you changed the labels. This makes the labels wrong.

EDIT: Okay, the fact that I some numbers are not being written is really weird. I can't wrap my head around it... There is no error, warning or anything. Now my old code is not working either....!?

Edit2: The code works for me if I put the following line after the scatterplot:

plot3(1,1,1000)

This additional plot prompts matlab to draw the text. Again, I cannot wrap my head around why this happens and I've tested the usual drawnow, refreshdata, pause... The text is there but invisible for some reason. I have checked the properties and everything looks normal. I have attached the m-file to this comment so that you can try it yourself.

jonas 2018-8-31

在 MATLAB Online 中打开

No problemo. What do you mean month name local?

This is where the month name is extracted and stored in a cell array:

Mname = unique(month(Date(G==ii),'name'));

you can change the 'name' to 'shortname' ( Jan,Feb ... ) or 'monthofyear' ( 1,2,...), but that's it. You can of course also remove this line and specify any strings manually:

Mname={'Janu','Febr'...,'Dece');

CdC 2021-4-26