use astropy unit for `flux_star_removal_threshold`
use astropy unit for `flux_star_removal_threshold`

Originally my choice is not to use astropy unit, but now that we are using it, we should consistently use it throughout this function.
sagasurvey/saga,sagasurvey/saga