library(leaflet)
leaflet()%>%addTiles()%>%setView(lng = 121.54, lat = 25.02, zoom = 14)
library(rgdal)
setwd("D:/1082SA/Data")
station=readOGR("EPA_STN1.shp",encoding = "utf8")
station=spTransform(station,CRS("+init=epsg:4326"))
leaflet()%>%addTiles()%>%addMarkers(data=station,popup =~SiteName)