Scatter plot over a map

New to Julia
,
1

Hello,

I am trying to use GMT.jl to plot some points overtop a map of a region of the U.S… I can get the map to show, but when I try to scatter points over it, I’m getting an error message. Any help here would be much appreciated – I am new to Julia, thank you!

function example()
    coast(region=[-123.4, -121, 36.7, 38.3], proj=:Mercator, frame=:g, area=10000, land=:gray, water=:white, show=true)

    longitudes = [-122.31441247906619, -122.09929770009289, -121.30891030707745]
    latitudes = [37.395990554112785, 37.37364347711264, 37.18752179269895]

    GMT.scatter!(
        longitudes, 
        latitudes, 
        fmt=:png, 
        marker=:circle,
        markeredgecolor=:black,
        size=0.1,  
        markerfacecolor=:cyan,
        show=true
    )
end

The error message is:
ERROR: Something went wrong when calling the module. GMT error number = 79

2

You are using the show option twice when you should use it only at the last command. The scatter! call means append to a previous plot but there is no previous one because you finished it when you used show in the coast call.

This works:

coast(region=[-123.4, -121, 36.7, 38.3], proj=:Mercator, frame=:g, area=10000,
      land=:gray, water=:white)
scatter!([-122.31, -122.1, -121.3], [37.396, 37.37, 37.18], marker=:circle,
         markeredgecolor=:black, markerfacecolor=:cyan, size=0.1, show=true)

GMTjl_j
GMTjl_j1790×1500 42.2 KB

3

Thank you so much!