if runtype not in ['raster', 'proximity']: print "type", runtype, "not recognized" sys.exit() startx = int(options.startx) starty = int(options.starty) place = MapGardening.get_place(options.place) if place is not None: print "starting with", options.place MapGardening.init_logging() MapGardening.init_db(place['dbname']) # analyze to nodes table, ignore others for now. nodetable = MapGardening.NodeTable() # TODO: test for spatial indexes before proceeding bstm = MapGardening.BlankSpotTableManager() # We assume the manager table exists. If not, create it. #if not bstm.table_exists(): # bstm.create_manager_table() params = { 'runtype': runtype,
usage = "usage: %prog [options]" p = optparse.OptionParser(usage) p.add_option('--place', '-p', default="all" ) options, arguments = p.parse_args() if options.place == "all": places = MapGardening.get_all_places() else: placename = options.place place = MapGardening.get_place(placename) places = {placename: place} MapGardening.init_logging() for placename in places.keys(): print "printing raster info for", placename MapGardening.init_db(places[placename]['dbname']) raster = MapGardening.Raster() raster.get_raster_stats() MapGardening.disconnect_db()
# TODO: query the BlankSpotTableManager to see which resolutions are present possible_resolutions = ["250", "500", "1000"] if options.resolution == "all": resolutions = possible_resolutions else: resolutions = [options.resolution] MapGardening.init_logging() for placename in places.keys(): # Each place exists in a different database, so we connect/disconnect inside the places loop (conn, cur) = MapGardening.init_db(places[placename]["dbname"]) bstm = MapGardening.BlankSpotTableManager() nt = MapGardening.NodeTable() # Using defaults for table name and projection for runtype in runtypes: for resolution in resolutions: st = time.time() lt = st print "starting", placename, resolution, runtype params = {"runtype": runtype, "resolution": resolution}