In master's thesis we discuss methods of optimization of geodetic networks. Optimization of geodetic networks allow us in advance determination of optimal geodetic datum, optimal positions of points in geodetic network, optimal weights of observations and optimal improvement of weak parts of network. The focus in this thesis is on second order design of optimization where optimal precision of geodetic observations is optimised according to a desired accurracy of points coordinates in geodetic network. For the automatization of the second order design, we made a computer program for second order design of geodetic network. Program allows importing data of geodetic network and criteria, which the program takes into account in the process of optimization. Program iteratively removes observations, which do not contribute significantly to the accuracy or reliability of geodetic network, until no observations can be removed, without breaking the criteria, which has been set at the beginning of the optimization. Using the assembled program, we solved few geodetic network optimization examples with different criteria to show the applicability of the computer program and methods of optimization for optimal geodetic network design.