Mixing UTM and geographic axes annotations
This is a question that comes up regularly such that GMT has a gallery example (famous example 28) to show how to do it. But even with that example it's not a trivial thing to do.
However it's trivial from the Julia wrapper under the condition that the grid or image to be displayed is referenced internally. If it is not,
grdedit can be used to assign a referencing system via its option
proj. The other condition is that grid or image to be displayed is already in memory (so that the internal magicks can work). Given that, example 28) boils down to this (with a little bit less fancy details).
using GMT G = gmtread("@Kilauea.utm.nc"); C = makecpt(cmap="copper", range=(0,1500)); imshow(G, cmap=C, shade=true, frame=(axes="WS", annot=true), coast=(shore=true, ocean=:lightblue, frame=(axes="EN", annot=true, grid=true)))
© GMT.jl. Last modified: January 20, 2023. Website built with Franklin.jl and the Julia programming language.
These docs were autogenerated using GMT: v0.44.4