Wau! Muchachos han sido de gran ayuida. Sirvió! Seguiré el el camino y espero devolver la mano! Muchas gracias.
El 21 de julio de 2016, 5:56, Francisco Rodriguez Sanchez < f.rodriguez.s...@gmail.com> escribió: > 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 > -- *Mauricio Mardones Inostroza* Investigador Departamento Evaluación de Recursos Instituto de Fomento Pesquero - IFOP Valparaíso - Chile +56-32-21514 <callto:+56-32-2151424>42 www.ifop.cl [[alternative HTML version deleted]] _______________________________________________ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es