Inaccurate distinction between ocean and land when plotting the Arctic

[xiuyongwu]

Now, I'm ready to plot the map for the Arctic, with land and ocean in different colors.
However, when I add_feature to the land and ocean, the color of ocean will overlap the land (or the continents near the Arctic). That is, the cartopy may be treat the ocean inaccurately in the Arctic.

The code for this map is shown as follows:

import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as cfeature
fig = plt.figure(figsize=[5,5])
ax1 = fig.add_subplot(1, 1, 1, projection=ccrs.NorthPolarStereo(central_longitude=0))
ax1.set_extent([-180, 180,60, 90], ccrs.PlateCarree())
ax1.gridlines()
ax1.add_feature(cfeature.LAND)
ax1.add_feature(cfeature.OCEAN)
plt.show()

However, if I exchange the order for the two sentences, write the "ax1.add_feature(cfeature.OCEAN)" first and then write "ax1.add_feature(cfeature.LAND)", it seems works.

Code to reproduce

Traceback

Full environment definition

Operating system

Cartopy version

conda list

pip list

[greglucas]

There is an interpolation that takes place when adding those features to a different projection than the one they are defined in. Some other options:

• Set the axes background color to the ocean color and overplot the land to make sure there are no gaps at all.
• The default OCEAN and LAND features are also very coarse resolution (110m), so you can try to use cfeature.LAND.with_scale('10m') to get higher resolution data to plot and look nicer too.