The present paper illustrates a methodology for the safety assessment of masonry domes. The dome is modelled as a thin shell made of a material satisfying Heyman’s hypotheses. Based on the static theorem of limit analysis, the method searches for statically admissible distributions of internal forces within the shell, suitably combining membrane forces and bending moments, by solving a convex optimisation problem. The solution is pursued numerically by means of an expressly developed collocation method that enables obtaining the analytical expressions for each internal force component. In its present formulation the method can be applied to domes of any shape, as well as to arbitrary load distributions. After validation against the benchmark case of the spherical dome under its self-weight, the paper illustrates application of the method to the dome of Pisa Cathedral under vertical loads as a first real case study.