利用mmap工具箱显示shp、tif和散点图（m_map）

close all, clear all, clc
fpni='E:\MeteoInfo\map\china.shp';%中国地图shp文件所在路径
China=shaperead(fpni);
boux=[China(:).X];bouy=[China(:).Y];%分别是获取经度X信息和纬度Y信息
set(gcf,'position',[0 0 1440 780]);%设置图形窗口位置和大小
m_proj('Mercator','lon',[-170 170],'lat',[-80 80]);%设置投影方式为：墨卡托，地图显示范围
m_plot(boux,bouy,'k');%最关键的一句，绘制地图
%下面这句设置图形横纵坐标为经纬度格式
m_grid('linestyle','none','linewidth',2,'tickdir','out','xaxisloc','bottom','yaxisloc','left','fontsize',12);

%还有一种简单的方式显示shp
%m_mapshow('F:\研究生\深圳碳中和会议\experiment\data\anshan.shp')

hold on;

% load corrcoef.mat; %加载散点数据
% lon=coce(:,1);lat=coce(:,2);dataco=coce(:,3);
% m_scatter(lon,lat,80, dataco,'filled', 'MarkerFaceColor', 'flat', 'MarkerEdgeColor', 'w','linewi',1) ;%画实心点图
% m_scatter(lon,lat,50, dataco, 'MarkerFaceColor', w','linewi',2) ;%画空心点图
% contourcmap( 'jet', [-1:0.1:1], 'colorbar', 'on','location', 'vertical', 'fontsize', 13);%显示图例
% title('***', 'Rotation', 0, 'FontSize', 14);%显示图标题

[info map] = imread('F:\研究生\贝叶斯大气反演\CT\flux_pri_bio.tif');
Plat=90-.5-[0:179];Plon=-180+.5+[0:359];
[Plg,Plt]=meshgrid(Plon-0.5,Plat+0.5);
m_pcolor(Plg,Plt,info);
m_plot(boux,bouy,'k');%最关键的一句，绘制地图