I find these contour plots very cute, and now finally implemented it in GMT.jl
using GMT
G = grdcut("@mars_relief_02m", R="220/233/13/25");
contourf(G, tanaka=true, title="Tanaka Olympus", show=1)
5 Likes