Cross sections are presented for various Dy, Tb and Gd radionuclides produced in the proton bombardment of 159 Tb as well as for the reactions 152 Gd(p,4n) 149 Tb and 155 Gd(p,4n) 152 Tb up to 66 MeV. The experimental excitation functions are compared with theoretical predictions by means of the geometrydependent hybrid (GDH) model as implemented in the code ALICE/ASH, as well as with values from the TENDL-2012 library and previous literature experimental data, where available. Physical yields have been derived for the production of some of the medically important radioterbiums, namely 149 Tb (radionuclide therapy), 152 Tb (PET) and 155 Tb (SPECT). The indirect production of high-purity 155 Tb via the decay of its precursor 155 Dy is reported. The possibility of a large-scale production facility based on a commercial 70 MeV cyclotron is also discussed.