This is to account for how I have handled multidim + non multidim input, see the following example unrotate scenario:

    # 1D input
    src_proj = ccrs.RotatedGeodetic(pole_longitude=178, pole_latitude=38)
    target_proj = ccrs.Geodetic()
    res = target_proj.transform_points(x=np.array([1, 2, 3, 4]),
                                                      y=np.array([1, 2, 3, 4]),
                                                      src_crs=src_proj)
    >>> res.shape
    (3, 1, 4)

    # Multidimensional input
    src_proj = ccrs.RotatedGeodetic(pole_longitude=178, pole_latitude=38)
    target_proj = ccrs.Geodetic()
    res = target_proj.transform_points(x=np.array([[1, 2, 3, 4],[5, 6, 7, 8]]),
                                                      y=np.array([[1, 2, 3, 4],[5, 6, 7, 8]]),
                                                      src_crs=src_proj)
    >>> res.shape
    (3, 2, 4)

Both scenarios return a shape that corresponds to if it was multidimensional (see [related cartopy issue](https://github.com/SciTools/cartopy/pull/192))