%% AMSR-E
% 读入SIC
sic = hdfread('POGOC_OUC_AMSR-E_DTASI_n_sic_YYYYMMDD.hdf','OUC_SIC');
sic=sic*100;
% 读入经纬度网格
flat = fopen('psn06lats_v3.dat','r');
polarlat = fread(flat,[608*2 896*2],'int32');
fclose(flat); clear flat;
lat = (polarlat/100000)';
flon = fopen('psn06lons_v3.dat','r');
polarlon = fread(flon,[608*2 896*2],'int32');
fclose(flon); clear flon;
lon = (polarlon/100000)';

% 画图(需提前导入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/AMSR (6.25 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')
%% AMSR-2
% 读入SIC(2019年之后文件名变化)
sic =hdfread('POGOC_OUC_AMSR2_DTASI_n_sic_YYYYMMDD.hdf','OUC_SIC'); %2012-2018
sic =hdfread('OUC_AMSR2_DTASI_n_sic_YYYYMMDD.hdf','OUC_SIC'); %2019-
sic=sic*100;
% 读入经纬度网格
 lat = hdfread('LongitudeLatitudeGrid-n6250-Arctic.hdf','Latitudes');
 lon = hdfread('LongitudeLatitudeGrid-n6250-Arctic.hdf','Longitudes');
% 画图(需提前导入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/AMSR (6.25 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')
