- 积分
- 1484
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-5-25
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 Wenfang 于 2017-1-6 12:36 编辑
程序:
clc;clear all;
currentFolder=pwd; %当前目录
work_rootdir=currentFolder(1:3);%工作文件夹所在磁盘 如j:\
addpath(genpath(work_rootdir));
%绘制网格和海岸线
clf; %清除所有绘图窗口
Lonlims=[70 140];
Latlims=[15 58];
m_proj('albers equal-area','lon',Lonlims,'lat',Latlims,'rect','on') ; %定义投影和矩形框(rect--on,否则默认为off,显示为扇形)
m_coast('color','k','linewidth',3); %加载海岸线
m_grid('linest','none','linewidth',2,'tickdir','out','xaxisloc','top','yaxisloc','right'); %加载经纬度网格
hold on;
%加载shapefile文件
%shapefile格式底图:国界、省界、市界、县界
%省界
shpdata_rootdir=[work_rootdir 'AtMyCodes\AtMyShpfiles\']; %shapefile文件所在文件夹
%国界
filepath_cnboundary=[shpdata_rootdir 'cn_boundary_wgs84\' 'cn_boundary_wgs84.shp'];
cnboundary=shaperead(filepath_cnboundary);
lonx_cnboud=[cnboundary(:).X]; %经度
laty_cnboud=[cnboundary(:).Y]; %纬度
m_plot(lonx_cnboud,laty_cnboud)
hold on;
%省界
filepath_provboundary=[shpdata_rootdir 'cn_prov_wgs84\' 'cn_prov_wgs84.shp'];
provboundary=shaperead(filepath_provboundary);
lonx_provboud=[provboundary(:).X]; %经度
laty_provboud=[provboundary(:).Y]; %纬度
m_plot(lonx_provboud,laty_provboud)
hold on;
(台风路径文件:
CH1949BST_hist.mat
(13.12 KB, 下载次数: 26, 售价: 5 贡献)
|
-
1949年CMA台风路径图
评分
-
查看全部评分
|