Problems with mask (raster)
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