为何CRS 4326的GIS Shapefile在NetLogo中无法正上方视角显示？

我将两个Shapefile加载至NetLogo中：一个原始投影为EPSG:4326（OSM数据），另一个原始为EPSG:28992（荷兰政府数据），已在R中通过st_transform(crs = 4326)转换投影。但在NetLogo中无法以正上方视角显示，R ggplot中却可以正常展示。

问题原因 #

• 投影处理逻辑差异：R的ggplot默认将EPSG:4326的经纬度数据当作平面笛卡尔坐标渲染，直接输出正上方视角；而NetLogo的GIS扩展默认会把WGS84（4326）数据识别为球面坐标，自动应用球面投影（如墨卡托），导致画面呈现倾斜的球面视角。
• 投影元数据残留：转换后的Shapefile仍保留EPSG:4326的投影元数据，NetLogo读取时会优先依据元数据启用球面投影，而非平面渲染。

解决办法 #

方法1：强制NetLogo使用平面坐标范围 #

加载数据前，手动设置世界范围为数据集的平面边界，强制关闭球面投影逻辑：

; 加载第一个数据集获取边界
gis:load-dataset "transformed-dataset-1.shp"
let dataset-bounds gis:envelope-of gis:feature-list-of gis:current-dataset
; 设置世界范围为数据集的平面边界
gis:set-world-envelope-ds dataset-bounds
; 加载第二个数据集
gis:load-dataset "transformed-dataset-2.shp"

方法2：在R中导出无投影元数据的Shapefile #

转换投影后移除文件的投影信息，让NetLogo将其视为纯平面坐标：

library(sf)
# 读取原始数据并转换至4326
nl_gov_data <- st_read("original-28992-data.shp") %>% st_transform(crs = 4326)
# 清除投影元数据
st_crs(nl_gov_data) <- NA
# 导出无投影的Shapefile
st_write(nl_gov_data, "flat-projection-data.shp")

将导出的文件导入NetLogo即可获得正上方视角。

方法3：显式设置NetLogo为平面投影 #

通过命令直接指定NetLogo使用平面投影模式：

gis:set-projection "flat"
; 依次加载两个数据集
gis:load-dataset "transformed-osm-data.shp"
gis:load-dataset "flat-projection-data.shp"

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