使用 tmap 绘制空间点数据的问题

Question

我已经从 csv 文件创建了 class 大型 SpatialPointDataFrame，并且我正在尝试使用 tmap 对数据点进行 plot。 当我使用 plot() 绘制 SpatialPointDataFrame 时，一切看起来都很好。 但是，当我尝试使用 tmap 对其进行 plot 时，它只显示一个点，而且它不在应有的位置，顺便说一句，显示的点是数据框中的最后一项。 我在想我的投影仪可能有问题？ SpatialPointDataFrame 的 CRS（使用 SVY21 或 EPSG:3414）属性：

crs: +init=epsg:3414 +proj=tmerc +lat_0=1.366666666666667 +lon_0=103.8333333333333 +k=1 +x_0=28001.642 +y_0=38744.572 +ellps=WGS84 +units=m +no_defs

这就是我将 dataframe 转换为空间点数据帧的方式，其中列表是从 csv 文件中读取的原始 dataframe。

coords <- listing[ , c("latitude", "longitude")]
crs <- CRS("+init=epsg:3414")
listing_sp <- SpatialPointsDataFrame(coords = coords, listing, proj4string = crs)

关于我应该如何调试的任何建议？

Answer 1

在您的代码中，您有crs <- CRS("+init=epsg:3414") ，您正在分配一个投影，但您的源数据尚未转换为该投影。 如果您的源数据采用纬度/经度，则应与 WGS84 投影（ epsg 代码 4326 ）一起使用。

使用与坐标匹配的投影创建 spatialPointDataFrame 后，您只需使用spTransform重新投影。 要重新投影到 SVY21，您仍然需要进行转换。 没有一些可重复的数据就无法测试，但它应该看起来像这样。

coords <- listing[ , c("latitude", "longitude")]
crs <- CRS("+init=epsg:4326")
listing_sp <- SpatialPointsDataFrame(coords = coords, listing, proj4string = crs)
listing_sp_svy21 <- spTransform(listing_sp, CRS("+init=epsg:3414"))