clc
figure('unit','centimeters','position',[10 1 17 17])
lon1=110;
lon2=145;
lat1=5;
lat2=25;
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\bwp212008\'; %指定批量数据所在的文件夹
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))
s = string(filename.textdata(:,9));
e=tfs(s).*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
数据来源:https://www.metoc.navy.mil/jtwc/jtwc.html?western-pacific
运行代码也不报错,但数据也出不来,不知道为什么,求大神解答
设几个断点调试一下,观察matlab工作区的变量值更新情况
我们的数据来源一样,可惜我不会,也在找解决办法