def add_geotiff_ax(ax, fp, band_ind, ix_t, crs):
    """
    Add a projected GeoTIFF raster to the Axes instance
    """
    import rasterio
    src = rasterio.open(fp)
    extent = (src.bounds.left, src.bounds.right, src.bounds.bottom, src.
        bounds.top)
    for ix_b in band_ind:
        ax.imshow(src.read(ix_b)[ix_t], transform=crs, extent=extent,
            origin='upper')
    return ax
