Dear list members,

I am trying to make a raster from a shapefile that I imported to r using
the function "readOGR".  When I use the function rasterize and try to plot
the object nothing is plotted.   Is this happening because the field is
categoric (factor) and not numeric?

> cr =
readOGR(dsn="/Users/manuelspinola/Documents/ProyectosRespacial/EvalHab_reporte",
layer="Cr_wgs84_meso")

> crproy = spTransform(cr, CRS("+proj=tmerc +lon_0=-84 +lat_0=0 +x_0=500000
+k=0.9999 +datum=WGS84"))

> lu =
readOGR(dsn="/Users/manuelspinola/Documents/ProyectosRespacial/Ecologia_del_paisaje",
layer="usosuelos2006_cr05")

> crr = raster(crproy)

> lur =  rasterize(lu, crr, field = "COB_AGROP")
Found 157480 region(s) and 227521 polygon(s)
Mensajes de aviso perdidos
1: In .polygonsToRaster(x, y, ...) : selected field is character type
2: In .polygonsToRaster(x, y, ...) : NAs introducidos por coerción

> lur
class       : RasterLayer
dimensions  : 10, 10, 100  (nrow, ncol, ncell)
resolution  : 37190.08, 35181.17  (x, y)
extent      : 286586.4, 658487.2, 889526.4, 1241338  (xmin, xmax, ymin,
ymax)
coord. ref. : +proj=tmerc +lon_0=-84 +lat_0=0 +x_0=500000 +k=0.9999
+datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
values      : in memory
min value   : NA
max value   : NA

plot(lur)

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@una.ac.cr
mspinola10@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]]

