NAME¶

ascii2ps - create a postscript file of the data in the 3 column ASCII input file

DESCRIPTION¶

ascii2ps v1.1 Doris software, TUDELFT, (c)1999-2009

USAGE:

ascii2ps

3col_asciifile

Create a postscript file of the data in the 3 column ASCII input file containing "lon lat hei" triplets (output of lonlathei2ascii program). GMT is used to interpolate these irregularly distributed data to a grid.

This script is intended to be transparent and understandable. Please adapt this script to your own needs. (lambda=longitude). For example, here one could easily set the correct approximate height of the terrain.

The output is a postscript file named "temp.ps".

The steps foreseen to be taken are:

[before running this, create a three column ascii file with lon lat hei

using philamh2ascii program. AND optionally convert this 3 column file to UTM coordinates using ascii2ascii_UTM

] 2) Call GMT program xyz2grd to perform interpolation of this data.

the number of points in the output grid is 300.

3) call GMT program cpt to create a color table 4) call GMT program to create a postscript file 5) call ghostview to view output.

This script can serve as an example.

If you want to have another resolution

output grid, for example exactly 100x100 meters, then you will have to change this script. With the public domain software PROJ, you can easily convert the WGS coordinates (Doris output) to X,Y values, and then use this script. See the script ascii2ascii_UTM for an example how to do this for WGS to UTM coordinates conversion.

DEPENDENCIES:

GMT tools (gmtset, minmax, xyz2grd, grd2cpt, psbasemap

psscale, grdgradient, psmask, grdimage, pscoast)

standard UNIX (echo, csh, rm, etc.) ghostview

EXAMPLE:

ascii2ps lonlathei.dat

ascii2ps xyz_UTM.dat

SEE ALSO:

cpxfiddle, http://gmt.soest.hawaii.edu/

USAGE:

ascii2ps

3col_asciifile

Create a postscript file of the data in the 3 column ASCII input file containing "lon lat hei" triplets (output of lonlathei2ascii program). GMT is used to interpolate these irregularly distributed data to a grid.

This script is intended to be transparent and understandable. Please adapt this script to your own needs. (lambda=longitude). For example, here one could easily set the correct approximate height of the terrain.

The output is a postscript file named "temp.ps".

The steps foreseen to be taken are:

[before running this, create a three column ascii file with lon lat hei

using philamh2ascii program. AND optionally convert this 3 column file to UTM coordinates using ascii2ascii_UTM

] 2) Call GMT program xyz2grd to perform interpolation of this data.

the number of points in the output grid is 300.

3) call GMT program cpt to create a color table 4) call GMT program to create a postscript file 5) call ghostview to view output.

This script can serve as an example.

If you want to have another resolution

output grid, for example exactly 100x100 meters, then you will have to change this script. With the public domain software PROJ, you can easily convert the WGS coordinates (Doris output) to X,Y values, and then use this script. See the script ascii2ascii_UTM for an example how to do this for WGS to UTM coordinates conversion.

DEPENDENCIES:

GMT tools (gmtset, minmax, xyz2grd, grd2cpt, psbasemap

psscale, grdgradient, psmask, grdimage, pscoast)

standard UNIX (echo, csh, rm, etc.) ghostview

EXAMPLE:

ascii2ps lonlathei.dat

ascii2ps xyz_UTM.dat

SEE ALSO:

cpxfiddle, http://gmt.soest.hawaii.edu/

SEE ALSO¶

The full documentation for ascii2ps is maintained as a Texinfo manual. If the info and ascii2ps programs are properly installed at your site, the command

info ascii2ps

should give you access to the complete manual.