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程序进行检查)<代码>无法打开文件
:请确保路径/文件夹和文件名正确