[QOCA学习笔记]利用mload模块计算地表质量负载

最新推荐文章于 2023-08-01 08:44:57 发布
最新推荐文章于 2023-08-01 08:44:57 发布
阅读量731 收藏 4
点赞数 1
分类专栏： 其他笔记 文章标签： ubuntu 地表质量负载 QOCA
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_42638388/article/details/108808800
版权

在季节性时间尺度上，地表质量负载引起的季节性场地位移是空间大地测量中站点坐标时间序列的主要特征。关于QOCA软件的mload简介可以移步下载《Briefing of mass load caused site displacements》
在处理数据时想按照文献中这样的方法计算安徽省区域的地表质量负载：（参考文献：武曙光. 区域CORS站坐标时间序列特征分析[D]. 2017.）
在这里插入图片描述
首先下载文中说到的需要的数据。其中，大气压数据、土壤湿度和积雪深度数据下载机构是NCEP，均采用NCEP-DOE Reanalysis 2的数据，网址为：NCEP-DOE Reanalysis 2: Summary，在该网站上可以找到数据下载的入口：
在这里插入图片描述
点击“Search all of the NCEP-DOE Reanalysis 2 available at PSL.”，进入根据你所需要的数据进行搜索并下载的界面，这里的分别下载大气压数据每6h、土壤湿度每24h和积雪深度每24h数据。
然后去NOPP-ECCO机构官网下载最新的非潮汐海洋负载的洋底压力模型数据，官网地址为：V4r4 User Guide这里的网址是最新的数据用户手册，下载数据之前可以先看清楚用户手册。下载这个数据之前首先得申请一个账户，申请了账户之后才会给用户分配密码，才可以进入数据下载网址下载数据。
在这里插入图片描述
进入这个数据下载中心，找到想要的年份，然后找到洋底压力模型文件（时间分辨率为12h）：
在这里插入图片描述
在这里插入图片描述
将相应年份的数据下载之后整理好，准备放入QOCA软件里进行解算。
【未完待续…】

