The global entropy solution to the radial motion of isothermal self-gravitating gases is studied. A transformation of variables is applied. The global approximate solutions are constructed by a generalized shock-capturing scheme. The source terms include an integral which is non-local and is carefully estimated. The uniform estimate of the approximate solutions is obtained. The convergence of the approximate solutions is proved by the compensated compactness method. The global existence of weak entropy solution with large $L^{\infty}$ initial data is established. The solution allows a possible vacuum state.