Color interpretation test failure with Rasterio 1.0 and GDAL 2.2

[sgillies]

This test fails, and seems to fail because the source file has no photometric interpretation explicitly set, and thus none is transferred via profile to the output file.

tmpdir = local('/private/var/folders/zq/fv6jds3j4wb_0z0zh35xky_c0000gn/T/pytest-of-sean/pytest-671/test_color_cli_16bit_photointe0')

    def test_color_cli_16bit_photointerp(tmpdir):
        output = str(tmpdir.join('color16color.tif'))
        runner = CliRunner()
        result = runner.invoke(
            color,
            ['-d', 'uint16',
             '-j', '1',
             'tests/rgb16.tif',
             output,
             "gamma 3 1.85",
             "gamma 1,2 1.95"])
        assert result.exit_code == 0

        with rasterio.open('tests/rgb16.tif') as src:
            with rasterio.open(output) as out:
                for b in src.indexes:
>                   assert out.colorinterp == src.colorinterp
E                   assert (<ColorInterp...undefined: 0>) == (<ColorInterp....terp.blue: 5>)
E                     At index 0 diff: <ColorInterp.grey: 1> != <ColorInterp.red: 3>
E                     Use -v to get the full diff

Otherwise, everything checks out with Rasterio 1.0rc5 and rio-mucho 1.0dev1.