作物模型狂奔 Apsim 脚本并行思路

2024-02-18 01:04

本文主要是介绍作物模型狂奔 Apsim 脚本并行思路,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

这里没图,去B站看看

整体思路:Python -》R -》Apsimx

0、准备工作

0.0 电脑环境

我用的Win10啦,Linux、Mac可能得自己再去微调一下。

PS:一直没上win11,因为win11 开始菜单用着实在难顶。
!在这里插入图片描述](https://img-blog.csdnimg.cn/direct/7758a7ad233946e89203f5fc5c9cc7b9.png)

0.1 安装Apsim 710

我这里ApsimNG 和 Apsim Classic 都有,用的多了发现NG好用,界面简洁,操作简单。

0.2 Python IDE

我用的Pycharm,个人感觉最好使的IDE,没有之一。

Python 解释器随便装个咯,我用的 Python 3.10 版本。

对于 Python,我其实也是个小白,只会写些小脚本,哈哈哈哈哈。

0.3 R IDE

继续用Pycharm,需要上 R环境插件。
PS:用过一天RStudio,感觉还是Pycharm顺手啊。尽管Pycharm R语言插件不支持绘图,这点坐等更新。对于 R,我其实还是个小白,这点代码全现百度现写的,哈哈哈哈哈。

R 解释器也随便装个咯,我用的 R 3.1.4 版本。

1、核心依赖包

1.1 rpy2(Python)

Python 环境下调用R包的依赖包。需要配一下R的环境变量,这里要注意哦,要不然会报错。

1.2 apsimx®

R环境下调用本地apsim的驱动包。(https://femiguez.github.io/apsimx-docs/apsimx.html)

2、Python部分编码思路

2.1 设置R环境变量

os.environ[‘R_HOME’] = ‘D:\application\R-4.1.3’

2.2 引包

import os
import pandas as pd
from pathlib import Path
from datetime import datetime
import rpy2.robjects as robjects

2.3 调用R部分编码

robjects.r’source’
apsimx_runner = robjects.r[‘apsimx_runner’]

result = apsimx_runner(2023, ‘mangshi’, ‘A’, param)

2.4 单线程

直接for,遍历就完了

for row in problems_parameters.values:
multi_run_func(row)

2.5 多进程

这里用的 进程池 ProcessPoolExecutor

pool = ProcessPoolExecutor(12)
for result in pool.map(multi_run_func, problems_parameters.values):
run_results.append(result)

3、R部分编码思路

3.1 函数封装

这里比较简单,封装了一个调用apsimx包的函数,留了接收作物模型参数的入口

apsimx_runner <- function (year, site.name, site.no, param) {}

3.2 函数体

apsimx 文件复制一个临时件,后续操作在这个临时文件上进行,不破坏原文件。

file.copy(ap_file.path, tmp_file.path)

设置气象数据

edit_apsimx(tmp_file.name, src.dir = tmp_file.dir, wrt.dir = tmp_file.dir, overwrite = TRUE, node = “Weather”, value = weather_file.path)

设置模型参数

edit_apsimx_replacement(file = tmp_file.name, src.dir = tmp_file.dir, wrt.dir = tmp_file.dir, overwrite = TRUE, node.string = p_path, parm = p_name, value = value, root = root, verbose = FALSE)

模型狂奔,这里用的Apsim NG

sim_result <- apsimx(tmp_file.name, src.dir = tmp_file.dir)

删除临时文件

file.remove(unlist(lapply(list.files(pattern = paste0(tmp_file.real_name, ‘.*’), path = tmp_file.dir), function(t) paste0(tmp_file.dir, ‘/’, t))))

返回结果

return(sim_result)
apsimx 包有使用demo,这里调用的几个函数都是根据demo调整的。

3.3 apsim classic

apsimx 包里 apsim 指代 Classic,apsimx 指代 NG

需要的话,自己研究下apsimx包里的apsim函数吧

3.4 R编码替代

上述代码可以在.R文件写好,然后在Python中调用;也可以直接写在Python代码里,这也是支持的。具体看rpy2的文档。

4、结果

4.1 格式

直接把结果接到list里

结果列表
数据项具体含义可以用 ApsimNg GUI 看看

4.2 耗时

跑一次大概10多秒,尽量开多进程或者多线程去跑,省点时间

欢迎各位道友关注、留言、私聊、交流病情。

如果需要的话去B站交流吧,平常不登CSDN

挂个小广告不会太过份吧~

需要指导的话,那就得让我挣点零花咯,嘿嘿。PS:不说虚的,程序狂奔才是最重要的。

Apsim classic 710、ApsimNg(ApsimX) 模型GUI运行、数据处理、脚本并行(R、Python)。

【闲鱼】https://m.tb.cn/h.5GmTkIi?tk=UcspW9B3gEY MF7997 「我在闲鱼发布了【apsim模型指导】」

这篇关于作物模型狂奔 Apsim 脚本并行思路的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



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

相关文章

Java的IO模型、Netty原理解析

《Java的IO模型、Netty原理解析》Java的I/O是以流的方式进行数据输入输出的,Java的类库涉及很多领域的IO内容:标准的输入输出,文件的操作、网络上的数据传输流、字符串流、对象流等,这篇... 目录1.什么是IO2.同步与异步、阻塞与非阻塞3.三种IO模型BIO(blocking I/O)NI

基于Flask框架添加多个AI模型的API并进行交互

《基于Flask框架添加多个AI模型的API并进行交互》:本文主要介绍如何基于Flask框架开发AI模型API管理系统,允许用户添加、删除不同AI模型的API密钥,感兴趣的可以了解下... 目录1. 概述2. 后端代码说明2.1 依赖库导入2.2 应用初始化2.3 API 存储字典2.4 路由函数2.5 应

通过Python脚本批量复制并规范命名视频文件

《通过Python脚本批量复制并规范命名视频文件》本文介绍了如何通过Python脚本批量复制并规范命名视频文件,实现自动补齐数字编号、保留原始文件、智能识别有效文件等功能,听过代码示例介绍的非常详细,... 目录一、问题场景:杂乱的视频文件名二、完整解决方案三、关键技术解析1. 智能路径处理2. 精准文件名

Java之并行流(Parallel Stream)使用详解

《Java之并行流(ParallelStream)使用详解》Java并行流(ParallelStream)通过多线程并行处理集合数据,利用Fork/Join框架加速计算,适用于大规模数据集和计算密集... 目录Java并行流(Parallel Stream)1. 核心概念与原理2. 创建并行流的方式3. 适

web网络安全之跨站脚本攻击(XSS)详解

《web网络安全之跨站脚本攻击(XSS)详解》:本文主要介绍web网络安全之跨站脚本攻击(XSS)的相关资料,跨站脚本攻击XSS是一种常见的Web安全漏洞,攻击者通过注入恶意脚本诱使用户执行,可能... 目录前言XSS 的类型1. 存储型 XSS(Stored XSS)示例:危害:2. 反射型 XSS(Re

C#集成DeepSeek模型实现AI私有化的流程步骤(本地部署与API调用教程)

《C#集成DeepSeek模型实现AI私有化的流程步骤(本地部署与API调用教程)》本文主要介绍了C#集成DeepSeek模型实现AI私有化的方法,包括搭建基础环境,如安装Ollama和下载DeepS... 目录前言搭建基础环境1、安装 Ollama2、下载 DeepSeek R1 模型客户端 ChatBo

SpringBoot快速接入OpenAI大模型的方法(JDK8)

《SpringBoot快速接入OpenAI大模型的方法(JDK8)》本文介绍了如何使用AI4J快速接入OpenAI大模型,并展示了如何实现流式与非流式的输出,以及对函数调用的使用,AI4J支持JDK8... 目录使用AI4J快速接入OpenAI大模型介绍AI4J-github快速使用创建SpringBoot

Python3脚本实现Excel与TXT的智能转换

《Python3脚本实现Excel与TXT的智能转换》在数据处理的日常工作中,我们经常需要将Excel中的结构化数据转换为其他格式,本文将使用Python3实现Excel与TXT的智能转换,需要的可以... 目录场景应用:为什么需要这种转换技术解析:代码实现详解核心代码展示改进点说明实战演练:从Excel到

Python脚本实现图片文件批量命名

《Python脚本实现图片文件批量命名》这篇文章主要为大家详细介绍了一个用python第三方库pillow写的批量处理图片命名的脚本,文中的示例代码讲解详细,感兴趣的小伙伴可以了解下... 目录前言源码批量处理图片尺寸脚本源码GUI界面源码打包成.exe可执行文件前言本文介绍一个用python第三方库pi

shell脚本自动删除30天以前的文件(最新推荐)

《shell脚本自动删除30天以前的文件(最新推荐)》该文章介绍了如何使用Shell脚本自动删除指定目录下30天以前的文件,并通过crontab设置定时任务,此外,还提供了如何使用Shell脚本删除E... 目录shell脚本自动删除30天以前的文件linux按照日期定时删除elasticsearch索引s