import cartopy.crs as ccrs
import cartopy.feature as cfeature
import matplotlib.pyplot as plt

crs_proj = ccrs.LambertAzimuthalEqualArea(central_latitude=-90)
fig, ax = plt.subplots(1, 1, subplot_kw={"projection": crs_proj})
ax.add_feature(cfeature.OCEAN)
ax.coastlines()
plt.show()

import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as cfeature


def ocean_fig():
    fig = plt.figure()
    ax = fig.add_subplot(111, projection=ccrs.NorthPolarStereo())

    ax.add_feature(cfeature.OCEAN.with_scale('10m'))
    
    return fig


import timeit

setup = """
from __main__ import ocean_fig
fig = ocean_fig()
"""

print(timeit.timeit('fig.draw_without_rendering()', setup=setup, number=5))