[R-sig-Geo] Problems with mask (raster)

Robert J. Hijmans r.hijmans at gmail.com
Thu Oct 14 22:44:28 CEST 2010


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
>
>



More information about the R-sig-Geo mailing list