M_Map Lambert Projection
9 次查看(过去 30 天)
显示 更早的评论
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 个评论
回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Climate Science and Analysis 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!