用matlab画台风轨迹图?

我的数据来源:https://www.metoc.navy.mil/jtwc/jtwc.html?western-pacific

代码如下:

clc
figure('unit','centimeters','position',[10 1 17 17])
lon1=-80;
lon2=85;
lat1=-80;
lat2=70;
load coast
m_proj('mercator','lon',[lon1 lon2],'lat',[lat1 lat2]);
axesm('MapProjection','mercator','MapLatLimit',[lat1 lat2],'MapLonLimit',[lon1 lon2])
framem
m_coast('patch',[0.7 0.7 0.6]);
m_grid('box','fancy','tickdir','in');
lat(end)=NaN;

datadir='D:\TYPHON\1\bwp282009.txt\'; %指定批量数据所在的文件夹
filelist=dir([datadir,'*.txt']); %指定批量数据的类型
k=length(filelist);
for s=1:k
    filename=importdata([datadir,filelist(s).name]);
  [m,n]=size(filename.textdata);
    a = string(filename.textdata(:,7));
    b = string(filename.textdata(:,8));
    tdata = string(filename.textdata(:,3)); 
    if n>=9
       if ~cellfun(@isempty,filename.textdata(:,9))
       g = string(filename.textdata(:,9)); 
       e=tfs(g).*1.852/3.6;
       end
    else
      e = filename.data.*1.852/3.6;
    end   
    tlat=tflatlon(a);
    tlon=tflatlon(b);
    ttdata=string(tfdate(tdata));
    t1data=char(tdata(1));
    textm(tlat(1)-0.7,tlon(1)-1,tdata(1),'fontsize',8,'SelectionHighlight','on')  
    geoshow(tlat,tlon, 'Marker','.','MarkerEdgeColor','k')  
    geoshow(tlat(1),tlon(1), 'Marker','pentagram','MarkerEdgeColor','k','MarkerSize',15)  
for tt = 8:4:length(a)
    geoshow(tlat(tt),tlon(tt), 'Marker','.','MarkerEdgeColor','r','MarkerSize',e(tt)*2)
    textm(tlat(tt)+0.2,tlon(tt)+0.2,tdata(tt),'fontsize',7,'SelectionHighlight','on')
end
end

 

不知道为什么画出来的只有底图,没有路径,求大神解答!!

 

您好,请问您的问题解决了吗?

 

经纬度的单位和脚本里的不一致