AVISO衛(wèi)星測(cè)高格網(wǎng)數(shù)據(jù)下載與讀取
法國(guó) AVISO(Archiving, Validation and Interpretation of Satellite Oceanographic)提供的格網(wǎng)化數(shù)據(jù)(MSLA).該數(shù)據(jù)融合了 T/P, Jason-1, ERS 和 ENVISAT 等多顆衛(wèi)星 的測(cè)高資料. 數(shù)據(jù)時(shí)間分辨率為 7 天, 空間分辨率為 1°×1°, 時(shí)間跨度為 1993 年 1 月至 2021年7月. 數(shù)據(jù)經(jīng)過了必要的標(biāo)準(zhǔn)改正, 如電離層延遲改正、對(duì)流層干濕分量改正、固體潮和海潮改正、海洋負(fù)荷潮汐 改正、極潮改正、電磁偏差改正、儀器改正和反變氣 壓計(jì)改正等。下面是鏈接的網(wǎng)址：
https://www.aviso.altimetry.fr/en/data/products/sea-surface-height-products/global/gridded-sea-level-anomalies-mean-and-climatology.html#c10358.
進(jìn)入官網(wǎng)后，下面的界面就是每月平均的格網(wǎng)數(shù)據(jù)（數(shù)據(jù)下載需要提前注冊(cè)賬號(hào)）
選擇第二個(gè)文件夾，即可進(jìn)入下載界面：
The set of variables are the following :
EKE (Eddy Kinetic Energy)
MSLA-H (Sea Level Anomaly H,?sea surface height above sea level)
MSLA-U? (Sea Level Anomaly U,?surface geostrophic eastward sea water velocity assuming sea level for geoid)
MSLA-V (Sea Level Anomaly V,?surface geostrophic northward sea water velocity assuming sea level for geoid)
MADT-H (Absolute Dynamics Topography H : sea surface height above geoid)
MADT-U (Absolute Dynamics Topography U, surface geostrophic eastward sea water velocity)
MADT-V (Absolute Dynamics Topography V, surface geostrophic northward sea water velocity)
這里我們下載了MSLA-H和MADT-H進(jìn)行對(duì)比：
文件內(nèi)部的信息：
Global Attributes:
? ? ? ? ? ?cdm_data_type? ? ? ? ? ? ? ? ? = 'Grid'
? ? ? ? ? ?date_issued? ? ? ? ? ? ? ? ? ? = '2022-02-25 16:57:35Z'
? ? ? ? ? ?time_coverage_resolution? ? ? ?= 'P1M'
? ? ? ? ? ?creator_email? ? ? ? ? ? ? ? ? = 'aviso@altimetry.fr'
? ? ? ? ? ?product_version? ? ? ? ? ? ? ? = '7.0'
? ? ? ? ? ?references? ? ? ? ? ? ? ? ? ? ?= 'www.aviso.altimetry.fr'
? ? ? ? ? ?Metadata_Conventions? ? ? ? ? ?= 'Unidata Dataset Discovery v1.0'
? ? ? ? ? ?creator_url? ? ? ? ? ? ? ? ? ? = 'http://www.aviso.altimetry.fr'
? ? ? ? ? ?keywords? ? ? ? ? ? ? ? ? ? ? ?= 'Oceans > Ocean Topography > Sea Surface Height'
? ? ? ? ? ?keywords_vocabulary? ? ? ? ? ? = 'NetCDF COARDS Climate and Forecast Standard Names'
? ? ? ? ? ?institution? ? ? ? ? ? ? ? ? ? = 'CNES, CLS'
? ? ? ? ? ?license? ? ? ? ? ? ? ? ? ? ? ? = 'http://www.aviso.altimetry.fr/fileadmin/documents/data/License_Aviso.pdf'
? ? ? ? ? ?geospatial_vertical_resolution = 'point'
? ? ? ? ? ?creator_name? ? ? ? ? ? ? ? ? ?= 'SSALTO/DUACS'
? ? ? ? ? ?ssalto_duacs_comment? ? ? ? ? ?= 'The reference mission used for the altimeter inter-calibration processing is Topex/Poseidon between 1993-01-01 and 2002-04-23, Jason-1 between 2002-04-24 and 2008-10-18, OSTM/Jason-2 between 2008-10-19 and 2016-06-25, Jason-3 since 2016-06-25.'
? ? ? ? ? ?standard_name_vocabulary? ? ? ?= 'NetCDF Climate and Forecast (CF) Metadata Convention Standard Name Table v37'
? ? ? ? ? ?date_modified? ? ? ? ? ? ? ? ? = '2022-02-25 16:57:35Z'
? ? ? ? ? ?summary? ? ? ? ? ? ? ? ? ? ? ? = 'Delayed Time Level-4 monthly mean of sea surface height above geoid products from multi-satellite observations over Global Ocean.'
? ? ? ? ? ?project? ? ? ? ? ? ? ? ? ? ? ? = 'SSALTO/DUACS: Data Unification and Altimeter Combination System'
? ? ? ? ? ?source? ? ? ? ? ? ? ? ? ? ? ? ?= 'Altimetry measurements'
? ? ? ? ? ?contact? ? ? ? ? ? ? ? ? ? ? ? = 'aviso@altimetry.fr'
? ? ? ? ? ?geospatial_vertical_positive? ?= 'down'
? ? ? ? ? ?title? ? ? ? ? ? ? ? ? ? ? ? ? = 'DT merged all satellites Global Ocean Ocean Gridded Monthly Mean of Absolute Dynamic Topography L4 product'
? ? ? ? ? ?geospatial_vertical_units? ? ? = 'm'
? ? ? ? ? ?processing_level? ? ? ? ? ? ? ?= 'L4'
? ? ? ? ? ?history? ? ? ? ? ? ? ? ? ? ? ? = '2022-02-25 16:57:35Z: Created by DUACS DT V7.0'
? ? ? ? ? ?date_created? ? ? ? ? ? ? ? ? ?= '2022-02-25 16:57:35'
? ? ? ? ? ?Conventions? ? ? ? ? ? ? ? ? ? = 'CF-1.6'
? ? ? ? ? ?geospatial_lat_min? ? ? ? ? ? ?= -89.875
? ? ? ? ? ?geospatial_lat_max? ? ? ? ? ? ?= 89.875
? ? ? ? ? ?geospatial_lon_min? ? ? ? ? ? ?= 0.125
? ? ? ? ? ?geospatial_lon_max? ? ? ? ? ? ?= 359.875
? ? ? ? ? ?geospatial_vertical_min? ? ? ? = '0.0'
? ? ? ? ? ?geospatial_vertical_max? ? ? ? = '0.0'
? ? ? ? ? ?geospatial_lat_units? ? ? ? ? ?= 'degrees_north'
? ? ? ? ? ?geospatial_lon_units? ? ? ? ? ?= 'degrees_east'
? ? ? ? ? ?geospatial_lat_resolution? ? ? = 0.25
? ? ? ? ? ?geospatial_lon_resolution? ? ? = 0.25
Dimensions:
? ? ? ? ? ?time? ? ? = 1
? ? ? ? ? ?latitude? = 720
? ? ? ? ? ?longitude = 1440
? ? ? ? ? ?nv? ? ? ? = 2
Variables:
? ? time? ? ? ? ? ??
? ? ? ? ? ?Size:? ? ? ?1x1
? ? ? ? ? ?Dimensions: time
? ? ? ? ? ?Datatype:? ?single
? ? ? ? ? ?Attributes:
? ? ? ? ? ? ? ? ? ? ? ?axis? ? ? ? ? = 'T'
? ? ? ? ? ? ? ? ? ? ? ?calendar? ? ? = 'gregorian'
? ? ? ? ? ? ? ? ? ? ? ?long_name? ? ?= 'Time'
? ? ? ? ? ? ? ? ? ? ? ?standard_name = 'time'
? ? ? ? ? ? ? ? ? ? ? ?units? ? ? ? ?= 'days since 1950-01-01 00:00:00'
? ? ? ? ? ? ? ? ? ? ? ?bounds? ? ? ? = 'climatology_bnds'
? ? climatology_bnds
? ? ? ? ? ?Size:? ? ? ?2x1
? ? ? ? ? ?Dimensions: nv,time
? ? ? ? ? ?Datatype:? ?single
? ? latitude? ? ? ??
? ? ? ? ? ?Size:? ? ? ?720x1
? ? ? ? ? ?Dimensions: latitude
? ? ? ? ? ?Datatype:? ?single
? ? ? ? ? ?Attributes:
? ? ? ? ? ? ? ? ? ? ? ?axis? ? ? ? ? = 'Y'
? ? ? ? ? ? ? ? ? ? ? ?bounds? ? ? ? = 'lat_bnds'
? ? ? ? ? ? ? ? ? ? ? ?long_name? ? ?= 'Latitude'
? ? ? ? ? ? ? ? ? ? ? ?standard_name = 'latitude'
? ? ? ? ? ? ? ? ? ? ? ?units? ? ? ? ?= 'degrees_north'
? ? ? ? ? ? ? ? ? ? ? ?valid_min? ? ?= -89.875
? ? ? ? ? ? ? ? ? ? ? ?valid_max? ? ?= 89.875
? ? lat_bnds? ? ? ??
? ? ? ? ? ?Size:? ? ? ?2x720
? ? ? ? ? ?Dimensions: nv,latitude
? ? ? ? ? ?Datatype:? ?single
? ? ? ? ? ?Attributes:
? ? ? ? ? ? ? ? ? ? ? ?comment = 'latitude values at the north and south bounds of each pixel.'
? ? ? ? ? ? ? ? ? ? ? ?units? ?= 'degrees_north'
? ? longitude? ? ? ?
? ? ? ? ? ?Size:? ? ? ?1440x1
? ? ? ? ? ?Dimensions: longitude
? ? ? ? ? ?Datatype:? ?single
? ? ? ? ? ?Attributes:
? ? ? ? ? ? ? ? ? ? ? ?axis? ? ? ? ? = 'X'
? ? ? ? ? ? ? ? ? ? ? ?bounds? ? ? ? = 'lon_bnds'
? ? ? ? ? ? ? ? ? ? ? ?long_name? ? ?= 'Longitude'
? ? ? ? ? ? ? ? ? ? ? ?standard_name = 'longitude'
? ? ? ? ? ? ? ? ? ? ? ?units? ? ? ? ?= 'degrees_east'
? ? ? ? ? ? ? ? ? ? ? ?valid_min? ? ?= 0.125
? ? ? ? ? ? ? ? ? ? ? ?valid_max? ? ?= 359.875
? ? lon_bnds? ? ? ??
? ? ? ? ? ?Size:? ? ? ?2x1440
? ? ? ? ? ?Dimensions: nv,longitude
? ? ? ? ? ?Datatype:? ?single
? ? ? ? ? ?Attributes:
? ? ? ? ? ? ? ? ? ? ? ?comment = 'longitude values at the west and east bounds of each pixel.'
? ? ? ? ? ? ? ? ? ? ? ?units? ?= 'degrees_east'
? ? nv? ? ? ? ? ? ??
? ? ? ? ? ?Size:? ? ? ?2x1
? ? ? ? ? ?Dimensions: nv
? ? ? ? ? ?Datatype:? ?int32
? ? ? ? ? ?Attributes:
? ? ? ? ? ? ? ? ? ? ? ?comment? ?= 'Vertex.'
? ? ? ? ? ? ? ? ? ? ? ?long_name = 'Number of cell vertices'
? ? ? ? ? ? ? ? ? ? ? ?units? ? ?= '1'
? ? crs? ? ? ? ? ? ?
? ? ? ? ? ?Size:? ? ? ?1x1
? ? ? ? ? ?Dimensions:?
? ? ? ? ? ?Datatype:? ?int32
? ? ? ? ? ?Attributes:
? ? ? ? ? ? ? ? ? ? ? ?comment? ? ? ? ? ? = 'This is a container variable that describes the grid_mapping used by the data in this file. This variable does not contain any data; only information about the geographic coordinate system.'
? ? ? ? ? ? ? ? ? ? ? ?grid_mapping_name? = 'latitude_longitude'
? ? ? ? ? ? ? ? ? ? ? ?inverse_flattening = 298.257
? ? ? ? ? ? ? ? ? ? ? ?semi_major_axis? ? = 6378136.3
? ? adt? ? ? ? ? ? ?
? ? ? ? ? ?Size:? ? ? ?1440x720x1
? ? ? ? ? ?Dimensions: longitude,latitude,time
? ? ? ? ? ?Datatype:? ?int32
? ? ? ? ? ?Attributes:
? ? ? ? ? ? ? ? ? ? ? ?_FillValue? ? = -2147483648
? ? ? ? ? ? ? ? ? ? ? ?long_name? ? ?= 'Averaged Absolute Dynamic Topography 2021/07'
? ? ? ? ? ? ? ? ? ? ? ?standard_name = 'sea_surface_height_above_geoid'
? ? ? ? ? ? ? ? ? ? ? ?cell_methods? = 'time: mean within years'
? ? ? ? ? ? ? ? ? ? ? ?scale_factor? = 0.0001
? ? ? ? ? ? ? ? ? ? ? ?units? ? ? ? ?= 'm'
? ? ? ? ? ? ? ? ? ? ? ?coordinates? ?= 'longitude latitude'
? ? ? ? ? ? ? ? ? ? ? ?grid_mapping? = 'crs'
讀取代碼：
madt = ncread('dt_global_allsat_madt_h_y2021_m07.nc','adt');
% ncdisp('dt_global_allsat_madt_h_y2021_m07.nc')
msla = ncread('dt_global_allsat_msla_h_y2021_m07.nc','sla');
lon = 0.25:0.25:360;
lat = -89.75:0.25:90;
[lon , lat] = meshgrid(lon ,lat);
LLZ.lon = lon';
LLZ.lat = lat';
LLZ.rg = madt;
subplot(3,1,1),rg_plot(LLZ),title('dt global allsat madt h y2021 m07')
h = colorbar;
set(get(h,'title'),'string','m')
LLZ.rg = msla;
subplot(3,1,2),rg_plot(LLZ),title('dt global allsat msla h y2021 m07')
h = colorbar;
set(get(h,'title'),'string','m')
LLZ.rg = msla-madt;
subplot(3,1,3),rg_plot(LLZ),title('difference between madt and msla')
h = colorbar;
set(get(h,'title'),'string','m')
中國(guó)海域部分結(jié)果：
歡迎交流學(xué)習(xí)！
參考文獻(xiàn)：聯(lián)合衛(wèi)星重力、衛(wèi)星測(cè)高和海洋資料研究中國(guó)南海 海平面變化，馮偉等，2012.
