M_Map Lambert Projection

12 次查看（过去 30 天）

Bianka Markovic 2021-6-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/864235-m_map-lambert-projection

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']);