I am trying to plot data for the region of American Samoa but I am having issues with the coastlines. I generally use the higher resolution data from the Global Self-consistent, Hierarchical, High-resolution Geography Database but the coastline appears shifted from the data I am plotting and from other datasets such as NaturalEarth and Google maps.

I try to plot a map with cartopy, and label latitudes with some digits.

I am trying to draw polygons on a map at arbitrary locations, including places that span the pole and the dateline. Conceptually, consider drawing instrument footprints for orbital measurements, where the corners of the footprints are known in lat/long (or cartesian). I have been able to get mid-latitude polygons to draw, but shapes spanning the pole come up blank.

consider the following script that gives grids of matplotlib axes with a cartopy projection