apply_by_band function

Apply numeric function by band