Ich versuche, mit QGIS einige Bilder von der Mars Reconnaissance Orbiter Context Camera anzuzeigen. Ich verwende ISIS3, um die Bilder herunterzuladen und in zugeordnete .cub-Dateien zu konvertieren. Befolgen Sie dazu die Anweisungen hier: http://isis.astrogeology.usgs.gov/IsisWorkshop/index.php/Working_with_Mars_Reconnaissance_Orbiter_CTX_Data
Jetzt versuche ich, diese ISIS-Cube-Dateien in GTiff-Dateien zu exportieren, um sie in QGIS anzuzeigen. Ich kann dies tun, aber die Georeferenzierung scheint nicht zu funktionieren. Ich habe dies auf drei verschiedene Arten versucht:
Verwenden Sie den ISIS-Befehl isis2std, um tif- (und tfw-) Dateien zu generieren.
Wenn ich versuche, die TIF-Datei in QGIS zu öffnen, werde ich aufgefordert, einen Projektionstyp auszuwählen (Nebenfrage, gibt es eine geeignete Projektion für den Mars? Ich wähle immer WGS 84, dumm, ich weiß, aber ich kenne keine vernünftigere Option). . Sobald das Raster geladen ist, wird es automatisch als Graustufenbild angezeigt, aber meine Koordinaten sind nicht in Grad angegeben. Sie scheinen in Pixel zu sein und haben einen Bereich von Hunderttausenden.
Diese Methode erzeugt ein TIF-Bild von ~ 100 MB aus einem ursprünglichen .cub von nahezu 1 GB.
Verwenden Sie gdal: gdal_translate -of GTiff input_mapped.cub output.tif.
Dieses Mal werde ich nicht gebeten, einen Projektionstyp auszuwählen. Die Projektion nimmt 'USER: 100002' an. Hoffentlich bedeutet dies, dass es gelingt, die entsprechende Mars-Projektion aus dem ISIS-Jungen zu lesen? Das Bild wird diesmal als graues Feld in Graustufen angezeigt, aber ich kann Details sehen, indem ich zu einem anderen Farbzuordnungsschema wechsle. Obwohl meine Koordinaten nicht in Grad angegeben sind, liegen sie wieder im Bereich von Hunderttausenden.
Diese Methode erzeugt einen GTiff mit der gleichen Größe von ~ 1 GB wie das ursprüngliche .cub-Bild.
Öffnen Sie die ISIS .cub-Datei direkt in QGIS.
Hat grundsätzlich den gleichen Effekt wie bei Methode 2.
Aufgrund der Dateigröße würde ich eine Variation von Methode 1 bevorzugen, jedoch mit einem geeigneten Projektions- und Georeferenzierungsschema für den Mars.
Bearbeiten: Hier ist die Georeferenz und die Ausgabe von gdalinfo in der ursprünglichen ISIS .cub-Datei:
Georeferenz (in 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 (vom Terminal):
$gdalinfo file.cub
Driver: ISIS3/USGS Astrogeology ISIS cube (Version 3)
Files: P01_001356_1747_XN_05S221W.map.cub
P01_001356_1747_XN_05S221W.map.cub.aux.xml
Size is 8526, 30111
Coordinate System is:
PROJCS["Sinusoidal Mars",
GEOGCS["GCS_Mars",
DATUM["D_Mars",
SPHEROID["Mars",3396190,0]],
PRIMEM["Reference_Meridian",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Sinusoidal"],
PARAMETER["longitude_of_center",138.45730511261],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0]]
Origin = (-22085.510544415999902,-234679.228851410007337)
Pixel Size = (5.169829247288500,-5.169829247288500)
Corner Coordinates:
Upper Left ( -22085.511, -234679.229) (138d 5' 1.74"E, 3d57'33.05"S)
Lower Left ( -22085.511, -390347.957) (138d 4'56.04"E, 6d35' 7.46"S)
Upper Right ( 21992.454, -234679.229) (138d49'45.19"E, 3d57'33.05"S)
Lower Right ( 21992.454, -390347.957) (138d49'50.86"E, 6d35' 7.46"S)
Center ( -46.528, -312513.593) (138d27'23.46"E, 5d16'20.26"S)
Band 1 Block=406x128 Type=Float32, ColorInterp=Undefined
Min=0.033 Max=0.150
Minimum=0.033, Maximum=0.150, Mean=0.080, StdDev=0.020
NoData Value=-3.40282265508890445e+38
Metadata:
STATISTICS_MAXIMUM=0.15037615597248
STATISTICS_MEAN=0.080181184414784
STATISTICS_MINIMUM=0.033307410776615
STATISTICS_STDDEV=0.020299974127104
