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 [[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo