想用matlab读CMIP5的NC文件

469436838

在CMIP5下载了一种模式下rcp降雨数据
看到一些人的帖子知道了怎么读取文件
接下来想做一下数据的筛选，想从全球的100年中挑出某一年感兴趣经纬度的数据，应该怎么用matlab实现ne ?
clc
clear
lon=ncread('pr_Amon_bcc-csm1-1_rcp45_r1i1p1_200601-209912.nc','lon');% 读取NC数据的经度数据
lat=ncread('pr_Amon_bcc-csm1-1_rcp45_r1i1p1_200601-209912.nc','lat');% 读取NC数据的纬度数据
pr=ncread('pr_Amon_bcc-csm1-1_rcp45_r1i1p1_200601-209912.nc','pr'); %读取NC数据的降水数据
time=ncread('pr_Amon_bcc-csm1-1_rcp45_r1i1p1_200601-209912.nc','time');% 读取NC数据中时间变量

接下来怎么提取年份和经纬度呢？请教大家。

sunzy025

你知道min函数吧~
用你想要的经纬度min这两个lon，lat[~ i]=min(abs(lat(:)-Want_lat));
[~ j]=min(abs(lon(:)-Want_lon));
然后就知道了位置（j,i）
之后你pr=pr(j,i,time,:);

469436838

sunzy025 发表于 2018-1-15 08:33
你知道min函数吧~
用你想要的经纬度min这两个lon，lat[~ i]=min(abs(lat(:)-Want_lat));
[~ j]=min(abs(l ...

你好
我用你这个方法编了一下，提示语法错误，能举个例子么？比如提取经度110-120，纬度30-40，应该怎么写？谢谢啦

sunzy025

469436838 发表于 2018-1-15 09:39
你好
我用你这个方法编了一下，提示语法错误，能举个例子么？比如提取经度110-120，纬度30-40，应该怎么 ...

刚才说的是提取单点的，其实你可以提取这个范围，然后 pr=pr（第一维范围,第二维范围,时间,:）

可能帮你编不了，因为没有安装matlab现在