Isolated neutron stars are highly magnetized, fast-rotating objects that form as an end point of stellar evolution. They are directly observable in X-ray emission, because of their high surface temperatures. Features in their X-ray spectra could in principle reveal the presence of atmospheres, or be used to estimate the strength of their magnetic fields through the cyclotron process, as is done for X-ray binaries. Almost all isolated neutron star spectra observed so far appear as featureless thermal continua. The only exception is 1E1207.4-5209 (refs 7-9), where two deep absorption features have been detected, but with insufficient definition to permit unambiguous interpretation. Here we report a long X-ray observation of the same object in which the star's spectrum shows three distinct features, regularly spaced at 0.7, 1.4 and 2.1 keV, plus a fourth feature of lower significance, at 2.8 keV. These features vary in phase with the star's rotation. The logical interpretation is that they are features from resonant cyclotron absorption, which allows us to calculate a magnetic field strength of 8 × 1010 G, assuming the absorption arises from electrons.