×

[PR]この広告は3ヶ月以上更新がないため表示されています。
ホームページを更新後24時間以内に表示されなくなります。

Bird's-eye view
Kinai(Kinki district) for example.


1. Prepare grid file.


In this case, "japan250m.nan.grd" is used.

2. Prepare colour palette table (cpt file).


Type,
% makecpt -Cglobe -Z > cpt
This is to use colour palette "globe" .
To type only `makecpt' show digest of usage.
To type `man makecpt' shows manual.
In this case, I modify this "globe" cpt file, namely cptg. cptg

3. Prepare Intensity file.


% grdcut /somewhere/japan250m.nan.grd -Gtopo.grd \
-R134.5/136.5/33/36 -V
% grdgradient topo.grd -A30 -Ghst.grd -V -Ne0.8
Note too intense shadow make topographic impression abrupt.

4. Write a shell script.


% cat > kinai.sh
#! /bin/sh
grdfile=topo.grd
output=kinai-cho.ps
range="134.5/136.5/33/36"
gmtset PAPER_MEDIA a3+
grdview $grdfile \
-JX10/12 -JZ10 -B5f1g20m -R$range \
-Qi -E225/60 -Ccptg -Ihst.grd -P -V -Z1000 \
> $output

You need not writing shell script because this is too simple command, only one. However to write shell script will help you for further work.

5. Run the shell script.


% sh kinai.sh
To open output file, type
% gv kinai-cho.ps
or open with Adobe Illustrator.
Figure

Problem in this case


Lake Biwa is not drawn. Note lakes and ponds can not be drawn by "grdview".
Marsh of Osaka (e.g. Fukushima-ku) sinks under water. This is derived from error values of used grd file.