Let $X$ be a real normed space equipped with $\lVert\cdot\rVert$ and $M$ be a convex set in $X$, such that $B(0,r) \subset M \subset B(0,R)$ for $r>0$. Here $B(a, t)$ stands for closed ball with center $a$ and radius $t$.
Further, let $\pi\colon S \rightarrow \partial M$ denote the central projection of the unit sphere onto boundary of $M$. We can show that the best Lipschitz constant $L(\pi)$ of the projection is at most $\frac{R(R+1)}{r}$$\frac{R(R+r)}{r}$. If $\mu$ is Minkowski functional of $M$, then $\pi(x) = \frac{x}{\mu(x)}$. Hence, for all $x,y \in S$ \begin{align*} \frac{\lVert \pi(x)-\pi(y) \rVert}{\lVert x-y \rVert} &= \frac{\lVert \mu(y)\,x-\mu(x)\,y \rVert}{\lVert x-y \rVert \,\mu(x)\mu(y)} \leqslant \frac{\lvert \mu(y)-\mu(x) \rvert \,\lVert x \rVert + \mu(x)\,\lVert x-y \rVert}{\lVert x-y \rVert \,\mu(x)\mu(y)} = \\ &= \frac{1}{\mu(x)\mu(y)} \frac{\lvert \mu(y)-\mu(x) \rvert}{\lVert x-y \rVert} + \frac{1}{\mu(y)} \leqslant \bigg( \frac{R}{r} \bigg)^2 \frac{1}{r} + \frac{R}{r} = \frac{R(R+1)}{r}, \end{align*}\begin{align*} \frac{\lVert \pi(x)-\pi(y) \rVert}{\lVert x-y \rVert} &= \frac{\lVert \mu(y)\,x-\mu(x)\,y \rVert}{\lVert x-y \rVert \,\mu(x)\mu(y)} \leqslant \frac{\lvert \mu(y)-\mu(x) \rvert \,\lVert x \rVert + \mu(x)\,\lVert x-y \rVert}{\lVert x-y \rVert \,\mu(x)\mu(y)} = \\ &= \frac{1}{\mu(x)\mu(y)} \frac{\lvert \mu(y)-\mu(x) \rvert}{\lVert x-y \rVert} + \frac{1}{\mu(y)} \leqslant \bigg( \frac{R}{r} \bigg)^2 \frac{1}{r} + R = \frac{R(R+r)}{r}, \end{align*} as for every $z \in X$ $$ \frac{\lVert z \rVert}{R} \leqslant \mu(z) \leqslant \frac{\lVert z \rVert}{r},$$ in particular, $\mu$ is $1/r$-Lipschitz.
So, we're getting closer to the point. Is this an exact bound of $L(\pi)$? No, at least in the case of X being inner product space. The exact bound then is $\frac{R^2}{r}$, it's proved in the paper "A note on starshaped sets" by Siniša Vrećica.
In fact, when dealing with $\frac{\lVert \pi(x)-\pi(y) \rVert}{\lVert x-y \rVert}$ we work inside $Y=\mathrm{span} \{x,y\}$ as $\pi(x),\pi(y) \in Y$. That means we can consider only two-dimensional normed spaces $X$. The proof of euclidean (planar) case of the estimate pretty uses trigonometry. (Siniša's proof looked too technical to me, so I proved it myself using just sines' theorem and elementary facts, like monotonicity of sine.)
It's also sufficient to consider convex bodies of the form $M = \mathrm{conv} (B(0,r) \cup \{z\})$, as that sets are the "worst" for $L(\pi)$. Also, WLOG, we can take $r=1$.
For example, due to my humble computations, for $X = \ell_\infty^2$ and $$M_1 = \mathrm{conv} (B(0,1) \cup \{(0,R)\}), \quad M_2 = \mathrm{conv} (B(0,1) \cup \{(R,R)\}), \quad R\geqslant 2$$ we have $L(\pi_1) = R(R-1)$ and $L(\pi_2) = R(R+1)/2$.
Now, I have no idea how to deal with non-euclidean case in general. It seems to me that the Lipschitz constant must be less than for euclidean case. Anyway, how can one obtain some estimations better than $\frac{R(R+1)}{r}$$\frac{R(R+r)}{r}$ (say, norm given by radial function of unit ball)? Or there is example when we can't do better?