Lectura de datos

Se cargan los datos y se crea la matriz de coordenadas.

ColCity<-read.csv2('ciudades_de_colombia.csv')
ColCityMat<-cbind(ColCity$long,ColCity$lat)
row.names(ColCityMat)<-ColCity$Nombre
head(ColCityMat)
##                      [,1]      [,2]
## Aguja, C. de la -74.20000 11.300000
## Antioquia       -75.91667  6.666667
## Apaporis        -69.41667 -1.383333
## Arauca          -70.66667  7.000000
## Arica           -71.83333 -2.000000
## Arjona          -75.36667 10.233333

Ahora se carga la librería y se construyen los objetos de la clase espaciales.

suppressPackageStartupMessages(require(sp))
SysProj<-CRS("+proj=longlat +ellps=WGS84")
ColCitySP<-SpatialPoints(ColCityMat,proj4string = SysProj)
summary(ColCitySP)
## Object of class SpatialPoints
## Coordinates:
##                 min       max
## coords.x1 -81.58333 -67.11667
## coords.x2  -4.15000  13.53333
## Is projected: FALSE 
## proj4string : [+proj=longlat +ellps=WGS84]
## Number of points: 90
bbox(ColCitySP)
##                 min       max
## coords.x1 -81.58333 -67.11667
## coords.x2  -4.15000  13.53333

Lectura del mapa en formato “Shape”

Para cargar los croquis de los municipios de Colombia, se utilizará la librería maptools.

suppressPackageStartupMessages(require(maptools))
ColShp<-readShapeSpatial('COL_adm2.shp')
bbox(ColShp)
##          min       max
## x -81.841530 -66.87033
## y  -4.228429  15.91247
proj4string(ColShp) <- CRS("+proj=longlat +ellps=WGS84")

par(mar=c(.5,.5,.5,.5))
plot(ColShp,ylim=c(-4,13),col='honeydew')
points(ColCitySP,col='blue',pch=19,cex=.7)