code: def test_transform_points_xyz(self):
      # Test geodetic transforms when using z value
      rx = np.array([2574.32516e3])
      ry = np.array([837.562e3])
      rz = np.array([5761.325e3])
  
      src_proj = ccrs.Geocentric()
      target_proj = ccrs.Geodetic()
  
      res = target_proj.transform_points(
          x=rx, y=ry, z=rz, src_crs=src_proj
      )
  
      glat = res[..., 0]
      glon = res[..., 1]
      galt = res[..., 2]
  
      # Solution generated by pyproj
      solx = np.array([18.0224043189])
      soly = np.array([64.9796515089])
      solz = np.array([5048.03893734])
  
      assert_arr_almost_eq(glat, solx)
      assert_arr_almost_eq(glon, soly)
      assert_arr_almost_eq(galt, solz)
  
imports: import copy
  from io import BytesIO
  import os
  from pathlib import Path
  import pickle
  import warnings
  import numpy
  from numpy.testing import assert_almost_equal, assert_array_equal
  from numpy.testing import assert_array_almost_equal
  import pyproj
  import pytest
  import shapely.geometry
  import cartopy.crs
label_desc: Use image node_id 10 for calling the TestCRS obj's underlying member method code with example usage: obj.test_transform_points_xyz() without return types