MATLAB工具包m_map最新版(2020)绘制中国地图

%% 参数
ColW = [0.6784    0.8510    0.9020];    % 海水颜色
ColL = [1 1 1];                         % 窗口颜色
ColB = [.6 .6 .6];                      % 中国底图边界
ColBk1 = [0.9294    0.6902    0.1294];  % 板块边界1颜色
ColBk2 = [0.0588    1.0000    1.0000];  % 板块边界2颜色
Lat = [15 53];                          % 纬度范围
Lon = [75 132];                         % 经度范围


m_proj('lambert', 'lat', Lat, 'long', Lon,'rectbox', 'on'); % 投影初始化
set(gcf, 'color', ColL);                                    % 图窗背景为白色
m_gshhs_i('patch',ColL,'edgecolor','none');                 % 高清海岸线
m_grid('box','fancy','backcolor',...
    ColW,'xtick',10,'ytick',6,...
    'linewidth',1,'Fontsize',12);                           % 图框参数

%% 中国边界
s = [pwd '\mapData\chinamap.txt']; 
d = load(s);
x = d(:, 1);
y = d(:, 2);
m_line(x, y, 'color', ColB, 'LineWidth', 1.2);

%% 版块边界
bk2 = load([pwd '\mapData\areaBankuai\2.txt']);
m_line(bk2(:,1),bk2(:,2),'color',ColBk2,'LineWidth',1.5);

bk1 = load([pwd '\mapData\areaBankuai\1.txt']);
m_line(bk1(:,1),bk1(:,2),'color',ColBk1,'LineWidth',1.5);

%% 画地震
str = 'D:\EqCatlog\CEIC_30d.EQT';
OpenEqt % 自定义函数
[~, T1] = DiscussFindTime(DataCell{7});
T0y = num2str(str2num(T1(1:4))-2);
T0 = [T0y '0101'];
EqT = DataCell{7};        % 时间
T_1 = datestr(datenum(T1,'yyyymmdd') ,'yyyymmdd');
T_2 = datestr(datenum(T1,'yyyymmdd') - 3 ,'yyyymmdd');
T_3 = datestr(datenum(T1,'yyyymmdd') - 7 ,'yyyymmdd');
color1 = [.2 .8 .1];
color2 = [.2 0 1];
color3 = [1 0 0];

xmin = Lon(1);
xmax = Lon(2);
ymin = Lat(1);
ymax = Lat(2);

MMplot(T0,T1,DataCell,color1,xmin,xmax,ymin,ymax,'1',0);   % 30天
MMplot(T_3,T_1,DataCell,color2,xmin,xmax,ymin,ymax,'1',0); % 15天
MMplot(T_2,T_1,DataCell,color3,xmin,xmax,ymin,ymax,'2',0); % 7天

%% ...以下代码省略