0

I want to plot some areas on a world map and contour them because some of them are overlapping each other.

First attempt: contour plot

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 )))

But I couldn't figure out how to plot a map over it.

Second attempt with GMT: GMT plot

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

But here i couldn´t figure out how to outline the areas or decrease the opacity. I also couldn't figure out how to plot multiple plots in one with grdimage.

vcv
  • 11
  • 2

0 Answers0