Shiny 形状文件:XY坐标和经纬度坐标

Shiny 形状文件:XY坐标和经纬度坐标,shiny,leaflet,shapefile,maptools,Shiny,Leaflet,Shapefile,Maptools,我有以下两个形状文件: > 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

我有以下两个形状文件:

> 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  

它们的坐标系是不同的。我可以在地图上用传单覆盖“选区1”的形状,但我不能用“bnd_nhd”覆盖。我正在使用shiny、maptools和传单。如何转换形状文件或更改地图上的设置,以便覆盖“bnd_nhd”的地图

这应该有效:

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
库(“rgdal”)
图书馆(“单张”)
这应该是可行的:

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
库(“rgdal”)
图书馆(“单张”)

你必须把这些坐标转换成一个通用的坐标系,这个坐标系也适用于传单,看看这个:我尝试了以下代码:>bnd_nhd proj4string(bnd_nhd)bnd_nhd_latlon虽然转换成功,但多边形位于错误的位置。我不知道如何正确地投影和转换文件。这是密苏里州圣路易斯市的ESRI形状文件,从该文件中获得的是名为“邻里与病房”的“GIS信息”。如果有人对解决这个问题有任何进一步的建议,请告诉我。谢谢。你必须将这些坐标转换成一个通用坐标系,该坐标系也适用于传单,看看这个:我尝试了以下代码:>bnd_nhd proj4string(bnd_nhd)bnd_nhd_latlon虽然转换成功,但多边形位于错误的位置。我不知道如何正确地投影和转换文件。这是密苏里州圣路易斯市的ESRI形状文件,从该文件中获得的是名为“邻里与病房”的“GIS信息”。如果有人对解决这个问题有任何进一步的建议,请告诉我。谢谢,非常感谢!!我试过了,但是我一直得到以下错误。您能在这里指出问题吗?>bnd_nhd>bnd_nhd
未找到图层
:图层名称区分大小写,“bnd_Nhd88_cw”与“bnd_Nhd88_cw”不同(您可以使用QGIS或ArcGIS Desktop之类的GIS程序进行检查)
无法打开文件
:请确保路径/文件夹和文件名正确无误,非常感谢!!我试过了,但是我一直得到以下错误。您能在这里指出问题吗?>bnd_nhd>bnd_nhd
未找到图层
:图层名称区分大小写,“bnd_Nhd88_cw”与“bnd_Nhd88_cw”不同(您可以使用QGIS或ArcGIS Desktop之类的GIS程序进行检查)<代码>无法打开文件
:请确保路径/文件夹和文件名正确