Ok. It’s not that hard to do with GeoMakie.jl (well I did it during dinner 
), it will be better soon, after the pending release is out. (xticks, yticks could be also added manually.)
using GeoMakie, GLMakie
using GraphMakie, Graphs
using Proj4
using Downloads
using GeoJSON, GeoInterface
states = Downloads.download("https://raw.githubusercontent.com/PublicaMundi/MappingAPI/master/data/geojson/us-states.json")
states_geo = GeoJSON.read(read(states, String))
n = length(GeoInterface.features(states_geo))
#g = wheel_graph(10)
gpos = Dict(
(-96.5967, 38.9617) => "KSCYng",
(-122.3, 47.6) => "STTLng",
(-105.0, 40.75) => "DNVRng",
(-122.026, 37.3858) => "SNVAng",
(-87.6167, 41.8333) => "CHINng",
(-85.5, 34.5) => "ATLAng",
(-77.0268, 38.8973) => "WASHng",
(-73.9667, 40.7833) => "NYCMng",
(-86.1595, 39.7806) => "IPLSng",
(-95.5174, 29.77) => "HSTNng",
(-118.25, 34.05) => "LOSAng",
(-84.3833, 33.75) => "ATLAM5")
g = complete_graph(length(keys(gpos)))
positions = Point2f.(collect(keys(gpos)))
fig = Figure(resolution = (1200, 800))
ax = Axis(fig[1, 1], aspect = DataAspect())
trans = Proj4.Transformation("+proj=longlat +datum=WGS84",
"+proj=lcc +lon_0=-100 +lat_1=33 +lat_2=45", always_xy = true)
ptrans = Makie.PointTrans{2}(trans)
ax.scene.transformation.transform_func[] = ptrans
lons = -130:-70
lats = 24:52
points = [Point2f0(lon, lat) for lon in lons, lat in lats]
rectLimits = FRect2D(Makie.apply_transform(ptrans, points))
limits!(ax, rectLimits)
lines!(ax, GeoMakie.coastlines(), color = :black)
poly!(ax, states_geo, color = 1:n, colormap = (:viridis, 0.25), strokecolor = :black,
strokewidth = 1)
graphplot!(ax, g; layout = _ -> positions, node_size = 20,
edge_color = to_colormap(:plasma, 66),
node_color = to_colormap(:plasma, length(keys(gpos))))
# drawing the correct lines and ticks is still a little bit cumbersome and manual
# you could ignore the rest and just hide everything
lonrange = range(-130, -70, length = 5)
latrange = range(24, 52, length = 5)
lonlines = [Point2f0(j, i) for i in lats, j in lonrange]
latlines = [Point2f0(j, i) for j in lons, i in latrange]
[lines!(ax, lonlines[:, i], color = (:black, 0.25),
linestyle = :dash, overdraw = true) for i in 1:size(lonlines)[2]]
[lines!(ax, latlines[:, i], color = (:black, 0.25), linestyle = :dash,
overdraw = true) for i in 1:size(latlines)[2]]
hidedecorations!(ax)
hidespines!(ax)
save("graphMap.png", fig)
fig
(proj4Geo) pkg> st
Status `~/Desktop/proj4Geo/Project.toml`
[13f3f980] CairoMakie v0.6.6
[e9467ef8] GLMakie v0.4.7
[ddc7317b] GeoDatasets v0.1.6
[cf35fbd7] GeoInterface v0.5.6
[61d90e0f] GeoJSON v0.5.1
[db073c08] GeoMakie v0.2.2
[5c1252a2] GeometryBasics v0.4.1
[1ecd5474] GraphMakie v0.3.0
[86223c79] Graphs v1.5.1
[ee78f7c6] Makie v0.15.3
[85f8d34a] NCDatasets v0.11.7
[9a7e659c] Proj4 v0.7.6