Hola Mauricio,
Como ya han indicado otros, spTransform y CRS sólo sirven para
reproyectar (de geográficas a UTM en este caso).
Para recortar a una zona determinada puedes usar por ejemplo la función
crop del paquete raster (entre otras opciones).
Algo así funciona:
chile <- readRDS("CHL_adm0.rds")
# project to UTM 19S
chile.UTM <- spTransform(chile,
CRS("+proj=utm +zone=19 +south
+ellps=WGS84 +datum=WGS84 +units=m +no_defs"))
bbox(chile.UTM) # bounding box (extent) of the entire country
# Now crop to selected region
library(raster)
area.extent <- extent(c(100000, 500000, 4500000, 5000000)) # define
limits (in UTM coordinates)
study.area <- crop(chile.UTM, area.extent) # crop
plot(study.area)
Suerte
Paco
El 20/07/2016 a las 18:43, Mauricio Mardones Inostroza escribió:
Hola a todos
Esta es mi primera pregunta en el grupo, y es sencilla pero me tiene
atascado. Estoy tratando de cortar mi mapa de (poner limites en UTM) en un
lugar definido como mi area de estudio (en este caso el sur de chile). Pero
creo no estar usando bien la función CRS ponendo bien los limites
requeridos.
study_area <- readRDS("CHL_adm0.rds")
study_area_UTM <- spTransform(study_area, CRS("+proj=utm +zone=19
+datum=WGS84"))
study_area_UTM <- spTransform(study_area_UTM, CRS(
+ paste("+x_0=-2000000.0 +y_0=-500000.0 +ellps=GRS80 +units=us-ft
+no_defs")))
Error in spTransform(study_area_UTM, CRS(paste("+x_0=-2000000.0
+y_0=-500000.0 +ellps=GRS80 +units=us-ft +no_defs"))) :
error in evaluating the argument 'CRSobj' in selecting a method for
function 'spTransform': Error in CRS(paste("+x_0=-2000000.0 +y_0=-500000.0
+ellps=GRS80 +units=us-ft +no_defs")) :
projection not named
Espero me puedan ayudar.
Saludos
--
Dr Francisco Rodriguez-Sanchez
Integrative Ecology Group
Estacion Biologica de Doñana - CSIC
Avda. Americo Vespucio s/n
41092 Sevilla (Spain)
http://bit.ly/frod_san
_______________________________________________
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
