The hybrid energy system of hydro-powers, pumped storages and renewable energies has become a new topic direction in modern power system developments. Consequently, it is essential to realize a rational and efficient allocation of different energy source capacities. Nevertheless, there is still a gap between the available studies and the requirement for further hybrid energy system development. This paper focuses on the optimal capacity configuration of a wind, photovoltaic, hydropower, and pumped storage power system. In this direction, a bi-level programming model for the optimal capacity configuration of wind, photovoltaic, hydropower, pumped storage power system is derived. To model the operating mode of a pumped storage power station, two 0-1 variables are introduced. To handle the nonlinear and nonconvex lower level programing problem caused by the two 0-1 variables, it is proposed that the 0-1 variables are treated as some uncertain parameters. Also, by treating the 0-1 variables as some uncertain parameters, a two-stage robust optimization problem to decompose the original bi-level programing one into a master problem and a subproblem is finally introduced. The Karush-Kuhn-Tucker (KKT) conditions are then applied to simplify and linearize the min-max problem and nonlinear terms in the master problem. This results in both the master problem and the subproblem being formulated as mixed integer linear programming (MILP) problems. By utilizing the powerful Column-and-Constraint Generation (C&CG) algorithm, the two-stage robust optimization model is decomposed into an iterative procedure of solving the master problem and the subproblem sequentially. This approach eliminates the need for intricate optimization algorithms as commonly used in existing bi-level planning problems in hybrid energy systems. Finally, the effectiveness and advantages of the proposed model is verified by the numerical results on a case study.