r - Plot series of filled hexagons in ggplot2 -


i trying create tesselation of filled hexagons (polygons centered around hexagonally-spaced lattice) in ggplot2. have accomplished using 'plot' command struggling transitioning ggplot.

here code set-up:

# generate lattice of points equally spaced in centers of hexagonal lattice dist = 1     # distance between centers of hexagons nx = dist*15 # horizontal extent ny = dist*15 # vertical extent        makehexlattice = function(nx, ny, dist, origin=c(0,0)) {   locations = cbind(location = 1:(nx*ny),                     x = sort(c(rep(seq(from=0, by=dist, length.out=nx),each=ceiling(ny/2)),                                rep(seq(from=dist/2, by=dist, length.out=nx),                                    each=floor(ny/2)))) + origin[1],                     y = rep(c(seq(from=0, = dist*sqrt(3), length.out=ceiling(ny/2)),                               seq(from=dist*sqrt(3)/2, by=dist*sqrt(3),                                   length.out=floor(ny/2))) + origin[2], times=nx))   class(locations) = c(class(locations), "lattice")   attr(locations, "gridsize") = dist   return(locations) }

here code creating image using 'plot', looks nice: 

landscape = makehexlattice(nx=nx,ny=ny,dist=dist,origin=c(0,0))  # plot hexagonal lattice points plot(x=landscape[,2],y=landscape[,3], pch=19, col="black", cex=0.5, asp=1/1)  # separate x , y coordinates lx = landscape[,2] # x-coordinates ly = landscape[,3] # y-coordinates        # plot hexagonal lattice filled hexagons hex.x = cbind(lx + 0, lx + 0.5, lx + 0.5, lx + 0, lx - 0.5, lx - 0.5)  hex.y = cbind(ly - 1/(sqrt(3)), ly - 1/(2*sqrt(3)), ly + 1/(2*sqrt(3)), ly + 1/(sqrt(3)), ly + 1/(2*sqrt(3)), ly - 1/(2*sqrt(3))) hex.vectors = cbind(hex.x, hex.y)  for(i in 1:(length(hex.vectors)/12)){   polygon(x=hex.vectors[i,1:6], y=hex.vectors[i,7:12], angle = 120, border=null, col="wheat",            lty = par("lty"), filloddeven = false) }

any tips on how accomplish same thing using ggplot2 (which transitioning using)? have tried using geom_polygon can't seem work out for-loop. (also, please don't tell me use 'hexbin' -- not goal trying accomplish!)

thank help!

as things in ggplot, plotting extremely straightforward, of work getting data in right shape makes sense. for loop entirely unnecessary, geom_polygon() needs dataframe x , y coordinates, , variable defining group belong to. data:

library(ggplot2) library(reshape2)  #get coordinates in long format id hexdat.x <- melt(cbind(id = 1:length(hex.x), as.data.frame(hex.x)), id.vars = "id", value.name = "x") hexdat.y <- melt(cbind(id = 1:length(hex.y), as.data.frame(hex.y)), id.vars = "id", value.name = "y")  #merge them same dataframe hexdat <- merge(hexdat.x, hexdat.y)  head(hexdat) #   id variable    x          y # 1  1       v1  0.0 -0.5773503 # 2  1       v2  0.5 -0.2886751 # 3  1       v3  0.5  0.2886751 # 4  1       v4  0.0  0.5773503 # 5  1       v5 -0.5  0.2886751 # 6  1       v6 -0.5 -0.2886751

now plot hexagons, need give ggplot x , y coordinates, , specify group each 1 belongs to:

ggplot(hexdat, aes(x, y)) +    geom_polygon(aes(group = id), fill = "wheat", colour = "black")

enter image description here


Comments

Post a Comment

Popular posts from this blog

commonjs - How to write a typescript definition file for a node module that exports a function? -

consider, toml node module can use: // toml.d.ts declare module toml { export function parse(value:string):any; } declare module "toml" { export = toml; } then: /// <reference path="../../../../../defs/toml/toml.d.ts"/> import toml = require('toml'); toml.parse(...); however, node module export single function, such 'glob' ( https://github.com/isaacs/node-glob ). the node usage of module be: var glob = require("glob") glob("*.js", {}, callback(err, files) { ... }); you'd naively expect this: // glob.d.ts declare function globs(paths:string, options:any, callback:{(err:any, files:string[]):void; ...but because typescripts 'import' semantics bit strange, seems can use 'import .. = require()' statement alias modules . trying call: /// <reference path="../../../../../defs/glob/glob.d.ts"/> import blog = require('glob'); results in: error ts2072: mod
Read more

openid - Okta: Failed to get authorization code through API call -

i'm integrating okta own idp server using okta's api. i'm implementing authorization code flow following steps below: in own server, use /api/v1/authn endpoint sessiontoken. use sessiontoken obtain authorization calling endpoint: /oauth2/v1/authorize?client_id=" + clientid + "&sessiontoken=" + sessiontoken + "&response_type=code&response_mode=query&scope=openid&redirect_uri=" + redirecturl + "&state=evanyang&nonce=" it's supposed return response status code 302 , location header containing redirect url code value. however, keep getting response status code 200 , without location header, html body saying "you using unsupported browser." , "javascript disabled on browser." according api documentation: http://developer.okta.com/docs/api/resources/oidc.html#authentication-request , sessiontoken parameter sufficient this: an okta one-time sessiontoken. allows api-based
Read more

thorough guide for profiling racket code -

did googling ("racket profiling", "racket measure performance"), didn't find , there no examples in docs. did google "profile" search on htdp - no luck. (profile (f ...)) output isn't obvious other small snippets. ideally want python's python -m cprofile usage examples. when search "python profiling", both duckduckgo , google give top result: 26.4. python profilers . (although scanned quickly, seems more of reference "thorough usage guide" examples. if had else in mind, perhaps please link that?) the equivalent racket documentation be: profile: statistical profiler . a usage example: #lang racket (require profile) (profile-thunk (thunk (function-to-profile arg0 arg1) )) here (thunk e) convenience (lambda () e) . bigger example: #lang racket (module mod racket (provide f) (define (f) (for/list ([i 10000]) i))) (require (prefix-in mod: 'mod)) (define (f) (for ([i 10000
Read more