Create a CEA map from a full disk HMI magnetogram

I am looking to make a CEA map of a full-disk magnetogram from HMI but I cannot figure out how to do it with sunpy. I think there was an old tutorial in the example gallery but I cannot find that either. I already have the magnetogram downloaded and saved as a .fits file, and the header is completed with jsoc keywords but I am not sure what to do after. Thanks in advance for the help!

You can adapt this example: Creating Carrington Maps — sunpy 8.0.0 documentation

To get the CEA projection, pass projection_code="CEA" when calling make_heliographic_header().