next up previous contents index
Next: 4.3.1 Exercises Up: 4. SESSION FOUR Previous: 4.2 Illumination and intensities   Contents   Index


4.3 Color images

Once a cpt file has been made it is relatively straightforward to generate a color image of a gridded data. Here, we will extract a subset of the global 30" DEM (data id 9) from USGS:

grdraster 9 -R-108/-103/35/40 -Gus.grd

Using grdinfo we find that the data ranges from $\sim$1000m to $\sim$4300m so we make a cpt file accordingly:

makecpt -Crainbow -T1000/5000/500 -Z >! topo.cpt

Color images are made with grdimage which takes the usual common command options (by default the -R is taken from the data set) and a cptfile; the main other options are


Table 4.3: The main options in grdimage.
Option Purpose
-Edpi Sets the desired resolution of the image [Default is data resolution]
-Iintenfile Use artificial illumination using intensities from intensfile
-M Force grayshade using the (television) YIQ conversion


We want to make a plain color map with a color bar superimposed above the plot. We try

grdimage us.grd -JM6i -P -B2 -Ctopo.cpt -V -K >! topo.ps
psscale -D3i/8.5i/5i/0.25ih -Ctopo.cpt -I0.4 -B/:m: -O >> topo.ps

The plain color map lacks detail and fails to reveal the topographic complexity of this Rocky Mountain region. What it needs is artificial illumination. We want to simulate shading by a sun source in the east, hence we derive the required intensities from the gradients of the topography in the N90$^{o}$E direction using grdgradient. Other than the required input and output filenames, the available options are


Table 4.4: The grdgradient options.
Option Purpose
-Aazimuth Azimuthal direction for gradients
-M Indicates that this is a geographic grid
-N[t$\vert$e][norm[/offset]] Normalize gradients by norm/offset [= 1/0 by default].
Insert t to normalize by the $tan^{-1}$ transformation.
Insert e to normalize by the cumulative Laplace distribution.


Figure 4.1 shows that raw slopes from bathymetry tend to be far from normally distributed (left). By using the inverse tangent transformation we can ensure a more uniform distribution (right). The inverse tangent transform simply takes the raw slope estimate (the x value at the arrow) and returns the corresponding inverse tangent value (normalized to fall in the $\pm1$ range; horizontal arrow pointing to the y-value).

Figure 4.1: How the inverse tangent operation works.
\includegraphics{eps/GMT_atan}

Both -Ne and -Nt yield well behaved gradients. Personally, we prefer to use the -Ne option; the value of norm is subjective and you may experiment somewhat in the 0.5-5 range. For our case we choose

grdgradient us.grd -Ne0.8 -A100 -M -Gus_i.grd

Given the cpt file and the two gridded data sets we can create the shaded relief image:

grdimage us.grd -Ius_i.grd -JM6i -P -B2 -Ctopo.cpt -K >! topo.ps
psscale -D3i/8.5i/5i/0.25ih -Ctopo.cpt -I0.4 -B/:m: -O >> topo.ps



Subsections
next up previous contents index
Next: 4.3.1 Exercises Up: 4. SESSION FOUR Previous: 4.2 Illumination and intensities   Contents   Index
Paul Wessel 2006-05-31