In this paper, a radial basis function method combined with the stochastic Galerkin method is considered to approximate elliptic optimal control problem with random coefficients. This radial basis function method is a meshfree approach for solving high dimensional random problem. Firstly, the optimality system of the model problem is derived and represented as a set of deterministic equations in high-dimensional parameter space by finite-dimensional noise assumption. Secondly, the approximation scheme is established by using finite element method for the physical space, and compactly supported radial basis functions for the parameter space. The radial basis functions lead to the sparsity of the stiff matrix with lower condition number. A residual type a posteriori error estimates with lower and upper bounds are derived for the state, co-state and control variables. An adaptive algorithm is developed to deal with the physical and parameter space, respectively. Numerical examples are presented to illustrate the theoretical results.