求助 ERA5数据怎么按照经纬度提取减小数据

ZOE

70多年的ERA5的NC数据，下载时候选了最大范围，现在MATLAB超过32G不能运行
想要把数据按照经纬度提取后数据变小，然后还是生成单个的NC数据 方便在电脑后续程序处理
单位不具备大型机计算能力 生娃后武功尽失 求支招求帮助一下

SunJiaming

%查看文件信息
finfo = ncinfo('E:\ecmwf_data\surface_pressure.nc');
%资料范围是0-60N, 70E-140E，水平分辨率0.1度
sp = ncread('E:\ecmwf_data\surface_pressure.nc', 'sp');

%%% sp大小701 x 601 x 8 ， 701是纬向数据个数，601是经向数据个数，8是层数
% 比如说要提取0-6N, 70E-77E的数据，就选择响应的纬向、经向数据，经纬度要和数据个数有个换算
% 法1，直接提取，经纬度对应的数据个数要换算一下

start = [1, 1, 1];
count = [60, 70, 1];
stride = [1, 1, 1];
sp_select1 = ncread('E:\ecmwf_data\surface_pressure.nc', 'sp',start, count, stride);

% 法2, 截取数组sp
sp_select2 = sp(1:60,1:70,1);

SunJiaming

https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-land-monthly-means?tab=form
以此为例，在网页下方Geographical Area ----> Sub-region extraction这里可以选取经纬度范围

ZOE

感谢回复 就是下载数据的时候没有选择经纬度才想着看看有啥办法处理下。
不想再重新下载数据了 ，依然感谢！

SunJiaming

不客气~~~~~~

伽蓝鸟

MATLAB可以只读取nc文件的部分数据啊，你挑你需要的经纬度范围读进来，就行了

五里

ncread里有选项可以选

ZOE

伽蓝鸟 发表于 2021-7-25 14:16
MATLAB可以只读取nc文件的部分数据啊，你挑你需要的经纬度范围读进来，就行了

麻烦能说的具体一点么 MATLAB才开始学 小白一枚

wjy_ecnu

查看ncread的帮助文档 里面可以读取你想要的部分的数据

ZOE

SunJiaming 发表于 2021-7-25 21:28
%查看文件信息
finfo = ncinfo('E:\ecmwf_data\surface_pressure.nc');
%资料范围是0-60N, 70E-140E，水 ...

去试试 非常感谢