gdallocationinfo

ラスタークエリツール

概要

Usage: gdallocationinfo [--help] [--help-general]
                        [-xml] [-lifonly] [-valonly]
                        [-E] [-field_sep <sep>] [-ignore_extra_input]
                        [-b <band>]... [-overview <overview_level>]
                        [-r {nearest|bilinear|cubic|cubicspline}]
                        [[-l_srs <srs_def>] | [-geoloc] | [-wgs84]]
                        [-oo <NAME>=<VALUE>]... <srcfile> [<x> <y>]

説明

gdallocationinfo ユーティリティは,さまざまな座標系の中での位置を指定することで,ピクセルに関する情報を問い合わせるメカニズムを提供します. いくつかのレポートオプションが提供されます.

--help: このヘルプメッセージを表示して終了します

--help-general: 一般的な GDAL コマンドラインオプションの使用方法を簡単に表示して終了します.

-xml: 便利な後処理のためのXML 形式で出力レポートが作成されます.

-lifonly: 出力は,データベースに対する LocationInfo リクエストから生成されたファイル名のみです (すなわち, VRT から影響を受けるファイルを識別するため).

-valonly: 選択されたピクセルの各選択されたバンドのピクセル値のみが出力されます. デフォルトでは,各バンドの値は別々の行に出力されますが, -field_sep が指定されている場合は別のフィールドに出力されます.

-b <band>: クエリするバンドを選択します. 複数のバンドをリストアップすることができます. デフォルトではすべてのバンドがクエリされます.

-r {nearest|bilinear|cubic|cubicspline}

Added in version 3.10.

サンプリングアルゴリズムを選択します. デフォルトは nearest です.

利用可能なメソッドは次のとおりです:

nearest で最近傍を利用できます.

bilinear でバイリニア畳み込みカーネルを利用できます.

cubic キュービック畳み込みカーネルを利用できます.

cubicspline B-Spline畳み込みカーネルを利用できます.

-overview <overview_level>: 基本バンドの代わりに (overview_level) 番目のオーバービュー (overview_level=1 は最初のオーバービュー) をクエリします. x,y 位置 (座標系がピクセル/ラインの場合) は基本バンドに対して引き続き指定する必要があります.

-l_srs <srs_def>: 入力 x, y 位置の座標系です.

-geoloc: 入力 x, y ポイントが画像のジオリファレンスシステムにあることを示します.

-wgs84: 入力 x, y ポイントが WGS84 経度, 緯度であることを示します.

-oo <NAME>=<VALUE>: データセットオープンオプション (フォーマット固有)

-ignore_extra_input: Added in version 3.9.

入力行の末尾に追加される余分な非数値コンテンツを出力行に追加しないようにするには,このフラグを設定します (-field_sep も定義されている必要があります). または,デフォルトまたは -xml モードで専用のフィールドとして使用できます.

-E: Added in version 3.9.

入力座標が -valonly モードで出力行に先行して表示されるエコーモードを有効にします.

-field_sep <sep>: Added in version 3.9.

異なる値を区切るために -valonly モードで使用されるフィールドセパレータを定義します. これはデフォルトで新しい行文字に設定されており, 複数のバンドを持つラスタをクエリすると, 出力には1行に1つの値が含まれるため, 複数の入力 x, y ポイントが提供された場合, どの値がどの入力 x, y ポイントのセットに属するかを認識するのが難しいかもしれません. フィールドセパレータを定義することも必要です.

<srcfile>: ソース GDAL ラスターデータソース名です.

<x>: ターゲットピクセルの X 位置です. デフォルトでは,座標系は -l_srs, -wgs84 または -geoloc が指定されていない限り,ピクセル/ラインです.

<y>: ターゲットピクセルの Y 位置です. デフォルトでは,座標系は -l_srs, -wgs84 または -geoloc が指定されていない限り,ピクセル/ラインです.

このユーティリティは,ピクセルに関するさまざまな情報を提供することを意図しています. 現在,次の情報を報告します:

ピクセルの位置をピクセル/ライン空間で指定します.
データソースに対する LocationInfo メタデータクエリの結果です. これは,そのピクセルに対するリクエストを満たすために使用されるファイルを報告する VRT ファイルと, MBTiles ドライバに実装されています.
すべてまたは一部のバンドのそのピクセルのラスタピクセル値です.
バンドにスケールと/またはオフセットが適用されている場合のスケーリングされていないピクセル値です.

選択されたピクセルは,コマンドラインで x/y 座標でリクエストされるか, stdin から読み取られます. stdin から座標を読み取る場合, 複数の座標ペアを提供することができます. デフォルトでは,整数ピクセル/ライン座標が想定されています. ただし, -geoloc, -wgs84, または -l_srs スイッチを使用すると, 他の座標系で位置を指定することができます.

デフォルトのレポートは人間が読めるテキスト形式です. 代わりに -xml スイッチで xml出力をリクエストすることも可能です.

スクリプト目的のために, -valonly と -lifonly スイッチが提供されており, 実際のピクセル値またはピクセルに識別された LocationInfo ファイルに出力を制限することができます.

将来, gdallocationinfo に追加のレポート機能が追加されることが予想されています.

例

ファイル utm.tif のピクセル (256,256) について報告する単純な例です.

$ gdallocationinfo utm.tif 256 256
Report:
Location: (256P,256L)
Band 1:
    Value: 115

WGS84 で位置を提供し, 結果を xml で取得する VRT ファイルをクエリします.

$ gdallocationinfo -xml -wgs84 utm.vrt -117.5 33.75
<Report pixel="217" line="282">
    <BandReport band="1">
        <LocationInfo>
        <File>utm.tif</File>
        </LocationInfo>
        <Value>16</Value>
    </BandReport>
</Report>

stdin から位置を読み取ります.

$ cat coordinates.txt
443020 3748359
441197 3749005
443852 3747743

$ cat coordinates.txt | gdallocationinfo -geoloc utmsmall.tif
Report:
  Location: (38P,49L)
  Band 1:
    Value: 214
Report:
  Location: (7P,38L)
  Band 1:
    Value: 107
Report:
  Location: (52P,59L)
  Band 1:
    Value: 148

$ cat coordinates.txt | gdallocationinfo -geoloc -valonly -E -field_sep , utmsmall.tif
443020,3748359,214
441197,3749005,107
443852,3747743,148