[R-es] un solo un favor

[rubenfcasal]

Hola Mauricio,

     No tengo muy claro que intentas hacer. En primer lugar la 
proyección/CRS no serviría para recortar los datos (simplemente sirve 
para poder interpretar las coordenadas), para esto en principio 
utilizaría la función 'over' del paquete sp (ejecuta 'vignette("over")').

     En segundo lugar tu objeto "CHL_adm0.rds" parece ser un polígono 
espacial correspondiente al contorno de Chile (descargado de 
http://www.gadm.org/country). Si simplemente lo vas a utilizar para 
representarlo bastaría con fijar los límites a la hora de hacer el 
gráfico (o pintarlo por encima...). Si necesitas algo más, recortarlo 
puede ser un problema. Puede que lo que te interese sea aumentar la 
resolución (empleando "CHL_admX.rds", i.e. áreas administrativas de 
nivel X) y seleccionar las áreas que te interesen. En cualquier caso te 
recomiendo que te documentes sobre como manipular este tipo de objetos 
(la referencia recomendada para empezar sería el libro Applied Spatial 
Data Analysis With R).

     Un saludo, Rubén.


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