yes you're correct.
note that row echelon form doesn't necessarily result in 'leading 1s'. it's 'reduced/canonical row echelon form' that requires that form.
having reduced your matrix to the set of the linearly independent rows/columns via the row transformations, you can choose either the new reduced vectors with leading pivots (1s or otherwise), or the corresponding vectors from the original matrix*. they are effectively 'the same'. i'd go with the reduced vectors however, as they make any further manipulation or plotting easier
*see caveat raised by user84413
Your procedures are rather non-standard in terms of finding basis for columnspaces and rowspaces. In particular your second procedure for finding the rowspace is not only unorthodox, but appears to be incorrect. To fully answer all of your questions takes a bit of time, and I apologize in advance for the length. The fact is, you can find both basis in a single shot by reducing the original matrix, say $A$, to its RREF, say $R$.
First of all, the RREF is fundamentally defined in terms of the rows. That's why it is the Reduced Row Echelon Form. Elementary row operations by design do not change the rowspace; each new row remains a linear combination of the old rows. That means you can row reduce your matrix all you want, but your rowspace remains the same. It follows that $R$ and $A$ share a rowspace.
But the rows of $R$ are linearly independent: each row has an entry, its leading pivot entry, that no other row has. The rows of $R$ clearly also span the rowspace of $R$, just by definition, and therefore it is a basis for the rowspace of $R$. Now since $\mathrm{row}(R) = \mathrm{row}(A)$, it follows that the rows of $R$ also form a basis for the rowspace of $A$. There is no need to go back to the original matrix because the rowspace has not changed between elementary row operations.
Finding a basis for the columnspace is a bit more complicated because elementary row operations does not preserve the columnspace. The columns of $R$ do not necessarily span the same space as the columns of $A$. Of course, the columnspace of $A$ is just the rowspace of $A^\mathrm{T}$ so we can always apply the previous procedure to $A^\mathrm{T}$. I believe this is what your second procedure tries to do, but it doesn't get it quite right. More on that later. In any case, reducing two matrices is hard work and luckily we can actually kill two birds with one stone.
If we encode the sequence of elementary row operations which take $A$ to $R$ as elementary matrices, then we can write
$$E_k\cdots E_1 A = R$$
where each $E_i$ is an elementary matrix. Let us just collectively write
$$EA = R$$
where $E$ is the product of all our elementary matrices. The important thing here is that $E$ is an invertible matrix since each $E_i$ is invertible. This fact actually characterizes row equivalent matrices: two (same sized) matrices $A$ and $B$ are row equivalent if and only if there exists an invertible matrix $E$ such that $A=EB$.
Let us write $\mathbf{a}_i$ as the $i$th column of $A$ and $\mathbf{r}_i$ as the $i$th column of $R$. By block matrix multiplication, it follows that
$$E\mathbf{a}_i = \mathbf{r}_i$$
My claim now is that linear relations are preserved between the columns of $A$ and $R$, i.e. for any set of scalar coefficients $\{c_i\}$, we have
$$\mathbf{0} = \sum_{i=1}^n c_i\mathbf{a}_i \iff \mathbf{0}=\sum_{i=1}^n c_i \mathbf{r}_i$$
The proof is simple:
$$\mathbf{0} = \sum_{i=1}^n c_i\mathbf{r}_i \iff \mathbf{0}=\sum_{i=1}^n c_iE\mathbf{a_i} \iff \mathbf{0}= E\left(\sum_{i=1}^n c_i\mathbf{a}_i\right) \iff \mathbf{0} = \sum_{i=1}^n c_i\mathbf{a}_i$$
The last line follows because $E$ is invertible and therefore its nullspace is just $\{\mathbf{0}\}$.
The above relations are very powerful and we have some immediate corollaries:
Let $\mathcal{A}=\{\mathbf{a}_{i_k}\}$ be a subset of $\{\mathbf{a}_i\}$ and let $\mathcal{R} = \{\mathbf{r}_{i_k}\}$ be the corresponding subset of $\{\mathbf{r}_i\}$. Then
$\mathcal{A}$ is linearly independent if and only if $\mathcal{R}$ is linearly independent.
$\mathcal{A}$ spans $\mathrm{col}(A)$ if and only if $\mathcal{R}$ spans $\mathrm{col}(R)$.
The two points above imply that $\mathcal{A}$ forms a basis for $\mathrm{col}(A)$ if and only if $\mathcal{R}$ forms a basis for $\mathrm{col}(R)$.
I will not prove the above propositions, the answer is getting a bit long and I think it would make a good exercise. I encourage you to attempt a proof yourself. The implications of the corollaries are straightforward: if you are given a basis for the columnspace of $R$, taking the corresponding vectors will give a basis for $A$. Taking the pivot columns of $R$ serves as an immediate basis for $\mathrm{col}(R)$ and the corresponding columns of $A$ serves as a basis for $\mathrm{col}(A)$. This is the underlying theory behind your first procedure.
To reiterate, take care to note that the pivot columns of $R$ do not form a basis for $\mathrm{col}(A)$ themselves, only the corresponding columns of $A$ do. Elementary row operations are designed to preserve the rowspace and not the columnspace. Oftentimes, $\mathrm{col}(A) \neq \mathrm{col}(R)$ and it makes no sense to take columns of $R$ to act as a basis for the columnspace of $A$. This is the fundamental reason why one of the procedures need to go back to the original matrix while the other one does not.
This now enables you to efficiently find the basis for the rowspace and columnspace all in one go. And we've explained how your first procedure works. With the benefit of hindsight, you should be able to see why your second procedure is a bit weird. You're row reducing the transpose, which means that you are essentially preforming elementary column operations on your original matrix. These operations do not preserve the rowspace and so the procedure as you've described it cannot be correct. At best, you are just preforming the procedure for finding the columnspace basis on $A^\mathrm{T}$, but that still requires you to go back to your original matrix. Even if your second procedure were correct, it turns out to be largely redundant as we've shown that simply reducing $A$ to its RREF $R$ is sufficient to find a basis for the columnspace and the rowspace (and the nullspace if you want).
Best Answer
[I’m sure that this must have been answered before, but a quick scan of the related questions at right didn’t turn up anything that addressed this particular point.]
Elementary row operations correspond to left-multiplication by particular invertible square matrices called elementary matrices. Thus, if $R$ is the rref of $A$, then $R=E_k\cdots E_2E_1A=SA$ for some invertible square matrix $S$. The column space of $R$ is obviously spanned by its pivot columns, which consist of a $1$ in some slot and zeros everywhere else. To find a basis for the column space of $A$, we need to let $S^{-1}$ act on the column space of $R$, but since each column of $R$ is the image under $S$ of the corresponding column of $A$, the images of $R$’s pivot columns under $S^{-1}$ are the corresponding columns of $A$. $S$ has full rank, so these columns are linearly independent and the rank-nullity theorem tells us that they form a basis for $C(A)$.
I’ll add here that finding a basis for the column space of $A$ by computing the rref of $A^T$ (eqv. column-reducing $A$) usually results in a “nicer” basis than the one you get by picking out a set of linearly-independent columns of $A$.