绘制IGS站点图

#绘制IGS站点图

下载m_map

进入官网，下载m_map存储包
https://www.eoas.ubc.ca/~rich/map.html

如图所示：

matlab添加路径:

点击“添加路径”
：

绘图

1.在不适用m_map包的情况下：转载于：https://blog.csdn.net/sinat_39238867/article/details/100177423

load('BL.txt');

h = geoshow('landareas.shp', 'FaceColor', 'c');
grid on
hold on

xlabel('Longitude');
ylabel('Latitude');%坐标轴标题
set(gca,'Ylim',[-90,90],'ytick',[-90:30:90]);
set(gca,'yticklabel',{'90°S','60°','30°','0°','30°','60°','90°N'});
set(gca,'Xlim',[-180,180],'xtick',[-180:30:180]);
set(gca,'xticklabel',{'180°W','150°','120°','90°','60°','30°','0°','30°','60°','90°','120°','150°','180°E'});
%坐标轴范围及刻度分划，坐标轴文字替代
set(gca,'Box','on');%坐标轴是否为四面
set(gca,'FontSize',10,'Fontname', 'Times New Roman','Fontweight', 'bold');%字号、字体、是否加粗
set(gca,'GridAlpha',1,'GridLineStyle','--');%格网透明度（0-1）及线型
title('MGEX Station Location','FontSize',14,'Fontweight', 'bold');%图标题

scatter(BL(:,2),BL(:,1),'MarkerEdgeColor','k','MarkerFaceColor','r');

print -f1 -r1200 -dpng MGEX;

result::

2.使用m_map转载于：https://blog.csdn.net/SmartTiger_CSL/article/details/103128917

%close all,clc, dbstop if error
load('station.txt');
rLon=station(:,1);
rLat=station(:,2);
%rsta=num2str(station(:,3));

r=size(rLon,1);
s=size(rLon,1);
 
m_proj('miller','lat',[-90 90],'long',[-180 180]);  %确定投影方式和绘图界线  矩形外框 
%m_proj('robinson','lat',[-90 90],'long',[-180 180]);  %确定投影方式和绘图界线  椭圆外框
m_coast('linewidth',1,'color','r');%画出海岸线,设置线宽和颜色
m_coast('patch',[0.7 1.0 0.7],'edgecolor','none'); %海岸线填充
m_grid('box','fancy','linestyle','-','gridcolor','w','backcolor',[0.2 0.65 1]);%添加格网
     
    for i=1:r
        if i<=s
            m_line(rLon(i),rLat(i),'Marker','o','LineWidth',0.75,'MarkerSize',8,'color','black');%绘制测站点，及其形状、大小和颜色
            hold on;
            m_line(rLon(i),rLat(i),'Marker','^','LineWidth',0.75,'color','r');
        end
        if i>s
            m_line(rLon(i),rLat(i),'Marker','v','LineWidth',20.755,'color','r');
        end
        %m_text(rLon(i),rLat(i),rsta(i,:),'vertical','top');%将测站号标在图中   
       % m_text(rLon(i),rLat(i),'vertical','top');%将测站号标在图中 
    end
      
    %set(gcf,'units','inches','position',[12 5 5 4]);%设置图像位置和大小 设置图像在scree中的位置[left bottom width height]

     title('GPS station map','fontsize',14,'fontweight','bold');
     
     set(gcf,'color','w');   % Need to do this otherwise 'print' turns the lakes black

在参考上述了上述两篇文章后：

%close all,clc, dbstop if error
m=load('station.txt');
rLon=m(:,1);
rLat=m(:,2);
%rsta=num2str(station(:,3));
 hold on;

r=size(rLon,1);
s=size(rLon,1);
disp(r);
disp(s);
%scatter(station(:,2),station(:,1),'MarkerEdgeColor','k','MarkerFaceColor','r');%绘制 
m_proj('miller','lat',[-90 90],'long',[-180 180]);  %确定投影方式和绘图界线  矩形外框 
%m_proj('robinson','lat',[-90 90],'long',[-180 180]);  %确定投影方式和绘图界线  椭圆外框
m_coast('linewidth',1,'color','r');%画出海岸线,设置线宽和颜色
m_coast('patch',[0.7 1.0 0.7],'edgecolor','none'); %海岸线填充
m_grid('box','fancy','linestyle','-','gridcolor','w','backcolor',[0.2 0.65 1]);%添加格网

 for i=1:r
   m_line(rLat(i),rLon(i),'Marker','o','MarkerFaceColor','g','LineWidth',0.75,'MarkerSize',5,'color','black');%绘制测站点，及其形状、大小和颜色       
 end     
title('GPS station map','fontsize',14,'fontweight','bold');
    
set(gcf,'color','w');   % Need to do this otherwise 'print' turns the lakes black