需求:在画好底图的图层上,根据显著性结果进行显著性打点
工具:matlab,函数:m_map
说明:使用打点函数时,需要注意单独构建一层 投影坐标,并且投影与原投影信息一致。
if strcmp(StipStr,'Yes')
% 建立新的坐标
axRi = axesm('miller','Frame','off','MapProjection',"miller","MapLatLimit",[-60 Lat(1,1)]);
hold on;
% 建立新的坐标
set(axRi,'Position',[fig.CurrentAxes.Position],'LineWidth',1.1,'Box','off','TickDir','in',...
'XMinorTick','off','YMinorTick','off','XGrid','off','YGrid','off','GridLineStyle','--',...
'FontName','Airl','FontSize',8,'FontWeight','bold','GridAlpha',.09) ;c
%底色调节为0
axRi.Color='none';
m_grid('linest','none','tickdir','in','fontsize',8,'fontname','Airl','FontWeight','bold','color','none','backgroundcolor','none');
% 选择投影地图与图片经纬度范围:与底图投影信息一致
M=m_proj('miller','lon',[Lon(1,1) Lon(1,end)],'lat',[-60 Lat(1,1)]);
hold on
% 这行的原理就是利用兴趣区mask和显著性点的数据进行掩膜
Region_PV=aimvalue .*LUCC_Mask.*pvalue;
Mask=Region_PV;
Mask(~isnan(Mask))=1;
%只显示感兴趣区的显著性点
scatterm(axRi,Lat.*Mask,Lon.*Mask,0.5,"filled","k");
end
注意:scatterm 函数使用时的投影