def statistics(xds, dv='IMAGE', mode='classic'):
    """
    Generate statistics on specified image data contents
    
    Resulting data is placed in the attributes section of the dataset

    Parameters
    ----------
    xds : xarray.core.dataset.Dataset
        input Image Dataset
    dv : str
        name of data_var in xds to compute statistics on. Default is 'IMAGE'
    mode : str
        algorithm mode to use ('classic', 'fit-half', 'hinges-fences', 'chauvenet', 'biweight'). Default is 'classic'
    
    Returns
    -------
    xarray.core.dataset.Dataset
        output Image
    """
    return xds
