A dataset of sea surface temperature (SST) estimates is generated from the temperature observations of surface drifting buoys of NOAA's Global Drifter Program. Estimates of SST at regular hourly time steps along drifter trajectories are obtained by fitting to observations a mathematical model representing simultaneously SST diurnal variability with three harmonics of the daily frequency, and SST low-frequency variability with a first degree polynomial. Subsequent estimates of non-diurnal SST, diurnal SST anomalies, and total SST as their sum, are provided with their respective standard uncertainties. This Lagrangian SST dataset has been developed to match the existing hourly dataset of position and velocity from the Global Drifter Program.