GESLA是英國海平面專家Woodworth等人收集了全球各個國家的驗潮站觀測構建的水位數(shù)據(jù)集，下載地址是?https://gesla.org/ 。如下圖所示， 從網(wǎng)上下載下來的數(shù)據(jù)都是文件格式，這種數(shù)據(jù)格式與我們常用的nc格式不太一樣。

下面我來介紹一下批量讀取GESLA數(shù)據(jù)的技巧

file_path='G:\數(shù)據(jù)\GESLA\';

img_path_list=dir(); ? %這個命令能獲取當前文件夾所有的文件名

k=0

for j=3:1278 ? ? ? %第一個和第二個似乎是隱藏文件

fid=fopen(img_path_list(j).name,'r');

?for i=1:4

? ?str=fgetl(fid); %讀取一行字符串

?end

?lat=fgetl(fid);lat=strtrim(lat(11:end));

?lon=fgetl(fid);lon=strtrim(lon(12:end));

?fclose(fid);k=k+1;

?lat1(k)=str2num(lat);lon1(k)=str2num(lon);

end

lon1(lon1>180)=lon1(lon1>180)-360;

m_proj('miller','lat',[-85 85]); ??

m_coast('patch',[.7 1 .7],'edgecolor','none');?

m_grid('box','fancy','linestyle','-','gridcolor','w','backcolor',[.2 .65 1]);

hold on;?m_plot(lon1,lat1,'r.');?set(gcf,'color','w');?