提问:ncl做svd分析时出现"SVD:Warning all X values in column are missing

suifeng1106

DevilMaySmile 发表于 2016-5-17 10:38
因为我后面没有用svd算了，所以就没管这个问题了，不过好像楼上说虽然有warning 算出来的结果也差不多， ...

谢谢楼主的回答，这个问题我差不多弄清楚了。现在情况是我的SST场里陆地上的确有缺测值，然后算的时候ncl一直报warming，而且计算特别慢，我差不多运行了40分钟连一个量都没算好，但是如果数据没有缺测值又一下就算出来了。

DevilMaySmile

suifeng1106 发表于 2016-5-17 16:07
谢谢楼主的回答，这个问题我差不多弄清楚了。现在情况是我的SST场里陆地上的确有缺测值，然后算的时候ncl ...

哦，是这样，弄清楚了就好!

qq469015280

请问 ak 和bk 是时间系数   官网上写的是homogeneous 同类系数，那么异类时间系数在哪的呢？？

小迪

请问，楼主进行SVD分解之后，进行检验了吗？想求教怎样进行检验？非常感谢！

DevilMaySmile

qq469015280 发表于 2016-9-26 11:21
请问 ak 和bk 是时间系数   官网上写的是homogeneous 同类系数，那么异类时间系数在哪的呢？？

可能之前看到忘记回复了，我记得，好像时间系数就是左场一个，右场一个，好像时间系数不分同类和异类吧，同类和异类应该主要指最后的空间型，同类相关图和异类相关图，你再看书确认一下，我很久没用svd这个方法了，记不太清楚了

DevilMaySmile

小迪 发表于 2016-10-13 23:43
请问，楼主进行SVD分解之后，进行检验了吗？想求教怎样进行检验？非常感谢！

我最后没用svd，所以也没做检验，svd检验好像比较麻烦，不知道ncl里面有没有一些现成的函数能进行检验，怎么检验建议你看一下书吧，比如那本叫《气候变率诊断和预测方法》里面就有svd的检验。

leoaries

应该是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？

DevilMaySmile

leoaries 发表于 2017-2-6 17:30
应该是constant zero values 或者  constant 1 values  导致的
我实际测试过HadISST 的SIC，constant zero ...

感谢你的回答，我因为很久没有再做svd方面的东西，所以记得不是很清楚了，homogenous 和 heterogeneous 应该是指 同类 和异类，因为svd 的空间型的表示，有左右 同类 异类相关图(左同类，左异类，右同类，右异类) 实际分析，比较常用异类的， 就是左异类 hetlft 和右异类 hetrgt

郝完美

DevilMaySmile 发表于 2016-4-3 11:43
你好，是这样的，我想你可能想复杂了，其实就是用reshape函数，将原来，lat lon二维转为大小是nlat*nlon ...

楼主，请问具体的语句是怎样的啊？官网的例子没看懂。谢谢！

DevilMaySmile

郝完美 发表于 2018-9-19 14:06
楼主，请问具体的语句是怎样的啊？官网的例子没看懂。谢谢！

你是指reshape的具体语句吗？官网例子感觉已经非常清楚了