Изображение спутника RGB в проекции Кассини-Сольдера на карте в Python

Question

У меня есть спутниковое изображение RGB с привязкой к географии с прибора MODIS в формате геотипов. Каков правильный способ построения карты на карте с использованием картографии и сохранения цветов RGB?

Основным препятствием я думаю, это проекция изображения, которое является Кассини-Soldner:

import numpy as np 
from osgeo import gdal, osr 

ds = gdal.Open('modis_201303261252_rgb.tif') 

print(ds.GetGeoTransform()) 
(-1669791.8857914428, 250.0, 0.0, 1669792.327327792, 0.0, -250.0) 

proj = ds.GetProjection() 

inproj = osr.SpatialReference() 
inproj.ImportFromWkt(proj) 
print(inproj) 

PROJCS["unnamed", 
    GEOGCS["unnamed ellipse", 
     DATUM["unknown", 
      SPHEROID["unnamed",6378137,0]], 
     PRIMEM["Greenwich",0], 
     UNIT["degree",0.0174532925199433]], 
    PROJECTION["Cassini_Soldner"], 
    PARAMETER["latitude_of_origin",72], 
    PARAMETER["central_meridian",-4], 
    PARAMETER["false_easting",0], 
    PARAMETER["false_northing",0], 
    UNIT["metre",1, 
     AUTHORITY["EPSG","9001"]]] 

Я пытался следовать этому примеру https://ocefpaf.github.io/python4oceanographers/blog/2015/03/02/geotiff/ и использовать cartopy определить проекцию из кода EPSG. Таким образом, я googled для кода Cassini-Soldner EPSG (9806), но cartopy's ccrs.epsg() не узнает его.

Я хочу использовать plt.imshow() метод, но я немного запутался, что использовать в качестве проекционного ключевого слова при создании оси и что передать в качестве аргумента в transform=imshow.

Answer 1

(я бы добавить это как комментарий, но мне нужно 50 респ.)

Из интереса, как вы пришли на проекции Кассини MODIS RGB?

Альтернативой является использование запроса на сборку Cartopy, который я сделал, что позволяет использовать измерение времени с помощью GeoAxes.add_wmts(), чтобы вы могли вытаскивать фрагменты изображений MODIS RGB для вашего времени из NASA GIBS прямо в Cartopy GeoAxes: https://github.com/SciTools/cartopy/pull/788

Edit: эта функция находится в V0.15