Py-Art progrma

122 views
Skip to first unread message

Yogesh Kolte

unread,
Oct 26, 2016, 6:12:04 AM10/26/16
to wradli...@googlegroups.com
Hello Sir,
Please find the attached file of radar data which i plotted  in pyart.


import numpy as np
import matplotlib.pyplot as plt
import matplotlib as mpl
import pyart
from netCDF4 import Dataset

filenamex = '/media/user/ExtHD2/20140731/MDV-20140731-155752-PPIVol.nc'
radarx = pyart.io.read(filenamex)
nc_fidx = Dataset(filenamex, 'r')
latsx = nc_fidx.variables['latitude'][:] 
lonsx = nc_fidx.variables['longitude'][:]
display = pyart.graph.RadarMapDisplay(radarx)
font = {'size': 16}
mpl.rc('font', **font)

fig1 = plt.figure(figsize=[30,30])
 ax1 = fig1.add_subplot(111 )
  
   
display.plot_ppi_map('DBZH',  vmin=-30, vmax=60,
                         min_lon=72.5, max_lon=75.2, min_lat=16.9, max_lat=19.2,
                         lon_lines=np.arange(72.5, 75.2, .5), projection='lcc',
                         lat_lines=np.arange(16.9, 19.2, .5), resolution='h',
                         lat_0=radarx.latitude['data'][0],
                         lon_0=radarx.longitude['data'][0])
   
    # plot range rings at 10, 20, 30 and 40km
display.plot_range_ring(30., line_style='r-')
display.plot_range_ring(50., line_style='r--')
display.plot_range_ring(100., line_style='r-')
display.plot_range_ring(125., line_style='r--')
   
# plots cross hairs
 display.plot_line_xy(np.array([-150000.0, 150000.0]), np.array([0.0, 0.0]),
                         line_style='k-')
display.plot_line_xy(np.array([0.0, 0.0]), np.array([-150000.0, 1500000.0]),
                         line_style='k-')
   ​
    Thank You

Yogesh Kolte

unread,
Dec 7, 2016, 1:34:45 PM12/7/16
to wradli...@googlegroups.com
Hello every one,

please help me to read attached file in wradlib.

https://drive.google.com/file/d/0B6FsIZaR20E8ZzhVTjZDcUhfYm8/view?usp=drive_web

Thank you,
yogesh kolte


Kai Muehlbauer

unread,
Dec 8, 2016, 3:20:45 AM12/8/16
to wradli...@googlegroups.com
Hi Yogesh,

attached is a simple proof of principle script to read the volume data
and display the first scan.

HTH,
Kai
>> MDV-20140731-155752-PPIVol.nc
>> <https://drive.google.com/file/d/0B6FsIZaR20E8ZzhVTjZDcUhfYm8/view?usp=drive_web>
>> ​
>> Thank You
>>
>

--
Kai Muehlbauer
Meteorological Institute University of Bonn
Auf dem Huegel 20 | +49 228 739083
D-53121 Bonn | kai.mue...@uni-bonn.de
netcdf_vol.py
ncvol.png
Reply all
Reply to author
Forward
0 new messages