||
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直接绘图能看到它是反过来的,所以要转置下;
下面是翻转前后对比;
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、完整代码截图
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-1 09:21
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社