matlab滴m_map学习心得6

乾坤毅志

已于 2024-03-29 17:02:44 修改

阅读量94

点赞数 1

分类专栏： m_map实战文章标签： 1024程序员节 matlab 算法学习方法学习经验分享

于 2023-10-24 16:10:32 首次发布

本文链接：https://blog.csdn.net/qq_58518586/article/details/134014945

版权

m_map实战专栏收录该内容

8 篇文章 1 订阅

订阅专栏

%% 读取文件�
LONG_MIN = 117;
LONG_MAX = 127;
LAT_MIN = 31;
LAT_MAX = 41;  

latstep=0.25*(1:720)-90.125;
lonstep=0.25*(1:1440)-0.125;    
lat_all = find(latstep>=LAT_MIN &latstep <=LAT_MAX);
lat_start = lat_all(1);
lat_end = lat_all(end);
lon_all = find(lonstep>=LONG_MIN & lonstep<=LONG_MAX);
lon_start = lon_all(1);
lon_end = lon_all(end);
lon = lonstep(lon_start:lon_end);
lat = latstep(lat_start:lat_end);    
[LON,LAT] = meshgrid(lon,lat);    
%% 循环读取
file=dir('ascatc_*');
nums = length(file);
for i=1:nums
    [windspd,winddir,scatflag,radrain,sosmap]=get_ascat_averaged(file(i).name);
    windspd_t(:,:,i) = windspd(lon_start:lon_end,lat_start:lat_end);
    winddir_t(:,:,i) = winddir(lon_start:lon_end,lat_start:lat_end);
end
index=find(windspd_t==-999);
windspd_t(index)=NaN;
index=find(winddir_t==-999);
winddir_t(index)=NaN;
%% 十二个月可视化
for i=1:12
    fi = figure;
    m_proj('miller','lon',[117 127],'lat',[31 41]); 
    m_pcolor(LON,LAT,windspd_t(:,:,i)');
    m_gshhs_l('color','k');
    m_gshhs_l('patch',[0.7 0.7 0.7]);
    m_grid('tickstyle','dd','tickdir','in','linestyle','none','fontsize',12,'backgroundcolor','none','fontname','Times New Roman','FontWeight','bold', 'linewidth',1,'ticklength',.01);    
    map=jet(199);
    colormap(map); 
    set(gca,'CLim',[2,10],'color','none');     
    ax=m_contfbar(1.05,[.15 .85],[0 30],[2:0.01:10],'edgecolor','none','endpiece','no','ticklength',[.01 0.05]);
    xlabel(ax,'(m/s)','color','k','fontweight','bold','fontname','Times New Roman'); 
    set(ax,'ytick',2:2:10,'yticklabel',{'2','4','6','8','10'})     
    m_ruler([0.02 0.22],.92,2,'fontsize',5,'ticklength',[.007 .007],'tickdir','out','fontname','Times New Roman');
    m_northarrow(126,40,1,'type',4,'linewi',0.5); 
    title(['2021/0' num2str(i)],'fontweight','bold','fontname','Times New Roman','fontsize',10);
    hold on
    u=(windspd_t(:,:,i).*(sin(winddir_t(:,:,i)*pi/180)))';
    v=(windspd_t(:,:,i).*(cos(winddir_t(:,:,i)*pi/180)))';
    m_quiver(LON,LAT,u/30,v/30,0,'k');
    m_quiver(124,40.5,10/30,0,0,'r','MaxHeadSize',5);
    m_text(124,40.6,'10m/s','color','r','fontsize',8);
end