單張上的R：如何在數據幀

爲每個因子水平層和顏色（順便說一下，有另一個包，會做我想要什麼？）單張上的R：如何在數據幀

這裏是我的數據集：

> head(df) 
    id groupID  lat  lon 
511 1 277706 -10.89834 -37.05893 
377 2 277706 -10.88870 -37.07695 
98 3 277705 -10.89050 -37.09406 
392 4 277697 -10.92131 -37.10525 
6 5 277705 -10.89050 -37.09406 
364 6 277697 -10.93730 -37.11600

我希望使用小冊子根據經緯度繪製地圖上的每一行。另外，每個標記（或彈出窗口或CircleMarker）應該具有基於groupID變量的不同顏色和圖層。

問題是數據每天都在變化，我不知道groupID變量有多少不同的獨特級別。在此示例數據集中，有12個級別，但範圍可能介於5到30之間。documentation中的示例使用預定義數量的級別。

這是我已經試過：

colorsmap = colors()[1:length(unique(df3$groupID))] 
groupColors = colorFactor(palette = "RdYlBu", domain = df3$groupID) 

leaflet(data = df3) %>% 
    addTiles() %>% 
    addCircleMarkers(lng = ~lon, lat = ~lat, color = ~groupColors(groupID), 
        group = ~groupID) %>% 
    #addLegend(position = "topright", pal = groupColors, values = ~groupID) %>% 
    addLayersControl(~groupID)

它確實提供了一個情節，但是當我選擇只有一個級別，其他的人不即消失，因爲他們應該：

圖像：

數據集本身：

> dput(df3) 
structure(list(id = 1:20, groupID = c(277698L, 277715L, 277704L, 
277706L, 277705L, 277705L, 277719L, 277705L, 277705L, 277709L, 
277706L, 277705L, 277704L, 277706L, 277715L, 277702L, 277719L, 
277706L, 277715L, 277706L), lat = c(-10.8172615660655, -10.8904055150991, 
-10.8887597563482, -10.90203509, -10.9001514, -10.8997748900025, 
-10.8960177351343, -10.8896179908615, -10.8991450456802, -10.9224848475651, 
-10.9000373151094, -10.8905013650562, -10.8889438100208, -10.90, 
-10.8861897462579, -10.9326053452642, -10.8916601751623, -10.902075281944, 
-10.8822231928033, -10.9079483812524), lon = c(-36.9248145687343, 
-37.0665064455395, -37.0921721937304, -37.05829295, -37.0969278, 
-37.0976847916125, -37.0840372102666, -37.0963566353117, -37.0945971936751, 
-37.0549293249471, -37.066113628594, -37.0940632483155, -37.095505683692, 
-37.0590422449149, -37.0782556623101, -37.0698746017798, -37.0841003949028, 
-37.0593363285999, -37.0724709841895, -37.0817244836096)), .Names = c("id", 
"groupID", "lat", "lon"), row.names = c(20L, 23L, 8L, 36L, 14L, 
13L, 16L, 2L, 11L, 1L, 26L, 6L, 5L, 31L, 22L, 50L, 17L, 34L, 
25L, 42L), class = "data.frame")

來源

2016-03-10 Sillas Teixeira Gonzaga

library(leaflet) 

groups = as.character(unique(df3$groupID)) 

map = leaflet(df3) %>% addTiles(group = "OpenStreetMap") 
for(g in groups){ 
    d = df3[df3$groupID == g, ] 
    map = map %>% addCircleMarkers(data = d, lng = ~lon, lat = ~lat, 
           color = ~groupColors(groupID), 
           group = g) 

} 
map %>% addLayersControl(overlayGroups = groups)

來源

2016-03-10 19:19:38 marbel

爲了完整，這是目前同時使用zcol和burst參數的MapView與輕鬆可行的：

library(mapview) 
library(sp) 

df3 <- structure(list(id = 1:20, 
         groupID = c(277698L, 277715L, 277704L, 
            277706L, 277705L, 277705L, 
            277719L, 277705L, 277705L, 277709L, 
            277706L, 277705L, 277704L, 277706L, 
            277715L, 277702L, 277719L, 
            277706L, 277715L, 277706L), 
         lat = c(-10.8172615660655, -10.8904055150991, 
           -10.8887597563482, -10.90203509, -10.9001514, 
           -10.8997748900025, -10.8960177351343, 
           -10.8896179908615, -10.8991450456802, 
           -10.9224848475651, -10.9000373151094, 
           -10.8905013650562, -10.8889438100208, 
           -10.90, -10.8861897462579, 
           -10.9326053452642, -10.8916601751623, 
           -10.902075281944, -10.8822231928033, 
           -10.9079483812524), 
         lon = c(-36.9248145687343, -37.0665064455395, 
           -37.0921721937304, -37.05829295, -37.0969278, 
           -37.0976847916125, -37.0840372102666, 
           -37.0963566353117, -37.0945971936751, 
           -37.0549293249471, -37.066113628594, 
           -37.0940632483155, -37.095505683692, 
           -37.0590422449149, -37.0782556623101, 
           -37.0698746017798, -37.0841003949028, 
           -37.0593363285999, -37.0724709841895, 
           -37.0817244836096)), 
       .Names = c("id", "groupID", "lat", "lon"), 
       row.names = c(20L, 23L, 8L, 36L, 14L, 13L, 16L, 2L, 11L, 
           1L, 26L, 6L, 5L, 31L, 22L, 50L, 17L, 34L, 25L, 42L), 
       class = "data.frame") 

## convert df3 to spatialPointsDataFrame for use with mapview 
coordinates(df3) <- ~ lon + lat 
proj4string(df3) <- "+init=epsg:4326" 

## now burst column "groupID" 
mapview(df3, zcol = "groupID", burst = TRUE)

注意：這是目前僅在開發版本mapview可以安裝

devtools::install_github("environmentalinformatics-marburg/mapview", ref = "develop")

來源

2016-05-24 15:47:36 TimSalabim

單張上的R：如何在數據幀

回答

相關問題