ERA-interim资料降水的格点值与站点观测值的对比差异太大

乔木

下载了1979.01.01-1981.12.31三年的ERA-interim资料tp（73*49*2196）值，三维分别表示经、纬度和时间，tp为每天两次的预报值，分辨率为0.125，地区在西藏；用资料算出1979年的站点附近的格点值大概是站点观测值的两倍多，比如一个站点1979年的年降水为500，该站点附近的格点值则大概为1000多一点，相差太大，个别站点附近的格点值达到该站点值的三倍，相差太大了；查找过一些可能的原因，比如数据类型，数据读取方式问题等，也可能是自己验证过程出错，但没能找到具体的原因；所以想请问下哪些大神了解这类资料？一般该情况的差异在哪个范围？年降水差异这么大的可能的原因是什么？做格点插值到站点与站点观测值的对比验证时，因为地形等原因误差应该也大，这如何是好？

andrewsoong

这个数据需要乘以scale，加上offset。也就是说需要将short型转为float才是真值。另外ERA的分析场是没有降水的。你的降水是预报场。预报的降水本来就差别大。还有就是青藏高原地区，卫星观测降水都不好，更何况是预报。

乔木

andrewsoong 发表于 2017-2-18 11:41
这个数据需要乘以scale，加上offset。也就是说需要将short型转为float才是真值。另外ERA的分析场是没有降水 ...

多谢回复！之前想过是否要做数值*scale+offset的转换，看帖子，提到用ncread打开不需要再进行转换，如果转换的话数值差别好像更大，用matlab转出来数据是double型，试过把它转成float型，看着没有多大变化；没有实测值，是预报值，高原地区预报值差异大怎么办？差别很大可以用做气候变化分析吗？主要是看趋势是否相同吗？

andrewsoong

乔木 发表于 2017-2-18 12:15
多谢回复！之前想过是否要做数值*scale+offset的转换，看帖子，提到用ncread打开不需要再进行转换，如果 ...

不转换肯定是不对的。至于能不能做气候变化分析，这个你要看数据精度。差别这么大，我看不大合适吧？你可以试试用基于观测的CRU

wuk

据我说知，对于netcdf数据，matlab自动处理*scale+offset，是不需要自己添加的，之前遇到过这个问题在论坛里问过。

乔木

wuk 发表于 2017-2-19 09:13
据我说知，对于netcdf数据，matlab自动处理*scale+offset，是不需要自己添加的，之前遇到过这个问题在论坛 ...

之前也是在论坛里看见过，应该不是scale与offset的问题，但是格点值与站点值差异太大了，太大，目前还不知道是什么原因.