M_map绘制台风路径的一个实例_新手入门笔记

Wenfang

程序：

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 贡献)

2017-1-6 12:29 上传

点击文件名下载附件
49年台风路径_strucutre array ）


year='1949';
test_tc_track_mat=['CH' year 'BST_hist.mat'];
load([work_rootdir 'AtMyCodes\CMABSTdata\CMABSTdata_hist_mat\' test_tc_track_mat]);

markersize=5; linewidth=3;

v_categories = [17.2 24.5 32.7 41.5 51.0 1000.0]; %(unit m/s)
colors_      = [ 99 184 255 ; %light blue  28 134 238 ; %blue  热带低压
%                   0 139   0 ; %green  
                255 215   0 ; %yellow   TS  热带风暴
                238 154   0 ; %dark yellow STS  强热带风暴
                238 118   0 ; %orange   TY   台风
                238  64   0 ; %red      STY  强台风
                139  37   0 ; %dark red  SuperTY 超强台风
                ]/255;


if isstruct(tc_track)
    track_count = length(tc_track);

    for track_i = 1:track_count %every track
        for node_i = 1:length(tc_track(track_i).MaxSustainedWind)-1 %every node for specific track
            if tc_track(track_i).MaxSustainedWind(node_i)>0

                v       = tc_track(track_i).MaxSustainedWind(node_i);
                v_color = find (v < v_categories);  %返回符合条件的所有数组下标
                if numel(v_color)>0
                    v_color = v_color(1); %强度类别归入第一个下标所属的类别，颜色为第一个
                    h(node_i) = m_plot(tc_track(track_i).lon(node_i:node_i+1), tc_track(track_i).lat(node_i:node_i+1),...
                            'color',colors_(v_color,:),'markersize',markersize,'linewidth',linewidth);
                end
            end
        end

    end
end

%绘制图例
    for cat_i = 1:length(v_categories)  %6类

        g(cat_i) = m_plot(200,40,'o','color',colors_(cat_i,:),'markerfacecolor',colors_(cat_i,:),'markersize',markersize);

    end

   legendstr = {'TD','TS','STS','TY','STY','SuperTY'};
   L = legend(g,legendstr,'location','southwest');

    legend('boxoff')
    set(L,'FontSize',8);

tianqh

没看明白，学习中

二爷名声在外

同时希望楼主把m文件当作附件发上来

aas314

好图 好看  学习 了

dx01024

数据要处理成什么样的啊？

2daidai

%加载shapefile文件
%shapefile格式底图：国界、省界、市界、县界
这个在哪里下？

妙龄童

非常感谢楼主的分享，仍在学习中

青之梅子

谢谢楼主分享

bingxinkeke

这个厉害了啊，我刚刚准备要画这样显示台风强度的路径图！

Yyyoungwu

xuexiyixia