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()
