Could you please provide a reproducible example with useable code? Cheers, Mike.
On Tue, Jun 22, 2010 at 5:36 PM, Roman Luštrik <roman.lust...@gmail.com> wrote: > Can someone explain to me why values of distances from a point are not in a > uniform, circle-like fashion, but rather form a hexagon shape (see > image<http://imagepaste.nullnetwork.net/viewimage.php?id=1075>)? > Here is the code I'm using to produce this plot (kudos to > Robert<http://r-sig-geo.2731867.n2.nabble.com/assigning-raster-cell-values-based-on-predefined-criteria-td4929775.html#a4929775>). > Raster projection is "+projs = NA". > > rst.poly <- polygonsToRaster(circle.polygon, rst) #circle.polygon is the > polygon I'm overlaying, rst is the raster > rst.lines <- linesToRaster(circle.polygon, rst) > rst.lines[rst.lines > 0] <- 0 > rst.point <- pointsToRaster(rst, xy[[i]][j,]) #xy is the point from where > the distances are calculated > rst.cover <- cover(rst.point, rst.lines) > rst.grid <- gridDistance(rst.cover) > rst.grid[!is.na(rst.lines)] <- NA > plot(rst.grid) > > > Cheers, > Roman > > > -- > In God we trust, all others bring data. > > [[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 > _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo
