15 - matlab m_map地学绘图工具基础函数 - 一些数据转换函数（二）

0. 引言
1. 关于m_geodesic
2. 关于mygrid_sand2
3. 结语

0. 引言

通过前面篇节已经将m_map绘图工具中大多绘图有关的函数进行过介绍，已经能够满足基本的绘图需求，本节介绍下m_map中关于提取数据的两个函数，在特定应用场景下也是极为重要的。首先介绍下面几种数据转换函数：。

1. 关于m_geodesic

m_geodesic 函数用于计算大圆航线或大圆弧的坐标点。大圆航线是地球表面上两点之间的最短路径，而大圆弧是大圆航线在地球表面上的投影。

m_geodesic函数的基本用法：

[lat, lon] = m_geodesic(lat1, lon1, lat2, lon2, npts)

其中：

lat1 和 lon1 是起始点的纬度和经度。
- lat2 和 lon2 是终点的纬度和经度。
npts 是用于近似大圆弧的点的数量。

这个函数返回一个列向量 lat 和一个列向量 lon，它们包含了大圆航线或大圆弧上的坐标点的纬度和经度。

例如，假设你想要在地图上绘制从纽约（40.7128° N，74.0060° W）到洛杉矶（34.0522° N，118.2437° W）的大圆弧，你可以使用 m_geodesic 函数来计算坐标点，然后使用 plot 函数将其绘制出来。

clc;clear;
lat1 = 40.7128; lon1 = -74.0060; % 纽约
lat2 = 34.0522; lon2 = -118.2437; % 洛杉矶
npts = 100; % 用于近似的点的数量

% 计算大圆弧上的坐标点
[lat, lon] = m_geodesic(lon1,lat1, lon2,lat2, npts);

% 绘制大圆弧
plot(lon, lat, 'r', 'LineWidth', 2);

这段代码计算了从纽约到洛杉矶的大圆弧上的坐标点，并使用 plot 函数将其绘制为一条红色的曲线。

2. 关于mygrid_sand2

mygrid_sand2函数用于提取经纬度范围内的水深数据。

[image_data,vlat,vlon] = mygrid_sand2(region,ssfname)
%  MYGRID_SAND2  Read bathymetry data from Sandwell Database
%    [Z,LAT,LON] = MYGRID_SAND2(REGION) extracts data from
%    the Sandwell and Smith bathymetry, which is now at 1-minute
%    resolution.
%
%
% WARNING: change ssfname and ssversion to the correct one for 
%          your machine
%
%                                               Catherine de Groot-Hedlin
%                                               modified Rich Pawlowicz
%
% latitudes must be between -80.738 and 80.738;
%       input:
%               REGION =[west east south north];
%       output:
%               Z - matrix of sandwell bathymetry/topography
%               LAT - vector of latitudes associated with image_data
%               LON - vector of longitudes

其中，region为研究区四至；ssfname数据源路径，需要注意数据的是否匹配。

mygrid_sand2函数使用示例：

Ulo1 = [-90,-68];                          % 经纬度范围大致是美国东海岸的地形水深范围
Ula1 = [30,42];
m_proj('mercator','long',Ulo1,'lat',Ula1); % 绘图前需指定投影方式与范围
 
[Sele,Sla,Slo] = mygrid_sand2([-90 -68 30 42]);
figure()
%m_contour(Slo-360,Sla,Sele); % 西半球270转化为-90度，以此类推
m_contourf(Slo-360,Sla,Sele); % 西半球270转化为-90度，以此类推{contourf效果更好}
m_grid('box','fancy','tickdir','in','gridlines','no','fontsize',12)
 
colorbar