M_map工具包函数研究

mytomorrow

- X& m/ [( w" ^- Z  U9 wm_map可以画出很漂亮的图

如果想在m_map的指定投影上画图必须使用m_text,m_quiver,m_pcolor等前面加了‘m_'的函数，那么这些函数与text,quiver,pcolor等普通matlab函数有什么区别呢？让我来分析一下m_plot的代码为大家解答吧。
. g! [, ]$ m5 ~0 r& {
m_plot代码如下：

function han = m_plot(varargin)
7 R9 o3 S" U# z7 w, i- R0 T
%将用户输入参数放在varargin中，这个varargin是个大小可变的cell数组，

% u  i' L" L: W%这样设计的好处是m_plot就和plot一样是一个输入参数可变的函数了
global MAP_PROJECTION MAP_VAR_LIST
$ u8 e3 J3 w6 S! U& [+ jif isempty(MAP_PROJECTION),  %判断有没有设定投影方式
  disp('No Map Projection initialized - call M_PROJ first!');
  return;
8 L* `9 S# Z2 O; U1 G3 uend;
if nargin < 2;  %判断输入参数的个数，小于2的话就好好查看一下m_plot函数的用法
  help m_plot
  return
2 ]6 l( n8 s# n2 _1 S, Qend
[x,y] = m_ll2xy(varargin{1},varargin{2});

6 R* e  C" ]. n8 K- u%m_ll2xy这个函数会根据当前投影的方式将varargin{1},varargin{2}也就是笛卡尔坐标系的横纵坐标进行转化
! C5 o6 c4 n8 A4 l) {varargin = varargin(：）; %将cell数组由一行变成一列
s = size(varargin,1); %s为cell数组的行数，也就是输入参数个数
/ ^, d  r. M& Wh=plot(x,y,varargin{3:s}); %调用plot函数，m_map里大部分函数都是如此，最终都要调用系统函数
if nargout == 1  %判断输出变量个数是否为1
  han = h;
0 M- r5 e; |1 d, {% ~# z' Jend

% }: _# W. X' J% _return
+ v3 {  E% w$ x. u5 @- I0 i3 ~# r0 Y3 R
- [" F  ~1 C  ^ 既然明白了m_map的基本原理，我们就可以自己为m_map添加函数了，比如m_map工具包没有m_scatter, 没关系，我们可以自己写啊，很简单的！

function han = m_scatter(varargin)
/ L# y: J, R5 [( M$ u/ V% written by Haidong Pan (Ocean University of China) 潘海东 （中国海洋大学）
% email: panhaidong_phd@qq.com
- q) }) m7 h$ H  E. J( u: s! p
global MAP_PROJECTION MAP_VAR_LIST

if isempty(MAP_PROJECTION),
6 H5 Z5 a8 K: ^( D2 C/ z  f0 Z  disp('No Map Projection initialized - call M_PROJ first!');
% _+ m8 ~0 L# k3 x  return;
end;

: X( b% n( F! X7 o6 v& H. l) I[x,y] = m_ll2xy(varargin{1},varargin{2});
. y' t. S- }, C$ T5 {varargin = varargin(：）;
s = size(varargin,1);
: |+ t  p" P9 P, S- `$ Y, Y  G: \h=scatter(x,y,varargin{3:s});
if nargout == 1
  han = h;
, l# f2 h2 i+ y$ H; ^8 gend
+ H0 X0 A: K& i3 S! ~
5 P2 G  _% x/ g& Q+ F9 ~return

NingW

M_map工具包函数研究