国土地理院標高タイルの数値データのダウンロードと画像表示の手順を説明する.
「レイヤ番号/X/Y」
最も荒いものがレベル0.レベルが増えることに,範囲は狭く,情報は精密になる. レイヤ番号=0のときは全世界,レイヤ番号が1増えると,面積は 1/4
※ 画像の地図の場合には,1個の地図タイル = 256×256の画像が標準
オンラインの地図システムでは,ダウンロードを早くするために,地図データは地図タイル化
※ゲームの「マップ」でもタイル地図と同じ技術が使われていることも
タイル地図形式になっているような標高データのこと
国土地理院が公開しているタイル形式の地図.種類はたくさんあり,次のWebページで確認できる.
標高データ(単位はm)が入った 256 × 256 サイズの数値データ.国土地理院が公開.
http://maps.gsi.go.jp/development/demtile.html
※ 著作権は,国土地理院にあります.特に,他の人に配ったり,インターネットで公開したい 合などは,念入りに,利用条件http://www.gsi.go.jp/LAW/2930-index.html#sec5 をよく確認すること.
※ この資料では,学習目的に,国土地理院タイルデータを使用する.
キーワード: 地図タイル,標高タイル,国土地理院タイル地図,国土地理院標高タイル,高さマップ(Height map),ダウンロード
http://maps.gsi.go.jp/development/tileCoordCheck.html
ズームイン,ズームアウトの操作は,Webブラウザの画面の左上の「+」, 「-」ボタンを利用するか,マウスのホイールを使う.
※ 「好きな場所」は,自宅に近くなど,慣れている場所を選んでみる.ある程度起伏がある場所を選ぶのが楽しい.
(例)14/14255/6519
※ Python でプロットを行いたい.
※ そのためには, PyCharmなどにある Python コンソールが便利である.
「E:/6519.txt」のところは,先ほどのファイル名に読み替える.
import pandas as pd %matplotlib inline import matplotlib.pyplot as plt import warnings warnings.filterwarnings('ignore') # Suppress Matplotlib warnings Z = pd.read_csv("E:/6519.txt", header=None) plt.imshow(Z)
次のように表示されたら OK. 標高が高いところが赤,低いところが青. になるように、自動的に色がつく.
※ 元データは 0 から 255 の範囲の数字だけど、表示されるときに、赤から青の色が勝手に付く、という意味です.元データはあくまでも 0 から 255 の範囲の数字.
※ Python でプロットを行いたい.
※ そのためには, PyCharmなどにある Python コンソールが便利である.
「E:/6519.txt」のところは,先ほどのファイル名に読み替える.
import pandas as pd %matplotlib inline import matplotlib.pyplot as plt import warnings warnings.filterwarnings('ignore') # Suppress Matplotlib warnings Z = pd.read_csv("E:/6519.txt", header=None) Z.shape
下の実行結果では,サイズは 256 x 256
上のプログラムに続いて,次の Python プログラムを実行
from mpl_toolkits.mplot3d.axes3d import * fig = plt.figure() ax = Axes3D(fig) X, Y = np.meshgrid(np.linspace(0,1,256), np.linspace(0,1,256)) ax.plot_surface(X, Y, Z) plt.show()