Je fais une copie de cette page en l'adaptant au cas de l'Europe.

Nous souhaitons placer quelques villes :

data = [
  {"nom":"Paris", "lat":48.866667, "lng": 2.333333},
  {"nom":"Londres", "lat":51.507222, "lng":0.1275},
  {"nom":"Berlin", "lat":52.5166667, "lng":13.383}
]

import matplotlib.pyplot as plt
import matplotlib.image as mpimg
img = mpimg.imread('europe.png')
plt.imshow(img)
plt.show()

La carte choisie a le format 857×712. Le coin supérieur gauche a les coordonnées (0;0) et le coint inférieur droit (856;711). Rien à voir avec latitude et longitude… Nous devons créer une fonction (lng;lat) → (x;y).

La carte n'est pas en projection plate, il faut faire quelques calculs. Je vous donne la fonction.

def latlng_to_xy(lat,lng):
    xA = 770
    xB = 179
    lngA = 8.226100
    lngB = -1.784251
    a = (xA - xB) / (lngA - lngB)
    b = xA - a*lngA
    return int(a*lng + b) # on approxime à l'entier

Le principe est exactement le même.

def lat_to_y(lat):
    yA = 190
    yB = 685
    latA = 48.968311
    latB = 43.372928
    a = (yA - yB) / (latA - latB)
    b = yA - a*latA
    return int(a*lat + b) # on approxime à l'entier

Il suffit de générer une liste pour les x et une autre pour les y :

data = [
  {"nom":"Paris", "lat":48.866667, "lng": 2.333333},
  {"nom":"Marseille", "lat":43.300000, "lng":5.400000},
  {"nom":"Berlin", "lat":52.520007, "lng":13.404954}
]

xs = [lng_to_x(item["lng"]) for item in data]
ys = [lat_to_y(item["lat"]) for item in data]

plt.scatter(xs,ys)

import matplotlib.pyplot as plt
import matplotlib.image as mpimg

def lng_to_x(lng):
    xA = 770
    xB = 179
    lngA = 8.226100
    lngB = -1.784251
    a = (xA - xB) / (lngA - lngB)
    b = xA - a*lngA
    return int(a*lng + b) # on approxime à l'entier

def lat_to_y(lat):
    yA = 190
    yB = 685
    latA = 48.968311
    latB = 43.372928
    a = (yA - yB) / (latA - latB)
    b = yA - a*latA
    return int(a*lat + b) # on approxime à l'entier

data = [
  {"nom":"Paris", "lat":48.866667, "lng": 2.333333},
  {"nom":"Marseille", "lat":43.300000, "lng":5.400000},
  {"nom":"Berlin", "lat":52.520007, "lng":13.404954}
]

xs = [lng_to_x(item["lng"]) for item in data]
ys = [lat_to_y(item["lat"]) for item in data]
img = mpimg.imread('france_vide.png')
plt.scatter(xs,ys)
plt.imshow(img)
plt.show()