def __init__(self, fileName, logLevel=30, cloudBits=None, debug=False, GCP_COUNT=10):
'''
Check if image is small or taken at night, add mask
'''
RasterYSize, RasterXSize, day_night_flag, mapperName = self.get_size_day_mapper(fileName)
isday = day_night_flag == 'Day'
isbig = (RasterXSize > 200) * (RasterYSize > 200)
if not isday:
print 'Night image'
time.sleep(0.5)
raise Exception('Night image!')
if not isbig:
print 'Too little image!'
time.sleep(0.5)
raise Exception('Too little image!')
Nansat.__init__(self, fileName, logLevel=logLevel, mapperName=mapperName, GCP_COUNT=GCP_COUNT)
# good bits for amazon plume studying
if cloudBits is None:
#cloudBits=[4, 5, 6, 10, 15, 20, 23, 30]
# good bits for NRT
cloudBits = [1, 6, 10, 16, 29]
#cloudBits=[1, 4, 5, 6, 9, 10, 13, 15, 20, 21, 23, 28, 29, 30]
# yellow sea:
#cloudBits = [1, 4, 5, 6, 9, 10, 15, 20, 29]
#cloudBits=[1, 4, 5, 6, 9, 10, 15, 20, 21, 29]
# for china
#cloudBits=[1, 5, 6, 9, 10, 13, 15, 20, 23, 28, 29, 30]
self.add_mask(cloudBits, debug=debug)
# replace cdom_index and sst with cdom_a and SST
self.rename_band('cdom_index', 'cdom_a')
self.rename_band('sst', 'SST')
