AVISO卫星测高格网数据下载与读取
法国 AVISO(Archiving, Validation and Interpretation of Satellite Oceanographic)提供的格网化数据(MSLA).该数据融合了 T/P, Jason-1, ERS 和 ENVISAT 等多颗卫星 的测高资料. 数据时间分辨率为 7 天, 空间分辨率为 1°×1°, 时间跨度为 1993 年 1 月至 2021年7月. 数据经过了必要的标准改正, 如电离层延迟改正、对流层干湿分量改正、固体潮和海潮改正、海洋负荷潮汐 改正、极潮改正、电磁偏差改正、仪器改正和反变气 压计改正等。下面是链接的网址:
https://www.aviso.altimetry.fr/en/data/products/sea-surface-height-products/global/gridded-sea-level-anomalies-mean-and-climatology.html#c10358.
进入官网后,下面的界面就是每月平均的格网数据(数据下载需要提前注册账号)
选择第二个文件夹,即可进入下载界面:
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进行对比:
文件内部的信息:
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')
全球的读取结果
中国海域部分结果:
欢迎交流学习!
参考文献:联合卫星重力、卫星测高和海洋资料研究中国南海 海平面变化,冯伟等,2012.
AVISO卫星测高格网数据下载与读取相关推荐
- MATLAB根据经纬度格网数据画地图
worldmap([20 60],[70 140])%描绘出中国地图的轮廓 %陆地水格网数据显示 geoshow(gca,data,refvec,'DisplayType','texturemap') ...
- Fledermaus和Surfer合作处理离散和格网数据
最近因为涉及到巨量的点数据处理,使用了一下Fledermaus,它的能力确实出众. Fledermaus是IVS 3D公司出品的一套功能强大的交互式3维数据可视化系统.它可以帮助用户完成包括海洋(海岸 ...
- OCO-2、OCO-3卫星数据下载、读取与可视化(适合初学者)
OCO-2.OCO-3卫星数据下载.读取与可视化(适合初学者) 有关OCO系列数据的下载方法已经有一些同行做了比较详细的解答了,这篇博文主要记录我读取及可视化此数据的过程.因为博主自身编程能力不强并且 ...
- CTU-UHB Intrapartum Cardiotocography Database 产时CTG信号数据下载与读取(2021.7.9 新增阿里云盘下载)
目录 CTU-CHB Abstract Data Description Acknowledgements 正文 数据读取 MATLAB代码 数据整理 注意! CTU-CHB 首先说明一下什么是CTU ...
- 点云数据在Matlab生成格网数据-2021-10-18
博主10-31日的更新链接如下:就是原模原样的代码,但是把R空间参考改了,可以储存地理坐标或者投影坐标了 点云生成栅格数据补充+带有地理坐标,2021-10-31_~追风筝的猫的博客-CSDN博客 点 ...
- c++实现快速气象格网数据成图: 风力图, 卫星云图,雷达反射图, 雨量图,温度图
输入数据格式说明: 以上DISK和REG开头都为气象云图为红外数据 PRCP为降雨量图 TMP为温度图 U和V开头的为水平和垂直风力分量数据 Z_OTH为雷达反射图 先放成图效果,后方成图算法,效果图 ...
- 香港cors网数据下载
香港卫星定位参考站网由18个平均分布于全港各处的 连续运行参考站 (CORS) 组成(包括16个参考站及2个完整性监测站),可免费下载使用,数据分辨率有1s.5s.30s三种. 1.下载网址 官方提供 ...
- 2019全国大众点评网数据下载获取
大众点评全网数据(324个国内城市,所有分类)美食数据超过1480万条(14,860,209条)全部分类4000千万+. 详细字段说明: "shop_id"(商户ID,唯一.固定不 ...
- Aviso 卫星测高数据处理(1)——数据下载和nc转TIF批处理
1.Aviso卫星测高数据下载: 在研究海平面变化时常用数据就是卫星测高数据,Aviso卫星测高提供的格网数据时间跨度为:1993至今.下载网站:Data: Aviso+ (altimetry.fr) ...
最新文章
- php 调用url field,使用 PHP curl 调用 API,如何传递如下格式的 POSTFIELD?
- 【MFC】工具栏按钮多选效果本
- 《SAS编程与数据挖掘商业案例》学习笔记之十
- vue element-ui 无法渲染表格数据解决
- 面试精讲之面试考点及大厂真题 - 分布式专栏 09 缓存必问:Reids持久化,高可用集群
- 华为留了一手!将继续发布P50、Mate50:搭载麒麟9000...
- 【机器学习】Softmax Regression算法原理与java实现
- iBarcoder如何在给定位置打印顺序标签?
- 技术书籍也香艳-《Head First设计模式》第二版即将出版
- 创客c语言程序设计软件下载,c语言编程软件
- 深入理解GBDT回归算法
- 动画 | 什么是红黑树?(与2-3树等价)
- Jzoj4722 跳楼机
- 遇到【java.lang.NoClassDefFoundError: javax/xml/bind/DatatypeConverter】的解决方案
- 又拍云叶靖:基于Docker的云处理服务平台
- 血氧心率测量仪(带温湿度测量功能)OLED显示
- Hexo个人博客搭建教程
- Qt图片绘图类QPixmap、QImage、QPicture
- 运维小白的python之路(一)
- 《数据结构》C语言版 链表的基本操作实现