def create_elev_circles(): from renderer.processing import MapRenderer BASENAME = "64883885" HOMEAREA = "./data/" + BASENAME gf = GridFloat(HOMEAREA, BASENAME) mr = MapRenderer("./renderer/application.linux/renderer") mr.start(gf.left, gf.bottom, gf.right, gf.top, 2000) mr.smooth() mr.fill(250, 230, 230) mr.background(255, 255, 255) mr.strokeWeight(0.00007) for y in floatrange(gf.bottom, gf.top, (gf.top - gf.bottom) / 50): for x in floatrange(gf.left, gf.right, (gf.right - gf.left) / 50): elev = gf.elevation(x, y) mr.ellipse(x, y, elev * 0.00001, elev * 0.00001) mr.saveLocal("elevs.png") mr.stop()
def create_elev_circles(): from renderer.processing import MapRenderer BASENAME = "64883885" HOMEAREA = "./data/"+BASENAME gf = GridFloat(HOMEAREA, BASENAME) mr = MapRenderer("./renderer/application.linux/renderer") mr.start( gf.left, gf.bottom, gf.right, gf.top, 2000 ) mr.smooth() mr.fill(250,230,230) mr.background(255,255,255) mr.strokeWeight(0.00007) for y in floatrange( gf.bottom, gf.top, (gf.top-gf.bottom)/50 ): for x in floatrange( gf.left, gf.right, (gf.right-gf.left)/50 ): elev = gf.elevation( x, y ) mr.ellipse( x, y, elev*0.00001, elev*0.00001 ) mr.saveLocal( "elevs.png" ) mr.stop()
for i in range( step, cutoff, step ): ret.append( sg.contour( i, closure_tolerance=0.05) ) ret.append( sg.contour( cutoff, closure_tolerance=0.05 ) ) return ret def travel_time_surface(points, cutoff=30*60, fudge=1.1, margin=2, closure_tolerance=0.05, cellsize=0.004): return points_to_surface_grid(points, cutoff, fudge, margin, closure_tolerance, cellsize).to_matrix() if __name__=='__main__': points = Times("portland.times").times l, b, r, t = bounds(points) ttc = travel_time_contour( points ) from renderer.processing import MapRenderer mr = MapRenderer("./renderer/application.linux/renderer") mr.start(l, b, r, t, 1000) mr.smooth() mr.background( 255, 255, 255 ) for c in ttc: coords = list(c) for (x1,y1), (x2,y2) in cons(coords): mr.line(x1,y1,x2,y2) mr.saveLocal( "contour.png" ) mr.stop()