[R-es] un solo un favor

Francisco Rodriguez Sanchez f.rodriguez.sanc en gmail.com
Jue Jul 21 11:56:12 CEST 2016


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



Más información sobre la lista de distribución R-help-es