BUG: to_pw ignores geometry specific rules ?

[neutrinoceros]

Bug report

Bug summary
In spherical geometry (and probably other non-cartesian geometries too), creating a plot via ds.slice(...).to_pw produces inconsistent results as compared by the supposedly equivalent plot done with yt.SlicePlot

Code for reproduction

import numpy as np
import yt
from yt.testing import fake_amr_ds

ds = fake_amr_ds(geometry="spherical")
field = ("stream", "Density")

yt.SlicePlot(ds, "phi", field).save("/tmp/method1.png")

slc = ds.slice("phi", np.pi)
plot = slc.to_pw(field)
plot.save("/tmp/method2.png")

Actual outcome

Note how the horizontal axis labels are incorrect in the second case
This was originally discovered and reported by @vdbma

[matthewturk]

Is this because it calls _setup_plots?

[neutrinoceros]

No, I think it's because some of the input sanitizing is only done in plot classes. More specifically, I'm thinking of this

yt/yt/visualization/plot_window.py

Lines 1950 to 1958 in 6061a1a

	if ds.geometry in (
	"spherical",
	"cylindrical",
	"geographic",
	"internal_geographic",
	"polar",
	):
	mylog.info("Setting origin='native' for %s geometry.", ds.geometry)
	origin = "native"