def METHOD_NAME(self):
    imagePath = datastore.fetch('landcover_sample_1000.tif')
    ts = self.basemodule.open(imagePath, projection='EPSG:3857', style='{"bands":[]}')
    # This gets the whole world
    region, _ = ts.getRegion(output=dict(maxWidth=1024, maxHeight=1024),
                             encoding='TILED')
    result = self.basemodule.open(str(region))
    tileMetadata = result.getMetadata()
    assert tileMetadata['bounds']['xmax'] == pytest.approx(20037508, 1)
    assert tileMetadata['bounds']['xmin'] == pytest.approx(-20037508, 1)
    assert tileMetadata['bounds']['ymax'] == pytest.approx(20037508, 1)
    assert tileMetadata['bounds']['ymin'] == pytest.approx(-20037508, 1)
    assert '+proj=merc' in tileMetadata['bounds']['srs']
    region.unlink()
    # Ask for a smaller part
    region, _ = ts.getRegion(
        output=dict(maxWidth=1024, maxHeight=1024),
        region=dict(left=-8622811, right=-8192317, bottom=5294998,
                    top=5477835, units='projection'),
        encoding='TILED')
    result = self.basemodule.open(str(region))
    tileMetadata = result.getMetadata()
    assert tileMetadata['bounds']['xmax'] == pytest.approx(-8192215, 1)
    assert tileMetadata['bounds']['xmin'] == pytest.approx(-8622708, 1)
    assert tileMetadata['bounds']['ymax'] == pytest.approx(5477783, 1)
    assert tileMetadata['bounds']['ymin'] == pytest.approx(5294946, 1)
    assert '+proj=merc' in tileMetadata['bounds']['srs']
    region.unlink()