Stack band files with rasterio

stack_band_files.py

def stack_band_files(infiles, outfile, creation_options={}):
    """Stack single-band files into one
    
    Parameters
    ----------
    infiles : list of str
        paths to input files
    outfile : str
        path to output file
    creation_options : dict
        destination file creation options
    """
    with rasterio.open(infiles[0]) as src:
        profile = src.profile.copy()
    profile.update(count=len(infiles))
    profile.update(creation_options)
    with rasterio.open(outfile, 'w', **profile) as dst:
        for i, infile in enumerate(infiles):
            dst.write_band(i+1, rasterio.band(src, 1))