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
