import matplotlib as mp
import numpy as np

ind='hoge.grd'
Dat = np.netcdf(ind,'r').variables['z'][::-1]
lon=np.netcdf(ind,'r').variables['x'][::-1]
lat=np.netcdf(ind,'r').variables['y'][::-1]
NX, NY=np.meshgrid(lon, lat)

#Die Mitte der Farbleiste unter Angabe der Minimal- und Maximalwerte der Farbleiste(z=0)Die Farbe von bleibt erhalten.
norm = mp.colors.TwoSlopeNorm(vcenter=0.0, vmin=-0.2, vmax=0.1)
#Schneiden Sie einen Teil der 2D-Daten aus[pixel]
extX0=950
extX1=1600
extY0=1100
extY1=1750

#Geben Sie auch die Figurengröße an
fig=plt.figure(figsize=(10,10))
#Jet Colormap
plt.pcolormesh(NX[extY0:extY1, extX0:extX1], NY[extY0:extY1, extX0:extX1], Dat[extY0:extY1, extX0:extX1], cmap='jet', norm=norm)　
plt.colorbar() #Farbbalkenanzeige
plt.xlabel('Longitude [deg.]')
plt.ylabel('Latitude [deg.]')
plt.show()

Recommended Posts