- 博客(59)
- 资源 (8)
- 收藏
- 关注
RSS订阅
原创 学术文章绘图常用颜色搭配(附RGB值)
参考博客：添加链接描述在其推荐的Color Wheel网站上（网址为：添加链接描述）找到了以下一些比较常用的颜色搭配，并附上其对应的rgb值。...
2020-02-22 10:35:24
142787
5
原创 [学习笔记]python绘制图中图（绘制站点分布图）
在绘制站点分布图时，有时需要采用图中图的方式，以便于在一张图中尽可能多的表达信息。此处记录一下利用python matplotlib绘制图中图的脚本，方便然后查询。该汇图脚本中包含以下数据：首先将CMONOC站点以及穿刺点分布画上，再通过plt.axes在图上再加一块画布绘制小范围的地图，接着在小范围的地图上标点。为了方便在大范围地图上找到图中图的位置，也要在相应位置上用红框标记下，这里采用了plt.vlines和plt.hlines两者结合的方法。注：此处仅提供绘图脚本，其中用到的站点坐标文件部分因一些原
2023-11-02 10:32:11
1193
6
原创 [电离层建模学习笔记]开源程序M_GIM学习记录
程序生成的观测值以mat格式保存在./OBS/regional/23001/目录下，精密轨道保存在./SP3，计算得到的P4平滑值保存在./P4/regional/GPS/23001/，最终最小二乘估计的结果保存在./M_Result/，接下来根据结果计算区域格网VTEC，并输出文件。该程序在M_DCB(Jin et al., 2012)的基础上发展而来，是便捷实用的电离层建模软件。本文记录了该程序的学习过程，该程序在运行前还有一些小的错误需要改正，在此也做一个记录，以备日后查询。
2023-06-16 16:11:36
3953
28
原创 [学习笔记]C/C++和python程序编译成Linux下可执行文件
利用cmake生成MakeFile将C/C++程序编译成可执行文件，记录一下编译过程
2022-06-16 10:56:26
2397
1
原创 [GAMP学习笔记]计算STEC程序中遇到的一些问题小结
计算STEC程序中遇到的一些问题小结文章目录计算STEC程序中遇到的一些问题小结1. 相关知识背景1.1 根据接收机和卫星zyx坐标穿刺点计算原理1.2 方法一计算穿刺点：1.3 根据三角关系计算穿刺点2. 关于两种不同方法计算穿刺点得到的高度角、穿刺点经纬度的差异3. 计算VTEC考虑经度改正对结果的影响3.1 根据穿刺点通过插值计算GIM 中VTEC方法3.2 两种方法计算VTEC的造成的差异1. 相关知识背景1.1 根据接收机和卫星zyx坐标穿刺点计算原理1.2 方法一计算穿刺点：首
2022-05-06 17:55:35
3348
3
原创 python3.8安装cartopy使用报错：DLL load failed while importing trace
背景python3.8环境下根据博客中的步骤安装cartopy后出现报错：安装cartopy参考博客：Cartopy安装教程其主要思路是在官方网站上下载自己所需要的whl文件到电脑，再利用pip install进行安装，本人在成功运行各pip install后，准备画图时，却发现不能用，报错为:ImportError: DLL load failed while importing trace: 找不到指定的模块。最后定位在这里：Traceback (most recent call last.
2022-04-22 20:24:48
6606
13
原创 [MATLAB学习笔记]采用快速傅里叶变换求时间序列的周期项
[MATLAB学习笔记]采用快速傅里叶变换求时间序列的周期项1. 背景现有长度为11年的5个时间序列，为某拟研究对象的5个参数。现计划先通过快速傅里叶变换求每个系数序列的显著周期项，再分别按照傅里叶序模型列进行拟合，得到系数和时间的模型，即达到一个效果：输入一个时间，可以得到对应的5个参数。2. 实现2.1 序列中含有缺失天数，即为非等间隔采样需先对序列进行插值计算，不同插值方法对结果的影响这里先不讨论。先用matlab中自带的spline函数对数据进行三次样条插值，得到完整序列。（实际上，本文
2022-04-13 20:41:49
5253
3
原创 [学习笔记]（科研作图）记录一个用于生成颜色表+颜色进制转换的网站
最近发现一个非常实用的网站，可以用来转换颜色RGB的表示方式，并且可以根据提供的两个颜色生成对应的颜色表，在科研绘图中非常使用，先看一下它的基本用法.网址：ColorHexa1. RGB颜色表示转换输入任何一种颜色，该网站可以给出该颜色的详细信息，甚至有和该颜色搭配的颜色，可以说非常详细了，比如我这里输入了RGB 131 60 64，点击搜索，就可以出现这个颜色的详细信息，以及各种进制的表示方法：下方给出了与此颜色搭配的其他颜色：2. 生成自定义的颜色表有时候我们需要生成颜色表来表示
2021-12-28 22:21:41
717
1
原创 [python学习笔记]利用python basemap绘制全球电子含量分布图
[python学习笔记]利用python basemap绘制全球电子含量分布图1. 绘图前准备1.1 需要安装的库numpymatplotlibmpl_toolkits1.2 其它函数在GitHub上找到一个GNSSpy库，里面可以实现RINEX以及IONEX等格式文件下载、读取功能，比较实用，下载地址： GNSSpy为了方便使用，我并没有直接安装这个库（一些函数还是需要修改的），只在这个库基础上做出一些修改方便自己用。在本次绘制电子含量图过程中，只用到了里面readFile.py里面的
2021-12-28 22:01:31
3128
22
原创 [学习笔记]天线相位中心改正
天线相位中心改正文章目录天线相位中心改正1. 相关概念：2. 天线相位中心改正2.1 卫星天线相位中心改正2.2 接收机天线相位中心改正3. IGS提供的天线相位改正文件3.1 绝对天线相位中心模型3.2 atx文件介绍4. GAMP中相关函数：4.1 读取atx文件：（RTKLIB/GAMP）4.2 天线相位改正4.2.1 接收机天线相位改正4.2.2 卫星天线相位改正1. 相关概念：天线参考点：Antenna Reference Point, ARP天线瞬时相位中心： Instantaneou
2021-12-14 15:44:51
6718
7
原创 [学习笔记]周跳探测
GNSS观测中的周跳探测文章目录GNSS观测中的周跳探测1. 周跳定义：2. 周跳探测与修复方法2.1 基于观测值随时间变化规律的方法：2.2 基于不同观测值组合的方法：2.3 基于观测值估值残差的方法3. 常用的周跳探测方法3.1 双频码相组合法3.2 双频电离层残差法3.3 GF组合和MW组合各自的探测盲区：4. GAMP中探测周跳1. 周跳定义：整周跳变是载波相位观测值特有的问题周跳：引起整周计数暂时中断，发生整周跳变。原因主要：卫星信号被某障碍物阻挡而无法到达接收机，或者由于外界干扰或接
2021-12-10 22:20:41
8753
6
原创 [学习笔记]GMT6.2在Win10下使用（Batch/Bash）
[学习笔记]GMT6.2在Win10下使用（Batch/Bash）1. GMT6.21.1 安装下载地址：[gmt-6.2.0-win64.exe][http://mirrors.ustc.edu.cn/gmt/bin/gmt-6.2.0-win64.exe]安装方法：请参考GMT中文社区在线手册中详细介绍[Windows 下安装 GMT][https://docs.gmt-china.org/6.2/install/windows/]1.2 配置Win系统下安装比较简单，点击ex
2021-10-21 17:24:34
2796
4
原创 [学习笔记]Rtklib单点定位中电离层相关函数
[学习笔记]Rtklib单点定位中电离层相关函数1. readtec函数所在文件：ionex.c主要参数：* read ionex ionospheric tec grid file* args : char *file I ionex tec grid file* (wind-card * is expanded)* nav_t *nav IO navigati
2021-10-15 22:19:16
1868
2
原创 [Matlab学习笔记]Matlab中删除矩阵指定行（指定列）
在用matlab编程进行GNSS粗差提取时，想实现对某矩阵删除指定行（指定列）的命令，下文中描述的均为行。通过百度，先明白对于一行数据而言，删除指定位置的元素的原理。主要思路是先提取要删除的位置，再对该位置数据进行删除。先通过以下范例说明上述思路：单行数据：删除指定位置的元素data='1999-08-25T12:00:00';OtherCharIndex=find( (data<'0'|data>'9') & data~='-');OtherCharIndex=[Other
2021-04-06 10:11:24
7865
原创 [GNSS时间序列分析笔记]Ubuntu下制作站点列表
在分析GNSS坐标时间序列时，如果研究区域为大范围的网络，那么站点往往会非常多。在分析时常常会用到站点名称列表，或者在写脚本进行批量处理时用到横向的站点名称，私以为若是站点个数在10个以下，可以采用手动输入的方法制作站点名称列表，但是超过这个数字甚至上百个上千个站点这样制作已经不太可能。这里给处本人在制作站点名称列表时总结出来的一些小Tips，以备日后查用。首先将所有站点的文件放进一个文件夹内，拷入ubuntu系统下。示例中，我用的是neu格式文件，其他xyz文件格式都行，只要文件名中包含了站点名称都可
2021-03-26 16:33:25
1063
2
原创 [Hector学习笔记]GNSS时间序列处理软件Hector使用备忘（批处理脚本）
由Machiel Bos 和Rui Fernandes共同研发的Hector软件可以利用时间相关噪声来估计时间序列线性趋势，是与优秀的GNSS时间序列处理软件之一。软件安装教程：Hector软件安装[附加软件环境变量设置]Hector软件用户手册中包含了教程，示例文件等等，并且有非常详细的说明。在教程中，大多采用“设置驱动文件+在终端中逐行输入代码”进行计算的方式进行解算，这适用于新手，或者是需要处理的站点较少的情况。对于超过5个站点的时间序列的分析，如果还用逐个修改驱动文件和终端逐个输入的方式，未免太
2021-03-25 17:04:20
3950
10
原创 win10下安装Java JDK1.8以及环境变量配置
JDK简介等详细信息可以直接百度，这里仅将安装过程做个记录，用以备忘。JAVA JDK1.8可以在官方网站上进行下载，贴上网址：Java SE Development Kit 15.0.1 ，我下载的文件为：此时点击运行该exe文件，可以在安装向导中设置自己的安装目录，也可以选择装在C盘的默认目录下，重要的是此时应该记录一下安装位置，这里分别对jdk和jre设置安装路径，我的安装位置在：安装向导里安装完成之后还需要进行环境变量配置，右键“此电脑”选择“属性”，选择“高级系统设置”：选择“环境变
2020-11-14 18:38:33
323
原创 [Ubuntu学习笔记]解决虚拟机ubuntu系统内存不足，无法正常开机进入系统的问题
在vm虚拟机上安装的ubuntu16.04系统，一直用来处理数据、画图等，导致内存不足，之前系统报告过这样的错误，但是没有清理过，以至于这次开机无法进入系统，开机后显示输入账号密码，然后就是黑屏，无法正常进入系统，终端输入reboot重启的时候发现开机过程中也有这样的报错，显示机器没有多余的内存。对于虚拟机而言，可以在关机之后编辑虚拟机的内存（图中我的已经修改好了）：有需要动手修改的小伙伴可以参考这篇博客：Vm虚拟机扩展Ubuntu系统磁盘空间之后需要从客户机操作系统内部对磁盘重新进行分区和扩展文件
2020-09-30 10:28:34
6038
原创 [Ubuntu学习笔记]由于操作失误造成ubuntu16.04下运行sudo 出现sudo: /usr/bin/sudo must be owned by uid 0 and have the s报错
最近在虚拟机上进行数据解算的时候，由于误操作，把/usr目录下的权限都改成777导致无法运行sudo -s以及一切带sudo的命令了，惨痛的教训再一次提醒自己和大家，作为一个linux小白，千万不要瞎改权限。报错的截屏刚才忘了保存，具体的就是：我在运行sudo -s 周后终端直接报错：sudo: /usr/bin/sudo must be owned by uid 0 and have the setuid bit set带着问题一顿百度发现还挺多人出现这样的失误，于是按照各方大神的博客，解决了这个问题。
2020-09-29 13:17:50
1772
1
原创 [QOCA学习笔记]利用mload模块计算地表质量负载
在季节性时间尺度上，地表质量负载引起的季节性场地位移是空间大地测量中站点坐标时间序列的主要特征。关于QOCA软件的mload简介可以移步下载《Briefing of mass load caused site displacements》在处理数据时想按照文献中这样的方法计算安徽省区域的地表质量负载：（参考文献：武曙光. 区域CORS站坐标时间序列特征分析[D]. 2017.）首先下载文中说到的需要的数据。其中，大气压数据、土壤湿度和积雪深度数据下载机构是NCEP，均采用NCEP-DOE Reanal
2020-09-26 10:52:23
809
4
原创 [Python画图笔记]利用Python画箱型图boxplot
最近在学习使用Python画图，想用subplot画两幅箱型图，分别用来表示三组数据在空间滤波前后WRMS值的变化情况。关于Pycharm的安装和配置可以自行百度，直接贴上代码，给出效果图。代码中给出标注以便日后查询。import pandas as pd #导入pandasimport matplotlib.pyplot as pltimport numpy as npfrom scipy import interpolate#fig,axes = plt.su
2020-08-27 19:59:05
5070
1
原创 [QOCA学习笔记]利用QOCA软件analyze_tseri模块进行时间序列分析
操作环境：win10+虚拟机ubuntu下QOCA（ver 1.33）数据：通过解算得到的安徽省CORS参考站坐标时间序列（neu格式）对于QOCA软件的介绍，在其官网https://qoca.jpl.nasa.gov/上可以查看相关的描述。其analyze_tseri模块可以用来分析GNSS坐标时间序列。对该模块的详细介绍可以在这里下载：analyze_tseri analysis.pdf与之前pca模块相似，在调用analyze_tseri模块之前亦需要仔细查阅其驱动文件的设置。下面简单介绍一次
2020-08-05 21:05:22
1437
2
原创 [QOCA学习笔记]利用QOCA软件PCA模块进行共模误差分析
操作环境：win10+虚拟机ubuntu下QOCA（ver 1.33）数据：通过解算得到的安徽省CORS参考站坐标时间序列（neu格式）QOCA软件是优秀的GNSS坐标时间序列分析软件，其简介及PCA模块的使用方法在QOCA online tutor上已经有详细介绍：下面简单介绍利用PCA模块对时间序列的共模误差进行分析的过程（具体PCA的介绍请下载相关文献）...
2020-08-03 11:22:16
1438
原创 Ubuntu 16.04 + gmt 5.4 画某省地图（以安徽省为例）
操作环境：win10虚拟机上64位的Ubuntu16.04；GMT版本为5.4GMT操作指导手册可以在这里下载：GMT中文社区想在论文中画出某一区域的地图，在之前一篇博文中参考GMT中文社区绘制区域CORS站的站点分布和速度场图（脚本）已经记录过，现在想画出更好看的图，且达到以下要求：带上地形起伏数据带省级区域划分只显示目标区域内的地形，以便突出研究区域内的地形变化类似这样的图：首先要明确要达到以上目标要准备的画图数据：1.地形起伏数据：http://mirrors.ustc.edu.c
2020-07-16 14:11:01
3658
4
原创 空间直角坐标转换成大地坐标matlab代码（附精度）
在实际生产中，通常需要进行空间直角坐标和大地坐标之间的相互转换。精度估计是大地测量数据处理的一项内容，而目前使用的的坐标转换软件中一般都没有提供空间直角坐标与大地坐标之间转换的精度。其中，空间直角坐标至大地坐标的转换常用的迭代法可以表示为：在文章《空间直角坐标与大地坐标的精度关系式_徐思达》中，作者论通过误差传播律推导了由空间直角坐标转到大地坐标的精度，分析了严密公式和简便公式的精度转换结果的差异，发现两者只有微小差异，所以可以用简便公式。根据其提供的公式，用matlab编写程序实现。function
2020-05-29 20:54:46
6035
2
原创 QOCA软件安装（附安装后出现缺少libg2c.so.0报错的解决办法）
QOCA（Quasi-observation Combination Analysis），由加州理工学院喷气推进实验室设计和开发的。其是将各种松散约束的大地测量站点坐标和速度解(作为准观测)相结合，获得地壳形变信息的软件包。QOCA对任何非商业用户都是免费的。许多大地测量数据分析研究人员用QOCA作后处理软件包。目前，它可以结合空间大地准观测(GPS, VLBI, SLR，…等)和大地测量准观测(EDM，三角测量，水准测量，…等等)。它有可能将SAR数据与重力、地震活动和地面运动数据结合起来。QOCA软件
2020-05-28 14:38:45
2707
5
原创 小木虫网站无法登陆（账号突然被封禁）的解决办法
去年下半年开始很多小木虫账号都被封禁了，而在用户这边会发现，首次登陆该论坛时，网页会一片空白，之后打开任何一个小木虫的网站都是空白，或者是显示找不到网页的404错误，。偶然一次在百度贴吧上看到才知道原来是账号被封了，贴吧里也有人提出了解决的办法：首先要借个账号，找同学帮忙申诉，结果他的号是新的还没办法发帖子…就很迷惑。后来看到贴吧里有人说可以给小木虫发邮件申诉账号，邮箱是：emuch2018@163.com经过实测这个方法效率挺高，我前一天晚上发的邮件第二天中午就回复了，只是理由有点无语。于是乎我
2020-05-26 13:45:17
8358
1
原创 [GAMIT学习笔记]报错：Error reading grid file ERROR 5002
在运用gamit解算时，有时会发生如下错误：其中错误：FATAL : 200426:1158:34.0 GRDTAB/get_otl_grid: Error reading grid file ERROR 5002可能是有以下两个原因导致的：工程目录/tables下的otl.grid的链接是无效链接。没有下载正确的海潮改正文件。正确的解决办法是，首先去官网上下载有效地海潮改正格...
2020-04-29 11:51:31
1325
原创 Median Absolute Deviation(MAD，绝对中位差)和3sigma准则剔除粗差方法对比（代码）
根据文献：Klos A, Bogusz J, Figurski M, Kosek W (2016) On the handling of outliers in the GNSS time series by means of the noise and probability analysis. In: Rizos C, Willis P (eds) IAG 150 Years. Spring...
2020-01-27 16:52:08
4457
2
原创 EMD端点效应MATLAB代码注释（G. Rilling版本）
以下是G. Rilling版本的EMD代码中，抑制端点效应的镜像延拓MATLAB代码注释：%BOUNDARY_CONDITIONS_EMD extends an extrema set to limit edge effects on the interpolations%延拓极值点% % [TMIN,TMAX,ZMIN,ZMAX] = BOUNDARY_CONDITIONS_EMD(I...
2019-12-31 17:21:48
2860
6
原创 参考GMT中文社区绘制区域CORS站的站点分布和速度场图（脚本）
软件：GMT5.4 操作环境：win10虚拟机下的ubuntu16.04系统在学习绘制区域CORS站的站点图以及速度场时，参考GMT中文社区里贴出来的一些脚本文件，结合自己的数据进行了一些修改。1.绘制站点分布图参考脚本网址（内含各省会的经纬度）：https://gmt-china.org/example/ex004/在绘制自己所需的区域站点分布图时，租要自己制作该区域内各个站点的经纬度...
2019-10-18 21:19:46
5703
2
原创 [GAMIT学习笔记]一些常用到的脚本合集(2020.0926更新)
操作环境：ubuntu16.04+gamit10.61学习GAMIT学习过程中，在Linux系统下经常用到的一些命令，利用脚本可以免去很多麻烦，现整理如下：解压*.Z文件：uncompress *.Z生成l文件：在 rinex文件夹下依次运行：grep POSITION *.19o > lfile.rnxrx2apr lfile.rnx 2019 216gapr_to_l l...
2019-09-11 13:50:42
1688
9
原创 [CATS软件学习笔记]CATS分析多个测站多种噪声模型组合的脚本
CATS（Create and Analyse Time Series）：cats is a program to use Maximum Likelihood Estimation to fit a multi-parameter model to a time series (such as continuous GPS). The routine solves for all paramet...
2019-09-06 21:05:33
1552
15
原创 CATS软件安装[2022.1.20更新]
官方网址：http://web.mit.edu/mfloyd/www/computing/cats/教学视频：https://ke.qq.com/course/2830241.Linux下解压文件与压缩文件：（1）压缩文件。命令格式： tar -zcvf 压缩文件名：.tar.gz（2）解压文件。命令格式： tar -zxvf 压缩文件名：.tar.gz 压缩后的文件只能放在当...
2019-09-04 10:28:59
4511
19
原创 [GAMIT/GLOBK学习笔记]制作station.info文件的方法（脚本）
在使用GAMIT/GLOBK软件解算GNSS数据时，首先需要运行GAMIT软件进行单日解解算，再将单日解的结算结果h文件作为GLOBK输入文件进行进一步解算，得到测站坐标、时间序列或者是速度场。在使用GAMIT软件时，无论是分部解算（刚学习时采用分布处理有助于了解软件的解算步骤，或者是在某天解算出问题的情况下通过分布解算来查找和排查错误）还是批处理（比较方便，适合时间跨度长的情况下），其间一个重要...
2019-08-14 17:13:00
2813
9
原创 ubuntu系统下gnuplot画时间序列图
在ubuntu系统下画时间序列图，源数据是QOCA的输入数据，格式如下：画出来的时间序列图如下：下面介绍如何用gnuplot软件的multiplot画出上述图。首先，在含有源数据的文件夹下 打开终端输入：出现上述界面即进入gnuplot绘图，再输入set multiplot layout 3,1表示我们要在一幅图上画三个图，排列是3*1：此时在终端输入后会弹出一块画布，现在还什么...
2019-07-17 10:02:46
939
原创 [QOCA学习笔记]analyze_tseri报错：get_st_lst :missing a station. nfnd nsit=1 22
第一次使用QOCA软件分析时间序列，查找原因是因为结算前没有提前准备先验坐标文件apr文件。用gamit/globk的apr文件，通过mk_net将其转换成QOCA需要用到的apr先验文件即可：在gamit解文件夹/rinex/下cat *.apr > globk.aprmk_net globk.apr qoca.apr globk即可生成名为qoca.apr的qoca需要用的先...
2019-07-16 20:45:46
490
原创 ubuntu报错：libg2c.so.0:cannot open shared object file: no such file or directory/wrong ELF class（已解决）
安装完QOCA软件运行analyze_tseri模块时遇到问题：百度过后发现可能是两个原因（参见：https://blog.csdn.net/wujiaqi168/article/details/17055543）经过检查发现在我的系统里没有libg2c.so.0这个库，所以要按照https://forum.ubuntu.org.cn/viewtopic.php?t=474325这里的方法进...
2019-07-16 13:32:42
2018
1
analyze_tseri analysis.pdf
2020-08-03
Second generation of QOCA st_filter .pdf
2020-08-03
Principal component analysis.pdf
2020-08-02
QOCA指南.zip
2019-07-16
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人