############################################
### CARGAMOS ESTOS PAQUETES
############################################
library(rgdal)
library(sp)
library(maptools)
library(raster)
library(scales)
library(foreign)
##############################################################################
# R te permite cargar información desde un URL
estacionesCNA<- read.csv(url("https://alanphd.com/archivos/estacionesLatLong.csv"))
# vemos el encabezado de la tabla
head(estacionesCNA)
# variables
names(estacionesCNA)
########### No hay metadatos, entonces NO podemos saber que nos dice cada # variable. Tampoco podemos saber el Datum de las coordenadas, asi que
### habrá que hacer una adivinanza educada.
# ¿Como podemos hacer un shapefile de puntos de esto?
# definimos X....Y
xy <- cbind(estacionesCNA$long,estacionesCNA$lat)
estaciones_geo<- SpatialPointsDataFrame(xy, estacionesCNA)
proj4string(estaciones_geo) <- CRS("+proj=longlat")
WGS84_Datum =CRS("+proj=longlat +datum=WGS84")
# precisamos el datum
projection(estaciones_geo) <- WGS84_Datum
### lo visualizamos
plot(estaciones_geo, axes = TRUE, pch=19)
##############################################################################
#### exportar como CSV, DEFINE TU DIRECTORIO DE TRABAJO
setwd('/Users/COLSON/Documents')
# write out a new shapefile (including .prj component)
writeOGR(estaciones_geo, ".", "estaciones_CNA_GCS_WGS84", driver="ESRI Shapefile")
##############################################################################
