Pythonでマップ上のCassini-Soldner投影のプロットRGB衛星画像

Question

私は地理座標系のgeos-reference RGB衛星画像をgeotiff形式で持っています。 カートリッジを使用してマップにプロットし、RGBカラーを保存する正しい方法は何ですか？私はこの例にhttps://ocefpaf.github.io/python4oceanographers/blog/2015/03/02/geotiff/ に従い、EP​​SGコードから投影を定義するcartopyを使用しよう

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"]]] 

：

らしい主な障害は、カッシーニ-Soldnerされた画像の投影です。だから私はカッシーナ - ソルナーのEPSGコード（9806）のためにgoogled、しかしcartopyのccrs.epsg()はそれを認識しません。

私はplt.imshow()メソッドを使いたいですが、imshowのtransform=引数として渡すべきものは、軸が作成されたときに投影キーワードとして何を使うのか混乱します。

Answer 1

（私はコメントとしてこれを追加することになりますが、私は50担当者を必要としています。）

の関心のうち、どのようにカッシーニ投影MODIS RGBで来ましたの？

代わりに、GeoAxes.add_wmts（）で時間ディメンションを使用できるように、私が作成したCartopyプルリクエストを使用して、MODIS RGB画像タイルをNASA GIBSからまっすぐに引き出すことができますCartopy GeoAxesへ：https://github.com/SciTools/cartopy/pull/788

編集：この機能は、私は駅、ダンディー大学、英国（http://www.sat.dundee.acを受信NEODAAS NERC衛星からそれを持っている