code: def test_natural_earth_custom():
      ax = plt.axes(projection=ccrs.PlateCarree())
      feature = cfeature.NaturalEarthFeature(
          "physical",
          "coastline",
          "50m",
          edgecolor="black",
          facecolor="none",
      )
      ax.add_feature(feature)
      ax.set_xlim((-26, -12))
      ax.set_ylim((58, 72))
      return ax.figure
  
file_path: cartopy/lib/cartopy/tests/mpl/test_features.py
imports: from xml.etree.ElementTree import ParseError
  import matplotlib.pyplot
  import pytest
  import cartopy.crs
  import cartopy.feature
  from cartopy.tests.conftest import _HAS_PYKDTREE_OR_SCIPY
label_desc: Use image node_id 2 for calling a global function with example usage: test_natural_earth_custom() and returns: ax