M_Map Lambert Projection

조회 수: 19 (최근 30일)
Bianka Markovic
Bianka Markovic 2021년 6월 24일
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개)

카테고리

Help CenterFile Exchange에서 Climate Science and Analysis에 대해 자세히 알아보기

제품


릴리스

R2019a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by