High angular resolution observations at optical wavelengths provide valuable insights into stellar astrophysics1,2, and enable direct measurements of fundamental stellar parameters3,4 and the probing of stellar atmospheres, circumstellar disks5, the elongation of rapidly rotating stars6 and the pulsations of Cepheid variable stars7. The angular size of most stars is of the order of one milliarcsecond or less, and to spatially resolve stellar disks and features at this scale requires an optical interferometer using an array of telescopes with baselines on the order of hundreds of metres. We report on the implementation of a stellar intensity interferometry system developed for the four VERITAS imaging atmospheric Cherenkov telescopes. The sy...