Se você tiver as latitudes e longitudes de todas as estações, pode colocá-las num mapa com o comando "points". Veja se te ajuda:

library(spdep)
library(akima) 
library(maps) 
library(geoR)

BR<-readShapePoly("BRASIL.shp")
CE<-BR[BR$UF == "CE",]

dados<-read.csv2("Ceará/estações.csv", h=T)

plot(CE, axes=T)
title("Ceará e Estações Pluviométricas")
points(dados$Long, dados$Lat,cex=0.75, pch=19)


--
YNWA