A 2D Heisenberg ferromagnet with exchange J and random magnetic anisotropy of strength has been studied. Analytical theory for the dependence of the average size of a pinned skyrmion on the magnetic field H, and for stability of such skyrmions on a lattice, has been developed. It has been complemented by numerical studies of 2D lattices containing up to 40 million spins. At low fields the average size of the skyrmion, λ, is determined by the average size of Imry–Ma domains. On increasing the field the skyrmions first shrink, with , and then collapse at fields distributed around . The concentration of the skyrmions goes down with the field as .