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
