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