KHScrashes

library(leaflet) pal1 <- colorFactor(c(“blue”, “red”), domain = KHS01\(A30) pal2 <- colorFactor(c("green", "red"), domain = KS01\)A30)

awesome <- makeAwesomeIcon( icon = “car”, iconColor = “black”, markerColor = “red”, library = “fa” ) awesome1 <- makeAwesomeIcon( icon = “car”, iconColor = “blue”, markerColor = “orange”, library = “fa” ) ## leaflet() %>% addTiles()%>%##urlTemplate = “http://mt0.google.com/vt/lyrs=m&hl=en&x={x}&y={y}&z={z}&s=Ga”, attribution = ‘Google’) %>%# use the default base map which is OpenStreetMap tile # addTiles() %>%# use the default base map which is OpenStreetMap tile24.958964, 121.222348 setView(lng = 120.3281 , lat = 22.60327, zoom = 10)%>% addMarkers(data=KS01,lng = ~Lon, lat = ~Lat, #color = ~pal2(A30),opacity = 0.9,fill = TRUE, label = ~as.character(C.CrashPattern), popup = ~paste0(“案件編號:”,CaseNo,“
座標：”,Lat,” “,Lon,”
日期時間:“,Day,” “,Time,”
地點:“,Location,”
分隔型態:“,R.Divison,”;
事故型態：“,S.CrashPattern ,”
涉入車種：“,InvPV,sep=”“), group =”全部事故”, clusterOptions = markerClusterOptions())%>% addCircleMarkers(data=KHS01s,lng = ~Lon, lat = ~Lat, color = ~pal1(A30),opacity = 0.9,fill = TRUE, label = ~as.character(InvPV), popup = ~paste0(“案件編號:”,CaseNo,“
座標：”,Lat,” “,Lon,”
日期時間:“,Day,” “,Time,”
地點:“,Location,”
分隔型態:“,R.Divison,”;
事故型態：“,S.CrashPattern ,”
涉入車種：“,InvPV,sep=”“), group =”自撞事故”, clusterOptions = markerClusterOptions())%>% addAwesomeMarkers(data=KHS01d,lng = ~Lon, lat = ~Lat, icon = awesome, label = ~as.character(InvPV), popup = ~paste0(“案件編號:”,CaseNo,“
座標：”,Lat,” “,Lon,”
日期時間:“,Day,” “,Time,”
地點:“,Location,”
分隔型態:“,R.Divison,”;
事故型態：“,S.CrashPattern ,”
涉入車種：“,InvPV,sep=”“), group =”自撞A30”)%>% addHeatmap(data=KHS01,lng=~Lon, lat=~Lat, group=“熱力圖”, blur = 15, max = 0.05, radius = 10 )%>%

addCircleMarkers( data = centroids, lng = ~Lon, lat = ~Lat, # 圓圈大小：依然隨數量變大 radius = ~sqrt(CrashCount) * 3, # 顏色：根據 CrashCount 賦值 fillColor = ~pal(CrashCount), color = “red”, # 外框線顏色 weight = 1, # 外框線粗細 label = ~as.character(CrashCount), fillOpacity = 0.8, popup = ~paste0(“熱點等級：”, ifelse(CrashCount > 20, “極高風險”, “高風險”), “
”, “事故總數：”, CrashCount, ” 件
“,”座標：“, round(Lat, 5),”, “, round(Lon, 5)), group =”熱點中心” ) %>% addLegend(pal = pal1, values = KHS01$A30)%>%

addLayersControl( overlayGroups =c(“全部事故”,“自撞事故”,“自撞A30”,“熱力圖”,“熱點中心”), options = layersControlOptions(collapsed=FALSE) )