如何在Matlab中基于地理图绘制矩阵形式的地形变形图？

一、别用geoplot，换对工具 #

geoplot仅用于绘制地理点、线，无法直接展示矩阵形式的变形数据。推荐以下实用工具组合：

• geoshow + surf/pcolor：可将地理底图与变形表面/伪彩图叠加，适合区域变形展示
• meshgrid + surf：若变形数据对应规则网格，先将经纬度转为网格坐标，再生成三维表面图
• contourf + geoshow：绘制变形数据的填充等高线并叠加地图，适合查看二维分布

二、一步步实现（以规则网格变形数据为例） #

假设你拥有：

• 目标区域经纬度范围：lon_min、lon_max、lat_min、lat_max
• 变形数据矩阵deform_data（M行N列，每个元素对应一个网格点的变形量）

1. 生成匹配的经纬度网格
   将经纬度转换为与变形矩阵尺寸一致的网格，代码如下：

   [Lon, Lat] = meshgrid(linspace(lon_min, lon_max, size(deform_data,2)), ...
                        linspace(lat_min, lat_max, size(deform_data,1)));
   

   注意：矩阵行对应纬度（上下方向），列对应经度（左右方向），不要搞反。

2. 两种可视化方案任选
   方案1：三维变形表面+地理底图

   figure
   % 加载Matlab自带陆地底图
   geoshow('landareas.shp','FaceColor',[0.5 0.7 0.5]);
   hold on
   % 绘制变形表面，设置透明度以显示底图
   surf(Lon, Lat, deform_data, 'FaceAlpha', 0.6, 'EdgeColor', 'none');
   colormap jet; % 用渐变色标区分变形大小
   colorbar; % 显示颜色对应的变形数值
   hold off
   view(3); % 切换到三维视角，需二维视图可改为view(2)
   

   方案2：二维填充等高线+地理底图

   figure
   geoshow('landareas.shp','FaceColor',[0.8 0.8 0.8]);
   hold on
   contourf(Lon, Lat, deform_data, 20); % 20条等高线，数量可自行调整
   colormap cool;
   colorbar;
   hold off
   

3. 踩坑提醒

   • 尺寸匹配：Lon、Lat的尺寸必须与deform_data完全一致，可用size()命令检查
   • 数值范围：若变形数据数值过小可能无法显示，用max(deform_data(:))和min(deform_data(:))查看范围，必要时用caxis([min_val max_val])手动调整颜色条范围
   • 坐标系：默认支持WGS84（GPS常用坐标系），若数据为其他坐标系需先转换为经纬度

内容的提问来源于stack exchange，提问作者Gero