파이썬의지도에서 Cassini-Soldner 투영의 RGB 위성 이미지를

Question

MODIS 기기의 geo-referenced RGB 위성 이미지가 geotiff 형식입니다. Cartopy를 사용하여 RGB 색상을 보존 한지도에 올바른 색상을 표시하는 방법은 무엇입니까? I는 예를 https://ocefpaf.github.io/python4oceanographers/blog/2015/03/02/geotiff/ 를 수행하고 EPSG 코드로부터 돌출부를 정의 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"]]] 

:

I 추측 주요 장애물 카시-Soldner되는 이미지의 투영이다. 그래서 나는 Cassini-Soldner EPSG 코드 (9806)를 찾아 봤지만 cartopy의 ccrs.epsg()은 그것을 인식하지 못했다.

plt.imshow() 메서드를 사용하고 싶지만 축이 만들어 질 때 transform= 인수로 전달할 때 투영 키워드로 사용할 항목을 imshow에 약간 혼동합니다.

Answer 1

(내가 주석으로이를 추가하지만 50 렙이 필요합니다.)

을 밖으로 관심, 당신은 어떻게 카시니 돌기 MODIS RGB에서 왔는가?

다른 방법은 GeoAxes.add_wmts()와 함께 사용할 시간 차원을 허용하는 Cartopy pull 요청을 사용하는 것이므로 NASA GIBS에서 MODIS RGB 이미지 타일을 가져올 수 있습니다. Cartopy GeoAxes로 : https://github.com/SciTools/cartopy/pull/788

편집 :이 기능은 내가 역, 던디 대학, 영국 (http://www.sat.dundee.ac 받기 NEODAAS NERC 위성에서있어