相关概念

• NSIDC 表示 “National Snow and Ice Data Center”，即美国国家雪和冰数据中心。NSIDC 是一个位于美国科罗拉多大学波尔得分校的研究中心，致力于收集、管理和分发全球雪和冰的科学数据。

Matlab绘图

clear;clc;close all
load('.\input\CDR\Daily\sic_nt_daily_7920.mat')
load('.\data\input\BFT sea\BFT_mask.mat')

data00=squeeze(nanmean(nanmean(nanmean(sic_nt_daily_7920(:,:,:,[20 30 34 38 41]).*mask_sea,1),2),4));
data01=smooth(data00,3);
m1=find(data01<0.8);
m2=find(data01<0.15);

x_0=0.1;
y_0=0.56;
len=0.85;
width=0.30;
d_x=0.43;
d_y=-0.30;
px=[0 0 0 1];
py=[0 1 1 1];
mon=[31 28 31 30 31 30 31 31 30 31 30 31];
mona=[1,32,60,91,121,152,182,213,244,274,305,335];
monb=[31,59,90,120,151,181,212,243,273,304,334,365];
for ii=1:length(mon)
    a(ii)=sum(mon(1:ii));
end
aa=a-mon/2;
month={'Jan.','Feb.','Mar.','Apr.','May','Jun.','Jul.','Aug.','Sep.','Oct.','Nov.','Dec.'};

siz=25;lind=1.5;lind1=3;
lw=2;
ms=18;

set(gcf,'color',[1 1 1],'position',[10 45 800 800*1.2]);%get(0,'screensize')
axes('position',[x_0+d_x*px(1), y_0+d_y*py(1), len, width]);

plot(data01,'k-','linewidth',lind+0.5);hold on
plot([1 365],[0.8 0.8],'color',[210 82 60]/255,'linestyle','--','linewidth',lind);hold on
plot([1 365],[0.15 0.15],'color',[36 53 245]/255,'linestyle','--','linewidth',lind);hold on
plot([m1(1) m1(1)],[0 data01(m1(1))],'color',[210 82 60]/255,'linestyle','--','linewidth',lind);hold on
plot([m1(end) m1(end)],[0 data01(m1(end))],'color',[210 82 60]/255,'linestyle','--','linewidth',lind);hold on

plot([m2(1) m2(1)],[0 data01(m2(1))],'color',[210 82 60]/255,'linestyle','--','linewidth',lind);hold on
plot([m2(end) m2(end)],[0 data01(m2(end))],'color',[210 82 60]/255,'linestyle','--','linewidth',lind);hold on

scatter(m1(1),data01(m1(1)),80,'filled','Marker','^','MarkerFaceColor',[210 82 60]/255);hold on
scatter(m1(end),data01(m1(end)),80,'filled','Marker','^','MarkerFaceColor',[210 82 60]/255);hold on
scatter(m2(1),data01(m2(1)),80,'filled','Marker','^','MarkerFaceColor',[36 53 245]/255);hold on
scatter(m2(end),data01(m2(end)),80,'filled','Marker','^','MarkerFaceColor',[36 53 245]/255);hold on

set(gca,'linewidth',1.5);grid on
set(gca,'xlim',[1 365],'xtick',a,'xticklabel',month,'fontname','Times New Roman','fontweight','bold','fontsize',siz-10);
set(gca,'ylim',[0 1],'ytick',[0:0.2:1.0],'yticklabel',num2str([0:20:100]'),'fontname','Times New Roman','fontweight','bold','fontsize',siz-10);
ylabel('Sea Ice Concentrration (%)','fontname','Times New Roman','fontweight','bold','fontsize',siz-10);

hold on
text(m1(1)-25,data01(m1(1))-0.05,'\it{DOO}','fontname','Times New Roman','fontweight','bold','fontsize',siz-10);
text(m2(1)-25,data01(m2(1))-0.05,'\it{DOR}','fontname','Times New Roman','fontweight','bold','fontsize',siz-10);
text(m2(end)+5,data01(m2(end))-0.05,'\it{DOA}','fontname','Times New Roman','fontweight','bold','fontsize',siz-10);
text(m1(end)+5,data01(m1(end))-0.05,'\it{DOC}','fontname','Times New Roman','fontweight','bold','fontsize',siz-10);

drawbrace([m1(1) 0.15],[m2(1) 0.15],10,'color',[160 32 240]/255,'linewidth',1.5);hold on
drawbrace([m2(end) 0.15],[m1(end) 0.15],10,'color',[160 32 240]/255,'linewidth',1.5);hold on
drawbrace([m2(1) 0.15],[m2(end) 0.15],10,'color',[30 144 255]/255,'linewidth',1.5);hold on
drawbrace([m1(1) 0.8],[m1(end) 0.8],10,'color',[30 144 255]/255,'linewidth',1.5);hold on

text([m1(1)+m2(1)]/2.2,0.25,'\it{SLIP}','color',[160 32 240]/255,'fontname','Times New Roman','fontweight','bold','fontsize',siz-10)
text([m1(end)+m2(end)]/2.1,0.25,'\it{SGIP}','color',[160 32 240]/255,'fontname','Times New Roman','fontweight','bold','fontsize',siz-10)
text([m2(1)+m2(end)]/2.1,0.25,'\it{IITP}','color',[30 144 255]/255,'fontname','Times New Roman','fontweight','bold','fontsize',siz-10)
text([m1(1)+m1(end)]/2.1,0.9,'\it{OIFP}','color',[30 144 255]/255,'fontname','Times New Roman','fontweight','bold','fontsize',siz-10)

text(1,-0.5,{'DOO(Day of opening): last day sea ice concentration drops below 80%',...
    'DOR(Day of retreat): last day sea ice concentration drops below 15%',...
    'DOA(First day of advance): first day sea ice concentration increases above 15%',...
    'DOC(First day of closing): first day sea ice concentration increases above 80%',...
    'OIFP(Outer ice-free period): DOC - DOO',...
    'IIFP(Inner ice-free period or open-water period): DOA - DOR',...
    'SLIP(Sensonal loss of ice period): DOR - DOO',...
    'SGIP(Seasonal gain of ice period): DOC - DOA'},...
    'color',[205 85 85]/255,'fontname','Times New Roman','fontweight','bold','fontsize',siz-10)

export_fig(['.\','NSIDC_海冰相关定义.png'],'-r300')
close all

结果展示