def genBilGLUCInputs():
grasssetup.grassConfig()
grass = grasssetup.grass
grass.run_command('r.out.gdal', input='boundary',
output='gluc/Data/boundary.bil', format='EHdr', type='Byte')
grass.run_command('r.out.gdal', input='landcover',
output='gluc/Data/landcover.bil', format='EHdr', type='Byte')
grass.run_command('r.out.gdal', input='nogrowth',
output='gluc/Data/nogrowth.bil', format='EHdr', type='Byte')
grass.run_command('r.out.gdal', input='probmap_res',
output='gluc/Data/growth_probmap_res.bil', format='EHdr', type='Float32')
grass.run_command('r.out.gdal', input='probmap_com',
output='gluc/Data/growth_probmap_com.bil', format='EHdr', type='Float32')
#runlog.p('***use GLUC default pop_density')
#runlog.p('***use GLUC default emp_density')
if grass.find_file('pop_density', element='cell'):
grass.run_command('r.out.gdal', input='pop_density',
output='gluc/Data/pop_density.bil', format='EHdr',type='Float32') # needs to be Float32
runlog.p('***use uploaded pop_density in GLUC')
else:
runlog.p('***use GLUC default pop_density')
if grass.find_file('emp_density', element='cell'):
grass.run_command('r.out.gdal', input='emp_density',
output='gluc/Data/emp_density.bil', format='EHdr', type='Float32') # needs to be Float32
runlog.p('***use uploaded emp_density in GLUC')
else:
runlog.p('***use GLUC default emp_density')
def genBilGLUCInputs():
grasssetup.grassConfig()
grass = grasssetup.grass
grass.run_command('r.out.gdal', input='boundary',
output='gluc/Data/boundary.bil', format='EHdr', type='Byte')
grass.run_command('r.out.gdal', input='landcover',
output='gluc/Data/landcover.bil', format='EHdr', type='Byte')
grass.run_command('r.out.gdal', input='nogrowth',
output='gluc/Data/nogrowth.bil', format='EHdr', type='Byte')
grass.run_command('r.out.gdal', input='probmap_res',
output='gluc/Data/growth_probmap_res.bil', format='EHdr', type='Float32')
grass.run_command('r.out.gdal', input='probmap_com',
output='gluc/Data/growth_probmap_com.bil', format='EHdr', type='Float32')
if grass.find_file('pop_density', element='cell'):
grass.run_command('r.out.gdal', input='pop_density',
output='gluc/Data/pop_density.bil', format='EHdr',type='Float32') # needs to be Float32
runlog.p('***use uploaded pop_density in GLUC')
else:
runlog.p('***use GLUC default pop_density')
if grass.find_file('emp_density', element='cell'):
grass.run_command('r.out.gdal', input='emp_density',
output='gluc/Data/emp_density.bil', format='EHdr', type='Float32') # needs to be Float32
runlog.p('***use uploaded emp_density in GLUC')
else:
runlog.p('***use GLUC default emp_density')
def main():
# The system stdout and stderr will be in /leam/scratch/<scenarioname>/<scenarioname>.log
jobstart = time.asctime()
os.environ['PATH'] = ':'.join(['./bin', os.environ.get('BIN', '.'),'/bin', '/usr/bin'])
configurl, user, password = parseFlags()
print "Parsing configuration file.............\n"
configfile = get_config(configurl, user, password)
luc = LUC(configfile)
print luc.growthmap[0]
# luc has sorted luc.growthmap(drivers) in year and luc.growth(projections)
print "Setting up GRASS environment.............\n"
get_grassfolder(luc.scenario['grass_loc'], user, password)
grasssetup.grassConfig()
global grass
grass = grasssetup.grass
print "Connect to the LEAM storage server............."
resultsdir = luc.scenario['results']
print resultsdir, '\n'
title = luc.scenario['title']
global site, runlog # run.log will be stored in Log repository
site = LEAMsite(resultsdir, user=user, passwd=password)
#create runlog file
site.createFolder("RunLog", resultsdir)
runlog = RunLog(resultsdir+"/runlog", site, initmsg='Scenario ' + title)
runlog.p('started at '+jobstart)
# creat two folders
site.createFolder("Details", resultsdir)
site.createFolder("Results", resultsdir)
global projectiontable
projectiontable = ProjTable()
growth = dict(deltapop=[0.0], deltaemp=[0.0])
if luc.growth: # note growthmap[0]...should change luc, using luc_new
runlog.h('Processing Growth driver set.............')
startyear, isprobmapcached = processDriverSets(luc.growthmap[0])
if not isprobmapcached:
runlog.h('Building Probability Maps..............')
runMulticostModel(resultsdir, site, runlog)
cacheProbmaps(luc.growthmap[0]['url'])
if luc.growthmap:
runlog.h('Processing Growth Projection set........')
demandstr = processProjectionSet(luc.growth[0])
genYearChangemap.executeGLUCModel(demandstr, title, runlog)
runlog.h('Publishing all results..............')
#wrap the simmap file to Results
publishResults(title, site, resultsdir+"/results")
def main():
# The system stdout and stderr will be in /leam/scratch/<scenarioname>/<scenarioname>.log
jobstart = time.asctime()
os.environ['PATH'] = ':'.join(
['./bin', os.environ.get('BIN', '.'), '/bin', '/usr/bin'])
configurl, user, password = parseFlags()
print "Parsing configuration file.............\n"
configfile = get_config(configurl, user, password)
luc = LUC(configfile)
print luc.growthmap[0]
# luc has sorted luc.growthmap(drivers) in year and luc.growth(projections)
print "Setting up GRASS environment.............\n"
get_grassfolder(luc.scenario['grass_loc'], user, password)
grasssetup.grassConfig()
global grass
grass = grasssetup.grass
print "Connect to the LEAM storage server............."
resultsdir = luc.scenario['results']
print resultsdir, '\n'
title = luc.scenario['title']
global site, runlog # run.log will be stored in Log repository
site = LEAMsite(resultsdir, user=user, passwd=password)
#create runlog file
site.createFolder("RunLog", resultsdir)
runlog = RunLog(resultsdir + "/runlog", site, initmsg='Scenario ' + title)
runlog.p('started at ' + jobstart)
# creat two folders
site.createFolder("Details", resultsdir)
site.createFolder("Results", resultsdir)
global projectiontable
projectiontable = ProjTable()
growth = dict(deltapop=[0.0], deltaemp=[0.0])
if luc.growth: # note growthmap[0]...should change luc, using luc_new
runlog.h('Processing Growth driver set.............')
startyear, isprobmapcached = processDriverSets(luc.growthmap[0])
if not isprobmapcached:
runlog.h('Building Probability Maps..............')
runMulticostModel(resultsdir, site, runlog)
#open cache function temporarily
cacheProbmaps(luc.growthmap[0]['url'])
if luc.growthmap:
runlog.h('Processing Growth Projection set........')
demandstr = processProjectionSet(luc.growth[0])
genYearChangemap.executeGLUCModel(demandstr, title, runlog)
runlog.h('Publishing all results..............')
#wrap the simmap file to Results
publishResults(title, site, resultsdir + "/results")
def writeConfig(confname='baseline.conf',
start='2010',
end='2040',
prefix='',
probmaps=False,
tmpl='gluc/Config/baseline.tmpl',
path='gluc/Config'):
""" write a GLUC config file from a template.
PREFIX,START_DATE, and END_DATE are replaced in the template. PREFIX
is applied to probability maps name but the START_DATE is not. The
GLUC model automatically appends the year during model execution.
@confname - name of config file to be written (.conf appended if needed)
@prefix - prefix applied to probmap names (growth, decline)
@start - fills the START_DATE parameter in the template
@end - fills the END_DATE parameter in the template
@probmaps - toggle initial and final probmaps writing
@tmpl - name of the temlate config file
@path - destination for
@Return - nothing
# function copied from original startup.py
"""
# runlog.debug('writeConfig: confname=%s, prefix=%s, start=%s, end=%s' % \
# (confname, prefix, start, end))
# cleanup confname and add path
grasssetup.grassConfig()
grass = grasssetup.grass
if not confname.endswith('.conf'):
confname += '.conf'
config = os.path.join(path, os.path.basename(confname))
# read the template file into string
with open(tmpl, 'r') as f:
template = f.read()
d = dict(PREFIX=prefix, START_DATE=start, END_DATE=end)
with open(config, 'w') as f:
f.write(template.format(**d))
if probmaps:
f.write('* INITIAL_PROB_RES_MAP M(M, 4, initial_probmap_res)\n')
f.write('* INITIAL_PROB_COM_MAP M(M, 4, initial_probmap_com)\n')
f.write('* FINAL_PROB_RES_MAP M(M, 4, final_probmap_res)\n')
f.write('* FINAL_PROB_COM_MAP M(M, 4, final_probmap_com)\n')
if grass.find_file('pop_density', element='cell'):
f.write('* DENSITY_MAP_RES d(M, pop_density.bil)\n')
if grass.find_file('emp_density', element='cell'):
f.write('* DENSITY_MAP_COM d(M, emp_density.bil)\n')
def writeConfig(confname='baseline.conf',
start='2010', end='2040',
prefix='', probmaps=False,
tmpl='gluc/Config/baseline.tmpl',
path='gluc/Config'):
""" write a GLUC config file from a template.
PREFIX,START_DATE, and END_DATE are replaced in the template. PREFIX
is applied to probability maps name but the START_DATE is not. The
GLUC model automatically appends the year during model execution.
@confname - name of config file to be written (.conf appended if needed)
@prefix - prefix applied to probmap names (growth, decline)
@start - fills the START_DATE parameter in the template
@end - fills the END_DATE parameter in the template
@probmaps - toggle initial and final probmaps writing
@tmpl - name of the temlate config file
@path - destination for
@Return - nothing
# function copied from original startup.py
"""
# runlog.debug('writeConfig: confname=%s, prefix=%s, start=%s, end=%s' % \
# (confname, prefix, start, end))
# cleanup confname and add path
grasssetup.grassConfig()
grass = grasssetup.grass
if not confname.endswith('.conf'):
confname += '.conf'
config = os.path.join(path, os.path.basename(confname))
# read the template file into string
with open(tmpl, 'r') as f:
template = f.read()
d = dict(PREFIX=prefix, START_DATE=start, END_DATE=end)
with open(config, 'w') as f:
f.write(template.format(**d))
if probmaps:
f.write('* INITIAL_PROB_RES_MAP M(M, 4, initial_probmap_res)\n')
f.write('* INITIAL_PROB_COM_MAP M(M, 4, initial_probmap_com)\n')
f.write('* FINAL_PROB_RES_MAP M(M, 4, final_probmap_res)\n')
f.write('* FINAL_PROB_COM_MAP M(M, 4, final_probmap_com)\n')
if grass.find_file('pop_density', element='cell'):
f.write('* DENSITY_MAP_RES d(M, pop_density.bil)\n')
if grass.find_file('emp_density', element='cell'):
f.write('* DENSITY_MAP_COM d(M, emp_density.bil)\n')