chaoshu的个人博客分享 http://blog.sciencenet.cn/u/chaoshu

博文

matlab读取GRACE_Mascon的NC格式文件

已有 1149 次阅读 2023-12-21 18:04 |个人分类:水文大地测量|系统分类:科研笔记

1、CSR官网下载NC文件:GRACE/GRACE-FO - Gravity Recovery and Climate Experiment (utexas.edu)

GLDAS数据也是同样的处理方式;

2、matlab读取NC内容

2.1 matlab将NC文件中的gregorian日期转换成年月日ymd

dt = datetime((time)*24*3600,'ConvertFrom','epochtime','Epoch','2002-01-01');

NC文件中时间time表示从2002年1月1日至目前所要读取文件时间为止累计的天数

% 提取给定时间的索引

t0=datetime(2020,06,15,00,00,00);

[~,ind1] = min(abs(datenum(dt)-datenum(t0)));

closest_time = dt(ind1);

2.2 NC读取数据需要翻转90°

rot90(csr_mascon(:,:,idx));

(不知道nc为什么是经纬度翻转了)matlab直接绘图能看到它是反过来的,所以要转置下;

下面是翻转前后对比;

image.pngimage.png

2.3 获取指定经纬度范围数据

% 经纬度范围

lat_range = [20, 30];

lon_range = [100, 110];

% 查找经纬度范围在nc文件中的索引

indices_lat = find(lat >= lat_range(1) & lat <= lat_range(2));

indices_lon = find(lon >= lon_range(1) & lon <= lon_range(2));

% 提取指定范围内的数据

data_subset = data(indices_lon, indices_lat, :);

3、完整代码截图

image.png



https://blog.sciencenet.cn/blog-858128-1414820.html

上一篇:地大测绘科学与技术和地球物理学期刊分类
下一篇:如何获得水文站经纬度
收藏 IP: 1.189.61.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-5-20 01:33

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部