Python3×地理空間データ　GDAL Python API 【未完】

概要

Datasetが含むバンドの数

出力

int

サンプル使用例

dataset.RasterCount

　1　

標高値の1バンドしか入ってません．

概要

幅（東西方向のピクセル数）

出力

int

サンプル使用例

dataset.RasterXSize

3032

概要

高さ（南北方向のピクセル数）

出力

int

サンプル使用例

dataset.RasterXSize 

4063

XSize(3032)×YSize(4063)=12,319,016個のセル数ってことですね．
ビックデータだあぁ

概要

空間参照系に関する情報を取得する．

出力

string（以下のフォーマットのいずれか）

• Well Known Text (WKT)
• Human-Readable OGC WKT
• Proj4
• OGC WKT
• JSON
• GML
• ESRI WKT
• .PRJ File
• USGS
• MapServer Mapfile | Python
• Mapnik XML | Python
• GeoServer
• PostGIS spatial_ref_sys INSERT statement
• Proj4js format

フォーマットについては「spatialreference.org」で登録・検索・閲覧ができます．

サンプル使用例

dataset.GetProjection()

'LOCAL_CS["JGD2000 / Japan Plane Rectangular CS V",
             GEOGCS["JGD2000", 
                     DATUM["unknown", SPHEROID["unretrievable - using WGS84", 6378137,298.257223563], TOWGS84[0,0,0,0,0,0,0]],
                    PRIMEM["Greenwich",0],
                      UNIT["degree",0.0174532925199433]], 
          AUTHORITY["EPSG","2447"],
               UNIT["metre",1, AUTHORITY["EPSG","9001"]]]'

上の結果は，入れ子構造が見やすいように少々加工してます．

ちなみに，このようなデータをWKT（well-known text）と呼びます．
WKTについては，このサイトなどに詳しく載っています．

入れ子構造の上位から見ていきましょう．

• LOCAL_CS
  ローカル測地系（Local Coordinate System）であることを示しています．
  地図作成や測量などにおいて共通的に使用できる経緯度や標高を定義するための基準を測地系といいます．
  この測地系には，特定地域のみで適用できるローカル測地系と，世界的範囲で適用できる世界測地系があります．
  サンプルデータは日本測地系であり，ローカル測地系のひとつです．
  この話は，さすがのEsriさんがわかりやすいので詳細はお譲りします．
  空間参照系全般については，ここがわかりやすいです．
  　
• JGD2000 / Japan Plane Rectangular CS V
  「JGD2000」は日本測地系を，「Japan Plane Rectangular CS V」は平面直角座標系の第５系であることを意味しています．
  詳細は，ここやそこが詳しいです．
  　
• GEOGCS
  地理座標であることを示しています．
  投影座標であれば「PROJCS」，地心から見た座標なら「GEOCCS」になります．
  ここに載っています．
  　
• DATUM
  水平方向の測地基準を示します．
  　
• SPHEROID
  選択されている回転楕円体に関する情報を示します．
  ここらへんの話です．
  　
• TOWGS84
  WGS84以外の測地系からWGS84に変換するためのパラメータです．
  ここに詳しく載っていました．
  　
• PRIMEM
  基準にする経線を定義しています．
  このデータでは「グリニッジ子午線（Greenwich）」です．
  また，数値が正の値ならグリニッジより東に，負の値なら西に基準があることを示します．
  ここを参考にしました．
  　
• AUTHORITY
  このサイトによるとオプション情報を格納できる場所とのことですが，ほとんどの場合はEPSGコードが入っています．
  EPSGコードとは，各座標参照系に割り振られたIDです．
  このサンプルの場合は「2447」 が「JGD2000 / Japan Plane Rectangular CS V」に対応しています．
  よく使う空間参照系のEPSGは覚えておくと便利です．
  日本だけならここが詳しいです．
  EPSGそのものについてはこのサイトに載ってます．
  　
• UNIT
  単位です．
  「metre」ならばメートル単位の距離，「degree」なら角度です．

概要

空間参照系に関する情報を設定します．

入力

String（以下のフォーマットのいずれか）

フォーマットの記述法の詳細はGetProjection()やspatialreference.orgを見てください．

サンプル使用例

# 更新モードで読み込んでいるものとする
wkt = 'GEOGCS["JGD2000",DATUM["unknown",SPHEROID["unnamed",6378137,298.2572221010042],TOWGS84[0,0,0,0,0,0,0]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4612"]]'
dataset.SetProjection(wkt)
dataset.FlushCache()       #変更の保存

概要

地図画像の情報（座標値など）を取得します．

出力

tuple
(num, num, num, num, num, num)

中身は以下の順番です．

1. 画像左上セルのX座標
2. 単一セルのXサイズ（東ー西）
3. 回転度（上が北なら0）
4. 画像左上セルのY座標
5. 回転度（上が北なら0）
6. 単一セルのYサイズ（北ー南）

画像は基本的に左上端を基準にします．

使用例

dataset.GetGeoTransform()

(-65112.075071631014,
 11.418434105844657,
 0.0,
 -129197.91163072696,
 0.0,
 -11.418434105844657)

概要

地図画像の情報（座標値など）を設定します．

入力

tuple
(num, num, num, num, num, num)

中身の詳細はGetGeoTransform()を見てください．

使用例

# 更新モード
gt = (-65112.075071631014, 11.418434105844657, 0.0, -129197.91163072696, 0.0, -11.418434105844657)
dataset.SetGeoTransform(gt)
dataset.FlushCache()

概要

データセットのファイル名（フルパス）を取得します．

出力

string

使用例

dataset.GetDescription()

'C:\\hogehoge\\raster.tif'

概要

メタデータ（補助的な情報）を取得します．

出力

dict

代表的なGDALメタデータは以下の通りです．
（ このページ より）

• AREA_OR_POINT：「Area」（デフォルト）または「Point」のいずれか．
• NODATA_VALUES：バンドごとのnodataの値のリスト．

また，Geotiffデータセットのメタデータはこのページに載っています．

サンプル使用例

dataset.GetMetadata()

{'AREA_OR_POINT': 'Area'}

概要

メタデータ（補助的な情報）を設定します．

入力

dict

詳細はGetMetaData()をご覧ください．

サンプル使用例

# 更新モード
dataset.SetMetadata({"TIFFTAG_DOCUMENTNAME": "metameta"})
dataset.FlushCache()

概要

データセットに含まれるバンドを取得します．

入力

int

取得したいバンドの番号（1~）です．

出力

bandオブジェクト

サンプル使用例

band1 = dataset.GetRasterBand(1)
type(band1)

osgeo.gdal.Band

概要

全てのバンドの配列部分をnumpy配列（ndarray）として取得します．

マルチバンドデータの場合，すべての配列を一度に取得すると非常に重くなるので，通常は「GetRasterBand()」で取得したバンドに対して「ReadAsArray()」を適用します．

出力

ndarray

サンプル使用例

dataset.ReadAsArray()

array([[0., 0., 0., ..., 0., 0., 0.],
       [0., 0., 0., ..., 0., 0., 0.],
       [0., 0., 0., ..., 0., 0., 0.],
       ...,
       [0., 0., 0., ..., 0., 0., 0.],
       [0., 0., 0., ..., 0., 0., 0.],
       [0., 0., 0., ..., 0., 0., 0.]], dtype=float32)

単バンドデータなので，1つの配列が取得できました．

概要

全てのバンドのラスタデータをバイト列で取得します．

出力

サンプル使用例

• 概　要：
• 返り値：
• サンプル使用例