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]