matlab读取gamit解算的met文件,获取ZTD、ZWD、PWV、ZHD等数据(附matlab代码)

2024-01-13 08:50

本文主要是介绍matlab读取gamit解算的met文件,获取ZTD、ZWD、PWV、ZHD等数据(附matlab代码),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

0. 码字不易,点赞加关注(公众号:WZZHHH,后期代码将在公众号可以下载),使用请注明出处(根据我的研究方向,我会不断更新代码)。

1. 读取gamit解算的met文件,根据前缀met_读取相应文件下met文件,met文件的列分别为1为year、2为doy(年积日)、3为小时、4为分钟、5为秒、6为ZTD、7为ZWD、8不知道、9为PWV、11为气压、12为温度、13为ZHD。

2. 下图为具体的met文件,对照上面列,找到自己需要的。

3. 我按照年积日存放的met文件,存放方式不同的,记得修改代码。我的存放方式如下图。

4. 下面是我的代码,可能没大佬写得好,凑合看吧,能用。

我用的ZTD数据是第六列,其他数据可在第84行更改列即可。

function load_ZTD_file
% 读取gamit解算的met文件
% 我按照年积日存放的met文件
% 存放方式不同的,记得修改
% 根据前缀met—读取相应文件下ztd文件清单,再根据需要更改列
% 我用的ZTD数据是第六列。其他数据可在第84行更改
% met文件的列分别为1为year、2为doy(年积日)、3为小时、4为分钟、5为秒、
% 6为ZTD、7为ZWD、8不知道、9为PWV、11为气压、12为温度、13为ZHD。clc,clear;
% met文件路径,2021年数据为例
ypath = 'D:\paper_write\ZTD\2021\';  
% 提取数据后,存放路径
save_path= 'D:\paper_write\paper_code\2\mat\ztd\2021ztddata.mat';
%% ---------------------找到存在met的站点名字---------------------------% 遍历文件夹,看总共有几个站点
ydata = dir(ypath);station = {};  % station 站点
for i=3:size(ydata,1)% 遍历年积日文件夹dpath = [ypath ydata(i).name '\'];ddata = dir(dpath);for ii=3:size(ddata,1)% i=3是第一天,第一天初步建立站点名称if i==3station = [station;{ddata(ii).name(5:8)}];else% 后面增加的站点加入到station中if ~contains(ddata(ii).name(5:8),station)station = [station;{ddata(ii).name(5:8)}];endendendend%% -----------------------------提取gamit解算的met文件------------------% 这个部分本可以和上面合在一起,为了便于理解,,所以分开写ztd_all = {};
parfor i = 3:size(ydata,1)dpath = [ypath ydata(i).name '\'];% 避免文件夹有其他文件,只读取met_开头List = dir(fullfile(dpath,'/','met_*'));for I = 1:size(List,1)% 站点名称stn = List(I).name(5:8);% 读取met文件filename = [dpath,List(I).name];fid = fopen(filename,'r');% 提取met文件的ztdztd = [];% 逐行读取while ~feof(fid)line = fgetl(fid);% 依次读取每行,met文件带*的删除,剩下就是想要的数据if ~contains(line,'*')ztd = [ztd;str2num(line)]; %#ok<ST2NM>endend% 我使用的是1小时时间分辨率% 避免个别时间可能会缺失,缺少数据为nanztd = hours(ztd);fclose(fid);% cellstr:char转cell、num2cell:double转cellstn = repmat(cellstr(stn),24,1);% 第六列是ztd,我只需要这个ztd_all = [ztd_all;[stn num2cell([ztd(:,1:3) ztd(:,6)])]];end
end
% ztd_2:把ztd相同测站放在一起,相同测站按照日期排列更好计算
% ztd_2就是排列好的ztd_all% 按照station名称排列站点数据
ztd = [];
for I =1:size(station)flag = contains(ztd_all(:,1),(station(I,:)));ztd_1 = ztd_all(flag,:);ztd = [ztd;ztd_1];
endsave(save_path,'ztd');
endfunction [ztd] = hours(ztd)
% 25行证明数据齐全,没有的为nan
if size(ztd,1)==25% 第25行为下一天的0时,对其删除ztd = ztd(1:24,:);
else% 年月日时,后面为0的数组ztdz = nan(24,14);yearz = repmat(ztd(1,1:2),24,1);hour = (0:1:23).';ztdz = [yearz hour ztdz];% 将存在的数据进行填充% ztd第3列为小时0-23,把存在数据的小时填充到ztdz中for i=1:size(ztd,1)ztdz(ztd(i,3)+1,:) = ztd(i,:);endztd = ztdz;
endend

这篇关于matlab读取gamit解算的met文件,获取ZTD、ZWD、PWV、ZHD等数据(附matlab代码)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/600894

相关文章

Python Transformers库(NLP处理库)案例代码讲解

《PythonTransformers库(NLP处理库)案例代码讲解》本文介绍transformers库的全面讲解,包含基础知识、高级用法、案例代码及学习路径,内容经过组织,适合不同阶段的学习者,对... 目录一、基础知识1. Transformers 库简介2. 安装与环境配置3. 快速上手示例二、核心模

Spring Boot读取配置文件的五种方式小结

《SpringBoot读取配置文件的五种方式小结》SpringBoot提供了灵活多样的方式来读取配置文件,这篇文章为大家介绍了5种常见的读取方式,文中的示例代码简洁易懂,大家可以根据自己的需要进... 目录1. 配置文件位置与加载顺序2. 读取配置文件的方式汇总方式一:使用 @Value 注解读取配置方式二

Spring 请求之传递 JSON 数据的操作方法

《Spring请求之传递JSON数据的操作方法》JSON就是一种数据格式,有自己的格式和语法,使用文本表示一个对象或数组的信息,因此JSON本质是字符串,主要负责在不同的语言中数据传递和交换,这... 目录jsON 概念JSON 语法JSON 的语法JSON 的两种结构JSON 字符串和 Java 对象互转

Java的栈与队列实现代码解析

《Java的栈与队列实现代码解析》栈是常见的线性数据结构,栈的特点是以先进后出的形式,后进先出,先进后出,分为栈底和栈顶,栈应用于内存的分配,表达式求值,存储临时的数据和方法的调用等,本文给大家介绍J... 目录栈的概念(Stack)栈的实现代码队列(Queue)模拟实现队列(双链表实现)循环队列(循环数组

C++如何通过Qt反射机制实现数据类序列化

《C++如何通过Qt反射机制实现数据类序列化》在C++工程中经常需要使用数据类,并对数据类进行存储、打印、调试等操作,所以本文就来聊聊C++如何通过Qt反射机制实现数据类序列化吧... 目录设计预期设计思路代码实现使用方法在 C++ 工程中经常需要使用数据类,并对数据类进行存储、打印、调试等操作。由于数据类

C# foreach 循环中获取索引的实现方式

《C#foreach循环中获取索引的实现方式》:本文主要介绍C#foreach循环中获取索引的实现方式,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考下吧... 目录一、手动维护索引变量二、LINQ Select + 元组解构三、扩展方法封装索引四、使用 for 循环替代

SpringBoot使用GZIP压缩反回数据问题

《SpringBoot使用GZIP压缩反回数据问题》:本文主要介绍SpringBoot使用GZIP压缩反回数据问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录SpringBoot使用GZIP压缩反回数据1、初识gzip2、gzip是什么,可以干什么?3、Spr

Linux下如何使用C++获取硬件信息

《Linux下如何使用C++获取硬件信息》这篇文章主要为大家详细介绍了如何使用C++实现获取CPU,主板,磁盘,BIOS信息等硬件信息,文中的示例代码讲解详细,感兴趣的小伙伴可以了解下... 目录方法获取CPU信息:读取"/proc/cpuinfo"文件获取磁盘信息:读取"/proc/diskstats"文

基于Python实现读取嵌套压缩包下文件的方法

《基于Python实现读取嵌套压缩包下文件的方法》工作中遇到的问题,需要用Python实现嵌套压缩包下文件读取,本文给大家介绍了详细的解决方法,并有相关的代码示例供大家参考,需要的朋友可以参考下... 目录思路完整代码代码优化思路打开外层zip压缩包并遍历文件:使用with zipfile.ZipFil

Vue3组件中getCurrentInstance()获取App实例,但是返回null的解决方案

《Vue3组件中getCurrentInstance()获取App实例,但是返回null的解决方案》:本文主要介绍Vue3组件中getCurrentInstance()获取App实例,但是返回nu... 目录vue3组件中getCurrentInstajavascriptnce()获取App实例,但是返回n