mapa

apiKey <- "AIzaSyBPxhDW9_bUJ1juxXEGvY_QcNIZNJPUy7A"

library(ggplot2)
library('ggmap')

## Warning: package 'ggmap' was built under R version 4.0.5

## Google's Terms of Service: https://cloud.google.com/maps-platform/terms/.

## Please cite ggmap if you use it! See citation("ggmap") for details.

register_google(apiKey)

qmap(c(-76.5, 3.43),zoom=12, maptype="hybrid")

## Source : https://maps.googleapis.com/maps/api/staticmap?center=3.43,-76.5&zoom=12&size=640x640&scale=2&maptype=hybrid&language=en-EN&key=xxx

#gasetero
#qmap("Colombia",zoom=12, maptype="hybrid")

#leer los datos de las estaciones de cali
paradasMio<- read.delim("C:/Users/Solutian/Desktop/Visualizacion/mapa/mapa/MIO_STOPS.txt")

#visualizar las estaciones de cali Mio en el mapa

#mapa<- ggplot(paradasMio, aes(DECIMALLONGITUDE, DECIMALLATITUDE))

## adiccionar columna de agrupamiento
k <- 7
nIteraciones <- 100


grupos <- kmeans(paradasMio[,c(7:8)], k , iter.max = nIteraciones, nstart = 1)

paradasMio$grupos <- grupos$cluster

# poligono

cluster <- subset(paradasMio, paradasMio$grupos == 1)

indicesC1 <- chull(cluster[,c(7,8)])

poligono <- paradasMio[indicesC1,]

city<- get_map(c(-76.5, 3.43),zoom=12)

## Source : https://maps.googleapis.com/maps/api/staticmap?center=3.43,-76.5&zoom=12&size=640x640&scale=2&maptype=terrain&language=en-EN&key=xxx

mapa <- ggmap(city)

#mapa <- mapa + geom_point(data=paradasMio, aes(x=DECIMALLONGITUDE, y=DECIMALLATITUDE),color=paradasMio$TYPE_BUS,size=paradasMio$PASSANGERS_NUM * 0.01)

mapa <- mapa + geom_point(data=paradasMio, aes(x=DECIMALLONGITUDE, y=DECIMALLATITUDE),color=paradasMio$grupos,size=paradasMio$PASSANGERS_NUM * 0.01)

mapa <- mapa + geom_polygon(data=poligono, aes(x=DECIMALLONGITUDE, y=DECIMALLATITUDE), color="black", alpha=0.2 )

mapa