在传单中绘制 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. 当我绘制 this is 时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 就是改变它的空间投影。一旦你知道了这一点，你就可以使用空间数据并在库中找到投影空间特征的正确方法。

在您的情况下， usingst_set_crs(IHS, 4326)仅表示投影是 Lon/lat，而不投影它​​。您需要使用st_transform：

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