read_netCDF.py
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)
#Le centre de la barre de couleur tout en spécifiant les valeurs minimum et maximum de la barre de couleur(z=0)La couleur de est conservée.
norm = mp.colors.TwoSlopeNorm(vcenter=0.0, vmin=-0.2, vmax=0.1)
#Découpez une partie des données 2D[pixel]
extX0=950
extX1=1600
extY0=1100
extY1=1750
#Spécifiez également la taille de la figure
fig=plt.figure(figsize=(10,10))
#couleur de jet
plt.pcolormesh(NX[extY0:extY1, extX0:extX1], NY[extY0:extY1, extX0:extX1], Dat[extY0:extY1, extX0:extX1], cmap='jet', norm=norm)
plt.colorbar() #Affichage de la barre de couleur
plt.xlabel('Longitude [deg.]')
plt.ylabel('Latitude [deg.]')
plt.show()
Un chiffre comme celui-ci s'affiche ↓
Recommended Posts