《m_map在matlab中使用及投影说明》由会员分享，可在线阅读，更多相关《m_map在matlab中使用及投影说明（4页珍藏版）》请在金锄头文库上搜索。
1、m_map 工具箱使用笔记matlab 中 m_map 地图工具包的使用m_map 工具包下载地址及英文使用说明和例子：http:/www.eos.ubc.ca/rich/map.html考虑到那些英文水平比我还低的人，故作简单的介绍如下：m_map 工具包的加载：下载 m-map 工具箱后，要想 matlab 能够调用，首先必须加载。可以先在 matlab 中查看当前目录，pwd。将 m_map 工具箱移动到此目录（方便以后操作） 。然后输入 addpath m_map工具包 （根据具体文件夹而定） 。此时 m_map 工具箱加载完毕，可以调用使用。可以 help m_coast 检验。下面
2、是 m_map 工具箱中常用的一些投影方式的说明：选择投影方式时定义区域命令（有些特别的例外）：m_proj(投影方式,long,lonmin lonmax,lat,latmin latmax);m_coast(patch,a,b,c);画海岸线。主要是配色。a,b,c 都在 0 1间，a,b,c 都乘以 255 后为RGB 配色方案。当然，边线还可以加载自己想要的地图，在http:/www.maproom.psu.edu/dcw/上选择区域边界（有一步要选“download points”我不知道原因） ，下载*2pts.txt 文件后：m_plotbndry(states/*,color,
3、r); # *为匹配下载的文件名。m_grid 格网命令。参数有：linest 设置格网线条，默认为虚线。 xticklabels（yticklabels）设置坐标上的标记方式，很少用。xtick（ytick ）格网范围（有些投影时表示格的份数） 。color 线条颜色。yaxislocaion（xaxislocation）坐标位置，选项是 rightleft（topbottom）middle。tickdir（outin）边框属性。 box（onofffancy ） 。Linewidth 线条宽。整个坐标位置：axes（position,a,b,c,d)。a 为 x 的最小值，b 为 y 上的
4、最小值，c 为宽，d 为高。m_elev 画海拔线。如： m_elev(contourf,a,b,c).海拔等高线。a 为起始高，b 为间距，c 为最高点。matlab 中导入数据命令：load data.*我们常用到的数据 data（经度纬度 值） ，从中提取某列如 x=data（:,1）提取第一列给 x。m_ungrid *能撤销之前的任何*操作。写到这里，完全没思路了，你们可以将说明书下下来。根据上面讲的应该看得懂了。下面呢将例子中比较难懂的语句解释：例一中：第一句中两个值表示视角的观测位置。patch(.55*-1 1 1 -1,.25*-1 -1 1 1-.55,b);第一个逗号前的
5、部分表示从左下角逆时针转的 x值，第二个逗号前中前表示从左下角逆时针点的 y 值，-.55 表示相对观测点 y 轴方向的偏移。w 为所画方块的填充颜色。例二中：m_elev(contourf,500:500:6000);从海拔 500 到 6000 以 500 为间距画等高线图。最后一句是配色方案。例三中：m_proj(stereographic,lat,90,long,30,radius,25);前两个数为视角位置。25 为纬度所选区域范围。例如画整个南半球：m_proj(stereographic,lat,-90,radius,90,rotangle,180); -90 为观测点纬度，90
6、 为纬度范围，180 为绕南北轴顺时针旋转的度数。m_grid(xtick,12,tickdir,out,ytick,70 80,linest,-)；12 为经度方向份数，这里强调下，份数不是严格的，份数的选择是 量子化的，不然不是你所期望的图。例四中的 set 一句是设置地图大小的。例 9 中选择投影方式时，最后的参数rect,on表示地图的正交性。为 on 时是块形图，为 off 时类似 lambert 所做的图。例 10 中dates=datenum(1997,10,23,15,1:41,zeros(1,41);设置时间属性。1997 年 10 月 23 日，15 时 1到 40 分，秒
7、全为 0.m_track(lons,lats,dates,ticks,0,times,4,dates,8,.clip,off,color,r,orient,upright)；画轨迹，时间间距为 4；数据等分为 8.最后两个参数是上面标记的属性，如方向、朝向。例 11 中：第一句中 clong 170表示经度的中心位置是东经 170 位置。m_line(100.5,13.5,marker,square,color,r);表示在（100.5，13.5）处加入方形标记，颜色为红色。m_range_ring(100.5,13.5,1000:1000:15000,color,b,linewi,2); 画
8、一系列间距线圈。离中心1000 开始，1000 为间距，15000 结束。颜色为蓝色，线宽 2.例 12 中：前两行是确定图中多边形的具体位置（各个顶点的坐标） 。m_hatch(bndry_lon,bndry_lat,single,30,5,color,k); % .with hatching added.填充多边形。填充线条倾角为 30 度，线条间距为 5 个单位。卫星数据还没具体看，自己动手了，只能。上面笔记忘对某些人有点用处，参考，仅供。Matlab M_map 工具箱 19 类投影(Project) 创建语句及图像生成一共有 19 种投影类型，分别作图，代码和图形在后面基本绘图三步走
9、1、m_proj(proj name,) 创建投影2、m_coast; 绘出海岸，可以用 patch 命令将陆地上色，也可用更精细的 m_gshhs_系列命令描绘海岸3、m_grid; 根据投影要求绘制网格不同投影的后面的设置内容略有不同，具体可以利用 m_proj(set,projection name)查阅可设置参数，或用 m_proj get 查询当前参数（必须已完成第一步创建投影）以下为代码，为保持格式，句首都添加了注释行。复制到 matlab m 文件里，全选并 ctrl+t即恢复。PS，要运行该代码，需添加 gshhs 海岸数据。否则，请将所有 m_gshhs 系列命令，改为m_c
10、oast 命令。% clear% clc% project=.% % 1-5 为方位图，圆形，lon,center long,lat,center lat,.% % rad,degree|long,lat % boundary% Stereographic,. % 1 方位图% Orthographic,. % 2 方位图% Azimuthal Equal-area,. % 3 方位图% Azimuthal Equidistant,. % 4 方位图% Gnomonic,. % 5 方位图% .% 6 透视图，方形，lon,center long,lat,center lat,.% .% al
11、t,altitude_fraction % 透视高度 % .a satellite in an orbit of radius 3 earth radii would have an altitude of 2% Satellite,. % 6 方位图% .% 7-8 锥顶投影图 适合中纬度地区，且东西延伸范围较大% Albers Equal-Area Conic,. % 7% Lambert Conformal Conic,. % 8% .% 9-11 lon (min max|center),lat(maxlat|min max)% .% lon 可设置中间位置，lat 若设置最大值，则表
12、示对称的 lat 范围，或给定值 % Mercator,. % 9% Miller Cylindrical,. %10 适合全球地图% Equidistant Cylindrical,. %11% .% 12 lon,G1 G1,lat,L1 L2,dir,(horizontal|vertical)% Oblique Mercator,. %12 适合细长图像， G1L1 G2L2 为两头中间坐标% Transverse Mercator,. %13 % Sinusoidal,. %14 两侧外凸，视觉鼓起，像地球仪的一片% Gall-Peters,. %15 经纬均直拉伸至正方，适合赤道区域
13、，% Hammer-Aitoff,. %16 全球% Mollweide,. %17 全球% Robinson,. %18 全球% UTM; %19% n=length(project) ;% n1=cell(1,n);% region=16 45 100 150;% region=south,north,east,west% for kk=1:19% clf% % =set projection=% if kk% m_proj(projectkk,lon,mean(region(3:4),lat,mean(region(1:2),rad,15)% elseif kk=6 % 6 alt 卫星轨道 圆的% m_proj(projectkk,lon,mean(region(3:4),lat,mean(region(1:2),.% rad,4.5,alt,2)% elseif kk=7|kk=8; % 7-8 Conic Projection 锥顶投影 方的% m_proj(projectkk,lon,region(3),region(4),lat,region(1),region(2);% elseif (kk=9 & kk=16 & kk=16 & kk =18% m_gshhs_c(patch,.8 .8 .8);% m_grid(box,on,xaxisloc
