I’m trying to map a rectangular lat/lon polygon on a non-rectangular map projection using GeoMakie,
When I use four corner points to define the region, GeoMakie draws straight lines between them. The example below shows what I am getting. I would like the polygon to follow the curvature of the meridians instead.
using Makie, CairoMakie
using GeoMakie
using Makie, CairoMakie
using GeoMakie
fig = Figure()
ax = GeoAxis(fig[1, 1], title="poly = [(60, -30), (180, -30), (180, 30), (60, 30)]")
poly = [(60, -30), (180, -30), (180, 30), (60, 30)]
poly!(ax, poly, color=:red)
fig
The default projection for GeoMakie is Equal Earth as I understand it, and the source default is Plate Caree, so I thought this would work. The corner points are obviously projected properly on the map.
What am I missing? Is there some way to accomplish what I want?