Hovmoller Diagram: longitude (x axis), time (y axis), and z(temperature values)

Question

Carlotta Dentico 2023-9-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2015806-hovmoller-diagram-longitude-x-axis-time-y-axis-and-z-temperature-values

评论： Mathieu NOE 2023-9-20

Hi,

I wuold like to extract temperature from a 3D matrix temp 609x881x372 (lonxlatxtime) and do an Hovmoller plot where I have time on the y-axis and longitude on the x-axis (see figure).

I created a logical matrix of 0 and 1 to identify where my coordinates meet my condition:

coord = (aa_lon>=-7 & aa_lon<=9 & aa_lat>=79 & aa_lat<=79.01);

coord is 609x881 logical but my matrix (temp) from which I want to extract the data is 609x881x372 double.

Hope my question is clear.

Please help me :)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Mathieu NOE 2023-9-13

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2015806-hovmoller-diagram-longitude-x-axis-time-y-axis-and-z-temperature-values#answer_1309251

移动：Mathieu NOE 2023-9-13

在 MATLAB Online 中打开

hello again

I modified my code so it should work now on your data

as I understand we have here only 12 months of data , so my display assumes we have here only the first 12 months (1 year) of the entire time vector (which starts in Jan 1991) - you may need to correct this

also I changed a bit the tolerance on the latitude : aa_lat>=79 & aa_lat<=79+0.1

so we pick more data points of the temp matrix ; otherwise your diagram will have a very coarse display

result so far :

code

%% load data 609x881x12  (lonxlatxtime)
load lat.mat
load lon.mat
load time.mat 
load temp.mat  
%% convert time serial number into Y/M/D format
[y,mo,d,h,mi,s] = datevec(TT);% samplin time = 1 month (data is stored every 15th of the month)
%% spatial filtering 
% coord = (aa_lon>=-7 & aa_lon<=9 & aa_lat>=79 & aa_lat<=79.01);
b = (aa_lat>=79 & aa_lat<=79+0.1); % play with tolerance (here 0.1) and see impact on how many points are selected
selected_lon = aa_lon(b);
ind_lon = (selected_lon>=-7 & selected_lon<=9);
final_lon = selected_lon(ind_lon);
%% main code
[m,n,p] = size(temp);
out = [];
for k = 1:p
    tmp = temp(:,:,k);
    selected_temp = tmp(b); % filtering for lat
    selected_temp = selected_temp(ind_lon); % filtering for lon
    % concatenation of the selected_temp data  
    out = [out; selected_temp']; % nb transpose of selected_temp to have data as row 
end
% plot 
% downsample Y tick spacing to display only years (tick position at january
% / 1st month of the year)
new_y = (1:12:p);
imagesc(final_lon,(1:p),out);
yticks(new_y)
yticklabels(num2str(y(new_y)))
xlabel('Longitude');
ylabel('Time (Years)');
% add horizontal line corresponding to january 15th (first data of the
% year)
hold on
for ci = 1:numel(new_y)
    plot(final_lon,new_y(ci)*ones(size(final_lon)),'k--');
end