weixin_43047707的博客

Matlab colorbar label标注

matlab 小时分辨率 X轴标注

中国省市区域 shapefile kml 文件格式转换

声明 微信公众号 专栏文集

GIA matlab GRACE aviso

gmt6 地形起伏 cpt 误差棒

印度洋偶极子 IOD json文件 matlab

matlab nc文件 批量处理

全球河流径流数据 GRDC 数据 逐日 逐月

matlab HYCOM 海洋温度 投影转换

matlab函数

GRACE 负荷位移 垂直位移 水平位移

负荷变形 弹性地球 负荷LOVE数 圆盘负荷 地表变形

湖泊水位测高产品、GNSS时间序列产品

Matlab nc文件 时间转换

matlab 科研数据 海平面变化 温度和盐度 比容海平面 GLORYS

matlab GFED数据 森林火灾 hdf5文件

GIA PGR matlab Caron Peltier Geruo level3数据

matlab EOF 经验正交分解 主成分 GRACE

全球降水数据产品 GPCC GPCP PREC/L CMAP

gmt 图框 colorbar

matlab 海平面 格网数据 scatteredInterpolant 北冰洋 重采样 投影转换

GNSS 陆地水储量反演 GRACE 开源代码 Github

LoadDef工具包 python 负荷love数 github

而为了将GLDAS提取的水储量与GRACE能够进行对比，需要进行同样的截断滤波处理。本专栏以GLDAS数据的土壤水为例，分析如何进行球谐展开并滤波。下面是截断滤波后的结果：我们看到截断滤波后，信号强度有所衰减，而随着截断阶数越大，信号越接近真实信号强度。而高斯平滑的半径越大，信号也衰减。

假设有一个均匀采样的时间序列，因此其轨迹矩阵为矩阵的列数为，Y的每个上升斜对角线都是相同的值。在之前的研究中，基于Broomhead和King(1986)提出的方法，使用矩阵Y形成滞后协方差矩阵，并采用主成分分析(PCA)程序进行处理。一种等效但更简洁的替代方法是通过奇异值分解(SVD)直接对Y进行分解:其中的下表代表矩阵大小，是对角矩阵，是正交的，也是滞后协方差矩阵的特征向量。主成分分析中经验正交函数(EOFs)和主成分(PCs)的输出可以表示如下：其中和的第列，以及的第。

目录一、如何在DEM中绘制特定测线的剖面二、如何在其他物理场中绘制特定测线的剖面（以负荷场为例）在GIS中，我们通常需要绘制特定剖面线上的值的变化情况，剖面可以清晰反映出某个物理场具体沿着特定方向的变化情况，比如以下的情况就非常的易于理解（图片来自于Pan et al, JGR）。接下来我们就针对于这样的需求来分析我们可以采用什么方法实现？我将尽可能采用多的方法和软件工具达到这一目的。首先我们先从直观的高程入手。【数据准备】我是直接利用GMT5工具将全球DEM提取出中国地区的DEM。具体的下载地址如下：

国际大地参考框架（ITRF）是国家大地参考系统的具体实现。联系不同的ITRF参考框架的转换关系是至关重要的，它们是基于转换参数，也称赫尔默特参数，它们可以用于将来自两个不同的ITRS系统的观测结果进行统一比较

Atmosphere and Ocean Level1B De-Aliasing (AOD1B)全称为大气和海洋去混叠level 1B产品。在GRACE数据处理时，由于卫星采样的限制，通常需要积累一定时间段的观测数据来进行全球重力场的解算，从而获得高空间分辨率的解，对GRACE而言，通常是7到30天。AOD1B产品旨在作为背景场模型移除重力场中的高频非潮汐变化信号。目前已经发布到AOD1B Release 06版本了。数据下载网站：https://podaac.jpl.nasa.gov/dataset/GR

是如何读取自带的RL05数据，目前数据已经发布到RL06了，因此我将在这里介绍本人改进的可以读取RL06数据的一套代码。我们下载冯伟老师的工具包，我在这里把我调好的可以读取RL05数据的程序放在以下的链接，感兴趣的可以下载调试，后续会上传可以读取最新的RL06数据的程序包。注意：控制文件中需要读取目录下所有GSM的文件名，手动修改是不可能的，可以采用以下的方法：在文件对应得目录下，新建txt文件，里面写入。（3）读取二阶项，C20没有改变，只是修改了C21S21C22S22的函数。2. 控制文件的修改。

与地理坐标不一样，长方体存在特殊的情况，比如需要把其分成上下、左右、前后6部分进行拼接。（1）我读取的文件是csv格式，因此选择文件后，Field Delimiter选择Comma,文件头没有则为0，下面选择的参数默认即可。具体的下载链接为：https://github.com/NWRichmond/xyz2stl.(2)下一步选择输出的路径和文件名，其余参数默认，点击运行即可。其中的 b01.txt为x,y,z三列数据。我在之前的一个专栏发过一个将地理坐标的数据。下面即为得到的stl文件。

此处我使用的是GMT5，对于GMT6，有另外一种方法。参考：https://docs.gmt-china.org/latest/examples/ex007/#gmtplot-7adb3518e066c53838e13ab7c9582743。这里记录一下如何在GMT中添加子图、指北针、比例尺和图例，以防每次全网搜索找不齐全。

都可以达到以下的效果，但是加上grid on命令也无法显示格网，原因是，如果采用上述的第二种方法，我们已经设置’linestyle‘为none，因此是没有办法绘制格网的。因此我们需要更改绘图的函数。如果我们使用pcolor绘制，则是以下的效果：pcolor中的每一个格子的边框都绘制了，因此对于小区域效果还行，但是要是范围过大，则是一片格网，黑黑一片。而且加上命令也没有效果。我们需要绘制一个这样的绘图效果：tick向外，加上格网虚线。但是经过我本人的试验，应该是要加上以下的命令，结果才是正常的。

假设在 10°W至10°E，10°S至10°N 的 5°×5°范围内存在一个 10 mm 的水层。尽管这里有两个分布完全不同的模 型，但是可以发现，在平滑和滤波后它们产生的信号几乎是一样的，也就是说细节基本上已 经丢失了，我们没有办法从观测信号中再去识别原始质量究竟是哪一种分布。那么，这种情 况下，我们不关心信号的细节分布，我们只关心总质量的大小。因而在GRACE观测的信号中，分辨率是大约300 km，在这个分辨率内部，存在众多的信号源，但是我们得到的是总的质量变化情况。因此需要进行精细的信号分离。

GRACE低阶项数据处理数据时间覆盖度

matlab编程时间日期转换格里高力历gregorian转年月日ymd年积日转年月日

GRACE数据处理 DDK1-8滤波 matlab

matlab读取CPC土壤水模型

MATLAB程序 GRACE数据处理 ICGEM网站 gfc数据处理

matlab 3D模型 stl文件 三维显示 地球重力场 陆地水储量变化EGM2008