def initVrt():
rasterpath = os.path.join(getEnv()['GISDBASE'], getEnv()['LOCATION_NAME'], getEnv()['MAPSET'], 'cellhd')
vrtdirR = os.path.join(getEnv()['GISDBASE'], getEnv()['LOCATION_NAME'], getEnv()['MAPSET'], 'vrt','raster/')
vrtdirV = os.path.join(getEnv()['GISDBASE'], getEnv()['LOCATION_NAME'], getEnv()['MAPSET'], 'vrt', 'vector/')
r = os.path.dirname(vrtdirR)
v = os.path.dirname(vrtdirV)
if not os.path.exists(r):
os.makedirs(r)
if not os.path.exists(v):
os.makedirs(v)
for i in RasterList():
mapfile = os.path.join(rasterpath, i)
vrtfilename = i + '.vrt'
output = os.path.join(vrtdirR, vrtfilename)
print output
try :
gdal.GetDriverByName('VRT').CreateCopy(output,gdal.Open(mapfile))
os.system('ossim-img2rr --create-histogram %s' % (output) )
except :
print 'vrt creation failed for map %s ' % (i)
for i in VectorList():
vectorpath = os.path.join(getEnv()['GISDBASE'], getEnv()['LOCATION_NAME'], getEnv()['MAPSET'], 'vector' )
mapfile = os.path.join(vectorpath, i, 'head')
vrtfilename = i + '.vrt'
outvrt = os.path.join(vrtdirV, vrtfilename)
try :
ogrvrt(mapfile,outvrt)
except :
print 'ogrvrt failed for map : %s' % (i)
def initVrt():
rasterpath = os.path.join(getEnv()['GISDBASE'],
getEnv()['LOCATION_NAME'],
getEnv()['MAPSET'], 'cellhd')
vrtdirR = os.path.join(getEnv()['GISDBASE'],
getEnv()['LOCATION_NAME'],
getEnv()['MAPSET'], 'vrt', 'raster/')
vrtdirV = os.path.join(getEnv()['GISDBASE'],
getEnv()['LOCATION_NAME'],
getEnv()['MAPSET'], 'vrt', 'vector/')
r = os.path.dirname(vrtdirR)
v = os.path.dirname(vrtdirV)
if not os.path.exists(r):
os.makedirs(r)
if not os.path.exists(v):
os.makedirs(v)
for i in RasterList():
mapfile = os.path.join(rasterpath, i)
vrtfilename = i + '.vrt'
output = os.path.join(vrtdirR, vrtfilename)
print output
try:
gdal.GetDriverByName('VRT').CreateCopy(output, gdal.Open(mapfile))
os.system('ossim-img2rr --create-histogram %s' % (output))
except:
print 'vrt creation failed for map %s ' % (i)
for i in VectorList():
vectorpath = os.path.join(getEnv()['GISDBASE'],
getEnv()['LOCATION_NAME'],
getEnv()['MAPSET'], 'vector')
mapfile = os.path.join(vectorpath, i, 'head')
vrtfilename = i + '.vrt'
outvrt = os.path.join(vrtdirV, vrtfilename)
try:
ogrvrt(mapfile, outvrt)
except:
print 'ogrvrt failed for map : %s' % (i)
def main():
add = flags['a']
remove = flags['r']
host = options['host']
dport = options['dport']
pport = options['pport']
grassenv = grass.gisenv()
mappa = options['map'].replace("@", " ")
mappa = mappa.split()
nflags = len(filter(None, [add, remove]))
if nflags > 1:
grass.run_command('g.message',
message='Cannot add & remove a map at same time.')
if nflags < 1:
grass.run_command(
'g.message',
message=
'No action requested , please choose one from "-a : add" or "-r : remove" flags.'
)
try:
vectorpath = os.path.join(grassenv['GISDBASE'],
grassenv['LOCATION_NAME'], mappa[1],
'vector', mappa[0])
mapfile = os.path.join(vectorpath, 'head')
vrtdir = os.path.join(grassenv['GISDBASE'], grassenv['LOCATION_NAME'],
mappa[1], 'vrt', 'vector/')
except:
vectorpath = os.path.join(grassenv['GISDBASE'],
grassenv['LOCATION_NAME'],
grassenv['MAPSET'], mappa[0])
mapfile = os.path.join(vectorpath, 'head')
vrtdir = os.path.join(grassenv['GISDBASE'], grassenv['LOCATION_NAME'],
grassenv['MAPSET'], 'vrt', 'vector/')
d = os.path.dirname(vrtdir)
if not os.path.exists(d):
os.makedirs(d)
vrtfilename = mappa[0] + '.vrt'
output = os.path.join(vrtdir, vrtfilename)
if not os.path.exists(output):
ogrvrt(mapfile, output)
print 'try to make omd'
makestile(output, options['brush'], options['pen'], options['size'],
options['fill'], options['thickness'])
proj_info = projinfo()
unit = proj_info['units']
if unit.lower() == 'meters':
zoom_position = setCPRJ(options['map'])
if unit.lower() == 'metres':
zoom_position = setCPRJ(options['map'])
if unit == 'Degrees':
zoom_position = setCLL(options['map'])
lat = zoom_position[0]
lon = zoom_position[1]
distance = zoom_position[2]
if add:
try:
addzoom(output, lon, lat, distance, host, dport, pport)
print 'Added vector file :', mappa[0]
print 'Camera positioned to : '
print 'Longitude = ', lon
print 'Latitude = ', lat
print 'Altitude = ', distance
except:
print "conecction error"
if remove:
removefile(output, host, dport)
print 'Removed vector file :', mappa[0]
def main():
add = flags['a']
remove = flags['r']
host = options['host']
dport = options['dport']
pport = options['pport']
grassenv = grass.gisenv()
mappa = options['map'].replace("@"," ")
mappa = mappa.split()
nflags = len(filter(None, [add, remove]))
if nflags > 1:
grass.run_command('g.message' , message = 'Cannot add & remove a map at same time.')
if nflags < 1:
grass.run_command('g.message' , message = 'No action requested , please choose one from "-a : add" or "-r : remove" flags.')
try :
vectorpath = os.path.join(grassenv['GISDBASE'], grassenv['LOCATION_NAME'], mappa[1], 'vector' , mappa[0] )
mapfile = os.path.join(vectorpath, 'head')
vrtdir = os.path.join(grassenv['GISDBASE'], grassenv['LOCATION_NAME'], mappa[1], 'vrt', 'vector/')
except :
vectorpath = os.path.join(grassenv['GISDBASE'], grassenv['LOCATION_NAME'], grassenv['MAPSET'], mappa[0] )
mapfile = os.path.join(vectorpath, 'head' )
vrtdir = os.path.join(grassenv['GISDBASE'], grassenv['LOCATION_NAME'], grassenv['MAPSET'], 'vrt', 'vector/')
d = os.path.dirname(vrtdir)
if not os.path.exists(d):
os.makedirs(d)
print d
vrtfilename = mappa[0] + '.vrt'
output = os.path.join(vrtdir, vrtfilename)
#file = open(output,'w')
#file.write('')
#file.close()
#print output
outfile = output.replace('.vrt','')
outfile = outfile+'.omd'
if os.path.isfile(outfile):
os.remove(outfile)
if not os.path.exists(output):
print 'ogvrt try'
#print mapfile, output
ogrvrt(mapfile,output)
print 'try to make omd'
makestile(output, options['brush'], options['pen'], options['size'], options['fill'], options['thickness'])
proj_info = projinfo()
unit = proj_info['units']
if unit.lower() == 'meters':
zoom_position = setCPRJ(options['map'])
if unit.lower() == 'metres':
zoom_position = setCPRJ(options['map'])
if unit == 'Degrees':
zoom_position = setCLL(options['map'])
lat = zoom_position[0]
lon = zoom_position[1]
distance = zoom_position[2]
print output,lon,lat,distance,host,dport,pport
if add :
try:
addzoom(output,lon,lat,distance,host,dport,pport)
print 'Added vector file :', mappa[0]
print 'Camera positioned to : '
print 'Longitude = ',lon
print 'Latitude = ', lat
print 'Altitude = ' , distance
except :
print "conecction error"
if remove :
removefile(output,host,dport)
print 'Removed vector file :', mappa[0]
