NSIDC定义的海冰相关概念

本文主要是介绍NSIDC定义的海冰相关概念，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

文章目录

- 相关概念
- Matlab绘图
- 结果展示

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 onplot([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 onscatter(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 onset(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 ontext([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