def test_testutils_geobox():
from datacube.testutils.io import dc_crs_from_rio, rio_geobox
from rasterio.crs import CRS
from affine import Affine
assert rio_geobox({}) is None
transform = Affine(10, 0, 4676,
0, -10, 171878)
shape = (100, 640)
h, w = shape
crs = CRS.from_epsg(3578)
meta = dict(width=w, height=h, transform=transform, crs=crs)
gbox = rio_geobox(meta)
assert gbox.shape == shape
assert gbox.crs.epsg == 3578
assert gbox.transform == transform
wkt = '''PROJCS["unnamed",
GEOGCS["NAD83",
DATUM["North_American_Datum_1983",
SPHEROID["GRS 1980",6378137,298.257222101, AUTHORITY["EPSG","7019"]],
TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6269"]],
PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],
],
PROJECTION["Albers_Conic_Equal_Area"],
PARAMETER["standard_parallel_1",61.66666666666666],
PARAMETER["standard_parallel_2",68],
PARAMETER["latitude_of_center",59],
PARAMETER["longitude_of_center",-132.5],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",500000],
UNIT["Meter",1]]
'''
crs_ = dc_crs_from_rio(CRS.from_wkt(wkt))
assert crs_.epsg is None
def test_testutils_geobox():
from datacube.testutils.io import dc_crs_from_rio, rio_geobox
from rasterio.crs import CRS
from affine import Affine
assert rio_geobox({}) is None
A = Affine(10, 0, 4676, 0, -10, 171878)
shape = (100, 640)
h, w = shape
crs = CRS.from_epsg(3578)
meta = dict(width=w, height=h, transform=A, crs=crs)
gbox = rio_geobox(meta)
assert gbox.shape == shape
assert gbox.crs.epsg == 3578
assert gbox.transform == A
crs_ = dc_crs_from_rio(CRS.from_wkt(crs.wkt))
assert crs_.epsg is None