In the near future, a much larger scale of low earth orbit (LEO) satellites will be deployed, and satellite selection is an effective way to overcome challenges to processing capability and channel limitation of receivers caused by superabundant satellites in view. However, existing algorithms do not take into account the fact that a large-scale LEO constellation satellite selection subset can be closer to the theoretical optimal configuration than the traditional brute force search and recursive algorithm computational overhead is too large for end devices. In this study, the solution of the minimum optimization problem of the GDOP (Geometric Dilution of Precision) under the condition of limited elevation angle is first analyzed to determine the optimal ratio configuration scheme of top and bottom satellites under the optimal geometric satellite positioning problem. Subsequently, a fast satellite selection algorithm is proposed, which selects a subset of satellites that meet the different numbers of observation satellites and is close to the optimal geometric configuration by geometric method. The comparative simulation shows that, in the 1586 LEO constellation system, the proposed algorithm has an average GDOP of 1.52 when the satellite subset size is 10 and 0.95 when the subset size is 30. Compared with the recursive, the Rotating Partition algorithm and the Quasi-optimal algorithm, the proposed algorithm achieves better accuracy when its size is not close to the total number of observable satellites.