% 读入SIC
sic = hdfread('OUC_FY3D_MWRI_DTASI_n_sic_YYYYMMDD.hdf','OUC SIC');
sic=sic*100;
% 读入经纬度网格
fid = fopen('psn12lats_v3.txt','r');
sizeA = [896 608];
lat = fscanf(fid,'%f', sizeA);
fclose(fid);
fid = fopen('psn12lons_v3.txt','r');
sizeA = [896 608];
lon = fscanf(fid,'%f', sizeA);
fclose(fid);    
% 画图(需提前导入m_map工具箱)
figure;
set(gcf,'Position',[60 60 800 600],'paperpositionmode','auto');
m_proj('Stereographic','lat',90,'lon',30,'radius',35);
m_pcolor(lon,lat,sic) ;
m_grid;
m_gshhs_l('save','mycoast');
m_usercoast('mycoast','patch',[.8 .8 .8]);
title('DT-ASI/FY (12.5 km)','FontSize',15);
caxis([0 100]) 
hc=colorbar;Pc=get(hc,'Position');
set(hc,'Position',[Pc(1)+0.088 Pc(2) Pc(3) Pc(4)],'FontSize',13,'FontWeight','b');
ylabel(hc,'SIC: %','FontSize',13,'FontWeight','b')
