PNG標高タイル ( PNG Elevation tile )

　PNG標高タイルは，PNGを利用した，Webブラウザで使用するための標高タイルフォーマットです．

バージョン

1.0.2 ( 2016-11-29 )

仕様

タイルファイルフォーマット

　PNG(Portable Network Graphics) 24ビットまたは32ビットRGB，8ビットインデックスカラー
　GIF(Graphic Interchange format

各セルの値が示す内容

　各セルの値はセル北西端の値を表します．
　※セルの代表値（平均値や中央値など）ではありません．

タイル東端・南端の扱いとタイプ

　タイル東端・南端の扱いによって通常タイプと拡大タイプの2種類を規定します．各セルの値がセル北西端の値を表すため，通常はタイル画像東端および南端の標高値を表現できません．これを通常タイプと呼びます．画素数を縦，横１ピクセルずつ増やして東端および南端の標高値を表現するものを拡大タイプとします．　

標高データの取得方法

　セルの画素から標高値を取得するには，そのRGB値(r,g,b ; それぞれ256階調で表現したもの)をこの順に並べて24bit符号付き整数として解釈して整数値を得え，この値にデータセット毎に定めらえた標高分解能u(単位m)をかけて標高値h(単位m)を取得します．
　具体的な式は以下のようになります．　

x = 2¹⁶r + 2⁸ g + b

h =
　xu,
　(x - 2²⁴)u,
( x < 2²³ )
( x ≧ 2²³ )

　ただし，不透明度が0の場合は無効値(NA)を示すものとします．

解説

　標高タイルとは，標高データを保持するタイル状に分割された画像またはテキストなどのファイルを指します．標高データを画像のピクセル情報として保持する方法自体は新しいものではなく，GeoTiff画像などとしてよく用いられます．また，多量の標高データを扱う場合，それを同じサイズの「タイル」に分割することは自然な方法であり，これもしばしば用いられます．一方，標高タイルにおける画像フォーマットやセルへの情報の保存方法は，その目的によって様々なものが用いられます．ここでは，特に，多量の標高データをWebブラウザで高速に扱うために適したフォーマットを規定しています．　

応用範囲

　国土地理院の標高データを扱うために考案されましたが，極地域を含む日本以外や，様々な精度に対応できるように設計されています．
　このため，本仕様の中では以下の内容は規定されていません．

地図投影法，タイル分割法，タイル命名規則
１枚のタイルのセルサイズ
データ精度(有効桁数）

地図投影法，タイル分割法，タイル命名規則

　仕様上は規定していませんが，Googleマップ互換のWebメルカトル図法や命名規則，OpenLayers.XYZなどが汎用性があり，扱いやすいでしょう．　
　一方，Webメルカトル図法では極が表現できませんので，用途によっては正距円筒図法などの投影法を採用することになります．　

セルサイズ

　一般に縦横それぞれ256が用いられますが，この仕様上は規定していません．

標高分解能

　地理院地図の標高タイルでは標高分解能0.01m(1cm)でデータが保存されています．この標高分解能の場合，本フォーマットでは-83,886.07m～+83,886.07mが表現可能で，界最高標高のエベレスト(8,844m)や世界最深のマリアナ海溝チャレンジャー海淵(10,911m)まで表現可能です．標高分解能0.001m(1mm)で表現する場合は-8,388.607m～+8,388.607mが表現可能で，日本最高標高の富士山（3,776m）を表現することができます．逆に精度を落とした場合は，PNGの圧縮効率が上がり，ファイルサイズが小さくなることが期待できます．　

Webブラウザでの情報の取得方法

　Webブラウザでは，HTML5のCanvasおよびJavaScriptを使って標高データを取り出すことができます．ただし，クロスオリジンによる制約を回避するため，imageオブジェクトのcrossOriginプロパティに'anonymous'を設定する等の対応が必要になります．
　以下に，タイル内座標(i, j)の位置の標高データを取得するJavaScriptのサンプルコードを示します（標高分解能は1mとします）．

var
    image = new Image(),
    filename = 'https://gsj-seamless.jp/labs/elev2/elev/9/201/454.png',
    i = 0, j = 0;    // ピクセル位置を指定する

image.crossOrigin = 'anonymous';
image.onload = function() {
    var
        r, g, b, x, h,
        u = 1,		// 標高分解能1m
        canvas = document.createElement('canvas'),
        context = canvas.getContext( '2d' );

    canvas.width = this.width;
    canvas.height = this.height;  
    context.drawImage( this, 0, 0 );
    data = context.getImageData( i, j, 1, 1 ).data;
    if ( data[ 3 ] == 0 ) ) {
        console.log( '無効値' );
    } else {
        r = data[ 0 ]
        g = data[ 1 ]
        b = data[ 2 ]
        x = Math.pow( 2, 16 ) * r + Math.pow( 2, 8 ) * g + b;
        h = ( x < Math.pow( 2, 23 )  ) ? x * u: ( x - Math.pow( 2, 24 ) ) * u;
        console.log( h );
    }
}
image.src = filename;

サンプル

　国土地理院の提供する標高タイルデータをもとにサンプルファイルを作成しました．
　Googleマップ互換のタイル分割方法で東京付近(z=9, y=201, x=454)のPNG標高タイルです(標高分解能は1m)．

履歴

日付	バージョン	概要
2016/11/29	1.0.2	「Webブラウザでの情報の取得方法」におけるタイルパスの誤りを修正
2016/8/30	1.0.1	無効値に関する不透明度の表記の誤りを修正
2016/5/20	1.0.0	データタイルマップ仕様(ver.1.0.0)に適合するようにPNG標高タイル(ver 0.6.3)をもとに作成
2016/4/20	0.6.3
2014/12/19	0.6.0
2014/6/25	0.5.1
2014/5/19	0.3.0
2014/4/25	0.2.4

x =	2¹⁶r + 2⁸ g + b
h =	xu, 　(x - 2²⁴)u,	( x < 2²³ ) ( x ≧ 2²³ )