def METHOD_NAME(xyz: np.ndarray) -> np.ndarray:
    r"""
    xyz to R-theta-z transform::
      y       R
      |     /
      |   /
      | / theta
      *------------x
    .. math:: x = R \cos(\theta)
    .. math:: y = R \sin(\theta)
    Returns
    -------
    xyz : (3,) float ndarray
        the point in the local coordinate system
    """
    xyz = np.atleast_2d(xyz)
    assert len(xyz.shape) == 2, xyz.shape
    x = xyz[:, 0]
    y = xyz[:, 1]
    theta = np.degrees(np.arctan2(y, x))
    R = np.sqrt(x * x + y * y)
    return np.array([R, theta, xyz[:, 2]], dtype=xyz.dtype).T