[英]Replicating leafletjs interactive choropleth map example using leaflet package in R
[英]Interactive Choropleth in R of Sweden
我正在尝试在R的Shiny应用程序中开发一个交互式choropleth。我曾尝试使用glot,gVis和rCharts,但仍然没有任何运气。 我现在需要在瑞典对其进行可视化,但以后可能还要在其他国家/地区使用它。 到目前为止,这是我对gvisGeoMap的了解:
polygons <- readOGR("/ggshape", layer="SWE_adm1")
polygons <- fortify(polygons, region="ID_1")
data.poly <- as.data.frame(polygons)
data.poly <- data.poly[,c(1,2)]
data.poly.final <- data.frame(locationvar = paste(data.poly[,2],data.poly[,1], sep = ":"),
numvar=1,
hovervar="test")
data.poly.final$locationvar <- as.character(data.poly.final$locationvar)
data.poly.final$hovervar <- as.character(data.poly.final$hovervar)
map <- gvisGeoMap(data=data.poly.final, locationvar = "locationvar",
options=list(width='800px',heigth='500px',colors="['0x0000ff', '0xff0000']",
dataMode = "markers"))
plot(map)
根据文档,我应该能够在这里尝试使用经度和纬度坐标,但是我还没有成功。 我正在使用的shapefile来自http://www.gadm.org/download
基本上,没有人知道如何从gadm.org获得交互式可视化文件以与shapefile一起使用吗?
这就是我要用ggplot做到的
SWE <- fortify(polygons, region="ID_1")
SWEplot <- merge(x=SWE, y=my_data, by="id")
p <- ggplot() +
geom_polygon(data = SWEplot , aes(x = long, y = lat, group = group, fill = Patients)) +
geom_path(color="black") +
theme(axis.ticks.y = element_blank(),axis.text.y = element_blank(), # get rid of x ticks/text
axis.ticks.x = element_blank(),axis.text.x = element_blank(), # get rid of y ticks/text
plot.title = element_text(lineheight=.8, face="bold", vjust=1),
panel.background = element_blank(), panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.text=element_text(size=14),
legend.title=element_text(size=16)) + # make title bold and add space
coord_equal(ratio=1)
哪个产生
根据需要,但没有交互性。 我本质上希望实现的是类似http://rcharts.io/viewer/?6735051#.V1px-7t97mE的内容,但是当然是针对瑞典。
给定
library(raster)
swe <- getData("GADM", country = "SWE", level = 1)
swe$Patients <- runif(1:nrow(swe))
你可以做例如
library(maptools)
library(rgeos)
library(broom)
library(ggplot2)
library(plotly)
swe_s <- gSimplify(swe, .01)
SWE <- fortify(swe_s, region="ID_1")
SWEplot <- merge(x=SWE, y=swe, by.x="id", by.y="ID_1")
ggplot() +
geom_polygon(data = SWEplot , aes(x = long, y = lat, group = group, fill = Patients)) +
coord_quickmap() +
ggthemes::theme_map() + theme(legend.position=c(.8, .2)) ->
p
ggplotly(p)
要么
library(leaflet)
pal <- scales::seq_gradient_pal(low = "#132B43", high = "#56B1F7", space = "Lab")(seq(0, 1, length.out = 255))
leaflet() %>%
addPolygons(
data = swe,
color = "#000", weight = 1, opacity = 0.5,
fillColor = ~colorNumeric(pal, swe$Patients)(Patients), fillOpacity = 1,
popup = with(swe@data, htmltools::htmlEscape(sprintf("%s: %s", NAME_1, Patients)))
)
给定来自getData
的SpatialPolygosDataFrame,您可以使用library(mapview)
为您完成其余工作:
library(raster)
library(mapview)
swe <- getData("GADM", country = "SWE", level = 1)
mapview(swe)
要控制要绘制的属性,请使用zcol
参数。
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.