如何在NetLogo模型运行时自动切换GIS数据文件

首先确保你已在NetLogo中启用gis扩展，所有GIS数据的坐标参考系、空间范围匹配。

1. 封装GIS数据加载函数 #

将两个季节的GIS数据加载逻辑拆分为独立可调用的过程，避免和全局初始化逻辑耦合：

extensions [ gis ]
globals [ gis-season1 gis-season2 current-season ]

; 加载第一个季节GIS数据并初始化景观
to load-season1-env
  set gis-season1 gis:load-dataset "替换为你的GIS数据1路径"
  gis:set-world-envelope gis:envelope-of gis-season1
  ask patches [
    ; 替换为你原来的景观属性赋值逻辑，比如资源量、植被类型等
    set resource gis:cell-value gis-season1 self
  ]
  set current-season "season1"
end

; 加载第二个季节GIS数据并初始化景观
to load-season2-env
  set gis-season2 gis:load-dataset "替换为你的GIS数据2路径"
  gis:set-world-envelope gis:envelope-of gis-season2
  ask patches [
    set resource gis:cell-value gis-season2 self
  ]
  set current-season "season2"
end

2. 封装单季节运行逻辑 #

将单轮迭代的运行、终止判断封装为可传入参数的过程：

; 运行单季节模拟，参数为该季节需要运行的总tick数
to run-season [ season-ticks ]
  reset-ticks
  while [ ticks < season-ticks ] [
    go ; 调用你现有编写的模型单步运行逻辑，无需修改原有代码
    tick
  ]
end

3. 编写串联双阶段的主调用函数 #

核心自动切换逻辑如下，可直接点击运行实现自动流转：

to run-full-cycle
  clear-all
  ; 运行第一个季节
  load-season1-env
  run-season 180 ; 替换为第一个季节对应的迭代次数
  
  ; 切换到第二个季节：仅重置景观属性，保留蜜蜂等agent的状态
  ask patches [ set resource 0 ]
  load-season2-env
  run-season 185 ; 替换为第二个季节对应的迭代次数
  
  print "双季节模拟全部运行完成"
end

常见适配调整点 #

• 如果需要保留第一个季节运行后的蜜蜂种群、位置、能量等状态，切换数据时不要调用clear-all，仅重置patches的环境属性即可
• 若两个GIS数据的栅格分辨率、行列数不一致，需先在GIS工具中做重采样对齐，避免patch赋值错位
• 若需要跑多年的双季节循环，可在run-full-cycle外层套循环逻辑即可

内容的提问来源于stack exchange，提问作者Raïssa Fernandes