GRACE数据处理和绘图
文章平均质量分 66
分享自己在GRACE数据处理和绘图中遇到的问题,并给出可行的解决方案
present1227
博士研究生在读,研究方向:GRACE时变重力场数据处理及其水文应用研究和智能测绘
本号用于总结和分享WX present1227 如有问题和需求,欢迎评论或私信交流
展开
-
ITSG、COST-G、Tongji和WHU Level-2数据产品读取绘图(Matlab)
Institute of Geodesy at Graz University of Technolog(格拉茨理工大学大地测量研究所) 2018版本,最高60阶球谐系数。代码中将所有球谐系数均人为截断至60阶,只对球谐系数做最基本的替换低阶项和扣除平均重力场处理,其余包括GIA改正、滤波处理和泄漏改正等均未涉及;:时变重力场解决方案组合服务解RL01版本,最高96阶球谐系数。:武汉大学 基于卫星间位势差的约束GRACE月重力场模型。:同济大学2022版本,最高96阶球谐系数。原创 2024-06-26 20:03:04 · 357 阅读 · 0 评论 -
绘制全球各大洲典型流域的时间序列图
以CSR RL05球谐系数为例,生成无约束格网解(grid_CSR)、高斯300km滤波解(grid300_CSR)、去相关滤波解(gridP4M6_CSR)和组合滤波解(grid300P4M6_CSR),并以CSR mascon解(C_mas)作为参考。原创 2024-06-25 23:19:08 · 249 阅读 · 0 评论 -
Matlab读取Swarm球谐系数,并绘制EWH全球格网图(存在疑问)
从结果可以看出40阶的Swarm球谐系数反演的噪声太大,即使经过强约束的高斯滤波后效果仍然不理想,不大清楚是Swarm数据本身所致还是代码实现有误?还请有懂的友友帮忙解答下,不胜感激!ICGEM官网下载 COST-G发布的40×40的球谐系数。原创 2024-05-28 22:02:12 · 462 阅读 · 0 评论 -
GRACE滤波数据处理之DDK系列滤波
Kusche 等(2007)从时变重力场模型球谐系数求解出发,对球谐系数的最小二乘解的法方程组应用贝叶斯估计,构造了 DDK 滤波。原创 2024-05-03 11:25:26 · 516 阅读 · 4 评论 -
利用Argo数据分别计算温度、盐度和温盐所造成的比容海平面变化
需下载计算海底压力、海水密度以及其他有关海水信息的“seawater”代码包,还需下载Yan-Ning Kuo编写的用于计算比容海平面变化的代码包(“steric_height_calculation.m”)按照上述盐容和热容海平面变化的计算公式,计算热容时,保持盐度为平均盐度不变;3、读取IPRC温盐数据,并绘制比容海平面变化全球趋势图与盐容、热容以及比容在太平洋某处的时间序列图。原创 2024-04-30 18:13:12 · 1016 阅读 · 0 评论 -
一个脚本文件读取CSR、JPL和GSFC的mascon产品,并绘制流域时序图
JPL:https://podaac.jpl.nasa.gov/dataset/TELLUS_GRAC-GRFO_MASCON_CRI_GRID_RL06.1_V3 (可能需要科学上网)由于下载的JPL mascon数据为“.nc4”文件,于是将其修改为“.nc”文件;并将文件名开头修改为“JPL_”,以与其它机构的mascon文件名对应起来。以2002年4月的全球EWH格网点为例,原创 2024-03-27 11:15:03 · 724 阅读 · 0 评论 -
GRACE/GRACE-FO泄漏误差数据处理之尺度因子法
先用GLDAS水文模型的格网产品经过球谐分析得到球谐系数,再经过截断和高斯滤波或者去相关滤波处理得到滤波后的格网,逐格网点求解滤波前后的差异因子。发现经过尺度因子改正之后,泄漏至海洋的信号明显减弱,但是以长江流域为例,求EWH时间序列发现,乘以GRACE/GRACE-FO经过截断和滤波后的格网数据得到尺度因子改正后的解。3、 将GLDAS滤波前和滤波后的格网按时间序列排列进行最小二乘求尺度因子。2、格网数据转换为球谐(质量),再经过组合滤波得到GLDAS滤波解。尺度因子法的基本原理,直接上图。原创 2024-01-13 21:49:17 · 1324 阅读 · 6 评论 -
结合GLDAS模型估计GRACE/GRACE-FO全球的RMS比率(信噪比)
上述公式出自 Chen J L, Wilson C R, Seo K W. Optimized smoothing of Gravity Recovery and Climate Experiment (GRACE) time‐variable gravity observations[J]. Journal of Geophysical Research: Solid Earth, 2006, 111(B6). 中的公式(6)在JGR: Solid Earth本篇文章网站显示所有引用这篇文章的学者论文。原创 2024-01-10 10:46:10 · 929 阅读 · 2 评论 -
GRACE泄露改正之全球正向建模恢复法(Forward-Modeling)
提前准备好:(1)CSR机构发布的Mascon产品,这里以CSR-Mascon RL06数据为例,准备了从2002年4月至2021年7月 共197个月的数据(实际试验的时候可以并不需要这么多数据)(2)全球陆地格网点的纬度加权值,即代码中的grid_weight。另一方面,南极区域外非冰盖质量变化信号泄漏至南极区域内,形成内泄漏,两者统称为泄漏误差。①M_obs为Mascon解经过球谐分析和500km高斯滤波后的结果,第一次循环开始前,假设M_tru=M_obs,每次循环的过程中,更新M_tru即可;原创 2023-10-23 20:50:00 · 376 阅读 · 6 评论 -
GRACE数据的滤波方法实现(Han、Fan和Wiener)
首先以等效水柱高形式绘制测量信号的阶方差图像,再运用上述前两个公式进行最小二乘拟合,求得拟合参数a、b、c、d,最后利用第三公式计算滤波系数h(注意:对于最大阶数60的球谐产品,上述第一个公式适用于0-20阶,第二个公式适用于31-60阶数)示例代码中滤波半斤r0=200km、r1=300km、选定的次m1=15,绘制滤波系数图如下。原创 2023-10-15 17:38:00 · 611 阅读 · 6 评论 -
绘制60阶无约束解和约束解的球谐系数格网图
①直接从机构数据读取,经过一阶项、二阶项和取平均等预处理操作后的61*121的球谐系数sc矩阵;②例如经过高斯300km之后的等效水柱高,再利用gmt_grid2cs函数转化为球谐系数。2、利用pcolor函数绘制球谐系数格网图。原创 2023-10-13 11:18:43 · 126 阅读 · 0 评论 -
GRACE数据绘制区域格网图
1、首先下载区域的边界经纬度坐标(如下图展示的为全球主要流域的边界经纬度),解压RBvect文件可以找到index用记事本打开,以编号为0011的Changjiang Yangtze流域边界为例进行说明,2、获取到研究区域边界后,即可参照下面程序进行区域格网图的绘制。出图效果为这样,需要人为手动修改标注条至合理的位置。原创 2023-10-12 09:47:52 · 624 阅读 · 20 评论 -
GRACE数据的广义三角帽(TCH)计算不确定度
1、首先整理好需要计算不确定度的GRACE数据,这里以CSR level-2截断至60阶的球谐系数、DDK3和4滤波解、CSR mascon以及JPL mascon数据为例(以JPL mascon数据为参照)[2]姚朝龙, 李琼, 罗志才,等. 利用广义三角帽方法评估GRACE反演中国大陆地区水储量变化的不确定性[J]. 地球物理学报, 2019, 62(3):15.2、根据TCH中的相关公式,利用迭代初值计算噪声协方差矩阵R(n*n)中的相关元素。函数求解在满足约束条件的情况下使得目标函数最小的参数。原创 2023-05-21 19:18:24 · 1411 阅读 · 10 评论