[R-es] un solo un favor
rubenfcasal
rubenfcasal en gmail.com
Jue Jul 21 11:04:53 CEST 2016
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
Más información sobre la lista de distribución R-help-es