spNetwork 的项目 SpatialLinesDataFrame答案

【问题标题】：Project SpatialLinesDataFrame for spNetworkspNetwork 的项目 SpatialLinesDataFrame
【发布时间】：2022-01-20 14:18:02
【问题描述】：

我正在尝试使用spNetworks 的nkde 函数在DC 沿其道路网络创建崩溃的KDE。为了准备函数，我正在为nkde 创建 lixels，但运行 lines_center(lixels) 总是给我一个错误。

lixels <- lixelize_lines(dc_lines,1000,mindist = 250)
samples <- lines_center(lixels)

每当我尝试运行 lines_center(lixels) 时，我都会收到以下错误：

Error in maptools::SpatialLinesMidPoints(with_length) : 
  is.projected(sldf) is not TRUE
In addition: Warning message:
In RGEOSMiscFunc(spgeom, byid, "rgeos_length") :
  Spatial object is not projected; GEOS expects planar coordinates

我尝试查找各种技术，例如将 SpatialLinesDataFrame 转换为普通数据帧 st_as_sf，然后对其进行投影，但到目前为止没有任何结果，而且我总是遇到同样的错误。

我正在加载这样的数据：

dc <- readOGR("assessment/test/Roads_2013", "Roads_2013")

由于上传的文件是大型 SpatialPolygonsDataFrame，我正在使用以下代码将其转换为 SpatialLines：

dc_lines <- as(dc, "SpatialLinesDataFrame")

知道我做错了什么或如何正确投影线条吗？

使用的 shapefile 在这里： https://opendata.dc.gov/datasets/roads

【问题讨论】：

欢迎来到 Stack Overflow！请收下tour 并阅读help center，如果您有任何问题，请在Meta Stack Overflow 上提问:)
@EkadhSingh-ReinstateMonica 你建议他们在 meta 上问什么？如果他们需要改进问题，他们应该在这里做，而不是通过询问这个问题或关于元的一般操作问题
@camille 如果他们对规则有任何疑问，他们应该询问 meta 是我的意思。不过我可能并不完全清楚。
@EkadhSingh-ReinstateMonica 我不建议盲目地在 meta 上发布有关网站规则的问题。阅读帮助中心和现有帖子以查看您的问题是否已经得到解答以及将一个好的问题放在一起的相同准则适用于元数据，就像他们在这里所做的那样，但在那里他们没有那么温和地执行。也就是说，这将是新用户获得大量反对票的简便方法
很好，您已经包含了指向数据的链接，但完整的 reproducible example 可能会包含您阅读它的方式，特别是因为问题似乎与您的 CRS 相关

标签： r spatial kde rgdal

【解决方案1】：

我是 spNetwork 开发人员。请考虑在 github 上发布您的问题以获得更快的响应。这里的问题是您的数据集不使用平面（以米为单位的 X/Y 坐标）CRS，而是使用地理坐标（以度为单位的经度/纬度）。您需要使用 sp 包中的函数 spTransform 将数据重新投影到适当的 CRS 中。这是一些示例的链接 (https://www.rdocumentation.org/packages/rgdal/versions/1.5-28/topics/spTransform-methods)。

也许 EPSG:2927 可能是要使用的 CRS（我不熟悉在美国使用的 EPSG）。

dc <- readOGR("assessment/test/Roads_2013", "Roads_2013")
dc_proj <- spTransform(dc, CRS("+init=epsg:2927"))

【讨论】：

这对我有用，非常感谢！！