Trabalhando com mapas no R

library(ggplot2)
library(ggmap)
library(maps)
library(mapdata)
library(dismo)
#Usando dados do pacote "mapdata"
usa <- map_data("usa") # Definindo qual dado queremos trabalhar
ggplot() + 
  geom_polygon(data = usa, aes(x=long, y = lat, group = group)) + # Aqui estamos usando o data frame "usa" para plotar-lo como poligono
  coord_fixed(1.3)

# Se quiser fazer o mesmo, mas com os limites dos estados:
states <- map_data("state")

ggplot(data = states) + 
  geom_polygon(aes(x = long, y = lat, group = group), color = "white") + 
  coord_fixed(1.3) +
  guides(fill=FALSE)
# para adicionar seus pontos (objeto chamado "species"), considerando que vc tem um dataframe com ao menos tres colunas (lon, lat e grupo resultate do K-means, que chaamrei de kmeans), ficaria algo assim:
mapsBase <- ggplot(data = states) + 
  geom_polygon(aes(x = long, y = lat, group = group), color = "white") + 
  coord_fixed(1.3) +
  guides(fill=FALSE)
mapBase + geom_point(data=species, aes(x=lon, y=lat),group = kmeans)