爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 6666|回复: 3

[程序设计] nc求降水阈值出错

[复制链接]
发表于 2018-6-11 06:28:13 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
想要利用nc文件算出某年降水95%的百分位值,用的prctile函数,
程序如下:fname = '/Users//Downloads/cpc-precip/precip.1979.nc';
ncid = netcdf.open(fname,'NOWRITE');
ncdisp(fname);
precip = ncread(fname,'precip');
lon = ncread(fname,'lon');
lat = ncread(fname,'lat');
time = ncread(fname,'time');
netcdf.close(ncid);
v = prctile(precip,95)

但得出来的结果不是0就是负数,想问一下这个是什么问题呢?
1.0e+36 *

  1 至 10 列

   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692

  11 至 20 列

   -9.9692   -9.9692   -9.9692         0         0         0         0         0         0         0

  21 至 30 列

         0         0    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  31 至 40 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  41 至 50 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  51 至 60 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  61 至 70 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  71 至 80 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  81 至 90 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  91 至 100 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  101 至 110 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  111 至 120 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  121 至 130 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  131 至 140 列

    0.0000    0.0000    0.0000    0.0000         0         0         0         0         0         0

  141 至 150 列

         0         0         0         0         0         0         0         0         0         0

  151 至 160 列

         0         0         0         0         0         0         0         0         0    0.0000

  161 至 170 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  171 至 180 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  181 至 190 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  191 至 200 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  201 至 210 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  211 至 220 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  221 至 230 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  231 至 240 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

  241 至 250 列

    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.0000    0.0000    0.0000

  251 至 260 列

         0         0         0         0         0         0         0         0   -9.9692   -4.9846

  261 至 270 列

   -9.9692         0         0   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692

  271 至 280 列

   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692

  281 至 290 列

   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692

  291 至 300 列

   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692

  301 至 310 列

   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692   -9.9692

  311 至 320 列

   -9.9692   -9.9692         0         0         0         0         0         0         0         0

  321 至 330 列

         0         0         0         0         0         0         0         0         0         0

  331 至 340 列

         0         0         0         0         0         0         0         0         0         0

  341 至 350 列

         0         0         0         0         0         0         0         0         0         0

  351 至 360 列

         0         0         0         0         0         0         0         0         0         0
密码修改失败请联系微信:mofangbao
发表于 2018-6-11 08:16:14 | 显示全部楼层
应该是默认填充数据的事情。
楼主可以看看量纲,很明显是不对的。
建议将默认填充数据改为nan再进行排序。
密码修改失败请联系微信:mofangbao
发表于 2019-3-14 09:28:17 | 显示全部楼层
请问楼主解决了吗?
密码修改失败请联系微信:mofangbao
发表于 2022-4-5 21:16:14 | 显示全部楼层
请问解决这个问题了吗

密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表