code: def test_sphere_transform(self):
      # USGS Professional Paper 1395, pg 337
      globe = ccrs.Globe(
          ellipse=None, semimajor_axis=3.0, semiminor_axis=3.0
      )
      aeqd = ccrs.AzimuthalEquidistant(
          central_latitude=40.0, central_longitude=-100.0, globe=globe
      )
      geodetic = aeqd.as_geodetic()
  
      other_args = {
          "a=3.0",
          "b=3.0",
          "lon_0=-100.0",
          "lat_0=40.0",
          "x_0=0.0",
          "y_0=0.0",
      }
      check_proj_params("aeqd", aeqd, other_args)
  
      assert_almost_equal(
          np.array(aeqd.x_limits), [-9.42477796, 9.42477796], decimal=6
      )
      assert_almost_equal(
          np.array(aeqd.y_limits), [-9.42477796, 9.42477796], decimal=6
      )
  
      result = aeqd.transform_point(100.0, -20.0, geodetic)
  
      assert_array_almost_equal(result, [-5.8311398, 5.5444634])
  
file_path: cartopy/lib/cartopy/tests/crs/test_azimuthal_equidistant.py
imports: import numpy
  from numpy.testing import assert_almost_equal, assert_array_almost_equal
  import cartopy.crs
  from .helpers import check_proj_params
label_desc: Use image node_id 5 for calling the TestAzimuthalEquidistant obj's underlying member method code with example usage: obj.test_sphere_transform() without return types