You can do it with GMT
using GMT
pts = rand(10,2);
c = convexhull(pts, gdataset=true); # gdataset option is optional but here avoids unnecessary conversions
gdalwrite(c, "lixo.shp") # Save it as a shape
# Confirm
imshow(gmtread("lixo.shp"), plot=(data=pts, marker=:star))
