And for those wanting color
G = gmtread("alps.grd", layout="TRB");
G2 = GMT.texture_img(G);
hill = gdaldem("alps.grd", "hillshade", zfactor=4);
color = gdaldem("alps.grd", "color-relief", color="faa.cpt");
G3 = blendimg!(G2, hill, new=true);
G4 = blendimg!(color, G3, new=true);
imshow(G4, fmt=:png)
