[step:Form a Gale transform and normalize its vectors onto the sphere]
Let $e_1,\dots,e_n$ denote the standard basis of $\mathbb{R}^n$. Define the homogenization map $A:\mathbb{R}^n \to \mathbb{R}^{r+1}$ by specifying
\begin{align*}
A e_i := (q_i,1)
\end{align*}
for each $i \in \{1,\dots,n\}$ and extending linearly. Since $q_1,\dots,q_n$ affinely span $\mathbb{R}^r$, the [linear map](/page/Linear%20Map) $A$ has rank $r+1$. Therefore its kernel
\begin{align*}
K := \ker A \subset \mathbb{R}^n
\end{align*}
has dimension
\begin{align*}
\dim K = n-(r+1)=d.
\end{align*}
Choose a linear isomorphism $\Phi:\mathbb{R}^d \to K$. For each $i \in \{1,\dots,n\}$, define the linear functional $\ell_i:\mathbb{R}^d \to \mathbb{R}$ by
\begin{align*}
\ell_i(a) := (\Phi(a))_i,
\end{align*}
where $(\Phi(a))_i$ denotes the $i$-th coordinate of the vector $\Phi(a) \in \mathbb{R}^n$. By the Euclidean representation of linear functionals on $\mathbb{R}^d$, there is a unique vector $g_i \in \mathbb{R}^d$ such that
\begin{align*}
\ell_i(a)=g_i \cdot a
\end{align*}
for every $a \in \mathbb{R}^d$. The ordered family $(g_1,\dots,g_n)$ is a Gale transform of the affine configuration $(q_1,\dots,q_n)$.
We claim that no Gale vector $g_i$ is zero. Fix $i \in \{1,\dots,n\}$. If $g_i=0$, then the coordinate functional $e_i^*:\mathbb{R}^n\to\mathbb{R}$, $x\mapsto x_i$, vanishes on $K=\ker A$. Thus every linear relation among the columns $(q_1,1),\dots,(q_n,1)$ has coefficient $0$ on the $i$-th column. Equivalently, the column $(q_i,1)$ is not in the linear span of the remaining columns, so deleting the $i$-th column lowers the rank of $A$.
On the other hand, since $n-1\geq 2k-1=r+1$, we may choose $r+1$ indices $m\neq i$. By affine general position of the moment-curve vertices, the corresponding points $q_m$ are affinely independent. Therefore their homogenizations $(q_m,1)$ are linearly independent in $\mathbb{R}^{r+1}$, so they span $\mathbb{R}^{r+1}$. Hence the columns with $m\neq i$ already have rank $r+1$, contradicting the conclusion that deleting the $i$-th column lowers the rank. Thus $g_i\neq 0$ for every $i$.
Define points $p_i \in S^{d-1}$ by
\begin{align*}
p_i := \frac{g_i}{|g_i|}
\end{align*}
for $i \in \{1,\dots,n\}$.
[/step]