ISIS to GDAL to QGIS: Georeferences Mars Images


8

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:

  1. 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.

  2. 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.

  3. 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

Chuỗi tham số Proj cho 'USER: 100002' là gì?
nhopton

@nhopton Sẽ là thế này: "+ proj = eqc + lat_ts = 0 + lat_0 = 0 + lon_0 = 180 + x_0 = 0 + y_0 = 0 + a = 3396190 + b = 3396190 + units = m + no_defs" Cài đặt> CRS tùy chỉnh và đây là trong hộp tham số.
EddyTheB

Tôi tin rằng điều này là do geotransform được tính bằng pixel, chứ không phải độ, đang điều tra.
EddyTheB

Tôi sẽ không nghĩ như vậy, không phải nếu tọa độ nằm trong hàng trăm ngàn. Có vẻ như nhiều khả năng raster được dự kiến, một người hy vọng vào CRS mà bạn đã cung cấp chi tiết ở trên. Bạn có thể gửi raster đến một nơi mà nó có thể được tải xuống không?
nhopton

@nhopton Chắc chắn, mặc dù nó rất lớn. Tôi cũng đã chỉnh sửa câu hỏi của mình để hiển thị đầu ra của gdalinfo và địa lý cho tệp. Xin lỗi vì chậm phản hồi ý kiến ​​của bạn, tôi đã bị cuốn hút vào những thứ khác. Tệp ISIS .cub có tại therm.boulder.swri.edu/images/iêu
EddyTheB 15/03/13

Câu trả lời:


2

Tạo CRS tùy chỉnh bằng chuỗi tham số sau:

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

Bạn có thể gọi CRS này (nói) "Mars_Geo". Sau đó tải raster * .cub, đặt CRS dự án thành "Mars_geo" và bật chuyển đổi CRS đang hoạt động. Làm điều này sẽ làm biến dạng raster đến một mức độ nào đó nhưng bạn sẽ nhận được tọa độ lon / lat khi bạn di chuyển con trỏ xung quanh màn hình.

Khi kiểm tra tọa độ góc của raster, hãy nhớ rằng các góc không phải là góc của dữ liệu. Hình ảnh là một hình chữ nhật thẳng đứng có các phạm vi được hiển thị trong hình ảnh đính kèm.

Xin lưu ý rằng tôi không biết gì về Sao Hỏa và bên cạnh không có gì về CRS hình sin, vì vậy hãy cẩn thận!

N

nhập mô tả hình ảnh ở đây

Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.