代码之家 › 专栏 › 技术社区 › Dhiraj

在传单中绘制sfc\u多边形

sf leaflet r

6

Dhiraj · 技术社区 · 7 年前

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 sf layer is not long-lat data 但事实并非如此,我得到了基本的世界地图。这是我用来绘图的代码 传单 .

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

st_set_crs(IHS, 4326) 这改变了我对从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

1 回复 | 直到 7 年前

1

7

Sébastien Rochette 7 年前

“转换”为Lon/Lat是指改变其 空间投影 . 一旦您知道了这一点,就可以使用空间数据并在库中找到正确的投影空间特征的方法。

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

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