r - 在传单中绘制 sfc_POLYGON

Question

我有一个 sfc_POLYGON，IHS使用 sf 作为几个形状的联合获得包裹。

Geometry set for 1 feature 
geometry type:  POLYGON
dimension:      XY
bbox:           xmin: 270194.7 ymin: 2255850 xmax: 279225 ymax: 2265469
epsg (SRID):    32643
proj4string:    +proj=utm +zone=43 +datum=WGS84 +units=m +no_defs
POLYGON((279224.998943279 2261084.1514869, 2792...

虽然我可以在 mapview 中绘制它，我在 leaflet 中绘制相同内容时遇到了一些问题.当我绘制这是 leaflet , 我收到警告 sf layer is not long-lat data事实并非如此，我得到了基础世界地图。这是我用于在 leaflet 中绘图的代码.

IHS%>%leaflet()%>%addTiles()%>%
  addProviderTiles(providers$Esri.WorldImagery, group ="ESRI")%>%
  addPolygons()

我怀疑这与CRS有关。我试过 st_set_crs(IHS, 4326)这改变了我对 IHS 的预测从 utm 到 longlat:

Geometry set for 1 feature 
geometry type:  POLYGON
dimension:      XY
bbox:           xmin: 270194.7 ymin: 2255850 xmax: 279225 ymax: 2265469
epsg (SRID):    4326
proj4string:    +proj=longlat +datum=WGS84 +no_defs

但是，传单警告和输出没有变化。我不确定这里，但它与改变尺寸和 bbox (似乎以米为单位)有关吗？

Answer 1

“转换”成 Lon/Lat 是改变它的 空间投影 .一旦你知道了，你就可以使用空间数据并在图书馆帮助中找到投影空间特征的正确方法。

在您的情况下，使用 st_set_crs(IHS, 4326)只说投影是 Lon/lat，没有投影它。
您需要使用 st_transform :

IHS_wgs84 <- st_transform(IHS, "+init=epsg:4326")