[R-sig-Geo] distances values from a point exhibit a hexagon-like distribution

[Michael Sumner]

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.lustrik at 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 at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>