GESLA数据读取技巧

GESLA是英国海平面专家Woodworth等人收集了全球各个国家的验潮站观测构建的水位数据集，下载地址是 https://gesla.org/ 。如下图所示， 从网上下载下来的数据都是文件格式，这种数据格式与我们常用的nc格式不太一样。

下面我来介绍一下批量读取GESLA数据的技巧

file_path='G:\数据\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');