python - Python の地図上に Cassini-Soldner 投影で RGB 衛星画像をプロットする

Question

私は MODIS 装置から geotiff 形式で地理参照された RGB 衛星画像を持っています。cartopy を使用してマップ上にプロットし、RGB カラーを保持する正しい方法は何ですか?

私が推測する主な障害は、Cassini-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 コードからの投影を定義しようとしました。そこで、Cassini-Soldner EPSG コード (9806) をググってみましたが、cartopyccrs.epsg()はそれを認識しません。

メソッドを使いたいのですが、軸を作成する際の射影キーワードとして何を使用し、 で引数plt.imshow()として何を渡すか、少し混乱しています。transform=imshow