GeoMakie.jl: Stereographic projection at pole did not work

Following @j_u’s suggestion on one of the gotchas. Playing with the xlimits! and ylimits! make the plot appear.


using GeoMakie, CairoMakie

lons = -180:180
lats = -90:90
field = [exp(cosd(l)) + 3(y/90) for l in lons, y in lats]

fig = Figure()
ga = GeoAxis(
    fig[1, 1],
    dest="+proj=sterea +lat_0=90",    # Only lined that was changed compared to example!
    lonlims = automatic,
    coastlines = true,
    title = "Stereographic projection"
)

xlims!(ga, -90, 90)
ylims!(ga, -50, 50)
sp = surface!(ga, lons, lats, field; shading = false, colormap = :rainbow_bgyrm_35_85_c69_n256)
cb = Colorbar(fig[1, 2], sp)
fig

Produces:

1 Like