如何在 Julia 中塑造网格单元的轮廓？

问：

我想在世界地图上绘制一些区域并绘制它们的轮廓，因为它们中的一些区域彼此重叠。

第一次尝试：等值线图

function scale4plot(A,f=100)
  resize(A,f) = [A[i,j] for i=1 .+(0:f*size(A,1)-1).÷f,j=1 .+(0:f*size(A,2)-1).÷f]
  B = resize(A,f)
  ax,ay = ( (axes(B,i)./f).+0.5 for i=1:2)
  ay,ax,B
end

plot(Plots.contour(scale4plot(dplot[:,:,i]), w=1, linecolor=:black, fill=true )))

但是我不知道如何在上面绘制地图。

第二次尝试使用 GMT：GMT 图

G2 = mat2grid(dplot)
grdimage(G2)    
coast!(region=[0 360 -90 90],show=true)  

但在这里，我无法弄清楚如何勾勒出区域轮廓或降低不透明度。我也无法弄清楚如何使用 grdimage 将多个图绘制在一个图中。