We show how to find upper and lower bounds to the median of a gamma distribution, over the entire range of shape parameter $k > 0$, that are the tightest possible bounds of the form $2^{-1/k} (A + Bk)$, with closed-form parameters $A$ and $B$. The lower bound of this form that is best at high $k$ stays between 48 and 50 percentile, while the uniquely best upper bound stays between 50 and 55 percentile. We show how to form even tighter bounds by interpolating between these bounds, yielding closed-form expressions that more tightly bound the median. Good closed-form approximations between the bounds are also found, including one that is exact at $k = 1$ and stays between 49.97 and 50.03 percentile.