def test_reduce_bands_comparison_ops(self):
input = self.create_spacetime_layer_singleband()
input = gps.Pyramid({0: input})
imagecollection = GeotrellisTimeSeriesImageCollection(
input, InMemoryServiceRegistry())
visitor = GeotrellisTileProcessGraphVisitor()
graph = {
"gt": {
"arguments": {
"x": {
"from_argument": "data"
},
"y": 6.0
},
"process_id": "gt",
"result": True
}
}
visitor.accept_process_graph(graph)
stitched = imagecollection.reduce_bands(
visitor).pyramid.levels[0].to_spatial_layer().stitch()
print(stitched)
self.assertEqual(1, stitched.cells[0][0][0])
def test_reduce_bands(self):
input = self.create_spacetime_layer()
input = gps.Pyramid({0: input})
imagecollection = GeotrellisTimeSeriesImageCollection(
input, InMemoryServiceRegistry())
visitor = GeotrellisTileProcessGraphVisitor()
graph = {
"sum": {
"arguments": {
"data": {
"from_argument": "dimension_data"
}
},
"process_id": "sum"
},
"subtract": {
"arguments": {
"data": {
"from_argument": "dimension_data"
}
},
"process_id": "subtract"
},
"divide": {
"arguments": {
"data": [{
"from_node": "sum"
}, {
"from_node": "subtract"
}]
},
"process_id": "divide",
"result": True
}
}
visitor.accept_process_graph(graph)
stitched = imagecollection.reduce_bands(
visitor).pyramid.levels[0].to_spatial_layer().stitch()
print(stitched)
self.assertEqual(3.0, stitched.cells[0][0][0])