Images clipped by coastlines

    geoid_cpt = grd2cpt("@india_geoid.nc", cmap=:rainbow)
    grdimage("@india_geoid.nc", I=:default, proj=:merc, figsize=16)

    # Then use gmt pscoast to initiate clip path for land
    coast!(region="@india_geoid.nc", res=:low, clip=:land)

    # Now generate topography image w/shading
    gray = makecpt(color=150, range=(-10000,10000), no_bg=:true)
    grdimage!("@india_topo.nc", I=:default, cmap=gray)

    # Finally undo clipping and overlay basemap
    coast!(clip=:end, frame=(annot=10, ticks=5, title="Clipping of Images"))

    # Put a color legend on top of the land mask
    colorbar!(cmap=geoid_cpt, pos=(anchor=:TR, inside=true, length=(10,0.5),
              offset=(0.8,0.25), horizontal=true), xaxis=(annot=5,ticks=1),
              yaxis=(label="m",), shade=true, Vd=1)

    # Add a text paragraph
    t = ["> 90 -10 12p 3i j"
        "@_@%5%Example 17.@%%@_  We first plot the color geoid image"
        "for the entire region, followed by a gray-shaded @#etopo5@#"
        "image that is clipped so it is only visible inside the coastlines."]
    text!(text_record(t), paragraph=true, fill=:white, pen=:thinner, clearance=(round=true,),
          offset=(-0.25,0.25), font=(12,"Times-Roman"), justify=:RB, show=true)