matlab滴m_map学习心得5

%% 读取文件�
LONG_MIN = 100;
LONG_MAX = 260;
LAT_MIN = 20;
LAT_MAX = 64;  

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);    
[row,col] = size(LON);

[sst,windLF,windMF,vapor,cloud,rain] = read_gmi_averaged_v8('f35_202101v8.2');    
sst_np = sst(lon_start:lon_end,lat_start:lat_end);
windlf_np = windLF(lon_start:lon_end,lat_start:lat_end);
windmf_np = windMF(lon_start:lon_end,lat_start:lat_end);
vapor_np = vapor(lon_start:lon_end,lat_start:lat_end);
cloud_np = cloud(lon_start:lon_end,lat_start:lat_end);
rain_np = rain(lon_start:lon_end,lat_start:lat_end);

% 创建文件
fid=fopen('文件.txt','w');   
fprintf(fid,'%4s  %7s  %7s  %7s  %7s  %7s  %7s\r\n','Lat','SST','WSPD_LF','WSPD_MF','VAPOR','CLOUD','RAIN');
fprintf(fid,'\n');
for i=1:4:row
    sst_txt = sst_np(:,i:i+3);
    windlf_txt = windlf_np(:,i:i+3);
    windmf_txt = windlf_np(:,i:i+3);
    vapor_txt = vapor_np(:,i:i+3);
    cloud_txt = cloud_np(:,i:i+3);
    rain_txt = rain_np(:,i:i+3);
    index = find(sst_txt<=250);        
    fprintf(fid,'%4d  %7.1f  %7.1f  %7.1f  %7.1f  %7.1f  %7.1f\n',round(lat(i)),mean(sst_txt(index)),mean(windlf_txt(index)),mean(windmf_txt(index)),mean(vapor_txt(index)),mean(cloud_txt(index)),mean(rain_txt(index)));
end
fclose(fid);

%% 求平均值
file=dir('f35_*');
nums = length(file);
for i=1:nums
    [sst,windLF,windMF,vapor,cloud,rain] = read_gmi_averaged_v8(file(i).name);
    sst_t(:,:,i) = sst(lon_start:lon_end,lat_start:lat_end);
end
for i=1:col
    for j=1:row
        if all(sst_t(i,j,:)<=250)   
            sst_annual(i,j) = mean(sst_t(i,j,:));
        else
            sst_annual(i,j) = NaN;
        end
    end
end
f1 = figure;
m_proj('miller','lon',[100 260],'lat',[20 64]); 
m_pcolor(LON,LAT,sst_annual');
m_coast('patch',[0.7 0.7 0.7]);
m_gshhs_l('color','k');
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',[0,25],'color','none');     
ax=m_contfbar(1.05,[.15 .85],[0 30],[0:0.01:30],'edgecolor','none','endpiece','no','ticklength',[.01 0.05]);
xlabel(ax,'(^oC)','color','k','fontweight','bold','fontname','Times New Roman'); 
set(ax,'ytick',0:5:30,'yticklabel',{'0','5','10','15','20','25','30'})   
title(ax,'SST','color','k','fontweight','bold','fontname','Times New Roman')    
m_ruler([0.02 0.22],.92,4,'fontsize',10,'ticklength',[.007 .007],'tickdir','out','fontname','Times New Roman');
m_northarrow(240,58,6,'type',4,'linewi',0.5); 
title('GMI Sea Surface Temperature: 2021/01 - 2021/12 - North Pacific','fontweight','bold','fontname','Times New Roman','fontsize',16)
    

sst_monthly = zeros(12,1);
for i=1:12
        index = find(sst_t(:,:,i)<251);
        sst_monthly(i) = mean(sst_t(index));
end
f2 = figure;
x=1:12;
plot(x,sst_monthly,'Marker','o','LineWidth',2,'Color',[1 0 0]); 
    set(gca,'Xtick',[1:1:12],'Ytick',[14:2:26])
    xlabel('月份','FontSize',15);
    ylabel('海水表面温度（℃）','FontSize',15);