matlab读取igs文件,Matlab 等間隔讀取IGS對流層天頂延遲（ZPD）文件

IGS(International GNSS Service)提供對流層天頂延遲(ZPD)文件，可從cddis分析中心下載(ftp://cddis.gsfc.nasa.gov/gps/products/trop_zpd)，其數據間隔為300s，使用MATLAB提取2小時間隔的對流層天頂延遲數據。

以DARW站2016年第1天至第366天數據為例 (注：ZPD文件自2017年第29天起記錄發生變化，使用該程序需做調整)

⑴數據提取

clear all;

clc;

day='缺失天數'; %定義變量名day，記錄缺失的文件年積日for n=1:366; %zpd文件天數

ns=num2str(n,'%03d'); % 輸值轉字符串，固定三位數字(001,010,100)

name=['darw',ns,'0.16zpd']; % 字符串拼接得到文件名稱

ffid=fopen(name,'rt'); % 打開zpd文件，從第一行開始讀取，數據存儲在cell里

if ffid>0 % 文件存在，讀取數據i = 1;

while feof(ffid) == 0

tline{i,1} = fgetl(ffid);

i = i+1;

end

% 提取間隔2h的數據，從0時至22時，zpd文件采樣間隔為300s，2h間隔為24行

for i=62:24:326

a((i-38)/24,:)=char(tline{i,1}); %從cell里提取數據轉成字符型

end

for i=1:12

zpd(i,:)=str2num(a(i,20:25)); %從字符串里提取第20至第25個字符並轉成數值

end

xlswrite(ns,zpd) %按天寫入excel表

else

day=strcat(day,'-',ns); %文件不存在，記錄缺失文件年積日

xlswrite(ns,nan); %生成空文件

end

⑵數據合並

clear;

clc;

for n=1:366

n=num2str(n,'%03d');

eval(['a',n,'=','xlsread(n)']);

end

n=2; ss=a001;

while(n<367)

n=num2str(n,'%03d');

ss=cat(1,ss,eval(['a',n]));

n=str2num(n);

n=n+1;

end

xlswrite('2016.xlsx',ss)

matlab读取igs文件,Matlab 等間隔讀取IGS對流層天頂延遲（ZPD）文件相关推荐

elsx文件转为csv并隔三取一
elsx文件转为csv并隔三取一 import csv import pandas as pd#转换xlsx文件为csv文件 def xlsx_to_csv_pd():data_xls = pd.re ...
matlab读取hdf显示,matlab读取.hdf文件
matlab读取.hdf文件打开hdf文件,以臭氧廓线数据为例 file = 'OMI-Aura_L2-OMO3PR_2009m0107t1808-o23846_v003-2009m1025t065 ...
matlab读取grib2数据,matlab读取grib2数据
2018 基于通用模型的 GRIB 格式数据读取技术王兵,李杰 (1. 南京航空航天大学民航学院,江苏南京 211106;2. 国家空管飞行流量技术重点实验室,江苏南京 211106) 摘要: ...
matlab读取二进制文件字符串,matlab读取内容为二进制的TXT文件
本方法同样适合读取十六进制和二进制以外的其他进制文件, txt使用一个最简单的命令就可以读取 textread 这是一个十分有用,简便的函数(对于fopen fscanf而言) 读取二进制txt文件: ...
matlab 读取通达信,Matlab 自动导入通达信板块文件
http://blog.sina.com.cn/s/blog_d6c58efd0102x0gx.html 通达信具有强大的板块处理能力,并且一些板块数据是动态更新的,做股票量化设计的时候,借助通达信的 ...
matlab读取nastran模型,Matlab 读Nastran结果文件.pch，并提取刚度阵和质量阵问题
我在前几天在CSDN上下载了一个提取Nastran结果.pch文件中的刚度阵和质量阵的m文件,代码如下: // function [K,M]=Get_K_M(N) %读取Nastran的pch文件,得 ...
matlab读取时间数据,Matlab有关数据库读取及时间项的设定
Matlab有关数据库的读取需下载mysql-connector-java-5.1.26-bin.jar文件放到指定位置,并修改classpath.txt相关文件,具体可网上自行查找. 看一下实例: ...
matlab读取图片亮度,Matlab中使用imagesc函数显示的亮度图像怎么和原图像的大小不一样？...
在matlab 2006上写的两个函数,这两个函数的调用格式为: 读DICOM图像数据: X = read_dicom(filename) 这个函数的作用是从filename所指定的文件中读取图像数据 ...
matlab读取trmm,基于Matlab的TRMM3B43数据处理的思维过程与技术流程
本文主要提供一套处理TRMM3B43数据的思维过程与技术流程,力求能够让读者在处理其他国际通用数据时也能够采用类似的方法来解决. 首先我们来看一下TRMM3B43数据的数据格式,是hdf格式的,相当于 ...

matlab读取igs文件,Matlab 等間隔讀取IGS對流層天頂延遲（ZPD）文件

matlab读取igs文件,Matlab 等間隔讀取IGS對流層天頂延遲（ZPD）文件相关推荐

最新文章

热门文章