Use single pass over datasets when computing bounds

Fixes issue #659
opendatacube · Feb 14, 2019 · 9bd8c07 · 9bd8c07
1 parent 8e63ab5
commit 9bd8c07
Showing 1 changed file with 2 additions and 5 deletions.
diff --git a/datacube/api/core.py b/datacube/api/core.py
@@ -659,11 +659,8 @@ def _fuse_measurement(dest, datasets, geobox, measurement,
 
 
 def get_bounds(datasets, crs):
-    left = min([d.extent.to_crs(crs).boundingbox.left for d in datasets])
-    right = max([d.extent.to_crs(crs).boundingbox.right for d in datasets])
-    top = max([d.extent.to_crs(crs).boundingbox.top for d in datasets])
-    bottom = min([d.extent.to_crs(crs).boundingbox.bottom for d in datasets])
-    return geometry.box(left, bottom, right, top, crs=crs)
+    bbox = geometry.bbox_union(ds.extent.to_crs(crs).boundingbox for ds in datasets)
+    return geometry.box(*bbox, crs=crs)
 
 
 def dataset_type_to_row(dt):