M_Map Lambert Projection
9 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
Hello everyone,
I'm currently trying to plot a lambert projection of my climate data. I have temperature values on a regular grid. The dimensions are 1920x374 (for Lon_oras,Lat_oras,an_temp). The dimension of longitude is 1920 and latitude 374.
I used the projection from M_Map which is available on this site: M_Map: A Mapping package for Matlab (ubc.ca)
Unfortunately, somehow the program won't run through. I'm not sure if this is because of the amount of data or is something wrong in my code. Woud love to hear some suggestions what could be the culprit of this problem. Thank you!
Bianka
[Lon_oras,Lat_oras] = meshgrid(lon_oras,lat_oras);
vq(:,:,cnt)=griddata(xsur,ysur,temp_layer,Lon_oras,Lat_oras);
%FESOM Mean temperature over the whole time period for reference
temp_mean_fesom=mean(vq_all,3);
%Annual mean
temp_annual_fesom=mean(vq,3);
%anomalie
an_temp_fesom=temp_mean_fesom-temp_annual_fesom;
m_proj('lambert','long',[-90 50],'lat',[45 85]);
h=figure;
m_pcolor(Lon_oras,Lat_oras,vq(:,:,1));
caxis([-2 10])
colormap;
cb=colorbar;
xlabel(cb,'Temperature');
m_grid('box','fancy','tickdir','in');
m_coast('patch',[.8 .8 .8],'edgecolor','k');
title(['year: ',num2str(year),'; depth: ',num2str(depth),' m'],'fontsize',16);
saveas(h, ['plots/',runid,'_temp_',int2str(year),'_depth_',int2str(depth),'.png']);
0 Kommentare
Antworten (0)
Siehe auch
Kategorien
Mehr zu Climate Science and Analysis finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!