Attenuation correction module
Attenuation due to precipitation plays a major role in underestimating the precipitation. Therefore, an attenuation correction according to the ZPHI method (Testud et. al, 2001) is performed.
Now we start step by step with the attenuation correction.
At first, we have to import following packages.
[1]:
import datetime as dt
import wradlib as wrl
import warnings
warnings.filterwarnings('ignore')
Note
It is recommended to perform a clutter correction first to remove all non-meteorological eochos from the data set.
The functions work also for a xarray dataset which is linking over the time.
Phase processing
First of all, a phase processing has to be done with the phase_zphi function. Therefore, a binary array (True = rain, False = not-rain) rolling a range-long sum from PHIDP is generated. Afterwards the first occurrence of the maximum (from front and back) of PHIDP binary array is searched for. This gives the indices of the centers of the rolling window. Then min, max, median and mean of these are calculated. At least the true start and end indices of the phase PHIDP are calculated (+- half
window length).
[2]:
cphase = wrf.attenuation_corr.phase_zphi(ds_clutter_corr.PHIDP,
rng=1000.0,
start_range=0.0)
display(cphase)
<xarray.Dataset>
Dimensions: (azimuth: 720, range: 936)
Coordinates:
* azimuth (azimuth) float64 0.25 0.75 1.25 1.75 ... 358.8 359.2 359.8
elevation (azimuth) float64 0.5 0.5 0.5 0.5 0.5 ... 0.5 0.5 0.5 0.5 0.5
rtime (azimuth) datetime64[ns] 2022-02-16T13:45:28.941817600 ... ...
* range (range) float32 37.5 112.5 187.5 ... 7.009e+04 7.016e+04
time datetime64[ns] 2022-02-16T13:45:01
sweep_mode <U20 'azimuth_surveillance'
longitude float64 13.24
latitude float64 53.55
altitude float64 38.0
Data variables: (12/15)
phib (azimuth, range) float64 nan nan nan nan ... nan nan nan nan
offset float64 450.0
offset_idx int64 6
start_range (azimuth) float32 7.988e+03 7.838e+03 ... 7.988e+03 7.988e+03
stop_range (azimuth) float32 2.651e+04 2.621e+04 ... 2.689e+04 2.689e+04
first_min (azimuth) float32 32.34 32.34 36.56 ... 36.56 30.94 33.75
... ...
first_idx (azimuth) int64 106 104 106 105 105 ... 124 123 123 106 106
last_min (azimuth) float32 32.34 30.94 29.53 ... 36.56 35.16 33.75
last_max (azimuth) float32 47.81 43.59 43.59 45.0 ... 46.41 45.0 49.22
last_mean (azimuth) float32 37.87 35.46 35.76 ... 40.08 39.38 38.97
last_median (azimuth) float32 35.86 34.45 36.56 ... 39.38 39.38 38.67
last_idx (azimuth) int64 352 348 346 325 325 ... 366 357 357 357 357xarray.Dataset
- azimuth: 720
- range: 936
- azimuth(azimuth)float640.25 0.75 1.25 ... 359.2 359.8
- standard_name :
- ray_azimuth_angle
- long_name :
- azimuth_angle_from_true_north
- units :
- degrees
- axis :
- radial_azimuth_coordinate
- a1gate :
- 1
- angle_res :
- 0.5
array([2.5000e-01, 7.5000e-01, 1.2500e+00, ..., 3.5875e+02, 3.5925e+02, 3.5975e+02]) - elevation(azimuth)float640.5 0.5 0.5 0.5 ... 0.5 0.5 0.5 0.5
- standard_name :
- ray_elevation_angle
- long_name :
- elevation_angle_from_horizontal_plane
- units :
- degrees
- axis :
- radial_elevation_coordinate
array([0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, ... 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5]) - rtime(azimuth)datetime64[ns]2022-02-16T13:45:28.941817600 .....
- standard_name :
- time
array(['2022-02-16T13:45:28.941817600', '2022-02-16T13:45:28.980652544', '2022-02-16T13:45:01.019417600', '2022-02-16T13:45:01.058252544', '2022-02-16T13:45:01.097087744', '2022-02-16T13:45:01.135922688', '2022-02-16T13:45:01.174757632', '2022-02-16T13:45:01.213592832', '2022-02-16T13:45:01.252427776', '2022-02-16T13:45:01.291262976', '2022-02-16T13:45:01.330097920', '2022-02-16T13:45:01.368932864', '2022-02-16T13:45:01.407768064', '2022-02-16T13:45:01.446603008', '2022-02-16T13:45:01.485438208', '2022-02-16T13:45:01.524273152', '2022-02-16T13:45:01.563108096', '2022-02-16T13:45:01.601943296', '2022-02-16T13:45:01.640778240', '2022-02-16T13:45:01.679613440', '2022-02-16T13:45:01.718448384', '2022-02-16T13:45:01.757283328', '2022-02-16T13:45:01.796118528', '2022-02-16T13:45:01.834953472', '2022-02-16T13:45:01.873788672', '2022-02-16T13:45:01.912623616', '2022-02-16T13:45:01.951458816', '2022-02-16T13:45:01.990293760', '2022-02-16T13:45:02.029128704', '2022-02-16T13:45:02.067963904', '2022-02-16T13:45:02.106798848', '2022-02-16T13:45:02.145634048', '2022-02-16T13:45:02.184468992', '2022-02-16T13:45:02.223303936', '2022-02-16T13:45:02.262139136', '2022-02-16T13:45:02.300974080', '2022-02-16T13:45:02.339809280', '2022-02-16T13:45:02.378644224', '2022-02-16T13:45:02.417479168', '2022-02-16T13:45:02.456314368', ... '2022-02-16T13:45:27.466085632', '2022-02-16T13:45:27.504920832', '2022-02-16T13:45:27.543755776', '2022-02-16T13:45:27.582590720', '2022-02-16T13:45:27.621425920', '2022-02-16T13:45:27.660260864', '2022-02-16T13:45:27.699096064', '2022-02-16T13:45:27.737931008', '2022-02-16T13:45:27.776765952', '2022-02-16T13:45:27.815601152', '2022-02-16T13:45:27.854436096', '2022-02-16T13:45:27.893271296', '2022-02-16T13:45:27.932106240', '2022-02-16T13:45:27.970941184', '2022-02-16T13:45:28.009776384', '2022-02-16T13:45:28.048611328', '2022-02-16T13:45:28.087446528', '2022-02-16T13:45:28.126281472', '2022-02-16T13:45:28.165116672', '2022-02-16T13:45:28.203951616', '2022-02-16T13:45:28.242786560', '2022-02-16T13:45:28.281621760', '2022-02-16T13:45:28.320456704', '2022-02-16T13:45:28.359291904', '2022-02-16T13:45:28.398126848', '2022-02-16T13:45:28.436961792', '2022-02-16T13:45:28.475796992', '2022-02-16T13:45:28.514631936', '2022-02-16T13:45:28.553467136', '2022-02-16T13:45:28.592302080', '2022-02-16T13:45:28.631137024', '2022-02-16T13:45:28.669972224', '2022-02-16T13:45:28.708807168', '2022-02-16T13:45:28.747642368', '2022-02-16T13:45:28.786477312', '2022-02-16T13:45:28.825312256', '2022-02-16T13:45:28.864147456', '2022-02-16T13:45:28.902982400'], dtype='datetime64[ns]') - range(range)float3237.5 112.5 ... 7.009e+04 7.016e+04
- units :
- meters
- standard_name :
- projection_range_coordinate
- long_name :
- range_to_measurement_volume
- spacing_is_constant :
- true
- axis :
- radial_range_coordinate
- meters_to_center_of_first_gate :
- 37.5
- meters_between_gates :
- 75.0
array([3.75000e+01, 1.12500e+02, 1.87500e+02, ..., 7.00125e+04, 7.00875e+04, 7.01625e+04], dtype=float32) - time()datetime64[ns]2022-02-16T13:45:01
- standard_name :
- time
array('2022-02-16T13:45:01.000000000', dtype='datetime64[ns]') - sweep_mode()<U20'azimuth_surveillance'
array('azimuth_surveillance', dtype='<U20') - longitude()float6413.24
- long_name :
- longitude
- units :
- degrees_east
- standard_name :
- longitude
array(13.24397)
- latitude()float6453.55
- long_name :
- latitude
- units :
- degrees_north
- positive :
- up
- standard_name :
- latitude
array(53.55478)
- altitude()float6438.0
- long_name :
- altitude
- units :
- meters
- standard_name :
- altitude
array(38.)
- phib(azimuth, range)float64nan nan nan nan ... nan nan nan nan
array([[nan, nan, nan, ..., nan, nan, nan], [nan, nan, nan, ..., nan, nan, nan], [nan, nan, nan, ..., nan, nan, nan], ..., [nan, nan, nan, ..., nan, nan, nan], [nan, nan, nan, ..., nan, nan, nan], [nan, nan, nan, ..., nan, nan, nan]]) - offset()float64450.0
array(450.)
- offset_idx()int646
array(6)
- start_range(azimuth)float327.988e+03 7.838e+03 ... 7.988e+03
array([ 7987.5, 7837.5, 7987.5, 7912.5, 7912.5, 7912.5, 7762.5, 7837.5, 7912.5, 7912.5, 7912.5, 8437.5, 8512.5, 8512.5, 8587.5, 8512.5, 8662.5, 8587.5, 8737.5, 8512.5, 8512.5, 8212.5, 9187.5, 9412.5, 9562.5, 9637.5, 9937.5, 10012.5, 10087.5, 10087.5, 10162.5, 10162.5, 10237.5, 10312.5, 10762.5, 10762.5, 10762.5, 10687.5, 10687.5, 10612.5, 10537.5, 10537.5, 10537.5, 13462.5, 13462.5, 10462.5, 10237.5, 10912.5, 11287.5, 14887.5, 14737.5, 14737.5, 10012.5, 9562.5, 10162.5, 10162.5, 9712.5, 9712.5, 10837.5, 14437.5, 14362.5, 14362.5, 14287.5, 14212.5, 11962.5, 11887.5, 11812.5, 11812.5, 12187.5, 12187.5, 10387.5, 10462.5, 9937.5, 9412.5, 9037.5, 9712.5, 9787.5, 9937.5, 9937.5, 9937.5, 9937.5, 17437.5, 17362.5, 10087.5, 10087.5, 9862.5, 9862.5, 9862.5, 9937.5, 9487.5, 9412.5, 9262.5, 9262.5, 9337.5, 9337.5, 9337.5, 9187.5, 9187.5, 9187.5, 9187.5, 8512.5, 8437.5, 8887.5, 9037.5, 9112.5, 9112.5, 9187.5, 9187.5, 9562.5, 9637.5, 11287.5, 11287.5, 11287.5, 11437.5, 10987.5, 9862.5, 9787.5, 10462.5, 11437.5, 9337.5, 9037.5, 10237.5, 11362.5, 9487.5, 11362.5, 9412.5, 10237.5, 10237.5, 9562.5, 9412.5, 7912.5, 7837.5, 7837.5, 7837.5, 7912.5, 7912.5, 7912.5, 7912.5, 8212.5, 8212.5, ... 7912.5, 7762.5, 7762.5, 7762.5, 7687.5, 7762.5, 7762.5, 7687.5, 7762.5, 8662.5, 8737.5, 8737.5, 8737.5, 8737.5, 8737.5, 8737.5, 8737.5, 8812.5, 9187.5, 9187.5, 9262.5, 9262.5, 9337.5, 9412.5, 9487.5, 9487.5, 9487.5, 9487.5, 9562.5, 9562.5, 9487.5, 9487.5, 9337.5, 9337.5, 9262.5, 9487.5, 9487.5, 8887.5, 8812.5, 8887.5, 9112.5, 9337.5, 10162.5, 10162.5, 10162.5, 10162.5, 9337.5, 9262.5, 9112.5, 8887.5, 8137.5, 8062.5, 8062.5, 7687.5, 7687.5, 7687.5, 7687.5, 7762.5, 7612.5, 8512.5, 8662.5, 8587.5, 8812.5, 8812.5, 8662.5, 8737.5, 9487.5, 13237.5, 13237.5, 13237.5, 13237.5, 12787.5, 12712.5, 12487.5, 12187.5, 12262.5, 12262.5, 12187.5, 12112.5, 12112.5, 12112.5, 12112.5, 12112.5, 12112.5, 12037.5, 12037.5, 12037.5, 11962.5, 11887.5, 11962.5, 12037.5, 12112.5, 12112.5, 12262.5, 12262.5, 12637.5, 14137.5, 14137.5, 14137.5, 14062.5, 13987.5, 14062.5, 13987.5, 14062.5, 14137.5, 12337.5, 12337.5, 12262.5, 11587.5, 11512.5, 11512.5, 11062.5, 10987.5, 9712.5, 9637.5, 9712.5, 9712.5, 9637.5, 9637.5, 9637.5, 9712.5, 10762.5, 10762.5, 10762.5, 9412.5, 9337.5, 9337.5, 9337.5, 9262.5, 9262.5, 7987.5, 7987.5], dtype=float32) - stop_range(azimuth)float322.651e+04 2.621e+04 ... 2.689e+04
array([26512.5, 26212.5, 26062.5, 24487.5, 24487.5, 24412.5, 69262.5, 24112.5, 23812.5, 23737.5, 23737.5, 23812.5, 30787.5, 31012.5, 31087.5, 31312.5, 30637.5, 30562.5, 30487.5, 30487.5, 30487.5, 30412.5, 32362.5, 32737.5, 28762.5, 32062.5, 31987.5, 31987.5, 31987.5, 31762.5, 41887.5, 41962.5, 50662.5, 50962.5, 50962.5, 50962.5, 50962.5, 50962.5, 51262.5, 51037.5, 51262.5, 52012.5, 52312.5, 52387.5, 53287.5, 53212.5, 53362.5, 53362.5, 53362.5, 52387.5, 51937.5, 51262.5, 49612.5, 49237.5, 49312.5, 48037.5, 58237.5, 58237.5, 58162.5, 57787.5, 57712.5, 57637.5, 33562.5, 33187.5, 33187.5, 33187.5, 33112.5, 32662.5, 31987.5, 31462.5, 31312.5, 31237.5, 64912.5, 65812.5, 65512.5, 65512.5, 65587.5, 64462.5, 64462.5, 63937.5, 63937.5, 63712.5, 63412.5, 63337.5, 63412.5, 63337.5, 62212.5, 58612.5, 60487.5, 60487.5, 60112.5, 60037.5, 57262.5, 55612.5, 55312.5, 54862.5, 53137.5, 52762.5, 51787.5, 51712.5, 51637.5, 49612.5, 49612.5, 49162.5, 49162.5, 48937.5, 48937.5, 49537.5, 48862.5, 48862.5, 48712.5, 48637.5, 48637.5, 48262.5, 48037.5, 47962.5, 47287.5, 46387.5, 46312.5, 46012.5, 40312.5, 39487.5, 38962.5, 38887.5, 38812.5, 38512.5, 32962.5, 32212.5, 32062.5, 32062.5, 32062.5, 31537.5, 31462.5, 31087.5, 30787.5, 30712.5, 30712.5, 30112.5, 30112.5, 30037.5, ... 15637.5, 15637.5, 16762.5, 16237.5, 16762.5, 16687.5, 16687.5, 12037.5, 12337.5, 14137.5, 14287.5, 12787.5, 13387.5, 16087.5, 15562.5, 15637.5, 13312.5, 13312.5, 15787.5, 16837.5, 13162.5, 29137.5, 15787.5, 13987.5, 15412.5, 15637.5, 15637.5, 12787.5, 14287.5, 26137.5, 28912.5, 30562.5, 30562.5, 30637.5, 30637.5, 30712.5, 30562.5, 30562.5, 30562.5, 30562.5, 30562.5, 30487.5, 29887.5, 29887.5, 29812.5, 29587.5, 15187.5, 18787.5, 18787.5, 19087.5, 19012.5, 19087.5, 19162.5, 19087.5, 10087.5, 10012.5, 8662.5, 8887.5, 8887.5, 9637.5, 10162.5, 9787.5, 10012.5, 14137.5, 16087.5, 16687.5, 17812.5, 18412.5, 18562.5, 20587.5, 20512.5, 20587.5, 19837.5, 20362.5, 21187.5, 21187.5, 21487.5, 21487.5, 21487.5, 21937.5, 21937.5, 21562.5, 21562.5, 21562.5, 21112.5, 20662.5, 20362.5, 19237.5, 19012.5, 19012.5, 18787.5, 67462.5, 67537.5, 67537.5, 66937.5, 68212.5, 67987.5, 68287.5, 68062.5, 67762.5, 36562.5, 36562.5, 36862.5, 37612.5, 41137.5, 42337.5, 44287.5, 44287.5, 43837.5, 43087.5, 43012.5, 42487.5, 41662.5, 41287.5, 41212.5, 40987.5, 58762.5, 58687.5, 39712.5, 58162.5, 36262.5, 36037.5, 35512.5, 35437.5, 35287.5, 31987.5, 30637.5, 27562.5, 26887.5, 26887.5, 26887.5, 26887.5], dtype=float32) - first_min(azimuth)float3232.34 32.34 36.56 ... 30.94 33.75
array([ 32.34375, 32.34375, 36.5625 , 36.5625 , 36.5625 , 33.75 , 32.34375, 33.75 , 33.75 , 33.75 , 33.75 , 33.75 , 30.9375 , 29.53125, 30.9375 , 32.34375, 32.34375, 33.75 , 35.15625, 35.15625, 36.5625 , 29.53125, 35.15625, 36.5625 , 36.5625 , 35.15625, 32.34375, 33.75 , 32.34375, 30.9375 , 30.9375 , 29.53125, 28.125 , 30.9375 , 32.34375, 30.9375 , 32.34375, 30.9375 , 30.9375 , 29.53125, 26.71875, 22.5 , 32.34375, 32.34375, 32.34375, 36.5625 , 39.375 , 40.78125, 33.75 , 32.34375, 37.96875, 33.75 , 37.96875, 37.96875, 29.53125, 33.75 , 32.34375, 26.71875, 32.34375, 30.9375 , 32.34375, 32.34375, 33.75 , 30.9375 , 28.125 , 29.53125, 28.125 , 25.3125 , 29.53125, 26.71875, 29.53125, 29.53125, 25.3125 , 26.71875, 29.53125, 33.75 , 32.34375, 25.3125 , 22.5 , 28.125 , 32.34375, 30.9375 , 32.34375, 26.71875, 26.71875, 30.9375 , 29.53125, 28.125 , 30.9375 , 30.9375 , 35.15625, 30.9375 , 32.34375, 33.75 , 32.34375, 28.125 , 30.9375 , 29.53125, 28.125 , 28.125 , 29.53125, 26.71875, 28.125 , 30.9375 , 28.125 , 28.125 , 30.9375 , 30.9375 , 28.125 , 19.6875 , 25.3125 , 21.09375, 23.90625, 25.3125 , 28.125 , 28.125 , 28.125 , 28.125 , 26.71875, 16.875 , ... 23.90625, 28.125 , 26.71875, 26.71875, 25.3125 , 25.3125 , 28.125 , 26.71875, 28.125 , 25.3125 , 22.5 , 22.5 , 23.90625, 25.3125 , 23.90625, 25.3125 , 25.3125 , 23.90625, 25.3125 , 25.3125 , 23.90625, 26.71875, 28.125 , 28.125 , 28.125 , 28.125 , 28.125 , 30.9375 , 23.90625, 22.5 , 26.71875, 28.125 , 25.3125 , 30.9375 , 32.34375, 23.90625, 29.53125, 30.9375 , 32.34375, 35.15625, 40.78125, 42.1875 , 42.1875 , 37.96875, 33.75 , 33.75 , 33.75 , 42.1875 , 33.75 , 35.15625, 30.9375 , 32.34375, 30.9375 , 32.34375, 32.34375, 32.34375, 32.34375, 33.75 , 33.75 , 33.75 , 32.34375, 33.75 , 33.75 , 32.34375, 32.34375, 32.34375, 32.34375, 33.75 , 33.75 , 32.34375, 30.9375 , 30.9375 , 33.75 , 33.75 , 33.75 , 33.75 , 33.75 , 29.53125, 37.96875, 36.5625 , 35.15625, 35.15625, 32.34375, 29.53125, 28.125 , 28.125 , 30.9375 , 32.34375, 32.34375, 33.75 , 32.34375, 32.34375, 36.5625 , 36.5625 , 33.75 , 36.5625 , 30.9375 , 30.9375 , 30.9375 , 29.53125, 32.34375, 33.75 , 35.15625, 35.15625, 35.15625, 33.75 , 33.75 , 35.15625, 35.15625, 35.15625, 37.96875, 36.5625 , 30.9375 , 33.75 ], dtype=float32) - first_max(azimuth)float3243.59 40.78 42.19 ... 43.59 43.59
array([ 43.59375, 40.78125, 42.1875 , 43.59375, 43.59375, 43.59375, 43.59375, 40.78125, 40.78125, 39.375 , 40.78125, 40.78125, 40.78125, 40.78125, 40.78125, 39.375 , 42.1875 , 43.59375, 43.59375, 42.1875 , 43.59375, 45. , 39.375 , 42.1875 , 42.1875 , 40.78125, 39.375 , 40.78125, 40.78125, 40.78125, 39.375 , 39.375 , 39.375 , 39.375 , 40.78125, 40.78125, 42.1875 , 39.375 , 37.96875, 35.15625, 39.375 , 37.96875, 40.78125, 40.78125, 40.78125, 47.8125 , 49.21875, 49.21875, 46.40625, 39.375 , 40.78125, 40.78125, 53.4375 , 47.8125 , 46.40625, 45. , 49.21875, 49.21875, 46.40625, 42.1875 , 40.78125, 39.375 , 37.96875, 40.78125, 43.59375, 43.59375, 42.1875 , 43.59375, 45. , 42.1875 , 40.78125, 39.375 , 40.78125, 46.40625, 43.59375, 42.1875 , 43.59375, 40.78125, 45. , 43.59375, 46.40625, 37.96875, 42.1875 , 46.40625, 43.59375, 42.1875 , 37.96875, 42.1875 , 43.59375, 47.8125 , 45. , 43.59375, 40.78125, 40.78125, 39.375 , 42.1875 , 45. , 46.40625, 40.78125, 40.78125, 40.78125, 39.375 , 39.375 , 40.78125, 40.78125, 40.78125, 37.96875, 39.375 , 39.375 , 42.1875 , 37.96875, 39.375 , 39.375 , 43.59375, 46.40625, 45. , 42.1875 , 46.40625, 37.96875, 37.96875, ... 35.15625, 35.15625, 35.15625, 36.5625 , 37.96875, 37.96875, 36.5625 , 37.96875, 39.375 , 37.96875, 37.96875, 36.5625 , 37.96875, 35.15625, 35.15625, 36.5625 , 33.75 , 30.9375 , 32.34375, 36.5625 , 39.375 , 37.96875, 39.375 , 37.96875, 39.375 , 39.375 , 39.375 , 39.375 , 39.375 , 36.5625 , 35.15625, 35.15625, 37.96875, 36.5625 , 39.375 , 39.375 , 40.78125, 42.1875 , 45. , 46.40625, 46.40625, 49.21875, 49.21875, 50.625 , 46.40625, 45. , 43.59375, 47.8125 , 52.03125, 45. , 39.375 , 40.78125, 39.375 , 39.375 , 39.375 , 37.96875, 39.375 , 39.375 , 39.375 , 42.1875 , 43.59375, 43.59375, 40.78125, 39.375 , 36.5625 , 39.375 , 36.5625 , 37.96875, 39.375 , 40.78125, 39.375 , 40.78125, 39.375 , 39.375 , 39.375 , 42.1875 , 39.375 , 42.1875 , 43.59375, 43.59375, 43.59375, 40.78125, 39.375 , 37.96875, 39.375 , 39.375 , 42.1875 , 40.78125, 42.1875 , 40.78125, 40.78125, 43.59375, 45. , 46.40625, 42.1875 , 42.1875 , 40.78125, 45. , 40.78125, 39.375 , 42.1875 , 42.1875 , 43.59375, 39.375 , 40.78125, 39.375 , 42.1875 , 45. , 43.59375, 45. , 45. , 43.59375, 43.59375, 43.59375], dtype=float32) - first_mean(azimuth)float3238.57 37.97 39.27 ... 37.37 38.87
array([ 38.57143 , 37.96875 , 39.274555, 39.274555, 39.776787, 38.973213, 37.566963, 37.06473 , 36.763393, 37.06473 , 38.069195, 36.964287, 36.160713, 36.160713, 36.86384 , 36.5625 , 39.07366 , 39.375 , 39.475445, 39.67634 , 40.78125 , 39.174107, 37.16518 , 37.96875 , 38.27009 , 37.767857, 36.662945, 36.5625 , 36.26116 , 35.75893 , 35.859375, 35.055805, 35.256695, 36.160713, 36.86384 , 36.5625 , 36.763393, 35.65848 , 34.453125, 33.317307, 33.44866 , 32.745537, 36.763393, 36.462055, 36.26116 , 42.68973 , 42.890625, 45.30134 , 40.982143, 37.16518 , 39.174107, 39.07366 , 44.799107, 44.19643 , 40.078125, 40.982143, 42.589287, 40.78125 , 39.274555, 38.069195, 37.265625, 36.06027 , 35.558037, 35.45759 , 35.256695, 36.129807, 34.151787, 35.264423, 36.160713, 35.357143, 35.45759 , 35.15625 , 36.26116 , 36.5625 , 37.36607 , 37.06473 , 37.06473 , 35.65848 , 35.15625 , 35.357143, 39.375 , 35.15625 , 37.66741 , 36.160713, 33.44866 , 36.462055, 33.317307, 34.754463, 36.86384 , 38.77232 , 39.274555, 38.47098 , 37.16518 , 37.66741 , 35.859375, 35.055805, 36.462055, 35.256695, 35.256695, 35.45759 , ... 29.732143, 29.83259 , 33.14732 , 33.75 , 32.243305, 33.75 , 34.65402 , 34.55357 , 34.85491 , 32.243305, 31.037947, 30.435268, 32.34375 , 33.348213, 33.549107, 36.26116 , 33.14732 , 35.357143, 36.86384 , 39.07366 , 41.420456, 44.4375 , 44.0625 , 45. , 42.728367, 40.379463, 38.973213, 39.475445, 45.100445, 41.58482 , 37.767857, 35.45759 , 36.361607, 35.055805, 35.558037, 35.75893 , 34.25223 , 35.558037, 35.65848 , 36.06027 , 37.265625, 36.26116 , 36.06027 , 35.65848 , 34.85491 , 35.055805, 34.65402 , 34.453125, 35.65848 , 36.160713, 37.06473 , 36.361607, 36.86384 , 36.06027 , 36.06027 , 36.06027 , 36.763393, 36.06027 , 36.462055, 40.17857 , 38.77232 , 38.370537, 36.964287, 35.45759 , 34.65402 , 34.05134 , 34.151787, 35.256695, 35.15625 , 35.558037, 36.5625 , 38.47098 , 38.87277 , 40.78125 , 41.58482 , 38.671875, 38.61779 , 36.462055, 34.85491 , 35.15625 , 34.85491 , 36.662945, 37.06473 , 38.47098 , 37.767857, 37.96875 , 37.46652 , 37.868305, 39.475445, 39.174107, 39.475445, 40.580357, 39.67634 , 37.36607 , 38.87277 ], dtype=float32) - first_median(azimuth)float3237.97 37.97 39.38 ... 37.97 38.67
array([ 37.96875 , 37.96875 , 39.375 , 39.375 , 39.375 , 39.375 , 37.265625, 36.5625 , 36.5625 , 37.265625, 39.375 , 36.5625 , 36.5625 , 35.859375, 37.96875 , 36.5625 , 39.375 , 40.78125 , 39.375 , 40.78125 , 40.78125 , 40.78125 , 37.265625, 37.265625, 37.96875 , 37.96875 , 36.5625 , 36.5625 , 36.5625 , 35.859375, 36.5625 , 35.859375, 35.15625 , 36.5625 , 37.265625, 36.5625 , 36.5625 , 35.15625 , 34.453125, 33.75 , 34.453125, 33.75 , 36.5625 , 36.5625 , 36.5625 , 43.59375 , 42.1875 , 45.703125, 41.484375, 37.96875 , 39.375 , 39.375 , 43.59375 , 45. , 40.78125 , 42.1875 , 42.1875 , 43.59375 , 40.078125, 39.375 , 37.96875 , 35.859375, 35.15625 , 35.15625 , 35.15625 , 35.15625 , 35.15625 , 35.15625 , 37.265625, 35.15625 , 35.859375, 35.15625 , 36.5625 , 37.265625, 37.265625, 37.265625, 37.265625, 37.265625, 35.859375, 33.75 , 38.671875, 35.15625 , 37.96875 , 35.859375, 33.046875, 36.5625 , 33.75 , 34.453125, 36.5625 , 39.375 , 39.375 , 38.671875, 37.96875 , 37.96875 , 35.15625 , 35.15625 , 35.859375, 33.75 , 36.5625 , 37.96875 , ... 28.125 , 29.53125 , 33.75 , 33.75 , 32.34375 , 33.75 , 36.5625 , 34.453125, 35.15625 , 31.640625, 32.34375 , 30.9375 , 32.34375 , 33.75 , 33.75 , 36.5625 , 35.15625 , 35.859375, 37.265625, 39.375 , 40.78125 , 45. , 42.890625, 44.296875, 43.59375 , 40.78125 , 37.96875 , 39.375 , 45.703125, 41.484375, 37.96875 , 35.15625 , 36.5625 , 34.453125, 35.15625 , 35.859375, 33.75 , 35.859375, 35.15625 , 36.5625 , 36.5625 , 35.859375, 35.15625 , 35.15625 , 34.453125, 35.15625 , 33.75 , 33.75 , 35.15625 , 36.5625 , 37.265625, 36.5625 , 37.265625, 35.859375, 36.5625 , 36.5625 , 36.5625 , 36.5625 , 37.265625, 39.375 , 37.96875 , 37.96875 , 36.5625 , 35.15625 , 33.75 , 33.75 , 34.453125, 35.15625 , 34.453125, 35.15625 , 35.859375, 39.375 , 39.375 , 40.78125 , 40.78125 , 39.375 , 37.96875 , 35.859375, 35.15625 , 34.453125, 34.453125, 36.5625 , 36.5625 , 37.96875 , 37.96875 , 38.671875, 37.96875 , 37.96875 , 38.671875, 38.671875, 38.671875, 39.375 , 39.375 , 37.96875 , 38.671875], dtype=float32) - first_idx(azimuth)int64106 104 106 105 ... 123 123 106 106
array([106, 104, 106, 105, 105, 105, 103, 104, 105, 105, 105, 112, 113, 113, 114, 113, 115, 114, 116, 113, 113, 109, 122, 125, 127, 128, 132, 133, 134, 134, 135, 135, 136, 137, 143, 143, 143, 142, 142, 141, 140, 140, 140, 179, 179, 139, 136, 145, 150, 198, 196, 196, 133, 127, 135, 135, 129, 129, 144, 192, 191, 191, 190, 189, 159, 158, 157, 157, 162, 162, 138, 139, 132, 125, 120, 129, 130, 132, 132, 132, 132, 232, 231, 134, 134, 131, 131, 131, 132, 126, 125, 123, 123, 124, 124, 124, 122, 122, 122, 122, 113, 112, 118, 120, 121, 121, 122, 122, 127, 128, 150, 150, 150, 152, 146, 131, 130, 139, 152, 124, 120, 136, 151, 126, 151, 125, 136, 136, 127, 125, 105, 104, 104, 104, 105, 105, 105, 105, 109, 109, 109, 109, 104, 105, 106, 107, 334, 112, 114, 113, 113, 113, 113, 115, 109, 115, 115, 99, 99, 115, 115, 106, 106, 106, 102, 102, 114, 11, 98, 98, 103, 115, 115, 105, 109, 107, 106, 106, 105, 104, 104, 104, 105, 105, 106, 106, 104, 105, 108, 105, 111, 111, 105, 105, 104, 104, 102, 102, 102, 108, 107, 106, 116, 106, 106, 105, 105, 104, 103, 103, 103, 109, 109, 112, 114, 114, 114, 111, 108, 107, 107, 106, 106, 108, 108, 108, 108, 106, 106, 104, 104, 104, 103, 104, 104, 104, 105, 105, 104, 104, 104, 102, 102, 102, 103, 103, 106, 105, 105, 103, 103, 103, 106, 106, 109, 109, 109, 106, 106, 106, ... 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 217, 216, 217, 215, 216, 215, 206, 200, 210, 210, 208, 206, 208, 211, 212, 207, 200, 199, 200, 194, 194, 200, 178, 176, 175, 171, 185, 156, 191, 190, 186, 176, 177, 158, 188, 188, 188, 152, 187, 188, 153, 181, 181, 181, 165, 165, 165, 109, 168, 167, 165, 159, 106, 169, 170, 169, 150, 146, 127, 135, 125, 120, 126, 125, 114, 114, 121, 102, 102, 102, 102, 105, 103, 104, 104, 105, 105, 109, 106, 109, 105, 105, 105, 105, 136, 107, 141, 141, 108, 133, 125, 139, 139, 139, 103, 103, 106, 106, 106, 105, 106, 106, 106, 106, 106, 105, 105, 103, 103, 103, 102, 103, 103, 102, 103, 115, 116, 116, 116, 116, 116, 116, 116, 117, 122, 122, 123, 123, 124, 125, 126, 126, 126, 126, 127, 127, 126, 126, 124, 124, 123, 126, 126, 118, 117, 118, 121, 124, 135, 135, 135, 135, 124, 123, 121, 118, 108, 107, 107, 102, 102, 102, 102, 103, 101, 113, 115, 114, 117, 117, 115, 116, 126, 176, 176, 176, 176, 170, 169, 166, 162, 163, 163, 162, 161, 161, 161, 161, 161, 161, 160, 160, 160, 159, 158, 159, 160, 161, 161, 163, 163, 168, 188, 188, 188, 187, 186, 187, 186, 187, 188, 164, 164, 163, 154, 153, 153, 147, 146, 129, 128, 129, 129, 128, 128, 128, 129, 143, 143, 143, 125, 124, 124, 124, 123, 123, 106, 106]) - last_min(azimuth)float3232.34 30.94 29.53 ... 35.16 33.75
array([ 32.34375, 30.9375 , 29.53125, 35.15625, 35.15625, 33.75 , 29.53125, 33.75 , 32.34375, 33.75 , 33.75 , 35.15625, 35.15625, 36.5625 , 37.96875, 35.15625, 36.5625 , 36.5625 , 37.96875, 36.5625 , 36.5625 , 35.15625, 32.34375, 30.9375 , 33.75 , 33.75 , 32.34375, 32.34375, 30.9375 , 32.34375, 33.75 , 33.75 , 26.71875, 29.53125, 33.75 , 35.15625, 37.96875, 33.75 , 33.75 , 35.15625, 39.375 , 40.78125, 36.5625 , 40.78125, 37.96875, 42.1875 , 40.78125, 42.1875 , 40.78125, 42.1875 , 39.375 , 43.59375, 39.375 , 37.96875, 39.375 , 42.1875 , 42.1875 , 40.78125, 39.375 , 36.5625 , 33.75 , 29.53125, 35.15625, 37.96875, 37.96875, 36.5625 , 36.5625 , 30.9375 , 33.75 , 29.53125, 32.34375, 36.5625 , 30.9375 , 35.15625, 35.15625, 35.15625, 35.15625, 36.5625 , 42.1875 , 39.375 , 39.375 , 37.96875, 37.96875, 36.5625 , 36.5625 , 39.375 , 35.15625, 36.5625 , 43.59375, 45. , 43.59375, 43.59375, 33.75 , 37.96875, 37.96875, 35.15625, 39.375 , 40.78125, 35.15625, 36.5625 , 36.5625 , 35.15625, 33.75 , 30.9375 , 32.34375, 33.75 , 33.75 , 36.5625 , 32.34375, 29.53125, 26.71875, 30.9375 , 30.9375 , 30.9375 , 33.75 , 33.75 , 32.34375, 32.34375, 28.125 , 21.09375, ... 22.5 , 16.875 , 28.125 , 19.6875 , 23.90625, 23.90625, 23.90625, 23.90625, 28.125 , 26.71875, 21.09375, 21.09375, 22.5 , 25.3125 , 26.71875, 29.53125, 28.125 , 25.3125 , 26.71875, 26.71875, 30.9375 , 30.9375 , 29.53125, 26.71875, 28.125 , 21.09375, 26.71875, 28.125 , 25.3125 , 18.28125, 23.90625, 21.09375, 22.5 , 23.90625, 28.125 , 30.9375 , 26.71875, 30.9375 , 32.34375, 35.15625, 43.59375, 42.1875 , 42.1875 , 37.96875, 33.75 , 29.53125, 30.9375 , 29.53125, 30.9375 , 28.125 , 32.34375, 25.3125 , 26.71875, 33.75 , 30.9375 , 35.15625, 25.3125 , 29.53125, 29.53125, 29.53125, 30.9375 , 28.125 , 28.125 , 29.53125, 29.53125, 28.125 , 30.9375 , 33.75 , 33.75 , 39.375 , 37.96875, 35.15625, 32.34375, 42.1875 , 45. , 42.1875 , 42.1875 , 43.59375, 40.78125, 40.78125, 40.78125, 39.375 , 33.75 , 33.75 , 33.75 , 29.53125, 32.34375, 39.375 , 29.53125, 32.34375, 29.53125, 35.15625, 29.53125, 35.15625, 32.34375, 32.34375, 29.53125, 35.15625, 33.75 , 30.9375 , 28.125 , 29.53125, 29.53125, 26.71875, 25.3125 , 30.9375 , 26.71875, 32.34375, 30.9375 , 35.15625, 36.5625 , 36.5625 , 35.15625, 33.75 ], dtype=float32) - last_max(azimuth)float3247.81 43.59 43.59 ... 45.0 49.22
array([ 47.8125 , 43.59375, 43.59375, 45. , 46.40625, 43.59375, 43.59375, 43.59375, 47.8125 , 54.84375, 50.625 , 52.03125, 61.875 , 54.84375, 47.8125 , 49.21875, 46.40625, 43.59375, 47.8125 , 42.1875 , 42.1875 , 40.78125, 46.40625, 40.78125, 46.40625, 45. , 45. , 42.1875 , 40.78125, 45. , 46.40625, 43.59375, 40.78125, 43.59375, 45. , 45. , 42.1875 , 42.1875 , 42.1875 , 46.40625, 50.625 , 56.25 , 57.65625, 52.03125, 52.03125, 52.03125, 46.40625, 46.40625, 49.21875, 49.21875, 52.03125, 56.25 , 52.03125, 47.8125 , 49.21875, 49.21875, 56.25 , 50.625 , 49.21875, 47.8125 , 46.40625, 45. , 42.1875 , 43.59375, 47.8125 , 50.625 , 54.84375, 45. , 42.1875 , 43.59375, 39.375 , 40.78125, 50.625 , 47.8125 , 46.40625, 45. , 46.40625, 49.21875, 53.4375 , 49.21875, 50.625 , 47.8125 , 46.40625, 47.8125 , 49.21875, 52.03125, 47.8125 , 49.21875, 54.84375, 56.25 , 54.84375, 56.25 , 50.625 , 57.65625, 59.0625 , 52.03125, 50.625 , 50.625 , 47.8125 , 46.40625, 46.40625, 45. , 43.59375, 39.375 , 40.78125, 43.59375, 43.59375, 47.8125 , 45. , 46.40625, 43.59375, 40.78125, 47.8125 , 42.1875 , 40.78125, 43.59375, 47.8125 , 42.1875 , 42.1875 , 45. , ... 36.5625 , 40.78125, 36.5625 , 33.75 , 35.15625, 37.96875, 35.15625, 42.1875 , 40.78125, 40.78125, 40.78125, 43.59375, 46.40625, 46.40625, 42.1875 , 37.96875, 36.5625 , 39.375 , 40.78125, 42.1875 , 42.1875 , 39.375 , 40.78125, 43.59375, 43.59375, 45. , 39.375 , 39.375 , 46.40625, 42.1875 , 37.96875, 39.375 , 39.375 , 42.1875 , 43.59375, 46.40625, 39.375 , 43.59375, 45. , 50.625 , 49.21875, 52.03125, 49.21875, 50.625 , 46.40625, 46.40625, 39.375 , 42.1875 , 45. , 42.1875 , 39.375 , 39.375 , 40.78125, 42.1875 , 45. , 49.21875, 45. , 45. , 42.1875 , 42.1875 , 40.78125, 43.59375, 47.8125 , 42.1875 , 40.78125, 42.1875 , 46.40625, 50.625 , 56.25 , 47.8125 , 49.21875, 53.4375 , 49.21875, 56.25 , 57.65625, 59.0625 , 54.84375, 59.0625 , 60.46875, 54.84375, 54.84375, 56.25 , 45. , 42.1875 , 43.59375, 47.8125 , 54.84375, 50.625 , 56.25 , 52.03125, 50.625 , 46.40625, 46.40625, 47.8125 , 43.59375, 46.40625, 42.1875 , 46.40625, 52.03125, 49.21875, 46.40625, 43.59375, 47.8125 , 45. , 42.1875 , 46.40625, 49.21875, 46.40625, 39.375 , 42.1875 , 40.78125, 46.40625, 45. , 49.21875], dtype=float32) - last_mean(azimuth)float3237.87 35.46 35.76 ... 39.38 38.97
array([ 37.868305, 35.45759 , 35.75893 , 40.982143, 40.078125, 38.370537, 37.566963, 39.97768 , 39.174107, 41.484375, 40.78125 , 40.47991 , 44.899555, 42.99107 , 43.09152 , 43.29241 , 40.881695, 40.580357, 41.28348 , 39.475445, 39.375 , 37.868305, 38.370537, 36.462055, 40.27902 , 39.97768 , 38.973213, 37.16518 , 37.36607 , 38.069195, 41.28348 , 38.973213, 34.182693, 38.77232 , 40.47991 , 40.379463, 39.97768 , 38.671875, 38.57143 , 40.580357, 42.890625, 48.91741 , 48.11384 , 47.31027 , 43.89509 , 44.899555, 43.392857, 44.799107, 44.296875, 46.004463, 44.09598 , 48.214287, 45.100445, 43.694195, 43.59375 , 46.70759 , 47.109375, 46.205357, 43.794643, 42.087055, 39.87723 , 38.671875, 39.375 , 40.78125 , 42.1875 , 43.59375 , 43.59375 , 36.763393, 36.662945, 36.964287, 37.265625, 38.77232 , 40.680805, 40.680805, 39.97768 , 40.379463, 41.58482 , 43.89509 , 46.808037, 44.49777 , 45.100445, 43.493305, 41.785713, 42.388393, 43.392857, 43.794643, 40.673077, 43.392857, 48.816963, 49.921875, 49.620537, 50.02232 , 44.598213, 46.70759 , 45.80357 , 41.986607, 46.506695, 45.200893, 42.287945, 42.1875 , ... 33.549107, 35.15625 , 35.357143, 36.26116 , 36.662945, 36.06027 , 36.26116 , 32.444195, 33.649555, 34.39904 , 33.549107, 33.14732 , 31.138393, 32.243305, 32.94643 , 35.45759 , 38.57143 , 35.45759 , 38.370537, 39.07366 , 42.826706, 46.125 , 45.200893, 45.9375 , 43.59375 , 41.38393 , 36.86384 , 36.361607, 38.069195, 38.671875, 36.160713, 36.06027 , 34.955357, 34.65402 , 37.566963, 37.46652 , 41.785713, 36.763393, 38.47098 , 36.160713, 36.662945, 36.964287, 34.35268 , 36.361607, 34.55357 , 34.453125, 35.055805, 38.973213, 41.08259 , 44.09598 , 43.392857, 45. , 44.799107, 41.484375, 49.118305, 50.32366 , 50.524555, 49.118305, 51.12723 , 49.759617, 48.013393, 46.10491 , 46.08173 , 38.87277 , 39.07366 , 38.169643, 38.973213, 43.493305, 44.69866 , 42.79018 , 40.881695, 41.88616 , 40.47991 , 38.973213, 40.78125 , 40.078125, 38.069195, 37.66741 , 39.274555, 43.29241 , 39.87723 , 37.536057, 36.662945, 39.87723 , 37.06473 , 35.055805, 37.566963, 36.995193, 39.87723 , 35.859375, 39.174107, 38.57143 , 40.078125, 39.375 , 38.973213], dtype=float32) - last_median(azimuth)float3235.86 34.45 36.56 ... 39.38 38.67
array([ 35.859375, 34.453125, 36.5625 , 41.484375, 40.078125, 38.671875, 37.96875 , 40.78125 , 38.671875, 40.078125, 40.078125, 38.671875, 42.1875 , 40.78125 , 43.59375 , 45. , 40.078125, 40.78125 , 40.78125 , 39.375 , 39.375 , 37.96875 , 37.265625, 37.265625, 40.078125, 40.078125, 39.375 , 36.5625 , 37.96875 , 37.96875 , 42.1875 , 38.671875, 35.15625 , 40.78125 , 42.1875 , 40.78125 , 39.375 , 39.375 , 38.671875, 40.78125 , 41.484375, 48.515625, 48.515625, 46.40625 , 43.59375 , 43.59375 , 43.59375 , 45. , 43.59375 , 46.40625 , 42.1875 , 47.8125 , 45.703125, 44.296875, 43.59375 , 46.40625 , 46.40625 , 46.40625 , 43.59375 , 42.1875 , 40.78125 , 38.671875, 39.375 , 40.78125 , 42.1875 , 43.59375 , 41.484375, 35.859375, 35.15625 , 36.5625 , 37.96875 , 37.96875 , 40.78125 , 40.078125, 39.375 , 40.78125 , 42.1875 , 43.59375 , 47.109375, 43.59375 , 45. , 43.59375 , 42.1875 , 42.1875 , 43.59375 , 42.1875 , 40.78125 , 45. , 49.21875 , 49.921875, 50.625 , 49.921875, 45. , 45. , 45. , 40.78125 , 47.109375, 45.703125, 42.890625, 42.890625, ... 33.046875, 33.75 , 35.15625 , 37.96875 , 37.265625, 36.5625 , 35.15625 , 33.046875, 33.75 , 32.34375 , 35.15625 , 33.75 , 32.34375 , 32.34375 , 33.75 , 35.859375, 38.671875, 36.5625 , 37.96875 , 39.375 , 40.78125 , 46.40625 , 43.59375 , 46.40625 , 43.59375 , 40.78125 , 37.265625, 36.5625 , 38.671875, 40.78125 , 35.859375, 35.859375, 35.15625 , 35.15625 , 39.375 , 36.5625 , 42.1875 , 35.859375, 39.375 , 36.5625 , 37.265625, 37.96875 , 33.046875, 35.15625 , 33.75 , 35.15625 , 34.453125, 40.078125, 40.78125 , 45. , 43.59375 , 46.40625 , 45. , 42.1875 , 49.21875 , 49.21875 , 50.625 , 49.921875, 50.625 , 49.21875 , 48.515625, 43.59375 , 45. , 39.375 , 39.375 , 39.375 , 39.375 , 42.890625, 45. , 42.890625, 40.78125 , 41.484375, 40.078125, 40.78125 , 40.78125 , 41.484375, 37.96875 , 39.375 , 38.671875, 43.59375 , 37.96875 , 37.96875 , 37.265625, 38.671875, 37.96875 , 35.859375, 37.265625, 35.15625 , 40.78125 , 36.5625 , 39.375 , 39.375 , 39.375 , 39.375 , 38.671875], dtype=float32) - last_idx(azimuth)int64352 348 346 325 ... 357 357 357 357
array([352, 348, 346, 325, 325, 324, 922, 320, 316, 315, 315, 316, 409, 412, 413, 416, 407, 406, 405, 405, 405, 404, 430, 435, 382, 426, 425, 425, 425, 422, 557, 558, 674, 678, 678, 678, 678, 678, 682, 679, 682, 692, 696, 697, 709, 708, 710, 710, 710, 697, 691, 682, 660, 655, 656, 639, 775, 775, 774, 769, 768, 767, 446, 441, 441, 441, 440, 434, 425, 418, 416, 415, 864, 876, 872, 872, 873, 858, 858, 851, 851, 848, 844, 843, 844, 843, 828, 780, 805, 805, 800, 799, 762, 740, 736, 730, 707, 702, 689, 688, 687, 660, 660, 654, 654, 651, 651, 659, 650, 650, 648, 647, 647, 642, 639, 638, 629, 617, 616, 612, 536, 525, 518, 517, 516, 512, 438, 428, 426, 426, 426, 419, 418, 413, 409, 408, 408, 400, 400, 399, 399, 398, 398, 397, 393, 392, 390, 390, 389, 128, 134, 130, 127, 133, 129, 129, 127, 120, 126, 126, 126, 127, 135, 123, 131, 125, 128, 118, 116, 116, 119, 128, 128, 128, 142, 142, 141, 142, 136, 137, 137, 134, 132, 132, 134, 135, 136, 137, 137, 137, 145, 157, 157, 157, 157, 157, 156, 156, 156, 160, 159, 233, 233, 237, 242, 242, 242, 250, 251, 251, 263, 268, 269, 272, 276, 277, 280, 280, 284, 287, 287, 288, 288, 289, 289, 288, 288, 206, 206, 209, 208, 208, 214, 213, 213, 224, 226, 234, 235, 235, 235, 247, 253, 253, 253, 253, 254, 254, 254, 187, 182, 182, 195, 198, 202, 198, 198, 206, 209, 209, ... 19, 19, 19, 19, 19, 19, 19, 19, 18, 18, 17, 17, 17, 17, 238, 234, 237, 227, 227, 227, 227, 227, 221, 222, 222, 221, 220, 223, 223, 230, 226, 211, 214, 206, 212, 212, 234, 189, 190, 193, 203, 202, 202, 202, 202, 188, 191, 200, 200, 201, 200, 202, 200, 201, 202, 193, 193, 196, 176, 199, 196, 196, 195, 195, 189, 219, 222, 209, 209, 209, 225, 225, 225, 311, 311, 313, 313, 319, 319, 319, 319, 319, 319, 242, 242, 242, 214, 212, 212, 212, 190, 209, 189, 188, 188, 195, 187, 195, 185, 188, 188, 188, 193, 190, 190, 182, 183, 163, 163, 180, 183, 187, 187, 188, 151, 160, 152, 163, 163, 170, 207, 207, 222, 215, 222, 221, 221, 159, 163, 187, 189, 169, 177, 213, 206, 207, 176, 176, 209, 223, 174, 387, 209, 185, 204, 207, 207, 169, 189, 347, 384, 406, 406, 407, 407, 408, 406, 406, 406, 406, 406, 405, 397, 397, 396, 393, 201, 249, 249, 253, 252, 253, 254, 253, 133, 132, 114, 117, 117, 127, 134, 129, 132, 187, 213, 221, 236, 244, 246, 273, 272, 273, 263, 270, 281, 281, 285, 285, 285, 291, 291, 286, 286, 286, 280, 274, 270, 255, 252, 252, 249, 898, 899, 899, 891, 908, 905, 909, 906, 902, 486, 486, 490, 500, 547, 563, 589, 589, 583, 573, 572, 565, 554, 549, 548, 545, 782, 781, 528, 774, 482, 479, 472, 471, 469, 425, 407, 366, 357, 357, 357, 357])
ZPHI - Method
In the zphi function the PHIDP is recalcuted and also the KDP. The specific attenuation (AH_ZPHI) is generated and can also be used for the precipitation estimation.
[3]:
ds_zphi=wrf.attenuation_corr.zphi(ds_clutter_corr.DBZH_no_clutter,
cphase=cphase,
alphax=0.28,
bx=0.78)
ds_zphi
[3]:
<xarray.Dataset>
Dimensions: (azimuth: 720, range: 936)
Coordinates:
* azimuth (azimuth) float64 0.25 0.75 1.25 1.75 ... 358.8 359.2 359.8
elevation (azimuth) float64 0.5 0.5 0.5 0.5 0.5 ... 0.5 0.5 0.5 0.5 0.5
rtime (azimuth) datetime64[ns] 2022-02-16T13:45:28.941817600 ... ...
* range (range) float32 37.5 112.5 187.5 ... 7.009e+04 7.016e+04
time datetime64[ns] 2022-02-16T13:45:01
sweep_mode <U20 'azimuth_surveillance'
longitude float64 13.24
latitude float64 53.55
altitude float64 38.0
Data variables:
PHIDP_RECALC (azimuth, range) float32 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
AH_ZPHI (azimuth, range) float32 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
KDP_RECALC (azimuth, range) float32 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0xarray.Dataset
- azimuth: 720
- range: 936
- azimuth(azimuth)float640.25 0.75 1.25 ... 359.2 359.8
- standard_name :
- ray_azimuth_angle
- long_name :
- azimuth_angle_from_true_north
- units :
- degrees
- axis :
- radial_azimuth_coordinate
- a1gate :
- 1
- angle_res :
- 0.5
array([2.5000e-01, 7.5000e-01, 1.2500e+00, ..., 3.5875e+02, 3.5925e+02, 3.5975e+02]) - elevation(azimuth)float640.5 0.5 0.5 0.5 ... 0.5 0.5 0.5 0.5
- standard_name :
- ray_elevation_angle
- long_name :
- elevation_angle_from_horizontal_plane
- units :
- degrees
- axis :
- radial_elevation_coordinate
array([0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, ... 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5]) - rtime(azimuth)datetime64[ns]2022-02-16T13:45:28.941817600 .....
- standard_name :
- time
array(['2022-02-16T13:45:28.941817600', '2022-02-16T13:45:28.980652544', '2022-02-16T13:45:01.019417600', '2022-02-16T13:45:01.058252544', '2022-02-16T13:45:01.097087744', '2022-02-16T13:45:01.135922688', '2022-02-16T13:45:01.174757632', '2022-02-16T13:45:01.213592832', '2022-02-16T13:45:01.252427776', '2022-02-16T13:45:01.291262976', '2022-02-16T13:45:01.330097920', '2022-02-16T13:45:01.368932864', '2022-02-16T13:45:01.407768064', '2022-02-16T13:45:01.446603008', '2022-02-16T13:45:01.485438208', '2022-02-16T13:45:01.524273152', '2022-02-16T13:45:01.563108096', '2022-02-16T13:45:01.601943296', '2022-02-16T13:45:01.640778240', '2022-02-16T13:45:01.679613440', '2022-02-16T13:45:01.718448384', '2022-02-16T13:45:01.757283328', '2022-02-16T13:45:01.796118528', '2022-02-16T13:45:01.834953472', '2022-02-16T13:45:01.873788672', '2022-02-16T13:45:01.912623616', '2022-02-16T13:45:01.951458816', '2022-02-16T13:45:01.990293760', '2022-02-16T13:45:02.029128704', '2022-02-16T13:45:02.067963904', '2022-02-16T13:45:02.106798848', '2022-02-16T13:45:02.145634048', '2022-02-16T13:45:02.184468992', '2022-02-16T13:45:02.223303936', '2022-02-16T13:45:02.262139136', '2022-02-16T13:45:02.300974080', '2022-02-16T13:45:02.339809280', '2022-02-16T13:45:02.378644224', '2022-02-16T13:45:02.417479168', '2022-02-16T13:45:02.456314368', ... '2022-02-16T13:45:27.466085632', '2022-02-16T13:45:27.504920832', '2022-02-16T13:45:27.543755776', '2022-02-16T13:45:27.582590720', '2022-02-16T13:45:27.621425920', '2022-02-16T13:45:27.660260864', '2022-02-16T13:45:27.699096064', '2022-02-16T13:45:27.737931008', '2022-02-16T13:45:27.776765952', '2022-02-16T13:45:27.815601152', '2022-02-16T13:45:27.854436096', '2022-02-16T13:45:27.893271296', '2022-02-16T13:45:27.932106240', '2022-02-16T13:45:27.970941184', '2022-02-16T13:45:28.009776384', '2022-02-16T13:45:28.048611328', '2022-02-16T13:45:28.087446528', '2022-02-16T13:45:28.126281472', '2022-02-16T13:45:28.165116672', '2022-02-16T13:45:28.203951616', '2022-02-16T13:45:28.242786560', '2022-02-16T13:45:28.281621760', '2022-02-16T13:45:28.320456704', '2022-02-16T13:45:28.359291904', '2022-02-16T13:45:28.398126848', '2022-02-16T13:45:28.436961792', '2022-02-16T13:45:28.475796992', '2022-02-16T13:45:28.514631936', '2022-02-16T13:45:28.553467136', '2022-02-16T13:45:28.592302080', '2022-02-16T13:45:28.631137024', '2022-02-16T13:45:28.669972224', '2022-02-16T13:45:28.708807168', '2022-02-16T13:45:28.747642368', '2022-02-16T13:45:28.786477312', '2022-02-16T13:45:28.825312256', '2022-02-16T13:45:28.864147456', '2022-02-16T13:45:28.902982400'], dtype='datetime64[ns]') - range(range)float3237.5 112.5 ... 7.009e+04 7.016e+04
- units :
- meters
- standard_name :
- projection_range_coordinate
- long_name :
- range_to_measurement_volume
- spacing_is_constant :
- true
- axis :
- radial_range_coordinate
- meters_to_center_of_first_gate :
- 37.5
- meters_between_gates :
- 75.0
array([3.75000e+01, 1.12500e+02, 1.87500e+02, ..., 7.00125e+04, 7.00875e+04, 7.01625e+04], dtype=float32) - time()datetime64[ns]2022-02-16T13:45:01
- standard_name :
- time
array('2022-02-16T13:45:01.000000000', dtype='datetime64[ns]') - sweep_mode()<U20'azimuth_surveillance'
array('azimuth_surveillance', dtype='<U20') - longitude()float6413.24
- long_name :
- longitude
- units :
- degrees_east
- standard_name :
- longitude
array(13.24397)
- latitude()float6453.55
- long_name :
- latitude
- units :
- degrees_north
- positive :
- up
- standard_name :
- latitude
array(53.55478)
- altitude()float6438.0
- long_name :
- altitude
- units :
- meters
- standard_name :
- altitude
array(38.)
- PHIDP_RECALC(azimuth, range)float320.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- standard_name :
- radar_differential_phase_hv
- long_name :
- Differential phase HV
- short_name :
- PHIDP
- units :
- degrees
array([[0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 0. , 0. , 0. ], ..., [0. , 0. , 0. , ..., 0. , 0. , 0. ], [0. , 0. , 0. , ..., 1.4429669, 1.4429669, 1.4429669], [0. , 0. , 0. , ..., 0. , 0. , 0. ]], dtype=float32) - AH_ZPHI(azimuth, range)float320.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- units :
- dB/km
- standard_name :
- specific_attenuation_h
- long_name :
- Specific attenuation H
- short_name :
- AH
array([[0., 0., 0., ..., 0., 0., 0.], [0., 0., 0., ..., 0., 0., 0.], [0., 0., 0., ..., 0., 0., 0.], ..., [0., 0., 0., ..., 0., 0., 0.], [0., 0., 0., ..., 0., 0., 0.], [0., 0., 0., ..., 0., 0., 0.]], dtype=float32) - KDP_RECALC(azimuth, range)float320.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- standard_name :
- radar_specific_differential_phase_hv
- long_name :
- Specific differential phase HV
- short_name :
- KDP
- units :
- degrees per kilometer
array([[0., 0., 0., ..., 0., 0., 0.], [0., 0., 0., ..., 0., 0., 0.], [0., 0., 0., ..., 0., 0., 0.], ..., [0., 0., 0., ..., 0., 0., 0.], [0., 0., 0., ..., 0., 0., 0.], [0., 0., 0., ..., 0., 0., 0.]], dtype=float32)
[4]:
import matplotlib.pyplot as plt
# plot PHIDP Recalc and PHIDP raw
# georefenced dataset for plot
ds_zphi_georef=ds_zphi.pipe(wrl.georef.georeference_dataset)
ds_clutter_corr_georef = ds_clutter_corr.pipe(wrl.georef.georeference_dataset)
# plot
fig = plt.figure(figsize=(14, 5))
# first subplot
ax1 = fig.add_subplot(121)
ds_clutter_corr_georef.PHIDP.plot(x="x", y="y", ax=ax1, vmin=0, vmax=40)
t = plt.title(r'Raw PHIDP')
t.set_y(1.1)
# second subplot
ax2 = fig.add_subplot(122)
ds_zphi_georef.PHIDP_RECALC.plot(x="x", y="y", ax=ax2)
t = plt.title(r'Recalculated PHIDP')
t.set_y(1.1)
Attenuation correction
The function attenuation_correction contains the phase processing and the zphi method.
[3]:
file_path="/tests/data/raw/2022/02/16/2006_20220216_134500.h5"
path_to_config_file="/tests/data/test_settings_wrainfo.json"
file_hdf5=wrf.reader.read_single_file(file_path,
path=path_to_config_file,
grp="sweep_0")
display(file_hdf5)
<xarray.Dataset>
Dimensions: (azimuth: 720, range: 936)
Coordinates:
* range (range) float32 37.5 112.5 187.5 ... 7.009e+04 7.016e+04
* azimuth (azimuth) float64 0.25 0.75 1.25 1.75 ... 358.8 359.2 359.8
elevation (azimuth) float64 0.5 0.5 0.5 0.5 0.5 ... 0.5 0.5 0.5 0.5 0.5
time (azimuth) datetime64[ns] ...
sweep_mode <U20 ...
longitude float64 ...
latitude float64 ...
altitude float64 ...
Data variables:
WRAD (azimuth, range) float32 ...
DBZH (azimuth, range) float32 ...
QC_INFO (azimuth, range) uint16 ...
PHIDP (azimuth, range) float32 ...
RATE (azimuth, range) float32 ...
VRAD (azimuth, range) float32 ...
KDP (azimuth, range) float32 ...
RHOHV (azimuth, range) float32 ...
ZDR (azimuth, range) float32 ...[4]:
file_clutter_corr=wrf.attenuation_corr.attenuation_correction(ds=file_hdf5,
moment="DBZH")
display(file_clutter_corr)
<xarray.Dataset>
Dimensions: (azimuth: 720, range: 936)
Coordinates:
* range (range) float32 37.5 112.5 187.5 ... 7.009e+04 7.016e+04
* azimuth (azimuth) float64 0.25 0.75 1.25 1.75 ... 358.8 359.2 359.8
elevation (azimuth) float64 0.5 0.5 0.5 0.5 0.5 ... 0.5 0.5 0.5 0.5 0.5
time (azimuth) datetime64[ns] ...
sweep_mode <U20 ...
longitude float64 ...
latitude float64 ...
altitude float64 ...
Data variables: (12/13)
WRAD (azimuth, range) float32 ...
DBZH (azimuth, range) float32 nan nan nan ... -32.0 -32.0 -32.0
QC_INFO (azimuth, range) uint16 ...
PHIDP (azimuth, range) float32 ...
RATE (azimuth, range) float32 ...
VRAD (azimuth, range) float32 ...
... ...
RHOHV (azimuth, range) float32 nan nan nan ... 0.1494 0.2291 0.4382
ZDR (azimuth, range) float32 ...
DBZH_CORR (azimuth, range) float32 nan nan nan ... -32.0 -32.0 -32.0
PHIDP_RECALC (azimuth, range) float32 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
AH_ZPHI (azimuth, range) float32 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
KDP_RECALC (azimuth, range) float32 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0[5]:
# smoothing of the dataset
ds_smooth = file_clutter_corr.pad(azimuth=(3, 3),
mode="wrap").rolling(azimuth=7,
min_periods=3,
center=True).mean().isel(azimuth=slice(3, -3))
ds_smooth
[5]:
<xarray.Dataset>
Dimensions: (range: 936, azimuth: 720)
Coordinates:
* range (range) float32 37.5 112.5 187.5 ... 7.009e+04 7.016e+04
* azimuth (azimuth) float64 0.25 0.75 1.25 1.75 ... 358.8 359.2 359.8
elevation (azimuth) float64 0.5 0.5 0.5 0.5 0.5 ... 0.5 0.5 0.5 0.5 0.5
time (azimuth) datetime64[ns] 2022-02-16T13:45:28.941817600 ... ...
sweep_mode <U20 ...
longitude float64 ...
latitude float64 ...
altitude float64 ...
Data variables: (12/13)
WRAD (azimuth, range) float32 nan nan nan 2.606 ... nan nan 3.581
DBZH (azimuth, range) float32 nan nan nan ... -16.57 -9.143 -16.86
QC_INFO (azimuth, range) float32 64.0 64.0 64.0 39.43 ... 0.0 0.0 0.0
PHIDP (azimuth, range) float32 nan nan nan ... 51.23 47.41 25.31
RATE (azimuth, range) float32 nan nan nan ... 0.2857 0.4286 0.2857
VRAD (azimuth, range) float32 nan nan nan -5.5 ... nan nan 4.833
... ...
RHOHV (azimuth, range) float32 nan nan nan ... 0.4318 0.4347 0.4816
ZDR (azimuth, range) float32 nan nan nan ... 0.1316 0.2059 -0.1426
DBZH_CORR (azimuth, range) float32 nan nan nan ... -16.57 -9.143 -16.86
PHIDP_RECALC (azimuth, range) float32 0.0 0.0 0.0 ... -1.788e-07 -3.151e-07
AH_ZPHI (azimuth, range) float32 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
KDP_RECALC (azimuth, range) float32 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0[9]:
# plot PHIDP Recalc and PHIDP raw
# georefenced dataset for plot
ds_smooth_georef=ds_smooth.pipe(wrl.georef.georeference)
# plot
fig = plt.figure(figsize=(14, 5))
# first subplot
ax1 = fig.add_subplot(121)
ds_smooth_georef.DBZH.plot(x="x", y="y", ax=ax1, cmap = "jet", vmin=0, vmax=40)
t = plt.title(r'Raw reflectivity')
t.set_y(1.1)
# second subplot
ax2 = fig.add_subplot(122)
ds_smooth_georef.DBZH_CORR.plot(x="x", y="y", ax=ax2, cmap="jet", vmin=0, vmax=40)
t = plt.title(r'Attenuation corrected reflecivity')
t.set_y(1.1)
Library Reference
Seealso
Get more information about the attenuation correction module in the library reference section.