Tôi đang cố gắng sử dụng QGIS để xem một số hình ảnh từ Máy ảnh bối cảnh tàu thám hiểm sao Hỏa. Tôi sử dụng ISIS3 tải những hình ảnh và chuyển đổi chúng sang file .cub ánh xạ, theo các hướng dẫn ở đây: http://isis.astrogeology.usgs.gov/IsisWorkshop/index.php/Working_with_Mars_Reconnaissance_Orbiter_CTX_Data
Bây giờ tôi đang cố gắng xuất các tệp khối ISIS này sang các tệp GTiff, để xem trong QGIS, tôi có thể làm điều này nhưng hội nghị địa lý dường như không hoạt động. Tôi đã thử cách này theo 3 cách khác nhau:
Sử dụng lệnh ISIS isis2std để tạo các tệp tif (và tfw).
Khi tôi cố gắng mở tệp tif trong QGIS, nó yêu cầu tôi chọn loại hình chiếu (câu hỏi phụ, có chiếu nào phù hợp cho Sao Hỏa không? Tôi luôn chọn WGS 84, tôi biết là ngu ngốc nhưng tôi không biết một tùy chọn hợp lý hơn) . Sau đó, khi raster được tải, nó sẽ tự động hiển thị dưới dạng hình ảnh thang độ xám nhưng tọa độ của tôi không tính bằng độ, chúng có vẻ như là pixel và có phạm vi 100 nghìn.
Phương pháp này tạo ra hình ảnh tif ~ 100 Mb từ .cub gốc gần 1 Gb.
Sử dụng gdal: gdal_translate -of GTiff input_mapped.cub output.tif.
Lần này tôi không được yêu cầu chọn loại hình chiếu. Phép chiếu diễn ra trên 'NGƯỜI DÙNG: 100002'. Hy vọng rằng điều này có nghĩa là nó đang quản lý để đọc phép chiếu sao Hỏa thích hợp từ khối ISIS? Hình ảnh hiển thị dưới dạng một hộp màu xám ở thang màu xám lần này, nhưng tôi có thể thấy chi tiết bằng cách chuyển sang sơ đồ màu khác. Một lần nữa, mặc dù tọa độ của tôi không tính theo độ, nhưng chúng lại nằm trong hàng trăm nghìn lần nữa.
Phương pháp này tạo ra một GTiff có cùng kích thước ~ 1 Gb như hình ảnh .cub ban đầu.
Mở tệp ISIS .cub trực tiếp trong QGIS.
Về cơ bản có tác dụng tương tự như đối với phương pháp 2.
Do kích thước tệp, tôi thích một biến thể của phương pháp 1, nhưng với sơ đồ chiếu và hội thảo địa lý phù hợp cho Sao Hỏa.
Chỉnh sửa: Đây là địa lý và đầu ra của gdalinfo trên tệp ISIS .cub ban đầu:
địa lý (trong 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 (từ thiết bị đầu cuối):
$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