code: def test_natural_earth():
      ax = plt.axes(projection=ccrs.PlateCarree())
      ax.add_feature(cfeature.LAND)
      ax.add_feature(cfeature.OCEAN)
      ax.coastlines()
      ax.add_feature(cfeature.BORDERS, linestyle=":")
      ax.add_feature(cfeature.LAKES, alpha=0.5)
      ax.add_feature(cfeature.RIVERS)
      ax.set_xlim((-20, 60))
      ax.set_ylim((-40, 40))
      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 1 for calling a global function with example usage: test_natural_earth() and returns: ax