[R-es] un solo un favor

[Mauricio Mardones Inostroza]

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.sanc en 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 en 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]]