GDAL3 support in ls_public_bucket.py script

[Kirill888]

GDAL3 swaps order of Lon/Lat and generated extents are not correct when running this script in the environment with GDAL3 installed.

cube-in-a-box/scripts/ls_public_bucket.py

Lines 80 to 85 in eb9ffa6

	def get_coords(geo_ref_points, spatial_ref):
	t = osr.CoordinateTransformation(spatial_ref, spatial_ref.CloneGeogCS())
	def transform(p):
	lon, lat, z = t.TransformPoint(p['x'], p['y'])
	return {'lon': lon, 'lat': lat}

In GDAL3 output of TransformPoint is lat,lon,z instead. But that code is wrong in other ways, since we need bounding box of the geometry in EPSG:4326 and not a mere transform or the 4 corners to ESPG:4326. This also doesn't work when dataset's footprint spans lon=180 line.
opendatacube/datacube-core#886

Also shouldn't that script be in datacube-dataset-config repo?
opendatacube/datacube-dataset-config#10

[alexgleith]

Yes, the script should be in datacube-dataset-config and it should be a git submodule. There's a fair bit of consolidation that should be undertaken in the indexing process in this repo.