code: def test_gshhs():
      ax = plt.axes(projection=ccrs.Mollweide())
      ax.set_extent([138, 142, 32, 42], ccrs.Geodetic())
  
      ax.stock_img()
      # Draw coastlines.
      ax.add_feature(cfeature.GSHHSFeature("coarse", edgecolor="red"))
      # Draw higher resolution lakes (and test overriding of kwargs)
      ax.add_feature(
          cfeature.GSHHSFeature("low", levels=[2], facecolor="green"),
          facecolor="blue",
      )
      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 3 for calling a global function with example usage: test_gshhs() and returns: ax