BUG: plot bounds are not computed correctly in cylindrical coordinates

[neutrinoceros]

Bug report

Bug summary

Code for reproduction

import yt
from yt.testing import fake_amr_ds, add_noise_fields

yt.load_uniform_grid(
    {},
    shape,
    bbox=np.array([[1, 5], [0.1, np.pi - 0.1], [0, 1]]),
    geometry="polar",
)
add_noise_fields(ds)
p = yt.SlicePlot(ds, "z", "noise1")
p.save("/tmp/")

Actual outcome

Expected outcome
We shouldn't waste half of the image space.
This is very similar to #3529

[neutrinoceros]

There's a similar issue for r-normal plots

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

shape = (32, 64, 128)
ds = yt.load_uniform_grid(
    {},
    shape,
    bbox=np.array([[1, 5], [0.1, np.pi / 2 - 0.1], [0, 1]]),
    geometry="polar",
)
add_noise_fields(ds)
p = yt.SlicePlot(ds, "r", "noise1")
p.save("/tmp/")