|
|
import numpy as np
|
|
|
import matplotlib.pyplot as pl
|
|
|
import wradlib
|
|
|
import warnings
|
|
|
#export WRADLIB_DATA=/path/to/wradlib-data
|
|
|
warnings.filterwarnings('ignore')
|
|
|
'''
|
|
|
try:
|
|
|
get_ipython().magic('matplotlib inline')
|
|
|
except:
|
|
|
pl.ion()
|
|
|
'''
|
|
|
filename = wradlib.util.get_wradlib_data_file("/home/soporte/Downloads/raa00-dx_10908-0806021735-fbg---bin.gz")
|
|
|
img, meta = wradlib.io.read_dx(filename)
|
|
|
print("Shape of polar array: %r\n" % (img.shape,))
|
|
|
print("Some meta data of the DX file:")
|
|
|
print("\tdatetime: %r" % (meta["datetime"],))
|
|
|
print("\tRadar ID: %s" % (meta["radarid"],))
|
|
|
|
|
|
img[200:250,:]= np.ones([50,img.shape[1]])*np.nan
|
|
|
|
|
|
img[300:360,:]= np.ones([60,img.shape[1]])*np.nan
|
|
|
|
|
|
cgax, pm= wradlib.vis.plot_ppi(img)
|
|
|
txt = pl.title('Simple PPI')
|
|
|
print("coordenada angular",img[:,0],len(img[:,0]))
|
|
|
print("COORDENADA 0",img[0],len(img[0]))
|
|
|
cbar = pl.gcf().colorbar(pm, pad=0.075)
|
|
|
|
|
|
#r = np.arange(40, 80)
|
|
|
#az = np.arange(200, 250)
|
|
|
#ax, pm = wradlib.vis.plot_ppi(img[200:250, 40:80], r, az, autoext=False)
|
|
|
#ax, pm = wradlib.vis.plot_ppi(img[200:250, 40:80], r, az)
|
|
|
|
|
|
#txt = pl.title('Sector PPI')
|
|
|
pl.show()
|
|
|
|
|
|
|
|
|
### PARA MULTIPLE PloteO REVISAR EL LINK
|
|
|
#https://docs.wradlib.org/en/stable/notebooks/visualisation/wradlib_plot_curvelinear_grids.html
|
|
|
'''
|
|
|
subplots = [221, 222, 223, 224]
|
|
|
fig = pl.figure(figsize=(10,8))
|
|
|
fig.subplots_adjust(wspace=0.2, hspace=0.35)
|
|
|
for sp in subplots:
|
|
|
cgax, pm = wrl.vis.plot_rhi(ma1, r, th, rf=1e3, ax=sp, proj='cg')
|
|
|
caax = cgax.parasites[0]
|
|
|
paax = cgax.parasites[1]
|
|
|
t = pl.title('CG RHI #%(sp)d' %locals(), y=1.1)
|
|
|
cgax.set_ylim(0, 15)
|
|
|
cbar = pl.gcf().colorbar(pm, pad=0.125, ax=paax)
|
|
|
caax.set_xlabel('range [km]')
|
|
|
caax.set_ylabel('height [km]')
|
|
|
gh = cgax.get_grid_helper()
|
|
|
# set theta to some nice values
|
|
|
locs = [0., 5., 10., 15., 20., 30., 40., 60., 90.]
|
|
|
gh.grid_finder.grid_locator1 = FixedLocator(locs)
|
|
|
gh.grid_finder.tick_formatter1 = DictFormatter(dict([(i, r"${0:.0f}^\circ$".format(i)) for i in locs]))
|
|
|
cbar.set_label('reflectivity [dBZ]')
|
|
|
'''
|
|
|
|
