The assimilation of satellite SST data into ocean models is complex. The temperature of the thin uppermost layer that is measured by satellites may differ from the much thicker upper layer used in numerical models, leading to biased results. This paper shows how canonical correlation analysis can be used to generate observation operators from existing datasets of model states and corresponding observation values. This type of operator can correct for near-surface effects when assimilating SST.