1.         %% 创建存放写出图像的文件夹
    2.         dirname =  'Result';
    3.         mkdir(dirname);
    6.     %% 绘制时间序列对比图
    7. %     fig = figure;
    8.     fig = figure('Visible','off'); % 不显示绘图图窗
    9.     set(gcf,'Position',get(0,'ScreenSize')); % 设置绘图为全屏绘制
    10.     set(gca,'FontSize',20); % 设置plot横纵轴数字字体大小，注意这是全局字体大小设置，需要在标题大小设置的前面，否则会覆盖
    12.     % 绘制 NDVI 的时间序列图
    13.     p1=plot(ts_16d,NDVI_16d_3_years(poi,:),':*y','MarkerSize',10,'MarkerFaceColor','w');
    14.     p1.LineWidth=2;
    16. %     % 标题和坐标轴区域
    17. %     title_name = ['Original vs Reconstruction(iHANTS, POI-',num2str(POI_FID(poi)+1),', ',LC_name,')'];
    18. %     title(title_name,'FontName','Times New Roman','FontSize',25);
    19. %         % 字体是Times New Roman，颜色是蓝色（'b'即blue），字体大小是20号。
    21.     xlim=[0,70]; ylim=[0,1];
    22.     set(gca,'xtick',0:2:70,'FontSize',15) % handles可以指定具体坐标轴的句柄
    23.     set(gca,'ytick',0:0.1:1,'FontSize',15) % handles可以指定具体坐标轴的句柄
    24.     axis([0 70 0 1]);%axis（[xmin xmax ymin ymax]）,设置绘制区域的范围
    25.     %% 这个放在后面设置字号，否则上面会被覆盖掉
    26.     xlabel('Time Series Index','FontSize',20,'FontWeight','bold');
    27.     ylabel('NDVI','FontSize',20,'FontWeight','bold'); %设置横纵轴名称及字体大小
    29.     % 绘制 NDVI 时间序列图
    30.       % 绘制 NDVI 时间序列图
    31.     hold on;
    32.     p2=plot(ts_16d,yr_16d_3_years1,'--.c','MarkerSize',8);
    33.     p2.LineWidth=2;
    34.     hold on;
    35.     p3=plot(ts_16d,yr_16d_3_years2,'--.g','MarkerSize',8);
    36.     p3.LineWidth=2;
    37.     hold on;
    38.     p4=plot(ts_16d,yr_16d_3_years3,'--.b','MarkerSize',8);
    39.     p4.LineWidth=2;
    40.     hold on;
    41.     p5=plot(ts_16d,yr_16d_3_years4,'--.m','MarkerSize',8);
    42.     p5.LineWidth=2;
    43.     hold on;
    44.     p6=plot(ts_16d,yr_16d_3_years5,'--.r','MarkerSize',8);
    45.     p6.LineWidth=2;
    46.     hold on;
    47. %     p7=plot(ts_16d,yr_16d_3_years6,'--.k','MarkerSize',8);
    48. %     p7.LineWidth=2;
    50.     hold on;
    51.     line([24 24],[0 0.9],'linestyle',':', 'Color','k', 'LineWidth', 3);
    53.     hold on;
    54.     line([46 46],[0 0.9],'linestyle',':', 'Color','k', 'LineWidth', 3);
    56.     %% 标注高质量点
    57.     for marker_index = 1:length(y_QC_3_years)
    58.         if y_QC_3_years(marker_index)==1
    59.             hold on;
    60.             p8=plot(ts_16d(marker_index),NDVI_16d_3_years(poi,marker_index),...
    61.                  'kp','MarkerSize',20,'LineWidth',2);
    62.         end
    63.     end
    65.     %% 绘制图例
    66. %     leg=legend({'16d (Original, 3 years)','16d (Reconstructed, 3 years)',...
    67. %         '16d (Original,1 year)','16d (Reconstructed, 1 year)'},'Location','northeast'); 
    68.     leg=legend([p1,p2,p3,p4,p5,p6,p8],{'Original','delta=0.1','delta=0.2',...
    69.          'delta=0.3', 'delta=0.4','delta=0.5','delta=0.6','High Quality Points'},...
    70.          'Location','northoutside','Orientation','horizontal');   
    71.     legend('boxoff'); % 不显示图例框
    72.     set(leg,'FontName','Times New Roman','FontSize',18,'FontWeight','normal')
    73.     %% 写出对比图像
    74.     img_name = ['.\',dirname,'\POI_',num2str(POI_index(poi)),'-',LC_name]; % 写出到dir_name文件夹中
    75.     print(img_name, '-dpng', '-r0'); % 保存到工作目录下，名字为img_name
    76. %     catch
    77. %         fprintf(['POI_',num2str(POI_FID(poi)+1),' has some problems!\n']);
    78. % end
    80. end
    81. %% 关闭绘图
    82. close all;