code: def test_ellipsoid_guam_transform(self):
      # USGS Professional Paper 1395, pp 339--340
      globe = ccrs.Globe(
          ellipse=None,
          semimajor_axis=6378206.4,
          flattening=1 - np.sqrt(1 - 0.00676866),
      )
      lat_0 = 13 + (28 + 20.87887 / 60) / 60
      lon_0 = 144 + (44 + 55.50254 / 60) / 60
      aeqd = ccrs.AzimuthalEquidistant(
          central_latitude=lat_0,
          central_longitude=lon_0,
          false_easting=50000.0,
          false_northing=50000.0,
          globe=globe,
      )
      geodetic = aeqd.as_geodetic()
  
      other_args = {
          "a=6378206.4",
          "f=0.003390076308689371",
          "lon_0=144.7487507055556",
          "lat_0=13.47246635277778",
          "x_0=50000.0",
          "y_0=50000.0",
      }
      check_proj_params("aeqd", aeqd, other_args)
  
      assert_almost_equal(
          np.array(aeqd.x_limits),
          [-19987726.36931940, 20087726.36931940],
          decimal=6,
      )
      assert_almost_equal(
          np.array(aeqd.y_limits),
          [-19919796.94787477, 20019796.94787477],
          decimal=6,
      )
  
      pt_lat = 13 + (20 + 20.53846 / 60) / 60
      pt_lon = 144 + (38 + 7.19265 / 60) / 60
      result = aeqd.transform_point(pt_lon, pt_lat, geodetic)
  
      # The paper uses an approximation, so we cannot match it exactly,
      # hence, decimal=1, not 2.
      assert_array_almost_equal(result, [37712.48, 35242.00], decimal=1)
  
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 7 for calling the TestAzimuthalEquidistant obj's underlying member method code with example usage: obj.test_ellipsoid_guam_transform() without return types