Large-scale multiobjective optimization problems have attracted increasing attention in both engineering applications and scientific research. Academically, large-scale multiobjective problems involve hundreds or thousands of decision variables. Due to the large decision space, the performance of traditional multiobjective evolutionary algorithms decreases dramatically when dealing with large-scale multiobjective problems, especially many-objective problems. With this in mind, a space sampling based large-scale many-objective evolutionary algorithm (LSMaOEA) is proposed in this article. Specifically, a space sampling method is developed that alternately performs upper/lower-linkage sampling and individual-linkage sampling to sample a set of individuals in the decision space. An environmental selection strategy based on nondominated sorting and reference vector association is proposed. Thus, the proposed LSMaOEA can alleviate excessively dense sampling at boundaries and improve the diversity of existing space sampling based algorithms for large-scale many-objective problems. In the experiments, the proposed algorithm is assessed by comparing it with eight state-of-the-art multi/many-objective evolutionary algorithms. The evaluation is conducted using two popular indicators across nine challenging multiobjective optimization benchmark problems with up to 2000 decision variables. The extensive experimental results consistently reveal that the proposed algorithm outperforms all the compared algorithms.