Re: [R-es] un solo un favor
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(10, 50, 450, 500)) # 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=-200.0 +y_0=-50.0 +ellps=GRS80 +units=us-ft
+no_defs")))
Error in spTransform(study_area_UTM, CRS(paste("+x_0=-200.0
+y_0=-50.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=-200.0 +y_0=-50.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
Re: [R-es] un solo un favor
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=-200.0 +y_0=-50.0 +ellps=GRS80 +units=us-ft
+no_defs")))
Error in spTransform(study_area_UTM, CRS(paste("+x_0=-200.0
+y_0=-50.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=-200.0 +y_0=-50.0
+ellps=GRS80 +units=us-ft +no_defs")) :
projection not named
Espero me puedan ayudar.
Saludos
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
Re: [R-es] un solo un favor
Hola Mauricio,
No entiendo 2 cosas en tu código:
1. Por qué proyectas el mismo objeto:
> study_area_UTM <- spTransform(study_area, CRS("+proj=utm +zone=19
+datum=WGS84"))
> study_area_UTM <- spTransform(study_area_UTM, CRS(
+ paste("+x_0=-200.0 +y_0=-50.0 +ellps=GRS80 +units=us-ft
+no_defs")))
2. Por qué usas "paste" (creo que no deberías usarlo).
Manuel Spínola
El 20 de julio de 2016, 10:43, Mauricio Mardones Inostroza <
mauricio.mardo...@ifop.cl> 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=-200.0 +y_0=-50.0 +ellps=GRS80 +units=us-ft
> +no_defs")))
> Error in spTransform(study_area_UTM, CRS(paste("+x_0=-200.0
> +y_0=-50.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=-200.0 +y_0=-50.0
> +ellps=GRS80 +units=us-ft +no_defs")) :
> projection not named
>
>
>
> Espero me puedan ayudar.
>
> Saludos
> --
>
> *Mauricio Mardones Inostroza*
>
> Investigador Departamento Evaluación de Recursos
> Instituto de Fomento Pesquero - IFOP
> Valparaíso - Chile
> +56-32-21514 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
>
--
*Manuel Spínola, Ph.D.*
Instituto Internacional en Conservación y Manejo de Vida Silvestre
Universidad Nacional
Apartado 1350-3000
Heredia
COSTA RICA
mspin...@una.cr
mspinol...@gmail.com
Teléfono: (506) 8706 - 4662
Personal website: Lobito de río
[R-es] un solo un favor
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=-200.0 +y_0=-50.0 +ellps=GRS80 +units=us-ft
+no_defs")))
Error in spTransform(study_area_UTM, CRS(paste("+x_0=-200.0
+y_0=-50.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=-200.0 +y_0=-50.0
+ellps=GRS80 +units=us-ft +no_defs")) :
projection not named
Espero me puedan ayudar.
Saludos
--
*Mauricio Mardones Inostroza*
Investigador Departamento Evaluación de Recursos
Instituto de Fomento Pesquero - IFOP
Valparaíso - Chile
+56-32-21514 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
