在R中通过reticulate调用Python脚本遇SystemExit:2错误求助

你的错误根源在于：你使用的Python脚本是依赖命令行参数运行的工具（从报错的usage信息能看出，它需要通过命令行传入path等必填参数），但reticulate::source_python()的作用是直接执行脚本内的代码，不会自动传递命令行参数，导致脚本启动时参数解析失败，触发退出错误。

下面提供两种可行的解决方式：

方式一：修改Python脚本，支持模块调用 #

把脚本的核心计算逻辑封装成函数，保留命令行解析的同时，允许从外部直接调用函数传入参数：

# raster_extract.py
import argparse

# 封装核心计算逻辑为函数
def calculate_precipitation(path, coords=None, input=None, odim=False, band=None, output=None, writer=None, decimals=None):
    # 这里放置原来脚本中处理降水量计算的核心代码
    # 示例：根据path和coords提取降水量并返回
    precipitation = 0.0  # 替换为实际计算结果
    return precipitation

# 保留原命令行执行逻辑，供直接运行脚本使用
if __name__ == "__main__":
    parser = argparse.ArgumentParser(description="Extract precipitation at specified coordinates")
    parser.add_argument('path', nargs='+', help="Path to raster file(s)")
    parser.add_argument('--coords', nargs='+', help="Coordinates (lon lat) to extract")
    # 添加其他原有的参数定义...
    args = parser.parse_args()
    
    # 调用封装的函数
    calculate_precipitation(**vars(args))

然后在R中调用这个封装好的函数：

setwd('C:/Users/username/Documents')
library(reticulate)
source_python("raster_extract.py")

# 传入实际参数调用函数，替换为你的文件路径和坐标
precip_result <- calculate_precipitation(
    path = "your_raster_data.tif",
    coords = c(120.1, 30.2)  # 经度、纬度
)

# 查看结果
print(precip_result)

方式二：在R中以命令行方式执行Python脚本 #

如果不想修改Python脚本，可以直接在R中通过命令行调用脚本并传入参数，用system()或processx包实现：

使用base R的system() #

setwd('C:/Users/username/Documents')
# 构造命令行命令，替换为你的实际参数
cmd <- 'python raster_extract.py "D:/data/raster.tif" --coords 120.1 30.2'
# 执行命令
system(cmd)

使用processx包（更安全，可捕获输出） #

setwd('C:/Users/username/Documents')
library(processx)

# 构造命令参数列表
args_list <- c(
    "raster_extract.py",
    "D:/data/raster.tif",  # path参数
    "--coords", "120.1", "30.2"  # 坐标参数
)

# 执行命令并获取结果
run_result <- run("python", args = args_list)

# 查看脚本输出
cat(run_result$stdout)
# 如果有错误输出，查看这里
cat(run_result$stderr)

内容的提问来源于stack exchange，提问作者Elke Debrie