利用matlab M_Map工具制作GPS测站分布图

原创最新推荐文章于 2026-09-22 12:17:22 发布 5.7k 阅读 AI写同款

收录于

MATLAB

SmartTiger_CSL

码龄8年墨衍会员 · AI 创作全网分发

217粉丝28篇原创

关注

matlab绘制站点分布图 author：wangjian time：2019年8月26号 function：读取地形数据，和是测站点数据的经纬度坐标并做图说明：程序参照了网上的m_map上的一些程序，本程序是以珠江河口为例写的 close all;clear;clc addpath m_map filename=‘E:\其他数据\水深\ETOPO1_Ice_g_gmt4.grd’; ncdisp(filename); %... 阅读详情

AI摘要

本文详细介绍了如何在MATLAB环境下安装M_Map工具箱，包括下载、解压及路径设置，并通过测试代码验证安装成功。同时，展示了如何利用M_Map绘制GPS测站分布图，包括数据读取、地图投影、测站点绘制及标注。

一、安装M_Map

首先去官网(https://www.eoas.ubc.ca/~rich/map.html)下载M_Map文件，进入M_Map下载页面，如下图：

M_Map官网下载界面
2. 点击gzipped tar-file下载m_map1.4.tar.gz压缩文件（大约8.28MB）。
3. 解压下载得到的m_map1.4.tar.gz压缩文件，将M_Map文件复制到MATLAB安装目录下的toolbox中，例如复制到D:\Program Files\MATLAB\R2014a\toolbox中。然后在MATLAB软件主页上找到“设置路径”按钮，点击“设置路径”按钮,出现“设置路径”对话框(如下图)，然后点击“添加文件夹”，找到某盘:\Program Files\MATLAB\R2014a\toolbox中m_map文件，然后选中m_map，点击“选择文件夹”。
在这里插入图片描述
M_Map文件路径被添加之后，点击“保存”，然后关闭即可。
4. 测试M_Map工具箱是否安装成功。在命令行窗口输入help m_map，如果出现M_Map版本、功能介绍等信息，就说明安装成功了。若显示m_map not found，说明没有安装成功，需要检查并重新配置。
在这里插入图片描述

二、测试m_map是否可用。

 打开MATLAB，新建一个脚本，然后输入如下代码：

close all, clear all, clc, dbstop if error

m_proj('ortho','lat',48','long',-123');
m_coast('patch','r');
m_grid('linest','-','xticklabels',[],'yticklabels',[]);

text(0,0,'M\_Map is installed successfully!','fontsize',25,'color','b',...
    'verticalalignment','middle','horizontalalignment','center');

保存、运行、得到下面图形：
在这里插入图片描述

三、绘制GPS测站分布图。

首先准备好GPS测站数据“station.txt“（经度、纬度和测站号）。按照以下代码，修改相应选项，将测站绘制在地图中：
station.txt中的数据格式：
119.75 49.22 50527
125.22 49.17 50557

close all, clear 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',[0 60],'long',[60 150]);  %确定投影方式和绘图界线  矩形外框 
%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');%将测站号标在图中     
    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