Python调用GRASS GIS导入grass_session遇RuntimeError报错求助

核心错误点 #

• 未设置GRASSBIN环境变量：grass_session模块依赖GRASSBIN环境变量定位GRASS启动脚本（即你指定的grass78.bat），你代码中仅定义了grass7bin变量，但未将其注入系统环境变量。
• GIS数据库路径无效：os.path.join(os.path.expanduser("~"), "C:/Users/user/Documents/grassdata")会生成冲突路径（如C:\Users\当前用户\C:/Users/user/Documents/grassdata），因为第二个参数是绝对路径，join操作会直接忽略前置路径。
• 过早删除配置文件：调用gsetup.init生成rcfile后立刻删除，导致后续导入grass_session时GRASS环境配置已失效。

修正后的代码 #

import os
import sys
import subprocess

# 1. 指定GRASS启动脚本并设置GRASSBIN环境变量
grass7bin = r'C:/Program Files/GRASS GIS 7.8/grass78.bat'
os.environ['GRASSBIN'] = grass7bin

# 验证GRASS路径（可选，提前排查启动脚本问题）
startcmd = [grass7bin, '--config', 'path']
try:
    p = subprocess.Popen(startcmd, shell=False,
                         stdout=subprocess.PIPE, stderr=subprocess.PIPE)
    out, err = p.communicate()
except OSError as error:
    sys.exit(f"ERROR: 找不到GRASS GIS启动脚本 {startcmd[0]}: {error}")
if p.returncode != 0:
    sys.exit(f"ERROR: 运行GRASS GIS启动脚本失败 {' '.join(startcmd)}: {err}")

gisbase = r'C:/Program Files/GRASS GIS 7.8'
os.environ['GISBASE'] = gisbase
grass_pydir = os.path.join(gisbase, "etc", "python")
sys.path.append(grass_pydir)

# 2. 修正GIS数据库路径（二选一）
gisdb = r"C:/Users/user/Documents/grassdata"
# 或者 gisdb = os.path.join(os.path.expanduser("~"), "Documents", "grassdata")

location = "Testlocation"
mapset = "test"

# 3. 先导入grass_session，用上下文管理器自动管理会话
import grass_session as session

with session.Session(gisdb=gisdb, location=location, mapset=mapset, grassbin=grass7bin):
    # 在此处调用GRASS函数
    print("GRASS会话已成功启动")

# 无需手动删除配置文件，上下文管理器会自动清理

额外说明 #

• 优先使用grass_session的上下文管理器（with语句），它会自动处理环境配置与会话清理，比手动调用gsetup.init更稳定。
• 确保Testlocation位置和test地图集已在指定的gisdb路径下存在，若不存在需提前通过GRASS GUI或grass.script命令创建。

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