Argo(Array for Real-time Geostrophic Oceanography)計劃，經(jīng)過多個國家的共同努力，到2005年基本顯示全球覆蓋。Argo觀測深度為2000m，每10天提供一組觀測數(shù)據(jù)。Argo是一個國際項目，它使用剖面浮標來觀察地球海洋的溫度、鹽度、洋流，以及最近的生物光學(xué)特性;它從21世紀初開始運行。它提供的實時數(shù)據(jù)用于氣候和海洋學(xué)研究。

下圖展示了Argo浮標的全球分布圖：

本節(jié)將顯示怎樣讀取Argo數(shù)據(jù)并繪圖顯示，代碼如下：

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?

歡迎交流??！