Hi Juliane,
may I also suggest having a look at the pick.from.ascii.grid function
(and the like) in the RSAGA package. This will use nearest neighbour
interpolation (or, for reasonably small grids, may also use an ad hoc
local kriging interpolation) to pick grid values (from grids with
numeric data).
http://cran.r-project.org/web/packages/RSAGA/RSAGA.pdf
There are some more suggestions in the R-sig-geo archives:
https://stat.ethz.ch/pipermail/r-sig-geo/2008-October/004469.html
http://www.mail-archive.com/r-sig-geo@stat.math.ethz.ch/msg04662.html
I hope this helps
Alex
Struve, Juliane wrote:
Dear list,
I am trying to incoporate spatial information (for example distance)
into a movement model by reading an ascii -b converted raster and
matching the raster with points of a movement path (for example
fishmove$coord[fishmove$day==1]. Below is a section of my code to
explain the problem. This code does not work because the coordinates
of the movement path are not the same as those in the raster. Rounding
the coordinates does not seem to solve the problem. Maybe the closest
point would have to be found, but i am not sure how to approach this.
Is there another perhaps easier way ?
Dist_DF=data.frame(readAsciiGrid("distance.txt"))
Dist_shoreDF$coord = paste(round(Dist_DF$s1),round(Dist_DF$s2),sep=",")
names(Dist_DF)[1]="Dist"
Dist_DF$Dist[Dist_shoreDF$coord==fishmove$coord[fishmove$day==1]]
I would be grateful for any help.
Many thanks,
Juliane
--
Alexander Brenning
brenn...@uwaterloo.ca - T +1-519-888-4567 ext 35783
Department of Geography and Environmental Management
University of Waterloo
200 University Ave. W - Waterloo, ON - Canada N2L 3G1
http://www.fes.uwaterloo.ca/geography/faculty/brenning/
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
