これは、Ubuntu Online、Fedora Online、Windows オンライン エミュレーター、MAC OS オンライン エミュレーターなど、複数の無料オンライン ワークステーションのいずれかを使用して、OnWorks 無料ホスティング プロバイダーで実行できるコマンド r.in.gdalgrass です。
プログラム:
NAME
リングダル - GDAL ライブラリを使用してラスター データを GRASS ラスター マップにインポートします。
KEYWORDS
ラスター、インポート
SYNOPSIS
リングダル
リングダル - 助けて
リングダル [-ojeflkcp] =名 出力=名 [バンド=整数[,整数、...]]
[メモリ=整数] [ターゲット=名] [タイトル=フレーズ] [オフセット=整数] [場所=名]
[テーブル=file] [-上書きする] [-助けます] [-詳細] [-静かな] [-ui]
フラグ:
-o
投影法をオーバーライド (場所の投影法を使用)
-j
映り込みチェックだけして終了
-e
新しいデータセットに基づいて領域範囲を拡張する
PERMANENT マップセット内の場合はデフォルトのリージョンも更新します
-f
サポートされている形式をリストして終了します
-l
緯度/経度マップを地理座標 (90N,S; 180E,W) に合わせます
-k
バンドの色名を使用する代わりにバンド番号を保持する
-c
「location」パラメータで指定された場所を作成して終了します。 輸入しないでください
ラスター ファイル。
-p
バンド数を表示して終了
-上書き
出力ファイルが既存のファイルを上書きできるようにする
- 助けて
使用状況の概要を印刷する
-詳細
冗長モジュール出力
- 静かな
静かなモジュール出力
--ui
GUIダイアログを強制的に起動する
パラメーター:
=名 [必要]
インポートするラスター ファイルの名前
出力=名 [必要]
出力ラスターマップの名前
バンド=integer [、integer、...]
選択するバンド (デフォルトはすべてのバンド)
メモリ=整数
使用する最大メモリ(MB単位)
ラスター行のキャッシュサイズ
オプション: 0-2047
デフォルト: 300
ターゲット=名
GCP ターゲット ロケーションの名前
GCP 変換のためにプロジェクションを作成または読み取る場所の名前
タイトル=フレーズ
結果として得られるラスター マップのタイトル
オフセット=整数
オフセットは、出力ラスター マップ名の作成中にバンド番号に追加されます
デフォルト: 0
場所=名
作成する新しい場所の名前
テーブル=file
ラスター属性テーブルのファイル プレフィックス
バンド番号と「.csv」がファイルのプレフィックスに追加されます
DESCRIPTION
リングダル ユーザーは (バイナリ) GRASS ラスター マップ レイヤーまたは画像グループを作成できます。
オプションのタイトルを付けて、GDAL がサポートするラスター マップ形式から。 インポートされたファイルは
オプションで、新しい場所を作成するためにも使用できます。
GDAL 。 ラスター フォーマット
GDAL でサポートされている形式の詳細については、次の URL を参照してください。
http://www.gdal.org/formats_list.html
サポートされている 80 を超える形式の選択された形式:
ロング フォーマット名コード作成ジオリファレンス最大ファイル サイズ
Arc/Info ASCII グリッド AAIGrid はい はい 制限なし
Arc/Info Binary Grid AIG いいえ はい --
AIRSAR ポラリメトリック AIRSAR いいえ いいえ --
Microsoft Windows Device Independent Bitmap (.bmp) BMP はい はい 4GiB
BSB 海図形式 (.kap) BSB いいえ はい --
VTP Binary Terrain Format (.bt) BT はい はい --
CEOS (例えばスポット) CEOS No No --
第 1 世代 USGS DOQ (.doq) DOQXNUMX いいえ はい --
新しいラベル付き USGS DOQ (.doq) DOQ2 いいえ はい --
データ (.dt0、.dt1) DTED いいえ はい --
ERMapper 圧縮ウェーブレット (.ecw) ECW はい はい
ESRI .hdr ラベル付き EHdr いいえ はい --
ENVI .hdr ラベル付きラスター ENVI はい はい 制限なし
Envisat イメージ プロダクト (.n1) Envisat いいえ いいえ --
EOSAT FAST フォーマット FAST いいえ はい --
FITS (.fits) FITS はい いいえ
Graphics Interchange Format (.gif) GIF はい いいえ 2GB
Arc/Info Binary Grid (.adf) GIO はい はい
GRASS ラスター GRASS いいえ はい --
TIFF / GeoTIFF (.tif) GTiff はい はい 4GiB
Hierarchical Data Format Release 4 (HDF4) HDF4 はい はい 2GiB
Erdas Imagine (.img) HFA はい はい 制限なし
Atlantis MFF2e HKV はい はい 制限なし
画像の表示と分析 (WinDisp) IDA あり あり 2GB
ILWIS ラスター マップ (.mpr、.mpl) ILWIS はい はい --
日本語 DEM (.mem) JDEM いいえ はい --
JPEG JFIF (.jpg) JPEG はい はい 4GiB (最大寸法 65500x65500)
JPEG2000 (.jp2、.j2k) JPEG2000 はい はい 2GiB
JPEG2000 (.jp2、.j2k) JP2KAK はい はい 制限なし
NOAA Polar Orbiter Level 1b データセット (AVHRR) L1B いいえ はい --
Erdas 7.x .LAN および .GIS LAN いいえ はい 2GB
インメモリ ラスター MEM はい はい 2GiB
アトランティス MFF MFF はい はい 制限なし
マルチ解像度シームレス画像データベース MrSID なし あり --
NDF NLAPS データ形式 いいえ はい 制限なし
NITF NITF はい はい
NetCDF netCDF はい はい 2GB
OGDI ブリッジ OGDI いいえ はい --
PCI .aux ラベル付き PAux はい いいえ 制限なし
PCI Geomatics データベース ファイル PCIDSK はい はい 制限なし
ポータブル ネットワーク グラフィックス (.png) PNG はい いいえ
PCRaster (.map) PCRaster はい いいえ
Netpbm (.ppm、.pgm) PNM はい いいえ 制限なし
RadarSat2 XML (product.xml) RS2 いいえ はい 4GB
USGS SDTS DEM (*CATD.DDF) SDTS いいえ はい --
SAR CEOS SAR_CEOS いいえ はい --
USGS ASCII DEM (.dem) USGSDEM いいえ はい --
X11 Pixmap (.xpm) XPM はい いいえ
会社名 創造
リングダル 場合、データセットをインポートするときに投影情報を保持しようとします。
ソース形式には投影情報が含まれており、GDAL ドライバーがそれをサポートしているかどうか。 もし
ソース データセットの投影法が現在の場所の投影法と一致しない
リングダル エラー メッセージが報告されます (データセットの射影が一致していないようです
現在の場所)、ソース データセットの PROJ_INFO パラメータをレポートします。
ユーザーが見かけの座標系の違いを無視したい場合は、
ソースデータと現在の場所、それらは -o 投影をオーバーライドするフラグ
チェック。
ユーザーが完全な投影定義でデータをインポートしたい場合は、可能です。
投影と範囲に基づいて r.ing.dal に新しい場所を自動的に作成させる
読み取られているファイルの。 これは、新しい
経由の場所 場所 パラメータ。 コマンドが完了すると、新しい場所が
(PERMANENT マップセットのみで) 作成され、ラスターがインポートされます。
示された 出力 名前を PERMANENT マップセットに追加します。
サポート for GCP
画像に GCP が含まれている場合、GCP は画像グループ内の POINTS ファイルに書き込まれます。
これらは i.rectify に直接使用できます。
この ターゲット オプションを使用すると、GCP を独自のものから自動的に再投影できます
ロケーション名の PROJ_INFO ファイルから読み取った別の投影への投影
ターゲット.
Status ターゲット 場所が存在しない場合、一致する新しい場所が作成されます
GCP の射影定義。 出力グループのターゲットは、新しい
場所、および i.rectify は、それ以上の準備なしで使用できるようになりました。
一部の衛星画像 (NOAA/AVHRR、ENVISAT など) には、数百または数千の画像が含まれる場合があります。
GCP。 これらの場合、薄板スプライン座標変換が推奨されます。
でインポートする前に グダルワープ -tps またはインポート後 i.修正する -t.
注意事項
設定すると、大きなファイルのインポートが大幅に高速化されます メモリ のサイズに
入力ファイル。
この リングダル コマンドは、基礎となるものである限り、次の機能をサポートします
フォーマットドライバーはそれをサポートしています:
カラーテーブル
カラーテーブルが関連付けられているバンドには、カラーテーブルが転送されます。 ご了承ください
ソースにカラーマップがない場合、GRASS 5.0 の r.in.gdal はカラーマップを出力しません。 使用する
r.colors map=... color=grey でグレースケール カラーマップを割り当てます。 の将来のバージョンでは
GRASS r.in.gdal は、グレースケール カラーマップを自動的に出力するようにアップグレードされる可能性があります。
データ型
ほとんどの GDAL データ型がサポートされています。 Float32 および Float64 タイプのバンドは次のように変換されます。
GRASS 浮動小数点セル (ただし、倍精度ではありません ... これは、次の場合に追加できます
必要)、他のほとんどのタイプは GRASS 整数セルとして変換されます。 これも
16 ビット整数データ ソース。 複雑な (一部の SAR 信号データ形式) データ バンドは、
XNUMX つの浮動小数点セル レイヤー (*.real および *.imaginary) に変換されます。
ジオリファレンス
データセットにアフィン地理参照情報がある場合、これを使用して
北端、南端、東端、西端。 回転係数は無視され、結果として
回転したデータセットの位置が正しくありません。
投影
データセットの投影法は、現在の場所と比較したり、
新しい場所。 内部的には、GDAL は OpenGIS Well Known Text の射影を表します
フォーマット。 GRASS プロジェクションの全セットの大きなサブセットがサポートされています。
Null値
null 値が GDAL によって認識されるラスター バンドには、null ピクセルがあります。
インポート中に GRASS スタイルの null に変換されます。 多くの一般的な形式 (および形式)
GDAL では十分にサポートされていません) には null ピクセルを認識する方法がありません。
r.null は、インポート後に使用する必要があります。
GCP
地上基準点を含むデータセットは、POINTS ファイルとしてインポートされます。
画像グループに関連付けられています。 バンドが XNUMX つだけのデータセット
単純なラスター マップとして変換され、関連付けられた画像グループも含まれます。
地上基準点がある場合。 地上管制の座標系
ポイントはr.in.gdalによって報告されますが、保存されません。 確認するかどうかはユーザー次第です
i.target で確立された位置が以前に互換性のある座標系を持っていること
i.rectifyでポイントを使用します。
計画された改善 リングダル 将来的には、すべてを報告するためのサポートが含まれます
データセットについて知られている場合、 出力 パラメータが設定されていません。
エラー メッセージ
"エラー: 入力 地図 is 回転させました - 輸入。"
この場合、保存されている回転情報を適用して、画像を最初に外部回転する必要があります。
ラスター イメージ ファイルのメタデータ フィールド。 たとえば、gdalwarp ソフトウェアは次のようになります。
マップをノースアップに変換するために使用されます (注: いくつかの gdalwarp パラメータがあり、
リサンプリング アルゴリズムを選択します):
gdalwarp 回転した.tif northup.tif
"エラー: 投影 of データセット ありません 現れる 〜へ match 現在 位置。"
インポートするデータと投影法が一致する場所を作成する必要があります。 試す
場所 の投影情報に基づいて新しい場所を作成するパラメータ
ファイル。 必要に応じて、別の場所に再投影できます r.プロジェクト.
または、次のコマンドを使用してこのエラーをオーバーライドできます。 -o フラグ。
"警告: G_set_window(): 違法 緯度 for 北"
GRASS の緯度/経度の場所は、北緯 90 度または北緯 XNUMX 度を超える地域を持つことはできません。
南。 地理参照されていない画像には、画像の数に基づく座標が含まれます。
ピクセル: 左下に 0,0。 cols,rows は右上にあります。 通常、画像は多くの
高さが 90 ピクセルを超えるため、GIS はインポートを拒否します。 確実な場合は、
データは緯度/経度の場所に適しており、マップの境界をリセットすることを目的としています
r.地域 モジュールをインポート直後に使用できます -l マップを制限するフラグ
有効な値に調整します。 結果の境界と解像度はおそらく
マップのデータが変更されず、安全になります。 既知にリセットした後
との境界 r.地域 あなたはそれらを再確認する必要があります r.info、特に注意してください
マップの解像度。 ほとんどの場合、データファイルのネイティブにインポートする必要があります
投影、または単純な XY 位置に投影し、地理補正ツール (i.修正する et
al.) ターゲットの場所に適切に投影します。 の -l フラグすべき の もしあなたが
投影が正しいことはわかっているが、内部のジオリファレンスが失われている。
マップの境界と解像度がどうあるべきかを事前に知っておいてください。
例
ECAD Rescale データ
European Climate Assessment and Dataset (ECAD) プロジェクトは、
1950 年から 2010 年までのヨーロッパ。
netCDF ファイルとしてプロジェクトを作成する場合、オフセット パラメータを使用して、毎日番号が付けられたラスターを取得できます。
1 年 1950 月 XNUMX 日以降のマップ。 LatLong の場所にいることを確認してください。
# 降水データのインポート
r.in.gdal -o 入力=rr_0.25deg_reg_1950-1964_v4.0.nc 出力=降水オフセット=0
r.in.gdal -o 入力=rr_0.25deg_reg_1965-1979_v4.0.nc 出力=降水オフセット=5479
r.in.gdal -o 入力=rr_0.25deg_reg_1980-1994_v4.0.nc 出力=降水オフセット=10957
r.in.gdal -o 入力=rr_0.25deg_reg_1995-2010_v4.0.nc 出力=降水オフセット=16436
# 空気圧データのインポート
r.in.gdal -o 入力=pp_0.25deg_reg_1950-1964_v4.0.nc 出力=air_pressure オフセット=0
r.in.gdal -o 入力=pp_0.25deg_reg_1965-1979_v4.0.nc 出力=air_pressure オフセット=5479
r.in.gdal -o 入力=pp_0.25deg_reg_1980-1994_v4.0.nc 出力=air_pressure オフセット=10957
r.in.gdal -o 入力=pp_0.25deg_reg_1995-2010_v4.0.nc 出力=air_pressure オフセット=16436
# 最低気温データをインポート
r.in.gdal -o 入力=tn_0.25deg_reg_1950-1964_v4.0.nc 出力=temperatur_min オフセット=0
r.in.gdal -o 入力=tn_0.25deg_reg_1965-1979_v4.0.nc 出力=temperatur_min オフセット=5479
r.in.gdal -o 入力=tn_0.25deg_reg_1980-1994_v4.0.nc 出力=temperatur_min オフセット=10957
r.in.gdal -o 入力=tn_0.25deg_reg_1995-2010_v4.0.nc 出力=temperatur_min オフセット=16436
# 最高気温データをインポート
r.in.gdal -o 入力=tx_0.25deg_reg_1950-1964_v4.0.nc 出力=temperatur_max オフセット=0
r.in.gdal -o 入力=tx_0.25deg_reg_1965-1979_v4.0.nc 出力=temperatur_max オフセット=5479
r.in.gdal -o 入力=tx_0.25deg_reg_1980-1994_v4.0.nc 出力=temperatur_max オフセット=10957
r.in.gdal -o 入力=tx_0.25deg_reg_1995-2010_v4.0.nc 出力=temperatur_max オフセット=16436
# 平均気温データをインポート
r.in.gdal -o 入力=tg_0.25deg_reg_1950-1964_v4.0.nc 出力=temperatur_mean オフセット=0
r.in.gdal -o 入力=tg_0.25deg_reg_1965-1979_v4.0.nc 出力=temperatur_mean オフセット=5479
r.in.gdal -o 入力=tg_0.25deg_reg_1980-1994_v4.0.nc 出力=temperatur_mean オフセット=10957
r.in.gdal -o 入力=tg_0.25deg_reg_1995-2010_v4.0.nc 出力=temperatur_mean オフセット=16436
GTOPO30 DEM
GTOPO30 データが正しく読み込まれないようにするには、新しい行「PIXELTYPE」を追加します。
SIGNEDINT" を .HDR に追加して、ファイルを無署名ではなく署名済みとして強制的に解釈します。
整数。 その後、.DEM ファイルをインポートできます。 最後に、たとえば「地形」カラー テーブルは次のことができます。
インポートされたマップに割り当てられます r.色.
GLOBE DEM
GLOBE DEM タイル (約 1km の解像度、GTOPO30 DEM データよりも優れている) をインポートするには、ユーザー
関連する HDR ファイルを追加でダウンロードする必要があります。 最後に、たとえば「地形」の色
テーブルはインポートされたマップに割り当てることができます r.色. 彼らの DEM ポータルも参照してください。
Worldclim.org
Worldclim データをインポートするには、次の行を各 .hdr ファイルに追加する必要があります。
ピクセルタイプ署名
HDF
HDF バンドのインポートには、個々のバンドの指定が必要です。
GDAL:
# MODIS FPAR の例
gdalinfo MOD15A2.A2003153.h18v04.004.2003171141042.hdf
...
サブデータセット:
SUBDATASET_1_NAME=HDF4_EOS:EOS_GRID:"MOD15A2.A2003153.h18v04.004.2003171141042.hdf":MOD_Grid_MOD15A2:Fpar_1km
SUBDATASET_1_DESC=[1200x1200] Fpar_1km MOD_Grid_MOD15A2 (8 ビット符号なし整数)
SUBDATASET_2_NAME=HDF4_EOS:EOS_GRID:"MOD15A2.A2003153.h18v04.004.2003171141042.hdf":MOD_Grid_MOD15A2:Lai_1km
SUBDATASET_2_DESC=[1200x1200] Lai_1km MOD_Grid_MOD15A2 (8 ビット符号なし整数)
...
# 最初のバンドのインポート、ここでは FPAR 1km:
r.in.gdal HDF4_EOS:EOS_GRID:"MOD15A2.A2003153.h18v04.004.2003171141042.hdf":MOD_Grid_MOD15A2:Fpar_1km \
out=fpar_1km_2003_06_02
# ... ファイル内の他の HDF バンドについても同様です。
onworks.net サービスを使用してオンラインで r.in.gdalgrass を使用する