将栅格 gis 数据应用于补丁 netlogo 的更快方法

Question

有没有办法缩短将栅格数据导入到面片时的计算时间？

这是我的代码。这通常需要大约 20 分钟才能完成。我的世界是 500x500。我发现使用 gis:intersect? 比 gis:apply-coverage 快得多，如果只考虑 true 或 false。干杯。

 to setup-gis
  set city gis:load-dataset "GSR_GIS.shp"
  set LGA_nodes gis:load-dataset "LGA_nodes.shp"
  set builtuparea gis:load-dataset "GSR_builtuparea.shp"
  set recreationalarea gis:load-dataset "GSR_recreationalareas.shp"
  set natural gis:load-dataset "GSR_natural.shp"
  set reserves gis:load-dataset "GSR_reserves.shp"
  set rail_network gis:load-dataset "Greater_Sydney_rail_networkt.shp"
  set roads gis:load-dataset "GSR_road_network.shp"

  gis:set-world-envelope gis:envelope-of city 
  gis:set-drawing-color white 
  gis:draw city 1 
  gis:apply-coverage city "LGA" zone 

  set patchesinlga (patch-set patches with [ zone > 0 ])

      ask patchesinlga [
    ifelse gis:intersects? builtuparea self [set builtuparea? true][set builtuparea? false]
    ifelse gis:intersects? recreationalarea self [set recreationalarea? true][set recreationalarea? false]
    ifelse gis:intersects? natural self [set natural? true][set natural? false]
    ifelse gis:intersects? reserves self [set reserves? true][set reserves? false]
    ifelse gis:intersects? roads self [set roads? true][set roads? false]
    ifelse gis:intersects? rail_network self [set rail? true][set rail? false]
    if gis:intersects? LGA_nodes self [sprout-LGAs 1 [set color red set shape "flag" set size 2 set LGAid [zone] of patch-here]]
    ]

    ask patches with [zone > 0 and (round([zone] of self) / [zone] of self) != 1][set zone 0] ;; fixed the zone floating point issue
end

*我的编辑注意事项：我找到了一种通过仅询问感兴趣的补丁来加快速度的方法。在我的例子中，它们是带有 LGA 的补丁（LGA 是一个分区 ID），因此我的 GIS 地图之外的补丁将被排除在底部的运行 ifelse 循环中。

但是有什么办法可以改善吗？

一些运行次结果

50*50用了23秒
100*100 用了 79.5 秒
400*400 花了 30 分钟

Answer 1

我遇到了同样的问题：我的代码需要 30 多分钟才能完成，而我的世界只有 500x700 的补丁。延迟不是加载数据，而是在我的例子中，显示数据。看看你上面的代码，我相信你也遇到了类似的现象。

在您的 askpatchesinlga[] 过程中，您要求每个补丁检查它是否与数据相交。我认为这就是需要这么长时间才能执行的原因。 2500 个代理，每一个仅需 45 分钟的处理时间。为了加快速度，事先创建光栅，并将它们分别加载到 netlogo 中。这会将 2500 个命令变成 7（ish），并且应该相应地减少你的处理时间。

总而言之：如果您正在加载数据，请不要要求每个补丁都做某事。首先进行处理，然后加载完成的数据集 - 在您的例子中，是建筑区域栅格、娱乐区域栅格、自然区域栅格等。任何地理处理工具包都应该能够为您制作这些。

将栅格 gis 数据应用于补丁 netlogo 的更快方法

Faster way to apply raster gis data to patches netlogo

gis

netlogo