GAMIT
achumoyangguang
GNSS气象学
展开
-
GAMIT基线解算_批处理
分步处理2017年001,002天(shao hkkt hkks,后面两个站点为香港CORS站)数据准备工作:1)mkdir demo2)cd demo3)mkdir rinex brdc igs4)cd brdc采用wget下载wget -c ftp://nfs.kasi.re.kr/gps/data/daily/2017/001/17n/brdc001...原创 2020-08-23 09:10:31 · 2114 阅读 · 2 评论 -
GAMIT基线解算_分部处理
分步处理2017年001天(chan shao bifs)数据准备工作:1)mkdir demo2)cd demo3)mkdir rinex brdc igs4)cd brdc两种下载方式:第一种方式:采用gamit自带的下载脚本sh_get_nav -archive sopac -yr 2017 -doy 001 -ndays 1第二种方式:采用...原创 2020-08-23 09:11:06 · 1816 阅读 · 1 评论 -
GAMIT10.61_处理BDS数据
数据下载广播星历:ftp://ftp.cddis.eosdis.nasa.gov/pub/gps/data/campaign/mgex/daily/rinex3/2018/001/18p/ (chmod +w brdm0010.18p)精密星历:http://www.igs.gnsswhu.cn/index.php/Home/DataProduct/mgex.html(文件只读改为可读写:...原创 2020-08-23 09:11:25 · 638 阅读 · 1 评论 -
批量提取met文件
find ./ -name "met_*" | xargs -i cp {} ./met-site/原创 2020-09-25 19:37:52 · 803 阅读 · 0 评论 -
在gamit中对rinex3的转换
下载的观测源文件:HKKT00HKG_R_20180010000_01D_30S_MO.crx.gz在gamit中采用如下方式1/解压:gunzip HKKT00HKG_R_20180010000_01D_30S_MO.crx.gz结果:HKKT00HKG_R_20180010000_01D_30S_MO.crx2/重命名:把.crx命名成.dsh_rename_rinex3 HKKT00HKG_R_20180010000_01D_30S_MO.crx结果:hkkt0010.18原创 2020-10-04 11:55:02 · 1215 阅读 · 1 评论 -
gamit10.72中对流层延迟解算参数
最大值只能为100原创 2020-12-01 11:20:58 · 1317 阅读 · 7 评论 -
查找空met
% 查找是否存在空文件clear;clc;Empty=[];site = ['bjfs';'alic';'darw';'daej';'hkoh';'hknp';'hkst';'hkpc';'hksc';'hksl';'hkkt';'hkss';'hklt';'hkws';'t430'];for j=1:15 for year=19 yr = num2str(year); for i=1:9 doy = num2str(i); .原创 2020-10-12 22:20:05 · 91 阅读 · 0 评论 -
更新gamit的tables表问题
网站下载说明:1/http://garner.ucsd.edu/pub/gamit/tables/(格网表文件不全)采用迅雷下载不了其他文件除了VMF文件但是可以采用filezilla下载,在虚拟机上和win10上均可以,亲测可用2/http://sopac.ucsd.edu/GAMIT.shtml直接点击目标文件即可下载但是发现有些文件不全,如Solar, Lunar, Nutation仅到2016年的3/ftp://everest.mit.edu//pub/GRIDS/(格网原创 2020-10-15 21:41:09 · 1898 阅读 · 1 评论 -
提取igs产品zpd
%提取zpdclear all;clc;day='缺失天数'; %定义变量名day,记录缺失的文件年积日ZPD = [];for n=1:31 %zpd文件天数 ns=num2str(n,'%03d'); % 输值转字符串,固定三位数字(001,010,100) name=['hksl',ns,'0.19zpd']; % 字符串拼接得到文件名称 ffid=fopen(na.原创 2020-10-22 22:13:16 · 739 阅读 · 0 评论 -
使用sh_metutil生成采样300秒的ztd
#!/bin/bashfor ((i=1;i<=365;i=i+1))dodoy=$(printf "%03d" "$i")cd ${doy}sh_metutil -f o*a.${doy} -z z*.${doy} -i 300cd ..done原创 2020-12-24 17:05:27 · 298 阅读 · 2 评论 -
gamit10.72处理BDS数据
1/下载数据p文件ftp://igs.gnsswhu.cn/pub/gnss/mgex/daily/rinex3/2020/001/20p/o文件ftp://ftp.geodetic.gov.hk/rinex3/2020/001/hkpc/30s/ftp://gssc.esa.int/gnss/data/daily/2020/001/sp3文件ftp://igs.gnsswhu.cn/pub/gnss/products/mgex/2086/ftp://gssc.esa.int/gnss/原创 2021-01-17 12:06:38 · 864 阅读 · 0 评论 -
下载电离层产品
#!/bin/bashyear=2019sitename=${site[j]}for ((i=1;i<=365;i=i+1))dodoy=$(printf "%03d" "$i")yr=${year:2:4}filename=igsg${doy}0.${yr}i.Zwget --auth-no-challenge https://cddis.nasa.gov/archive/gnss/products/ionex/${year}/${doy}/${filename}gunzip $.原创 2021-04-21 19:36:41 · 824 阅读 · 0 评论