def sample_statistics(image, geom, numPixels, scale):
    """ Sample within a region of a single band float image,
    returning a dictionary of statistics:
    accessed as ["values"]["<stat>"]
    Note this is an estimate"""
    bn = image.bandNames().getInfo()[0]
    fc = image.sample(region=geom, numPixels=numPixels, scale=scale, tileScale=16, dropNulls=True)
    return fc.aggregate_stats(bn)