help="comma separated list of bands to download (all: B01,B02,B03,B04,B05,B06,B07,B08,B8A,B09,B10,B11,B12)", ) parser.add_argument( "--days", dest="days", default=DEFAULT_DAYS, help="how many days to look back for new images (default %s)" % DEFAULT_DAYS, ) parser.add_argument( "--rebuild", dest="rebuild", default=False, action="store_true", help="force downloading of images" ) parser.add_argument("--quiet", dest="quiet", default=False, action="store_true", help="suppress debug messages") args = parser.parse_args() # init utility library util = Util() download_images = {} entries = [] # set geometry for regions geom = {"atwest": "POLYGON((9.53076 46.6515, 9.53076 47.743, 12.9659 47.743, 12.9659 46.6515, 9.53076 46.6515))"} # bail out if no geometry is available if not args.region in geom: print "ERROR: No geometry to query found for region %s." % args.region # set query query = { "q": 'ingestionDate:[NOW-%sDAYS TO NOW] AND footprint:"Intersects(%s)"' % (args.days, geom[args.region]), "start": 0,
if __name__ == '__main__': """ create RGB and NDVI images from Sentinel-2 bands in JPEG2000 format """ import argparse parser = argparse.ArgumentParser(description='create RGB and NDVI images from Sentinel-2 bands in JPEG2000 format') parser.add_argument('--region', dest='region', help='set region abreviation (e.g. atwest)') parser.add_argument('--tiles', dest='tiles', help='comma separated list of tiles to process (e.g. 32TNT,32TPT)') parser.add_argument('--rgb', dest='rgb', action="store_true", help='create RGB-image') parser.add_argument('--ndvi', dest='ndvi', action="store_true", help='create NDVI-image') parser.add_argument('--rebuild', dest='rebuild', default=DEFAULT_REBUILD, action="store_true", help='force downloading of images') parser.add_argument('--quiet', dest='quiet', default=DEFAULT_QUIET, action="store_true", help='suppress debug messages') args = parser.parse_args() # init utility library util = Util() # find JPEG files to process jpegs = {} tiles = {} for dirpath, dirnames, filenames in os.walk(util.get_satdir()): if args.region and not re.search(args.region,dirpath): continue for fname in sorted(filenames): if fname[-3:] == 'jp2': # take tiles argument into account if args.tiles and not filename_matches_tile(fname,args.tiles): continue