1. %% 绘图
    2. % fig = figure;% 新建一个figure,并将图像句柄保存到fig
    3. fig = figure('Visible','off'); % 不显示figure
    4. set(gcf,'Position',get(0,'ScreenSize')); % 设置绘图为全屏绘制
    5. set(gca,'FontName','Times New Roman','FontSize',15); % 设置坐标轴的字体和大小
    6. p1=plot(ts,y,':xk','MarkerSize',12,'MarkerFaceColor','k');
    7. p1.LineWidth=2;
    8. % title_name=['Remove One or Multi points in one iteration','(POI ',num2str(i),')'];
    9. % title(title_name);
    10. xlabel('Day Of Year (2019)','FontSize',20);
    11. ylabel('NDVI','FontSize',20);
    12. xlim=[0,366];ylim=[0,1];
    13. set(gca,'xtick',1:16:353,'FontSize',15) % handles可以指定具体坐标轴的句柄
    14. set(gca,'ytick',0:0.1:1,'FontSize',15) % handles可以指定具体坐标轴的句柄
    15. axis([0 366 0 1]);%axis([xmin xmax ymin ymax]),设置绘制区域的范围
    16. %% HANTS重建
    17. hold on;
    18. p2=plot(ts,yr1,'--*g','MarkerSize',14,'MarkerFaceColor','g');
    19. p2.LineWidth=3;
    20. %% iHANTS重建
    21. hold on;
    22. p3=plot(ts,yr2,'-om','MarkerSize',10,'MarkerFaceColor','m');
    23. p3.LineWidth=1.5;
    24. %% 标注参与拟合的点(即权重为1的点)
    25. for marker_index = 1:ni
    26. if p_final_1(marker_index)==1
    27. hold on;
    28. p4=plot(ts(marker_index),y(marker_index),...
    29. 'bs','MarkerSize',18,'LineWidth',2,'MarkerEdgeColor','b');
    30. end
    31. if p_final_2(marker_index)==1
    32. hold on;
    33. p5=plot(ts(marker_index),y(marker_index),...
    34. 'rp','MarkerSize',18,'LineWidth',2,'MarkerEdgeColor','r');
    35. end
    36. end
    37. %% 绘制图例
    38. leg = legend([p1,p2,p3,p4,p5],{'Original',...
    39. 'Reconstructed (HANTS-revised)','Reconstructed (iHANTS-revised)',...
    40. 'Fitting points (HANTS-revised)','Fitting points (iHANTS-revised)'},...
    41. 'Location','best');
    42. set(leg,'FontName','Times New Roman','FontSize',18,'FontWeight','normal')
    43. %% 写出对比图像
    44. img_name = ['.\',dirname,'\HANTS-revised vs iHANTS-revised','(', LC_type, 'POI ',num2str(i),')'];% 按照POI的序号写出图像
    45. print(img_name,'-dpng','-r0'); % '-r0':以屏幕分辨率保存
    46. %% 关闭绘图
    47. % close all;