Hi

I spend the whole day with it, trying grass, gdal, R, but was not successful
(and I know it can be done...).

I have a shape file and a raster file, both in the same projection and in
thew same area (shape file completely overlapping the raster) and I want to
"burn" the values from one column in the attribute table of the shape file
(DNSTY) into the raster.

I am doing the following:

####
library(sp)
library(maptools)
library(raster) # <- this might work as well. I also have the shape file as
a .csv file, including coordinates and DNSTY

s <- readShapePoints("./plants.shp")
r <- readAsciiGrid("./plants.asc")

## and now?
####

I seem to remember that it works by using overlay(), but I don't remember
how.\
Also, it should be working with pointsToRaster() from the package raster,
but I also was not successful.

Cheers,

Rainer

-- 
Rainer M. Krug, Centre of Excellence for Invasion Biology, Stellenbosch
University, South Africa

	[[alternative HTML version deleted]]

