对pcse\wofost官方代码的一些改进:根据时期变化的变量的设置

2024-02-17 09:30

本文主要是介绍对pcse\wofost官方代码的一些改进:根据时期变化的变量的设置,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

这几天在学习WOFOST模型的过程中遇到了一个问题:对于根据时期变化的变量该怎么设置才可以成功进行 SOBOL 敏感性分析和 SUBPLEX 参数优化。苦于国内关于 PCSE 的资料甚少,大部分都是直接照搬官方文档,而官方文档并没有这方面的示例,困扰了挺久后问了一位大佬让我有了新思路,多亏了他才有我这篇文章。

由于本人资质尚浅,某些想法和文章中的话语难免有些幼稚,这个望各位不啬指正。

先上大佬的文章,他是利用simlab软件进行敏感性分析,然后用图片中的方法进行模型参数替换。

 但是我想用代码的方式实现 SOBOL 方法进行敏感性分析。还有就是本人想用面向对象的写法实现模型,尽量的精简代码,避免遇到不同的作物、不同的情况又得频繁的修改代码。所以在这里我详细地为大家介绍我的方法:

首先对于根据时期变化的变量我们可以分成多个变量

以 AMAXTB 为例,下图为官方示例的小麦作物参数文件

 可以看见这有四个数值,分别对应发育期为0.00、1.00、1.30、2.00的 AMAXTB 值。那我们在设置参数范围时,可以这样:

 分成四个变量,后面的数字代表发育阶段。每个变量对应四个范围(在这里我是四个[20,70])。

注意这个变量后的数字要和作物参数文件上的发育阶段表示方法一致,不然我后面修改的代码会无法成功匹配和替换。

设置完成后读取这个参数范围yaml文件,在敏感性分析代码(我将敏感性分析写成了一个单独的模块,各位根据自己的情况修改)中添加一个判断:

def Sensitivity(self,sampledata: str):# 使用 Saltelli 采样器采集生成参数集calc_second_order = Truensamples = ((65*(self.problem['num_vars'])//500)+1)*500print("样本数为 {}".format(nsamples))paramsets = saltelli.sample(self.problem, nsamples, calc_second_order=calc_second_order)print("We are going to do %s simulations" % len(paramsets))nsets = len(paramsets)printProgressBar(0, nsets, prefix = 'Progress:', suffix = 'Complete', length = 50)# 循环遍历参数值的样本,并为每个样本运行 WOFOSTtarget_results = []for i, paramset in enumerate(paramsets):self.params.clear_override()                 # 清除任何现有的覆盖 for name, value in zip(self.problem["names"], paramset):pername = re.findall(r"\d+\.\d+|[A-Z]+",name)if pername[0] in ["SLATB","KDIFTB","EFFTB","AMAXTB","TMPFTB"]:index = self.params[pername[0]].index(float(pername[1]))+1self.params[pername[0]][index] = valueelse:self.params.set_override(name, value)    # 更改每个 WOFOST 模拟的模型参数wofost = Wofost72_PP(self.params, self.wdp, self.agro) # 导入模型wofost.run_till_terminate()             # 运行模型r = wofost.get_summary_output()         # 返回已在模拟期间存储的汇总变量。target_result = r[0][self.target_variable]if target_result is None:print("Target variable is not available in summary output!")breaktarget_results.append(target_result)printProgressBar(i+1, nsets, prefix = 'Progress:', suffix = 'Complete', length = 50)target_results = np.array(target_results)# 计算灵敏度指数,使用 Sobol 方法来评估模型对参数变化的响应。Si = sobol.analyze(self.problem, target_results, calc_second_order=calc_second_order)

这句代码是利用正则(re库)对变量进行分割变成字符串和数字(包括小数)。这样做的目的是分割出来的字符串模型是可以正常识别的,而数字对应变量的发育阶段,之后可以依据这个数字进行值替换。

pername = re.findall(r"\d+\.\d+|[A-Z]+",name)

这句代码是进行判断,是否是根据时期变化的变量。列表中的变量是一些根据时期变化的变量,如果有其他可以自己增加。

if pername[0] in ["SLATB","KDIFTB","EFFTB","AMAXTB","TMPFTB"]:

这句代码是根据发育阶段得到该进行值替换的索引

index = self.params[pername[0]].index(float(pername[1]))+1

最后就可以成功分情况进行值替换啦。

self.params[pername[0]][index] = value

最后再次感谢大佬的启发!希望其他在学习WOFOST的小伙伴看到我这篇文章可以有所收获!

这篇关于对pcse\wofost官方代码的一些改进:根据时期变化的变量的设置的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/717358

相关文章

变量与命名

