登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 lynpatty 于 2018-6-25 19:41 编辑
想向各位老师请教,在计算某区域全年降水百分位值时需要先把missing value 设置为NaN后进行计算,代码如下。但我查了一下最后得出来的百分位值要比实际低很多。应该是在处理missing value这一块时在数据量级上出了问题,想请各位老师看看代码哪里出了错。
clear all; clc
fname='/User/Downloads/cpc-precip/precip.1982.nc'; ncid = netcdf.open(fname,'NOWRITE'); %?¨°??nc???? ncdisp(fname); %% mm/day precip=ncread(fname,'precip'); lon=ncread(fname,'lon'); lat=ncread(fname,'lat'); time=ncread(fname,'time'); netcdf.close(ncid);
data=precip(472:491,105:114,1:365); data(data==-9.969209968386869e+36)=NaN; p1=prctile(data(:),95);
在工作区出现的data值都偏小,有的甚至为零。。。 0 NaN NaN NaN NaN
0 0 0 NaN NaN
0 0 0 0 NaN
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0.2710 0 0 0 0
0.5730 0.0262 0 0 0
0.0277 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
|