J'essaie d'utiliser QGIS pour afficher certaines images de la caméra de contexte Mars Reconnaissance Orbiter. J'utilise ISIS3 pour télécharger les images et les convertir en fichiers .cub mappés, en suivant les instructions ici: http://isis.astrogeology.usgs.gov/IsisWorkshop/index.php/Working_with_Mars_Reconnaissance_Orbiter_CTX_Data
Maintenant, j'essaie d'exporter ces fichiers de cube ISIS vers des fichiers GTiff, pour les afficher dans QGIS, je peux le faire mais le géoréférencement ne semble pas fonctionner. J'ai essayé cela de 3 manières différentes:
Utilisez la commande ISIS isis2std pour générer des fichiers tif (et tfw).
Lorsque j'essaie d'ouvrir le fichier tif dans QGIS, il me demande de choisir un type de projection (question secondaire, y a-t-il une projection appropriée pour Mars? Je choisis toujours WGS 84, stupide je sais mais je ne connais pas d'option plus sensée) . Ensuite, une fois le raster chargé, il s'affiche automatiquement sous la forme d'une image en niveaux de gris, mais mes coordonnées ne sont pas en degrés, elles semblent être en pixels et avoir une plage de 100s de milliers.
Cette méthode produit une image tif d'environ 100 Mo à partir d'un fichier original de près de 1 Go.
Utilisez gdal: gdal_translate -of GTiff input_mapped.cub output.tif.
Cette fois, on ne me demande pas de choisir un type de projection. La projection prend «USER: 100002». Espérons que cela signifie qu'il parvient à lire la projection appropriée de Mars à partir du petit ISIS? Cette fois, l'image s'affiche sous la forme d'un cadre gris uni dans une échelle de gris, mais je peux voir les détails en passant à un schéma de mappage de couleurs différent. Encore une fois, bien que mes coordonnées ne soient pas en degrés, elles se situent à nouveau par centaines de milliers.
Cette méthode produit un GTiff avec la même taille ~ 1 Go que l'image .cub d'origine.
Ouvrez le fichier ISIS .cub directement dans QGIS.
A essentiellement le même effet que pour la méthode 2.
En raison de la taille des fichiers, je préférerais une variante de la méthode 1, mais avec un schéma de projection et de géoréférencement approprié pour Mars.
Edit: Voici la géoréférence et la sortie de gdalinfo sur le fichier ISIS .cub d'origine:
géoréférence (en 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 (depuis le 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