def METHOD_NAME(self):
    wcs = WCS(self.msx_header)
    with galactocentric_frame_defaults.set("latest"):
        coord = SkyCoord(0 * u.kpc, 0 * u.kpc, 0 * u.kpc, frame="galactocentric")
    fig = plt.figure()
    ax = fig.add_subplot(1, 1, 1, projection=wcs)
    (point,) = ax.plot_coord(coord, "ro")
    np.testing.assert_allclose(point.get_xydata()[0], [0, 0], atol=1e-4)