M_map工具包函数研究

mytomorrow

m_map可以画出很漂亮的图
. @$ t6 E* G/ v3 Z0 _; p8 j
如果想在m_map的指定投影上画图必须使用m_text,m_quiver,m_pcolor等前面加了‘m_'的函数，那么这些函数与text,quiver,pcolor等普通matlab函数有什么区别呢？让我来分析一下m_plot的代码为大家解答吧。
6 K8 L: f; g/ k8 r! `
m_plot代码如下：
% S8 }- R; U5 P  O
" t+ M& l- \& @" Pfunction han = m_plot(varargin)

0 b& p& P0 A+ s7 \%将用户输入参数放在varargin中，这个varargin是个大小可变的cell数组，

4 _; x9 d+ X  T& O5 ?. j%这样设计的好处是m_plot就和plot一样是一个输入参数可变的函数了
2 }( q! [/ }/ f6 E5 M: @: T0 F! hglobal MAP_PROJECTION MAP_VAR_LIST
if isempty(MAP_PROJECTION),  %判断有没有设定投影方式
  disp('No Map Projection initialized - call M_PROJ first!');
  return;
. }8 L+ q% m( V% n# mend;
if nargin < 2;  %判断输入参数的个数，小于2的话就好好查看一下m_plot函数的用法
3 C# e! c6 s' Y2 n; i' [  help m_plot
! j  s- p0 C! _; M9 ~& M. |  return
, {8 F! U4 Y/ vend
[x,y] = m_ll2xy(varargin{1},varargin{2});

%m_ll2xy这个函数会根据当前投影的方式将varargin{1},varargin{2}也就是笛卡尔坐标系的横纵坐标进行转化
+ C/ p4 f2 s2 I2 }5 Jvarargin = varargin(：）; %将cell数组由一行变成一列
+ M5 X8 h" T9 j8 U! _s = size(varargin,1); %s为cell数组的行数，也就是输入参数个数
h=plot(x,y,varargin{3:s}); %调用plot函数，m_map里大部分函数都是如此，最终都要调用系统函数
if nargout == 1  %判断输出变量个数是否为1
4 S# G4 B5 U3 I5 k8 Z' x  han = h;
" F% D2 L: `* q" Pend

return

' A% x* f1 }9 S2 k: [7 I% | 既然明白了m_map的基本原理，我们就可以自己为m_map添加函数了，比如m_map工具包没有m_scatter, 没关系，我们可以自己写啊，很简单的！

function han = m_scatter(varargin)
6 `2 T8 w1 T  n+ ?9 R, b' ~% written by Haidong Pan (Ocean University of China) 潘海东 （中国海洋大学）
7 K) h! ~( N& @' i) S% email: panhaidong_phd@qq.com
) ~& b  A& G# o# N% _% d1 k7 W% n4 L$ {
global MAP_PROJECTION MAP_VAR_LIST

if isempty(MAP_PROJECTION),
  disp('No Map Projection initialized - call M_PROJ first!');
% H) j- t& c6 D  return;
end;
1 I/ i" P( C6 b
[x,y] = m_ll2xy(varargin{1},varargin{2});
6 p# c/ C( k2 z+ L& u$ Pvarargin = varargin(：）;
+ b  X2 R& Y, V3 S( ss = size(varargin,1);
h=scatter(x,y,varargin{3:s});
9 Q1 k+ y9 e4 b- n% o# H7 Y5 Zif nargout == 1
/ K" W: L1 W% V( O  han = h;
1 Y+ M% c. R' W& l: Lend

return

NingW

M_map工具包函数研究