Argo数据的读取与显示

Argo(Array for Real-time Geostrophic Oceanography)计划，经过多个国家的共同努力，到2005年基本显示全球覆盖。Argo观测深度为2000m，每10天提供一组观测数据。Argo是一个国际项目，它使用剖面浮标来观察地球海洋的温度、盐度、洋流，以及最近的生物光学特性;它从21世纪初开始运行。它提供的实时数据用于气候和海洋学研究。

Argo浮标的观测流程示意图

下图展示了Argo浮标的全球分布图：

图片来源：https://argo.ucsd.edu/

本节将显示怎样读取Argo数据并绘图显示，代码如下：

address = 'ArgoData2005.nc';

lon = ncread(address,'LONGITUDE');

lat = ncread(address,'LATITUDE');

temp = ncread(address,'TEMP');

salt = ncread(address,'SALT');

%-------------grid--------------

[lon, lat] = meshgrid(lon,lat);

LLZ.lon = lon';

LLZ.lat = lat';

LLZ.rg = temp(:,:,1);

% rg_plot(LLZ)

subplot(121)

m_proj('miller','lon',[0,360],'lat',[-90,90]); 

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

m_pcolor(lon,lat,temp(:,:,2)');

m_grid('tickdir','out','linewi',2); 

As = shaperead('operational_floatst.shp')

subplot(122)

m_proj('miller','lon',[-180,180],'lat',[-90,90]); 

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

% m_pcolor(lon,lat,temp(:,:,2)');

m_grid('tickdir','out','linewi',2); 

A = [As.X;As.Y]';

for i=1:length(A)

    [X,Y]=m_ll2xy(A(i,1),A(i,2));

    line(X,Y,'marker','.','markersize',10,'color','b');

end 

运行结果

欢迎交流！！【内容同步见本人的B站账号：我是水怪的哥】