ISISからGDAL、QGIS：火星画像の地理参照
QGISを使用して火星偵察オービターコンテキストカメラからの画像を表示しようとしています。ISIS3を使用して画像をダウンロードし、マップされた.cubファイルに変換します。手順は次のとおりです。http：//isis.astrogeology.usgs.gov/IsisWorkshop/index.php/Working_with_Mars_Reconnaissance_Orbiter_CTX_Data 現在、これらのISISキューブファイルをGTiffファイルにエクスポートしようとしています。QGISで表示するには、これを実行できますが、ジオリファレンスが機能していないようです。私はこれを3つの異なる方法で試しました： ISISコマンドisis2stdを使用して、tif（およびtfw）ファイルを生成します。 QGISでtifファイルを開こうとすると、投影タイプを選択するように求められます（質問、火星に適切な投影はありますか？常にWGS 84を選択しますが、私は愚かですが、もっと賢明なオプションを知りません）。 。次に、ラスターが読み込まれると、自動的にグレースケールイメージとして表示されますが、座標は度単位ではなく、ピクセル単位で表示され、範囲は数十万です。 この方法では、1 Gbに近い元の.cubから〜100 Mbのtifイメージが生成されます。 gdalを使用します：gdal_translate -of GTiff input_mapped.cub output.tif。 今回は、投影タイプを選択するように求められません。投影は「USER：100002」をとります。うまくいけば、これは、ISISカブから適切な火星投影を読み取ることができていることを意味しますか？今回は画像がグレースケールで塗りつぶされた灰色のボックスとして表示されますが、別のカラーマッピングスキームに切り替えると詳細を確認できます。繰り返しになりますが、座標は度単位ではありませんが、数十万の範囲です。 このメソッドは、元の.cubイメージと同じ〜1 GbのサイズのGTiffを生成します。 QGISで直接ISIS .cubファイルを開きます。 基本的に方法2と同じ効果があります。 ファイルサイズのため、方法1のバリエーションをお勧めしますが、火星に適した投影および地理参照スキームを使用します。 編集：ジオリファレンスと、元のISIS .cubファイルでのgdalinfoの出力を次に示します。 ジオリファレンス（Python）： >>>import gdal >>>from gdalconst import * >>>fn = 'P01_001356_1747_XN_05S221W.map.cub' >>>ds = gdal.Open(fn, GA_ReadOnly) >>>gt = ds.GetGeoTransform() >>>gt (-22085.510544416, 5.1698292472885, 0.0, -234679.22885141, 0.0, -5.1698292472885) gdalinfo（ターミナルから）： $gdalinfo file.cub …