next up previous contents index
Next: Linear Projection with Polar Up: Non-map Projections Previous: Power projection

      
Geographical linear projection


  
Figure 5.4: Linear transformation of map coordinates
\begin{figure}\centering\epsfig{figure=eps/GMT_linear_d.eps}\end{figure}

While these linear projections are primarily designed for generic x,y data, it is sometimes necessary to plot geographical data in a linear projection. This poses a problem since longitudes have a 360o periodicity. GMT therefore needs to be informed that it has been given geographical data although a linear projection has been chosen. We do so by appending a d (for degrees) to the end of the -–Jx (or -–JX) option. As an example, we want to plot a crude world map centered on 125oE. Our command will be





gmtset GRID_CROSS_SIZE 0.1i BASEMAP_TYPE FANCY DEGREE_FORMAT 3
pscoast -R-55/305/-90/90 -Jx0.014id -B60g30f15/30g30f15WSen -Dc -A1000 -G200 -W0.25p -P \
    >! GMT_linear_d.ps
gmtset GRID_CROSS_SIZE 0





with the result reproduced in Figure 5.4.

     



Paul Wessel
1999-07-01