关于如何读取某边界范围以及某段时间的NC数据的问题。

tonghua111

请问楼主，你读取以黑河流域为边界提取数据时，按照你所编写的程序，不是提取出来一个方形区域的数据吗？就是说提取的不是完全符合黑河流域边界的数据，只是提取的数据包含黑河流域。

DoneThat

学习一下                             

wangdapang

子达 发表于 2016-4-6 18:41
air_temperature= ncread('tas_Amon_bcc-csm1-1_historical_r1i1p1_185001-201212.nc','tas',[1,1,1321],[1 ...

在吗，我想知道你是怎么确定起始月份的位置的，我现在就遇到这个问题啦

ychun

谢谢，学习到了

zqzq

lon=ncread('C:\Users\mmm\Desktop\data.nc','X');
变量里自己就含带了时间信息
可以参考
lat=ncread('C:\Users\mmm\Desktop\data.nc','Y');
i_inx=find(lon>=100&lon<=180);
j_inx=find(lat>=0&lat<=40);
nx=length(i_inx);
ny=length(j_inx);
ssta=ncread('C:\Users\mmm\Desktop\data.nc','ssta',[i_inx(1) j_inx(1) 1],[nx ny 456]);

zqzq

lon=ncread('G:\DHW\ct5km_dhw-max_v3.1_201904.nc','lon');
lat=ncread('G:\DHW\ct5km_dhw-max_v3.1_201904.nc','lat');
i_inx=find(lon>=104&lon<=122);
j_inx=find(lat>=2&lat<=24);
nx=length(i_inx);
ny=length(j_inx);
HS(:,:,s)  = ncread(filename,'hotspot',[i_inx(1) j_inx(1) 1],[nx ny 1]);; %读入变量lat
这样应该也是可以的

kaiofchina

学习了，谢谢分享

zhx1111

tonghua111 发表于 2017-8-20 08:10
请问楼主，你读取以黑河流域为边界提取数据时，按照你所编写的程序，不是提取出来一个方形区域的数据吗？就 ...

请问您理解了吗，我刚做这一块，也有这个疑问

zhx1111

zhx1111 发表于 2020-8-4 16:57
请问您理解了吗，我刚做这一块，也有这个疑问

已解决，我用的是arcgis中的剪切功能

tonghua111

可以用inpolygone函数判断在边界内的数据，然后提取flag==1的值就OK了