ISIS a GDAL para QGIS: Marte de georreferenciación de imágenes

Question

Estoy tratando de usar QGIS para ver algunas imágenes de la Mars Reconnaissance Orbiter de la Cámara de Contexto. Yo uso ISIS3 para descargar las imágenes y convertirlas a asignar .cachorro de archivos, siguiendo las instrucciones aquí: http://isis.astrogeology.usgs.gov/IsisWorkshop/index.php/Working_with_Mars_Reconnaissance_Orbiter_CTX_Data

Ahora estoy tratando de exportar estos ISIS archivos de cubo para GTiff archivos, para ver en QGIS, me puede hacer esto, pero la georeferenciación no parece funcionar. He intentado esto en 3 formas diferentes:

1. Utilice el ISIS comando isis2std para generar tif (y tfw) archivos.

   Cuando intento abrir el archivo tif en QGIS me pide que elija un tipo de proyección (de lado la pregunta, ¿existe una adecuada proyección de Marte? Yo siempre elegir WGS 84, estúpido lo sé, pero no sé más sensible opción). A continuación, una vez que la trama se carga automáticamente se muestra como una imagen en escala de grises, pero mi coordenadas no son en grados, que aparecen en píxeles y tiene un rango de 100s de miles de personas.

   Este método produce una imagen tif de ~100 Mb de un original .cachorros de cerca de 1 Gb.

2. El uso de gdal: gdal_translate de GTiff input_mapped.cachorros de salida.tif.

   Esta vez yo no le pida que elija un tipo de proyección. La proyección toma en "USUARIO:100002'. Esperemos que esto significa que se está gestionando para ver la correspondiente Marte proyección de la ISIS cub? La imagen muestra como un sólido cuadro gris en escala de grises de este tiempo, pero puedo ver los detalles de cambiar a un color diferente esquema de asignación. De nuevo a pesar de mi coordenadas no son en grados, rango en el que cientos de miles de nuevo.

   Este método produce un GTiff con el mismo ~1 Gb de tamaño que el original .cub imagen.

3. Abra el ISIS .cub archivo directamente en QGIS.

   Básicamente tiene el mismo efecto que para el método 2.

Debido a los tamaños de archivo prefiero una variación del método 1, pero con una adecuada proyección y georeferenciación de un esquema de Marte.

Edit: Aquí está la georeferenciación, y la salida de gdalinfo en el original de ISIS .cub archivo:

georeferencia (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 (desde el 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

Answer 1

Hacer un custom CRS mediante el parámetro siguiente cadena:

+proj=longlat +a=3396190 +b=3396190 +no_defs

Puede llamar a esta CRS (decir) "Mars_Geo". A continuación, carga el *.cub trama, establece el proyecto de CRS para "Mars_geo" y activar la opción on-the-fly CRS transformación. Haciendo esto distorsionará la trama, en cierta medida, pero el lon/lat coordenadas reportadas como mover el puntero por la pantalla.

Cuando la comprobación de la esquina de las coordenadas de la trama tener en cuenta que las esquinas no son las esquinas de los datos. La imagen es vertical rectángulo con las dimensiones que se muestra en la imagen adjunta.

Por favor, tenga en cuenta que yo no sé nada acerca de Marte y nada acerca de los Sinusoidal Sir, así que ten cuidado!

N.