i wanna plot rainfall on Nile river basin
from netCDF4 import Dataset
import numpy as np
from mpl_toolkits.basemap import Basemap
import warnings
import climdex.precipitation as pdex
import matplotlib as mpl
import matplotlib.pyplot as plt
from matplotlib import cm
import geopandas
warnings.filterwarnings('ignore')
nc_file=r"E:\climate change\code\clipped03.nc"
data=Dataset(nc_file)
lat=data.variables['lat'][:]
print("Extracting Latitdues")
lon=data.variables['lon'][:]
print("Extracting Longitude")
time=data.variables['time'][:]
print("Extracting time")
rainfall=data.variables['precip'][:]
print("Extracting rainfall")
map = Basemap(projection='merc',llcrnrlon=22.,llcrnrlat=-4.,urcrnrlon=40.,urcrnrlat=32.,resolution='i')
lons,lats=np.meshgrid(lon,lat)
x,y=map(lons,lats)
c=map.pcolor(x,y,np.squeeze(rainfall[0,:,:]),cmap='jet')
map.drawcountries
cbar=map.colorbar(c,location='right',pad='10%')
plt.show()
- Badges
- Science topic
- Similar topics
- Social Psychology
- Attitude