

def test_warp_reproject_like(runner, tmpdir):
    likename = str(tmpdir.join('like.tif'))
    kwargs = {
        'crs': {
            'init': 'epsg:4326',
        },
        'transform': ((- 106.523), 0.001, 0, 39.6395, 0, (- 0.001)),
        'count': 1,
        'dtype': rasterio.uint8,
        'driver': 'GTiff',
        'width': 10,
        'height': 10,
        'nodata': 0,
    }
    with rasterio.drivers():
        with rasterio.open(likename, 'w', **kwargs) as dst:
            data = numpy.zeros((10, 10), dtype=rasterio.uint8)
            dst.write(data, indexes=1)
    srcname = 'tests/data/shade.tif'
    outputname = str(tmpdir.join('test.tif'))
    result = runner.invoke(warp.warp, [srcname, outputname, '--like', likename])
    assert (result.exit_code == 0)
    assert os.path.exists(outputname)
    with rasterio.open(outputname) as output:
        assert (output.crs == {
            'init': 'epsg:4326',
        })
        assert numpy.allclose([0.001, 0.001], [output.affine.a, (- output.affine.e)])
        assert (output.width == 10)
        assert (output.height == 10)

dataset/ETHPy150Open mapbox/rasterio/tests/test_rio_warp.py test_warp_reproject_like/original