使用HWSD数据集，点击 Download database 和 HWSD Raster，下载得到原始全球的栅格数据及其对应的参数库。数据下载后即可在ArcGIS中打开查看。

绘制的分面图每个分面的标题标签包含特殊符号和下标，没办法正确显示，和chatgpt反复拉扯一个多小时，反复提问反复修改就是不能正确显示公式，网站搜索也无果搜到一个帖子，感觉好像看到了希望，结果一运行提示错误，是我的ggplot2版本更新的比较高了，如果有比较低的朋友就可以参考这个方法来进行修改如何在分面中添加数学表达式标签?_ggplot修改分面标签-CSDN博客Error:!variablean

在阅读文献时，看到了一个作者提供的，需要从github安装的包运行 devtools::install_github()开始安装，但出现问题，安装不成功，反复研究，终于找到问题所在，最终成功安装并library，非常激动！！

但需要注意的是，导入软件以及自己计算天气发生器所需参数时，不能包含缺失值，因此缺失值也得先简单处理一下。除太阳辐射外的气象数据直接整理完毕即可，太阳辐射没有直接的数据，需要基于日照时数进行计算，参考相关教程（原链接提供的excel表非常清楚，也很简洁明了，可以直接使用），由于数据较多，在R中完成批量处理（代码如下），最终得到对应的太阳辐射数据。当所需的六类数据中，有缺失项，或者是某一项中有缺失值时

使用HWSD数据集，点击 Download database 和 HWSD Raster，下载得到原始全球的栅格数据及其对应的参数库。数据下载后即可在ArcGIS中打开查看。

但需要注意的是，导入软件以及自己计算天气发生器所需参数时，不能包含缺失值，因此缺失值也得先简单处理一下。除太阳辐射外的气象数据直接整理完毕即可，太阳辐射没有直接的数据，需要基于日照时数进行计算，参考相关教程（原链接提供的excel表非常清楚，也很简洁明了，可以直接使用），由于数据较多，在R中完成批量处理（代码如下），最终得到对应的太阳辐射数据。当所需的六类数据中，有缺失项，或者是某一项中有缺失值时