def test(self):
ndvi_i = indices.ndvi('B8', 'B4', addBand=False)(image)
evi_i = indices.evi('B8', 'B4', 'B2', addBand=False)(image)
nbr_i = indices.nbr('B8', 'B12', addBand=False)(image)
nbr2_i = indices.nbr2('B11', 'B12', addBand=False)(image)
values = get_value(image, p, 10, side='client')
NIR = float(values['B8'])
SWIR1 = float(values['B11'])
SWIR2 = float(values['B12'])
RED = float(values['B4'])
BLUE = float(values['B2'])
GREEN = float(values['B3'])
ndvi = (NIR - RED) / (NIR + RED)
evi = (2.5) * ((NIR - RED) / (NIR + (6 * RED) - ((7.5) * BLUE) + 1))
nbr = (NIR - SWIR2) / (NIR + SWIR2)
nbr2 = (SWIR1 - SWIR2) / (SWIR1 + SWIR2)
ndvi_from_i = get_value(ndvi_i, p, 10, 'client')['NDVI']
evi_from_i = get_value(evi_i, p, 10, 'client')['EVI']
nbr_from_i = get_value(nbr_i, p, 10, 'client')['NBR']
nbr2_from_i = get_value(nbr2_i, p, 10, 'client')['NBR2']
self.assertEqual(ndvi, ndvi_from_i)
self.assertEqual(evi, evi_from_i)
self.assertEqual(nbr, nbr_from_i)
self.assertEqual(nbr2, nbr2_from_i)
def test_ndvi(self):
ndvi_col = col.map(indices.ndvi('B8', 'B4'))
ndvi_i = ee.Image(ndvi_col.first())
values = get_value(ndvi_i, p, 10, side='client')
NIR = float(values['B8'])
RED = float(values['B4'])
ndvi = (NIR - RED) / (NIR + RED)
ndvi_from_i = get_value(ndvi_i, p, 10, 'client')['NDVI']
self.assertEqual(ndvi, ndvi_from_i)
def test_expressions(self):
from geetools import expressions
generator = expressions.ExpGen()
exp_max = generator.max("b('B1')", "b('B2')")
exp_min = generator.min("b('B1')", "b('B2')")
img_max = self.l8SR.expression(exp_max).select([0], ["max"])
img_min = self.l8SR.expression(exp_min).select([0], ["min"])
vals_max = get_value(img_max, self.p_l8SR_no_cloud, 30, 'client')
vals_min = get_value(img_min, self.p_l8SR_no_cloud, 30, 'client')
self.assertEqual(vals_max["max"], 580)
self.assertEqual(vals_min["min"], 517)
def test_clouds(self):
masked = cloud_mask.landsat457TOA_BQA(['cloud'])(self.image)
vals = get_value(masked, self.p_cloud, 30, 'client')
show(self.iid, self.image, masked, bands457, 0, 0.5)
self.assertEqual(vals["B1"], None)
def test(self):
iid = 'LANDSAT/LT05/C01/T1_TOA'
collection = ee.ImageCollection(iid)
date = '2010-03-21'
p_cloud = ee.Geometry.Point([-65.2835, -24.3422])
p_clear = ee.Geometry.Point([-64.6106, -24.7469])
image = getimage(collection, date, area)
masked = cloud_mask.landsat457TOA_BQA()(image)
vals = get_value(masked, p_cloud, 30, 'client')
show(iid, image, masked, bands457, 0, 0.5)
self.assertEqual(vals["B1"], None)
def test(self):
iid = 'COPERNICUS/S2'
collection = ee.ImageCollection(iid)
date = '2017-03-07'
p_cloud = ee.Geometry.Point([-65.84304, -24.82382])
p_clear = ee.Geometry.Point([-65.88415, -24.82608])
image = getimage(collection, date, area)
masked = cloud_mask.hollstein_S2(['cloud','snow','shadow', 'cirrus'])(image)
vals = get_value(masked, p_cloud, 30, 'client')
show(iid, image, masked, bandsS2, 0, 5000)
self.assertEqual(vals["B1"], None)
def test(self):
iid = 'LANDSAT/LC08/C01/T1_SR'
collection = ee.ImageCollection(iid)
date = '2017-03-08'
p_cloud = ee.Geometry.Point([-65.7539, -25.0327])
p_clear = ee.Geometry.Point([-66.034, -24.9512])
image = getimage(collection, date, area)
# masked = cloud_mask.cfmask_bits(image)
masked = cloud_mask.landsatSR()(image)
vals = get_value(masked, p_cloud, 30, 'client')
show(iid, image, masked, bands8, 0, 5000)
self.assertEqual(vals["B1"], None)
def test(self):
iid = 'LANDSAT/LE07/C01/T1_SR'
collection = ee.ImageCollection(iid)
date = '2010-03-13'
p_cloud = ee.Geometry.Point([-65.5784, -24.773])
p_clear = ee.Geometry.Point([-65.599, -24.9313])
image = getimage(collection, date, area)
# masked = cloud_mask.cfmask_bits(image)
masked = cloud_mask.landsatSR()(image)
vals = get_value(masked, p_cloud, 30, 'client')
show(iid, image, masked, bands457, 0, 5000)
self.assertEqual(vals["B1"], None)
def test(self):
iid = 'LANDSAT/LT04/C01/T1_SR'
collection = ee.ImageCollection(iid)
date = '1989-07-25'
p_cloud = ee.Geometry.Point([-65.506, -24.9263])
p_clear = ee.Geometry.Point([-64.6133, -24.4522])
image = getimage(collection, date, area)
# masked = cloud_mask.cfmask_bits(image)
masked = cloud_mask.landsatSR()(image)
vals = get_value(masked, p_cloud, 30, 'client')
show(iid, image, masked, bands457, 0, 5000)
self.assertEqual(vals["B1"], None)
def test_snow(self):
masked = cloud_mask.landsat457TOA_BQA(['snow'])(self.image)
vals = get_value(masked, self.p_cloud, 30, 'client')
show(self.iid, self.image, masked, bands457, 0, 0.5)