Make a plt from a 4-D multi array variable from nc file
显示 更早的评论
I have a nc file with:
Dimensions:
time = 8
depth = 48
latitude = 109
longitude = 25
and make this:
u='global-reanalysis-phy-001-030-daily_1591804259623.nc'
ncdisp(u)
east = ncread(u,'uo');
lat = ncread(u,'latitude');
lon = ncread(u,'longitude');
time=ncread(u, 'time');
depth=ncread(u, 'depth');
dtime = datetime(1950, 1, 1, time, 0, 0);
jd=juliandate(dtime);
plot(lat, east(:,:,15,1));
east is a 4-D 8371200 double
But I want a plot (surfc) in time, this 8 days, latitude, one depth (lets say 48, the last) and the east variable that is a velocity m/s.
Any idea? What I should do with east variable? squeeze, quiver, etc realy dont know or find solution
Thanks.
1 个评论
KSSV
2020-6-17
OP commented:
the plot(lat, east(:,:,15,1)) is working but only for one day if I change 1 for : then the error
Data cannot have more than 2 dimensions.
采纳的回答
ncfile = "global-reanalysis-phy-001-030-daily_1591804259623.nc";
time = ncread(ncfile,'time') ;
depth = ncread(ncfile,'depth') ;
lat = ncread(ncfile,'latitude') ;
lon = ncread(ncfile,'longitude') ;
u = ncread(ncfile,'uo') ;
% plots
for i = 1:length(time) % time loop
for j = 1:length(depth) % depth loop
z = u(:,:,j,i)' ; % extract the matrix here
pcolor(lon,lat,z) ;
colorbar
shading interp
xlabel('lon')
ylabel('lat')
title(sprintf('time = %d, depth = %d',i,j))
drawnow
end
end
17 个评论
Pepe Grillo
2020-6-17
KSSV
2020-6-17
try transposing the matrix..
surf(lon,lat,east(:,:,1)')
If not working share the dimensions of each variable here..
Pepe Grillo
2020-6-17
KSSV
2020-6-17
Edited the answer.
Pepe Grillo
2020-6-17
Pepe Grillo
2020-6-17
KSSV
2020-6-17
All the matrices are not NaN....when depth is increasing....data is having NaN's....see the plot, it is in a loop..it is plotting all the values..
Pepe Grillo
2020-6-17
KSSV
2020-6-17
Okay..I got it..you want to plot for specific depth, all the time data...
For this case, you need to fix the range of lat and lon...
Pepe Grillo
2020-6-17
KSSV
2020-6-17
ncfile = "global-reanalysis-phy-001-030-daily_1591804259623.nc";
time = ncread(ncfile,'time') ;
depth = ncread(ncfile,'depth') ;
lat = ncread(ncfile,'latitude') ;
lon = ncread(ncfile,'longitude') ;
u = ncread(ncfile,'uo') ;
% plots
for i = 1:length(lon)
ud = squeeze(u(:,i,1,:)) ;
pcolor(lon,time,ud')
shading interp ;
colorbar
title(sprintf("depth = 1, lat rows = %d",i)) ;
drawnow
pause
end
KSSV
2020-6-17
This is called Hovmoller diagrams ..am I right?
Pepe Grillo
2020-6-17
Pepe Grillo
2020-6-17
Pepe Grillo
2020-6-17
KSSV
2020-6-17
You can try it your self......check the dimensions before and after squeeze.
squeeze : it will remove the unwanted extra dimensions
Pepe Grillo
2020-6-17
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Graphics Object Properties 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
