- 积分
- 18
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-8-21
- 最后登录
- 1970-1-1
|
发表于 2017-2-6 17:30:07
|
显示全部楼层
应该是constant zero values 或者 constant 1 values 导致的
我实际测试过HadISST 的SIC,constant zero values 的格点还是有不少的,但1的基本没有
建议你 对原始数据的时间维求均值好了,基本上均值为零的都是所有值为零,然后可以在挑选数据的时候加上这个条件
比如 对 sst
ls = dim_avg(sst(latitude|:,longitude|:,time|:))
nsst = num(ls.ne.sst@_FillValue)
print("total numbers of ocean grids in the chosen area : " + nsst)
ls1d = ndtooned(ls)
lsind = ind(ls1d.ne.sst@_FillValue.and.ls1d.ne.0.and.ls1d.ne.1)
nsst = num(ls1d.ne.sst@_FillValue.and.ls1d.ne.0.and.ls1d.ne.1)
print("after removing constant 0/1 value," + nsst + " grids left" )
lsind就是后面提取去算SVD的关键指针了
不过我也有问题
NCL SVD函数,实际分析的时候应该用hetlft和hetrgt吧,homogeneous array 和 heterogeneous array应该怎么理解和翻译? 国内的书基本都说 左右 奇异向量,是分别对应hetlft和hetrgt? |
|