Spatial analysis for a list of coordinates

Q&A
2

Hi, this is how I would probably do it with GMT.jl
It takes a polygon (the circle in your case) and uses gmtselect (or similar) to find the x,y,z grid points that are inside that polygon and than do some stats with them.The following example computes the average height of Swisserland as given by the 6m global grid.

using GMT

julia> G = gmtread("@earth_relief_10m");

julia> suiss= coast(DCW="CH", dump=true, minpts=300);

julia> zonal_statistics(G, suiss, mean)

BoundingBox: [1207.8861083984375, 1207.8861083984375]

1×1 GMTdataset{Float64, 2}
 Row │       X
─────┼─────────
   1 │ 1207.89