# 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
granule = util.parse_granule(fname[:-4])
img_dir = re.sub(r'jpeg$','img',dirpath)
img_rgb = os.path.join(img_dir, "rgb_%s.img" % granule['tilename'])
img_ndvi = os.path.join(img_dir, "ndvi_%s.img" % granule['tilename'])
if args.rgb and (not os.path.exists(img_rgb) or args.rebuild):
#print fname,granule,granule['band']
if granule['band'] in ('B02','B03','B04'):
if not img_dir in jpegs:
jpegs[img_dir] = {
'tiles' : {}
}
if not granule['tilename'] in jpegs[img_dir]['tiles']:
jpegs[img_dir]['tiles'][granule['tilename']] = {
'epsg' : util.get_epsg_code(granule['zone']),
}
