GNSS
achumoyangguang
GNSS气象学
展开
-
提取VMF1中的Tm
clc;clear;name=['y2018','.vmf1_g'];fid = fopen(name);a=[];a1=[];A1=[];if fid>0 i=0; while feof(fid) == 0 i = i+1; tline{i} = fgetl(fid); a=tline(i); if i>22 a1=cell2mat(a); if (a1(1.原创 2020-05-09 16:08:24 · 237 阅读 · 1 评论 -
计算VMF1和探空Tm的精度
clc;clear;name1=['chan','.txt'];[site,mjd,ah,aw,zhd,zwd,tm]=textread(name1,'%s%s%s%s%s%s%s%*[^\n]','delimiter',' ');%提取前5列数据n=length(tm);for i=1:n%在1-n行数据内判断数据是否为空,如果为空则跳出本次循环 Tm(i,:)=str2num(tm{i}); s=i/4+1; j=floor(s); if mod(i,4)==.原创 2020-05-09 16:08:57 · 353 阅读 · 0 评论 -
批量提取met文件
find ./ -name "met_*" | xargs -i cp {} ./met-site/原创 2020-09-25 19:37:52 · 743 阅读 · 0 评论 -
提取探空站pwv并存储
% ①数据提取pwv并存储clear;clc;year=2018;month=10;b=[];PWV=[];pwv=[];month_day=7;for day=1:month_day for time=[00,12];%以下将年月日一天的次数数字转为字符串 years=num2str(year); months=num2str(month); days=num2str(day); times=num2str(time).原创 2020-10-05 15:11:37 · 1091 阅读 · 10 评论 -
查找空rinex
% 查找是否存在空文件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:19:27 · 80 阅读 · 0 评论 -
提取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 · 719 阅读 · 0 评论 -
对流层延迟产品网站
igs产品:https://cddis.nasa.gov/archive/gnss/products/troposphere/zpd/ftp://cddis.gsfc.nasa.gov/gnss/products/troposphere/zpd/(Anonymous ftp service will be discontinued on October 31, 2020.)ftp://gssc.esa.int/gnss/products/troposphere_zpd/%以下内容取自igs官原创 2020-10-05 20:43:43 · 1854 阅读 · 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 · 772 阅读 · 0 评论