0

我有以下两个 shapefile：

> summary(precincts1)

Object of class SpatialPolygonsDataFrame

Coordinates:
        min       max
x -74.25545 -73.70002
y  40.49613  40.91540

    Precinct        Shape_Leng       Shape_Area       
 Min.   :  1.00   Min.   : 17083   Min.   : 15286897  
 1st Qu.: 31.50   1st Qu.: 29900   1st Qu.: 37593804  
 Median : 64.50   Median : 46887   Median : 65891025  
 Mean   : 62.57   Mean   : 65720   Mean   :111231564  
 3rd Qu.: 95.50   3rd Qu.: 76375   3rd Qu.:133644443  
 Max.   :123.00   Max.   :309518   Max.   :781725787

和

> summary(bnd_nhd)

Object of class SpatialPolygonsDataFrame

Coordinates:
       min       max
x 871512.3  912850.5
y 982994.4 1070956.9

   SHAPE_area         SHAPE_len    
 Min.   : 3173813   Min.   : 7879  
 1st Qu.: 9687122   1st Qu.:13514  
 Median :14363449   Median :17044  
 Mean   :19674314   Mean   :19516  
 3rd Qu.:27161251   3rd Qu.:23821  
 Max.   :68101106   Max.   :49269

它们的坐标系不同。我可以用传单在地图上覆盖“区域 1”的形状，但我不能对“bnd_nhd”做同样的事情。我正在使用闪亮、地图工具和传单。如何转换 shapefile 或更改地图上的设置，以便为“bnd_nhd”覆盖地图？

4

1 回答 1

4

这应该有效：

library("rgdal")
library("leaflet")

bnd_nhd <- readOGR("C:/data/BND_Nhd88_cw.shp",
                  layer="BND_Nhd88_cw")
pol_wrd <- readOGR("C:/data/POL_WRD_2010_Prec.shp",
                  layer="POL_WRD_2010_Prec")

bnd_nhd4326 <- spTransform(bnd_nhd, CRS("+init=epsg:4326"))
pol_wrd4326 <- spTransform(pol_wrd, CRS("+init=epsg:4326"))

m <- leaflet() %>% 
  addTiles() %>%
  addPolygons(data=bnd_nhd4326, weight=2, color="red", group="bnd_nhd") %>%
  addPolygons(data=pol_wrd4326, weight=2, color="blue", group="pol_wrd") %>%
  addLayersControl(
    overlayGroups = c("bnd_nhd", "pol_wrd"),
    options = layersControlOptions(collapsed = FALSE)
  )
m

R 传单结果

于 2016-01-24T10:39:49.157 回答