1

我想在 NetLogo 中创建一个人口。因此，我想根据该地区的人口在地区内创建海龟。但是，我不完全确定如何做到这一点。

我得到了该地区的人口作为补丁值，如下所示：

gis:apply-coverage lor-dataset "GESBEV" population

但是当我用这个创建人口时，我得到的是每个区域内一个区域内的人口数量，而不是每个区域内的人口数量。

是否有可能每个地区只获得一次人口值？我也将感谢任何我可以做进一步阅读的资源。我自己没有找到任何东西。

4

1 回答 1

3

可能有很多方法可以解决这个问题，但这里有两种选择。首先，我在 Esri Open 数据中找到的一些示例数据——海地 2015 年人口数据。我将 shapefile 和相关文件提取到一个名为“gis”的文件夹中。我们将使用在 `gis:property-name': "POPULATION" 中找到的人口值。使用此设置：

extensions [ gis ]

globals [
  state-layer 
]

to setup
  ca
  resize-world 0 110 0 80
  set-patch-size 6.5
  set state-layer gis:load-dataset "gis/Population_2015.shp"
  gis:set-world-envelope gis:envelope-of state-layer
  gis:set-drawing-color white
  gis:draw state-layer 1
  reset-ticks
end

第一种选择是在每个特征的质心处只萌发整个种群（在这种情况下，除以 1000 以不产生太多海龟）。

to sprout-at-centroid
  foreach gis:feature-list-of state-layer [
    state ->

    ; Get the population for this state, divided by 1000
    ; to get a reasonable number
    let pop ( gis:property-value state "POPULATION" ) / 1000

    ; get the 'gis:location-of', an x-y list, for the 
    ; centroid of the current state / district
    let center gis:location-of gis:centroid-of state

    ; get the patch with the center xy values to sprout pop
    ask patch first center last center [
      sprout pop
    ]    
  ] 
end

输出如下所示： 在此处输入图像描述

看起来还不错！所有海龟都在每个要素的地理中心发芽。但是，根据您的数据集，您可能会遇到问题。请注意，中心的岛屿实际上是多部分要素的一部分，因此该多部分要素的人口已经在边界之外产生。这可能不会对您造成问题，具体取决于您所在地区的形状。

选项 2 稍微复杂一些，您可能会引入一些舍入误差。它也慢得多，并且根据您的世界和人口的大小可能需要很长时间/您可能会耗尽内存。首先，获取您的人口数量和您所在地区的补丁数量。然后，将人口除以该地区的斑块，得到每个斑块的平均人口。然后，让每个包含的补丁发芽该数量的海龟：

to apply-evenly
  foreach gis:feature-list-of state-layer [
    state ->

    ; Get the population for this state, divided by 10000
    ; to get a reasonable number
    let pop ( gis:property-value state "POPULATION" ) / 1000

    ; Get the patches contained by the state. This is slow!
    ; Using 'gis:intersecting' alone is much faster, but results
    ; in overlaps as geographic boundaries don't align with  patch boundaries
    let target-patches ( patches gis:intersecting state ) with [ gis:contained-by? self state ] 

    if any? target-patches [

      ; Get the number of turtles that should be in each target-patch:
      let avg-turtles round ( pop / count target-patches )

      ; Get the contained patches to sprout the appropriate number of turtles
      ask target-patches [
        sprout avg-turtles
      ]    
    ]
  ] 
end

这个输出看起来像：

在此处输入图像描述

但请注意，与我在这里作为示例使用的人口相比，它确实需要更长的时间，特别是因为我除以 1000。还要注意，如果您的补丁太大而无法被某个区域包含，您将不会有任何海龟在你的区域内产卵（例如南海岸的小岛）。

希望这能让你指出正确的方向！

于 2018-07-12T18:47:48.947 回答