Two caveats are noteworthy.
In both steps 2 and 3, datum shifts are ignored because (i) the relevant empirical parameters are not always available for converting all CRSs to WGS84 coordinates and (ii) datum information is usually missing from CF NetCDF files.
This can lead to geolocation errors of up to a few hundred metres. Secondly, this scheme does not implement any kind of interpolation: the resulting map images reproduce the shape of the source grid cells.