Abstract The surface temperature distributions of central compact objects (CCOs) are powerful probes of their crustal magnetic field strengths and geometries. Here we model the surface temperature distribution of RX J0822−4300, the CCO in the Puppis A supernova remnant, using 471 ks of XMM-Newton data. We compute the energy-dependent pulse profiles in 16 energy bands, fully including the general relativistic effects of gravitational redshift and light bending, to accurately model the two heated surface regions of different temperatures and areas, in addition to constraining the viewing geometry. This results in precise measurements of the two temperatures: kT warm = ( 1 + z ) × 0.222 − 0.019 + 0.018 keV and kT _hot = (1 + z) × 0.411 ± 0.011 keV. The two heated surface regions are likely located very close to the rotational poles, with the most probable position of the hotter component ≈ 6° from the rotational pole. For the first time, we are able to measure a deviation from a pure antipodal hot-spot geometry, with a longitudinal offset δ γ = 11 .° 7 − 2 .° 5 + 2 .° 6 . The discovery of this asymmetry, along with the factor of ≈2 temperature difference between the two emitting regions, may indicate that RX J0822−4300 was born with a strong, tangled crustal magnetic field.