color对照 python,python3 matplotlib多个子图对应同一个colorbar

# -*- coding: utf-8 -*-

"""

Created on Tue Nov 12 19:31:07 2019

@author: BAI Depei

"""

from osgeo import gdal

import numpy as np

import sympy as sp

import matplotlib.pyplot as plt

from matplotlib.colors import Normalize

#-------------------------读取2007、2008、2009年TRMM降水遥感图像-----------------------

path = 'E:\\课程PPT\\Python空间数据处理\\python空间数据处理-期中大作业02\\卫星降雨数据TRMM3B43_分辨率0.25度\\'

#有时候路径必须双反斜杠

dataset = gdal.Open(path + '3B43_2007.TIF')

dataset2 = gdal.Open(path + '3B43_2008.TIF')

dataset3 = gdal.Open(path + '3B43_2009.TIF')

samples = dataset.RasterXSize

lines = dataset.RasterYSize

bands = dataset.RasterCount

img_geotrans = dataset.GetGeoTransform()

#Out[10]: (-180.0, 0.25, 0.0, 50.0, 0.0, -0.25)

img_proj = dataset.GetProjection()

#Out[12]: 'GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84"...]]]

#可以看到此遥感图像未进行投影定义，因此还不是投影坐标，只是地理坐标

#在求解这些实测站点位于图像中哪个行列位置时不用将其转为投影坐标

im_data1 = dataset.ReadAsArray(0,0,samples,lines)

im_data2 = dataset2.ReadAsArray(0,0,samples,lines)

im_data3 = dataset3.ReadAsArray(0,0,samples,lines)

del dataset

del dataset2</