Message-ID: <AANLkTik5H2f4t3nVR0LuVBgmaQrk-vanExGZQGWEcOT1@mail.gmail.com>
Date: 2010-10-14T20:44:28Z
From: Robert J. Hijmans
Subject: Problems with mask (raster)
In-Reply-To: <4CB7667A.1010200@gmail.com>
I think it should be:
res <- mask(bio1cr, crmask)
On Thu, Oct 14, 2010 at 1:22 PM, Manuel Sp?nola <mspinola10 at gmail.com> wrote:
> ?Dear list members,
>
> I am working with bioclim variables to model species distribution.
> The files downloaded from Worldclim.org are generic raster files (.bil)
> I transform this files to raster:
>
> bio1 <- raster("bio1_23.bil")
>
> cr = readOGR(dsn="C:/ProyectoRespacial/EvalHab_reporte",
> layer="Cr_wgs84_meso")
>
> bio1cr <- crop(bio1, cr)
> crmask <- polygonsToRaster(cr, bio1cr)
> res <- mask(bio1, crmask)
> Error en compare(x, mask) : Different bounding box
>
> plot(bio1cr) the scale legend shows numbers from 0 to 55000
>
> Best,
>
> Manuel
>
>
> --
> *Manuel Sp?nola, Ph.D.*
> Instituto Internacional en Conservaci?n y Manejo de Vida Silvestre
> Universidad Nacional
> Apartado 1350-3000
> Heredia
> COSTA RICA
> mspinola at una.ac.cr
> mspinola10 at gmail.com
> Tel?fono: (506) 2277-3598
> Fax: (506) 2237-7036
> Personal website: Lobito de r?o
> <https://sites.google.com/site/lobitoderio/>
> Institutional website: ICOMVIS <http://www.icomvis.una.ac.cr/>
>
> ? ? ? ?[[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
>