我的数据来源: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
不知道为什么画出来的只有底图,没有路径,求大神解答!!
您好,请问您的问题解决了吗?
经纬度的单位和脚本里的不一致