hr11asdfg
关注
  • 1
    点赞
  • 踩
  • 4
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
专栏目录
智能合约_智能合约编写之Solidity运行原理
weixin_39700625的博客
12-24 887
3月25日，BSN第二次开发者大赛正式启动，本次大赛以“编写基于多种底层框架的智能合约”为主题，开发者可基于FISCO BCOS等主流底层框架，结合业务场景设计、开发并部署智能合约。为了让大家更好上手智能合约开发，区块链服务网络发展联盟与FISCO BCOS开源社区共同推出“7个课时速成智能合约全能型开发”系列教程，助力开发者学习并熟悉合约开发，轻松应对此次大赛并拔得头筹。系列专题 | 超话区块链...
zde:ZDE 1.6的重构源代码（Z系统显示编辑器）-Source code editor
03-24
ZDE 1.6（Z系统显示编辑器）重构源 以的源代码为指南，我重新构造了的源代码。 该源已与以下组件组装在一起： Al Hawley的ZMAC：按原样组装。 MicroSoft的M80：重命名为ZDE16.MAC，取消注释前两行并进行组装。 使用RELHEX创建ZDE16.HEX。 ZASM（Cromemco的ASMB）：重命名为ZDE16.Z80并进行组装。 使用RELHEX创建ZDE16.HEX。 使用MLOAD创建ZDE16.COM。 对于ZDE，还有两个我还没有弄清楚的例程（尚未）。 但是大多数都被怀疑了。 ZDE 1.7 我修复了v1.6中存在的“在ZSDOS下不保留大于单个范围的文件的时间戳”的错误。 现有的ZDENST16.COM程序将与1.7一起使用，以为您的终端机和打印机设置程序。
Mass load caused site displacements.pdf
08-03
该文件中包含GNSS坐标时间序列QOCA（Quasi-observation Combination Analysis）软件mload模块的简介（英文版）包括驱动文件中各参数的释义from QOCA online tutor
QOCA指南.zip
07-16
共包含5篇文档，分别介绍了QOCA中analyze_tseri、Analyze_strain、st_filter、mstoke等模块的一般应用。其中analyze_tseri模块应用较多，文档介绍比较详细，有各个参数的设置。
[QOCA学习笔记]利用QOCA软件analyze_tseri模块进行时间序列分析
无
08-05 1327
操作环境：win10+虚拟机ubuntu下QOCA（ver 1.33） 数据：通过解算得到的安徽省CORS参考站坐标时间序列（neu格式） 对于QOCA软件的介绍，在其官网https://qoca.jpl.nasa.gov/上可以查看相关的描述。其analyze_tseri模块可以用来分析GNSS坐标时间序列。对该模块的详细介绍可以在这里下载：analyze_tseri analysis.pdf 与之前pca模块相似，在调用analyze_tseri模块之前亦需要仔细查阅其驱动文件的设置。下面简单介绍一次
利用GRACE和地表质量模型计算地球弹性负荷变形实验笔记
最新发布
weixin_43047707的博客
08-01 935
地球负荷变形 love数 格林函数 python GRACE mascon 负荷源
Principal component analysis.pdf
08-02
该文件中包含GNSS坐标时间序列QOCA（Quasi-observation Combination Analysis）软件PCA模块的简介（英文版）包括驱动文件中各参数的释义from QOCA online tutor
QOCA软件大地测量专业、地球物理专业必备的专业软件
04-17
地学必备，大地测量专业、地球物理专业必备的专业软件
Google Earth Engine——NCEP-DOE Reanalysis 2项目正在使用最先进的分析/预报系统，利用1979年至前一年的过去数据进行数据同化。
qq_31988139的博客
10-04 717
NCEP-DOE Reanalysis 2 project is using a state-of-the-art analysis/forecast system to perform data assimilation using past data from 1979 through the previous year. NCEP-DOE Reanalysis 2项目正在使用最先进的分析/预报系统，利用1979年至前一年的过去数据进行数据同化。 Dataset Availability 1979
QOCA软件安装（附安装后出现缺少libg2c.so.0报错的解决办法）
无
05-28 2552
QOCA（Quasi-observation Combination Analysis），由加州理工学院喷气推进实验室设计和开发的。其是将各种松散约束的大地测量站点坐标和速度解(作为准观测)相结合，获得地壳形变信息的软件包。QOCA对任何非商业用户都是免费的。许多大地测量数据分析研究人员用QOCA作后处理软件包。目前，它可以结合空间大地准观测(GPS, VLBI, SLR，…等)和大地测量准观测(EDM，三角测量，水准测量，…等等)。它有可能将SAR数据与重力、地震活动和地面运动数据结合起来。 QOCA软件
[QOCA学习笔记]利用QOCA软件PCA模块进行共模误差分析
无
08-03 1302
操作环境：win10+虚拟机ubuntu下QOCA（ver 1.33） 数据：通过解算得到的安徽省CORS参考站坐标时间序列（neu格式） QOCA软件是优秀的GNSS坐标时间序列分析软件，其简介及PCA模块的使用方法在QOCA online tutor上已经有详细介绍： 下面简单介绍利用PCA模块对时间序列的共模误差进行分析的过程（具体PCA的介绍请下载相关文献） ...
matlab代码绿颜色表示-Registration:登记
06-03
matlab代码绿颜色表示登记 用于大脑发育中尺度研究的时空 3D 图像配准 用法： 将包含脚本的文件夹添加到 Matlab 路径。 为此，打开 Matlab 并在 Matlab 中导航到存储文件夹Registration位置。 右键单击文件夹Registration ，然后从下拉菜单中选择Add to Path ，然后Selected Folders and Subfolders 。 导入大脑图像 在 Matlab 中，导航到包含单个 3D 大脑样本的 z 堆栈图像的文件夹。 进入文件夹后，在 Matlab 命令窗口中执行以下命令： IM = mload( ' *.tif ' , 3 , 1 / 4 ); 这里的第一个参数是一个正则表达式，它允许从当前文件夹中选择文件。 此处'*.tif'表示我们选择具有任何名称 ( * ) 的 TIFF 格式 ( .tif ) 文件。 我们选择所有文件，因为我们当前所在的文件夹仅包含单个大脑样本的文件。 有关正则表达式的更多信息，请访问 此处的第二个参数 ( 3 ) 表示与 TIFF 图像的原始分辨率相比，我们将图像分辨率粗化了多少。 例如，如果
智能合约源代码TUSD.pdf
02-27
智能合约源代码 智能合约源代码TUSD /** *Submitted for verification at Etherscan.io on 2018-12-18 */ pragma solidity ^0.4.23; // This is the proxy contract for the TrustToken Registry // File: contracts/Proxy/Proxy.sol /** * @title Proxy * @dev Gives the possibility to delegate any call to a foreign implementation. */ contract Proxy { /** * @dev Tells the address of the implementation where every call will be delegated. * @return address of the implementation to which it will be delegated */ function implementation() public view returns (address); /** * @dev Fallback function allowing to perform a delegatecall to the given implementation. * This function will return whatever the implementation call returns */ function() external payable { address _impl = implementation(); require(_impl != address(0), "implementation contract not set"); assembly { let ptr := mload(0x40) calldatacopy(ptr, 0, calldatasize) let result := delegatecall(gas, _impl, ptr, calldatasize, 0, 0) let size := returndatasize returndatacopy(ptr, 0, size) switch result case 0 { revert(ptr, size) } default { return(ptr, size) } } } } // File: contracts/Proxy/UpgradeabilityProxy.sol /** * @title UpgradeabilityProxy * @dev This contract represents a proxy where the implementation address to which it will delegate can be upgraded */ contract UpgradeabilityProxy is Proxy { /** * @dev This event will be emitted every time the implementation gets upgraded * @param implementation representing the address of the upgraded implementation */ event Upgraded(address indexed implementation); // Storage position of the address of the current implementation bytes32 private constant implementationPosition = keccak256("trueUSD.proxy.implementation"); /** * @dev Tells the address of the current implementation * @return address of the current implementation */ function implementation() public view returns (address impl) { bytes32 position = implementationPosition; assembly { impl := sload(position) } } /** * @dev Sets the address of the
学术文章绘图常用颜色搭配(附RGB值)
热门推荐
无
02-22 13万+
参考博客： 添加链接描述 在其推荐的Color Wheel网站上（网址为：添加链接描述）找到了以下一些比较常用的颜色搭配，并附上其对应的rgb值。 ...
小木虫网站无法登陆（账号突然被封禁）的解决办法
无
05-26 8044
去年下半年开始很多小木虫账号都被封禁了，而在用户这边会发现，首次登陆该论坛时，网页会一片空白，之后打开任何一个小木虫的网站都是空白，或者是显示找不到网页的404错误，。偶然一次在百度贴吧上看到才知道原来是账号被封了，贴吧里也有人提出了解决的办法： 首先要借个账号，找同学帮忙申诉，结果他的号是新的还没办法发帖子…就很迷惑。后来看到贴吧里有人说可以给小木虫发邮件申诉账号，邮箱是：[email protected] 经过实测这个方法效率挺高，我前一天晚上发的邮件第二天中午就回复了，只是理由有点无语。 于是乎我
[学习笔记]周跳探测
无
12-10 7091
GNSS观测中的周跳探测 文章目录GNSS观测中的周跳探测1. 周跳定义：2. 周跳探测与修复方法2.1 基于观测值随时间变化规律的方法：2.2 基于不同观测值组合的方法：2.3 基于观测值估值残差的方法3. 常用的周跳探测方法3.1 双频码相组合法3.2 双频电离层残差法3.3 GF组合和MW组合各自的探测盲区：4. GAMP中探测周跳 1. 周跳定义： 整周跳变是载波相位观测值特有的问题 周跳：引起整周计数暂时中断，发生整周跳变。 原因主要：卫星信号被某障碍物阻挡而无法到达接收机，或者由于外界干扰或接
python3.8安装cartopy使用报错：DLL load failed while importing trace
无
04-22 6309
背景 python3.8环境下根据博客中的步骤安装cartopy后出现报错： 安装cartopy参考博客： Cartopy安装教程 其主要思路是在官方网站上下载自己所需要的whl文件到电脑，再利用pip install进行安装，本人在成功运行各pip install后，准备画图时，却发现不能用，报错为:ImportError: DLL load failed while importing trace: 找不到指定的模块。 最后定位在这里： Traceback (most recent call last.
空间直角坐标转换成大地坐标matlab代码（附精度）
无
05-29 5624
在实际生产中，通常需要进行空间直角坐标和大地坐标之间的相互转换。精度估计是大地测量数据处理的一项内容，而目前使用的的坐标转换软件中一般都没有提供空间直角坐标与大地坐标之间转换的精度。其中，空间直角坐标至大地坐标的转换常用的迭代法可以表示为： 在文章《空间直角坐标与大地坐标的精度关系式_徐思达》中，作者论通过误差传播律推导了由空间直角坐标转到大地坐标的精度，分析了严密公式和简便公式的精度转换结果的差异，发现两者只有微小差异，所以可以用简便公式。根据其提供的公式，用matlab编写程序实现。 function
参考GMT中文社区绘制区域CORS站的站点分布和速度场图（脚本）
无
10-18 5396
软件：GMT5.4 操作环境：win10虚拟机下的ubuntu16.04系统 在学习绘制区域CORS站的站点图以及速度场时，参考GMT中文社区里贴出来的一些脚本文件，结合自己的数据进行了一些修改。 1.绘制站点分布图 参考脚本网址（内含各省会的经纬度）：https://gmt-china.org/example/ex004/ 在绘制自己所需的区域站点分布图时，租要自己制作该区域内各个站点的经纬度...
erc20的签名值为629680c82f82ad9ea02704a988d70565741def2d93ec95efe07d8bd6976340e15b87a837d96406516436eef6d88491598e9c4f1c15aaf7260db8ef25af5b5851c，求erc20的permit的r，s，v的值
06-03
根据您提供的签名值，可以得到以下信息： ``` address owner = 0x1234567890123456789012345678901234567890; address spender = 0x0987654321098765432109876543210987654321; uint256 value = 100; uint256 nonce = 1; uint256 deadline = 1634028589; ``` 下面是使用Solidity代码计算r，s，v的值的示例： ```solidity function splitSignature(bytes memory signature) private pure returns (uint8 v, bytes32 r, bytes32 s) { require(signature.length == 65, "Invalid signature length"); assembly { r := mload(add(signature, 32)) s := mload(add(signature, 64)) v := byte(0, mload(add(signature, 96))) } if (v < 27) { v += 27; } require(v == 27 || v == 28, "Invalid signature value"); return (v, r, s); } function getPermitRSV(bytes32 domainSeparator, bytes32 structHash, bytes memory signature) private pure returns (uint256 r, uint256 s, uint8 v) { (v, bytes32 rBytes, bytes32 sBytes) = splitSignature(signature); r = uint256(rBytes); s = uint256(sBytes); bytes32 digest = keccak256(abi.encodePacked("\x19\x01", domainSeparator, structHash)); address signer = ecrecover(digest, v, rBytes, sBytes); require(signer != address(0) && signer == owner, "Invalid signature"); return (r, s, v); } function getPermitSignature() public view returns (uint256 r, uint256 s, uint8 v) { bytes32 domainSeparator = keccak256(abi.encode( EIP712DOMAIN_TYPEHASH, keccak256(bytes("My Token")), keccak256(bytes("1.0")), chainId, address(this) )); bytes32 structHash = keccak256(abi.encode( PERMIT_TYPEHASH, owner, spender, value, nonces[owner]++, deadline )); bytes memory signature = hex"629680c82f82ad9ea02704a988d70565741def2d93ec95efe07d8bd6976340e15b87a837d96406516436eef6d88491598e9c4f1c15aaf7260db8ef25af5b5851c"; return getPermitRSV(domainSeparator, structHash, signature); } ``` 其中，splitSignature函数用于将签名值拆分成v，r，s三个部分；getPermitRSV函数用于从签名值中获取r，s，v的值，并对签名值进行验证；getPermitSignature函数则是调用getPermitRSV函数，返回r，s，v的值。 注意，上述代码中的EIP712DOMAIN_TYPEHASH和PERMIT_TYPEHASH需要根据具体的ERC20代币合约进行修改。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 4
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值