1.Generating the data for the cities and their populatión
lat<-c(4.598889,4.702354,4.861785,4.720331,4.297113,5.632386,4.434044,5.205835,4.592852)
lon<-c(-74.080833,-74.231873,-74.056778, -73.968356,-74.800415,-73.520508,-75.234375,-74.739990,-74.207153)
type<-c("city","city","city","city","city","city","city","city","city")
nameci<-c("Bogotá","Funza","Chía","La Calera","Girardot","Villa de leyva","Ibague","Honda","Soacha")
pop<-c(8080734,75350,29652,27527,97834,16984,553526,26873,522442 )
dfcities<-data.frame(matrix(nrow = 9,ncol=4))
dfcities[,1]<-lat
dfcities[,2]<-lon
dfcities[,3]<-type
dfcities[,4]<-nameci
dfcities[,5]<-pop
dfcities<-rename(dfcities,Lat=X1,Lng=X2,Type=X3,Nombre=X4,population=V5)
mymap<-leaflet(dfcities)%>%addTiles()%>%addCircleMarkers(lng=dfcities$Lng,lat=dfcities$Lat,weight=1,radius=sqrt(sqrt(dfcities$population)))2.Generating the points of interes in Bogotá
lat<-c(4.696391,4.813530,4.711613,4.707115,4.668868,4.628092,4.734820,4.647807,4.781130,4.741770)
lon<-c(-74.141387,-74.064895,-74.029798,-74.051898,-74.056735, -74.064058,-74.050635,-74.054328,-74.054393,-74.025913)
type<-c("Airport","Airport","Hospital","Hospital","Hospital","Hospital","School","School","School","School")
nameci<-c("<a href='https://eldorado.aero/en/'>El Dorado</a>","<a href='https://en.wikipedia.org/wiki/Guaymaral_Airport'>Guaymaral</a>","<a href='http://www.clinicaelbosque.com.co'>El Bosque</a>","<a href='https://portal.colsanitas.com/portal/web/clinica-reina-sofia'>La Reina Sofia</a>","<a href='http://www.clinicadelcountry.com'>El Country</a>","<a href='http://www.husi.org.co'>San Ignacio</a>","<a href='http://www.anglocolombiano.edu.co/index.php/es/'>Colegio Anglo Colombiano</a","<a href='http://www.cng.edu'>CNG</a>","<a href='http://www.nogales.edu.co/index.htm'>Los Nogales</a>","<a href='http://www.campestre.edu.co'>Campestre</a>")
colores<-c("red","red","green","green","green","green","purple","purple","purple","purple")
mymap<-addMarkers(mymap,lng = lon,lat = lat,group = type,popup = nameci)
mymap<-addCircleMarkers(mymap,lng=lon,lat=lat,radius=2,weight = 10,color=colores)%>%addLegend(labels=unique(type),colors=(unique(colores)))