导入的GIS文件和网络之间未对齐节点具有X和Y坐标-NetLogo

如何解决导入的GIS文件和网络之间未对齐节点具有X和Y坐标-NetLogo

我，再次有一个问题。我以前对问题的解决方案未按预期工作。提醒您，我导入了一些多边形，并在适当地设置NetLogo世界中的某些海龟（城市）时遇到了问题（要从非gis文件中读取GIS坐标并设置其在世界中的位置，请参见Netlogo doesen't recognize the GIS coordinates in the defined envelope of the world,it treats them as netologo world coordinates）。因此，最后，我想出了一种将NetGIS坐标转换为NetLogo坐标的解决方案（nl-x和nl-y过程）。这些城市实际上是带有x和y属性的graphml格式。所以，我的代码是这样的：

extensions [nw gis]
directed-link-breed [highways highway]
breed [cities city]
highways-own [ name ]
cities-own [ x y ]
globals [ paldrino ]

to setup
  ca
  ;; set the world envelope
  gis:load-coordinate-system "wgs84.prj"
  set paldrino gis:load-dataset "paldrino.shp"               
  let world ( gis:envelope-of paldrino )
  gis:set-world-envelope (world)

  ;; Make them visible
  foreach gis:feature-list-of paldrino  [ ;for each polygon
  polygon ->
   let temp-color one-of base-colors
      ask patches gis:intersecting polygon [
      set pcolor  temp-color
   ]]

  ;; load network
  nw:set-context cities highways
  nw:load-graphml "highway-network.graphml"
  ask cities[
    set xcor nl-x(read-from-string x)  ;; if I put set xcor read-from-string  x,then it will put all the nodes in one point in Netlogo world,same for setxy fix-x(read-from-string x) fix-y (read-from-string y)
    set ycor nl-y (read-from-string y)
  ]

end

to-report nl-x [#x] 
   let world gis:envelope-of paldrino
   let minx item 0 world 
   let maxx item 1 world 
   report ((#x - minx) / (maxx - minx)) * (max-pxcor - min-pxcor) + min-pxcor 
end 

 to-report nl-y [#y] 
   let world gis:envelope-of paldrino 
   let miny item 2 world 
   let maxy item 3 world 
   report ((#y - miny) / (maxy - miny)) * (max-pycor - min-pycor) + min-pycor 
 end  

最后，问题在于NetLogo不能完美地对齐这两个文件：多边形（paldrino）和网络（城市）。大概让他们觉得还可以，但我需要让他们在世界上完美设置。例如，某些城市本应位于一个多边形中，但它们位于相邻的多边形中。 因此，例如，这是NetLogo放置它们的方式：

这就是应该如何使用它们（在QGIS中阅读）：

如果有人可以帮助我或将我指向另一个方向。我开始失去理智。

谢谢！