近期问题汇总（五）

做这个系列的目的是记录，好记性不如烂笔头。时不时回头看也许能获取些灵感。

1. 话说有知道这个报错原因是什么的吗

库没装或者该更新了

2. 请教各位大佬，这种RMSE的图怎么画

rmse = np.sqrt(np.mean((np.array(actual_values) - np.array(predicted_values))**2)) 将对应的数据填入括号即可

3. 请问下各位老师，对全球大西洋海温研究中，大西洋的研究区域一般是怎么选取的呢？

这个推荐去查文献，可能有多种定义，还需自身决定用哪种 4. 请问大家nc数据带bnds维度存和直接用经纬度点存有什么不一样呀 存储空间占用差异巨大 5. 大家有啥利用本地文献pdf批量生成文献引用的工具 zotero可以自动识别外文期刊PDF

6. 最近，有老师用过CDS下载ERA5数据吗，我今天下载数据一直在排队中，以前从来没有遇到过排这么久的情况 欧洲中心转移服务器中，还需静待 7. 请问大家谁做过WRF的obsnuding吗？如何将地面站的观测数据csv转化为nuding所需要的little_R格式呀？ 去github上搜一下，我记得有对应的python脚本，直接little_r 8. 请问python调用的 matplotlib.cm下的colormap 怎么下载其rgb呢 ？

在Python中，你可以使用matplotlib库的`cm`模块来获取colormap。这个模块包含了许多预定义的colormap，你可以通过名字来获取。

例如，如果你想获取名为"viridis"的colormap，你可以这样做：

import matplotlib.cm as cm

cmap = cm.get_cmap('viridis')

然后，你可以使用`cmap(i)`来获取特定索引的颜色。例如，要获取索引为0.5的颜色，你可以这样做：

rgb = cmap(0.5)

这将返回一个RGBA元组，其中A（alpha）是透明度，R、G和B分别代表红色、绿色和蓝色的强度。

如果你只关心RGB值，你可以忽略alpha值：

rgb = cmap(0.5)[:3]

这将返回一个包含R、G和B值的元组

9. 请问大家可以读取已知经纬度的组合反射率吗？

使用pycinrad或pycwr读取获得的雷达变量已经是xarray格式，就当普通的数据去索引即可 不知怎么用xarray可查看往期推文如何快速熟悉一个陌生的nc格式数据

10. 大佬们，之前都是z-x、z-y垂直剖面，有没有画过z-沿着想要的水平线做切割的垂直剖面呀，有相关代码吗~ 可参考metpy的官方示例

https://unidata.github.io/MetPy/latest/examples/cross_section.html#sphx-glr-examples-cross-section-py

11. 请问一下wrf中的诊断输出变量pressure（full model pressure）是表示大气总气压（干燥气压加水汽压）吗？

no

12. 请问各位老师，如何将tif格式遥感数据从unit16转换为float32. 参考往期如何转换tif格式遥感数据的数据类型

13. 哪些是pycwr有而pycinrad没有的呢？

比较直观来说，是CAPPI算法。不过这个实现起来比较容易，就是对三维数据进行插值，考虑放到有空去搞搞 14. 最近要从wrfout文件里提取风剖面，求问各位大佬有什么推荐的帖子吗画

https://www.heywhale.com/mw/project/618e8ca486227300178d33df

15. 请问各位大佬，卫星降水数据可以代替站点数据比较嘛

结论，一般而言和实况差距是比较大的，毕竟降水数据是反演的。

16. 想问一下大家降水的站点数据能在哪里找到

答案是ISD (Integrated Surface Dataset，全球地面观测数据集）

网址是

https://www.ncei.noaa.gov/access/search/data-search/global-hourly

但是数据的处理可能较为困难，如果你是小白，推荐去和鲸社区的这个训练营练练。

17. 想请问下用python的 metpy.calc.gradient 梯度函数做计算的时候，二维数组（lat,lon）带入，对应这里的(y,x)。那么这里所得的 【'Gradient in y direction: ', grad[0]】，对应的是不是 V（经向风方向）？还是说U（纬向风方向）？

是V

18. 请问怎么计算能把气压值转换成高度值呢

metpy有函数

19. 请问各位老铁，有人知道这种图叫什么吗

shap值图

20. 请问python可以用shp截取数据吗，好奇雷达数据core.dataset.Dataset格式的可以剪裁么

震惊！雷达数据绘图居然能白化

好，下期攒够二十个问题，咱们再见