cartopy_example.py

import numpy as np
import matplotlib.pyplot as plt

import cartopy.crs as ccrs
import cartopy.feature as cfeature


lon = 10 + 15 * np.random.random(30)
lat = 55 + 15 * np.random.random(30)
data = np.random.randn(30)

land_10m = cfeature.NaturalEarthFeature('physical', 'land', '10m',
                                        edgecolor='k', facecolor='none')
borders_10m = cfeature.NaturalEarthFeature('cultural', 'admin_0_boundary_lines_land',
                                           '10m', edgecolor='grey', facecolor='none')
subborders_10m = cfeature.NaturalEarthFeature('cultural', 'admin_1_states_provinces',
                                              '10m', edgecolor='grey', alpha=0.3, facecolor='none')
                                              
                                              
projection = ccrs.TransverseMercator(central_longitude=15, scale_factor=0.9996) #SWEREF89
plate = ccrs.PlateCarree()

fig, ax = plt.subplots(1, 1, constrained_layout=True,
                     subplot_kw=dict(projection=projection),
                     figsize=(4, 8))


ax.add_feature(land_10m)
ax.add_feature(borders_10m)
ax.add_feature(subborders_10m)
c = ax.scatter(lon, lat, c=data, transform=plate,
           s=5, vmin=-2.5, vmax=2.5, cmap='bwr')
ax.set_title('Data')
plt.colorbar(c, ax=ax)

plt.show()