引言         在前两个课时中,我们已经了解了 Python 程序的基本结构,学习了如何正确地使用缩进来组织代码,并且知道了注释的重要性。现在我们将进一步深入到 Python 编程的核心——变量与命名。变量是我们存储数据的主要方式,而合理的命名则有助于提高代码的可读性和可维护性。 变量的概念与使用         在 Python 中,变量是一种用来存储数据值的标识符。创建变量很简单,

活用c4d官方开发文档查询代码

当你问AI助手比如豆包,如何用python禁止掉xpresso标签时候,它会提示到 这时候要用到两个东西。https://developers.maxon.net/论坛搜索和开发文档 比如这里我就在官方找到正确的id描述 然后我就把参数标签换过来

Android实现任意版本设置默认的锁屏壁纸和桌面壁纸(两张壁纸可不一致)

客户有些需求需要设置默认壁纸和锁屏壁纸  在默认情况下 这两个壁纸是相同的  如果需要默认的锁屏壁纸和桌面壁纸不一样 需要额外修改 Android13实现 替换默认桌面壁纸: 将图片文件替换frameworks/base/core/res/res/drawable-nodpi/default_wallpaper.*  (注意不能是bmp格式) 替换默认锁屏壁纸: 将图片资源放入vendo

poj 1258 Agri-Net(最小生成树模板代码)

感觉用这题来当模板更适合。 题意就是给你邻接矩阵求最小生成树啦。~ prim代码:效率很高。172k...0ms。 #include<stdio.h>#include<algorithm>using namespace std;const int MaxN = 101;const int INF = 0x3f3f3f3f;int g[MaxN][MaxN];int n

计算机毕业设计 大学志愿填报系统 Java+SpringBoot+Vue 前后端分离 文档报告 代码讲解 安装调试

🍊作者:计算机编程-吉哥 🍊简介:专业从事JavaWeb程序开发,微信小程序开发,定制化项目、 源码、代码讲解、文档撰写、ppt制作。做自己喜欢的事,生活就是快乐的。 🍊心愿:点赞 👍 收藏 ⭐评论 📝 🍅 文末获取源码联系 👇🏻 精彩专栏推荐订阅 👇🏻 不然下次找不到哟~Java毕业设计项目~热门选题推荐《1000套》 目录 1.技术选型 2.开发工具 3.功能

SWAP作物生长模型安装教程、数据制备、敏感性分析、气候变化影响、R模型敏感性分析与贝叶斯优化、Fortran源代码分析、气候数据降尺度与变化影响分析

查看原文>>>全流程SWAP农业模型数据制备、敏感性分析及气候变化影响实践技术应用 SWAP模型是由荷兰瓦赫宁根大学开发的先进农作物模型,它综合考虑了土壤-水分-大气以及植被间的相互作用;是一种描述作物生长过程的一种机理性作物生长模型。它不但运用Richard方程,使其能够精确的模拟土壤中水分的运动,而且耦合了WOFOST作物模型使作物的生长描述更为科学。 本文让更多的科研人员和农业工作者

代码随想录冲冲冲 Day39 动态规划Part7

198. 打家劫舍 dp数组的意义是在第i位的时候偷的最大钱数是多少 如果nums的size为0 总价值当然就是0 如果nums的size为1 总价值是nums[0] 遍历顺序就是从小到大遍历 之后是递推公式 对于dp[i]的最大价值来说有两种可能 1.偷第i个 那么最大价值就是dp[i-2]+nums[i] 2.不偷第i个 那么价值就是dp[i-1] 之后取这两个的最大值就是d

pip-tools:打造可重复、可控的 Python 开发环境,解决依赖关系,让代码更稳定

在 Python 开发中,管理依赖关系是一项繁琐且容易出错的任务。手动更新依赖版本、处理冲突、确保一致性等等,都可能让开发者感到头疼。而 pip-tools 为开发者提供了一套稳定可靠的解决方案。 什么是 pip-tools? pip-tools 是一组命令行工具,旨在简化 Python 依赖关系的管理,确保项目环境的稳定性和可重复性。它主要包含两个核心工具:pip-compile 和 pip

D4代码AC集

贪心问题解决的步骤: (局部贪心能导致全局贪心)    1.确定贪心策略    2.验证贪心策略是否正确 排队接水 #include<bits/stdc++.h>using namespace std;int main(){int w,n,a[32000];cin>>w>>n;for(int i=1;i<=n;i++){cin>>a[i];}sort(a+1,a+n+1);int i=1

uniapp设置微信小程序的交互反馈

链接:uni.showToast(OBJECT) | uni-app官网 (dcloud.net.cn) 设置操作成功的弹窗: title是我们弹窗提示的文字 showToast是我们在加载的时候进入就会弹出的提示。 2.设置失败的提示窗口和标签 icon:'error'是设置我们失败的logo 设置的文字上限是7个文字,如果需要设置的提示文字过长就需要设置icon并给