def makeStateCounties(): db = pg.Database( database = 'usageo_500k' ) for level in ( '512', ): table = 'county_' + level fulltable = schema + '.' + table shpname = 'us2012-county-500k-%s' % level shpfile = '%(path)s/us2012/%(shpname)s/%(shpname)s.shp' %({ 'path': private.OUTPUT_SHAPEFILE_PATH, 'shpname': shpname, }) db.dropTable( fulltable ) db.loadShapefile( shpfile, private.TEMP_PATH, fulltable, googGeom, '3857', 'LATIN1', True ) #mergeStates( db, table, level ) #writeEachState( db, table, level ) writeEachStateCounties( db, fulltable, level ) #db.connection.commit() #db.connection.close() #level = '512' #if level is not None: # addLevel( db, table, level ) #mergeCounties( db, level ) #writeEachState( db, table, level ) #writeAllStates( db, table, level ) db.connection.commit() db.connection.close()
def makeDepartments(): db = pg.Database(database='france2012') table = 'departement' #for level in ( None, 512, 1024, 2048, 4096, 8192, ): for level in ( None, 4096, ): if level is not None: addLevel(db, table, 'code_dept', level) mergeGeometries( db, 'france.departement', 'france.reg2012', level, ''' WHERE france.departement.code_reg = france.reg2012.region AND substring(france.departement.code_reg from 1 for 1) != '0' GROUP BY france.departement.code_reg ''') mergeGeometries(db, 'france.reg2012', 'france.france2012', level, ''' WHERE true ''') #mergeGeometries( db, table, level ) #writeEachDepartment( db, table, level ) writeAllDepartments(db, table, level) db.connection.commit() db.connection.close()
def main(): db = pg.Database(database=database) loadNameTables(db) loadContinents(db) makeNation(db) db.connection.commit() db.connection.close()
def makeLegislative(): db = pg.Database(database='france2012') table = 'legislative' #for level in ( None, 512, 1024, ): for level in (1024, ): if level is not None: addLevel(db, table, 'id_circo', level) mergeGeometries( db, 'france.legislative', 'france.departement', level, ''' WHERE france.legislative.id_dep = france.departement.code_dept GROUP BY france.legislative.id_dep ''') mergeGeometries( db, 'france.departement', 'france.reg2012', level, ''' WHERE france.departement.code_reg = france.reg2012.region AND substring(france.departement.code_reg from 1 for 1) != '0' GROUP BY france.departement.code_reg ''') mergeGeometries(db, 'france.reg2012', 'france.france2012', level, ''' WHERE true ''') #mergeGeometries( db, table, level ) #writeEachLegislative( db, table, level ) writeAllLegislative(db, table, level) db.connection.commit() db.connection.close()
def main(): global db db = pg.Database(host=opt.hostname, database=opt.database, user=opt.username, password=opt.password) process() db.connection.close()
def makeStatesNational(): db = pg.Database( database = 'usageo_500k' ) table = 'state' level = '4096' if level is not None: addLevel( db, table, level ) mergeStates( db, table, level ) #writeEachState( db, table, level ) writeAllStates( db, table, level ) db.connection.commit() db.connection.close()
def makeGopDetail(): db = pg.Database( database = 'usageo_500k' ) table = 'gop2012loc' level = '512' if level is not None: addLevel( db, table, level ) mergeStates( db, table, level ) writeEachState( db, table, level ) writeAllStates( db, table, level ) db.connection.commit() db.connection.close()
def makeMunis(): db = pg.Database( database = database ) createSchema( db, schema ) for level in ( #'full', '512', '1024', ): munitable = schema + '.muni' + level provtable = schema + '.prov' + level nltable = schema + '.nl' + level db.dropTable( munitable ) db.dropTable( provtable ) db.dropTable( nltable ) srcfile = 'nl2012-%s' % level filename = '../shapes/shp/%s/%s.shp' %( srcfile, srcfile ) print 'Loading %s' % filename db.loadShapefile( filename, private.TEMP_PATH, munitable, geom, '3857', 'LATIN1', True ) mergeGeometries( db, munitable, provtable, 'prov', ''' prov varchar(2) ''', ''' WHERE %(munitable)s.prov IS NOT NULL GROUP BY %(munitable)s.prov ''' % { 'munitable': munitable, }) mergeGeometries( db, munitable, nltable, "'NL'", ''' nation varchar(2) ''', ''' WHERE %(munitable)s.prov IS NOT NULL ''' % { 'munitable': munitable, }) writeAllMunis( db, munitable, level ) db.connection.commit() db.connection.close()
def makeLegimunes(): db = pg.Database(database='france2012') table = 'legimune' for level in (None, ): if level is not None: addLevel(db, table, 'code_leg,code_comm', level) mergeGeometries( db, 'france.legimune', 'france.legislative', level, ''' WHERE france.legimune.code_dept = france.legislative.id_dep AND france.legimune.code_leg = france.legislative.id_circo GROUP BY france.legimune.code_dept, france.legimune.code_leg ''') writeEachLegislative(db, table, None) db.connection.commit() db.connection.close()
def makeCommunes(): db = pg.Database(database='france2012') table = 'commune' for level in ( None, 128, ): if level is not None: addLevel(db, table, 'code_dept,code_comm', level) mergeGeometries( db, 'france.commune', 'france.departement', level, ''' WHERE france.commune.code_dept = france.departement.code_dept GROUP BY france.commune.code_dept ''') writeEachDepartment(db, table, None) db.connection.commit() db.connection.close()
def process(): #db = pg.Database( database='postgres' ) #db.createGeoDatabase( 'census' ) #db.connection.close() db = pg.Database(database='census') #for year in '2000', '2010': for year in '2010', : schema = 'c' + year db.createSchema(schema) for kind in 'state', 'county', 'tract', 'bg': table = schema + '.' + kind create = True for state in states.array: fips = state['fips'] print 'Loading', table, fips, state['name'] zipfile = censusFileName(kind, year, fips) db.loadShapeZip(zipfile, private.TEMP_PATH, table, create) create = False db.connection.close()
def makeHouseNational(): db = pg.Database( database = 'usageo_500k' ) #level = '4096' #for level in '512', '4096': for level in ( '2048', ): table = 'house2012_' + level fulltable = schema + '.' + table shpname = 'us2012-house2012-%s' % level shpfile = '%(path)s/shp/us2012/%(shpname)s/%(shpname)s.shp' %({ 'path': private.OUTPUT_SHAPEFILE_PATH, 'shpname': shpname, }) db.dropTable( fulltable ) db.loadShapefile( shpfile, private.TEMP_PATH, fulltable, googGeom, '3857', 'LATIN1', True #, None, tweakHouseSQL ) #mergeStates( db, table, level ) #writeEachState( db, table, level ) writeAllStatesHouse( db, fulltable, level )
def process(): #db = pg.Database( database='postgres' ) #db.createGeoDatabase( 'turkey' ) #db.connection.close() db = pg.Database(database='turkey') #db.createGeoDatabase( 'turkey' ) schema = 't2011' db.createSchema(schema) db.connection.commit() loadProvinceFT(db, schema) loadDistrictFT(db, schema) #for level in '00 50 60 70 80 90'.split(' '): for level in '90'.split(' '): loadSHP(db, schema, 'districts', level) for tbl in ( 'provinces', 'districts', ): #for tbl in ( 'provinces', ): name = tbl # 'Turkey' table = '%s.%s' % (schema, tbl) gid = '-1' targetGeom = 'geom_' + level #boxGeom = 'geom_00' boxGeom = targetGeom #db.addGoogleGeometry( table, geom, googGeom ) filename = '%s/turkey-%s-%s.jsonp' % (private.GEOJSON_PATH, tbl, targetGeom) db.makeGeoJSON(filename, table, boxGeom, targetGeom, tbl, name, gid, 'loadGeoJSON') db.connection.commit() db.connection.close()
def createDatabase( database ): print 'Creating database %s' % database db = pg.Database( database='postgres' ) db.createGeoDatabase( database ) db.connection.close()
def main(): global db db = pg.Database( database='turkey' ) process() db.connection.close() print 'Done!'
def main(): db = pg.Database(database=database) loadContinents(db) makeRegions(db) db.connection.commit() db.connection.close()
def openDatabase(): print 'Opening database %s' % database return pg.Database(database=database)
def main(): global db db = pg.Database(database='census') process() db.connection.close()