Next: 7.3 Spectral estimation and Up: 7. Creating GMT Graphics Previous: 7.1 The making of Contents Index

#
7.2 Image presentations

As our second example we will demonstrate how to make color images from gridded data sets (again, we will defer the actual making of grid files to later examples). We will use the supplemental program grdraster to extract 2-D grid files of bathymetry and Geosat geoid heights and put the two images on the same page. The region of interest is the Hawaiian islands, and due to the oblique trend of the island chain we prefer to rotate our geographical data sets using an oblique Mercator projection defined by the hotspot pole at (68°W, 69°N). We choose the point (190°, 25.5°) to be the center of our projection (e.g., the local origin), and we want to image a rectangular region defined by the longitudes and latitudes of the lower left and upper right corner of region. In our case we choose (160°, 20°) and (220°, 30°) as the corners. We use grdimage to make the illustration:

    #!/bin/sh
    #               GMT EXAMPLE 02
    #
    # Purpose:      Make two color images based gridded data
    # GMT progs:    gmtset, grd2cpt, grdgradient, grdimage, makecpt, psscale, pstext
    # Unix progs:   rm
    #
    ps=example_02.ps
    gmtset HEADER_FONT_SIZE 30 OBLIQUE_ANNOTATION 0
    makecpt -Crainbow -T-2/14/2 > g.cpt
    grdimage HI_geoid2.nc -R160/20/220/30r -JOc190/25.5/292/69/4.5i -E50 -K -P \
            -U/-1.25i/-1i/"Example 2 in Cookbook" -B10 -Cg.cpt -X1.5i -Y1.25i > $ps
    psscale -Cg.cpt -D5.1i/1.35i/2.88i/0.4i -O -K -Ac -B2:GEOID:/:m: -E >> $ps
    grd2cpt HI_topo2.nc -Crelief -Z > t.cpt
    grdgradient HI_topo2.nc -A0 -Nt -GHI_topo2_int.grd
    grdimage HI_topo2.nc -IHI_topo2_int.grd -R -J -E50 -B10:."H@#awaiian@# T@#opo and @#G@#eoid:" -O -K \
            -Ct.cpt -Y4.5i --HEADER_OFFSET=0.5i >> $ps
    psscale -Ct.cpt -D5.1i/1.35i/2.88i/0.4i -O -K -I0.3 -Ac -B2:TOPO:/:km: >> $ps
    pstext -R0/8.5/0/11 -Jx1i -O -N -Y-4.5i >> $ps << END
    -0.4 7.5 30 0.0 1 CB a)
    -0.4 3.0 30 0.0 1 CB b)
    END
    rm -f HI_topo2_int.grd ?.cpt .gmt*

The first step extracts the 2-D data sets from the local data base using grdraster, which is a supplemental utility program (see Appendix A) that may be adapted to reflect the nature of your data base format. It automatically figures out the required extent of the region given the two corners points and the projection. The extreme meridians and parallels enclosing the oblique region is -R159:50/220:10/3:10/47:35. This is the area extracted by grdraster. For your convenience we have commented out those lines and provided the two extracted files so you do not need grdraster to try this example. By using the embedded grid file format mechanism we saved the topography using kilometers as the data unit. We now have two grid files with bathymetry and geoid heights, respectively. We use makecpt to generate a linear color palette file geoid.cpt for the geoid and use grd2cpt to get a histogram-equalized cpt file topo.cpt for the topography data. To emphasize the structures in the data we calculate the slopes in the north-south direction using grdgradient; these will be used to modulate the color image. Next we run grdimage to create a color-code image of the Geosat geoid heights, and draw a color legend to the right of the image with psscale. Similarly, we run grdimage but specify -Y4.5i to plot above the previous image. Adding scale and label the two plots a) and b) completes the illustration (Figure 7.2).

Figure 7.2: Color images from gridded data.

Next: 7.3 Spectral estimation and Up: 7. Creating GMT Graphics Previous: 7.1 The making of Contents Index