[R-sig-Geo] distances values from a point exhibit a hexagon-like distribution
Michael Sumner
mdsumner at gmail.com
Tue Jun 22 09:55:57 CEST 2010
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
>
More information about the R-sig-Geo
mailing list