I think it should be: res <- mask(bio1cr, crmask)
On Thu, Oct 14, 2010 at 1:22 PM, Manuel Spínola <[email protected]> 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 > [email protected] > [email protected] > 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 > [email protected] > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > > _______________________________________________ R-sig-Geo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
