【gnuplot5】色の付いた世界地図を書いてみたい

ホント感動した・・・
こんな図がgnuplotでも作れる
(正確にはシェルコマンドで綺麗に整形したデータをプロットしてるだけだけど・・・そこは野暮か・・・)

■ 参考 : Plotting raster data from Natural Earth

このデータは convert_natural_earth.sh という(おそらく)自作のスクリプトで作られている
使い方も書いてあるけど、元の世界地図のTIFデータ(？)がないと使えない・・・・
■ 参考 : convert_natural_earth

上記のデータよりももうちょい解像度が高いデータはgithubにあった

■ 参考 : もうちょい解像度が高いデータ

これで書いたのが以下、

関係ありそうな記事リンク

■ 参考 : Plotting the world

■ 過去記事 : 【gnuplot】塗りつぶしプロット

■ 過去記事 : 【gnuplot】world.datを活用したい

(2017/05/22 追記)

TIFファイルという地図データは

■ 過去記事 : 1:50m Raster Data - Natural Earth

からダウンロードできるっぽい

上にリンクを貼ったデータはそれぞれ
・Cross-blended Hypsometric Tints
・Natural Earth 2
からダウンロードされたデータを加工したもの

加工するためのシェルスクリプト convert_natural_earth.sh をダウンロードしてきたTIFファイルにかけても、エラーが出ておかしなデータが出てきてしまう・・・
原因はMacのsedとLINUXのsedが違う挙動をすることで、

extra characters at the end of N command

みたいなエラーが出てくるが、
これはMacのsedで-iとするときはバックアップのファイルのための拡張子を指定しないといけないらしい
なので以下ではsed -i '.bak'と置き換えた

あとはpm3dのために、適時空行をいれないといけないが、元々のソースコードにあったsedではうまく動かなかったのでawkで代用

途中で257で割っているのは、RGBの値がきちんと取得できていないので、よくわからんが257で割っている
時間をかけたらうまくできそうだけど、めんどくさいからいいや・・・・

#!/bin/sh

# data : http://www.naturalearthdata.com/downloads/50m-raster-data/

if [ $# -ne 2 ]; then
echo "usage : `basename $0` [resolution] [inputfile]"
exit
fi

RES="$1"
IN="$2"
SCALED="${IN%.tif}_${RES}px.tif"
OUT="${SCALED%.tif}.txt"

# resize tiff image to desired resolution
convert "$IN" -resize $RES "$SCALED"

# convert to txt file
convert "$SCALED" "$OUT"

# remove unneeded header
sed -i '.bak' '1d' "$OUT"

# replace ": (" with ","
sed -i '.bak' 's/: (/\,/' "$OUT"

# remove all unneeded entries at end of line
sed -i '.bak' 's/).*$//' "$OUT"

# convert longitude and latitude to the right values
# 0..RES => -180..180
# 0..RES/2 => -90..89
awk -F, -v res=${RES} '{$1=$1*360/res-180; $2=-1*($2*360/res-90); }1' OFS=, "$OUT" > tmp.txt
awk -F, '{print $1, $2, $3/257, $4/257, $5/257}' OFS=, tmp.txt > tmp2.txt
mv tmp2.txt ${OUT}
rm tmp.txt

# add blank lines for pm3d
awk -v n=${RES} '{print $0} NR>2 && 0==NR % n{printf "\n"}' ${OUT} > tmp.txt
mv tmp.txt ${OUT}

# add new header
awk 'BEGIN{print "# long, lat, r, g, b"} {print $0}' ${OUT} > tmp.txt
mv tmp.txt ${OUT}

rm *.bak

これならMacでも動くはず・・・・

n=2000 でプロットしたのが以下の図
この図の上に色々なものを書いていきたい